两点经纬度，计算距离
这种公式我必然是不知道的，谷歌翻了翻，有人（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;
??? }
两点经纬度，计算方位角
最后在http://bbs.godeyes.cn/Upload/2006/12/06/205631.rar下了一个excel，用它的计算公式写了一个。
计算方位角pab，
其中lat_a, lng_a是A的纬度和经度； lat_b, lng_b是B的纬度和经度。代码如下：
??? // 计算方位角pab。
??? privatedouble gps2d(double lat_a, double lng_a, double lat_b, double lng_b) {
?????? 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
}