数据提取

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

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 = gdb
    feature_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”关键字的要素类，并将其保存到统一的地理数据库中。