01.简介
偶然间获得了Open Street Map的2021年全国水体分布Shp矢量数据,大约400M左右。
当我加载到arcgis中,感觉到了明显的卡顿。由于之前写过分省、市、县裁剪矢量的代码,顺手就将水体数据按照2021年的行政边界裁剪了。
该数据分为3个层次结构,分别为省、地级市、县。其中,省级水体31个(不包含HongKong、TaiWan等), 地级市水体366个,县级市水体2847个。每一个文件夹中都包含该省、地级市或者县的水体矢量数据、行政边界以及数据说明。
02.制作方法
2.1制作方法流程图
2.2Python裁剪脚本
主要是用Geopandas读取水体矢量,并按照行政边界的字段名称进行行政边界的裁剪。
def clip_parking(Mask_PATH,Data_path,OutPut_path):
tifPaths_folder_SHENG = os.listdir(Mask_PATH)
for mask_path in tq