java已知两点经纬度，两点连线某一点到任一点的距离，计算该点的经纬度

秒速五厘米 2024-02-05 19:49 49阅读 0赞

计算逻辑大概是这样（我们把两个点分为 点1 和 点2； 计算的点设为中间点）：

1.根据两点的经纬度， 计算点2相对于点1的方位角度

2. 根据点1 的经纬度， 点1到中间点的距离，上面计算的方位角度。来计算中间点的经纬度。

具体如下：

准备实体类Point（也可以直接传经纬度，需要自己稍作修改）

@AllArgsConstructor
    @NoArgsConstructor
    @Data
    public class Point {
    
        private Double lon;
        private Double lat;
    
    
    }

util：

import java.text.DecimalFormat;
    import java.util.ArrayList;
    import java.util.Arrays;
    import java.util.List;
    import java.util.regex.Matcher;
    import java.util.regex.Pattern;
    
    import static java.lang.Math.cos;
    import static java.lang.Math.sin;
    
    @Component
    @Slf4j
    public class JWDUtil {
    
       
    
        /**
         * 默认地球半径
         */
        private static final double EARTH_RADIUS_KM = 6378.2;//赤道半径(单位Km)
    
    
    
    
        /**
         *
         * @param point1  第一个点
         * @param point2  第二个点
         * @param distance  第一个点到中间点的距离
         * @return  double[]  point[0]  计算得到的中间点的纬度    point[1]; // 计算得到的中间点的经度
         */
        public static double[] calculateIntermediatePoint(Point point1, Point point2, double distance) {
    
            double lon1 = point1.getLon();
            double lat1 = point1.getLat();
            double lon2 = point2.getLon();
            double lat2 = point2.getLat();
            double bearing = bearingDistOffset(lon1, lat1, lon2, lat2);
    
    
            double[] doubles = LongLatOffset(lon1, lat1, bearing, distance);
    
            return doubles;
    
        }
    
    
        /**
         * 根据一点的经纬度, 到另外一点的距离, 方位角, 计算另外一点的经纬度
         * @param lon1   经度 以度为单位
         * @param lat1   纬度 以度为单位
         * @param bearing 方位角（以度为单位）
         * @param distance   移动距离 千米
         * @return	返回double数组 double[经度， 纬度]
         */
        public static double[] LongLatOffset(double lon1, double lat1, double bearing, double distance) {
    //
            //  方位角 角度转换到弧度
            lon1 = Math.toRadians(lon1);
            lat1 = Math.toRadians(lat1);
            bearing = Math.toRadians(bearing);
            double lat2 = Math.asin(Math.sin(lat1) * Math.cos(distance / EARTH_RADIUS_KM) +
                    Math.cos(lat1) * Math.sin(distance / EARTH_RADIUS_KM) * Math.cos(bearing));
    
            double lon2 = lon1 + Math.atan2(Math.sin(bearing) * Math.sin(distance / EARTH_RADIUS_KM) * Math.cos(lat1),
                    Math.cos(distance / EARTH_RADIUS_KM) - Math.sin(lat1) * Math.sin(lat2));
    
            return new double[]{Math.toDegrees(lon2), Math.toDegrees(lat2)};
        }
    
        /**
         * 求坐标2相对于坐标1的方位角度和距离
         * 单位：度
         * @param lon1
         * @param lat1
         * @param lon2
         * @param lat2
         * @return	返回 double数组[方位角]，方位角单位: 度 
         */
        public static double bearingDistOffset(double lon1, double lat1, double lon2, double lat2) {
            lon1 = Math.toRadians(lon1);
            lat1 = Math.toRadians(lat1);
            lon2 = Math.toRadians(lon2);
            lat2 = Math.toRadians(lat2);
            // 是坐标1的纬度和经度（以弧度为单位）
            double bearing = Math.atan2(Math.sin(lon2 - lon1) * Math.cos(lat2),
                    Math.cos(lat1) * Math.sin(lat2) - Math.sin(lat1) * Math.cos(lat2) * Math.cos(lon2 - lon1));
            bearing = Math.toDegrees(bearing);
            log.info("计算后方位角：{}°", bearing);
    
            return bearing;
        }
    
    
    
    
    
    
    }

测试调用：

public static void main(String[] args) throws IOException {
    
        double lon1 = 105.3438608712768;
        double lat1 = 28.73352486773717;
    
        double lon2 = 105.3438818347129;
        double lat2 = 28.74322995888552;
    
        Point point1 = new Point(lon1, lat1);
        Point point2 = new Point(lon2, lat2);
        double distance = 0.56;        // 中间点到第一个点的距离（单位km）
        // 计算中间点经纬度
        double[] point = JWDUtil.calculateIntermediatePoint(point1, point2, distance);
    
    
        System.out.println("中间点的经度：" + point[0]);
        System.out.println("中间点的维度：" + point[1]);
    
    
    
    }

亲测有效， 使用方位角计算，偏差可以几乎忽略不计。。。。。。

又是摸鱼的一天。。。