地统计私分析方法.pptVIP

例1 从上图可以看出，空间上有些点，由于某种原因没有采集到。如果没有缺失值，可直接对正方形网格数据结构计算变异函数；在有缺失值的情况下，也可以计算变异函数。只要“跳过”缺失点位置即可。 首先计算南北方向上的变异函数值，由变异函数的计算公式可得 同样计算出： 同理可以计算东西方向和西北—东南方向上的变异函数。最后，得到南北方向和西北—东南方向上的变异函数计算结果见下表： ①纯块金效应模型:其一般公式为 式中：c00，为先验方差。该模型相当于区域化变量为随机分布，样本点间的协方差函数对于所有距离h均等于0，变量的空间相关不存在。 ②球状模型:其一般公式为 式中：c0为块金（效应）常数;c为拱高;c0+c为基台值;a为变程。当c0=0，c=1时，称为标准球状模型。球状模型是地统计分析中应用最广泛的理论模型，许多区域化变量的理论模型都可以用该模型去拟合。 ③指数模型:其一般公式为 式中：c0和c意义与前相同，但a不是变程。当h=3α时， ，即 ，从而指数模型的变程ɑ约为3ɑ。当c0=0，c=1时，称为标准指数模型。 ④高斯模型:其一般公式为 式中：c0和c意义与前相同，a也不是变程。当 时 ，即 ，因此高斯模型的变程ɑ约为 .当c0=0,c=1时，称为标准高斯函数模型。 ?幂函数模型:其一般公式为 式中：θ为幂指数。当θ变化时，这种模型可以反映在原点附近的各种性状。但是θ必须小于2，若θ≧2，则函数r(-h)就不再是一个条件非负定函数了，也就是说它已经不能成为变异函数了。 ⑥对数模型:其一般公式为 显然，当h→0,logh→-∞，这与变异函数的性质 r(h)≧0不符。因此，对数模型不能描述点支撑上的区域化变量的结构。 ?线性有基台值模型:其一般公式为 式中:该模型的变程为a，基台值为c0+ c。 ⑧线性无基台值模型:其一般公式为 从式中可以看出，该模型没有基台值，也没有变程。 例2:某地区降水量是一个区域化变量，其变异函数 的实测值及距离h的关系见下表，下面我们试用回归分析方法建立其球状变异函数模型。 如果记 ，则可以得到线性模型 根据表中的数据，对上式进行最小二乘拟合，得到 计算可知，上式的显著性检验参数F=114.054，R2=0.962，可见模型的拟合效果是很好的。 3)普通克立格法 首先假设区域化变量Z(x)满足二阶平稳和本征假设，其数学期望为m，协方差函数c(h)及变异函数r(h)存在。即 假设在待估块段V的邻域 内，有一组n个已知样本 v(xi)(i=1,2,…,n)，其实测 值为Z(xi)(i=1,2,…,n) 克立格法的目标就是求一组权重系数λi(i=1,2,…,n),使得加权平均值 成为待估地段V的平均值ZV（x0）的线性、无偏最优估计量，即克立格估计值。为此，需要满足以下两个条件： 在变异函数存在的条件下，根据协方差与变异函数的关系：c(h)=c(0)-r(h)或r(h)=c(0)-c(h)，也可以用变异函数表示普通克立格方程组和克立格估计方差，即 用变异函数表示普通克立格方程组和克立格方差 例3 在例1中，假设降水量的变异函数为例2中的函数，它是一个各向同性的二维球状模型，已知四个观测点x1，x2，x3，x4的观测值分别为Z(x1)=37(mm)、Z(x2)=42(mm)、Z(x3)=36(mm)、Z(x4)=35(mm),试用普通克立格法内插估计观测点x0的降水量值Z(x0)。 根据普通克立格法的基本原理，我们知道，Z(x0)估计的基本公式应该是 因为i=0,1,2,3,4;j=1,2,3,4,故 当i=j时，c11=c22=c33=c44=c(0)=c0+c=2.048+1.154=3.202 当i≠j时，由 ，及根据克立格矩阵的对称性，得 观测点x0的降水量的克立格估计值为： 克立格估计方差为 1）克立格法概述 (1) 适用条件 变异函数和相关分析的结果表明区域化变量存在空间相关性。其实质是利用区域化变量的原始数据和变异函数的结构特点，对未采样点的区域化变量的取值进行线性无偏、最优估计。 (2) 克立格法的类型 ?普通克立格法； ?泛克立格法； ?协同克立格法； ?对数正态克立格法； ?指示克立格法； ?折取克立格法等。 克立格法是一簇空间局部插值模型的总称。 2）克立格估计量 对于研究区域内任一点x的测量值Z(x)，其估计值 的估