野外空间采样的渐进式策略.pdfVIP

第33卷第5期 武 汉 大学 学 报 ·信 息 科 学 版 Vo1．33 NO．5 2008年 5月 Geomatics and Information Science of Wuhan University May 2008 文章编号：1671—8860(2008)05—0441—05 文献标志码：A 野外空问采样的渐进式策略 张景雄 Michael F Goodchild (1 武汉大学遥感信息工程学院，武汉市珞喻路 129号，430079) (2 加利福尼亚大学圣芭芭拉分校地理学系，美国CA 93106—4060，UsA) 摘 要：介绍了有关地统计学(geostatistics)的基本理论，论述了顾及样本点位之间空间依赖性的空间采样方 案的设计。基于总体上极大限度地降低克里金方差的思想，阐述了一个快速确立采样点位的序贯算法和块段 克里金方法。由于野外作业的计算环境所限，数值方法和启发式搜索相结合的策略才是开发实用系统的明智 选择。以假想数据为例，验证了所提算法的有效性。 关键词：空间采样；协变量；块段克里金法；协同克里金法 中图法分类号：P237．3 地统计学在空间信息科学中具有重要地位，也 越广泛的应用。常用于空间预测的一个重要方法 给空间采样方案设计提供了理论基础。传统的采 就是克里金法，该方法在给出最优线性无偏估值 样是非空间的，即不考虑样本所在的空间位置；空 的同时，还可计算出估计值的量化评价指标—— 间采样设计时，应考虑采样点位的分布。在实际工 估计方差，即克里金方差(Kriging variance)。 作中，采集数据的间隔不能太大，也不能太小，间距 克里金法是一种基于最优权重的空间估值方 越小，采样成本越高；反之，采样间距越大，则成本 法，位于z处的变量z的克里金估值是由周围一 越小，但数据不一定能反映研究区的真实特征。因 定区域内采样点的已有数据经线性加权组合得到 此，野外采样通常是以动态方式进行的，在保证成 的。如图1所示，假设有一样本由”个采样点组 本尽量小的前提下，降低估计量的方差。对于某一 成，各采样点的观测值为z(z )(i一1，2，…，”)， 个空间区域，分散在一定距离内的点位的某些特性 观测值的均值为m ，则未采样的位置z的克里金 通常会呈现空间相关性，即空间依赖性¨1]。也就是 估计值z(z)为： 说，在一定邻域范围内的采样点是相关的，它们所 z(z)一m + (z(z )一m ) (1) 提供的信息有冗余。因此，在空间采样设计时，须 i=l 考虑采样点间的相关性，使不确定性达到最小。 式中， (i一1，2，…，”)对应邻近搜索区域内位于 在已有一组采样点的前提下，如何确定后续 z ( 一1，2，…，”)点位的样本权重系数。将 (i 采样点的位置，并且分析整个问题域的克里金方 一 1，2，…，”)组成的列向量记为w，则w可由下 差和实施效率／成本随采样点累积的动态效应，是 式解算： 本文的主旨所在。在野外采样实施过程中，由于 COV (z 32 )W — COV (xx ) (2) 计算机软、硬件和网络传输等方面都受到很大的 估计值z(z)的不确定性可用克里金方差来 限制，因而需要研究一种高效、优化的算法。 衡量：