arcpy数据分析自动化

数据提取

假设土地利用数据存储在多个地理数据库中，我们需要将这些数据提取到一个统一的地理数据库中。

import arcpy# 设置工作空间
arcpy.env.workspace = r"C:\GIS_Data\Land_Use"# 获取所有地理数据库中的土地利用数据
gdb_list = arcpy.ListWorkspaces("*", "FileGDB")
output_gdb = r"C:\GIS_Data\Land_Use\Combined_Land_Use.gdb"# 创建输出地理数据库
if not arcpy.Exists(output_gdb):arcpy.CreateFileGDB_management(r"C:\GIS_Data\Land_Use", "Combined_Land_Use.gdb")# 遍历每个地理数据库，提取土地利用数据
for gdb in gdb_list:arcpy.env.workspace = gdbfeature_classes = arcpy.ListFeatureClasses("*Land_Use*")for fc in feature_classes:output_fc = arcpy.FeatureClassToFeatureClass_conversion(fc, output_gdb, fc.split("_")[-1])print(f"数据 {fc} 已提取到 {output_fc}")

这段代码会从每个地理数据库中提取包含“Land_Use”关键字的要素类，并将其保存到统一的地理数据库中。