什么是经验贝叶斯克里金插值法?

原创 2016年08月29日 17:17:36

一.简介

经验贝叶斯克里金法 (EBK) 是一种地统计插值方法[还有确定性插值方法],可自动执行构建有效克里金模型过程中的那些最困难的步骤。Geostatistical Analyst 中的其他克里金方法需要您手动调整参数来接收准确的结果,而 EBK 可通过构造子集和模拟的过程来自动计算这些参数。

经验贝叶斯克里金法与其他克里金方法也有所不同,它通过估计基础半变异函数来说明所引入的误差。其他克里金方法通过已知的数据位置计算半变异函数,并使用此单一半变异函数在未知位置进行预测;此过程隐式假定估计的半变异函数是插值区域的真实半变异函数。由于不考虑半变异函数估计的不确定性,其他克里金方法都低估了预测的标准误差。经验贝叶斯克里金法在地统计向导中以地理处理工具的形式提供。

二.优点和缺点

优点

  • 需要极少的交互式建模。
  • 预测标准误差比其他克里金方法更准确。
  • 可准确预测一般程度上不稳定的数据。
  • 对于小型数据集,比其他克里金法更准确。

缺点

  • 处理时间会随着输入点数、子集大小或重叠系数的增加而快速增加。应用变换还会增加处理时间,尤其是选择 K-Bessel去除趋势的 K-Bessel 作为半变异模型类型时。这些参数将在本主题的后续章节中进行介绍。
  • 处理速度比其他克里金方法慢,尤其是输出为栅格时。
  • 协同克里金法和各向异性校正不可用。
  • 对数经验变换对异常值尤其敏感。如果将该变换用于含有异常值的数据,则可能会得到大于或小于输入点值若干个数量级的预测结果。该参数将在下面的“变换”部分进行介绍。

三.半变异函数估计

与其他克里金法(使用加权最小二乘)不同,EBK 中的半变异函数参数是使用受限最大似然法 (REML) 估计的。由于 REML 对大型数据集有计算限制,输入数据首先被分为多个特定大小的重叠子集(默认为每子集 100 个点)。在每个子集中,按以下方式估计半变异函数:

  1. 通过子集中的数据估计半变异函数。
  2. 将此半变异函数用作模型,新数据会在子集的每个输入位置进行无条件模拟。
  3. 通过已模拟的数据估计新的半变异函数。
  4. 将步骤 2 和步骤 3 重复执行指定次数。在每次重复中,步骤 1 中估计的半变异函数用于模拟输入位置的一组新数据,已模拟的数据用于估计新的半变异函数。

此过程将为每个子集创建大量半变异函数,并且在将它们绘制在一起时,结果是按密度着色的半变异函数分布(蓝色越深,通过该区域的半变异函数就越多)。经验半方差用蓝色十字表示。此外,分布的中值用红色实线表示,25% 和 75% 百分数值用红色虚线表示,如下图所示。


每个子集中模拟的半变异函数数量默认为 100,其中每一个半变异函数都是子集的真实半变异函数的估计。

对于每个预测位置,都使用新的半变异函数分布计算预测,该分布是通过对点邻域内半变异函数光谱中的各个半变异函数进行基于可能性的采样而生成的。例如,如果预测位置在三个不同的子集(由搜索邻域指定)中均有邻域,则将使用每个子集的某些模拟的半变异函数计算预测;这些半变异函数根据其可能性值以概率的方式选择。

四.克里金模型

经验贝叶斯克里金法与 Geostatistical Analyst 中的其他克里金方法不同,它使用固有的随机函数作为克里金模型。

其他克里金模型假定过程遵循一个总体平均值(或指定趋势),并且各种变化均围绕该平均值。较大的偏差将向平均值拉回,因此值不会偏差过大。但是,EBK 不会呈现出趋于总体平均值的趋势,因此较大偏差变大变小的可能性相同。因此,固有的随机函数本身就会对数据的趋势进行校正。

半变异函数模型

对于给定距离 h,经验贝叶斯克里金法支持以下半变异函数:

  • 幂函数
    • γ(h)= Nugget + b|h|α
  • 线性函数
    • γ(h)= Nugget + b|h|
  • 薄板样条函数
    • γ(h)= Nugget + b|h2|*ln(|h|)

块金值和 b(坡度)必须为正值,而 α(幂)必须介于 0.25 和 1.75 之间。在这些限制下,使用 REML 估计参数。这些半变异函数模型没有变程或基台参数,因为函数没有上限。

在 EBK 中,可以分析参数估计的经验分布,因为在每个位置都估计了多个半变异函数。单击块金值坡度选项卡可显示关联参数的分布。下图显示了前一图片中显示的模拟半变异函数的半变异函数参数分布:



半变异函数

所有的地统计方法都假设空间自相关(即距离较近的事物比距离较远的事物更相似),半变异函数则定义这种相似性如何随距离的增加而减少。某些半变异函数(如指数函数)假设相似性快速减少。另一方面,消减半变异函数模型则假设相似性缓慢减少。即使包含相同的块金、变程、基台,这两种半变异函数也会以截然不同的方式定义相似性的减少。要获得可靠的结果,关键在于选择与您的现象的行为方式最为匹配的半变异函数。半变异函数模型的可用情况取决于您如何设置变换。

如果将变换设置为,以下半变异函数模型将可用:

  • 幂函数(默认)
  • 线性函数
  • 薄板样条函数

如果将变换设置为经验对数经验,则以下半变异函数模型将可用:

  • 指数函数(默认)
  • 去除趋势的指数函数
  • 消减函数
  • 去除趋势的消减函数
  • K-Bessel
  • 去除趋势的 K-Bessel

除应用一阶趋势移外,这三种去除趋势的半变异函数模型与未去除趋势的对应模型是相同的。移除趋势对计算速度的影响可以忽略不计。可使用趋势分析 ESDA 工具检查趋势是否存在。

每种模型的优点和缺点

每种半变异函数都具有其各自的优缺点。选择半变异函数时,应考虑函数模型的计算时间和灵活性(精确覆盖广大范围数据集的能力):

  • 幂函数
    • 优点:相对较快和灵活。通常是平衡性能和精度的安全选择。
    • 缺点:较其他选择而言灵活性较差,速度也较慢。
  • 线性函数
    • 优点:速度非常快。
    • 缺点:灵活性最差的模型。
  • 薄板样条函数
    • 优点:速度非常快。在存在强大趋势时最有效。
    • 缺点:灵活性较差,尤其在不存在任何趋势时。
  • 指数函数:
    • 优点:可灵活进行变换。速度比 K-Bessel 和去除趋势的 K-Bessel 更快。
    • 缺点:半变异函数的形状不灵活。速度比幂函数、线性函数和薄板样条函数更慢。
  • 去除趋势的指数函数
    • 优点:可灵活进行变换。速度比 K-Bessel 和去除趋势的 K-Bessel 更快。移除一阶趋势。
    • 缺点:半变异函数的形状不灵活。速度比幂函数、线性函数和薄板样条函数更慢。
  • 消减函数
    • 优点:可灵活进行变换。速度比 K-Bessel 和去除趋势的 K-Bessel 更快。
    • 缺点:半变异函数的形状不灵活。速度比幂函数、线性函数和薄板样条函数更慢。
  • 去除趋势的消减函数
    • 优点:可灵活进行变换。速度比 K-Bessel 和去除趋势的 K-Bessel 更快。移除一阶趋势。
    • 缺点:半变异函数的形状不灵活。速度比幂函数、线性函数和薄板样条函数更慢。
  • K-Bessel
    • 优点:最灵活、最精确。
    • 缺点:计算时间最长。
  • 去除趋势的 K-Bessel
    • 优点:最灵活、最精确。移除一阶趋势。
    • 缺点:计算时间最长。

选择半变异函数

大多数情况下,应根据下列条件明确选择半变异函数:

  • 如果您愿意耐心等待获取最准确的结果,应选择 K-Bessel 或去除趋势的 K-Bessel。应视是否存在趋势来选择相应的函数。
  • 如果您需要快速获得结果,并愿意牺牲一定的精度,应选择线性函数或薄板样条函数。如果不存在趋势或趋势很弱,则更适合选择线性函数。
  • 如果需要平衡精度和速度,幂函数是个不错的选择。
  • 如果需要进行变换,但又无法长时间等待结果,应选择指数函数或消减函数(或未去除趋势的对应函数)。您应选择与地统计向导中的经验半方差最为匹配的函数(如下所述)。另外还应考虑进行交叉验证

如果试图在指数函数、消减函数及其去除趋势的对应函数之间进行选择,应选择在视觉上与经验半方差最拟合的半变异函数(见下图中的蓝色十字)。理想情况下,经验半方差应落在半变异函数光谱的中间。例如,在下图中,蓝色十字没有落在半变异函数光谱的中间(大多都落到光谱的顶部):







版权声明:本文为生命奇迹泉原创文章,转载请注明出处生命奇迹泉http://blog.csdn.net/shengmingqijiquan

经验贝叶斯克里金——基于ArcGIS 10.1Geostatistical Analyst

经验贝叶斯克里金 基于ArcGIS 10.1Geostatistical Analyst 作者:康斯坦丁 克里沃切科(Konstantin Krivoruchko)  Esri软件研...
  • Esri_GovIndTec_HB
  • Esri_GovIndTec_HB
  • 2013年04月17日 10:31
  • 3969

python的克里金算法,使用pyKriging模块

python版本:2.7.11pykriging网址:http://www.pykriging.com/ 安装pykriging的其他必备库:numpy,cython,matplotlib,pykri...
  • wildder
  • wildder
  • 2016年07月15日 16:40
  • 2174

ArcGIS教程:什么是经验贝叶斯克里金法?

  • 2014年12月16日 10:15
  • 137KB
  • 下载

ArcGIS教程:不同的克里金模型

克里金方法依赖于数学模型和统计模型。通过添加包含概率的统计模型,可将克里金方法从空间插值的确定性方法中描述的确定性方法中分离出来。对于克里金法,您会将某种概率与预测值相关联;也就是说,这些值不能完全基...
  • u010687924
  • u010687924
  • 2015年01月28日 15:19
  • 891

数值分析-拉格朗日插值法python3实现绘图

import matplotlib.pyplot as plot import numpy as np print('拉格朗日插值法:') """ 10 拉格朗日插值法:(曲线差值法) 1.定义:找到...
  • zzy123zzy123_
  • zzy123zzy123_
  • 2017年10月26日 21:01
  • 281

克里金(Kriging)插值的原理与公式推导

转自:http://xg1990.com/blog/archives/222 学过空间插值的人都知道克里金插值,但是它的变种繁多、公式复杂,还有个半方差函数让人不知所云 本文讲简单介绍基本克里...
  • tsroad
  • tsroad
  • 2015年12月07日 14:35
  • 14745

ArcGIS教程:经验半变异函数

  • 2014年09月25日 11:50
  • 140KB
  • 下载

克里金模型

1.模型建立 [dmodel,perf]=dacefit(S,Y,regr,corr,theta0) [dmodel,perf]=dacefit(S,Y,regr,corr,theta0,lob,...
  • m62260849
  • m62260849
  • 2017年05月11日 11:29
  • 779

通过C#代码实现空间离散点的克里金(kriging)插值(二) 代码实现

程序下载本文的程序下载:克里金插值及DEM等高线生成回顾算法流程 求取已知点的距离以及点对的半方差 筛选第一步求取的结果,计算出几个均值点,用于拟合 选定拟合模型,为了方便代码实现,我选择了指数模型 ...
  • zhangkechun
  • zhangkechun
  • 2016年01月06日 21:09
  • 5066

通过C#代码实现空间离散点的克里金(kriging)插值(一) 计算原理

克里金插值较为复杂,但效果也是比较好的。为了能够通过代码实现克里金插值的过程,首先需要阅读ArcGIS中的帮助文档,以及在ArcGIS中实际操作克里金插值,以了解其详细的计算过程。...
  • zhangkechun
  • zhangkechun
  • 2015年12月30日 23:45
  • 5450
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:什么是经验贝叶斯克里金插值法?
举报原因:
原因补充:

(最多只允许输入30个字)