我目前正在使用 Python/Numpy 处理地理/GPS 数据(喜欢它!),并且我面临着计算由坐标对定义的地理点之间的距离的重复任务pn = [lon, lat]
.
我有一个这样使用的函数:dist = geodistance(p1, p2)
它类似于线性代数中的欧几里德距离(向量减法/差),但发生在测地(球形)空间而不是矩形欧几里德空间中。
从编程上讲,欧氏距离由下式给出
dist = ((p2[0] - p1[0])**2 + (p2[1] - p1[1])**2)**0.5
从数学上讲,这相当于“惯用”(因为没有更好的词)句子
dist = p1 - p1 # the "norm" of the vector difference, subtraction.
目前,我得到的距离是这样的:
p1 = [-51.598354,-29.953363]
p2 = [-51.598701,-29.953045]
dist = geodistance(p1, p2)
print dist
>> 44.3904032407
我想这样做:
print p2 - p1 # these points now are from some fancy datatype
>> 44.3904032407
最终目标:
track = numpy.array([[-51.203018 -29.996149]
[-51.203018 -29.99625 ]
[-51.20266 -29.996229]
[-51.20229 -29.996309]
[-51.201519 -29.99416 ]], dtype=fancy) # (**) or something like
print numpy.diff(track)
>> ndarray([[ 0. ]
[ 7.03531252]
[ 39.82663316]
[ 41.50958596]
[ 172.49825765]])
类似的事情是:如果你拿两个datetime
对象并减去它们,该操作返回一个timedelta
目的。我想减去两个坐标并得到测地距离作为结果。
我想知道一个类是否可以工作,但是 dtype (例如,float32 的“子类型”)对于从列表创建数组(** 这是我从 xml 文件中读取内容的方式)会有很大帮助。
多谢!