Python ogr融合两个线性矢量后，多段线转换单线

在 OGR 中，可以通过将多个线性矢量合并成一个多段线矢量，然后再将其转换为单线矢量来实现。下面是一个示例代码：

from osgeo import ogr

# 打开第一个线性矢量
ds1 = ogr.Open("path/to/vector1.shp")
lyr1 = ds1.GetLayer()

# 打开第二个线性矢量
ds2 = ogr.Open("path/to/vector2.shp")
lyr2 = ds2.GetLayer()

# 创建一个多段线矢量
multi_line = ogr.Geometry(ogr.wkbMultiLineString)

# 将第一个线性矢量中的所有要素添加到多段线矢量中
for feat in lyr1:
    geom = feat.GetGeometryRef()
    multi_line.AddGeometry(geom)

# 将第二个线性矢量中的所有要素添加到多段线矢量中
for feat in lyr2:
    geom = feat.GetGeometryRef()
    multi_line.AddGeometry(geom)

# 将多段线矢量转换为单线矢量
line = multi_line.GetGeometryRef(0)

# 创建一个新的矢量文件
driver = ogr.GetDriverByName("ESRI Shapefile")
ds_out = driver.CreateDataSource("path/to/output.shp")
lyr_out = ds_out.CreateLayer("output", geom_type=ogr.wkbLineString)

# 将单线矢量写入输出图层中
feat_out = ogr.Feature(lyr_out.GetLayerDefn())
feat_out.SetGeometry(line)
lyr_out.CreateFeature(feat_out)

# 关闭数据源
ds_out = None
ds1 = None
ds2 = None

在上面的代码中，我们首先打开两个线性矢量，并将它们的要素逐个添加到一个多段线矢量中。然后，我们使用 `GetGeometryRef(0)` 方法从多段线矢量中提取单线矢量，并将其写入一个新的矢量文件中。

请注意，如果你想将多个线性矢量合并成一个多边形矢量，也可以使用类似的方法。只需要将 `geom_type` 参数设置为 `ogr.wkbPolygon`，并在创建 `multi_line` 和 `line` 变量时分别改为 `ogr.Geometry(ogr.wkbMultiPolygon)` 和 `ogr.Geometry(ogr.wkbPolygon)` 即可。