应用最大熵方法改进克里金负权系数.pdf

第３５卷 第２期 桂 林 理 工 大 学 学 报 Ｖｏｌ３５Ｎｏ２ ２０１５年５月 ＪｏｕｒｎａｌｏｆＧｕｉｌｉｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ Ｍａｙ　２０１５  文章编号：１６７４－９０５７（２０１５）０２－０３４０－０３　　　　　 ｄｏｉ：１０３９６９／ｊ．ｉｓｓｎ１６７４－９０５７２０１５０２０２０ 应用最大熵方法改进克里金负权系数 １ １ ２ 陈天伟 ，卢献健 ，江海东 （１桂林理工大学 测绘地理信息学院，广西 桂林　５４１００４； ２中南大学 ａ地球科学与信息物理学院；ｂ有色金属找矿预测教育部重点实验室，长沙　４１００８３） 摘　要：针对基于地统计学的克里金数学模型，利用最大熵原理分析模型系统对空间权重系数的依赖度增 设非负约束条件，研究消除构建格网ＤＥＭ时的负权问题。通过 ＭＡＴＬＡＢ编程验证算法的正确性、准确 性，并与拟牛顿法、内点算法解算二次规划的结果进行了比较。对比结果显示，最大熵法得到的权系数大 小比例与点位关系相适应，且其估值精度优于后两者。 关键词：克里金法；负权；最大熵 中图分类号：Ｐ２０８；Ｐ２０９　　　　 　 　　　　　　　文献标志码：Ａ 数字高程模型 （ＤＥＭ）作为空间信息系统的 即模型误差最小，那么熵值应该最大，使得插值 一个重要组成部分，蕴含着大量的、各种各样的 模型更能反映空间数据的离散性，则模型会最稳 地形结构和特征信息，是国家空间数据基础设施 健，由此可以修正权估值系数的解算，避免负权 的基本产品之一和赖以进行地形分析的基础数据。 系数的产生。 但目前网格ＤＥＭ建模的质量仍然不尽如意，其原 因主要与建模方法的局限性有关。 １　原　理 构建格网ＤＥＭ，普遍采用克里金法、杨赤中 最大熵理论的数学模型为 法等内插方法进行网格点的高程估计。由于上述 ｎ 插值数学模型中都含有权重函数的约束条件，即 　　ｍａｘＳ［Ｐ］＝－ ＰｌｏｇＰ， （１） ∑ ｉ ２ ｉ ｉ＝１ Ｗ ＝１，在解算过程中，权重值 Ｗ时常会出现 ｎ ∑ ｉ ｉ ｓ．ｔ． Ｐｇ（ｘ）＝Ｅ（ｇ），ｊ＝１，２，…，ｍ； 负值。负权重值在地学上没有任何意义，并且可 ∑ ｉｊ ｉ ｊ ｉ＝１ 　　　　Ｐ ０，ｉ＝１，２，… ，ｎ。 能会随之产生负的估计值，模型也变得高度不稳 ｉ≥ 定，这不仅给工程带来较大误差甚至错误，也影 其中：ｇ（ｊ＝１，２，…，ｍ）表示各阶统计矩函数；