标签:python latitude-longitude gps numpy user-defined-types
我目前正在使用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文件中读取内容的方式).
非常感谢!
解决方法:
您可以通过创建类并编写__add__或__sub__方法来定义自己的类型.
例如:
class P(object):
def __init__(self, lon, lat):
self.lon = lon
self.lat = lat
def __sub__(self, other):
dist = ((other.lon - self.lon)**2 + (other.lat - self.lat)**2)**0.5
return dist
鉴于您目前正在使用列表索引语法获取点的坐标,您还可以实现这些:
class P(object):
def __init__(self, lon, lat):
self.lon = lon
self.lat = lat
def __sub__(self, other):
dist = ((other[0] - self[0])**2 + (other[1] - self[1])**2)**0.5
return dist
def __getitem__(self, key):
if key == 0:
return self.lon
elif key == 1:
return self.lat
else:
raise IndexError
def __setitem__(self, key, value):
if key == 0:
self.lon = value
elif key == 1:
self.lat = value
else:
raise IndexError
(我意识到上述可能不是最优雅的方式).
这样,您的新课程就是您当前使用的列表的替代品.
Python documentation包含有关您需要编写的双下划线方法的更多信息,以便创建用户定义的类型. (您正在寻找的信息大约在页面的一半处开始)
标签:python,latitude-longitude,gps,numpy,user-defined-types 来源: https://codeday.me/bug/20190620/1247299.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。