计算地理坐标间的距离

只要开始永远不晚

于 2021-12-02 13:17:50 发布

阅读量1k

点赞数

分类专栏：基于位置的服务LBS 文章标签： LBS map python

本文链接：https://blog.csdn.net/haohaizijhz/article/details/121675054

版权

基于位置的服务LBS 专栏收录该内容

4 篇文章 4 订阅

订阅专栏

基础知识

示例代码

tips

计算经纬度间的距离
计算地理坐标间的距离
计算两个经纬度对间的距离
计算两个地理坐标间的距离
计算两点间的距离

基础知识

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

大圆距离（ great_circle ）使用地球的球形模型，使用国际大地测量学和地球物理学联合会定义的平均地球半径，（2）a + b) /3=6371.0087714150598公里，约6371.009公里（WGS-84），误差高达0.5%。半径值存储在 distance.EARTH_RADIUS ，因此可以对其进行自定义（但是，它应该始终以公里为单位）。

The geodesic distance is the shortest distance on the surface of an ellipsoidal model of the earth. The default algorithm uses the method is given by Karney (2013) (geodesic); this is accurate to round-off and always converges.

geopy.distance.distance 目前使用 geodesic .

有许多流行的椭球体模型，哪一个模型最精确取决于你的点在地球上的位置。默认值是WGS-84椭球体，它是全局最精确的。geopy还包括 distance.ELLIPSOIDS 字典：

model major (km) minor (km) flattening
ELLIPSOIDS = {'WGS-84': (6378.137, 6356.7523142, 1 / 298.257223563),
'GRS-80': (6378.137, 6356.7523141, 1 / 298.257222101),
'Airy (1830)': (6377.563396, 6356.256909, 1 / 299.3249646),
'Intl 1924': (6378.388, 6356.911946, 1 / 297.0),
'Clarke (1880)': (6378.249145, 6356.51486955, 1 / 293.465),
'GRS-67': (6378.1600, 6356.774719, 1 / 298.25),
}

geopy文档简介
https://www.osgeo.cn/geopy/#module-geopy.distance

示例代码

使用geopy
代码

# 计算两点距离
from geopy.distance import geodesic
#25.102350,121.548432 # 台北故宫
#25.07639,121.22389 # 台北机场
apoint = (25.102350,121.548432) # 维度在前，经度在后
bpoint = (25.07639,121.22389)
print(geodesic(apoint, bpoint).km) # 默认使用GPS坐标系
print(geodesic(apoint, bpoint).m) # 默认使用GPS坐标系
apoint = (52.2322,21.0083) #华沙
bpoint = (52.5186,13.4081) #柏林
print(geodesic(apoint, bpoint).km)