在去年cosbeta曾经发布了一个网页计算工具，这个作用就是根据地球上两点之间的经纬度计算两点之间的直线距离。经纬度到距离的计算在通信工程中应用比较广泛，所以cosbeta通过搜索找到了一个js的计算脚本（其实是google map的计算脚本，应该算是比较准确了），做成了这个经纬度算距离的工具。
今天有人给cosbeta发邮件，询问计算的公式是什么样的。其实，若是把地球当作一个正常的球体（其实它是椭球）来说，球面两点之间的距离计算并不复杂，运用球坐标很容易就能计算出两点之间的弧长。当然这都是高中的知识，我和你一样，也没有那个耐心来将其推导，所以我就利用google map的经纬度到距离计算的js脚本，将球面弧长的公式给还原出来（估计这个公式是经过部分修正的），还原出来的公式如下：
对上面的公式解释如下：
公式中经纬度均用弧度表示，角度到弧度的转化应该是很简单的了吧，若不会，依然请参考这个这个经纬度算距离的工具；
Lat1 Lung1 表示A点纬度和经度，Lat2 Lung2 表示B点纬度和经度；
a=Lat1 – Lat2 为两点纬度之差 b=Lung1 -Lung2 为两点经度之差；
6378.137为地球半径，单位为公里；
计算出来的结果单位为公里；
相应的js脚本代码为：
//进行经纬度转换为距离的计算
function Rad(d){
return d * Math.PI / 180.0;//经纬度转换成三角函数中度分表形式。
}
//计算距离，参数分别为第一点的纬度，经度；第二点的纬度，经度
function GetDistance(lat1,lng1,lat2,lng2){
var radLat1 = Rad(lat1);
var radLat2 = Rad(lat2);
var a = radLat1 - radLat2;
var b = Rad(lng1) - Rad(lng2);
var 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 *6378.137 ;// EARTH_RADIUS;
s = Math.round(s * 10000) / 10000; //输出为公里
//s=s.toFixed(4);
return s;
}
下面这个网址写的也比较好：
http://www.360doc.com/content/10/1228/10/284485_81937345.shtml