2021年全国水体分布（按省、市、县）矢量数据

1264

收藏 2021-12-20

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 intqdm(tifPaths_folder_SHENG):
      try:
         #省目录
         DEM_SHENG_PATH =os.path.join(mask_path, Mask_PATH)
         #获得文件中的名字
         SHENG_PATH=DEM_SHENG_PATH+"\\"+mask_path+"\\"+mask_path+".shp"
         #裁剪函数
         City =gpd.read_file(SHENG_PATH)
         Parking_shp =gpd.read_file(Data_path)
         City=City.to_crs(Parking_shp.crs)
         Beijing_Parking = gpd.clip(Parking_shp,City)
         mkdir(OutPut_path+mask_path+"_2021年全省水体分布_矢量")
         #保存文件SHP
         save_SHP_path = OutPut_path +mask_path + "\\"+ mask_path +"水体分布"+ ".shp"
         Beijing_Parking.to_file(save_SHP_path,driver="ESRI Shapefile",encoding="utf-8")
      except:
         pass

其中，裁剪之后的结果文件保存为shp格式，由于字段包含中文，编码格式一定要选择UTF8。

2.3数据后处理

这一步主要是在已经处理好的水体数据文件夹中，添加水体效果图以及各个行政区的矢量数据。其中，添加水体效果图到每一个文件夹的脚本为：

# 第一部分，准备工作，拼接出要存放的文件夹的路径
file_path = r"待添加文件的路径"
# 文件列表
current_folder = os.listdir(file_path)
# 添加文件的位置
add_file_path1 = "需要添加的文件"

#将名称为file的文件复制到名为file_dir的文件夹中
for list intqdm(current_folder):
# 拼接出要存放的文件夹的路径
dir_path =file_path + list
# 将指定的文件file移动到file_dir的文件夹里面
shutil.copy(add_file_path1, dir_path)

之后，使用一个脚本，将每一个水体文件对应的行政边界添加到对应文件中，下面以添加对应的县行政边界为例：

# 第一部分，准备工作，拼接出要存放的文件夹的路径
file_path = r"待添加文件的路径"
#矢量边界的目录
gsons_path = r"矢量边界的路径"

# 文件列表
current_sheng_folder = os.listdir(file_path)
# 循环省
for sheng intqdm(current_sheng_folder):
# 省份的目录
path_sheng =file_path + sheng+"\\"
#循环市
current_shi_folder= os.listdir(path_sheng)
for shi in current_shi_folder:
      # 市的目录
      path_shi =path_sheng + shi+"\\"
      # 循环县
   current_xian_folder = os.listdir(path_shi)
      for xian incurrent_xian_folder:
         try:
            #行政边界矢量位置
            gson_name=sheng+"_"+shi+"_"+xian+".gson"
            gson_path=gsons_path+gson_name
            #该县的路径
            xian_path=path_shi+xian
            shutil.copy(gson_path,xian_path)
         except:
            pass