1:100万公众版基础地理信息数据(2021)—python数据预处理

1:100万公众版基础地理信息数据(2021)—python数据预处理

0 前言

关于1:100万公众版基础地理信息数据(2021)全网的介绍、下载和处理教程有很多,有些是用ArcMap建模处理,有些是用bash和python分步骤处理,本文只需要一个python代码+arcgis pro的融合dissolve即可把所有的GDB文件处理成合并好的shp。

Python代码


# Weiguo Xie YITC 20250101  gisxieweiguo@qq.com

import arcpy
import os

# 输入参数
gdb_folder = r"D:\BaiduNetdiskDownload\2021公众版国家基础信息系统1:100万行政区划数据\GDB"  # GDB文件夹路径
output_folder = r"D:\Output"  # 输出SHP文件夹路径
feature_list = ["CPTP", "CPTL", "HYDA", "HYDL", "HYDP", "HFCA", "HFCL", "HFCP", 
                "RESA", "RESL", "RESP", "RFCA", "RFCL", "RFCP", "LRRL", "LRDL", 
                "LFCL", "LFCP", "PIPL", "PIPP", "BOUA", "BOUL", "BOUP", "BRGA", 
                "BRGL", "BRGP", "TERA", "TERL", "TERP", "VEGA", "VEGL", "VEGP", 
                "AGNP", "AANP"]  # 需要合并的要素名称列表,有些GDB里面没有所有的类别,所以会显示某要素不存在某GDB中,这个是正常情况

# 创建输出文件夹
if not os.path.exists(output_folder):
    os.makedirs(output_folder)

# 创建一个字典,用于存储每个要素的临时合并结果
merged_features = {feature: [] for feature in feature_list}

# 遍历文件夹中的所有GDB
for gdb in os.listdir(gdb_folder):
    gdb_path = os.path.join(gdb_folder, gdb)
    if os.path.isdir(gdb_path) and gdb_path.endswith(".gdb"):
        print(f"正在处理GDB: {gdb_path}")
        arcpy.env.workspace = gdb_path

        # 遍历需要合并的要素
        for feature in feature_list:
            if arcpy.Exists(feature):
                print(f"  找到要素: {feature}")
                merged_features[feature].append(os.path.join(gdb_path, feature))
            else:
                print(f"  要素 {feature} 不存在于 {gdb_path}")

# 合并并导出为SHP
for feature, paths in merged_features.items():
    if paths:
        print(f"正在合并要素: {feature}")
        output_shp = os.path.join(output_folder, f"{feature}.shp")

        # 合并要素
        arcpy.management.Merge(paths, output_shp)
        print(f"  已导出合并结果: {output_shp}")
    else:
        print(f"要素 {feature} 在所有GDB中均不存在,跳过合并。")

print("所有要素处理完成!")

将代码放在Arcgis Pro的python模块中运行即可,整个数据的处理过程比较漫长,我的电脑内存24G,用了一个小时左右处理完。

在这里插入图片描述

合并完成后,有些线和面还需要根据相同的属性进行融合,在Arcgis Pro中打开融合工具,输入shp和融合字段
在这里插入图片描述

  • 融合前:
    在这里插入图片描述

  • 融合后:
    在这里插入图片描述


  • 成果数据组织
    全国1:100万公众版地形数据(2021)内容包括水系、居民地及设施、交通、管线、境界与政区、地貌与土质、植被、地名及注记9个数据集。

数据分层的命名采用四个字符,第一个字符代表数据分类,第二三个字符是数据内容的缩写,第四个字符代表几何类型。
在这里插入图片描述
在这里插入图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Vigo*GIS_RS

来瓶可乐~

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值