之前发过使用批量合并矢量数据的文章:【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) # 调用函数
3.总结
融合、合并、创建、修改属性表等都是矢量数据的常见操作,同时也是每个GISer必备的技能,我会不定期地在博客上分享一些自己在进行RS、GIS工作时使用到的代码以及学习经验。如果大家感兴趣可以点个关注,有什么问题可以评论或者私信!