在知识传播途中,向涉及到的相关著作权人谨致谢意!
南师国家精品课程《地理信息系统》——主讲人:韦玉春老师

文章目录

  • 1 空间插值
    • 1.1 插值的方法
      • 1.1.1 趋势面分析
      • 1.1.2 IDW插值
      • 1.1.3 [克里格插值](https://blog.csdn.net/summer_dew/article/details/83072804)
      • 1.1.4 几何方法
    • 1.2 空间插值的工作流程
    • 1.3 思考的问题
    • 1.4 原则

1 空间插值

【空间自相关】距离越近越相似
【什么叫变异】距离越近,变异越小

【插值的道理】近朱者赤,近墨者黑
【解释】某种属性的空间分布是受制于某个过程,了解了过程,才可以找到好的插值方法。过程控制格局,而数据是格局的一个表现。所以我们要找的是那个特定的条件是什么

【插值时为什么要考虑方向?】

  1. 不考虑方向,那就是说这个数据在四周是均匀分布的,是各向同性的
  2. 但实际不是这样的,一个东北风往这吹,那东北方向的东西和你就不是一致关系–>那插值时就应该按东北的方向来做

【如何考虑插值的方向?】是考虑了一个空间作用的关系,和研究的问题有关

  1. 若研究的对象是空间同性的,不考虑方向
  2. 若研究的对象是空间异性的,而且是沿着某个方向变化的,插值的时候要考虑方向,要有个约束条件

1.1 插值的方法

分类的角度 方法
数据分布规律 1. 基于规则分布数据的内插方法
2. 基于不规则分布的内插方法
3. 适合于等高线数据的内插方法等
内插函数与参考点的关系方面 1. 曲面通过所有采样点的纯二维插值方法
2. 曲面不通过参考点的曲面拟合插值方法
内插曲面的数学性质 1. 多项式内插
2. 样条内插
3. 最小二乘配置内插等内插函数;
对曲面理解的角度 1. 克立金法;2. 多层曲面叠加法;3. 加权平均法;4. 分形内插
内插点的分布范围 1. 整体内插;2. 局部内插和;3. 点内插法

1.1.1 趋势面分析

【趋势面】反映数据在空间上分布规律的曲面
【趋势面分析】空间的变化规律->找趋势->原因分析
【注意】不注重精度和检验
【解释】

  1. 找的一个大的规律、趋势,看这个规律是存在还是不存在的?
  2. 需不需要用其他的数据去检测?不需要。
    这个方法本身是不可以去检验的。并不是说高度拟合,它的精度就高;不拟合精度就低。只是某种趋势,某种启发

【假设】观测面=区域趋势+局部异常+随机干扰
【进一步假设】随机干扰不考虑

  • 观测面=趋势面+局部异常:判断水流、养分–>对趋势面感兴趣,对异常感兴趣
  • 局部异常=观测面-趋势面:找矿->对趋势面不感兴趣,对异常感兴趣

【趋势】变化比较缓慢、影响遍及整个研究区的区域成分组成
【局部异常】变化比较快,其影响在区内并非处处可见的成分

【怎么获取趋势面?】做一个回归,回归分析
【多项式回归分析】趋势面分析的实质是进行数据的拟合,对因变量无特别的要求,自变量是地理坐标

公式:Y(x,y)- Y:属性- x,y:坐标- Y=c0+a1*x1+b1*y1+c1x1*y1+… 用多项式来拟合,看看Y和x,y有什么关系
说明
数据 ①离散的点;②最好是均匀分布:各个方位的点数基本一致;分布范围规则
步骤
结果表示 等值线图;等值面图;3D图
趋势面分析的要点 离散数据的获取;多项式拟合的阶数选择
【解释】取几阶多项式,是对问题的理解有关的
1. 先做一阶,看能不能解释这个结果
2. 如果一阶不能解释,再做二阶,看看能不能解释这个结果,能解释就没有必要做下去了
3. 并不是说某一阶的结果才是正确的,要看对问题的理解是否有关,能不能解释这种现象

1.1.2 IDW插值

【IDW】距离倒数加权插值
【解释】距离越近影响越大,距离越远影响越小

1.1.3 克里格插值

【需要知道的东西】

  1. 线性无偏是BLUE的(最好的方法)
  2. 怎么才能BLUE呢?要达到它的假设
  3. 它的假设是什么呢?空间平稳性假设
  4. 插值之后有一个东西,方差变异函数

1.1.4 几何方法

  1. 整体内插法
  2. 分块内插法
  3. 逐点内插法

1.2 空间插值的工作流程

【工作流程】

  1. 插值方法(模型)的选择
  2. 空间数据的探索分析,包括对数据的均值、方差、协方差、独立性和变异函数的估计等
  3. 插值方法评价
  4. 重新选择内插方法,直到合理
  5. 空间内插

【要点】

  1. 内插点邻域范围的确定
  2. 权值确定方法-自相关程度
  3. 内插函数的选择

1.3 思考的问题

  1. 任何的插值方法都建立在一定的假设基础上,这些假设都有一定的局限性
    在这些假设中,我们往往回避的是空间不连续性和空间异质两个问题
    【解释】为什么能插值呢?我们相信它在空间上是连续的
    【举例】温度受一些大的空间现象影响(寒流),认为它在空间上是连续的
  2. 如何将地理空间分解为不同区域?分区后是否存在界面处的不连续现象?
    【解释】受边界效应的影响,分区域的时候应该考虑哪些边界?
  3. 能够进行智能化、通用化的插值吗?

1.4 原则

  1. 对于众多的空间内插方法而言,没有绝对最优的空间内插方法,只有特定条件下的最优方法
  2. 必须依据数据的内在特征,依据对数 据的空间探索分析,经过反复实验,选择最优的空间内插方法。同时,应对内插结果做严格的检验

[GIS原理] 10.2 空间插值相关推荐

  1. [GIS原理] 10 空间统计分析

    在知识传播途中,向涉及到的相关著作权人谨致谢意! 南师国家精品课程<地理信息系统>--主讲人:韦玉春老师 文章目录 背景 空间统计分析 1 研究目的 2 研究内容 3 空间统计分析的基本流 ...

  2. [GIS原理] 8.1 空间分析-网络分析

    在知识传播途中,向涉及到的相关著作权人谨致谢意! 文章目录 1 网络分析(Network Analysis) 1.1 矢量数据的网络分析 1.1.1 网络组成和属性 1.1.2 几何网络的建立 1.1 ...

  3. python2 使用arcpy自定义GIS工具箱,实现excel读取、空间插值

    一.实验目的 1.掌握Arcpy基本语法.功能函数 2.学会使用Arcpy自定义GIS工具箱 3.使用arcpy进行空间插值 4.自定义创建的Arcmap制图模板,使用arcpy调节图例.注记.比例尺 ...

  4. gstat | 空间插值(二)——克里金插值之普通克里金

    说明:昨天的推文误把可吸入颗粒物当作PM2.5,实应该为PM10,这里修正后重发. 从本篇开始计划分三篇介绍克里金插值.与反距离权重插值不同,克里金插值是无偏估计,其中也涉及到模型估计.本篇先对普通克 ...

  5. gstat | 空间插值(一)——反距离权重插值;使用ggplot2绘制地图

    本篇既是空间插值系列的第一篇推文,也是ggplot2工具包系列推文中的一篇.空间插值使用的工具包是gstat,该工具包主要用于地统计分析. library(gstat) 示例数据来自HSAR工具包: ...

  6. 关于GIS原理的实际分析应用题的一些解法

    话不多说,看题. 01 公园选址问题1 题目 请写出利用GIS技术进行公园选址的空间操作步骤.其中公园选址条件: 1)为了安静舒适,要求该园区离主要公路1公里以外,且交通方便,离主要公路3公里以内. ...

  7. 《Web GIS原理与应用开发》读书笔记(5)

    第五章 切片地图 以下均掺杂个人理解 书目:Web GIS 原理与应用开发 -- 刘光 曾敬文 曾庆丰 著 如果要提高Web地图的访问速度,使用地图切片的非常有效的方法. 地图切片就是在多个比例尺下配 ...

  8. 《Web GIS原理与应用开发》读书笔记(1)

    第一章 Web GIS 概述 以下均掺杂个人理解 书目:Web GIS 原理与应用开发 -- 刘光 曾敬文 曾庆丰 著 1.1 GIS的发展 20世纪60年代,世界上第一个地信系统:加拿大地理信息系统 ...

  9. 完美解决ArcGIS10.x栅格空间插值报错无法进行和导出插值栅格结果出错的问题

    不少GISer,尤其是初学者,对栅格数据的操作,主要是栅格插值无法进行,或者是插值后的栅格结果无法导出,或者导出结果为空,或者是栅格格式无法转换,今天刘一哥GIS就帮您排忧解难吧. 其实,关于栅格(尤 ...

  10. gstat | 空间插值(三)——克里金插值之泛克里金和简单克里金

    本篇接着上篇继续介绍克里金插值.首先加载相关工具包和上篇使用的示例数据: library(gstat) library(sf) library(tidyverse) library(readxl) l ...

最新文章

  1. IBM Thinkpad T43-44U 升级到 2G 内存后少了 66M
  2. linux ftp服务
  3. UIScrollView用法
  4. 你是“细”精你就赢了 游戏开发者怎样做好玩家细分
  5. 小米路由器r3gv2/r4a Lean的OpenWrt固件
  6. 2022华中杯C题矿井提升机钢丝绳的缺陷分析思路讲解
  7. MapABC地图API
  8. 周杰伦要出新专辑了?上 Instagram 看看
  9. 开源 iOS 项目分类索引大全
  10. 【TF2】Eager Execution机制
  11. 如何定位crach问题
  12. windows 凭据管理
  13. 弥补自己的无知: GPL与LGPL
  14. PS2019零基础教程分享给大家
  15. 对于家政020 APP平台如何走出资本寒冬?
  16. 【暴力枚举】速算游戏 fun.pas/c/cpp
  17. 载药脂质体并表面修饰各种分子定制-带有多肽的脂质分子
  18. PySide6中QlineEdit设置自动所有字母大写以及首字母大写的案例
  19. 简述UIView的属性和用法
  20. AirPods国行和海外版有什么区别

热门文章

  1. 【毕业求职季】-听说你想去大厂看学妹,带你看看字节广告运营岗面试长啥样?
  2. mysql 批量录入 id_Mysql 批量插入数据 提前获得主键id
  3. matlab的基本语法规则_Matlab基本语法与指令
  4. MySQL数据库环境变量设置
  5. pdf.js、pdfobject.js实现pdf文件预览
  6. android soundpool 封装,android使用SoundPool 播放短音效
  7. 微信音频通话数据保存服务器,微信语音通话怎么录MP3音频文件
  8. java画地图热力图_高德地图+热力图+AJAX(SSM)
  9. AD快捷键无法使用的解决方法
  10. zabbix3.2短信告警脚本