qq_20373723的博客

前面已经有一篇实现分区统计功能的文章了（https://blog.csdn.net/qq_20373723/article/details/111350741），但是核心计算部分是用的numpy，比较难实现majority这个属性，这里用了geopandas，rasterio，gdal地理处理模块组合实现了自动分区统计到表，希望对大家有所帮助python 3.7.3，geopandas 0.8.2 ，GDAL 3.2.1，rasterio 1.2.0，其中的gdal和rasterio是严格对应的，如果不对

也不知道为啥，矢量裁剪矢量这么普通的东西这么难找，赶紧放出来让大家用用import osimport numpy as npimport geopandas as gpdimport warningswarnings.filterwarnings('ignore', 'GeoSeries.notna', UserWarning)shp_path = './osm/gis_osm_roads_free_1.shp' #被裁剪的对象mask = 'D:/ToolSet/Test_data