地学统计_第四章.pptVIP

第四章 变异函数和拟合模型 资源与环境学院 杨勇 目录 半方差函数 影响半方差函数的因素 变异函数的理论拟合模型 变异函数（半方差函数）r(h)是地统计学的基石 4.1 半方差函数—半方差云图 4.1 半方差函数—经验半方差函数 在实践中，一般是将滞后距（h）划分为不同的级别，计算每个滞后级别上的平均半方差值，进而得到经验半方差值和经验半方差图来分析数据的空间关系。 一般地，在实际计算时，假设步长为lag，当前滞后级别为n(n为正整数)，则h=n*lag，应该这样处理： （1）研究区所有点，找到点对（Pi,Pj），其符合条件：(n-1)*lagdis (Pi,Pj)=n*lag，它们之间的距离记为DISi （2）计算[z(pi)-z(pj)]2,记为Si. （3）设找到N(h)个这样的点对，计算平均距离 （4）计算 ，为n滞后级别上的经验半方差值。 （5）将各个级别的（havg，r*(havg)），绘制在图上，形成经验半方差图 课堂练习 假设给出某地土壤元素含量采样数据，其格式为（x,y,k），共N个采样点，其中，x和y分别表示采样点的横纵坐标值，k为实验室测得的元素含量值，步长为h，试给出相应的程序步骤，计算该区域内土壤元素的经验半方差图。（需计算的最大滞后据为样点间最大距离的一半） 变异函数的一些现象 1、跃迁现象 变异函数的一些现象 变异函数的一些现象 3、各向异性 4.2 影响半方差函数的因素 一、统计方面 1、统计分布：特异值的影响，极大或极小值越靠近区域中间，影响越大。采样数据的不均匀分布（h-散点图）。 2、屏蔽：特异值可能是错误的值，也可能是真实的值，可区分对待。如果检查发现时错误的值，必须去除。但如果是真实的值，可以去除，也可以保留（如污染区的污染物含量值）。 3、偏斜：数据的偏态分布，也会造成方差的不稳定 4.2 影响半方差函数的因素 二、样品大小和设计 1、更多的样品能增强经验半方差值的稳定性，但考虑到工作量和资金，实际取样不可能无限取点，一般要求在变程a以内，各距离上的点对数据部应小于20对，有的学者认为不应小于30对。 2、在采样的设计上，一般来讲，规则格网采样更有助于后续的分析。 4.2 影响半方差函数的因素 三、采样间距：应大致地了解变异的空间尺度以选择一个合适的采样密度。 （1）对于完全未知的领域，不知道变异的尺度或模式，应预采样，找到合适的采样间距 （2）若先前的采样间距造成经验半方差函数呈纯块金效应，没有明显的空间自相关性，则先前的采样间距过大，应缩小。 （3）若依先前的采样间距得到的半方差函数有明显的结构，但研究区内样点分布不均匀，则可进一步优化其分布。 对于第（1）（2）情况，可以采取套合采样解决。 4.3 变异函数的理论拟合模型 理论变异函数用来拟合一些列经验变异函数值，供后续进行插值估计时使用。 选用理论变异函数模型是，要根据经验半方差图的性状来选取合适的模型 4.3 变异函数的理论拟合模型 变异函数在原点处的性状 4.3 变异函数的理论拟合模型 变异函数的理论模型： 有基台值模型 无基台值模型 有基台值模型—球状模型 有基台值模型—指数模型 有基台值模型—高斯模型 三种常用模型比较 有基台值模型—线性有基台值模型 有基台值模型—纯块金效应模型 无基台值模型——线性无基台值模型 无基台值模型——幂函数值模型 无基台值模型——对数值模型 套合模型 在实际中，有时区域化随机变量Z(x)的变化相当复杂，往往包含各种尺度及各种层次的变化，反映在变异函数r(h)上，就是单一的模型结构不能将其合理表达，而是多层次的结构相互叠加在一起，地统计学上称为套合。所谓套合结构，就是把分别出现在不同距离h上或不同方向上同时起作用的变异性组合起来，对全部有效的结构信息，作定量化的概括，以表示区域化变量的主要特征。 套合模型 土壤是一个不均与、具有高度空间异质性的复合体，它与土壤母质、气候、水文、地形和生物等因素有关，分析土壤空间变异的因素，可将其变异