如何在Python中将浙江省河流矢量图层的SHP数据转换为GeoJSON格式？请提供详细的转换步骤和代码示例。

为了有效地处理地理信息系统（GIS）中的数据，转换矢量图层格式是一项基本技能。特别是当你拥有详细的浙江省河流水系矢量图层数据时，能够将SHP数据转换为GeoJSON格式将大大扩展你的数据可用性和兼容性。以下是在Python中进行该转换的详细步骤和代码示例：

参考资源链接：[2024浙江省河流矢量数据及转换代码免费获取](https://wenku.csdn.net/doc/5h9ovwhznr?spm=1055.2569.3001.10343)

首先，确保你已经安装了GDAL库，这是一个在GIS数据处理中广泛使用的库，特别是用于读取和写入不同格式的矢量数据。你可以使用pip安装GDAL库：

pip install GDAL

接下来，使用Python的GDAL/OGR接口进行数据转换。以下是一个示例代码，展示了如何加载一个SHP文件并将其保存为GeoJSON格式：

from osgeo import ogr, osr

# 打开SHP文件
driver = ogr.GetDriverByName('ESRI Shapefile')
dataSource = driver.Open('zhejiang_rivers.shp', 0)
if dataSource is None:
    raise RuntimeError('无法打开文件')

# 获取图层
layer = dataSource.GetLayer()

# 创建GeoJSON驱动
geojson_driver = ogr.GetDriverByName('GeoJSON')
geojson_out = geojson_driver.CreateDataSource('zhejiang_rivers.geojson')

# 由于创建一个新的数据源，需要先创建图层
geojson_layer = geojson_out.CreateLayer(layer.GetName(), layer.GetSpatialRef(), layer.GetGeomType())

# 复制字段信息
layer_defn = layer.GetLayerDefn()
for i in range(layer_defn.GetFieldCount()):
    field_defn = layer_defn.GetFieldDefn(i)
    geojson_layer.CreateField(field_defn)

# 将特征复制到新的GeoJSON文件
feature = layer.GetNextFeature()
while feature:
    geojson_layer.CreateFeature(feature)
    feature = layer.GetNextFeature()

# 清理
geojson_out.FlushCache()
geojson_out = None
dataSource = None

print(

参考资源链接：[2024浙江省河流矢量数据及转换代码免费获取](https://wenku.csdn.net/doc/5h9ovwhznr?spm=1055.2569.3001.10343)