写在前面的话:

我知道很多同学都在等克里金,但是空间插值这个系列我已经预定好了撰写的思路,所以暂时没办法直接跳过前面的内容直奔克里金,所以大家只能将就一下先把前面的内容看完,但是应该快了。


前文再续,书接上一回。

上一章节我们讲了IDW插值法的计算,核心就在于根据距离把权重给算出来,然后带着权重用值预测出需要插值的位置的值是多少。而且也演示了ArcGIS如何使用这个工具来进行插值。

本章我们先来讲讲ArcGIS工具中IDW的一些参数设置。

首先,IDW插值,默认的情况是一种区域局部性的插值,无论是选择N个观测点的子集,还是选择固定距离的方式,都是利用了局部数据来进行预测,从理论上来看,这是符合地理空间相关性和异质性理论的。

但是有一种情况,就是认为,如果我的插值区域很小,例如:只有一个县,几十平方公里,而且观测站点非常少(比如少于12个),如果还是按照12个点来进行插值,就等于每个位置都利用了所有的观测点了——这种插值思路,实际上是一种全局插值的思想,有时候还挺有用的哦。

局部可变插值出来的结果,与采样点的分布有很大的关系,比如还是按照12个点来进行查找,那么观测点密集的地方,很近的区域就能找到12个点;而观测点稀疏地方,就需要很大的范围才能找到12个点,这样插值的情况,可能是这样的:

那么我们如果换成固定模式呢:

可以看见,如果采用固定模式,一般来说寻找到的子集大多不会发生剧烈的跳变(当然也有可能跳变,视观点分布的情况而定),那么如果说,我就是需要做到最平滑,忽略你的各种地理原理,最简单的方法,就是采用全局的模式,认为整个区域中,所有的观测点都参与到插值中:

在这种模式下,所有的点均参与插值计算,每个位置的计算,都不会因为出现子集不统一而产生的跳变。所以这种全局搜索模式下,可以生成最平滑的结果。

那么还有的同学发现,如果我要的是覆盖北京的全部的一个插值结果,比如你说的这个ArcGIS的“问题”:

前面我们已经说过,这是默认的方法是内插,也就是范围内的插值,但是现在我们要覆盖全图的结果,而不仅仅是内插得来,应该怎么去设置呢?

在ArcGIS里面,把内插改为外插,主要通过环境变量里面的处理范围这个选项来设定,如下:

我们修改处理范围,为我们需要覆盖的整个区域,比如设置为北京市的行政区划:

注意,我这里的观测点数据是投影坐标系——插值因为涉及到距离运算,所以最好不要采用地理坐标系,如果是地理坐标系的,最好先做投影。

运行插值,结果如下:

这样就生成了一个能够覆盖整个区域的插值栅格。

然后有同学又产生了一个问题,如果我仅仅需要北京市范围(严丝合缝的北京行政区划,多一点都不要,少一点都不行)咋办?

遇上这种情况,只需要在用一个工具就可以了。在ArcGIS的toolbox里面,有一个工具,叫做按掩膜提取,位置在:Spatial Analyst Tools —— Extraction —— Extract by Mask:

点击打开工具,设置如下:

然后运行工具,就可以得到你需要的区域了:

当然,采用数据管理工具箱里面的栅格裁剪也能实现,大家自己去测试好了,我个人习惯喜欢用这个提取工具,方便省事参数少。

到此,IDW插值,打完收工。

下期预告:插值中仅次于IDW第二简单的插值算法:径向基函数插值(RBF),有时候也叫样条函数插值。它可以生成最平滑的曲线,是强迫症患者的福音哦。

待续未完。

白话空间统计二十九:空间插值(六)IDW部分完结篇相关推荐

  1. 白话空间统计二十九:空间插值(二)

    前文再续,书接上一回--上回书说到,空间插值可以对数据进行估算,这一张我们来具体说说空间插值的一些概念. 首先,来说说插值的输入与输出问题. 首先输入的肯定是用来进行插值的观测点数据,一般来说都是点要 ...

  2. 白话空间统计二十四:地理加权回归(八)结果解读(一)

    地理加权回归分析完成之后,与OLS不同的是会默认生成一张可视化图,像下面这张一样的: 这种图里面数值和颜色,主要是系数的标准误差.主要用来衡量每个系数估计值的可靠性.标准误差与实际系数值相比较小时,这 ...

  3. 白话空间统计二十四:地理加权回归(九)结果解读(二)

    实际上,除了辅助表以外,GWR还会生成一份全要素的表.对回归的每一个样本都给出相应的信息,今天就来看看这些信息代表了什么内容. 生成的新的要素类字段信息如下: 实际上,Coeffcient(系数)和S ...

  4. 白话空间统计二十一:密度分析(一)

    白话空间统计二十一:密度分析(一) 密度分析这个概念其实很早就想写了,也有无数同学都问过我,虾神你能不能讲讲那些漂亮的热度图是怎么做的啊?比如下面这种: 如果说,这是互联网地图里面,最让人喜闻乐见的一 ...

  5. 白话空间统计二十一:密度分析(五)带宽与核表面曲率的关系

    白话空间统计二十一:密度分析(五) 上次讲密度分析的时候,有同学问道带宽的问题,实际上我翻 了一下以前写的文章,在密度分析一.二里面,都对这个有过描述,详细的可以回去翻一下(很老的文章了--可能要翻到 ...

  6. 白话空间统计二十三:回归分析(四)

    白话空间统计二十三:回归分析(四) 今天把回归最后几个概念介绍一下,后面就不再说了--开始写回归的时候,我发现我掉入了一个巨大的坑里面.就说回归分析这种东东,汗牛充栋啊,在任何大学里面,妥妥一个学期的 ...

  7. 2021年大数据Hadoop(二十九):​​​​​​​关于YARN常用参数设置

    全网最详细的Hadoop文章系列,强烈建议收藏加关注! 后面更新文章都会列出历史文章目录,帮助大家回顾知识重点. 目录 本系列历史文章 前言 关于yarn常用参数设置 设置container分配最小内 ...

  8. android游戏性能测试,评测 篇二十九:用游戏数据来彰显实力,华为nova5i Pro游戏性能篇...

    评测 篇二十九:用游戏数据来彰显实力,华为nova5i Pro游戏性能篇 2020-01-02 00:20:00 4点赞 1收藏 2评论 笔者带着大家对华为nova5i Pro做了开箱测试,了解外观与 ...

  9. CSDN 编程竞赛二十九期题解

    竞赛总览 CSDN 编程竞赛二十九期:比赛详情 (csdn.net) 竞赛题解 题目1.订班服 小A班级订班服了!可是小A是个小糊涂鬼,整错了好多人的衣服的大小.小A只能自己掏钱包来补钱了.小A想知道 ...

最新文章

  1. pip、NumPy、Matplotlib在Windows上的安装过程
  2. COM组件转换成.NET组件[转]
  3. EOS共识机制——DPoS代理权益证明
  4. C++ explicit的作用
  5. 如何通过parent object ID获得underlying object structure IBASE ID
  6. java ftp 中文上传_java实现ftp文件上传下载,解决慢,中文乱码,多个文件下载等问题...
  7. 互联网级监控系统必备-时序数据库之Influxdb集群及踩过的坑
  8. Incorrect string value: '/xE7/xA8/x8B/xE5/xBA/x8F...' for column 'course' at row 1
  9. dedecms m view.php,dedecms在php7下的使用方法 织梦dedecsm后台一片空白的解决方法
  10. java-第十章-类和对象-创建管理员对象
  11. [原] 求薪水第三高的雇员的SQL题(图解)
  12. 设计模式---桥接模式(C++实现)
  13. Springboot 统一maper 和逆向工程
  14. Android 友盟分享(截图指定的View分享)
  15. 宇宙最强API接口调试工具Apipost
  16. AHP(层次分析法)学习笔记及多层权重Python实践
  17. 增量式编码器工作原理以及使用
  18. 奇偶校验c语言ascii,奇偶校验(parity check)
  19. 矩阵的Kronecker积的相关结论
  20. 华硕打开桌面计算机没有反应,华硕window10电脑打开就这样,进不去桌面,怎么处理...

热门文章

  1. 北大青鸟消防设备类型编码_北大青鸟JBF4372E编码型火灾声光警报器
  2. LeeCode 1579 并查集
  3. Windows上Android手机屏幕共享指南
  4. URL转二维码,图片
  5. STM32F103学习笔记(1)掌握GPIO操作—输出模式(推挽/开漏)
  6. 谈信息化监理实施的方法
  7. 什么是延时电路?6种延时电路原理讲解
  8. 三种方式实现双飞翼布局
  9. java基于springboot+vue的社区居民健康管理系统 nodejs 前后端分离
  10. [BZOJ4455][UOJ185][Zjoi2016]小星星(树形DP+容斥)