反距离权重法插值
即IDW(Inverse Distance Weight) ,也可以称为距离倒数乘方法。计算一个格网结点时给予一个特定数据点的权值与指定方次的从结点到观测点的该结点被赋予距离倒数成比例。
1、选取气象站点的点矢量数据进行插值
2、在ArcGIS中“Spatial Analyst工具”——插值分析——反距离权重法填写以下信息,z值字段为需要进行计算的字段,幂默认为2,搜索半径根据点矢量进行填写,注意最大距离需要进行转换。(假设在ArcGIS中测量距离最远的两个点的距离为1689千米,1689/111=15.2,则插值半径最大距离可设置为16)
3、结果如下
Arcpy批处理方法
# 批处理遍历所有点矢量文件
for pointshp_filename in pointshp_path:
if pointshp_filename.split(".")[1]=="shp":
print pointshp_filename
real_pointshp_filename = " "+pointshp_filename
# 暂存路径
temp_raster = " "+pointshp_filename.split(".")[0]+".tif"
# 输出路径
outraster = " "+pointshp_filename.split(".")[0]+".tif"
# 调用工具集
arcpy.CheckOutExtension("Spatial")
# 反距离权重插值(点矢量文件,字段名,插值像元大小,反距离幂值默认取2,矢量点个数和搜索最大半径
outIDW = Idw(real_pointshp_filename, "V12001_701", 0.00026949459, 2, RadiusVariable(39, 16))
outIDW.save(temp_raster)
# 掩膜矢量范围
outExtractByMask = ExtractByMask(temp_raster, " ")
outExtractByMask.save(outraster)
完整脚本可直接使用
反距离权重法插值批处理——IDW.py