综述指示克里格方法的原理与应用.pdfVIP

综述指示克里格方法的原理与应用 学号：2013301610318 姓名： 在土壤质量评价、大气污染物浓度分布评估等研究中，克里格 空间插值方法是一个有力的工具。但是观测数据中往往存在一些特异 值，观测数据不成(对数)正态分布，影响了变异函数的稳健性。如果用 参数地质统计学方法, 必须剔除这些特异值或者对观测值进行非 线性转换，以使观测值的概率分布满足正态，但这会影响变量空间变 异的真实信息。非参数地质统计学中的指示克里格法是处理有偏数据 的有效方法，它能在不必去除重要而实际存在的特异值的条件下处理 不同的现象，并能抑制特异值对变异函数的稳健型的影响。指示克立 格法是一种最常用的非参数地质统计学方法,它是因把对区域化变量 的研究转换为对其指示函数的研究而得名.有关指示克里格方法的研 究与应用，国内外学者已经做了很多工作，但是大部分研究都是单元 指示克里格在单一尺度下的应用。本文将主要讨论指示克里格方法的 基本原理，指示克里格的应用方法(比如多尺度指示克里格、多元指示 克里格等)以及指示克里格法的限制和不足。 1 基本原理 克里格(Kriging)插值法是空间统计分析方法的重要内容之一， 1 它是建立在半变异函数理论分析基础上的，是对有限区域内的区域化 变量取值进行无偏最优估计的一种方法。基于这种方法进行插值时， 不仅考虑了待预测点与邻近样点数据的空间距离关系，还考虑了各参 与预测的样点之间的位置关系，充分利用了各样点数据的空间分布结 构特征，使其估计结果比传统方法更精确，更符合实际，更有效避免 了系统误差的出现。 在空间统计分析方法中，可以通过选择阈值，将一个连续的变 量转换成一个值为 0 或 1 的二进制变量。比如在研究区域 D 内,Z(X)表 示采样点 X 上的采样值,设 Z 为研究区域 D 上的一个临界值(阈值), 在 D 上的每点 X ∈D 上定义一个 Z 的指示函数如下: I(Z;X) = { 1 Z(X) ≤ Z ( ) 0 Z X Z 对于指示函数I(Z;X),可以用条件概率来描述:当Xa, a=1,2, ⋯,n 为 采样点时, ( ) ( ) I X ;Z = P{Z(X ) ≤ Z|Z X = Z } a a a a 这时,某待估点 X 的指示函数估计值I * (X;Z)可以表示为： * ( ) ( ) I X;Z = P{Z(X) ≤ Z|Z X = Z , a = 1, 2,⋯,n} a a 对于采样点来说,指示值可解释为已知该点的实测值为 Za 时,该 点的真实值小于等于阈值的概率,而对于待估点,其指示函数估计值可 解释为已知待估点周围信息(样本的实测值) 时,该点的真实值小于等 于阈值的概率。 2 克里格插值法的研究工具是半变异函数，指示克里格插值法的 研究工具是指示半变异函数 (h)： 1 ( ) ( ) 2 h = E{[I X;Z - I(X + h;Z)]