osr.CoordinateTransformation坐标转换结果错误求助（EPSG4326转32630）

你的转换结果错误的核心原因是坐标顺序传入错误：GDAL/OSR的CoordinateTransformation.TransformPoint方法默认遵循EPSG 4326的官方定义，要求传入纬度在前、经度在后的坐标，但你传入的是经度在前、纬度在后，直接导致坐标错位，转换出完全偏离预期的结果。

解决方案 #

有两种方式可以修正这个问题：

方式1：调整传入TransformPoint的参数顺序 #

直接把纬度放在第一个参数位置，经度放在第二个位置，匹配OSR的默认要求：

from osgeo import ogr, osr, gdal

bbox=[-7.84831, 39.71484, -7.73302, 39.81967]

sourceSR = osr.SpatialReference()
sourceSR.ImportFromEPSG(4326)
targetSR = osr.SpatialReference()
targetSR.ImportFromEPSG(32630)
coordTrans = osr.CoordinateTransformation(sourceSR, targetSR)

min_lon, min_lat, max_lon, max_lat = bbox
# 调整为纬度在前、经度在后
minX, minY, _ = coordTrans.TransformPoint(min_lat, min_lon)
maxX, maxY, _ = coordTrans.TransformPoint(max_lat, max_lon)

trans_bbox = [minX, minY, maxX, maxY]
print(trans_bbox)

方式2：强制设置源坐标系的轴顺序为经度优先 #

如果你习惯使用经度在前、纬度在后的顺序，可以显式修改源坐标系的轴映射策略，让它遵循传统GIS的坐标顺序：

from osgeo import ogr, osr, gdal

bbox=[-7.84831, 39.71484, -7.73302, 39.81967]

sourceSR = osr.SpatialReference()
sourceSR.ImportFromEPSG(4326)
# 强制设置轴顺序为经度（X）在前，纬度（Y）在后
sourceSR.SetAxisMappingStrategy(osr.OAMS_TRADITIONAL_GIS_ORDER)

targetSR = osr.SpatialReference()
targetSR.ImportFromEPSG(32630)
coordTrans = osr.CoordinateTransformation(sourceSR, targetSR)

min_lon, min_lat, max_lon, max_lat = bbox
# 按经度、纬度顺序直接传入
minX, minY, _ = coordTrans.TransformPoint(min_lon, min_lat)
maxX, maxY, _ = coordTrans.TransformPoint(max_lon, max_lat)

trans_bbox = [minX, minY, maxX, maxY]
print(trans_bbox)

问题原理说明 #

EPSG 4326（WGS84）的官方坐标定义是纬度（Y轴）在前，经度（X轴）在后，但日常GIS应用和数据使用中，多数场景习惯采用经度在前、纬度在后的顺序。GDAL/OSR库默认严格遵循EPSG的官方标准，因此直接传入经纬度顺序会导致坐标被错误解析，转换出的UTM坐标自然完全偏离预期范围。

内容的提问来源于stack exchange，提问作者Andrea Viñuales Navarro