【数据分享】2021年全国水体分布(按省、市、县)矢量数据
文章来源:锐多宝的地理空间
偶然间获得了Open Street Map的2021年全国水体分布Shp矢量数据,大约400M左右。
当我加载到arcgis中,感觉到了明显的卡顿。由于之前写过分省、市、县裁剪矢量的代码,顺手就将水体数据按照2021年的行政边界裁剪了。
该数据分为3个层次结构,分别为省、地级市、县。其中,省级水体31个(不包含HongKong、TaiWan等), 地级市水体366个,县级市水体2847个。每一个文件夹中都包含该省、地级市或者县的水体矢量数据、行政边界以及数据说明。
2.1、制作方法流程图
2.2、Python裁剪脚本
主要是用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年全省水体分布_矢量")
#保存文件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 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+xian
shutil.copy(gson_path,xian_path)
except:
pass
通过上述步骤,我们获得了全国各省、地级市、县等行政单位的水体数据。首先展示省级水体数据:
广东省水体分布图
如果你要找县级市的水体分布图,进入"分县"文件夹中,选择对应省份-地级市-县进行选择。我们以四川省资阳市乐至县的水体数据为例, 直接选择分县-四川省-资阳市-乐至县,并将其加到Google影像图中,查看水体矢量效果。
本次资料需参与转发活动领取,Ai尚研修会员无需转发可直接领取
长按下方二维码关注"Ai尚研修"公众号
回复关键词“154” 查看领取方式!
1
教学方式
INTRODUCE
线上学习+导师群辅助+不定期线上答疑
2
主讲专家
INTRODUCE
刘瑞民(副教授),北京师范大学副教授、博士生导师。主要关注流域综合管理,数字流域建模,流域水土过程模拟,洪水淹没风险评估,流域面源模拟与控制,流域水资源优化配置等,为流域管理决策提供重要支撑。在Environmental Science &Technology、Water Research、Journal of Hazardous Materials,Journal of Cleaner Production等杂志发表论文100余篇。
3
教学特色
INTRODUCE
1、原理深入浅出的讲解;
2、技巧方法讲解,提供所有案例数据及代码;
3、与项目案例相结合讲解实现方法,对接实际工作应用 ;
4、跟学上机操作、独立完成案例操作练习、全程问题跟踪解析;
5、建立导师助学交流群,长期进行答疑及经验分享,辅助学习及应用,不定期召开线上答疑;
4
发票及报销文件
INTRODUCE
发票开具:培训费、会议费、资料费、技术咨询费等,配有盖章文件等,用于参会人员报销使用。
5
ArcGIS洪水灾害普查、风险评估及淹没制图
INTRODUCE
直播时间:2021年12月11日—13日
上下滑动查看更多
6
ArcGIS 在水土流失模拟及分析中的实践技术应用
INTRODUCE
直播时间:2021年12月24日—26日
上下滑动查看更多
7
SWAT模型在水文水资源、面源污染模拟
INTRODUCE
直播时间:2022年1月14日—1月17日
教学方式:线上学习+导师群辅助+不定期线上答疑
上下滑动查看更多
1
FVCOM三维水动力与水交换、污染物扩散数值模拟
INTRODUCE
直播时间:2021年12月18日—19日、25日
教学方式:线上学习+导师群辅助+不定期线上答疑
上下滑动查看更多
2
HYPE分布式水文模型建模
INTRODUCE
直播时间:2021年12月24日—12月26日
教学方式:线上学习+导师群辅助+不定期线上答疑
上下滑动查看更多
3
Python在气象与海洋中的实践技术应用
INTRODUCE
直播时间:2022年1月8日-9日 15日-16日
教学方式:线上学习+导师群辅助+不定期线上答疑
上下滑动查看更多
4
FVCOM溢油及物质输运数值模拟
INTRODUCE
直播时间:2022年1月14日—1月16日 (共计3天学习)
教学方式:线上学习+导师群辅助+不定期线上答疑
上下滑动查看更多
5
Delft3D建模、水动力模拟方法及在地表水环境
INTRODUCE
直播时间:2022年1月14日-16日(三全天)
教学方式:线上学习+导师群辅助+不定期线上答疑
上下滑动查看更多
6
联系我们
INTRODUCE
视频课程最高享受7.5折优惠
贾莲:185-1037-1192(同微)
免费数据资源
Ai尚研修海量资源(数据、课件、书籍、视频教程)一键免费领取
(点击数据标题,直接进入)
土利用数据
上下滑动查看更多
超级会员专享
END
声明: 本号旨在传播、传递、交流,对相关文章内容观点保持中立态度。涉及内容如有侵权或其他问题,请与本号联系,第一时间做出撤回。
END
Ai尚研修丨专注科研领域
技术推广,人才招聘推荐,科研活动服务
科研技术云导师,Easy Scientific Research