根据经纬度计算两点之间的距离

package com.thinkgem.jeesite.modules.business.utils;
/**
 * @author 罗超
 * @create 2018-03-20 上午7:13
 **/
public class Distance {
    //地球平均半径
    private static final double EARTH_RADIUS = 6376984;
    //把经纬度转为度（°）
    private static double rad(double d){
        return d * Math.PI / 180.0;
    }

    /**
     * 根据两点间经纬度坐标（double值），计算两点间距离，单位为米
     * @param lng1
     * @param lat1
     * @param lng2
     * @param lat2
     * @return
     */
    public static double getDistance(double lng1, double lat1, double lng2, double lat2){
        double radLat1 = rad(lat1);
        double radLat2 = rad(lat2);
        double a = radLat1 - radLat2;
        double b = rad(lng1) - rad(lng2);
        double s = 2 * Math.asin(
                Math.sqrt(
                        Math.pow(Math.sin(a/2),2)
                                + Math.cos(radLat1)*Math.cos(radLat2)*Math.pow(Math.sin(b/2),2)
                )
        );
        s = s * EARTH_RADIUS;
        s = Math.round(s * 10000) / 10000;
        return s;
    }

    /**
     * 中国大地原点纬度
     */
    private static double chinaGeodeticOriginLat = 34.3227;
    /**
     * 用椭圆的参数方程
     * x=acosθ=6378*cosθ
     * y=bsinθ=6357*sinθ
     * 将θ=37°带入
     * x=5093.70
     * y=3825.74
     * 点(5093.70 , 3825.74)到原点(0,0)的距离d=√(x^2+y^2) =6370.41 km
     */
    private static double getRadius() {
        double x = 6378140*Math.cos(chinaGeodeticOriginLat);
        double y = 6356755*Math.sin(chinaGeodeticOriginLat);
        return Math.sqrt(x*x+y*y);
    }

    /**
     * test
     * @param args
     */
    public static void main(String[] args){
        double distance = getDistance(121.491909,31.233234,121.411994,31.206134);
        System.out.println("Distance is:"+distance);
        double radius = getRadius();
        System.out.println("Radius is:"+radius);
    }
}