<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) {
        // 定義兩個地理位置 (經(jīng)度, 緯度)
        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("新地點的經(jīng)緯度: ");
        System.out.println("緯度: " + destination.getLatitude());
        System.out.println("經(jīng)度: " + 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("經(jīng)度: " + boundaryPoint.getLongitude());
        }
    }
}