如何使用Python代码将shapefile（shp）格式的河流水系数据转换为geojson格式？请提供转换过程中的关键步骤和代码示例。

在地理信息系统（GIS）数据处理中，将shapefile（shp）格式转换为geojson格式是一个常见需求，尤其当需要在Web环境中使用这些地理数据时。根据提供的资源《四川河流矢量图层2024版及shp转geojson代码》，这里将介绍使用Python进行shp到geojson转换的关键步骤和示例代码。

参考资源链接：[四川河流矢量图层2024版及shp转geojson代码](https://wenku.csdn.net/doc/4act433y4y?spm=1055.2569.3001.10343)

首先，为了完成这个任务，需要安装GDAL/OGR库，这是一个用于读写多种格式矢量空间数据的开源工具包。你可以通过pip安装GDAL的Python绑定：

pip install GDAL

接着，使用GDAL/OGR库中的功能，我们可以创建一个Python脚本来读取shp文件并将其转换为geojson格式。以下是一个简单的代码示例：

from osgeo import ogr

def shp_to_geojson(input_shp, output_geojson):
    driver = ogr.GetDriverByName('ESRI Shapefile')
    source_ds = driver.Open(input_shp, 0)  # 0 means read-only
    if source_ds is None:
        raise RuntimeError('无法打开文件: ' + input_shp)
    
    source_layer = source_ds.GetLayer()
    geojson_driver = ogr.GetDriverByName('GeoJSON')
    if geojson_driver:
        if os.path.exists(output_geojson):
            os.remove(output_geojson)
        out_ds = geojson_driver.CreateDataSource(output_geojson)
        out_layer = out_ds.CreateLayer('water_system', geom_type=ogr.wkbMultiPolygon)
        source_layer_defn = source_layer.GetLayerDefn()
        for i in range(source_layer_defn.GetFieldCount()):
            field_defn = source_layer_defn.GetFieldDefn(i)
            out_layer.CreateField(field_defn)
        feature = source_layer.GetNextFeature()
        while feature:
            out_layer.CreateFeature(feature)
            feature = source_layer.GetNextFeature()
        out_ds = None
    
    source_ds = None

# 使用示例
shp_to_geojson('path_to_your_shp_file/river_system.shp', 'path_to_save_geojson/river_system.geojson')

在这段代码中，我们首先使用GDAL/OGR的驱动程序打开shp文件，然后创建一个新的GeoJSON数据源，并定义一个新图层。之后，我们遍历shp文件中的每一个要素，并将其添加到新的GeoJSON图层中。

完成转换后，你可以使用这个geojson文件在Web GIS应用或者数据分析中，如使用Leaflet或Mapbox等Web地图服务。

根据你的具体需求，还可以对代码进行修改和扩展，比如添加坐标转换、属性过滤等功能。如果你需要更深入地了解shp数据处理和转换的具体细节，可以参考《四川河流矢量图层2024版及shp转geojson代码》中的完整示例和解释，这将有助于你更好地利用GIS数据进行空间分析和开发工作。

参考资源链接：[四川河流矢量图层2024版及shp转geojson代码](https://wenku.csdn.net/doc/4act433y4y?spm=1055.2569.3001.10343)