罗湖区做网站的公司,浏览wap网站,兄弟们来个能看的,亚马逊公司网站建设的目的之前发过使用批量合并矢量数据的文章#xff1a;【PythonGIS】基于Python批量合并矢量数据#xff0c;正好前段时间有需求把矢量数据进行融合#xff0c;然后就编了一段融合矢量数据的代码。今天就和大家分享一下如何使用Python对矢量数据实现融合的操作。
1.定义 首先… 之前发过使用批量合并矢量数据的文章【PythonGIS】基于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_56729804import geopandas as gpd
from shapely.geometry import Polygon, MultiPolygondef integrate_shp(file1r彭俊喜/1.shp, file2r彭俊喜/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, driverESRI Shapefile)if __name__ __main__:path1 # 输入矢量文件path2 # 输出融合后的路径integrate_shp(path1, path2) # 调用函数
3.总结 融合、合并、创建、修改属性表等都是矢量数据的常见操作同时也是每个GISer必备的技能我会不定期地在博客上分享一些自己在进行RS、GIS工作时使用到的代码以及学习经验。如果大家感兴趣可以点个关注有什么问题可以评论或者私信
