Python ogr融合

Python的`ogr`库可以用于矢量数据的读取、写入、转换和操作。如果要将两个矢量数据集（例如shapefile）融合在一起，可以使用`ogr`库中的`Union()`函数。以下是一个示例代码：

from osgeo import ogr

# 打开第一个shapefile
shp1 = ogr.Open("path/to/shapefile1.shp")
layer1 = shp1.GetLayer()

# 打开第二个shapefile
shp2 = ogr.Open("path/to/shapefile2.shp")
layer2 = shp2.GetLayer()

# 创建新的shapefile
driver = ogr.GetDriverByName("ESRI Shapefile")
new_shp = driver.CreateDataSource("path/to/new_shapefile.shp")

# 创建新的图层
new_layer = new_shp.CreateLayer("new_layer", layer1.GetSpatialRef(), layer1.GetGeomType())

# 添加字段
layer_def = layer1.GetLayerDefn()
for i in range(layer_def.GetFieldCount()):
    field_def = layer_def.GetFieldDefn(i)
    new_layer.CreateField(field_def)

# 将第一个图层的要素添加到新图层
for feature in layer1:
    new_layer.CreateFeature(feature)

# 将第二个图层的要素添加到新图层
for feature in layer2:
    new_layer.CreateFeature(feature)

# 对新图层进行融合（Union）
new_layer = None
new_shp = None

union_shp = driver.Open("path/to/new_shapefile.shp")
union_layer = union_shp.GetLayer()
union_layer.SetSpatialFilter(None)

results = []
for feature in union_layer:
    results.append(feature.Clone())

union_layer = None
union_shp = None

这段代码首先打开两个shapefile，然后创建一个新的shapefile和图层。它复制了第一个图层的字段，并将其添加到新图层中。接下来，它将第一个图层和第二个图层的要素添加到新图层中。最后，它使用`Union()`函数将新图层中的要素融合在一起，并将结果存储在一个列表中。