Mikel
两点经纬度，计算距离
这种公式我必然是不知道的，谷歌翻了翻，有人（http://xxyyyboy.blog.163.com/blog/static/765832620110410457662/）说是
1.Lat1 Lung1 表示A点经纬度，Lat2 Lung2 表示B点经纬度；
2.a=Lat1 – Lat2 为两点纬度之差 b=Lung1 -Lung2 为两点经度之差；
3.6378.137为地球半径，单位为千米；
计算出来的结果单位为千米。
也有人http://panyee.cnblogs.com/archive/2006/07/04/442771.html）说直接从google maps的脚本里扒了段代码,maps的代码：计算的结果是米为单位。
// 计算两点距离
privatefinaldoubleEARTH_RADIUS = 6378137.0;
privatedouble gps2m(double lat_a, double lng_a, double lat_b, double lng_b) {
double radLat1 = (lat_a * Math.PI / 180.0);
double radLat2 = (lat_b * Math.PI / 180.0);
double a = radLat1 – radLat2;
double b = (lng_a – lng_b) * Math.PI / 180.0;
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;
}
两点经纬度，计算方位角
计算方位角pab，
其中lat_a, lng_a是A的纬度和经度； lat_b, lng_b是B的纬度和经度。代码如下：
// 计算方位角pab。
double d = 0;
lat_a=lat_a*Math.PI/180;
lng_a=lng_a*Math.PI/180;
lat_b=lat_b*Math.PI/180;
lng_b=lng_b*Math.PI/180;
d=Math.sin(lat_a)*Math.sin(lat_b)+Math.cos(lat_a)*Math.cos(lat_b)*Math.cos(lng_b-lng_a);
d=Math.sqrt(1-d*d);
d=Math.cos(lat_b)*Math.sin(lng_b-lng_a)/d;
d=Math.asin(d)*180/Math.PI;
// d = Math.round(d*10000);
return d
}
