### 将ECEF坐标转换为地理坐标（经纬度和高度） 为了实现从地心地球固定(ECEF)坐标到地理坐标（纬度、经度和高度）的转换，在 Python 中可以采用多种方式。一种常见的方式是利用已有的库，如 `pyproj` 或者编写自定义函数来完成这一过程。 #### 使用 pyproj 库 `pyproj` 是一个强大的投影变换库，支持多种坐标系统的相互转换。对于 ECEF 到 LLA 的转换，可以通过如下代码片段实现： ```python import pyproj def ecef_to_lla(x, y, z): transformer = pyproj.Transformer.from_crs( {"proj": "geocent", "ellps": "WGS84", "datum": "WGS84"}, {"proj": "latlong", "ellps": "WGS84", "datum": "WGS84"}, ) lon, lat, alt = transformer.transform(x, y, z, radians=False) return lat, lon, alt ``` 此段代码创建了一个从 ECEF (地心直角坐标系) 转换至 LLA (地理坐标系) 的转换器，并返回给定 ECEF 坐标的相应纬度、经度以及海拔高度[^2]。 #### 自定义实现算法 如果不希望依赖外部库，则可以根据 WGS84 椭球模型的相关公式自行编码实现转换逻辑。下面是一个简单的例子: ```python from math import pi, sqrt, atan, pow a = 6378137.0 # 半径单位米 f = 1 / 298.257223563 # 扁率 b = a * (1 - f) def ecef_to_geodetic(x, y, z): eps = b*b/(a*a) p = sqrt(x*x+y*y) q = atan((z*eps)/(p)) sinQ = sin(q) cosQ = cos(q) Ro = sqrt(1-(2*f-f*f)*sinQ*sinQ) while abs(z-b*sinQ)>1e-6: sinQOld = sinQ N = a/Ro h = (p/cosQ)-N phi = atan((z+(N*(f+f*f*sinQ*sinQ))*sinQ)/p) sinPhi = sin(phi) cosPhi = cos(phi) Ro = sqrt(1-(2*f-f*f)*sinPhi*sinPhi) N = a/Ro h = (p/cosPhi)-N q = atan((z*eps)/(p)) sinQ = sin(q) if(abs(sinQ-sinQOld)<1e-6): break lambda_ = atan(y/x) latitude = degrees(phi) longitude = degrees(lambda_) altitude = h return latitude, longitude, altitude ``` 这段代码实现了基于迭代法求解大地高程的方法，最终得到纬度、经度及高度值[^1]。