在Python中使用osgeo.osr模块读取和导出Shapefile文件的地理参考信息

在Python中，可以使用osgeo.osr模块读取和导出Shapefile文件的地理参考信息。osr模块是GDAL（Geospatial Data Abstraction Library）的一部分，用于处理和转换地理坐标系统。

首先，需要安装GDAL库和相关Python绑定。可以使用pip命令安装GDAL库：

pip install gdal

代码示例：

import osgeo.osr as osr

# 读取Shapefile文件
def read_shapefile_metadata(filepath):
    # 创建空间参考对象
    spatial_ref = osr.SpatialReference()
    
    # 打开Shapefile文件
    datasource = osr.Open(filepath)
    
    # 获取      个图层
    layer = datasource.GetLayerByIndex(0)
    
    # 获取图层的空间参考
    spatial_ref.ImportFromWkt(layer.GetSpatialRef().ExportToWkt())
    
    # 打印空间参考信息
    print("空间参考信息：")
    print("坐标系名称：", spatial_ref.GetAttrValue("AUTHORITY", 1))
    print("地理坐标系：", spatial_ref.GetAttrValue("GEOGCS"))
    print("投影坐标系：", spatial_ref.GetAttrValue("PROJCS"))
    

# 导出Shapefile文件的地理参考信息
def export_shapefile_metadata(filepath, dest_filepath):
    # 创建空间参考对象
    spatial_ref = osr.SpatialReference()
    
    # 打开Shapefile文件
    datasource = osr.Open(filepath)
    
    # 获取      层图层
    layer = datasource.GetLayerByIndex(0)
    
    # 获取图层的空间参考
    spatial_ref.ImportFromWkt(layer.GetSpatialRef().ExportToWkt())
    
    # 创建目标文件，写入地理参考信息
    with open(dest_filepath, 'w') as file:
        file.write(spatial_ref.ExportToPrettyWkt())
    
    print("地理参考信息导出成功。")

使用示例：

# 读取Shapefile文件的地理参考信息
read_shapefile_metadata("shapefile.shp")

# 导出Shapefile文件的地理参考信息
export_shapefile_metadata("shapefile.shp", "metadata.txt")

在上面的示例中，read_shapefile_metadata函数用于读取Shapefile文件的地理参考信息并打印到控制台。export_shapefile_metadata函数用于导出Shapefile文件的地理参考信息到指定的文件。需要将"shapefile.shp"替换为实际的Shapefile文件路径，"metadata.txt"替换为地理参考信息导出的目标文件路径。

使用osgeo.osr模块，可以轻松读取和导出Shapefile文件的地理参考信息，以便进行地理数据的处理和分析。