GIS应用技巧之空间插值理论介绍
一、空间插值概论
空间插值常用于将离散点的测量数据转换为连续的数据曲面,以便探究空 间现象的分布模式,该方法通常用来分析地区环境污染、地区降水量、地区气候 变化、资源利用程度、公共基础设施影响效应等。空间插值方法分为两类:一类 是确定性方法,另一类是地质统计学方法。确定性插值方法是基于信息点之间 的相似程度或者整个曲面的光滑性来创建一个拟合曲面,比如反距离加权平均插值法(IDW)、趋势面法、样条函数法等。地质统计学插值方法是利用样本点 的统计规律,使样本点之间的空间自相关性定量化,从而在待预测的点周围构建样本点的空间结构模型,比如克立格(Kriging)插值法。确定性插值方法的 特点是在样本点处的插值结果和原样本点实际值基本一致,若是利用非确定性 插值方法的话,在样本处的插值结果与样本实测值就不一定一致了,有的相差 甚远。在实际的GIS 应用过程中,往往要根据实际的需要,综合考虑插值的精 度、效率等因素,来选择不同的插值方法。这里着重介绍下反距离加权,样条函 数和克里金插值法。
(1)反距离加权法(IDW)
GIS 中最常用的空间内插方法之一是反距离加权法,是以插值点与样 本点之间的距离为权重的插值方法,插值点越近的样本点赋予的权重越大,其 权重贡献与距离成反比。可表示为:
其中 Z 是插值点估计值,Zi(i =1、2、3、、、、)是实测样本值,n为参与计算 的实测样本数,Di为插值点与第 i 个站点间的距离,p是距离的幂,它显著影响 内插的结果,它的选择标准是最小平均绝对误差。
(2)样条函数内插法
样条函数是一个分段函数,进行一次拟合只有少数点拟合,同时保证曲线 段连接处连续,这就意味着样条函数可以修改少数数据点配准而不必重新计算 整条曲线。样条函数的一些缺点是:样条内插的误差不能直接估算,同时在实 践中要解决的问题是样条块的定义以及如何在三维空间中将这些“块”拼成复 杂曲面,又不引入原始曲面中所没有的异常现象等问题。
(3)克里格插值法
克里格法是 GIS 软件地理统计插值的重要组成部分。这种方法充分吸收 了地理统计的思想,认为任何在空间连续性变化的属性是非常不规则的,不能 用简单的平滑数学函数进行模拟,可以用随机表面给予较恰当的描述。这种连 续性变化的空间属性称为“区域性变量”,可以描述像气压、高程及其他连续性 变化的描述指标变量。地理统计方法为空间插值提供了一种优化策略,即在插 值过程中根据某种优化准则函数动态的决定变量的数值。Kriging 插值方法着 重于权重系数的确定,从而使内插函数处于最佳状态,即对给定点上的变量值 提供最好的线性无偏估计。
GIS 中有七类克里格法,下表是这七种方法的名称和适用范围:
克里格法的优点是以空间统计学作为其坚实的理论基础,可以克服内插中误差难以分析的问题,能够对误差做出逐点的理论估计;不但能估计测定参数 的空间变异分布,而且还可以估算估计参数的方差分布。其缺点是计算步骤较 繁琐,计算量大,且变异函数有时需要根据经验人为选定。
GIS应用技巧之空间插值理论介绍相关推荐
- 一、DMSP/OLS、NPP/VIIRS等夜间灯光数据之GDP空间化的理论介绍
一.前言 自上世纪90年代以来,DMSP/OLS夜间灯光数据问世,由于和其他遥感数据不同,该数据可以探测到地球夜间微弱的灯光.随着数据的积累,可以获取长时间序列的数据之后,随后常用于人类社会经济指标的 ...
- python2 使用arcpy自定义GIS工具箱,实现excel读取、空间插值
一.实验目的 1.掌握Arcpy基本语法.功能函数 2.学会使用Arcpy自定义GIS工具箱 3.使用arcpy进行空间插值 4.自定义创建的Arcmap制图模板,使用arcpy调节图例.注记.比例尺 ...
- 开源的前端GIS空间分析库介绍 (三)turf与ol结合
前言 turf是mapbox出品的前端空间分析库,官网:http://turfjs.org/ turf库中包含的空间分析计算功能比较多,也非常简单易用.相比于jsts,turf的官方文档维护的非常好, ...
- 开源的前端GIS空间分析库介绍 (一)jsts与turf
文章目录 1 前言 2 JSTS 3 turf 4 安装使用 4.1 jsts 4.1.1 直接引入 4.1.2 NPM 4.2 turf 4.1.1 直接引入 4.1.2 NPM 5 空间分析 5. ...
- GEE:克里金 Kriging 空间插值(以陕西省2013年生物量为例)
作者:CSDN @ _养乐多_ 本文记录了在Google Earth Engine(GEE)平台上进行 Kriging 插值的介绍和代码案例.本文通过选取的2013年陕西省生物量样本点数据为例,利用 ...
- 完美解决ArcGIS10.x栅格空间插值报错无法进行和导出插值栅格结果出错的问题
不少GISer,尤其是初学者,对栅格数据的操作,主要是栅格插值无法进行,或者是插值后的栅格结果无法导出,或者导出结果为空,或者是栅格格式无法转换,今天刘一哥GIS就帮您排忧解难吧. 其实,关于栅格(尤 ...
- gstat | 空间插值(三)——克里金插值之泛克里金和简单克里金
本篇接着上篇继续介绍克里金插值.首先加载相关工具包和上篇使用的示例数据: library(gstat) library(sf) library(tidyverse) library(readxl) l ...
- gstat | 空间插值(二)——克里金插值之普通克里金
说明:昨天的推文误把可吸入颗粒物当作PM2.5,实应该为PM10,这里修正后重发. 从本篇开始计划分三篇介绍克里金插值.与反距离权重插值不同,克里金插值是无偏估计,其中也涉及到模型估计.本篇先对普通克 ...
- gstat | 空间插值(一)——反距离权重插值;使用ggplot2绘制地图
本篇既是空间插值系列的第一篇推文,也是ggplot2工具包系列推文中的一篇.空间插值使用的工具包是gstat,该工具包主要用于地统计分析. library(gstat) 示例数据来自HSAR工具包: ...
- 白话空间统计二十九:空间插值(二)
前文再续,书接上一回--上回书说到,空间插值可以对数据进行估算,这一张我们来具体说说空间插值的一些概念. 首先,来说说插值的输入与输出问题. 首先输入的肯定是用来进行插值的观测点数据,一般来说都是点要 ...
最新文章
- 第九章-构建图形用户界面
- 使用 jQuery 查询属性不包含 disabled 的 input radio
- LeetCode 684. 冗余连接(并查集)
- 全开源纵横支付多平台多功能支付系统源码
- 层间距离对ct图像的影响_CT专业知识掌握基础、了解原理
- 支付功能如何软件测试
- 视觉平台搭建——光源选型介绍
- Docker安装迅雷下载工具实现远程下载
- QGIS教程01:为什么要用QGIS?
- 大型网站技术架构 读书笔记2 核心架构要素
- RNA剪接体 Spliceosome | 冷冻电镜 | 结构生物学
- mysql排序之if(isnull(字段名),0,1),fild 或者 if(isnull(字段名),1,0),fild
- PaddlePaddle证件照换底换大小工具
- 希望各位大神能帮我指导一下!!!感激不尽!!!
- 简单的几句话让你理解”什么是备份、容灾、集群、负载均衡”
- 淘宝开店首页装修教程分享
- LeetCode-86
- 【Chrome 浏览器自带谷歌翻译用不了】
- libreCAD源码阅读笔记2
- 基于imx8的成熟360全景方案 ,NVP6324+IMX8,N4+IMX8