ArcGIS教程:半变异函数与协方差函数

  半变异函数和协方差函数将邻近事物比远处事物更相似这一假设加以量化。半变异函数和协方差都将统计相关性的强度作为距离函数来测量。

  对半变异函数和协方差函数建模的过程将半变异函数或协方差曲线与经验数据拟合。目标是达到最佳拟合,并将对现象的认知纳入模型。之后模型便可用于预测。

  在拟合模型时,浏览数据中的方向自相关。基台、变程和块金是模型的重要特征。如果数据中有测量误差,请使用测量值误差模型。跟踪这一链接来了解如何将模型与经验半变异函数拟合。

  半变异函数

  半变异函数定义为

  γ(si,sj) = ? var(Z(si) - Z(sj)),

  其中 var 是方差。

  如果两个位置 si 和 sj,在 d(si, sj) 的距离测量上彼此相近,那么您会希望这两个位置相似,以便缩小两个位置的差值 Z(si) - Z(sj) 的大小。当 si 和 sj 距离逐渐增大时,它们变得越来越不相似,它们的值 Z(si) - Z(sj) 的差异也会增大。在下图中可以看到这一情况,其中显示了典型半变异函数的解析图。

  

  请注意,差值的方差会随距离的增大而增加,因此可以将半变异函数视为相异度函数。与这一函数经常关联的术语也可用在 Geostatistical Analyst 中。半变异函数在其呈平稳状态时所达到的高度称为基台。它通常由两部分组成:原点处不连续(称为块金效应)和偏基台;二者一起形成基台。块金效应可以细分为测量误差和微刻度变化。块金效应就是测量误差和微尺度变化的和,由于任一组件都可为零,因此块金效应可以完全由一个组件或另一个组件形成。变程是半变异函数达到平稳基台处的距离。

  协方差函数

  协方差函数定义为

  C(si, sj) = cov(Z(si), Z(sj)),

  其中 cov 是协方差。

  协方差是相关性的缩放版。因此当两个位置,si 和 sj 彼此相近时,您会希望这两个位置相似,而他们的协方差(相关性)会变大。当 si 和 sj 距离逐渐增大时,它们变得越来越不相似,并且它们的协方差会变为零。在下图中可以看到这一情况,下图显示典型协方差函数的解析图。

  

  请注意,协方差函数随距离的增大而减小,因此可将其视为一种相似度函数。

  半变异函数和协方差函数之间的关系

  在半变异函数和协方差函数之间存在以下关系:

  γ(si, sj) = sill - C(si, sj),

  从图中可看出该关系。由于这一相等关系,您可以在 Geostatistical Analyst 中使用两种函数中的任一种来执行预测。(Geostatistical Analyst 中所有半变异函数都拥有基台。)

  半变异函数和协方差不是任意函数皆可。为使预测具有非负的克里金标准误差,只有部分函数可以用作半变异函数和协方差。Geostatistical Analyst 提供了多种可接受的选项,您可以为数据尝试不同的选项。您也可以通过同时添加多个模型的方式获得模型 - 此构造提供有效的模型,可以在 Geostatistical Analyst 中添加其中的最多四个模型。有一些当半变异函数存在时,协方差函数却不存在的实例。例如,有一个线性半变异函数,但它没有基台,并且没有相对应的协方差函数。Geostatistical
Analyst 中仅使用带有基台的模型。在选择“最佳”半变异函数模型时,没有必须遵守的规则。您可以查看经验半变异函数或协方差函数并选择看起来适合的模型。你也可以使用验证和交叉验证作为指南。

时间: 2024-10-09 06:29:48

ArcGIS教程:半变异函数与协方差函数的相关文章

ArcGIS教程:使用方向半变异函数和协方差函数考虑各向异性

由于您是在二维空间内进行操作,因此可能希望半变异函数和协方差函数不仅随距离发生改变,还会随方向发生改变.这被称为各向异性. 考虑两个点 si 和 sj,以及分开这两个点的矢量,这个矢量记为 si - sj.该矢量在 x 轴和 y 轴上都具有距离.或者,在极坐标中可将该矢量视为具有距离和角度.此处,各向异性用于描述半变异函数;思路与协方差函数类似. 各向同性模型在所有方向上都是相同的,而各向异性模型在某些方向上达到基台的速度高于另一些方向.达到基台的较长轴的长度称为主变程,达到基台的较短轴的长度称

ArcGIS教程:根据经验半变异函数拟合模型

半变异函数/协方差建模是空间描述和空间预测之间的关键步骤.地统计的主要应用是预测未采样位置处的属性值(克里金法). 经验半变异函数和协方差可提供有关数据集的空间自相关的信息.但是,不提供所有可能方向和距离的信息.因此,为确保克里金法预测的克里金法方差为正值,根据经验半变异函数/协方差拟合模型(即连续函数或曲线)是很有必要的. 经验半变异函数/协方差值的不同视图 地统计向导可提供经验半变异函数值的三种不同视图.可以使用任意数量(一个.两个或全部三个)的视图来帮助您根据数据拟合模型.默认视图显示了已

ArcGIS教程:不同方向的经验半变异函数

有时,测量位置的值会包含一种方向影响,这种方向影响可以通过统计方式进行量化但可能无法通过已知的可识别过程予以说明.此方向影响称为各向异性.各向异性用于分析样本是否在不同方向上呈现出不同的范围.容差角度用于确定在哪个角度下包括或排除接近点,直至达到带宽.带宽用于指定当确定将在半变异函数中绘制哪些点对时应采用的搜索宽度. 条柱中的点是分别位于特定距离和方向内的位置对.从概念上讲,可以将方向分组视为对将在分组过程中绘制的点对的限制,或者视为绘制所有点对并仅考虑表示特定方向的图形部分.以下场景描述了带宽

ArcGIS教程:半变异函数建模

当进行半变异函数建模时,可对自相关性进行检查和量化.在地统计中,这称为空间建模,也称为结构分析或变异分析.在半变异函数的空间建模中,可以从经验半变异函数图开始,计算为, Semivariogram(distance h) = 0.5 * average [ (value at location i– value at location j)2] (所有成对位置的相隔距离为 h).该公式涉及到计算配对位置的差值平方的一半.快速绘制所有配对则变得难以处理.并不绘制每个配对,而是将配对分组为各个步长条

了解半变异函数:变程、基台和块金

变程和基台 块金 半变异函数显示测量采样点的空间自相关.绘制每对位置,然后根据这些位置拟合模型.通常使用某几个特征来描述这些模型. 变程和基台 查看半变异函数的模型时,您将注意到模型会在特定距离处呈现水平状态.模型首次呈现水平状态的距离称为变程.比该变程近的距离分隔的样本位置与空间自相关,而距离远于该变程的样本位置不与空间自相关. 半变异函数模型在变程处所获得的值(y 轴上的值)称为基台.偏基台等于基台减去块金. 半变异函数示例 块金 从理论上讲,在零间距(步长 = 0)处,半变异函数的值是 0

ArcGIS教程:地统计工作流

这一主题将介绍地统计研究的概化工作流以及主要步骤.正如什么是地统计中所述,地统计是用于分析和预测与空间现象或时空现象相关联的值的统计数据类.ArcGIS Geostatistical Analyst 提供一套工具,该工具允许构建使用空间(和时间)坐标的模型.这些模型可以应用于各种情况并通常用于生成未采样位置的预测,也可用于生成这些预测的不确定性的度量值. 与几乎所有数据驱动研究相同,第一步是仔细检查数据.这通常从映射数据集.使用可使数据集可能显示的重要特征清晰可见的分类和颜色方案开始,例如,从北

ArcGIS教程:Geostatistical Analyst 应用示例

Geostatistical Analyst 满足多种不同应用的需求.以下是 Geostatistical Analyst 的一小部分应用示例. 探索性空间数据分析 Geostatistical Analyst 用于使用研究区中已测量的采样点为同一区域内其他未测量位置创建准确预测.Geostatistical Analyst 中包括的探索性空间数据分析工具用于评估数据的统计属性,比如空间数据变异性.空间数据相关性和全局趋势. 以下示例中使用了多个探索性空间数据分析工具来研究喀尔巴阡山中监测站处获

ArcGIS教程:合并半变异函数模型

通常有两个或更多处理过程将决定某一现象的空间分布.例如植被量(生物量)可能与高程和土壤湿度相关.如果这种关系已知,则能够使用协同克里金法预测生物量.可将生物量的测量值用作数据集 1,高程用作数据集 2,土壤湿度用作数据集 3.可能要根据每个数据集拟合不同的半变异函数模型,因为每个数据集的空间结构各不相同.也就是说,球面模型可能最适合高程,指数模型可能最适合土壤湿度,而这两种模型的组合可能最适合生物量.然后可通过一种最适合数据结构的方式将模型合并. 但有时并不知道在某一现象中将确定空间结构的因素的

ArcGIS教程:区域插值参数

构建有效模型 与所有地统计插值方法一样,区域插值中的预测准确性取决于模型的准确性.了解此事项后,在地统计向导 中构建有效模型时应多加注意. 由于 ArcGIS Geostatistical Analyst 扩展模块中的区域插值是通过克里金框架实现的,因此交互式变异分析是构建模型的重要步骤.通常很难从视觉上判断协方差曲线的质量,因此为每个经验协方差(下图中的蓝十字)提供了置信区间(下图中的红色垂直线段).如果正确指定了协方差模型,预计有 90% 的经验协方差落在置信区间内.在下图中,12 个经验协