Java使用Geodesy进行地理计算的技术指南

<dependency>
    <groupId>org.gavaghan</groupId>
    <artifactId>geodesy</artifactId>
    <version>1.1.3</version>
</dependency>

import org.gavaghan.geodesy.*;

public class GeodesyExample {
    public static void main(String[] args) {
        // 定义两个地理位置 (经度, 纬度)
        GlobalCoordinates pointA = new GlobalCoordinates(34.052235, -118.243683); // 洛杉矶
        GlobalCoordinates pointB = new GlobalCoordinates(40.712776, -74.005974); // 纽约

        // 使用 WGS84 椭球体
        GeodeticCalculator calculator = new GeodeticCalculator();
        Ellipsoid reference = Ellipsoid.WGS84;

        // 计算距离（米）
        GeodeticCurve curve = calculator.calculateGeodeticCurve(reference, pointA, pointB);
        double distance = curve.getEllipsoidalDistance();

        System.out.println("洛杉矶到纽约的距离: " + distance / 1000 + " 公里");
    }
}

public class BearingExample {
    public static void main(String[] args) {
        GlobalCoordinates pointA = new GlobalCoordinates(34.052235, -118.243683); // 洛杉矶
        GlobalCoordinates pointB = new GlobalCoordinates(40.712776, -74.005974); // 纽约

        GeodeticCalculator calculator = new GeodeticCalculator();
        Ellipsoid reference = Ellipsoid.WGS84;

        // 计算方位角
        GeodeticCurve curve = calculator.calculateGeodeticCurve(reference, pointA, pointB);
        double azimuth = curve.getAzimuth();

        System.out.println("洛杉矶到纽约的方位角: " + azimuth + " 度");
    }
}

public class DestinationExample {
    public static void main(String[] args) {
        GlobalCoordinates startPoint = new GlobalCoordinates(34.052235, -118.243683); // 洛杉矶

        // 距离（米）和方向（角度）
        double distance = 100000; // 100 公里
        double azimuth = 45; // 东北方向

        GeodeticCalculator calculator = new GeodeticCalculator();
        Ellipsoid reference = Ellipsoid.WGS84;

        GlobalCoordinates destination = calculator.calculateEndingGlobalCoordinates(
                reference, startPoint, azimuth, distance);

        System.out.println("新地点的经纬度: ");
        System.out.println("纬度: " + destination.getLatitude());
        System.out.println("经度: " + destination.getLongitude());
    }
}

public class BoundaryExample {
    public static void main(String[] args) {
        GlobalCoordinates center = new GlobalCoordinates(34.052235, -118.243683); // 洛杉矶

        // 距离范围（米）
        double radius = 50000; // 50 公里

        GeodeticCalculator calculator = new GeodeticCalculator();
        Ellipsoid reference = Ellipsoid.WGS84;

        // 计算4个方向的边界点
        double[] azimuths = {0, 90, 180, 270}; // 北、东、南、西
        for (double azimuth : azimuths) {
            GlobalCoordinates boundaryPoint = calculator.calculateEndingGlobalCoordinates(
                    reference, center, azimuth, radius);

            System.out.println("方位角 " + azimuth + " 的边界点:");
            System.out.println("纬度: " + boundaryPoint.getLatitude());
            System.out.println("经度: " + boundaryPoint.getLongitude());
        }
    }
}