# -*- 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("生成文件已添加至当前图层！")
Arcpy 提取质心、添加属性值、插值、批量出图脚本文件汇总 （244个子文件） 
扬州市.cpg 5B 苏州市.cpg 5B 淮安市.cpg 5B 南通市.cpg 5B 徐州市.cpg 5B 泰州市.cpg 5B JS_county_zx.cpg 5B 常州市.cpg 5B 镇江市.cpg 5B 宿迁市.cpg 5B 南京市.cpg 5B 连云港市.cpg 5B 无锡市.cpg 5B 盐城市.cpg 5B ZX2.cpg 5B ZX1.cpg 5B province.cpg 5B county.cpg 5B city.cpg 5B JS_county_zx.dbf 30KB ZX2.dbf 28KB ZX1.dbf 28KB county.dbf 23KB city.dbf 2KB 宿迁市.dbf 365B 镇江市.dbf 365B 徐州市.dbf 365B 连云港市.dbf 365B 苏州市.dbf 365B 南通市.dbf 365B 淮安市.dbf 365B 常州市.dbf 365B 扬州市.dbf 365B 盐城市.dbf 365B 南京市.dbf 365B 无锡市.dbf 365B 泰州市.dbf 365B province.dbf 241B .gitignore 50B Geo_Task.iml 357B
2.jpg 232KB 4.jpg 230KB 5.jpg 229KB 8.jpg 228KB 1.jpg 228KB 13.jpg 228KB 7.jpg 221KB 11.jpg 217KB 9.jpg 215KB 12.jpg 212KB 10.jpg 212KB 3.jpg 206KB 6.jpg 205KB city.shp.YUNJIA.18164.14540.sr.lock 0B Population_style.lyr 14KB Jiangsu_city_population.mxd 1.82MB 人口分布图_绘图模板.mxd 1.6MB
盐城市.png 61KB 淮安市.png 59KB 宿迁市.png 55KB 徐州市.png 55KB 扬州市.png 54KB 苏州市.png 54KB 南通市.png 54KB 连云港.png 54KB 泰州市.png 53KB 南京市.png 52KB 无锡市.png 52KB 常州市.png 51KB 镇江市.png 49KB 徐州市.prj 140B 镇江市.prj 140B 盐城市.prj 140B 南京市.prj 140B 扬州市.prj 140B 连云港市.prj 140B 无锡市.prj 140B 常州市.prj 140B 苏州市.prj 140B 南通市.prj 140B 泰州市.prj 140B 淮安市.prj 140B 宿迁市.prj 140B province.prj 140B city.prj 140B county.prj 140B Geo_Task_Interpolation.py 3KB Geo_Task_Centroid.py 2KB Geo_Task_Batch_Export.py 2KB Geo_Task_Add_Attribute.py 2KB Geo_Task_Batch_Clip.py 1KB Geo_Task_Batch_ExportToPNG.py 1KB county.sbn 1KB city.sbn 260B 扬州市.sbn 132B 宿迁市.sbn 132B 连云港市.sbn 132B 南京市.sbn 132B 徐州市.sbn 132B 盐城市.sbn 132B共 244 条
- 1
- 2
- 3
资源评论
