gdal 图层内相交要素合并

一、整体思路

将图层内所有要素合并要一个集合类型中，集合类型进行内部合并，遍历集合类型找原图层内查询属性信息，进行属性信息补充。

二、代码实现

# pgurl pg数据库连接地址
# lyrname 合并图层名
from osgeo import gdal, ogr, osr
pg_driver = ogr.GetDriverByName("PostgreSQL")
ds = pg_driver.Open(pgurl,gdal.GA_Update)
sr = osr.SpatialReference()
sr.ImportFromEPSG(4490)
newlyr = ds.CreateLayer('unionLayer', sr, ogr.wkbPolygon, ["LAUNDER=NO","GEOMETRY_NAME=geom"])
searchsql = 'select * from "'+ lyrname +'"'
layer = ds.ExecuteSQL(searchsql)
def_feature = layer.GetLayerDefn()
multi = ogr.Geometry(ogr.wkbMultiPolygon)
for i, feature in enumerate(layer):
    feature.geometry().CloseRings() 
    wkt = feature.geometry().ExportToWkt()
    multi.AddGeometryDirectly(ogr.CreateGeometryFromWkt(wkt))
current_union = multi.UnionCascaded()
for geom in current_union:
    new_feature = ogr.Feature(def_feature)
    layer.SetSpatialFilter(geom)
    fieldvalDic = layer.GetNextFeature().items()
    for field in fieldvalDic:
        new_feature.SetField(field,fieldvalDic[field])
    new_feature.SetGeometry(geom)
    newlyr.CreateFeature(new_feature)
ds.FlushCache()

三、效果图

合并前

合并后