Python中使用osr模块读取和导出地理坐标系统信息的示例

osr模块（OGR Spatial Reference System）是Python中用于读取和导出地理坐标系统（Coordinate Reference System）信息的模块。地理坐标系统是用来描述地球表面上的点的坐标的一种系统，它通过一组参数来定义坐标的表示方法。

首先，我们需要安装osgeo库，它包含了osr模块。可以使用pip命令来安装osgeo库：

pip install osgeo

下面是一个使用osr模块读取和导出地理坐标系统信息的示例：

from osgeo import osr

# 创建一个SpatialReference对象
srs = osr.SpatialReference()

# 设置地理坐标系统的EPSG代码
srs.ImportFromEPSG(4326)

# 输出地理坐标系统的描述信息
print("地理坐标系统描述信息：", srs.ExportToPrettyWkt())

# 输出地理坐标系统的EPSG代码
print("地理坐标系统EPSG代码：", srs.GetAuthorityCode(None))

# 输出地理坐标系统的名称
print("地理坐标系统名称：", srs.GetAttrValue('AUTHORITY', 0))

# 修改坐标系名称
srs.SetAttrValue('AUTHORITY', 0, 'EPSG')
srs.SetAttrValue('AUTHORITY', 1, '4326')

# 输出修改后的地理坐标系统描述信息
print("修改后的地理坐标系统描述信息：", srs.ExportToPrettyWkt())

在上面的示例中，首先创建了一个SpatialReference对象，然后通过ImportFromEPSG方法设置了地理坐标系统的EPSG代码为4326，即WGS 84经纬度坐标系。

然后使用ExportToPrettyWkt方法输出地理坐标系统的描述信息，它返回一个WKT（Well-Known Text）格式的字符串，用来描述地理坐标系统的详细信息。

使用GetAuthorityCode方法可以获取地理坐标系统的EPSG代码，并使用GetAttrValue方法获取地理坐标系统的名称。

接着使用SetAttrValue方法修改地理坐标系统的名称，将其设置为"EPSG:4326"。

最后使用ExportToPrettyWkt方法输出修改后的地理坐标系统的描述信息。

需要注意的是，osr模块还提供了其他方法用于读取和修改地理坐标系统的信息，如ImportFromProj4、ImportFromWkt、ExportToProj4等方法。可以根据不同的需求选择合适的方法来读取和导出地理坐标系统信息。

总结来说，osr模块是Python中用于读取和导出地理坐标系统信息的功能模块，通过它可以方便地获取和修改地理坐标系统的EPSG代码、描述信息和名称等。