Python中OSGeo的osr模块实现投影坐标转换

在Python中，可以使用OSGeo的osr模块来实现投影坐标转换。osr模块提供了一种简单而方便的方式来处理不同坐标系之间的转换。下面是一个使用例子，展示了如何使用osr模块进行投影坐标转换。

# 导入osr模块
from osgeo import osr

# 创建源坐标系EPSG:4326
wgs84 = osr.SpatialReference()
wgs84.ImportFromEPSG(4326)

# 创建目标坐标系EPSG:3857（Web Mercator）
web_mercator = osr.SpatialReference()
web_mercator.ImportFromEPSG(3857)

# 创建投影转换实例
transform = osr.CoordinateTransformation(wgs84, web_mercator)

# 定义源坐标
lon = 120.10
lat = 30.20

# 将源坐标转换为目标坐标
transformed = transform.TransformPoint(lon, lat)

# 打印结果
print("源坐标：", lon, lat)
print("目标坐标：", transformed[0], transformed[1])

在上面的例子中，首先导入osr模块。然后，创建源坐标系和目标坐标系的SpatialReference实例。这里使用的是EPSG代码来定义坐标系，4326代表WGS84坐标系，3857代表Web Mercator投影坐标系。

然后，创建CoordinateTransformation实例，它接受源坐标系和目标坐标系作为参数。在这里，我们将WGS84坐标系转换为Web Mercator投影坐标系。

接下来，我们定义源坐标lon和lat，并使用TransformPoint方法将其转换为目标坐标系中的点。TransformPoint方法接受两个参数，表示经度和纬度。转换后的结果以元组的形式返回， 个元素表示转换后的x坐标，第二个元素表示转换后的y坐标。

最后，打印转换前的源坐标和转换后的目标坐标。

这只是osr模块的一个简单示例，osr模块还提供了更多的方法和功能，用于处理不同的空间参考系统和坐标转换需求。可以参考OSGeo官方文档中的详细说明来了解更多信息。