四种方法获取地球上两点的真实地面距离

叫我锐多宝

已于 2022-02-09 14:51:18 修改

阅读量1.2k

点赞数 1

分类专栏：锐多宝的地理空间文章标签：地理信息系统投影地面距离两点间距离

于 2021-09-26 23:29:38 首次发布

本文为博主原创文章，未经博主允许不得转载。

本文链接：https://blog.csdn.net/weixin_36396470/article/details/120499361

版权

锐多宝的地理空间专栏收录该内容

120 篇文章

订阅专栏

1.背景

使用非等距离投影，无法获取较大空间距离的两点真实地面距离。
特别是我们经常使用的墨卡托投影是等角投影，特性是水平方向形变必须与竖直方向形变一致。

而如果没有投影坐标，使用arcgis等地理软件，无法计算真实地面距离。

2.四种处理方法

2.1 公式计算法

d＝111.12*cos{1/[sin ΦAsin ΦB 十 cos ΦAcos ΦBcos( λB—λA)]}

其中A点经度，纬度分别为λA和ΦA，B点的经度、纬度分别为λB和ΦB，d为距离。
111.12的来源是赤道周长4万公里，则 40000 公里/360 度＝111.11公里/度(地球是非常近似球体的)。
附上一段JavaScript的代码:

function degreesToRadians(degrees) {
  return degrees * Math.PI / 180;
}

function distanceInKmBetweenEarthCoordinates(lat1, lon1, lat2, lon2) {
  var earthRadiusKm = 6371;

  var dLat = degreesToRadians(lat2-lat1);
  var dLon = degreesToRadians(lon2-lon1);

  lat1 = degreesToRadians(lat1);
  lat2 = degreesToRadians(lat2);

  var a = Math.sin(dLat/2) * Math.sin(dLat/2) +
          Math.sin(dLon/2) * Math.sin(dLon/2) * Math.cos(lat1) * Math.cos(lat2); 
  var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); 
  return earthRadiusKm * c;
}