ArcGIS Pro for Python批量处理实战10大高频脚本+逐行注释GISer直接抄作业.

# 导入ArcGIS核心库import arcpy# 导入路径处理库import os# 允许覆盖输出文件arcpy.env.overwriteOutput = True# 设置输入栅格所在文件夹（改成你自己的栅格目录）arcpy.env.workspace = r"D:\data\input_rasters"# 裁剪边界文件（shp或gdb要素类）clip_mask = r"D:\data\mask.shp"# 裁剪结果输出文件夹out_folder = r"D:\data\clip_result"# 自动列出工作空间下所有栅格文件（tif、img等）rasters = arcpy.ListRasters()# 开始循环，逐个处理每一张栅格for ras in rasters:    # 拼接输出路径：输出文件夹 + 原文件名，保证名字不变    out_ras = os.path.join(out_folder, ras)# 执行栅格裁剪工具arcpy.management.Clip(in_raster=ras,rectangle="",out_raster=out_ras,in_template_dataset=clip_mask,nodata_value="",clipping_geometry="ClippingGeometry")# 输入当前栅格# 范围留空，使用边界范围# 输出栅格路径# 裁剪边界# 空值不设置 # 按边界形状裁剪    # 控制台打印进度，方便查看运行到哪一个    print(f"已完成裁剪：{ras}")

import arcpyimport osarcpy.env.overwriteOutput = True# 设置输入矢量所在的地理数据库（gdb路径）arcpy.env.workspace = r"D:\GIS项目\原始矢量数据.gdb"# 裁剪边界路径clip_mask = r"D:\GIS项目\项目边界.shp"# 输出成果地理数据库（提前新建空gdb）out_gdb = r"D:\GIS项目\裁剪后矢量成果.gdb"# 自动遍历工作空间内所有矢量要素类fcs = arcpy.ListFeatureClasses()# 循环批量裁剪for fc in fcs:    # 输出要素命名：原名称+_clip，方便区分    out_fc = os.path.join(out_gdb, fc + "_clip")    # 执行矢量裁剪工具    arcpy.analysis.Clip(fc, clip_mask, out_fc)    # 打印处理日志    print(f"矢量裁剪完成：{fc}")

import arcpyimport osarcpy.env.overwriteOutput = True# ====================== 【这里改成你自己路径】 ======================arcpy.env.workspace = r"E:\测试数据\待投影数据.gdb"  # 原始GDBout_gdb = r"E:\测试数据\CGCS2000投影成果.gdb"  # 输出GDB# ！！重要：如果你的数据原本是某已知坐标系，在这里指定原始WKIDsource_wkid = 4537  # 举例：比如原本是西安80 3度带，改成你实际的。# =================================================================# 创建输出GDB（不存在就新建）if not arcpy.Exists(out_gdb):    arcpy.management.CreateFileGDB(os.path.dirname(out_gdb), os.path.basename(out_gdb))# 目标坐标系 CGCS2000 地理坐标系(4490)out_sr = arcpy.SpatialReference(4490)# 遍历所有要素类fcs = arcpy.ListFeatureClasses()for fc in fcs:    # 获取当前要素空间参考    desc = arcpy.Describe(fc)    in_sr = desc.spatialReference    # 处理：无坐标系/未知坐标系    if in_sr.name == "Unknown":        print(f"【警告】{fc} 无定义坐标系，临时指定预设原始坐标系WKID:{source_wkid}")        in_sr = arcpy.SpatialReference(source_wkid)        # 先给图层定义投影        arcpy.management.DefineProjection(fc, in_sr)    out_fc = os.path.join(out_gdb, f"{fc}_CGCS2000")    try:        arcpy.management.Project(fc, out_fc, out_sr)        print(f"✅ 投影成功：{fc} -> {os.path.basename(out_fc)}")    except Exception as e:        print(f"❌ 投影失败：{fc}，错误信息：{str(e)}")

import arcpyimport osarcpy.env.overwriteOutput = True# 输入需要处理的矢量数据路径arcpy.env.workspace = r"E:\测试数据\测试数据1.gdb"# 遍历所有矢量要素fcs = arcpy.ListFeatureClasses()for fc in fcs:    # 1. 添加面积字段，字段类型为双精度（DOUBLE），存储小数    arcpy.management.AddField(fc, "mj", "DOUBLE", field_alias="实测面积(㎡)")    # 2. 计算面积，单位：平方米，如需公顷改成"SQUARE_HECTARES"    arcpy.management.CalculateGeometryAttributes(fc, [["mj", "AREA"]], area_unit="SQUARE_METERS")    # 3. 添加x字段    arcpy.management.AddField(fc, "CENTROID_X", "DOUBLE", field_alias="X坐标")    # 4. 添加y字段    arcpy.management.AddField(fc, "CENTROID_Y", "DOUBLE", field_alias="y坐标")    # 4. 游标更新坐标    with arcpy.da.UpdateCursor(fc,["SHAPE@","CENTROID_X","CENTROID_Y"]) as cursor:        for r in cursor:            c = r[0].centroid            r[1] = c.X            r[2] = c.Y            cursor.updateRow(r)    print(f"字段添加+几何计算完成：{fc}")

#添加长度字段arcpy.management.AddField(in_table=fc,field_name="Shape_Length",field_type="DOUBLE",field_alias="长度(m)")#计算长度arcpy.management.CalculateGeometryAttributes(in_features=fc,geometry_property=[["Shape_Length","LENGTH"]],length_unit="METERS")print(f"字段与几何计算完成：{fc}")

import arcpyarcpy.env.overwriteOutput = True# 输入需要分割的完整矢量数据路径input_fc = r"E:\测试数据\测试数据1.gdb\test"# 分割后成果输出文件夹（空文件夹）out_workspace = r"E:\测试数据\按乡镇分割成果"# 执行按属性分割工具，最后一个参数是分割字段（改成你的字段名：乡镇/区县/图幅号）arcpy.analysis.SplitByAttributes(input_fc, out_workspace, "CITY_NAME")

import arcpyarcpy.env.overwriteOutput = True# 输入待合并的矢量数据所在gdbarcpy.env.workspace = r"E:\测试数据\处理结果.gdb"# 筛选同类要素：比如只合并前缀为road_的路网数据，按需修改匹配规则fcs = arcpy.ListFeatureClasses("test_*")# 合并后成果输出路径out_fc = r"E:\测试数据\合并成果.gdb\megre"# 执行合并工具，自动拼接所有符合条件的要素arcpy.management.Merge(fcs, out_fc)print(f"批量合并完成，共合并{len(fcs)}个图层")

import arcpyimport osarcpy.env.overwriteOutput = True# 输入原始栅格数据文件夹arcpy.env.workspace = r"E:\测试数据\原始栅格dem"# 输出重采样后栅格文件夹out_folder = r"E:\测试数据\成果栅格dom"# 遍历所有栅格rasters = arcpy.ListRasters()for ras in rasters:    out_ras = os.path.join(out_folder, ras)    # 重采样工具：输入、输出、分辨率（1.0=1米）、重采样方法（双线性插值，适合栅格）    arcpy.management.Resample(ras, out_ras, 1.0, "BILINEAR")    print(f"栅格重采样完成：{ras}")

import arcpyimport osarcpy.env.overwriteOutput=True#输入gdb数据路径arcpy.env.workspace=r"D:\GIS项目\成果数据.gdb"#导出SHP输出文件夹out_folder=r"D:\GIS项目\SHP格式成果"#遍历所有矢量，批量导出SHPfcs=arcpy.ListFeatureClasses()for fc in fcs:    arcpy.conversion.FeatureClassToShapefile(fc,out_folder)    print(f"批量导出SHP完成：{fc}")

importarcpyimportosarcpy.env.overwriteOutput=True#输入gdb数据路径arcpy.env.workspace=r"D:\GIS项目\成果数据.gdb"#导出SHP输出文件夹out_folder=r"D:\GIS项目\CAD格式成果"#遍历所有矢量，批量导出CAD#【CAD导出替换代码】，把上面导出SHP的循环换成这段即可for fc in fcs:    out_cad=os.path.join(out_folder,fc+".dwg")    arcpy.conversion.ExportCAD(fc,"DWG_R2018",out_cad)    print(f"批量导出CAD完成：{fc}")

import arcpyarcpy.env.overwriteOutput = True##输入待质检数据路径arcpy.env.workspace = r"E:\测试数据\测试数据1.gdb"# 遍历数据库下所有图层fcs = arcpy.ListFeatureClasses()for fc in fcs:    # 检查几何错误    arcpy.management.CheckGeometry(fc)    # 自动修复几何错误    arcpy.management.RepairGeometry(fc)    # 选择没有几何图形的空要素    arcpy.management.SelectLayerByAttribute(in_layer_or_view=fc,selection_type="NEW_SELECTION",where_clause="Shape IS NULL")    # 删除选中的空要素    arcpy.management.DeleteRows(fc)    print(f"数据质检与修复完成：{fc}")

import arcpyimport osarcpy.env.overwriteOutput=True#输入待拼接栅格文件夹arcpy.env.workspace=r"E:\测试数据\原始栅格路径"#拼接后成果输出路径out_mosaic=r"E:\测试数据\成果栅格dom\成果DOM.tif"#遍历所有栅格rasters=arcpy.ListRasters()#执行栅格拼接工具arcpy.management.MosaicToNewRaster(rasters,os.path.dirname(out_mosaic),os.path.basename(out_mosaic),pixel_type="8_BIT_UNSIGNED",number_of_bands=3)print("栅格批量拼接完成！")

import arcpyaoi_fc = r"E:\测试数据\原始数据\乡界.shp"        # 乡镇边界out_pdf_folder = r"E:\测试数据\成果数据"         # 导出PDF目录name_field = "XZQMC"             # 名称字段# 获取当前打开的ArcGIS工程aprx = arcpy.mp.ArcGISProject("CURRENT")# 获取第一个布局（你自己做好的制图模板）layout = aprx.listLayouts()[0]# 获取布局中的地图框元素map_frame = layout.listElements("MAPFRAME_ELEMENT")[0]# 遍历每个乡镇，自动定位并导出with arcpy.da.SearchCursor(aoi_fc, ["SHAPE@", name_field]) as cursor:    for row in cursor:        # 获取当前乡镇范围        ext = row[0].extent        # 地图框缩放到该范围        map_frame.camera.setExtent(ext)        # 拼接PDF路径        out_pdf = os.path.join(out_pdf_folder, f"{row[1]}.pdf")        # 导出高清PDF        layout.exportToPDF(out_pdf, resolution=300)        print(f"已导出地图：{row[1]}")