如何基于(x,y)获取所有点经纬度？pyproj支持该逆向转换吗？

一、能否基于现有(x,y)点获取所有点的经纬度？ #

完全可以，但需要明确两个核心前提：

• 你的(x,y)属于平面投影坐标系（比如UTM、矿区独立坐标系、自定义局部平面坐标系），经纬度属于地理坐标系（比如WGS84），两者存在确定的数学转换关系
• 现有1/7带经纬度的点是关键控制点，可用于确定转换参数或验证转换模型

具体分两种场景：

1. 已知标准坐标系：如果明确(x,y)和经纬度对应的标准坐标系（比如(x,y)是UTM 18N，经纬度是WGS84），直接用标准转换参数就能批量转换所有点，无需依赖控制点
2. 自定义/未知坐标系：如果(x,y)是自定义局部坐标系，可利用这批控制点拟合转换模型（如四参数、七参数、仿射变换等），拟合完成后即可将剩余(x,y)点转换为经纬度

二、pyproj是否支持此类逆向经纬度转换？ #

pyproj是Python生态中专业的坐标转换库，完全支持这类需求，覆盖上述两种场景：

1. 标准坐标系间的直接转换 #

如果已知源（x,y的平面坐标系）和目标（经纬度的地理坐标系）的EPSG代码或参数，可直接用Transformer完成正向/逆向转换：

from pyproj import Transformer

# 示例：假设(x,y)是UTM 18N（EPSG:32618），目标为WGS84经纬度（EPSG:4326）
transformer = Transformer.from_crs("EPSG:32618", "EPSG:4326", always_xy=True)

# 转换示例中的点
lon, lat = transformer.transform(60.682697, 35.01935)
# 输出结果应接近示例给出的-66.230895, 18.346775

2. 自定义转换参数的拟合与应用 #

如果是自定义坐标系，可通过控制点计算最优转换参数，再应用到批量转换：

from pyproj import Transformer
from pyproj.transformer import TransformerGroup
from pyproj.aoi import AreaOfInterest

# 整理你的控制点列表：(x, y, 经度, 纬度)
control_points = [
    (60.682697, 35.01935, -66.230895, 18.346775),
    # 补充更多控制点以提升拟合精度
]

# 定义源（自定义平面坐标系，可先假设相近的标准CRS）和目标坐标系
src_crs = "EPSG:32618"  # 示例，若未知可通过AOI和控制点自动匹配
tgt_crs = "EPSG:4326"

# 基于控制点和区域范围获取最优转换器
aoi = AreaOfInterest(west_lon=-67, east_lon=-66, south_lat=18, north_lat=19)
tg = TransformerGroup(src_crs, tgt_crs)
best_transformer = tgt_crs.best_transformer(aoi, control_points=control_points)

# 转换任意(x,y)点
x, y = 120.5, 40.2
lon, lat = best_transformer.transform(x, y)

注意事项 #

• 控制点数量越多、分布越均匀，拟合的转换模型精度越高
• 如果转换结果误差较大，需排查坐标系是否匹配，或尝试更复杂的转换模型（如七参数）

内容的提问来源于stack exchange，提问作者Minh Huynh