【Python&GIS】基于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, MultiPolygon


def 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)  # 调用函数

3.总结

        融合、合并、创建、修改属性表等都是矢量数据的常见操作，同时也是每个GISer必备的技能，我会不定期地在博客上分享一些自己在进行RS、GIS工作时使用到的代码以及学习经验。如果大家感兴趣可以点个关注，有什么问题可以评论或者私信！