连续小波变换定量反演土壤有机质含量
王延仓1,3, 张兰1,3, 王欢1,3, 顾晓鹤2,4,*, 庄连英1,3, 段龙方1,3, 李佳俊1,3, 林靖1,3
1. 北华航天工业学院, 河北 廊坊 065000
2. 北京农业信息技术研究中心, 北京 100097
3. 河北省航天遥感信息处理与应用协同创新中心, 河北 廊坊 065000
4. 农业部农业信息技术重点实验室, 北京 100097
*通讯联系人 e-mail: guxh@nercita.org.cn

作者简介: 王延仓, 1986年生, 北华航天工业学院讲师 e-mail: yancangwang@163.com

摘要

以北京市东部地区96个潮土土样的土壤参数及对应光谱数据为数据源, 采用连续小波多尺度分析处理与分析。 首先将土壤光谱进行初步处理, 生成小波系数, 其次将土样的有机质含量与小波分解系数开展相关性分析, 提取特征波段, 最后采用特征波段建立预测耕层有机质含量的模型。 结果表明: 经连续小波处理后, 光谱对耕层有机质含量的预测能力明显优于传统光谱变换技术; 经连续小波分解后, 对土壤有机质含量的预测能力随光谱分辨率降低呈先降后升再降的趋势; 连续小波分析算法可提升土壤光谱对有机质含量的估测能力, 与土壤高光谱反射率相比, 基于连续小波变换的土壤有机含量最佳的精度提高19%; 由于光谱分辨率为80 nm建立的模型精度较高, 其 R2达到0.632, 这表明在连续小波算法下, 光谱分辨率较低的宽波段数据可用于土壤有机质含量的监测。

关键词: 土壤有机质; 连续小波变换; 高光谱
中图分类号:S153 文献标志码:A
Quantitative Inversion of Soil Organic Matter Content Based on Continuous Wavelet Transform
WANG Yan-cang1,3, ZHANG Lan1,3, WANG Huan1,3, GU Xiao-he2,4,*, ZHUANG Lian-ying1,3, DUAN Long-fang1,3, LI Jia-jun1,3, LIN Jing1,3
1. Institute of Computer and Remote Sensing Information Technology, North China Institute of Aerospace Engineering, Langfang 065000, China
2. National Engineering Research Center for Information Technology in Agriculture, Beijing 100097, China
3. Aerospace Remote Sensing Information Processing and Application Collaborative Innovation Center of Hebei Province, Langfang 065000, China
4. Key Laboratory of Information Technology in Agriculture, Ministry of Agriculture, Beijing 100097, China
Abstract

In this study, the data sourced from hyperspectral data of 96 tidal soil samples in Miyun, Tongzhou and Shunyi Districts of Beijing are processed and analyzed by means of continuous wavelet multiscale analysis technique. Firstly, the hyperspectral data are decomposed to generate wavelet coefficients and the correlation between the coefficients and soil organic matter content is analyzed, and the characteristic band is selected. Finally, the model to estimate soil organic matter content is constructed by using the characteristic band. The research results show that the estimation of soil organic matter by the reflectivity of soil spectrum is better than that of the traditional spectral transformation technology after continuous wavelet transformation. The ability of estimating soil organic matter by continuous wavelet decomposition decreases first and then increases with the reduction of spectral resolution. The results of continuous wavelet analysis can improve the ability to estimate the content of organic matter by the soil spectrum. Compared with the high spectral reflectivity of soil, the accuracy of soil organic content based on continuous wavelet is improved by 19%. Since the model accuracy is higher when built with the spectral resolution of 80 nm, its R2 reaches 0.632, which indicates that the wide band data can be used for the monitoring of soil organic matter content by using the continuous wavelet technique.

Keyword: Soil organic matter; Continuous wavelet transformation; Hyperspectral
引 言

精准农业即精准农作, 其将信息与知识相结合实现农作的知识化与信息化的现代农耕系统。 而精准施肥是精准农业的重要构成部分。 精准施肥, 即以作物产量为目标, 将作物生长规律与耕层土壤养分空间变异状况相结合, 定时、 定量的施肥技术[1]。 故, 实现快速、 准确、 无损地对大尺度耕层土壤实施监测, 具有迫切而重要的现实意义。 遥感技术凭其大尺度、 快速、 精准、 无损的优势, 为实现耕层土壤养分的监测, 提供了新的思路。 目前, 将遥感技术应用于耕层土壤养分监测的研究已成为热点[2, 3, 4, 5, 6, 7]

有机质是耕层土壤重要构成成分, 为农作物提供必需的营养成分, 具有涵养水源, 改良土壤等作用。 因此, 开展耕层土壤有机质含量的监测的研究具有重要意义[8]。 借助光谱技术监测耕层有机质含量已有大量工作, 且已取得一定的研究成果: 位于VIS— NIR区域内的光谱基本涵盖估测有机质含量的有益信息[9, 10, 11, 12, 13, 14]; 耕层有机质含量的监测精度一定程度上受其含量高低的影响, 但2%并非耕层土壤有机质含量监测的临界点; 此外, 还利用航空影像开展针对耕层土壤有机质含量监测的研究, 并取得一定研究结果[15, 16]。 但针对基于宽波段数据开展土壤有机质含量监测的研究相对较少, 且对最适宜光谱分辨率仍存在异议。

受时空变异与气候影响, 土壤类型多且构成成分复杂。 利用光谱对耕层土壤养分的监测, 须依据土壤类型构建相应监测模型。 但目前研究多偏于单一土壤类型(如东北黑土), 而对华北地区潮土的研究偏少, 且未取得较满意的研究成果。 近年来, 小波分析算法得到了长足发展, 且在信号分析领域得到了广泛应用, 如何将小波分析算法用于土壤养分监测的研究, 已成为当前研究的热点[17, 18, 19, 20]。 本研究采用连续小波算法对光谱开展分析, 提取有益信息, 建立土壤有机质含量的预测模型, 并探讨在连续小波算法下, 不同光谱分辨率的对土壤有机质含量预测精度的影响。

1 实验部分
1.1 土样采集

实验区位于北京市的密云区、 通州区及顺义区, 土壤类型为潮土, 土样采集均在2011年与2012年10月份开展, 共采集96份土样; 土样点均匀分布于实验区内(图1)。 土壤有机质含量测定方法与土样处理流程、 方法参照本刊34卷第7期1922页。

图1 研究区域与样点空间分布Fig.1 Thespatial distribution of regional and sample points

表1 土壤有机质含量的统计性描述 Table 1 The statistical description of soil organic matter content
1.2 光谱测量

土壤高光谱数据采用ASD(analytical spectral devices)生产的地物光谱仪采集, 光谱覆盖范围为350~2 500 nm, 分辨率为1 nm。 土壤光谱测量在暗室内开展, 其测量方式如图2所示, 将土壤置于黑色容器内, 并用直尺刮平土壤表面, 土样厚度为2 cm以上; 光源采用可提供平行光束的1 000 W的卤光灯, 其光线与垂线夹角为15° ; 采用5° 视场角的光纤探头进行收集, 并将探头置于土壤上垂直距离为30 cm处。 测量前用白板定标。 每个土样均获取10条光谱, 并将其均值作为该土样的光谱。

图2 实验室内土壤光谱测量示意图Fig.2 The sketch map of soil spectrum in laboratory

1.3 光谱处理

由于环境不可控因素、 仪器暗电流及不同传感器响应精度存在差异, 致使测量的高光谱数据混入噪声信息; 噪声信息的加入导致光谱边缘有效信息受损, 因此删除350399和2 4512 500 nm波段。 为削弱噪声引起的负面影响, 提升光谱数据的信噪比, 采用海明窗低通滤波器去噪处理[21], 并将光谱数据重采样至5 nm。

1.4 连续小波分析

小波分析是源于傅里叶算法的信号处理技术, 可在频率和时间同时开展分析, 能从信号中有效提取可用信息, 而非成像高光谱数据分析与电子信号分析具有较高类似性, 因此, 小波分析技术为光谱数据的处理与分析创建了新的思路。 小波分析分为连续小波和离散小波, 由于随离散小波分解尺度的增加, 土壤光谱维数减半, 未能对固定土壤光谱数据开展较好的分析, 故本研究采用连续小波分析(continuous wavelet transformation, CWT)将土壤高光谱数据进行不同尺度分解, 生成一系列小波系数。

Ψa, b=1aΨλ-ba(1)

式(1)中, a为伸缩因子, b为平移因子, λ 为土壤高光谱数据的波段数。 利用连续小波变换将土壤高光谱数据分解后, 可生成在不同分解尺度与不同波长上的小波系数。

Wf(a, b)=(f, Ψa, b)=-+f(λ)Ψa, b(λ)dλ(2)

式(2)中, f(λ )为土壤光谱反射率, 小波系数Wf(a, b)包含二维, 分别为波长(3502 500)与分解尺度(1, 2, …, 10), 故小波系数行为尺度数, 列为波长数的矩阵。

1.5 模型精度验证

光谱经连续小波变换处理后, 利用偏最小二乘法(partial least squares regression, PLS)构建有机质含量的估测模型。 采用随机抽样法对土样分为建模组、 验证组, 其中建模组用于构建估测模型共64个土样, 验证组用于检验预测模型的精度与稳定性共32个土样; 模型精度采用决定系数(R2)与均方根误差(RMSE)共同评价。

2 结果与讨论
2.1 光谱分析

图3为土壤光谱经连续小波处理后的结果, 其中图3(a) 为5 nm, 图3(b)为160 nm。 由图3(a)可知, 1~6尺度的土壤光谱起伏变化较大且无明显规律, 它们均在原始土壤光谱吸收处有明显下凹现象; 7~8尺度土壤光谱呈上凸的光滑曲线, 9~10尺度的光谱呈锯齿状。 由图3(b)知, 在2~3尺度的土壤光谱为光滑的上凸曲线; 3~4尺度呈分别呈曲线与直线状态, 6尺度呈现直线, 究其原因为光谱分辨率降低所致。 图3(a)与(b)对比分析可知, 随着光谱分辨率的降低, 经连续小波分解后的光谱曲线形状逐步呈反转状态。

图3 土壤光谱特性Fig.3 The soil spectral characteristics

2.2 相关性分析

图4为不同光谱分辨率的光谱经连续小波分解后, 其与有机质含量的相关系数矩阵。 在相关系数矩阵内总存在相关系数变化不大的区域, 该区域简称“ 信息平原” 。 由图知, 有效光谱信息主要集中于低尺度分解, 而在高尺度分解较少; 在光谱分辨率方面, 随光谱分辨率的降低, 光谱对土壤有机质的敏感性也逐渐降低的趋势; 在分解尺度上, 随着光谱分辨率的降低, 信息平原呈向低尺度移动的现象。 表明连续小波算法可深度挖掘宽波段内的有效信息。 光谱分辨率为5~20 nm时, 有效信息主要集结于1~5分解尺度, 而6~10分解尺度内的信息相对均一; 光谱分辨率为40~160 nm时, “ 信息平原” 逐步下移, 有效信息在高尺度逐步凸显。 究其原因是由宽波段内的有效光谱信息被隐藏, 故针对宽波段数据适度增大解尺度可有效发现其内的隐含信息。 由5, 10和20 nm相关性矩阵可知, 部分1和2级相关性较低的区域在3, 4和6级的相关性得到了较大提升, 这表明连续小波变换可深度挖掘有效信息。

图4 土壤有机质含量与小波系数的相关性Fig.4 The correlation between soil organic matter content and wavelet coefficients

表2 光谱反射率与土壤有机质含量的偏最小二乘回归模型列表 Table 2 The list of partial least squares regression models of spectral reflectance and soil organic matter content
2.2 反射率建模

采用11类传统光谱处理技术分析土壤光谱特性, 采用偏最小二乘算法建立土壤有机质含量预测模型, 并筛选6个较优模型作为参考。

2.3 连续小波变换建模

光谱信息经连续小波变换处理后, 采用偏最小二乘算法建立有机质含量预测模型结果如表3所示。 从表可知, 经连续小波处理的光谱构建的模型, 预测精度均较高, 且预测精度随光谱分辨率的降低而呈现先降后升的规律, 5 nm分辨率的预测精度最高, 分辨率为80 nm宽波段取得较为理想的精度; 由模型稳定性与精度综合分析, 可知基于5 nm构建的模型最优[散点图如图5(a)所示], 其次是光谱分辨率为80 nm的宽波段数据的模型[散点图如图5(b)所示], 其R2与RMSE分别为: 0.685, 0.068, 0.632, 0.058。 从整体分析, 基于连续小波变换分解的光谱信息建立的模型精度与稳定性均较高, 且随光谱分辨率的降低其对应模型的稳定性逐步增加, 故宽波段数据可有效压制高频噪声的负面影响。 由于80和160 nm构建的模型均取得了较高的精度, 因此, 连续小波变换可有效促进宽波段光谱数据用于监测土壤有机质含量。

表3 基于连续小波构建的土壤有机质含量估测模型列表 Table 3 The list of soil organic matter estimation models based on continuous wavelet

图5 真实值与估测值的散点图Fig.5 The scatter plots of real values and estimated values

2.4 连续小波变换与光谱反射率的对比分析

表2表3对比分析可知, 基于连续小波处理的光谱信息构建的模型的精度与稳定性明显较优, 究其源由为土壤构成成分复杂致使土壤有机质信息在土壤光谱信息中较弱; 而连续小波变换可有效的从土壤光谱信息中将有益信息分离, 从而提升可用信息的利用率。 与基于高光谱数据建立的模型对比分析可知, 经连续小波分解后, 光谱分辨率为80 nm的宽波段建立的模型精度与稳定性均较好, 表明连续小波变换适宜于宽波段数据的提取有效信息, 可用于分析宽波段数据; 与基于高光谱数据构建的模型相比, 经连续小波变换分解后, 估测土壤有机含量的最佳精度可提高19%, 宽波段最佳精度提高了10%。

在前人研究基础上, 以常用11种光谱变换技术方法为参照, 开展基于连续小波算法土壤有机质含量估测的研究。 研究表明: 与传统光谱变换技术相比, 该算法可较高的提升土壤光谱对土壤有机质含量的敏感性, 且模型精度与稳定性均较高。 也研究了光谱分辨率的降低对土壤光谱预测能力的影响; 研究表明在连续小波算法下, 5 nm分辨率预测精度最高, 其次为80 nm的宽波段。

本研究结果与其他相关研究存在差异, 究其原因, 主要为不同研究地理环境对土壤组分与土壤光谱均有较大影响, 且不同算法对土壤光谱内有效信息的整合与处理算法不同所致。 本研究结果可为利用遥感技术快速、 实时、 精准、 大尺度监测耕层土壤有机质含量的应用提供技术支撑, 但受实验条件制约与土壤有机质本身复杂性, 需在以下几方面作改进:

(1)土壤光谱决定于土壤各构成成分的共同作用, 故耕层土壤内其他成分必然对土壤有机质含量的光谱存在一定影响或干扰, 该干扰的影响需进一步的深入研究。

(2)土壤光谱是在相对可控的环境条件下获取的, 但自然界不可控因素较多, 复杂多变, 如何将实验室研究成果转为现实可用技术, 仍需开展深入研究。

3 结 论

采用传统光谱变换技术和连续小波技术对光谱分别开展处理和分析, 并建立预测有机质含量的模型, 旨在研究基于连续小波技术预测土壤有机质含量的最优宽波段, 结论如下:

(1)利用连续小波分析技术估测土壤有机质含量, 其最佳估测模型对应的分辨率是5 nm, 其R2=0.685, RMSE=0.300; 其宽波段的最佳模型对应的分辨率为80 nm, R2=0.632, RMSE=0.320。

(2)光谱经连续小波变换后, 其光谱信息对土壤有机质的估测能力随光谱分辨率的降低呈先降后增再降的趋势。

(3)连续小波分析算法可渗入挖掘土壤光谱内的有益信息, 提升对有机质含量的估测能力; 与土壤高光谱反射率相比, 经连续小波技术处理后, 光谱对土壤有机含量的最优估测精度提升19%, 宽波段的最优精度提升了10%。

(4)在连续小波算法下, 光谱分辨率为80 nm是预测土壤有机质含量的最优宽波段。

(5)基于光谱分辨率为80 nm建立的模型预测精度达到0.632, 这表明基于宽波段数据开展土壤有机质含量的监测是可行的。

The authors have declared that no competing interests exist.

参考文献
[1] CHEN Gui-fen, MA Li, CHEN Hang(陈桂芬, 马丽, 陈航). Journal of Jilin Agricultural University(吉林农业大学学报), 2013, 35(3): 253. [本文引用:1]
[2] MENG Ji-hua, WU Bing-fang, DU Xin, et al(蒙继华, 吴炳方, 杜鑫, ). Remote Sensing for Land & Resources(国土资源遥感), 2011, 90(3): 1. [本文引用:1]
[3] YANG Bang-jie, PEI Zhi-yuan, ZHOU Qing-bo, et al(杨邦杰, 裴志远, 周清波, ). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2002, 18(3): 191. [本文引用:1]
[4] ZHAO Ming-song, ZHANG Gan-lin, LI De-cheng, et al(赵明松, 张甘霖, 李德成, ). Acta Ecologica Sinica(生态学报), 2013, 33(16): 5058. [本文引用:1]
[5] CHEN Feng-rui, QIN Feng, LI Xi, et al(陈锋锐, 秦奋, 李熙, ). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2012, 28(20): 188. [本文引用:1]
[6] YANG Jian-feng, MA Jun-cheng, WANG Ling-chao(杨建锋, 马军成, 王令超). Geospatial Information(地理空间信息), 2015, 13(2): 47. [本文引用:1]
[7] BAI Ling-yu, ZENG Xi-bai, LI Lian-fang, et al(白玲玉, 曾希柏, 李连芳, ). Agricultural Sciences in China(中国农业科学), 2010, 9(11): 1650. [本文引用:1]
[8] Masserschmidt I, Cuelbas C J, Poppi R J, et al. Journal of Chemometrics, 1999, Chemometrics(13): 265. [本文引用:1]
[9] Moslem Ladoni, Hosein Ali Bahrami, Sayed Kazem Alavipanah, et al. Precision Agric, 2010, 11: 82. [本文引用:1]
[10] Krishnan P, Alexand er J D, Butler B J, et al. Soil Science Society of America Journal, 1980, 44: 1282. [本文引用:1]
[11] Alicia Palacios-Orueta, Jorge E Pinzon, Susan L Ustin, et al. Remote Sensing of Environment, 1998, 68: 138. [本文引用:1]
[12] Liu H J, Zhang Y Z, Zhang B, et al. Environmental Monitoring and Assessment, 2009, 154: 147. [本文引用:1]
[13] JI Wen-jun, SHI Zhou, ZHOU Qing, et al(纪文君, 史舟, 周清, ). Journal of Infrared and Millimeter Waves(红外与毫米波学报), 2012, 31(3): 277. [本文引用:1]
[14] HE Ting, WANG Jing, et al(何挺, 王静, ). Geomatics and Information Science of Wuhan University(武汉大学学报·信息科学版), 2006, 31(11): 975. [本文引用:1]
[15] Viscarra Rossela R A, Chappella A, de Caritatb P, et al. European Journal of Soil Science, 2011, 62: 442. [本文引用:1]
[16] Roberts D F, Adamchuk V I, Shanahan J F, et al. Precision Agricg, 2011, 12: 82. [本文引用:1]
[17] CHEN Hong-yan, ZHAO Geng-xing, LI Xi-can, et al(陈红艳, 赵庚星, 李希灿, ). Scientia Agricultura Sinica(中国农业科学), 2012, 45(7): 1425. [本文引用:1]
[18] Jahn B R, Brooksby P A, Upadhyaya S K. Transactions of the American Society of Agricultural Engineers, 2005, 48(6): 2065. [本文引用:1]
[19] LI Rui-ping, SHI Hai-bin, ZHANG Xiao-hong, et al(李瑞平, 史海滨, 张晓红, ). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2012, 6(3): 82. [本文引用:1]
[20] ZHAO Ze-hai, ZU Yuan-gang, CONG Pei-tong(赵则海, 祖元刚, 丛沛桐). Acta Ecologica Sinica(生态学报), 2002, 22(10): 1660. [本文引用:1]
[21] XU Yong-ming, LIN Qi-zhong, WANG Lu, et al(徐永明, 蔺启忠, 王璐, ). Acta Pedologica Sinica(土壤学报), 2006, 43(5): 709. [本文引用:1]
[22] Burrus C S, Gopinath R A, Guo H. Upper Saddle River, NJ: Prentice Hall, 1997. [本文引用:1]