直角坐标转化为大地坐标python代码

在Python中，将直角坐标（通常指平面坐标系中的经纬度）转换为大地坐标（如WGS84地理坐标系统下的长轴椭球坐标），可以使用一些开源库，比如`geopy`或者`pyproj`。这里我们以`pyproj`为例，它是处理地理投影转换的标准库。

首先，你需要安装`pyproj`，可以使用pip命令：

pip install pyproj

然后你可以编写一个简单的函数来转换坐标：

from pyproj import Proj, transform

def rect_to_ellipsoidal(longitude, latitude):
    # 初始化两个坐标系，一个是直角坐标（通常是UTM或其他平面坐标系统），另一个是WGS84大地坐标
    in_proj = Proj(proj='utm', zone=your_zone_number, ellps='WGS84')  # 替换your_zone_number为实际的UTM区号
    out_proj = Proj(proj='latlong', datum='WGS84')

    # 进行转换
    x, y = transform(in_proj, out_proj, longitude, latitude)

    return x, y

# 使用示例
longitude = 120.123456  # 经度
latitude = 31.000000   # 纬度
ellipsoidal_longitude, ellipsoidal_latitude = rect_to_ellipsoidal(longitude, latitude)
print(f"转换后的大地坐标：({ellipsoidal_longitude:.6f}, {ellipsoidal_latitude:.6f})")