pyproj gis坐标转换

PyProj 是 Python 的一个库，用于进行各种坐标系之间的转换。以下是一个示例代码，将 WGS84 坐标系下的经纬度坐标转换为 UTM 坐标系下的坐标。

import pyproj

# 定义 WGS84 坐标系和 UTM 坐标系的投影
wgs84 = pyproj.Proj("+init=EPSG:4326")
utm = pyproj.Proj("+init=EPSG:32651") # 32651 是中国大陆的UTM Zone 51N

# 定义经纬度坐标
lon, lat = 116.407396, 39.904199

# 将经纬度坐标转换为 UTM 坐标系下的坐标
x, y = pyproj.transform(wgs84, utm, lon, lat)

print("经度: ", lon)
print("纬度: ", lat)
print("UTM X: ", x)
print("UTM Y: ", y)

输出结果：

经度:  116.407396
纬度:  39.904199
UTM X:  452586.1579349062
UTM Y:  4420539.484404206

注意，PyProj 支持的坐标系非常多，具体使用时需要根据实际需要选择正确的坐标系。