python-2.7 - 在 Python 中将纬度、经度和海拔高度转换为本地 ENU 坐标

Question

如何使用 Python 将大地测量(纬度、经度、高度)坐标转换为本地切线平面 ENU(东、北、上)坐标？
pyproj 包似乎没有正确的功能......

Answer 1

您可以使用 pymap3d 包裹:
安装

pip install pymap3d

简单的例子
让我们以 this page 上显示的示例值为例.

import pymap3d as pm

# The local coordinate origin (Zermatt, Switzerland)
lat0 = 46.017 # deg
lon0 = 7.750  # deg
h0 = 1673     # meters

# The point of interest
lat = 45.976  # deg
lon = 7.658   # deg
h = 4531      # meters

pm.geodetic2enu(lat, lon, h, lat0, lon0, h0)

产生

(-7134.757195979863, -4556.321513844541, 2852.3904239436915)

分别是东、北和上分量。
使用的椭球默认为 WGS84。所有可用的椭球模型(可以作为 geodetic2enu 的参数)都可以看到 here .以下是使用 WGS72 引用椭球计算相同 ENU 坐标的方法:

pm.geodetic2enu(lat, lon, h, lat0, lon0, h0, ell=pm.utils.Ellipsoid('wgs72'))
# (-7134.754845247729, -4556.320150825548, 2852.3904257449926)