1、12第十一讲第十一讲 空间插值空间插值地统计学基础地统计学基础3主要内容：主要内容：1、引言、引言2、空间插值基础、空间插值基础3、kriging空间插值空间插值4一、地统计学发展概述一、地统计学发展概述引言引言5地统计学是以具有空间分布特点的区域化变量理论为基础，研究自然现象地统计学是以具有空间分布特点的区域化变量理论为基础，研究自然现象的空间变异与空间结构的一门学科。的空间变异与空间结构的一门学科。它针对像矿产、资源、生物群落、地貌等有着特定的地域分布特征而发展它针对像矿产、资源、生物群落、地貌等有着特定的地域分布特征而发展的统计学。由于最先在地学领应用，故称为地统计学。地统计学的主要理的
2、统计学。由于最先在地学领应用，故称为地统计学。地统计学的主要理论是法国统计学家论是法国统计学家G.Matheron 1963年创立的，经过不断完善和改进，目年创立的，经过不断完善和改进，目前已成为具有坚实理论基础和实用价值的数学工具。前已成为具有坚实理论基础和实用价值的数学工具。地统计学的应用范围十分广泛，不仅可以研究空间分布数据的结构性和随地统计学的应用范围十分广泛，不仅可以研究空间分布数据的结构性和随机性、空间相关性和依赖性、空间格局与变异，还可以对空间数据进行最机性、空间相关性和依赖性、空间格局与变异，还可以对空间数据进行最优无偏内插，以及模拟空间数据的离散性及波动性。优无偏内插，以及模
3、拟空间数据的离散性及波动性。地统计学由分析空间变异与结构的变异函数及其参数和空间局部估计的地统计学由分析空间变异与结构的变异函数及其参数和空间局部估计的Kriging插值法两个主要部分组成，目前已在地球物理、地质、生态、土插值法两个主要部分组成，目前已在地球物理、地质、生态、土壤等领域应用。气象领域的应用目前还不多见，主要使用壤等领域应用。气象领域的应用目前还不多见，主要使用Kriging法进行法进行降水、温度等要素的最优内插的研究及气候对农业影响方面的研究。降水、温度等要素的最优内插的研究及气候对农业影响方面的研究。引言引言引言引言6引言引言引言引言 地统计学是以区域化变量理论为基础，以变异
4、地统计学是以区域化变量理论为基础，以变异函数为主要工具，研究那些在空间分布上既有随机函数为主要工具，研究那些在空间分布上既有随机性又有结构性，或空间相关和依赖性的自然现象的性又有结构性，或空间相关和依赖性的自然现象的科学。科学。协方差函数和变异函数是以区域化变量理论为协方差函数和变异函数是以区域化变量理论为基础建立起来的地统计学的两个最基本的函数。地基础建立起来的地统计学的两个最基本的函数。地统计学的主要方法之一，克立格法就是建立在变异统计学的主要方法之一，克立格法就是建立在变异函数理论和结构分析基础之上的。函数理论和结构分析基础之上的。7引言引言引言引言 当一个变量呈现为空间分布时，就称之为
5、区域化变量当一个变量呈现为空间分布时，就称之为区域化变量（regionalized variable）。这种变量常常反映某种空间现象的特征，）。这种变量常常反映某种空间现象的特征，用区域化变量来描述的现象称之为区域化现象。用区域化变量来描述的现象称之为区域化现象。区域化变量，亦称区域化随机变量，区域化变量，亦称区域化随机变量，G.Matheron（1963）将它定）将它定义为以空间点义为以空间点x的三个直角坐标为自变量的随机场的三个直角坐标为自变量的随机场 。区域化变量具有两个最显著，而且也是最重要的特征，即随机性和区域化变量具有两个最显著，而且也是最重要的特征，即随机性和结构性。结构性。区域
6、化变量区域化变量),(wvxxxZZux8引言引言引言引言9引言引言引言引言10超出超出标准的标准的重金重金属属土壤土壤分布范围分布范围11引言引言制图和地理可视化制图和地理可视化 区域化变量分析过程区域化变量分析过程12引言引言经过链接和级联更新的多变量探究性数据分析经过链接和级联更新的多变量探究性数据分析 13引言引言空间自相关分析空间自相关分析-LISA聚集图和重要性图 14151617小结小结地统计学是以具有空间分布特点的区域化变量理论为基础，研究自然现地统计学是以具有空间分布特点的区域化变量理论为基础，研究自然现象的空间变异与空间结构的一门学科。象的空间变异与空间结构的一门学科。理论
7、核心理论核心-区域化变量、变异函数区域化变量、变异函数地统计学的主要用途，是研究对象空间自相关结构地统计学的主要用途，是研究对象空间自相关结构(或空间变异结构或空间变异结构)的探的探测以及变量值的估计和模拟。测以及变量值的估计和模拟。不管哪一种用途，地统计学分析的核心是根据样本点来确定研究对象不管哪一种用途，地统计学分析的核心是根据样本点来确定研究对象(某某一变量一变量)随空间位置而变化的规律，以此去推算未知点的属性值。这个规随空间位置而变化的规律，以此去推算未知点的属性值。这个规律，可用变异函数来表达。律，可用变异函数来表达。应用范围在不断扩大应用范围在不断扩大18地统计学基础地统计学基础1
8、9地统计学简介地统计学简介20地统计学简介地统计学简介21地统计学简介地统计学简介22地统计学简介地统计学简介23地统计学简介地统计学简介24地统计学简介地统计学简介25地统计学简介地统计学简介26地统计学简介地统计学简介27地统计学简介地统计学简介28地统计学简介地统计学简介29地统计学简介地统计学简介变异曲线变异曲线30地统计学简介地统计学简介31地统计学简介地统计学简介32地统计学简介地统计学简介33地统计学简介地统计学简介34地统计学简介地统计学简介35地统计学简介地统计学简介36地统计学简介地统计学简介37地统计学简介地统计学简介l空间统计分析，即地统计学(geostatistics
9、)，亦称地质统计学，于20世纪5O年代初开始形成，60年代在法国著名统计学家Matheron的大量理论研究工作基础上形成一门新的统计学分支(Matheron，1963)。地统计学是以区域化变量理论(theory of regionalized variable)为基础，以变异函数(variogram)为基本工具来研究分布于空间并呈现出一定的随机性和结构性的自然现象的科学，38l它的应用已被扩展到分析各种自然现象的空间异质性和空间格局。而森林受物理环境、自然干扰、人为干扰、树种特性、树种对干扰的反应等因素及其相互作用的影响，它具有高度空间异质性和时间异质性。地统计学简介地统计学简介39地统计学核
10、心地统计学核心l不管哪一种用途，地统计学分析的核心是根据样本点来确定研究对象(某一变量)随空间位置而变化的规律，以此去推算未知点的属性值。运用地统计学进行空间分析基本包括以下几个步骤，即数据探索性分析，空间连续性的量化模型，未知点属性值的估计，对未知点局部及空间整体不确定性的预测。40地质统计学主要软件工具地质统计学主要软件工具地质统计学大型软件地质统计学大型软件S-GeMS(开放代码软件开放代码软件)The Stanford Geostatistical Modeling Software(S-GeMS)大型3D地统计学建模软件，包含绝大多数经典地统计学算法。开放源代码，开放插件再开发标准。
11、作者：Dr.Nicolas RemySgems官方网站：http:/ekofisk.stanford.edu/SCRFweb/sgems/程序下载源代码:http:/ekofisk.stanford.edu/SCRFweb/sgems/.gems-1.2.tar.gz安装程序+源代码:http:/ekofisk.stanford.edu/SCRFweb/sgems/.staller-1.2.exe软件文档:用户手册:http:/ekofisk.stanford.edu/SCRFweb/sgems/.gems_manual.pdf程序文档:http:/ekofisk.stanford.edu/S
12、CRFweb/sgems/.html/index.html插件标准:http:/ekofisk.stanford.edu/SCRFweb/sgems/.ugin_howto.html SGeMS 发布最新版本发布最新版本2.0，新增加的特性包括：，新增加的特性包括：1）新增请多地质统计算法）新增请多地质统计算法2）改进多点地质统计算法）改进多点地质统计算法3）volume rendering to hide selected grid cells(?)4）新的网格文件格式）新的网格文件格式41地质统计学主要软件工具地质统计学主要软件工具Goldensoft Surfer 9.1 42地质统计学
13、主要软件工具地质统计学主要软件工具Goldensoft Surfer 9.1 滤滤波波功功能能43地质统计学主要软件工具地质统计学主要软件工具Goldensoft Surfer 9.1 滤滤波波功功能能44地质统计学主要软件工具地质统计学主要软件工具Goldensoft Surfer 9.1 45地质统计学主要软件工具地质统计学主要软件工具Golden Software Surfer(简称简称Surfer)是一款画三维图（等高线是一款画三维图（等高线,image map,3d surface）的软件，该软件简单易学，且其自带的英文帮助）的软件，该软件简单易学，且其自带的英文帮助（help菜单）
14、对如何使用菜单）对如何使用surfer解释的很详细，其中的解释的很详细，其中的tutorial教程更教程更是清晰的介绍了是清晰的介绍了surfer的简单应用，应该说的简单应用，应该说surfer软件自带的帮助文件软件自带的帮助文件是相当完美且容易阅读的是相当完美且容易阅读的,只要学过英语的人都可以很快上手。只要学过英语的人都可以很快上手。Surfer是具有插值功能的绘图软件，因此，即使你的数据是不等间距的，是具有插值功能的绘图软件，因此，即使你的数据是不等间距的，依然可以用它作图。但依据作者的经验，最好不使用依然可以用它作图。但依据作者的经验，最好不使用Surfer自带的插值自带的插值功能，尤
15、其是要精确确定等高线时。由于功能，尤其是要精确确定等高线时。由于surfer是美国的一个软件，它是美国的一个软件，它不提供对中文的支持，这可以算的上一个小的遗憾。不提供对中文的支持，这可以算的上一个小的遗憾。Surfer的主要功能是绘制等值线图（的主要功能是绘制等值线图（contour map），此外它还可以绘），此外它还可以绘制制post map，classed post map，vector map，image map，wireframe map，3d surface map，等形式的图形。其功能是比较强，等形式的图形。其功能是比较强的，但没有各种投影变化是它的一大缺点。尤其是在等值线领域
16、，这不的，但没有各种投影变化是它的一大缺点。尤其是在等值线领域，这不能不说是它的应用受到限制的地方。能不说是它的应用受到限制的地方。46主要工具主要工具-GeoDA47主要工具主要工具-GS+4849数据转换505152半方差函数分析53半方差函数模型l球状模型，指数模型-聚集分布l线性模型-随机分布54图表设置55转换成excel图表56查看GS+计算值列表57复制到excel中58在excel中制图59分维数60Moran s index 61Kriging 准备62Kriging 632D Map-等值线图643D Map65Arc/info Geostatistic模块66Arc/in
17、fo Geostatistic模块l将excel表格存储为数据库文件.xls .dbf 方法：excel：文件菜单 另存为 文件格式选择：dbase4（3）l 在arcview中打开数据库文件，将其转化为shapefile文件 .shp格式 打开tables，使用add命令添加数据库文件，选择view为当前操作界面，点击view菜单中的add event theme命令，将数据添加主题，选中后点击theme菜单中的convert to shapefile命令，转化为.shp文件67l将.shp文件转化为coverage 打开arc/info中的arctoolbox模块，选择import to
18、coverage下的shapefile to coverage命令，将.shp文件转化为coverage68l使用geostatistic模块分析 打开arcmap，使用add命令添加coverage，选中后选择geostatistic analysis命令下各选项 explore data semivariogram/conviarance 变异函数的分析 geostatistic wizard kriging制图 直方图69l区域化变量：小结小结70l变异函数与变异曲线小结小结71小结小结平稳假设与内蕴假设平稳假设与内蕴假设72小结小结平稳假设与内蕴假设平稳假设与内蕴假设73小结小结平稳假设与内蕴假设平稳假设与内蕴假设74
