我特别想将LLA(纬度,经度,海拔)数据转换为ECEF(以地球为中心,固定在地球上)。 (LLA数据为WGS84。)
anaconda的安装程序附带的软件包中是否有此功能?
[如果此包能做的不仅仅是LLA到ECEF,例如反向转换和其他坐标帧(例如ENU),那就很好了。
剩下的只是澄清:
换句话说,假设我在计算机上安装了香草水蟒并且无法安装新软件包。 (Linux,64位,以确保问题的完整性。)
我知道pyproj软件包具有功能,但是根据anaconda网站上的list,它不包含在安装程序中。安装程序中包含很多软件包,我怀疑其中的一个功能存在,但到目前为止,在Stack Overflow,Google和anaconda文档上进行搜索都无法产生答案。
编辑:此外,我知道在LLA和ECEF之间执行单向转换的函数仅需十几行代码,并且该数学是公开可用的。但是,我想发现一个包具有用于坐标转换的功能。 (三个主要原因:(1)进行反向转换以及到其他系统的转换将很有用。(2)我喜欢社区支持的软件,而不必键入和集成自己的程序包。(3)我很好奇。)
from astropy.coordinates import EarthLocation
longitude = -76.899626 # degrees, West (so negative), AKA 283.10037 degrees East
latitude = 39.163002 # degrees
x = EarthLocation.from_geodetic(longitude, latitude)
x.geocentric
Out:
(<Quantity 1122390.43640202 m>,
<Quantity -4823042.93548404 m>,
<Quantity 4006364.03705792 m>)
x.geocentric[0].value # the ECEF x-value
Out: 1122390.4364020245