1、第第5讲讲 空间统计分析空间统计分析授课人:王 杰Email: 地理信息系统科研方法课程本讲内容本讲内容探索性空间统计分析探索性空间统计分析地统计分析方法地统计分析方法 空间统计分析空间统计分析, ,即空间数据(即空间数据(spatial data)的统)的统计分析,是现代计量地理学中一个快速发展的方向计分析,是现代计量地理学中一个快速发展的方向和领域。和领域。空间统计分析,其核心就是认识与地理位置相关的空间统计分析,其核心就是认识与地理位置相关的数据间的数据间的空间依赖、空间关联或空间自相关空间依赖、空间关联或空间自相关,通过,通过空间位置建立数据间的统计关系。空间位置建立数据间的统计关系。
2、空间统计分析的任务,就是运用有关统计方法,建空间统计分析的任务,就是运用有关统计方法,建立立空间统计模型空间统计模型,从凌乱的数据中,从凌乱的数据中挖掘空间自相关挖掘空间自相关与空间变异规律与空间变异规律。 空间统计分析空间统计分析空间数据分析与传统统计分析主要有空间数据分析与传统统计分析主要有两大差异:两大差异:(1)空间数据间并非空间数据间并非独立独立,而是在维空间中具有某种空间相关,而是在维空间中具有某种空间相关性,且在不同的空间分辨率下呈现不同之相关程度;性,且在不同的空间分辨率下呈现不同之相关程度;(2)地球只有一个,大多数空间问题仅有一组(空间分布不规地球只有一个,大多数空间问题仅
3、有一组(空间分布不规则的)观测值,而则的)观测值,而无重复观测数据无重复观测数据。因此,空间现象的了。因此,空间现象的了解与描述是极为复杂的,而传统方法,尤其是建立在独立解与描述是极为复杂的,而传统方法,尤其是建立在独立样本上的统计方法,不适合分析空间数据。样本上的统计方法,不适合分析空间数据。空间统计空间统计 VS. 经典统计经典统计n 经典统计:独立性、随机性假设经典统计:独立性、随机性假设n 空间统计:自相关、依赖性、异质性空间统计:自相关、依赖性、异质性地理学第一定律(地理学第一定律(FLGFLG): :空间统计的基本思想:Waldo Tobler(born in 1930) rece
4、iving a plaque for his contributions to geography. On the event of his November 2000 birthday. Tobler, W. R. (1970). A computer movie simulating urban growth in the Detroit region. Economic Geography, 46(2): 234-240.FLGFLG的一般性的一般性: : 自然地理、人文地理、社会经济n 空间自空间自相关相关是普遍存在的是普遍存在的,否则地理分析便没有多大意义。 经典统计:独立 空间自相
5、关的存在,使得经典统计学所要求的样本独立性假设不满足。n 如果地理学从根本上值得研究,必然是因为地理现象在空间上的变化地理现象在空间上的变化不是随机不是随机的的。 经典统计:随机 可以借助空间统计更好地理解地理现象。可以借助空间统计更好地理解地理现象。 或许学习空间统计最重要的原因是我们不仅仅想知道问题或许学习空间统计最重要的原因是我们不仅仅想知道问题“怎么怎么样样”,更想知道,更想知道“哪里怎么样哪里怎么样” 空间统计学可以帮助我们准确地判断具体地理模空间统计学可以帮助我们准确地判断具体地理模式的原因。式的原因。 John Snow的霍乱地图的霍乱地图 当发现某种病仅仅发生在靠近河流的村庄时
6、,河流中的寄生物可当发现某种病仅仅发生在靠近河流的村庄时,河流中的寄生物可能是病源。能是病源。 空间统计学可以帮助我们处理大的复杂数据集空间统计学可以帮助我们处理大的复杂数据集, 这是这是GIS经常面对的事情。经常面对的事情。为什么要用空间统计为什么要用空间统计霍乱病死者居住分布图(John Snow, 1854)v 18541854年年8 8月到月到9 9月月英国伦敦霍乱流英国伦敦霍乱流行时,当局始终找不到发病的原行时,当局始终找不到发病的原因,后来医生约翰因,后来医生约翰斯诺斯诺(John (John Snow) Snow) 参与调查。参与调查。v 他在绘有霍乱流行地区所有道路、他在绘有霍
7、乱流行地区所有道路、房屋、饮用水机井等内容的房屋、饮用水机井等内容的1 1:65006500比例尺地图上,标出了每个比例尺地图上,标出了每个霍乱病死者的霍乱病死者的居住位置居住位置,得到了,得到了霍乱病死者居住分布图霍乱病死者居住分布图。一一. 探索性空间统计分析探索性空间统计分析 基本原理与方法 应用实例 探索性空间数据分析探索性空间数据分析(ESDA)ESDA是指利用统计学原理统计学原理和图形图表图形图表相结合对空间信息的性质进行分析、鉴别,用以引导确定性模型的结构和解法。 ESDA与EDA区别在于它考虑了数据的空间特性,它考虑了数据的空间特性,在方法上它将数据分解为一般趋势和叠加于其上的
8、在方法上它将数据分解为一般趋势和叠加于其上的局部变化两部分局部变化两部分。然后用一定的数学函数去拟合由样本点产生的经验变率函数,进行诸如克立格内插等空间操作。 通常定义一个二元对称空间权重矩阵W,来表达n个位置的空间区域的邻近关系,其形式如下式中:Wij表示区域i与j的临近关系,它可以根据邻接标准或距离标准来度量。 1. 基本原理与方法基本原理与方法 (一)空间权重矩阵 nnnnnnwwwwwwwwwW212222111211 简单的二进制邻接矩阵简单的二进制邻接矩阵 基于距离的二进制空间权重矩阵基于距离的二进制空间权重矩阵 n两种最常用的确定空间权重矩阵的规则 其他相邻接和当区域01jiwi
9、j其他时的距离小于和当区域01djiwij(二)全局空间自相关(二)全局空间自相关 Moran指数反映的是空间邻接或空间邻近的区域单元属性值的相似程度。 Geary 系数与Moran指数存在负相关关系。 Patrick A.P.Moran (1917-1988) Moran指数和Geary系数是两个用来度量空间自相关的全局指标。 全局空间自相关概括了在一个总的空间范围内空间依赖的程度。 如果是位置(区域)的观测值,则该变量的全局Moran指数I,用如下公式计算 式中: I 为Moran指数; ninjniiijninjjiijxxwxxxxwnI111211ninijijninijjiijwS
10、xxxxw121)(22)(1iixxnSniixnx11;。 Geary 系数C计算公式如下 式中:C为Geary系数;其他变量同上式。 如果引入记号 ninjniiijninjjiijxxwxxwnC111211221ninjijwS110)(xxzii)(xxzjj,21nTzzzz 则全局Moran指数I的计算公式也可以进一步写成 Moran指数指数I的取值一般在的取值一般在-1,1之间之间,小于小于0表示负相表示负相关,等于关,等于0表示不相关,大于表示不相关,大于0表示正相关;表示正相关; Geary系数系数C的取值一般在的取值一般在0,2之间,大于之间,大于1表示负表示负相关,等
11、于相关,等于1表示不相关,而小于表示不相关,而小于1表示正相关。表示正相关。 niininjjiijxxxxxxwSnI12110)()(zzWzzSnzzzwSnTTniininjjiij012110 (三)局部空间自相关(三)局部空间自相关 描述一个空间单元与其领域的相似程度,能够表示每个局部单元服从全局总趋势的程度(包括方向和量级),反映了空间异质性,说明空间依赖是如何随位置变化的。 局部空间自相关分析方法包括局部空间自相关分析方法包括3种种:p 空间联系的局部指标空间联系的局部指标(LISA); p G统计量统计量;p Moran散点图散点图n空间联系的局部指标空间联系的局部指标(LI
12、SA)(LISA) 空间联系的局部指标(空间联系的局部指标(local indicators of spatial association ,缩写为,缩写为LISA)满足下列两个条件:)满足下列两个条件: (1)每个区域单元的)每个区域单元的LISA,是,是描述该区域单元周围描述该区域单元周围显著的相似值区域单元之间空间集聚程度的指标显著的相似值区域单元之间空间集聚程度的指标; (2)所有区域单元)所有区域单元LISA的总和与全局的空间联系指的总和与全局的空间联系指标成比例。标成比例。 LISA包括包括局部局部Moran指数指数(local Moran index)和和局部局部Geary指数指
13、数(local Geary index),下面重),下面重点介绍和讨论局部点介绍和讨论局部Moran指数。指数。 局部局部Moran指数被定义为指数被定义为可进一步写成可进一步写成 式中:式中: 和和 是经过标准差标准化的观测值。是经过标准差标准化的观测值。 局部局部Moran指数检验的标准化统计量为指数检验的标准化统计量为 jjijiixxwSxxI)()(2iijjijiixxxxwxxnI2)()()(jjijiTjjijizwzzzzwnzizjz)()()(iiiiIVARIEIIZn G统计量统计量 v 全局全局G统计量的计算公式为统计量的计算公式为v 对每一个区域单元的统计量为对
14、每一个区域单元的统计量为 ijijjijiijxxxxwG/ijjjijixxwG/ 探测区域单元是属于高值集聚还是低值集聚的空探测区域单元是属于高值集聚还是低值集聚的空间分布模式间分布模式. . 对统计量的检验与局部对统计量的检验与局部Moran指数相似,其检验值为指数相似,其检验值为 显著的正值表示在该区域单元周围,高观测值的区域单显著的正值表示在该区域单元周围,高观测值的区域单元趋于空间集聚,而显著的负值表示低观测值的区域单元趋元趋于空间集聚,而显著的负值表示低观测值的区域单元趋于空间集聚于空间集聚, ,与与Moran指数只能发现相似值指数只能发现相似值( (正关联正关联) )或非相似或
15、非相似性观测值性观测值( (负关联负关联) )的空间集聚模式相比,的空间集聚模式相比,具有能够探测出区具有能够探测出区域单元属于高值集聚还是低值集聚的空间分布模式域单元属于高值集聚还是低值集聚的空间分布模式。)()()(iiiiGVARGEGGZn Moran散点图散点图 以(以(Wz,z)为坐标点的为坐标点的Moran散点图,常来研究局部的散点图,常来研究局部的空间不稳定性,它对空间滞后因子空间不稳定性,它对空间滞后因子Wz和和z数据对进行了可视数据对进行了可视化的二维图示。化的二维图示。 全局全局Moran指数,可以看作是指数,可以看作是Wz对于对于z的线性回归系数,的线性回归系数,对界外
16、值以及对对界外值以及对Moran指数具有强烈影响的区域单元,可通指数具有强烈影响的区域单元,可通过标准回归来诊断出。过标准回归来诊断出。 由于数据对(由于数据对(Wz,z)经过了标准化,因此界外值可易由)经过了标准化,因此界外值可易由2sigma规则可视化地识别出来。规则可视化地识别出来。 Moran散点图的散点图的4 4个象限个象限,分别对应于区域单元与其邻居分别对应于区域单元与其邻居之间之间4 4种类型的局部空间联系种类型的局部空间联系形式:形式:l 第第1 1象限代表了象限代表了高观测值高观测值的的区域单元被区域单元被高值高值的区域所包围的区域所包围的空间联系形式;的空间联系形式;l 第
17、第2 2象限代表了象限代表了低观测值低观测值的的区域单元被区域单元被高值高值的区域所包围的区域所包围的空间联系形式;的空间联系形式; l 第第3 3象限代表了象限代表了低观测值低观测值的区域单元被的区域单元被低值低值的区域所的区域所包围的空间联系形式;包围的空间联系形式;l 第第4 4象限代表了象限代表了高观测值高观测值的区域单元被的区域单元被低值低值的区域所的区域所包围的空间联系形式。包围的空间联系形式。2. 应用实例应用实例 中国大陆30个省级行政区人均GDP的空间关联分析。根据各省(直辖市、自治区)之间的邻接关系,采用二进制邻接权重矩阵,选取各省(直辖市、自治区)19982002年人均G
18、DP的自然对数,依照公式计算全局Moran指数I,计算其检验的标准化统计量Z(I),结果如下表所示。年份年份IZP19980.50014.503 50.000 019990.506 94.555 10.000 020000.511 24.597 80.000 020010.505 94.553 20.000 020020.501 34.532 60.000 0 从表从表中中可以看出,在可以看出,在19982002年期间,中国大陆年期间,中国大陆30个个省级行政区人均省级行政区人均GDP的全局的全局Moran指数均为正值;在正态指数均为正值;在正态分布假设之上,对分布假设之上,对Moran指数检
19、验的结果也高度显著。这就指数检验的结果也高度显著。这就是说,在是说,在19982002年期间,中国大陆年期间,中国大陆30个省级行政区人均个省级行政区人均GDP存在着显著的、正的空间自相关,也就是说各省存在着显著的、正的空间自相关,也就是说各省级行级行政区政区人均人均GDP水平的空间分布表现出相似值之间的空间集水平的空间分布表现出相似值之间的空间集聚,其空间联系的特征是:较高人均聚,其空间联系的特征是:较高人均GDP水平的省级行政水平的省级行政区相对地趋于和较高人均区相对地趋于和较高人均GDP水平的省水平的省级行政区级行政区相邻,或相邻,或者较低人均者较低人均GDP水平的水平的省级行政区省级行
20、政区相对地趋于和较低人均相对地趋于和较低人均GDP水平的水平的省级行政区省级行政区相邻。相邻。 选取2001年我国30个省级行政区人均GDP数据,计算局部Gi统计量和局部Gi统计量的检验值Z(Gi),并绘制统计地图如下。 检验结果表明,贵州、四川、云南西部3省的Z值在0.05的显著性水平下显著,重庆的Z值在0.1的显著性水平下显著,该4省市在空间上相连成片分布,而且从统计学意义上来说,与该区域相邻的省区,其人均GDP趋于为同样是人均GDP低值的省区所包围。由此形成人均GDP低值与低值的空间集聚,据此可认识到西部落后省区趋于空间集聚的分布特征。 东部的江苏、上海、浙江三省市的Z值在0.05的显著
21、性水平下显著,天津的Z值在0.1的显著性水平下显著。而东部上海、江浙等发达省市趋于为一些相邻经济发展水平相对较高的省份所包围,东部发达地区的空间集聚分布特征也显现出来。 以(Wz,z)为坐标,进一步绘制Moran散点图 可以发现,多数省(直辖市、自治区)位于第1和第3象限内,为正的空间联系,属于低低集聚和高高集聚类型,而且位于第3象限内的低低集聚类型的省(直辖市、自治区)比位于第1象限内的高高集聚类型的省(直辖市、自治区)更多一些。 上图进一步显示了30个省级行政区人均GDP局部集聚的空间结构。可以看出,从人均GDP水平相对地来看: 高值被高值包围的高高集聚省(直辖市)有:北京、天津、河南、安
22、徽、湖北、江西、海南、广东、福建、浙江、山东、上海、江苏; 低值被低值包围的低低集聚省(自治区)有:黑龙江、内蒙古、新疆、吉林、甘肃、山西、陕西、青海、西藏、四川、云南、辽宁、贵州; 被低值包围的高值省(直辖市)有:重庆、广西、河北;被高值包围的低值省份只有湖南。1978、1990、2001、2007年全国各省市人均GDP的LISA显著水平图二、地统计分析方法二、地统计分析方法地统计方法的基本原理 应用实例 地统计学地统计学(Geostatistics), 又称地质统计学,是法国著名统计学又称地质统计学,是法国著名统计学家家G.Matheron在大量理论研究基础上提出的。在大量理论研究基础上提
23、出的。 地统计学是以地统计学是以区域化变量理论区域化变量理论为基础,以变异函数为主要工为基础,以变异函数为主要工具,研究那些在空间分布上既有随机性又有结构性,或空间具,研究那些在空间分布上既有随机性又有结构性,或空间相关和依赖性,或空间格局与变异,并对这些数据进行最优相关和依赖性,或空间格局与变异,并对这些数据进行最优无偏内插估计,或模拟这些数据的离散性、波动性。无偏内插估计,或模拟这些数据的离散性、波动性。 协方差函数协方差函数和和变异函数变异函数是以区域化变量理论为基础建立起来是以区域化变量理论为基础建立起来的地统计学的两个最基本的函数。地统计学的主要方法之一,的地统计学的两个最基本的函数
24、。地统计学的主要方法之一,克立格法就是建立在变异函数理论和结构分析基础之上的。克立格法就是建立在变异函数理论和结构分析基础之上的。 当一个变量当一个变量呈现为空间分布呈现为空间分布时,就称之为区域化变量。这种变时,就称之为区域化变量。这种变量常常反映某种空间现象的特征,用区域化变量来描述的现象称之量常常反映某种空间现象的特征,用区域化变量来描述的现象称之为为区域化现象区域化现象。 区域化变量区域化变量,亦称,亦称区域化随机变量区域化随机变量,G. Matheron(1963)将)将它定义为以空间点它定义为以空间点x的三个直角坐标为自变量的随机场的三个直角坐标为自变量的随机场 区域化变量具有区域
25、化变量具有两个最显著两个最显著,而且也是最重要的特征,即,而且也是最重要的特征,即随机随机性和结构性性和结构性。1. 地统计方法的基本原理地统计方法的基本原理 (一)区域化变量 ),(wvxxxZZux区域化变量的功能:区域化变量的功能:由于区域化变量是一种随机函数,因而能同时反映空间变量的结构性和随机性。l一方面,当空间点 x 固定后,Z(x)就是一个随机变量,这体现了其随机性随机性。l另一方面,在空间两个不同点 x 与 x+h 处的区域化变量值具有某种程度的相关性,这体现了其结构性结构性。区域化变量的组成部分区域化变量的组成部分l 数据点 结构性 可以用均值和常数趋势表示 空间相关 数据通
26、常呈现正空间相关性 随机性 测量误差,其他误差 distance elevation结构性结构性随机性随机性实际值实际值(二)协方差函数(二)协方差函数 协方差函数的概念协方差函数的概念 区域化随机变量之间的差异,可以用空间协方差来表示。区域化随机变量之间的差异,可以用空间协方差来表示。 在概率论中在概率论中,随机向量随机向量X与与Y的协方差被定义为的协方差被定义为 区域化变量区域化变量 在空间点在空间点x和和x+h处的两个随机处的两个随机变量和的二阶混合中心矩定义为变量和的二阶混合中心矩定义为Z(x)的自协方差函数,即的自协方差函数,即),()(wvuxxxZxZ)()()()()(),(h
27、xZExZEhxZxZEhxZxZCov(4.2.24.2.2))(),(EyyExxEyxcov(4.2.1)(4.2.1)n协方差函数的计算公式协方差函数的计算公式 式中:式中:h h为两样本点空间分隔距离或距离滞后;为两样本点空间分隔距离或距离滞后; 为为 在空间位置在空间位置 处的实测值;处的实测值; 是是 在在 处距离偏离处距离偏离h的实测值的实测值 i=1=1,2 2, , 是分隔距离为是分隔距离为h时的样本点对(时的样本点对(parisparis)总数,)总数, 和和 分别为分别为 和和 的样本平均数的样本平均数, ,即即 ix)(1)()()()()(1)(hNiiiiihxZ
28、hxZxZxZhNhc)(xZ)(hxZi)(xZ)(hN)(ixZ)(hxZi)(ixZ)(hxZi(4.2.3)(4.2.3)NiiixZNxZ1)(1)()(ixZNiiihxZNhxZ1)(1)(4.2.4)(4.2.4)(4.2.5)(4.2.5) 若若 = =m(常数),则上式可以改写为(常数),则上式可以改写为 式中:式中:m为样本平均数,可由一般算术平均数公式求得,即为样本平均数,可由一般算术平均数公式求得,即 )(ixZ)(hxZi)(12)()()(1)(hNiiimhxZxZhNhcniixZNm1)(1(4.2.6)(4.2.6)(三)变异函数(三)变异函数 变异函数的
29、概念变异函数的概念 变异函数变异函数(variograms),又称变差函数、变异矩,是地统计,又称变差函数、变异矩,是地统计分析所特有的基本工具。分析所特有的基本工具。 在一维条件下变异函数定义为,在一维条件下变异函数定义为,当空间点当空间点x在一维在一维x轴上变轴上变化时,区域化变量化时,区域化变量Z(x)在点在点x和和x+h处的值处的值Z(x)与与Z(x+h)差的方差的方差的一半为区域化变量差的一半为区域化变量Z(x)在在x轴方向上的变异函数轴方向上的变异函数,记为,记为(h),即即 )()(21),(hxZxZVarhx22)()(21)()(21hxZExZEhxZxZE(4.2.7)
30、(4.2.7) 在在二阶平稳假设二阶平稳假设条件下,对任意的条件下,对任意的h有有 因此,公式可以改写为因此,公式可以改写为 从上式可知,变异函数依赖于两个自变量从上式可知,变异函数依赖于两个自变量 x和和h,当变异函,当变异函数数 仅仅依赖于距离仅仅依赖于距离h h而与位置而与位置x无关时,无关时, 可改写可改写成成 ,即,即 )()(xZEhxZE2)()(21),(hxZxZEhx),(hx)(h),(hx2)()(21)(hxZxZEh(4.2.9)(4.2.9)(4.2.8)(4.2.8)变异函数的性质变异函数的性质 设设Z(x)是区域化变量,在满足二阶平稳假设条件下,变异函是区域化
31、变量,在满足二阶平稳假设条件下,变异函数式具有如下性质:数式具有如下性质: (1)(1) =0 =0,即在,即在h=0=0处,变异函数为处,变异函数为0 0; (2(2) = ) = ,即,即 关于直线关于直线h=0=0是对称的,它是一个偶函是对称的,它是一个偶函数;数; (3) 0(3) 0,即,即 只能大于或等于只能大于或等于0 0;)0()( h)(h)(h)( h)(hn变异函数的计算公式 v 设设 是系统某属性是系统某属性Z在空间位置在空间位置x处的值,处的值, 为一区为一区域化随机变量,并满足二阶平稳假设,域化随机变量,并满足二阶平稳假设,h为两样本点空间为两样本点空间分隔距离,分
32、隔距离,v 和和 分别是区域化变量分别是区域化变量 在空间位置在空间位置 和和 处的实测值处的实测值i=1,2,N(h),那么,变异函数,那么,变异函数 的离散计的离散计算公式为算公式为)(xZ)(xZ)(ixZ)(hxZiix)(xZhxi)(h)(12)()()(21)(hNiiihxZxZhNh(4.2.10)(4.2.10) 这样对不同的空间分隔距离这样对不同的空间分隔距离h,计算出相应的,计算出相应的 和和 值。值。如果分别以如果分别以h为横坐标,为横坐标, 或或 为纵坐标,画出协方差函数和为纵坐标,画出协方差函数和变异函数曲线图,就可以直接展示区域化变量变异函数曲线图,就可以直接展
33、示区域化变量Z(x)的空间变异特的空间变异特点。可见,变异函数能同时描述区域化变量的随机性和结构性,点。可见,变异函数能同时描述区域化变量的随机性和结构性,从而在数学上对区域化变量进行严格分析,是空间变异规律分析从而在数学上对区域化变量进行严格分析,是空间变异规律分析和空间结构分析的有效工具。和空间结构分析的有效工具。)(hc)(h)(hc)(h例如:假设某地区降水量Z(x)(单位:mm)是二维区域化随机变量,满足二阶平稳假设,其观测值的空间正方形网格数据如图4.2.1所示(点与点之间的距离为h=1 km)。试计算其南北方向及西北和东南方向的变异函数。图图4.2.1 4.2.1 空间正方形网格
34、数据(点间距空间正方形网格数据(点间距h h=1 km=1 km) 从图4.2.1可以看出,空间上有些点,由于某种原因没有采集到。如果没有缺失值,可直接对正方形网格数据结构计算变异函数;在有缺失值的情况下,也可以计算变异函数。只要“跳过”缺失点位置即可(图4.2.2)。图图4.2.2 4.2.2 缺失值情况下样本数对的组成和计算过程缺失值情况下样本数对的组成和计算过程为缺失值 图图1 空间正方形网格数据(点间距空间正方形网格数据(点间距h=1km) 图图2 空间正方形网格数据(点间距空间正方形网格数据(点间距h=2km) 首先计算南北方向上的变异函数值,由变异函数的计算公式可得 =385/72
35、=5.35 图图4.2.2 4.2.2 缺失值情况下样本数对的组成和计算过程缺失值情况下样本数对的组成和计算过程为缺失值 222222)3837()3836()3635()3537()3742()4240(3621) 1 (222222)4242()3536()3743()4340()3735()3538(222222)3738()3839()3940(3535)3535()3542()(222222)3336()3637()3739()3939()3034()3437(222222)4036()2932()3236()3637()3741()4137()3228()3438()3029()2
36、935()3533()3340(222222 同样计算出 最后,得到南北方向和西北东南方向上的变异函数计算结果见下表。同样可以计算东西方向上的变异函数。 9.26)2(17.55)3(25.69)4(22.90)5( 方向 南北 方向 西北东南 h12345h1.412.824.245.657.07N(h) 36 27 21 13 5 N(h) 322113825.359.2617.5525.6922.907.0612.9530.8558.1350.00)(h)(hn变异函数的参数 变异函数有变异函数有4个非常重要的参数,即个非常重要的参数,即基台值(基台值(sill)、变变程(程(range
37、)或称空间依赖范围(或称空间依赖范围(range of spatial dependence)、)、块金值(块金值(nugget)或称区域不连续性值或称区域不连续性值(localized discontinuity)和)和分维数(分维数(fractal dimension)。 前前3个参数可以直接从变异函数图中得到。它们决定变个参数可以直接从变异函数图中得到。它们决定变异函数的形状与结构。异函数的形状与结构。 变异函数的形状反映自然现象空间分布结构或空间相变异函数的形状反映自然现象空间分布结构或空间相关的类型,同时还能给出这种空间相关的范围。关的类型,同时还能给出这种空间相关的范围。n 当变异
38、函数随着间隔距离当变异函数随着间隔距离h的增大,从非零值达到一个相对的增大,从非零值达到一个相对稳定的常数时,该常数称为稳定的常数时,该常数称为基台值基台值C0+C。n 当间隔距离当间隔距离h=0时,时,(0)= C0,该值称为,该值称为块金值块金值或块金方差或块金方差(n nugget variance)。)。n 基台值是系统或系统属性中最大的变异,变异函数达到基台基台值是系统或系统属性中最大的变异,变异函数达到基台值时的间隔距离值时的间隔距离a a称为称为变程变程。变程表示在。变程表示在ha以后,区域化变以后,区域化变量量Z Z( (x) )空间相关性消失。空间相关性消失。 第第4个参数,
39、即个参数,即分维数用于表示变异函数的特性分维数用于表示变异函数的特性,由变异函数由变异函数 和间隔距离和间隔距离h之间的关系确定之间的关系确定 其中,分维数其中,分维数D为双对数直线回归方程中的斜率,它是一个无量纲为双对数直线回归方程中的斜率,它是一个无量纲数。分维数数。分维数D的大小,表示变异函数曲线的曲率,可以作为随机变异的大小,表示变异函数曲线的曲率,可以作为随机变异的量度。的量度。)24()(2Dhh)(h理论变异函数模型理论变异函数模型实践中,常用的是变异函数图:偏基台值:C(partial sill)块金值:C0(nugget)变程: a(range)h( )h基台值(sill)n
40、ot related anymore变程范围内才有结构性变程范围内才有结构性变化(有规律的变化)变化(有规律的变化)反映随机性大小:反映随机性大小:主要来源于区域化变量主要来源于区域化变量Z(x)在小于抽样尺度在小于抽样尺度h时所具有的内时所具有的内部变异;另外还有抽样分析误差。部变异;另外还有抽样分析误差。变异函数是一个单变异函数是一个单调不减函数。当调不减函数。当h超过某一个范围,超过某一个范围,例如变程,变异函例如变程,变异函数不再增大,而是数不再增大,而是趋于一个极限值,趋于一个极限值,即为基台值。实际即为基台值。实际上等于区域化变量上等于区域化变量的先验方差。即,的先验方差。即,即基
41、台值与块金值之即基台值与块金值之差,表示数据中存在差,表示数据中存在空间相关性引起的方空间相关性引起的方差变化范围。差变化范围。( ) ( )(0)Var Z xC 00lim ( )hChn变异函数的理论模型 地统计学将变异函数理论模型分为地统计学将变异函数理论模型分为3 3大类:大类: 第第1 1类是类是有基台值模型有基台值模型,包括球状模型、指数模型、高,包括球状模型、指数模型、高斯模型、线性有基台值模型和纯块金效应模型;斯模型、线性有基台值模型和纯块金效应模型; 第第2 2类是类是无基台值模型无基台值模型,包括幂函数模型、线性无基台,包括幂函数模型、线性无基台值模型、抛物线模型;值模型
42、、抛物线模型; 第第3 3类是类是孔穴效应模型孔穴效应模型。 下面有代表性地介绍几种常见的变异函数理论模型。下面有代表性地介绍几种常见的变异函数理论模型。 纯块金效应模型纯块金效应模型: :其一般公式为其一般公式为 式中:式中:c0 000,为先验方差。该模型相当于区域化变量,为先验方差。该模型相当于区域化变量为随机分布,样本点间的协方差函数对于所有距离为随机分布,样本点间的协方差函数对于所有距离h均均等于等于0 0,变量的空间相关不存在。,变量的空间相关不存在。 000)(0hchh(4.2.11)(4.2.11) 球状模型球状模型: :其一般公式为其一般公式为 式中:式中:c0为块金(效应
43、)常数为块金(效应)常数; ;c为拱高为拱高; ;c0+c为基为基台值台值; ;a为变程。当为变程。当c0=0,c=1时,称为标准球状模型。时,称为标准球状模型。球状模型是地统计分析中应用最广泛的理论模型,球状模型是地统计分析中应用最广泛的理论模型,许多区域化变量的理论模型都可以用该模型去拟合。许多区域化变量的理论模型都可以用该模型去拟合。 ahccahahahcchh03300)223(00)(4.2.12)(4.2.12) 指数模型指数模型: :其一般公式为其一般公式为 式中:式中:c c0 0和和c c意义与前相同,但意义与前相同,但a a不是变程。当不是变程。当h h=3=3时,时,
44、,即,即 ,从而指数模型的变程,从而指数模型的变程 约为约为 约为约为3a。当。当c c0 0=0=0,c c=1=1时,称为标准指数模型。时,称为标准指数模型。0)1 (00)(0hecchhah195. 0113eeahcca0)3(a(4.2.13)(4.2.13) 高斯模型高斯模型:其一般公式为其一般公式为 式中:式中:c0 0和和c意义与前相同,意义与前相同,a也不是变程。当也不是变程。当 时,时, ,即,即 ,因此高斯模型的,因此高斯模型的变程变程 约为约为 。当。当 时,称为标准高斯函数时,称为标准高斯函数模型。模型。0)1 (00)(220hecchhahah3195. 011
45、322eeahcca0)3(1, 00ccaa3(4.2.14)(4.2.14)幂函数模型幂函数模型: :其一般公式为其一般公式为 式中:式中:为幂指数。当为幂指数。当变化时,这种模型可以反变化时,这种模型可以反映在原点附近的各种性状。但是映在原点附近的各种性状。但是必须小于必须小于2,若,若 ,则函数则函数 就不再是一个条件非负定函数了,也就是就不再是一个条件非负定函数了,也就是说它已经不能成为变异函数了。说它已经不能成为变异函数了。 20 ,)(Ahh2(-h)(4.2.15)(4.2.15) 对数模型对数模型: :其一般公式为其一般公式为 显然,当显然,当 ,这与变异函,这与变异函数的性
46、质数的性质 不符不符 。因此,对数模型不能。因此,对数模型不能描述点支撑上的区域化变量的结构。描述点支撑上的区域化变量的结构。hAhlg)(hhlog, 00)(h(4.2.16)(4.2.16) 线性有基台值模型线性有基台值模型: :其一般公式为其一般公式为 式中式中:该模型的变程为该模型的变程为a,基台值为,基台值为 。 线性无基台值模型线性无基台值模型: :其一般公式为其一般公式为 从式中可以看出,该模型没有基台值,也没有从式中可以看出,该模型没有基台值,也没有变程。变程。 ahcahAhhh00c00c)(cc 000)(0hAhhch(4.2.18)(4.2.18)(4.2.17)(
47、4.2.17)例如:某地区降水量是一个区域化变量,其变异函数 的实测值及距离h的关系见下表,下面我们试用回归分析方法建立其球状变异函数模型。)(h实测值(h)距离h实测值(h)距离h2.10.69.24.94.31.110.35.15.72.210.56.26.52.510.97.57.83.111.29.58.83.812.49.8 从上面的介绍和讨论,我们知道,球状变异函数的一般形式为 当 时,有ahccahahahcchh03300)223(00)(ah 0330)2()23()(hachacch 如果记 ,则可以得到线性模型 根据表中的数据,对上式进行最小二乘拟合,得到 (4.2.20
48、) 计算可知,上式的显著性检验参数F=114.054,R2=0.962,可见模型的拟合效果是很好的。32132100,21,23,),(hxhxacbacbcbhy22110 xbxbby2192007. 0731. 1048. 2xxy(4.2.19) 比较(4.2.20)式与(4.2.19)式,并做简单计算可知:c0=2.048,c=1.154,a=8.353,所以,球状变异函数模型为535. 8202. 3535. 80)535. 821535. 823(154. 1048. 200)(33*hhhhhh(4.2.21) ( (四四) )克立格插值方法克立格插值方法 克立格(克立格(Kr
49、iging)插值法)插值法,又称空间局部估计或空间局部插值法。它建立在变异函数理论及结构分析变异函数理论及结构分析基础之上的,在有限区域内对区域化变量的取值进行无偏最优估计的一种方法。 克立格法适用的条件是,如果变异函数和相关分析的结果表明区域化变量存区域化变量存在空间相关性在空间相关性。 其实质是利用区域化变量的原始数据和变异函数的结构特点,对未采样点的区域化变量的取值进行线性无偏、最优估计线性无偏、最优估计。具体来说,它是根据待估样点(或待估块段)有限邻域内若干已测定的样点数据,在认真考虑样点的形状、大小和空间相互位置关系,它们与待估样点相互空间位置关系,以及变异函数提供的结构信息之后,对
50、该待估样点值进行的一种线性无偏最优估计。 克立格插值(克立格插值(kriging interpolation)是根据变异函数模型而是根据变异函数模型而发展起来的一系列地统计的空间插值方法,包括:发展起来的一系列地统计的空间插值方法,包括:p 普通克立格法(普通克立格法(ordinary kriging);p 泛克立格法(泛克立格法(universal kriging);p 指示克立格法(指示克立格法(indicator kriging);p 析取克立格法(析取克立格法(disjunctive kriging);p 协同克立格法(协同克立格法(cokriging)等。)等。 下面仅对普通克立格法