使用osgeo.osr模块在Python中设置和变换地理坐标系统

osgeo.osr（Open Source Geospatial Foundation's Spatial Reference System）是一个用于地理坐标系统定义、转化和变换的Python模块。该模块提供了一系列函数和方法，可以方便地操作和处理地理坐标系统。

在Python中使用osgeo.osr模块设置和变换地理坐标系统主要包括以下几个步骤：

1. 导入模块：

from osgeo import osr

2. 创建SpatialReference对象：

srs = osr.SpatialReference()

此时，可以使用SpatialReference对象的方法设置和获取地理坐标系统的属性，如坐标系名称、地理坐标系类型、投影坐标系类型等。

srs.SetWellKnownGeogCS("WGS84") # 设置地理坐标系为WGS84

3. 设置或获取当前地理坐标系统的属性：

srs.SetProjCS("UTM Zone 50N") # 设置投影坐标系为UTM Zone 50N
print(srs.GetAttrValue("PROJCS")) # 获取投影坐标系名称

4. 坐标转化：

source_srs = osr.SpatialReference()
target_srs = osr.SpatialReference()
source_srs.ImportFromEPSG(4326) # 设置源地理坐标系统为WGS84
target_srs.ImportFromEPSG(3857) # 设置目标地理坐标系统为Web Mercator

transform = osr.CoordinateTransformation(source_srs, target_srs)
x = 113.55  #源坐标的经度
y = 22.20   #源坐标的纬度
transformed_x, transformed_y, _ = transform.TransformPoint(x, y) #坐标转换
print(transformed_x, transformed_y) #输出坐标转换后的结果

在上述例子中，首先使用ImportFromEPSG方法分别导入了源坐标系（WGS84，EPSG编号4326）和目标坐标系（Web Mercator，EPSG编号3857）。然后创建CoordinateTransformation对象，使用TransformPoint方法进行坐标转换，并输出转换后的结果。

通过以上步骤，就可以实现在Python中使用osgeo.osr模块设置和变换地理坐标系统。该模块提供了丰富的功能和方法，可以满足地理坐标系统的各种需求。