【PythonGIS】基于Python融合矢量数据（多面合一）

之前发过使用批量合并矢量数据的文章：【Python&GIS】基于Python批量合并矢量数据，正好前段时间有需求把矢量数据进行融合，然后就编了一段融合矢量数据的代码。今天就和大家分享一下如何使用Python对矢量数据实现融合的操作。

1.定义

首先大家要明白矢量数据的合并和融合有什么区别！合并是指将多个shp合并成1个shp，而融合是指将1个shp中的多个要素融合成1个要素。如下所示：

2.代码部分

我们这里使用geopandas和shapely.geometry两个库，使用geopandas内置的函数实现矢量数据的融合。

# -*- coding: utf-8 -*-
"""
@Time ： 2024/1/22 16:09
@Auth ： RS迷途小书童
@File ：Vector Data Integrate.py
@IDE ：PyCharm
@Purpose：融合矢量数据,多面成单面
@Web：博客地址:https://blog.csdn.net/m0_56729804
"""
import geopandas as gpd
from shapely.geometry import Polygon, MultiPolygondef integrate_shp(file1=r"彭俊喜/1.shp", file2=r"彭俊喜/2.shp"):  # shp将多个图层融合成一个图层gdf = gpd.read_file(file1)crs_original = gdf.crs  # 获取坐标系信息# 确保所有的几何对象都是Polygon类型gdf['geometry'] = gdf['geometry'].apply(lambda x: x if isinstance(x, Polygon) else x.convex_hull)# 将所有面合并为一个面merged_polygon = gdf['geometry'].unary_union# 将结果保存到新的shp文件中output_gdf = gpd.GeoDataFrame({'geometry': merged_polygon}, index=[0]).set_crs(crs_original)output_gdf.to_file(file2, driver='ESRI Shapefile')if __name__ == "__main__":path1 = ""  # 输入矢量文件path2 = ""  # 输出融合后的路径integrate_shp(path1, path2)  # 调用函数