【ArcGIS】插值验证

1.抽取测试数据和验证数据

在被插值的数据中提取测试数据和验证数据
从地统计分析工具条（Geostatistical Analyst）中选择子集要素

选择验证数据和插值数据存储位置，推荐放在一个数据库中
填写抽取百分比数

以上数据准备完毕

2.插值

使用某种插值方法对上一步形成的插值数据进行插值，形成插值预测图，右键点击插值预测图，生成栅格文件

建议导出至上一步数据库中，像元大小可选为20

3.插值验证

（1）图斑赋值

将上一步生成的栅格数据，赋值到待插值图斑上

在工具箱中选择 3D Analyst Tools.tbx →功能性表面→添加表面信息

以上图斑赋值完毕

（2）验证点赋值

通过点位置挂接面，来给验证点赋值

打开测试点的属性表，在属性表中进行连接操作，选择基于位置关联


挂接成功后，验证点位显示插值数据

（2）插值赋值

对验证数据新建一个验证字段RE，类型为int，计算原始值和插值值相对误差，相对误差小于或等于0.3，赋1，否则赋0
最后统计1的数量占比，是否大于70%，若大于则插值有效，否则则无效。
以下为计算RE代码

def re(x,y):
    if((math.fabs(x-y)/x)<=0.3):
        return 1
    else:
       return 0