文档介绍：第 27卷第 1期自然资源学报 Vol27No 1
2012年 1月 JOURNALOFNATURALRESOURCES Jan.,201 2
土壤连续属性空间插值方法及其精度的研究进展
史文娇1,岳天祥1,石晓丽2,3,宋伟1
(,北京 100101;,石家庄 050016;
,石家庄 050016)
摘要:作为土壤变化的时空定量监测方法,土壤属性空间插值方法及其精度是计量土壤学和
“数字土壤”领域的重要研究内容。文章首先介绍了土壤属性空间插值的常用方法,包括克立格
插值法(Kriging)、反距离权重法(IDW)、样条插值法(Spline)、贝叶斯最大熵(BME)、地理加权
回归(GWR)以及高精度曲面建模方法(HASM);其次阐述了土壤属性空间插值精度验证的方法
和指标;再次总结了能够提高土壤属性插值精度的 6种途径,包括合理选择插值方法、准确设定
插值方法参数、合理设计采样数目和密度、注意空间自相关程度和范围的影响、科学安排实验分
析顺序以及结合辅助信息进行插值;最后从插值方法的选择、验证指标的选取以及辅助信息的
结合三个方面指出了土壤属性空间插值方法及其精度的未来研究方向。
关键词:计量土壤学;地统计学;土壤空间插值方法;高精度曲面建模(HASM);插值精度
中图分类号:S159 文献标志码:A 文章编号:1000-3037(2012)01-0163-13
作为土壤变化的时空定量监测方法,土壤属性空间插值方法及其精度是计量土壤学
(Pedometrics)[1]和“数字土壤”(DigitalSoil)[2]领域的重要研究与方向[3]。由于土壤母质、
土壤类型、地形、气候、生物等环境条件的差异及人类活动的影响,土壤属性存在空间异质
性[4]。地理信息系统(GIS)、地统计学和曲面建模逐渐成为适合农业、环境管理等研究领域
的常用工具。近些年来,由于精准农业的实施和土壤污染管理的需要使得这些技术进一步
发展。例如,土壤肥力制图的质量影响定位土壤肥力管理的效果[5],土壤重金属分布图的
质量影响区域环境保护措施的具体实施等[6]。土壤属性空间插值能够为土壤的可持续利
用及环境保护的宏观决策提供技术支持,但在一定程度上依赖于土壤属性空间插值的精度。
因此,土壤属性的空间插值方法、精度及其验证是土壤属性空间插值中必须考虑的因
素[5,79]。本文系统总结了土壤连续属性空间插值及其精度的研究进展,并阐述了提高土壤
属性插值精度的方法,最后探讨了未来的研究趋势。
1 土壤属性空间插值方法
11 克立格插值法(Kriging)
Kriging插值法是一种无偏线性最优估值方法。对变量在点 x处的估计值 Z(x),可以
通过该点影响范围内的 n个有效观测值 Z(xi)的线性组合得到,即:
收稿日期:2011-05-10;修订日期:2011-07-20。
基金项目:国家自然科学基金青年基金项目(41001057);国家杰出青年科学基金项目(40825003);地表过程与资源
生态国家重点实验室开放基金项目(2011-KF-06)。
第一作者简介:史文娇(1982-),女,辽宁葫芦岛人,博士,研究方向为土壤属性空间模拟、全球变化与区域农业。
Email:******@lreis.
书
1 64 自然资源学报 27卷
n

Z (x)=∑λiZ(xi) (1)
i=1

式中:λi是赋予观测值 Z(xi)的权重,表示各观测值对估计值 Z (x)的贡献,在保证估计值
无偏性(即估值偏差的平均值为 0)和最优性(即估值方差最小)条件下,可由变量的半方差
函数计算得到[10]。其半方差函数可用下式表示:
N(h)
1 2
γ(h)= N(h)×∑[Z(xi)-Z(xi+h)] (2)
2 i=1
其中:γ(h)是半方差,N(h)是距离等于 h时的点对数,Z(xi)是在位置 xi处的数值,Z(xi+
h)是在距离(xi+h)处的数值。
在进行半方差函数的拟合之前,需对数据进行正态分布检验,因为数据的非正态分布会
影响半方差函数的精度,即当样品的平均值增加时,样品的方差也增加,使实际变异函数产
生畸变、基台值和块金值增大,导致某些结构特征不明显[10]。如不呈正态分布,需对其进行
变换。正态变换的常用方法主要包括对数变换、BoxCox变换和正态计分变换(NormalScore
Transformation)等。
Kriging方法包括单变量 Kriging方法和与