- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我正在尝试检查在使用欧氏距离而不是大圆距离 (gcd) 计算地球上两点的距离时引入的错误。我有两个由纬度和经度定义的点。我为 great circle distance 使用了 python geopy 框架.这里是 gcd 的代码:
def measure(self, a, b):
a, b = Point(a), Point(b)
lat1, lng1 = radians(degrees=a.latitude), radians(degrees=a.longitude)
lat2, lng2 = radians(degrees=b.latitude), radians(degrees=b.longitude)
sin_lat1, cos_lat1 = sin(lat1), cos(lat1)
sin_lat2, cos_lat2 = sin(lat2), cos(lat2)
delta_lng = lng2 - lng1
cos_delta_lng, sin_delta_lng = cos(delta_lng), sin(delta_lng)
d = atan2(sqrt((cos_lat2 * sin_delta_lng) ** 2 +
(cos_lat1 * sin_lat2 -
sin_lat1 * cos_lat2 * cos_delta_lng) ** 2),
sin_lat1 * sin_lat2 + cos_lat1 * cos_lat2 * cos_delta_lng)
return self.RADIUS * d
所以还是两点:
p1=[39.8616,-75.0748], p2=[-7.30933,112.76]
gcd = 78.8433004543197 klm
使用来自 geopy 的 great_circle(p1,p2).kilometers
函数
然后我使用以下公式将这两个点转换为笛卡尔坐标:
def spherical_to_cartesian(r,la,lo):
x=r*np.sin(90-la)*np.cos(lo)
y=r*np.sin(90-la)*np.sin(lo)
z=r*np.cos(90-la)
return (x,y,z)
其中 r=6372.795
,这导致以下笛卡尔坐标
p1=[ -765.81579368, -256.69640558, 6321.40405587],
p2=[480.8302149,-168.64726394,-6352.39140142]
然后输入:np.linalg.norm(p2-p1)
我得到 1103.4963114787836
作为他们的欧几里德范数,与来自的 ~78klm 相比似乎不合理gcd。我是不是推断错了?
最佳答案
Python 在 math 包中包含两个函数; radians 将度数转换为弧度,degrees 将弧度转换为度数。
sin() 方法返回 x 的正弦值,以弧度为单位。
import math
def spherical_to_cartesian(r,la,lo):
rlo = math.radians(lo)
rla = math.radians(90-la)
x=r*np.sin(rla)*np.cos(rlo)
y=r*np.sin(rla)*np.sin(rlo)
z=r*np.cos(rla)
return (x,y,z)
关于python - 如何使用python比较大圆距离和两个球面点的欧氏距离?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/22562739/
我在PostGIS数据库(-4326)中使用经纬度/经度SRID。我想以一种有效的方式找到最接近给定点的点。我试图做一个 ORDER BY ST_Distance(point, ST_GeomF
我想在几对上创建一个 Voronoi 图 纬度/经度,但想使用大圆距离 它们之间,而不是(不准确的)毕达哥拉斯距离。 我可以让 qhull/qvoronoi 或其他一些 Linux 程序做到这一点吗?
我是一名优秀的程序员,十分优秀!