计算已知经纬度两点的距离

2023-12-27 12:09| 来源: 网络整理| 查看: 265

计算经纬度点之间的距离有三种方法，一种是利用geopy计算；一种是利用haversine公式计算；还可以直接在Excel中计算。

1.直接利用Python第三方库：geopy

代码：

from geopy.distance import geodesic dist = geodesic((37.193462, 100.1347), (37.194363, 100.1348)).km print(dist)

更新：利用geopy计算距离

官网资料：https://www.osgeo.cn/geopy/#module-geopy.distance

Geopy可以用 geodesic distance 或 great-circle distance ，默认测地线距离作为函数 geopy.distance.distance .

classgeopy.distance.geodesic(*args, **kwargs)

用于计算两点间的测地线距离。

参考资料： https://blog.csdn.net/study_in/article/details/86136504

2.半正矢公式(haversine equation)用于计算两经纬度点的距离在这里插入图片描述 https://www.geeksforgeeks.org/haversine-formula-to-find-distance-between-two-points-on-a-sphere/

具体代码：

import math # Python 3 program for the # haversine formula def haversine(lat1, lon1, lat2, lon2): # distance between latitudes # and longitudes dLat = (lat2 - lat1) * math.pi / 180.0 dLon = (lon2 - lon1) * math.pi / 180.0 # convert to radians lat1 = (lat1) * math.pi / 180.0 lat2 = (lat2) * math.pi / 180.0 # apply formulae a = (pow(math.sin(dLat / 2), 2) + pow(math.sin(dLon / 2), 2) * math.cos(lat1) * math.cos(lat2)); rad = 6371 c = 2 * math.asin(math.sqrt(a)) return rad * c # Driver code if __name__ == "__main__": lat1 = 35.289448 lon1 = 79.79672 lat2 = 35.290348 lon2 = 79.796616 print(haversine(lat1, lon1,lat2, lon2), "K.M.")

3.在Excel中计算经纬度两点的距离 http://blog.sciencenet.cn/blog-1094854-1025358.html

计算A（lonA，latA），B（lonB,latB）两点的距离：公式1：

输入公式=6371004*ACOS(1-(POWER((SIN((90-latA)*PI()/180)*COS(lonA*PI()/180)-SIN((90-latB)*PI()/180)*COS(lonB*PI()/180)),2)+POWER((SIN((90-latA)*PI()/180)*SIN(lonA*PI()/180)-SIN((90-latB)*PI()/180)*SIN(lonB*PI()/180)),2)+POWER((COS((90-latA)*PI()/180)-COS((90-latB)*PI()/180)),2))/2)

:公式2：

输入公式=6371004*ACOS((SIN(RADIANS(latA))*SIN(RADIANS(latB))+COS(RADIANS(latA))*COS(RADIANS(latB))*COS(RADIANS(lonB-lonA))))

4.三种方法计算结果对比在这里插入图片描述 Excel中计算的结果（公式2）： 100.5196888m

更新于2021-06-02

【本文地址】

公司简介

联系我们