package xxTest;


public class Test {

    // 地球半徑
    public static final double EARTH_RADIUS = 6378137.0; // 單位M

    private static double getRad(double d) {
        return d * Math.PI / 180.0;
    }

    /**
     * 兩座標（經緯度）之間距離計算
     * GD/BD地圖採用此演算法
     * @param startLng 起點座標經度
     * @param startLat 起點座標緯度
     * @param endLng 終點座標經度
     * @param endLat 終點座標緯度
     * @return 起點座標與終點座標的距離 單位:m
     */
    public static double getGreatCircleDistance(double startLng, double startLat, double endLng, double endLat) {
        double radLat1 = getRad(startLat);
        double radLat2 = getRad(endLat);
        double dy = radLat1 - radLat2; // a
        double dx = getRad(startLng) - getRad(endLng); // b
        double s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(dy / 2), 2) + Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(Math.sin(dx / 2), 2)));
        s = s * EARTH_RADIUS;
        s = Math.round(s * 10000) / 10000.0;
        return s;
    }

    public static void main(String[] args) {
        
        // 東經106.486654度 北緯23.490295度 向東移動4000米 向北移動3000米
        double lng = getMoveLng(23.490295, 4000) + 106.486654;
        double lat = getMoveLat(3000) + 23.490295;
        System.out.println(lat);
        System.out.println(lng);
        System.out.println("移動後兩點距離（米）："+getGreatCircleDistance(106.486654, 23.490295,lng,lat));
    }

    /**
     * 獲取緯度移動latMovement米，緯度移動的度數
     * @param latMovement 單位：米
     * @return 緯度移動的度數
     */
    public static double getMoveLat(double latMovement) {
        return (180/(Math.PI * EARTH_RADIUS))*latMovement;
    }
    
    /**
     * 獲取經度移動lngMovement米，經度移動的度數
     * @param lat 緯度
     * @param lngMovement 單位：米
     * @return 經度移動的度數
     */
    public static double getMoveLng(double lat, double lngMovement) {
        return (180/(Math.PI * EARTH_RADIUS*Math.cos(Math.toRadians(lat))))*lngMovement;
    }
}