来自于： http://shuhui018125.blog.163.com/blog/static/771105020118263164820/
算法可参考：http://blog.csdn.net/zhyl8157121/article/details/7961315
有两种公式可以计算距离：
第一种公式：
地球是一个近乎标准的椭球体，它的赤道半径为6378.140千米，极半径为 6356.755千米，平均半径6371.004千米。如果我们假设地球是一个完美的球体，那么它的半径就是地球的平均半径，记为R。如果以0度经线为基 准，那么根据地球表面任意两点的经纬度就可以计算出这两点间的地表距离（这里忽略地球表面地形对计算带来的误差，仅仅是理论上的估算值）。设第一点A的经 纬度为(LonA, LatA)，第二点B的经纬度为(LonB, LatB)，按照0度经线的基准，东经取经度的正值(Longitude)，西经取经度负值(-Longitude)，北纬取90-纬度值(90- Latitude)，南纬取90+纬度值(90+Latitude)，则经过上述处理过后的两点被计为(MLonA, MLatA)和(MLonB, MLatB)。那么根据三角推导，可以得到计算两点距离的如下公式：
C = sin(MLatA)*sin(MLatB)*cos(MLonA-MLonB) + cos(MLatA)*cos(MLatB)（这是标准的形式，在什么地方都可以适用。这里适用的经纬度为度数，即°在excel中用三角函数的时候都是弧度，需要换算）
Distance = R*Arccos(C)*Pi/180 （换算成弧度）
这里，R和Distance单位是相同，如果是采用6371.004千米作为半径，那么Distance就是千米为单位。
如果仅对经度作正负的处理，而不对纬度作90-Latitude(假设都是北半球，南半球只有澳洲具有应用意义)的处理，那么公式将是：
C = sin(LatA)*sin(LatB) + cos(LatA)*cos(LatB)*cos(MLonA-MLonB)（同样这里的都是°）
Distance = R*Arccos(C)*Pi/180 （换算成弧度）
以上通过简单的三角变换就可以推出。
如果三角函数的输入和输出都采用弧度值，那么公式还可以写作：
C = sin(LatA*Pi/180)*sin(LatB*Pi/180) + cos(LatA*Pi/180)*cos(LatB*Pi/180)*cos((MLonA-MLonB)*Pi/180) （这个公式可以直接在excel中计算）
Distance = R*Arccos(C) （就不需要再换算成弧度）
第二种公式：
（在excel中还是要换算成弧度）
Lat1 Lung1 表示A点经纬度，Lat2 Lung2 表示B点经纬度；
a=Lat1 – Lat2 为两点纬度之差 b=Lung1 -Lung2 为两点经度之差；
6378.137为地球半径，单位为公里；
计算出来的结果单位为公里。
//地球半径
public static final double EARTH_RADIUS = 6378.137d;
public double getDistance(double lng1,double lat1,double lng2,double lat2){
//将角度差转变为弧度
double dLat = Math.toRadians(lat2 - lat1);
double dLng = Math.toRadians(lng2 - lng1);
//计算两点间的距离
double a = Math.sin(dLat /2) * Math.sin(dLat /2) + Math.cos(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2))
* Math.sin(dLng /2 ) * Math.sin(dLng /2);
double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
double distance = EARTH_RADIUS * c * 1000.00d;
return distance ;
}
两种方式计算出来的结果是一样的。
两种公式的原理是一样的。