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 tqdm(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年全省水体分布_矢量")#保存文件SHPsave_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 in tqdm(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 in tqdm(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 in current_xian_folder:try:#行政边界矢量位置gson_name=sheng+"_"+shi+"_"+xian+".gson"gson_path=gsons_path+gson_name#该县的路径xian_path=path_shi+xianshutil.copy(gson_path,xian_path)except:pass

03.数据展示

通过上述步骤，我们获得了全国各省、地级市、县等行政单位的水体数据。首先展示省级水体数据：

广东省水体分布图

四川省水体分布图

如果你要找县级市的水体分布图，进入"分县"文件夹中，选择对应省份-地级市-县进行选择。我们以四川省资阳市乐至县的水体数据为例，直接选择分县-四川省-资阳市-乐至县，并将其加到Google影像图中，查看水体矢量效果。

04.数据分享

直接公众号回复，"水体"，获得百度云链接（无套路）。

05.写在后面

1.水体数据是我进行全国数据大整合的一小步。

2.之前我提到的分享100期各省、地级市、县的土地利用数据，目前已经处理好了20期500m分辨率的MODIS土地利用数据、30期30m分辨率的武汉大学土地利用数据、1期10m分辨率的清华大学数据、1期Esri的10m分辨率数据。另外，DEM数据目前已经分县扩展到了12.5m分辨率。

3.土壤数据、温度数据、降水数据等也在路上。

4.关于数据的获取，前往gzh，回复“水体”

2021年全国水体分布（按省、市、县）矢量数据的制作与分享相关推荐

全国各省10米分辨率的土地利用数据的制作与分享
1.简介该数据来源是清华大学的宫鹏老师团队做的全球土地利用数据,但该数据只能分影像块进行下载,没有制作按区划数据进行归纳,本人使用python爬虫下载了全球共计7000景影像,并利用2021年的区划 ...
数字-模拟信号混合传输收发机（E 题）--2021 年全国大学生电子设计竞赛
数字-模拟信号混合传输收发机(E 题)–2021 年全国大学生电子设计竞赛一任务设计并制作在同一信道进行数字-模拟信号混合传输的无线收发机.其中,数字信号由 4 个 0~9 的一组数字构成:模拟 ...
2021 年全国职业院校技能大赛网络搭建与应用赛项
2021 年全国职业院校技能大赛网络搭建与应用赛项正式赛卷 (总分 1000 分) ZZ-2021030 网络搭建与应用赛项执委会及专家组 2021 年 05 月 20 日 2021 年全国职业院 ...
2021年全国大学生数据统计与分析竞赛赛题B—基于机器学习的用户消费行为预测（上）
2021年全国大学生数据统计与分析竞赛赛题B-基于机器学习的用户消费行为预测(上) 一.赛题B:用户消费行为价值分析二.基于机器学习的用户消费行为预测 1.数据简介及清洗 (一)数据简介 (二)数据 ...
2021 年全国职业院校技能大赛网络搭建与应用赛项正式赛卷
2021 年全国职业院校技能大赛网络搭建与应用赛项正式赛卷技能要求 (总分1000 分) ZZ-2021030 网络搭建与应用赛项执委会及专家组2021 年 05 月 20 日竞赛说明一.竞赛 ...
2021年全国职业院校技能大赛：网络系统管理项目-真题-模块C-2全套视频讲解含无线地勘
2021年全国职业院校技能大赛高职组网络系统管理赛项模块C竞赛试题(C-2) 目录考试说明.................................................... ...
2021年全国在映影院数量、观影人次及平均票价走势分析[图]
影院规模可分为五类,从特大型到VIP型不等,而现今影院发展也已经经历四代,从老式电影院发展到了现在的商业综合体.在中国,影院需要加盟特定院线才能够放映. 2015-2021年我国影院分级城市分布无显著 ...
2021年全国大学生电子设计竞赛F题中数字识别这一技术分支实践与学习笔记
文章目录 Chapter 1 简介 Chapter 2 制作思路 Section 1 环境与成像效果及思路 Section 2 制作并训练数据集 Section 3 电脑环境配置 Section 4 ...
湖北省计算机四级报名考试,2021上半年全国大学英语四级考试(湖北考区)报名须知...
2021上半年全国大学英语四级考试(湖北考区)报名须知英语四级考试目的是推动大学英语教学大纲的贯彻执行,对大学生的英语能力进行客观.准确的'测量,为提高我国大学英语课程的教学质量服务.以下是小编为大 ...
2021 年全国硕士研究生入学统一考试管理类专业学位联考逻辑试题
2021 年全国硕士研究生入学统一考试管理类专业学位联考逻辑试题一. 逻辑推理:第 26~55 小题,每小题 2 分,共 60 分.下列每题给出的 A.B.C.D.E 五个选项中,只有一项是符合试题 ...

2021年全国水体分布（按省、市、县）矢量数据的制作与分享