arcgis空间内插值教程

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

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

3、横线所示）。第二需要注意的是怎样选择控制点。一种方法是用离估算点最近的点作为控制点；另 一种方法是通过半径来选择控制点，半径的大小必须根据控制点的分布来调整。S6、按照不同方法进行空间插值，并比较各自优劣打开ArcToolboxSpatial Analyst工具插值，打开插值方法列表，如下图:匚小描值I- 克里金去i 港去- 含隆需的样条函数，地fi缀栅格卜*.样条函数法 *气自然邻域法-5趋势酉法A采用反距离权重法（IDW0对降水量数据进行插值:反距离权重法的特点是按照距离待插值点的远近核定已知数据点的权重，从而对待插值点进行插值的过程。一个已知数据点距离待插值点越远， 权重就越低，它的值对

4、待插值点 的影响就越小。影响的程度用点之间距离乘方的倒数表示，通过“ power”设置乘方。乘方 为1意味着点之间数值变化率为恒定， 称为线性插值法；乘方为2或更高则意味着越靠近已 知点，数值的变化率越大。这种插值方法的优点是对于数据分布均匀的区域，插值效果好；缺点是在数据分布不均地区插值容易出现小的封闭等值线（“球状突起”）和因数据缺乏而产生的不规则等值线。双击ArcToolbox里面的“反距离权重法”，输入点要素选择"prec”,Z值字段选择"prec”, 输出像元大小选择1000。点击确定，效果如下图：由图像可以看出，山东省内陆地区由于布点均匀，插值效果比较好；但是沿

5、海地区出现两条“球状弧线”，如下图:如，TI-iW：MaitIJMI 此:画电鼻U rHU ©Li D ES胃1r:审Fn>-Ataw同Km MLLLU312 M期期询解 2回口皿-KNr1苓匚 A V 1二W口乂：*SiHuZ r L，* n0r：i L Kre»im ；、B Q 叫常«- «*,«*«A - iHMtB、采用克里金法对降水量数据进行插值:克里金插值与IDW插值的区别在于权重的选择。IDW仅仅是将距离的倒数作为权重，而克里金考虑了空间相关性的问题。克里金法首先将每两个点进行配对，这样就能产生一个以两点之间距离为

6、自变量的函数。对于这种方法，原始的输入点可能会发生变化。在数据点多时，内插的结果可信度较高。1、双击ArcToolbox里面的"克里金法"，输入点要素选择"prec”, Z值字段选择"prec”, 半变异函数属性中，设置克里金方法为普通克里金( Ordinary Krigging , OK),半变异模 型选择球面(Sphere);输出栅格命名为"ok_sphere"，输出像元大小选择 1000。点击确 定，效果如下图：发现，局部存在比较明显的“锯齿现象”，如图:2、同样是克里金插值方法，本次采用线形普通克里金插值，效果如图:TU却%1

7、用蛀if U w仲iMh的F“阳m ” '2 niHUjja n FIHW" TB MUTJVH IL i 回阳fn & FWE.+J效果比较理想。对于各种方法，依次展示如下:高斯:口也* r ；卡，叶 3加鬲jsi鼻却白 >崎.鼠丁 ;.： *，酢，：* 田二力/1闻.«: kAhui -!?/*- «Mk； - - * ttdiR-、，":*»力-nr勺 PMWfWU图 * o里 口Q ifh.nAn,白«K"vda 口 胆， ui.<* o A2k»4MEHWl亚耳BMdJ >

8、;Vc -1匚W” 而初0!9*2*同H P« fb G! "*."*7口 口"£nm井lXM雷z同. msg»*，H ofc_kdb dKjUuMl - BdLiriwq P QlLj*yMjiM nw K* 2d.L一 eu n 2工a8aE3 Hiwflff圆:-工I SE nfli ag Hflui E9Q-B1mi bhq R孙 »>u+ 11 FW.EMWr a '.a*B辛 K法叵奇昌it3H矽口(THli =< r ： ： ： * 事 r It O 日 MJHLEA £ 机口 I

9、a - AB' , ae«s a*S!ID4 AkJUb.El &m.r niEuxa *.*加KJ1 iZ£- B«T.n E5J »- 3IS 134d473M：75泛克里金一一与一次漂移函数成线性关系:口一总i 息JM." Zra-WFIM- JQ 以 3L7.1US1 Xi 晦/闺问 = 占的MJ4 *5 -412.4UMK- -<T .43K71M -M SMJMM4 * M7.HZUU - 3C SK-mmS - £j- 总g Hi MlnrLrSB E ELaTWE - ri .7评4餐|一口-77

10、1 "g2LH4 . ?lt IIALHMW - Si- IHUniili M 抽冲aiAl贺骷 -VSL3MZZ1J - W «KAff7«W - LJC L-” LX O &_DMk* E日与W菱MILLE灯也落泛克里金一一与二次漂移函数成线性关系: dllim E 岬 / I -=0二 A *?事，Q :工，4 事. * MF 1工 鼠J."r 1 t * * WHI *3*h! -； liWi- 4 U A'网 -a J 龊 4 4 上二, B0 中网 I - -5 WdWVTVJ QI JJHmiJM WJKLKTTH . Mt

11、fmMTH - gi>mM 叩mm - UFn-jnno? -IMJ-HJCfl - -L LttlM*I 1 -317MS3H3-> TLVM»HI I i+ri-!H77M - H3、依次尝试各种克里金插值法，对比得出:克里金方 法半变异模型插值效果普通克里 金球面在数据点分布稀疏地区，锯齿现象和撕裂现象明显线性结合了球面方法和指数方法的优点，插值效果比较 理想数据范围被剧烈拉伸，不能很好地反映降水量分布指数斑块状撕裂现象基本避免了，但是锯齿现象仍旧明 显圆插值效果和线性类似，比较理想泛克里金与一次漂移函数成线性 关系中部插值效果比较理想，但边缘地带出现严重的数 据误

12、差与二次漂移函数呈线性 关系数据拉伸剧烈；区分度小C采用薄板样条函数法对降水量数据进行插值:薄板样条函数法使用样条函数来对空间点进行插值，它有两个基本条件：a.表面必须完全通过控制点（样本点）b.使所有点的坡度变化最小，换句话说，表面的二阶曲率是最小的。该方法适用于高程面、水位面、气候数据（如平均降水量）的插值。有一个缺点是，在数据贫乏的地区坡度较大，经常涉及如同过伸的情况。样条函数插值分为两种，规则样条函数和张力样条函数。下面分别进行操作：1、双击ArcToolbox里面的“样条函数法”，输入点要素选择"prec”,Z值字段选择"prec”, 样条函数类型选择：REGUL

13、ARIZED （规则样条函数），输出像元大小选择1000。点击确定，效果如下图：丽w lfMLZZm-74 0星u 7WWDM F nLLLNnai k KrLUUAU iH+，国工常EiM聿口 A”. y ；： ： ： * M 用用 4 =力力果 耳金RlMd> - . r/ Jk 'Mbtaia - - J Mto-、 I - A - : ：. r B 4电上,上，.心嘉於ita I 4观察发现，规则样条函数也是对数据进行了一定程度的拉伸，但是拉伸幅度不大，拉伸后的像元值介于 991108之间（相比之下，普通克里金的高斯方法和泛克里金的与二次漂移函数呈线性关系方法拉伸幅度分别

14、达到了-597 4903、-53474292）,可以说是拉伸范围适宜；在山东省内部区域插值分布比上述所有方法都要更加均匀一些；虽然在边缘地带出现一定的数据突兀（如下图，色调偏白的区域年平均降雨量只有200mm以下），但总体效果基本上能令人满意。2、双击ArcToolbox里面的“样条函数法”，输入点要素选择"prec”,Z值字段选择"prec”, 样条函数类型选择：TENSION （张力样条函数），输出像元大小选择1000。点击确定， 效果如下图：I<皿毡LU7NM .IM*5M - JIMia期川期巾9啪新 曲J电t片 Hu 2潮娴 薛团涉喻UM做 曲*1rt4 f

15、! H*"«E心口 明金小,tad - %, A -可见得，张力样条插值不存在数据过拉伸现象，插值平滑且无突兀地带， 插值效果非常好。D采用自然邻域法对降水量数据进行插值：这种方法的原理是构建voronoi多边形，也就是泰森多边形。首先将所有的已知数据点构建一组voronoi多边形，然后将待插值点也构建一组voronoi多边形，这样后一组多边形与原多边形有很多相交的地方。接下来，根据每一块的面积按比例设置权重，这样就能求得待插值点的值了。双击ArcToolbox里面的“自然邻域法”，输入点要素选择“prec”，Z值字段选择“prec”, 输出像元大小选择1000。点击确定，

16、效果如下图:吗乂产k a 二 II Mr. A黑犀岫-，匚 / * WHlt -0 _WWV：FJld而, 3 E 3M1) 3g 与成雄 mid -=W1 E他 i珏网哂0问 HTi：R *(*!m-+ I I ai. HA；MlUl- Rr ' c - A. - ."3 id_EW*ftr 口L_1 &iGISajg-i GImm !.C 口 PTK口得0 577J77777I：a sn.mrFMi .空巧珏5骐 HB15S5&55?.制43Hmm 【能&科上於M gLLLIHll 66131111L2 4mm够> 也喻粉"聚班版

17、11反领茄JjUH3a 必4乜*5 77-2223 TTL21H223 - KOK O yJrjrn S*TE 睢口安可见，该方法在有数据的区域,插值效果非常好,优于上述各种方法;但是由于自然邻域法是基于邻近点的算法，所以图上山东省轮廓边缘（如鲁西北、威海荣成成山头附近）现数据空白区，不符合最后成图要求。S7、对插值效果最好的一种方法，调整插值参数接下来调整搜索半径,经过上述比较，得出插值效果最理想的是薄板张力样条插值方法,比较在不同的插值半径下不同的效果。1、双击ArcToolbox里面的“样条函数法”,输入点要素选择"precZ值字段选择“ prec”,样条函数类型选择:TENS

18、ION （张力样条函数）,输出像元大小选择1000,权重值填0.1,，4UMd 口 i Q 1口 ” . Mkaiiit只改变权重值，不改变其他参数，重新尝试。权重值填2,点数12。确定，效果如下图:点数填12。文件名设为ten_A_12。点击确定，效果如下图:. 4TliD*WM*FL lEcjy mriaor - ac.1 unnu攸 Mwwuara M&UTHAI - EJMHN HLJMdG - NflMa&ia?讣鼻4 .,IB勺 就留、倒 sa*Ed¥35*W!.H P* EnKIWi0+15 ML4MC- - -4U-KHH3I ll.nTMK- 1WH

19、4U4谀电网小即蝌 3S4R0J - WUHZin Mi&ium n»j-xa?aL TKUXm TMJ-aSN?m EflW* - g*g Ef UW幻，刷岫MlJl H iw.U3_UJE1 - RLtHUi/- -Ma.nimr串.|1 Ci | 0 - -再次改变权重值，不改变其他参数。权重值填 12,点数填12。点击确定，效果如下图:通过对比，可得：对薄板张力样条插值方法来说，在选取同样多的控制点的前提下，权重越小，生成的表面越光滑。 例如，权重0.1的效果权重2的效果优于权重 12的效果。2、同理，对比在同一权重值下，控制点数目不同所产生的效果差异：在权重为0.1

20、的前提下，我们设置控制点数量为5个、30个，跟控制点12个的作比较： 5个控制点：可见得，控制点5个和30个的情况均不如12个的时候效果好，也就是说，控制点不一定要太多，也不一定要少，要结合已知数据点的个数和分布来判定最佳方案。因此我们最终选择权重0.1、12个控制点的图像进行插值。S8、制作山东省县域年平均降水量地图J3WI3即 -4MIU.：33 M.1、首先制作一个山东省轮廓图。保存一份山东省县界图的副本，命名为 编辑器一一开始编辑，将所有要素选中，点击编辑器下拉菜单中的就是山东省轮廓范围了；sd.shp,鼠标点击“合并”，生成的文件0 史 lOII.Ed 5 Mep旧 呆jaw *Ei

21、i|HI*川 HTAjOQ #(*! IRSH m n j* u a r _j t 二1的，/Hjln -c *<Li25iW - 4tMy dMJH 箝1 3ITaMiai|HT.M9M，5i . SJDJJiMae- M4-L3WH4 65-T.47LLM3 T¥7.3JM4« - tmPWi1IUM377TI 皿 30个控制点:吗?x 0 FZ , *：曰 *®irq， I WH- * U* A -1图*it怯训a-口 BE. o次i5-1*9/忖0Hli已产图 4bJHH33 - 3n.?HM3TJM3H 3 +Kai7Ul JCi9 - ：m 1

22、； - Ui 电 F 声!HBXIC1 ,>如 » 六b g«4nM"，团E5U 4WJ71ZM3 Wjj猊必喀帕困一 - wRDTnUIXIDlliiL编*J*；!出 MhMmjnvau骑 1M1-T3KM R1_5«nH U1JK3IL5 - H71.IM4SI5i DK M3X mmJlRhBg4ui鼠B中百看口 3-f尸 r-狐. 建W启.口:正认*£5 dd-UHfs. =1 M 由 口2、利用山东省轮廓范围矢量图对ten_A_12栅格图进行裁剪。打开数据管理工具一一栅格一一栅格处理一一裁剪，弹出“裁剪”对话框。输入栅格选择te

23、n_A_12,输出范围选择sd,在“将输入要素用于裁剪几何”处打对勾，目的是为了按照山东省的边界来裁剪降水量图（否则，默认是按照矢量多变性的最小外接矩形来裁剪栅格图像）。如下图:3、现在生成了一幅灰度图像。我们还是按照颜色对其进行分类:-鼠标右键点击sd_prec,选择属性一一符号系统一一已分类，弹出计算直方图对话框，询 问是否计算直方图，选择“是”。目前最大值是 835,最小值是490,数据后夸度为345, 因此我们按照每 25mm为一个降水量带，对图像进行分类制图。345+ 25+1 = 15,因此我们将分类设为15档，色带选择蓝色调，如图：3 0 sdjprec岳04897696533 - 50050