ArcGIS教程:地统计模拟的重要概念

  模拟概念

  模拟在广义上是指使用模型复制现实的过程。在地统计中,模拟是随机函数(表面)的实现,其与生成该模拟的样本数据拥有相同的地统计要素(使用均值、方差和半变异函数来度量)。更具体地说,高斯地统计模拟 (GGS) 适用于连续数据,并假设数据或数据的变换具有正态(高斯)分布。GGS 所依托的主要假设是数据是静态的 - 均值、方差和空间结构(半变异函数)在数据空间域上不发生改变。GGS 的另一个主要假设是建模的随机函数为多元高斯随机函数。

  同克里金法相比,GGS 具有优势。由于克里金法是基于数据的局部平均值的,因此,其可生成平滑的输出。另一方面,GGS 生成的局部变异性的制图表达比较好,因为 GGS 将克里金法中丢失的局部变异性重新添加到了其生成的表面中。对于由 GGS 实现添加到特定位置的预测值中的变异性,其平均值为零,这样,很多 GGS 实现的平均值会趋向于克里金预测。下图对此概念进行了说明。各种实现以一组堆叠输出图层的形式表示出来,并且特定坐标位置的值服从高斯分布,其平均值等于该位置的克里金估计值,而扩散程度则由该位置上的克里金法方差给出。

  

  提取值到表工具可以用来为上图中的图形生成数据,在对 GGS 生成的输出进行后处理时该工具也很有用。

  对 GGS 的使用在地统计实际操作中日益呈现出一种趋势,它不是追求获得每个未采样位置的最佳无偏预测结果(正如克里金法所体现的),而是强调对决策分析和风险分析的不确定性的特证描述,这样更适合于呈现数据中的全局趋势 (Deutsch and Journel 1998, Goovaerts 1997)。模拟还会克服克里金估计值中的条件偏差带来的问题(高值区域预测值通常偏低,而低值区域预测值通常偏高)。

  对于所研究属性的空间分布,地统计模拟可为其生成多个具有同等可能性的制图表达。可基于这些制图表达来测量未采样位置的不确定性,这些未采样位置在空间上被一起选取,而不是逐个被选取(如同通过克里金法方差进行测量一样)。此外,克里金法方差通常独立于数据值,且通常不能用作估计精度的测量值。另一方面,可以通过使用多个模拟实现(该实现用呈正态分布的输入数据通过简单克里金模型进行构建,即,数据呈正态分布或已使用常态得分变换或其他类型的变换对数据进行了变换)为未采样位置的估计值构建分布来测量估计精度。对于使用估计数据值的风险评估和决策分析而言,这些不确定性的分布很关键。

  GGS 假设数据呈正态分布,但在实际中,很少会出现这种情况。对数据执行常态得分变换,使得数据符合标准正态分布(均值 = 0,方差 = 1)。然后对此正态分布数据进行模拟,并对结果做反向变换,以便以原始单位获得模拟输出。对正态分布数据使用简单克里金法时,该克里金法所提供的克里金估计值和方差可完全定义研究区域中每个位置的条件分布。这样,您可以在只知道每个位置的这两个参数的情况下绘制随机函数(未知采样表面)的模拟实现,这也是GGS 基于简单克里金模型和正态分布数据的原因。

  “高斯地统计模拟”工具支持两种类型的模拟:

  1、条件模拟遵循数据值(除非克里金模型中包含测量误差)。由于模拟会在格网像元中心生成值,因此,如果此值与采样点的位置不完全对应,则采样位置的测量值与模拟值可能会不同。条件模拟也将以平均方式(即,在很多实现上平均)复制数据的均值、方差和半变异函数。模拟表面看起来很像克里金预测地图,但其将显示更多的空间变异性。

  2、非条件模拟不遵循数据值,但会以平均方式复制数据的均值、方差和半变异函数。模拟表面所显示的空间结构类似于克里金地图,但输入数据中存在高值或低值的地方不一定会出现高值和低值。

时间: 2024-10-09 22:49:49

ArcGIS教程:地统计模拟的重要概念的相关文章

ArcGIS教程:“高斯地统计模拟”的工作原理(一)

高斯地统计模拟的一般工作流包括:准备数据.创建实现.在后台将结果传输到原始单元中,以及后处理结果和/或将结果用作传输函数(模型)的输入以评估模型输出中的变异性. 实现的生成方法 使用"高斯地统计模拟"时,首先要创建一个基于标准正态分布(平均值 = 0 且 方差 = 1)绘制的随机分配值的格网.然后,将协方差模型(基于在需要作为模拟的输入的简单克里金图层中指定的半变异函数)应用到栅格.这样可以确保栅格值遵循输入数据集中的空间结构.生成的栅格构成一个非条件实现,而且通过每次使用不同的包含正

ArcGIS教程:“高斯地统计模拟”的工作原理(二)

用于模拟的简单克里金模型 "高斯地统计模拟"工具可接受任何的简单克里金模型.不过,仅当输入数据(用于拟合半变异函数以及条件化实现)呈正态分布时模拟结果才有效.这可通过应用常态得分变换来实现.还应进行去聚处理,以便从已聚类的数据中获得典型直方图,而且应移除趋势,以确保平均值在整个空间域是稳定的.这些步骤如以下各图所示.图 (a) 显示简单克里金法,其中趋势移除.去聚和常态得分变换选项被选中.图 (b) 显示在进行去聚.常态得分变换和变异分析之前要从数据集中移除的趋势.图 (c) 按像元显

ArcGIS教程:3D Analyst基础知识

创建3D视图 以三维形式查看数据能为您提供一个全新的认识.通过三维视图可以深入了解通过相同数据的平面地图不易察觉的内容.例如,您不必根据配置等值线来推断是否存在山谷,您能够实际看到山谷和感到谷底和谷脊的高度差异. ArcGlobe 和 ArcScene 可用于构建多图层 3D 环境,并控制如何对各图层进行符号化.渲染各图层和在 3D 空间中定位各图层.还可以控制 3D 视图的全局属性,如照明度或垂直夸大.可以通过以下方式选择要素:使用要素的属性或要素相对于其他要素的位置,或者在场景或地球中单击各

ArcGIS教程:区域插值参数

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

ArcGIS教程:绘制数据

任何分析的第一步都是绘制和检查数据.这可以提供数据集的空间组成部分的第一印象,而且可能给出异常值和错误数据值.全球趋势和其他系数间的空间自相关的主导方向的指示,所有这些在开发正确反映感兴趣的现象的插值模型的过程中都非常重要. ArcGIS 提供很多方法来可视化数据:ArcMap 可以访问用于高亮显示数据不同方面的很多分类方案和色带,而ArcScene 可以在 3D 空间渲染数据,这在查找局部异常值和全球趋势时非常有用.尽管没有正确的方法来显示数据,下图显示了相同数据的不同渲染,从中可以看出不同方

ArcGIS教程:检查数据的分布

虽然预测图在数据不服从正态分布的情况下可能不是最佳的,但 ArcGIS Geostatistical Analyst 提供的大多数插值方法都不要求数据服从正态分布.即,更改数据形状(分布)的数据变换不要求是插值模型的一部分.不过,有些克里金方法要求数据近似服从正态分布(接近钟形曲线).特别地,利用普通克里金法.简单克里金法或泛克里金法创建的分位数图和概率图都假定数据服从多元正态分布.而且,用作地统计模拟的基础的简单克里金模型应使用服从正态分布的数据或者将常态得分变换作为模型的一部分以确保数据服从

ArcGIS教程:ArcGIS中矢量裁剪栅格图像

(1)是否需要裁剪栅格图象区域通过一个面状的shapefile表达出来? 如果可以,那么就很简单了. 在ArcMap中,调用空间分析扩展模块,将你感兴趣区的shapefile多边形图层设置为掩膜,然后在栅格计算器中重新计算一下你的图象,它就会沿掩膜裁出. 设置掩膜:空间分析工具条的下拉菜单>option里面设置 (2)用任意多边形剪切栅格数据(矢量数据转换为栅格数据) 2.1在ArcCatlog下新建一个要素类(要素类型为:多边形),命名为:ClipPoly.shp 2.2在ArcMap中,加载

ArcGIS教程:ArcGIS栅格数据的合并和剪切

1.合并:ArcToolBox->DataManagement->Raster->Mosaic. 2.剪切:在ArcMap中,调用空间分析扩展模块,将你感兴趣区的shapefile多边形图层设置为掩膜(空间分析工具条的下拉菜单>option里面设置),然后在栅格计算器中重新计算一下你的图象,它就会沿掩膜裁出. 以上教程来源地理国情监测云平台,更多ArcGIS教程及空间地理信息数据请咨询本平台,电话:010-84896208转898.官方微信号DLGQJC,为您提供3S行业知识及每日

ArcGIS教程:ArcGIS中图层添加投影及投影转换

由于GIS描述的是位于地球表面的空间信息,所以在表示时必须嵌入到一个空间参照系中,这个参照系就是坐标系它是根据椭球体等参数建立的.我国目前所用到的三种主要坐标系是:北京1954.西安1980及WGS84. 有了坐标系,即椭球体,数据图层就有了地理坐标,地理坐标是用经纬度表示球面的位置的.为了能够将地图从三维地理坐标通过投影转换成二位平面坐标,还要进行投影,这样的坐标系叫做投影坐标系.常用的投影有高斯-克吕格(Gauss-Kruger)投影,阿尔博斯(Albers Equal-Area Conic