Arcpy 提取质心、添加属性值、插值、批量出图脚本文件汇总

# -*- coding: GB2312 -*- # (4)编写脚本工具，对江苏2021年常住人口数据进行插值，插值方法任选； import arcpy # 定义添加文件到当前图层的函数add_currentMxd() def add_currentMxd(fileDir): mxd = arcpy.mapping.MapDocument("CURRENT") df = arcpy.mapping.ListDataFrames(mxd)[0] addLayer = arcpy.mapping.Layer(fileDir) arcpy.mapping.AddLayer(df, addLayer, "AUTO_ARRANGE") # 图层刷新 arcpy.RefreshActiveView() arcpy.RefreshTOC() del mxd arcpy.env.overwriteOutput = True # 文件路径设置 pointFile = arcpy.GetParameterAsText(0) # 插值点文件 outFolder = arcpy.GetParameterAsText(1) # 输出文件夹 outFile = outFolder + "/JS_County_Point3_" # 环境设置 maskFile = arcpy.GetParameterAsText(2) # 掩膜文件 arcpy.env.workspace = "D:/Geo_Task/WorkSpace" # 设置工作空间 arcpy.env.mask = maskFile # 设置掩膜文件 arcpy.env.extent = maskFile # 设置处理范围 arcpy.CheckOutExtension("Spatial") IDW_Flag = arcpy.GetParameterAsText(3) Kriging_Flag = arcpy.GetParameterAsText(4) Spline_Flag = arcpy.GetParameterAsText(5) # 插值参数设置 zField = "Num" # 插值字段：Z值字段 # 输出文件名称设置 IDW_outFileName = "IDW_1.tif" Kriging_outFileName = "Kriging_1.tif" Spline_outFileName = "Spline_1.tif" # 输出像元大小设置 cellSize = "1.74676370500002E-02" # IDW插值-反距离权重法 if IDW_Flag == 'true': arcpy.AddWarning("IDW插值中……") arcpy.gp.Idw_sa(pointFile, # 输入点要素 zField, # Z值字段 # 输出栅格路径 outFile + IDW_outFileName, cellSize, # 输出像元大小 "2", "VARIABLE 12", # 搜索半径点数 "") arcpy.AddMessage("IDW插值完成！") add_currentMxd(outFile + IDW_outFileName) # Kriging插值-克里金法 if Kriging_Flag == 'true': arcpy.AddWarning("Kriging插值中……") arcpy.gp.Kriging_sa(pointFile, zField, outFile + Kriging_outFileName, # 半变异模型-球面函数，步长大小 "Spherical 0.015495", cellSize, "VARIABLE 12", # 搜索半径点数 "") arcpy.AddMessage("Kriging插值完成！") add_currentMxd(outFile + Kriging_outFileName) # Spline插值-样条函数法 if Spline_Flag == 'true': arcpy.AddWarning("Spline插值中……") arcpy.gp.Spline_sa(pointFile, zField, outFile + Spline_outFileName, cellSize, "TENSION", # 样条函数类型 "0.1", # 权重 "12" # 点数 ) arcpy.AddMessage("Spline插值完成！") add_currentMxd(outFile + Spline_outFileName) arcpy.AddMessage("插值结束！") arcpy.AddWarning("生成文件已添加至当前图层！")