基于克里金的水深插值模型比较.pdfVIP

第 15卷 第4期 2015年2月 科 学 技 术 与 工 程 Vol_15 No．4 Feb．2015 1671～ 1815(2015)04—0135—05 ScienceTechnologyandEngineering ⑥ 2015 Sci．Tech．Engrg． 基于克里金的水深插值模型比较 陆婷婷 林 敏 (天津大学建筑工程学院，天津 300072) 摘 要 对于河道及近海区域的水动力计算和工程规划，水深的确定和地形的模拟有着重要的作用，而对于某些水域，只存 在离散的水深值，则可用Kriging方法对离散值进行插值。因此，将水深作为区域化变量，根据普通Kriging插值方法，选定不 同的模型进行比较，包括高斯模型，球面模型，五球模型，孔洞模型以及指数模型。依据渤海湾的250个水深值，并借助ArcGIS 软件进行模型参数分析。结果显示球面模型的精度最高，表明在水深插值中运用球面模型具有更好的效果。 关键词 普通Kriging 水深值 插值模型 中图法分类号 TP391．75； 文献标志码 A 随着科技的进步，船舶航行、港 口建设、海洋工 的规律进行预测。地统计学上假设所有样本随机产 程等重要活动对海底三维地形精度都具有一定的要 生，样本之间具有内在的联系。 求。但是由于技术 以及特殊地理环境的影响，实测 水深可以表示为空问点函数 海域的水深值都为离散的。又因为海底地形具有平 h=l厂(，Y) (1) 滑性，与普通克里金插值方法的特点相符合，所以， 其呈现一定的空间分布，为明显的区域化变量， 基本采用克里金插值方法构造出相应的三维高程模 具有随机性和结构性。 型，为需求提供精确的数据参考。 1．1 变异分析 目前，地统计分析被运用于许多领域，在现代空 1．1．1 半变异函数 间统计学中扮演重要的角色。而地统计分析又能与 半变异函数为地统计分析中的函数，又称半变 GIS分析联系在一起，使得复杂的地统计过程能够 异矩 ，半变差函数。定义为 在ArcGIS中轻易的实现 lJ，为可视化提供 了理论 (，h)=@Var[Z()一z(+h)] (2) 基础。运用ArcGIS软件 ，可以对插值模型精度进行 量化 。 即 克里金法是基于统计模型和数学模型的一种空 y(，h)=÷E[z()一z(+h)]一 间插值_3，4J，将克里金法运用到水深插值时，为了得 到精确度较高的插值结果，就必须针对现有的数据 1{z(x)一E[Z(x)一z(x+h)]} (3) 进行分析，选取精度较高的模型，在一般 的克里金插 又因为Z()满足内蕴平稳假设，因此 值方法上提出更适合水深插值的模型。 目前国内外关于克里金插值方法应用进行很多 (，h)=÷E[Z(x)一z(x+h)] (4) 研究 。本文利用 ArcGIS软件 ，基于渤海湾水 1．1．2 协方差函数 深的250个水深数据(离散点)，在普通克里金方法 概率论 中定义随机变量 和 y的协方差为 下，利用不同的模型进行数据分析，并检验模型的合 Cov(X，Y)=E{[ 一E(y)][y—E(y)]}(5) 理性，找寻在水深插值中最优的克里金