根据两点经纬度计算地球表面两点间的距离

时间:2022-03-04 19:52:22

现在基于地理信息的应用很火,无论是从google map上获取的地理位置信息,还是通过iphone或Android这类拥有GPS模块的移动终端获取的地理位置信息,都是十进制的经纬度信息。在我们做实际应用时可能需要计算在地球表面某两点间的距离,这里以JavaScript为示例语言,提供计算方法,这也是google map所采用的计算方法。

1. 十进制经纬度转度、分、秒

view plaincopy to clipboardprint?
  1. function dec2deg(dec){   
  2.   var dec = Math.abs(dec) + "";   
  3.   dec = dec.split(".");   
  4.   
  5.   var deg = dec[0];   
  6.   
  7.   dec[1] = "0."+ dec[1];   
  8.   var min_sec = dec[1] * 3600;   
  9.   var min = Math.floor( min_sec / 60 );   
  10.   var sec =  ( min_sec - ( min * 60 ) ) ;   
  11.   
  12.   return [deg, min, sec];   
  13. }  

2. 度、分、秒转十进制经纬度

view plaincopy to clipboardprint?
  1. function deg2dec(deg, min, sec){   
  2.   var deg = Math.abs(deg);   
  3.   var min = Math.abs(min);   
  4.   var sec = Math.abs(sec);   
  5.   return deg * 1 +  ( sec  * 1 + min * 60 )/3600 ;   
  6. }  

3. 计算两点间距离

view plaincopy to clipboardprint?
  1. // 弧度转换   
  2. function rad(d){   
  3.   return d * Math.PI / 180.0;   
  4. }   
  5.   
  6. // 计算距离,结果的单位为千米(km)   
  7. function calDistance(lat1, lng1, lat2, lng2){   
  8.   if( ( Math.abs( lat1 ) > 90  ) ||(  Math.abs( lat2 ) > 90 ) )   
  9.     return false;   
  10.   
  11.   if( ( Math.abs( lng1 ) > 180  ) ||(  Math.abs( lng2 ) > 180 ) )   
  12.     return false;   
  13.   
  14.   var radLat1 = rad(lat1);   
  15.   var radLat2 = rad(lat2);   
  16.   var a = radLat1 - radLat2;   
  17.   var b = rad(lng1) - rad(lng2);   
  18.   var s = 2 * Math.asin(   
  19.       Math.sqrt(   
  20.         Math.powMath.sin( a/2 ), 2 ) + Math.cos( radLat1 ) * Math.cos( radLat2 ) *   
  21.         Math.powMath.sin( b/2 ), 2 )   
  22.       )   
  23.   );   
  24.   s = s * 6378.137 ; // 地球半径 6378.137   
  25.   s = Math.round(s * 10000) / 10000;   
  26.   return s;   
  27. }