当前位置:首页 >> 信息与通信 >>

arcgis空间内插值教程_图文

GIS 空间插值 (局部插值方法) 实习记录
一、空间插值的概念和原理
当我们需要做一幅某个区域的专题地图, 或是对该区域进行详细研究的时候, 必须具备 研究区任一点的属性值,也就是连续的属性值。但是,由于各种属性数据(如降水量、气温 等)很难实施地面无缝观测,所以,我们能获取的往往是离散的属性数据。例如本例,我们 现有一幅山东省等降雨量图,但是最终目标是得到山东省降水量专题图(覆盖全省,统计完 成后,各地均具有自己的降雨量属性) 。 空间插值是指利用研究区已知数据来估算未知数据的过程, 即将离散点的测量数据转换 为连续的数据曲面。利用空间插值,我们就可以通过离散的等降雨量线,来推算出山东省各 地的降雨量了。

二、空间插值的几种方法及本次实习采用的原理和方法
整体插值方法 ? 边界内插方法 ? 趋势面分析 ? 变换函数插值 – 局部分块插值方法 ? 自然邻域法 ? 移动平均插值方法:反距离权重插值 ? 样条函数插值法(薄板样条和张力样条法) ? 空间自协方差最佳插值方法:克里金插值 ■局部插值方法的控制点个数与控制点选择问题 局部插值方法用一组已知数据点(我们将其称为控制点)样本来估算待插值点(未知 点)的值,因此控制点对该方法十分重要。 为此,第一要注意的是控制点的个数。控制点的个数与估算结果精确程度的关系取决 于控制点的分布与待插值点的关系以及控制点的空间自相关程度。 为了获取更精确的插值结 果,我们需要着重考虑上述两点因素(横线所示) 。 第二需要注意的是怎样选择控制点。一种方法是用离估算点最近的点作为控制点;另 一种方法是通过半径来选择控制点,半径的大小必须根据控制点的分布来调整。 –

S6、按照不同方法进行空间插值,并比较各自优劣
打开 ArcToolbox——Spatial Analyst 工具——插值,打开插值方法列表,如下图:

A、采用反距离权重法(IDW)对降水量数据进行插值:
反距离权重法的特点是按照距离待插值点的远近核定已知数据点的权重,从而对待插 值点进行插值的过程。一个已知数据点距离待插值点越远,权重就越低,它的值对待插值点 的影响就越小。影响的程度用点之间距离乘方的倒数表示,通过“power”设置乘方。乘方 为 1 意味着点之间数值变化率为恒定, 称为线性插值法; 乘方为 2 或更高则意味着越靠近已 知点,数值的变化率越大。 这种插值方法的优点是对于数据分布均匀的区域,插值效果好;缺点是在数据分布不 均地区插值容易出现小的封闭等值线( “球状突起” )和因数据缺乏而产生的不规则等值线。 双击 ArcToolbox 里面的 “反距离权重法” , 输入点要素选择 “prec” , Z 值字段选择 “prec” , 输出像元大小选择 1000。点击确定,效果如下图:

由图像可以看出,山东省内陆地区由于布点均匀,插值效果比较好;但是沿海地区出 现两条“球状弧线” ,如下图:

B、采用克里金法对降水量数据进行插值:
克里金插值与 IDW 插值的区别在于权重的选择。IDW 仅仅是将距离的倒数作为权重, 而克里金考虑了空间相关性的问题。 克里金法首先将每两个点进行配对, 这样就能产生一个 以两点之间距离为自变量的函数。对于这种方法,原始的输入点可能会发生变化。在数据点 多时,内插的结果可信度较高。 1、双击 ArcToolbox 里面的“克里金法” ,输入点要素选择“prec” ,Z 值字段选择“prec” , 半变异函数属性中,设置克里金方法为普通克里金(Ordinary Krigging,OK) ,半变异模 型选择球面(Sphere) ;输出栅格命名为“ok_sphere” ,输出像元大小选择 1000。点击确 定,效果如下图:

发现,局部存在比较明显的“锯齿现象” ,如图:

2、同样是克里金插值方法,本次采用线形普通克里金插值,效果如图:

效果比较理想。 对于各种方法,依次展示如下: 高斯:

指数:

圆:

泛克里金——与一次漂移函数成线性关系:

泛克里金——与二次漂移函数成线性关系:

3、依次尝试各种克里金插值法,对比得出: 克里金方 法 半变异模型 球面 线性 普通克里 金 高斯 指数 圆 与一次漂移函数成线性 关系 与二次漂移函数呈线性 关系 插值效果 在数据点分布稀疏地区,锯齿现象和撕裂现象明显 结合了球面方法和指数方法的优点,插值效果比较 理想 数据范围被剧烈拉伸,不能很好地反映降水量分布 斑块状撕裂现象基本避免了,但是锯齿现象仍旧明 显 插值效果和线性类似,比较理想 中部插值效果比较理想,但边缘地带出现严重的数 据误差 数据拉伸剧烈;区分度小

泛克里金

C、采用薄板样条函数法对降水量数据进行插值:
薄板样条函数法使用样条函数来对空间点进行插值,它有两个基本条件: a.表面必须完全通过控制点(样本点) b.使所有点的坡度变化最小,换句话说,表面的二阶曲率是最小的。 该方法适用于高程面、水位面、气候数据(如平均降水量)的插值。有一个缺点是, 在数据贫乏的地区坡度较大,经常涉及如同过伸的情况。 样条函数插值分为两种,规则样条函数和张力样条函数。下面分别进行操作: 1、 双击 ArcToolbox 里面的 “样条函数法” , 输入点要素选择 “prec” , Z 值字段选择 “prec” , 样条函数类型选择:REGULARIZED(规则样条函数) ,输出像元大小选择 1000。点 击确定,效果如下图:

观察发现,规则样条函数也是对数据进行了一定程度的拉伸,但是拉伸幅度不大, 拉伸后的像元值介于 99—1108 之间 (相比之下, 普通克里金的高斯方法和泛克里金的与 二次漂移函数呈线性关系方法拉伸幅度分别达到了 -597—4903、-5347—4292) ,可以说 是拉伸范围适宜;在山东省内部区域插值分布比上述所有方法都要更加均匀一些;虽然

在边缘地带出现一定的数据突兀 (如下图, 色调偏白的区域年平均降雨量只有 200mm 以 下) ,但总体效果基本上能令人满意。

2、 双击 ArcToolbox 里面的 “样条函数法” , 输入点要素选择 “prec” , Z 值字段选择 “prec” , 样条函数类型选择:TENSION(张力样条函数) ,输出像元大小选择 1000。点击确定, 效果如下图:

可见得,张力样条插值不存在数据过拉伸现象,插值平滑且无突兀地带,插值效果 非常好。

D、采用自然邻域法对降水量数据进行插值:
这种方法的原理是构建 voronoi 多边形,也就是泰森多边形。首先将所有的已知数据点 构建一组 voronoi 多边形,然后将待插值点也构建一组 voronoi 多边形,这样后一组多边形 与原多边形有很多相交的地方。接下来,根据每一块的面积按比例设置权重,这样就能求得 待插值点的值了。 双击 ArcToolbox 里面的 “自然邻域法” , 输入点要素选择 “prec” , Z 值字段选择 “prec” , 输出像元大小选择 1000。点击确定,效果如下图:

可见,该方法在有数据的区域,插值效果非常好,优于上述各种方法;但是由于自然邻 域法是基于邻近点的算法,所以图上山东省轮廓边缘(如鲁西北、威海荣成成山头附近)出 现数据空白区,不符合最后成图要求。

S7、对插值效果最好的一种方法,调整插值参数
经过上述比较, 得出插值效果最理想的是薄板张力样条插值方法, 接下来调整搜索半径, 比较在不同的插值半径下不同的效果。 1、双击 ArcToolbox 里面的“样条函数法” ,输入点要素选择“prec” ,Z 值字段选择“prec” , 样条函数类型选择: TENSION (张力样条函数) , 输出像元大小选择 1000, 权重值填 0.1, 点数填 12。文件名设为 ten_A_12。点击确定,效果如下图:

只改变权重值,不改变其他参数,重新尝试。权重值填 2,点数 12。确定,效果如下图:

再次改变权重值,不改变其他参数。权重值填 12,点数填 12。点击确定,效果如下图:

通过对比,可得:对薄板张力样条插值方法来说,在选取同样多的控制点的前提下,权 重越小,生成的表面越光滑。例如,权重 0.1 的效果>权重 2 的效果>优于权重 12 的效果。 2、同理,对比在同一权重值下,控制点数目不同所产生的效果差异: 在权重为 0.1 的前提下,我们设置控制点数量为 5 个、30 个,跟控制点 12 个的作比较: ■5 个控制点:

■30 个控制点:

可见得,控制点 5 个和 30 个的情况均不如 12 个的时候效果好,也就是说,控制点不一 定要太多,也不一定要少,要结合已知数据点的个数和分布来判定最佳方案。因此我们最终 选择权重 0.1、12 个控制点的图像进行插值。

S8、制作山东省县域年平均降水量地图
1、首先制作一个山东省轮廓图。保存一份山东省县界图的副本,命名为 sd.shp,鼠标点击 编辑器——开始编辑,将所有要素选中,点击编辑器下拉菜单中的“合并” ,生成的文件 就是山东省轮廓范围了;

2、利用山东省轮廓范围矢量图对 ten_A_12 栅格图进行裁剪。 打开数据管理工具——栅格— —栅格处理——裁剪,弹出“裁剪”对话框。输入栅格选择 ten_A_12,输出范围选择 sd, 在“将输入要素用于裁剪几何”处打对勾,目的是为了按照山东省的边界来裁剪降水量 图(否则,默认是按照矢量多变性的最小外接矩形来裁剪栅格图像) 。如下图:

3、现在生成了一幅灰度图像。我们还是按照颜色对其进行分类:

鼠标右键点击 sd_prec,选择属性——符号系统——已分类,弹出计算直方图对话框,询 问是否计算直方图,选择“是” 。目前最大值是 835,最小值是 490,数据跨度为 345, 因此我们按照每 25mm 为一个降水量带,对图像进行分类制图。345÷25+1≈15,因此 我们将分类设为 15 档,色带选择蓝色调,如图:

4、鼠标右键点击 sd_county——属性——标注,勾选“标注此图层中的要素” ,使县市的名 称显示在地图中; 5、按照地图符号及视觉特征实现的要求,对地图进行设计和整饰,最终成图如下:


相关文章:
arcgis空间内插值教程_图文.doc
arcgis空间内插值教程 - GIS 空间插值 (局部插值方法) 实习记录 一、空间插值的概念和原理 当我们需要做一幅某个区域的专题地图, 或是对该区域进行详细研究的...
ArcGIS 克里金空间插值方法介绍_图文.ppt
ArcGIS 克里金空间插值方法介绍_计算机软件及应用_IT/计算机_专业资料。克里金插值法 龚资林 2014.10.10 基本思想: 克里金插值与IDW(反距离权重法)插值的区别在于...
ArcGIS 自然邻域空间插值方法介绍2_图文.ppt
ArcGIS几种空间插值方法介绍 自然临域插值法 自然临域插值法 基本思想: 自然邻域法插值工具使用的算法可找到距查询点 最近的输入样本子集,并基于区域大小按比例对...
ARCGIS插值操作_图文.ppt
ARCGIS插值操作_理学_高等教育_教育专区。实验4.1 空间数据统计、插值 1.空间数据统计 GIS/LIS数据库中的专题数据进行统计分析包括 频数和频率 属性数据的集中特征 ...
ARCGIS实习课程之空间内插实习_图文.ppt
ARCGIS实习课程之空间内插实习_天文/地理_自然科学_专业资料。ARCGIS实习课程之...介于最 大值和最小值之间 23 4、对AMT进行空间插值- Spline(样条)插值 ?...
12 ArcGIS中几种空间内插方法的比较_图文.pdf
12 ArcGIS中几种空间内插方法的比较_数学_自然科学_专业资料。第 29 卷 第 ...1. 1 反距离加权法 反距离加权插值 ( Inverse D istance W eighting, 简称...
ArcGIS中空间数据统计、插值_图文.ppt
ArcGIS中空间数据统计、插值_工学_高等教育_教育专区。ArcGIS中空间数据统计、插值 实验4.1 空间数据统计、插值 1.空间数据统计 GIS/LIS数据库中的专题数据进行统计...
ArcGIS中几种空间内插方法的比较_图文.pdf
ArcGIS中几种空间内插方法的比较_互联网_IT/计算机_专业资料。第 29 卷 第 ...1. 1 反距离加权法 反距离加权插值 ( Inverse D istance W eighting, 简称...
ArcGIS中几种空间插值方法.doc
反距离加权法(IDW) ArcGIS 中最常用的空间内插方法之一,反距离加权法是以插值点与样本点 之间的距离为权重的插值方法, 插值点越近的样本点赋予的权重越大,其...
ArcGIS零基础者快速实现空间插值.doc
ArcGIS 零基础者快速实现空间插值以太湖贡湖湾水质 数据 Spline 插值为例 ...教程所用软件:ArcGIS10.1 英文版,其他版本部分具体操作稍有差异,步骤基本一致。...
ARCGIS中几种空间插值简单比较.doc
ARCGIS中几种空间插值简单比较_设计/艺术_人文社科_专业资料。ARCGIS 几种空间插值简单比较 (2012-01-10 22:09:14) 1.IDW。基本思想是目标离观察点越近则...
ArcGIS实验操作(十四)---采样数据的空间内插.pdf
ArcGIS实验操作(十四)---采样数据的空间内插_专业资料。ArcGIS教程ArcGIS 实验操作(十四)采样数据的空间内插(Interpolate) 空间插值常用于将离散点的测量数据转换为连...
ARCGIS插值操作_图文.pdf
ARCGIS插值操作_数学_自然科学_专业资料。 文档贡献者 ljb0034908 贡献于2015-03-21 1 /2 相关文档推荐 ARCGIS插值方法、原理 62页 1下载券 ArcGIS中几种空间...
arcgis教程_ArcGIS栅格数据的空间分析_图文.ppt
arcgis教程_ArcGIS栅格数据的空间分析 - 第13章 ArcGIS栅格数据的空间分析 基于栅格数据的空间分析是GIS空间分析的重要组成部分。 其中,缓冲区分析在GIS中用得较多,...
ArcGIS中几种空间内插方法的比较_图文.pdf
( 2007) 01) 0072) 04 ArcGIS中几种空间内插方法的比较秦 涛, 付宗堂 (...反距离加 权插值 ( Inverse D istance W eight ing , 简称 ID ) 作为 一...
arcgis栅格数据空间分析实验报告课案_图文.doc
arcgis栅格数据空间分析实验报告课案_中职中专_职业教育_教育专区。arcgis 实验五 一、 实验目的 栅格数据的空间分析 理解空间插值的原理,掌握几种常用的空间差值分析...
ArcGIS地理信息系统空间分析实验教程PPT-第10章 地统计....ppt
ArcGIS地理信息系统空间分析实验教程PPT-第10章 地统计分析解读_中职中专_职业...? ? ? ? ? 10.1 地统计基础 10.2探索性数据分析 10.3空间确定性插值 ...
arcgis4 空间基本操作_图文.ppt
采样点数据空间插值问题: ? 根据气象观测站点数据,生成某一区域内降雨 量等值线...Arcgis8-9空间基本操作 78页 2下载券 ARCGIS教程 实验五、空间... 27页 ...
第13章 ArcGIS栅格数据的空间分析_图文.ppt
13.1 设置分析环境在ArcGIS中,使用栅格数据进行空间分析前,先要设 置分析环境...1.核函数密度制图 2.简单密度制图 13.4.1 反距离加权插值反距离加权插值(...
基于ArcGIS两种空间插值方法的比较.doc
基于ArcGIS两种空间插值方法的比较_计算机硬件及网络_IT/计算机_专业资料。2008 ...2004,28(2):184 [4] 汤国安, 杨昕.ArcGIS 地理信息 系统空间分析实验教程[M...
更多相关标签: