李鑫星， 曹闪闪， 白雪冰， 李 辉

1. 中国农业大学信息与电气工程学院食品质量与安全北京实验室， 北京 100083 2. 中国农业大学信息与电气工程学院， 北京 100083

引 言

土壤是农业生产的基础， 不同的土壤成分含量对作物的生长有不同的影响， 对于作物而言， 土壤中的营养成分供应的丰缺情况至关重要。 不同地区的土壤由于环境不同导致土壤成分含量有所不同， 土壤成分含量的不均会造成土壤退化、 土壤污染， 影响作物生长等问题[1-2]， 如： 土壤含水量的多少影响作物的生长； 土壤有机质直接影响土壤肥力； 土壤重金属很难被分解， 长期积累转化为有毒化合物， 会造成土壤污染， 农作物自身的自净能力有限， 对农作物的生长有一定影响， 农作物转化为粮食可能间接危害人体健康； 土壤盐分含量过多会造成土壤盐渍化， 引起土壤板结、 肥力下降； 土壤酸碱度影响土壤中的化学反应， 土壤过酸或过碱都造成植物体内元素失衡， 影响植物生长； 土壤缺少氮、 磷、 钾使作物光合作用减弱， 造成叶片发黄等问题制约着我国农业的可持续发展。

多光谱技术是近年来发展较为迅猛的高新分析技术， 国内外很多学者将它应用于土壤学研究领域， 利用无人机搭载多光谱传感器采集土壤样本， 对样本数据进行实时动态监测和分析。 多光谱成像技术能够在宽光谱范围内的窄光谱带上进行检查， 从而可以直接评估具有诊断重要性的不可见或低对比度特征。 传统的土壤成分检测多采用化学分析方法， 存在人为误差， 因此将多光谱技术应用到土壤成分含量的检测中是一个新的发展方向。

本文综述了近6年来国内外多光谱技术在土壤成分含量检测领域应用的相关文献， 分析了多光谱成像技术的快速性、 准确性和便捷性等特点， 简述了多光谱技术对土壤成分含量的检测过程， 即： 土壤样品制备、 采集多光谱图像、 图像预处理、 特征提取、 建立模型， 重点阐述了多光谱技术在土壤成分含量检测中的研究进展， 分析多光谱在土壤的水分、 有机质、 氮磷钾及其他成分(土壤盐分和重金属)中的研究进展及存在的问题， 并展望了多光谱技术在土壤成分含量检测中的未来发展趋势， 为土壤成分含量检测研究提供一定的参考。

1 土壤成分多光谱成像检测技术概述

土壤中营养成分含量的变化直接影响作物的生长， 因此需要对土壤成分含量进行实时检测， 并能预测土壤成分含量的变化趋势， 及时采取措施对土壤成分含量进行调整。 传统的数据采集基于物理手段和化学方法， 非常费时、 费力， 为了高效采集土壤成分含量的数据， 采用多光谱成像的方法获取样本图像信息和光谱信息。

1.1 多光谱技术特点

多光谱成像技术是一种结合了光谱学和成像技术的方法， 能够有效的检测土壤成分含量， 与传统的土壤成分含量检测方法相比有如下特点：

(1)快速性： 多光谱技术能够同时获取土壤样本的空间信息和光谱信息， 不需要物理或化学处理， 利用多光谱成像仪直接获取多光谱图像， 仅需10～20 s。

(2)便捷性： 利用无人机搭载多光谱相机直接照射需要检测地区， 通过分析直接获取该地区土壤含量信息， 具有便捷性。

(3)准确性： 多光谱技术通常采用10～20个工作波段， 根据土壤表面的光谱反射率曲线， 使用多光谱成像技术重建光谱反射率， 从而获取土壤成分含量的准确信息。 例如， 陈硕博等[3]根据光谱反射率曲线建立多元模型， 对模型进行分析得到相关系数达到0.9以上。

1.2 多光谱技术检测过程

多光谱成像检测技术是一种间接检测方法， 应用多光谱对土壤成分定性或定量分析， 检测过程可划分为5个阶段， 即： 土壤样品制备、 采集多光谱图像、 图像预处理、 特征提取、 建立模型[4]。

图1 土壤成分多光谱技术检测流程图

(1)土壤样品集选择： 土壤样品集选择是土壤分析的关键环节， 直接关系到分析的结果[5]： ①选择制备覆盖度较低的区域， 对多光谱反演影响较小； ②根据地形随机和均匀地选择采样点， 保证研究结果的准确性。

(2)采集多光谱图像： 使用无人机搭载多光谱相机， 获取土壤的光谱图像， 为了使检测结果更加准确与稳定， 必须严格控制测试时间、 路径、 高度、 光照、 仪器参数等因素。

(3)数据预处理： 多光谱建模过程中存在许多不稳定因素， 常常会受到一些与待测样品性质无关的信息的干扰， 为了建立一个更加稳定、 可靠的土壤多光谱模型， 需要对多光谱数据进行预处理[6]。 数据预处理方法中效果较好的有： 多元散射校正、 Savitzky-Golay平滑算法、 小波分析等。

(4)特征波段选择： 在土壤成分含量的多光谱反演过程中， 特征性波段选择即是找出地表反射率与土壤成分含量相关性高的波段[7]。

(5)建立模型： 使用不同的建模方法对采集的多光谱图像建模， 分析得出用于检测土壤成分含量的最合适模型。 目前较为常用数学建模方法包括非线性建模方法， 如循环神经网络(recurrent neural network, RNN)和线性建模方法如多元线性回归(multivariable linear regression, MLR)， 偏最小二乘法(partial least squares, PLS)和主成分分析(principal component analysis, PCA)。

2 土壤成分多光谱成像检测

2.1 水分检测

土壤含水量是影响农作物生长的一个重要指标， 土壤含水量偏高或偏低都会影响农作物的正常生长。 土壤含水量的快速估算对于干旱半干旱地区的精准农业具有重要意义[8-10]， 因此研究一种简单、 快速、 准确的方法测定土壤含水量非常重要。

国内外学者利用多光谱成像对土壤含水量进行实时、 准确的检测， 并已取得显著成果。 高瞻等[11]以关中塿土为研究对象， 制备含水量为1%～25%的土壤样本， 采集反射光谱， 建立模型， 结果表明通过偏最小二乘法和RBF神经网络建立模型可以更好地检测土壤水分。 Hassan-Esfahani等[12]利用无人机载高光谱相机获取图像数据建立人工神经网络模型， 估计表面土壤水分的有效性， 相关系数达到88%。 Helene等[13-14]利用光谱成像仪， 获取田间土壤表面绿度及冠层温度， 利用水分胁迫图反演土壤含水率， 显示整个季节所需灌溉量， 以确保作物健康生长。 Petropoulos等[15]利用遥感数据大面积测定土壤水分。 张智韬等[16]利用多光谱遥感技术对土壤含水率大范围快速检测， 结果表明土壤表层约1 cm处是检测土壤含水率的最佳检测深度。 Wang等[17]采用相关系数法对表层土壤含水率的敏感带进行筛选， 对单个敏感带的光谱数据构建一元回归模型并分析定量关系， 预测相关系数达0.92以上， 为多光谱技术监测土壤表层含水率提供了一条新途径。

每年总用水量的80%都用于作物灌溉， 但实际利用率仅40%， 由于这些不当的灌溉方式导致水资源严重浪费, 多光谱技术能够有效的大范围动态监测土壤含水量。 但是， 采集土壤多光谱数据时容易受到土壤地貌、 植被覆盖率和当日气候等因素的影响， 并且多光谱在可见光和近红外光谱区域内波段较少， 因此可以与高光谱数据融合提高反演精度， 进行互补式土壤水分反演。

2.2 有机质检测

土壤有机质(soil organic matter， SOM)是土壤成分中供农作物生长必不可少的营养成分之一, SOM在很大程度上影响大气中的CO2浓度甚至全球碳循环[18-21]， 其含量的多少与土壤肥力和植被状况直接相关。

Liu等[22-24]提出基于相似性与深度函数相结合的方法， 预测土壤中有机质浓度的三维变化。 袁征等[25]利用光谱对土壤有机质含量定量估测， 采用不同的建模方法分析， 其中模糊识别模型决定系数达到0.973， 优于其他模型精度。 Anne等[26]采用光谱与偏最小二乘回归建模方法， 预测佛罗里达州西海岸湿地土壤有机质的含量， 相关系数达到0.81。 陈思明等[27]为了提高湿地土壤有机质含量的预测精度， 分别用不同多光谱的不同波段反射率测量土壤有机质的含量， 通过3种不同的数据模型(多元逐步回归方法、 BP神经网络和支持向量机)建立土壤有机质反演模型， 通过实验数据分析得出基于支持向量机模型为光谱检测土壤成分含量的最佳模型。 该研究为湿地土壤有机质含量的预测精度提供了技术可行性。 王锐等[28]利用LandSat8多光谱遥感技术分析地表反射率与土壤有机质的相关性， 建立遥感反演模型计算土壤有机质含量分布， 结果表明， 地表反射率与土壤有机质含量呈负相关， 并且呈现不规律分布。

多光谱技术对土壤有机质光谱反演机理研究的检测准确度较高[6]， 但是由于土壤的复杂性， 光照强弱的影响及植被的覆盖率和种类对多光谱反演土壤有机质数学模型有一定的影响， 因此在今后的研究中数学建模技术的进一步创新仍占有很大的发展空间。

目前较为常用数学建模方法包括非线性建模方法RNN和线性建模方法MLR， PLS和PCA。 它们在建模中有各自的特点及适用范围， 如表1所示。

表1 四种建模方法的对比分析

2.3 氮、 磷、 钾检测

土壤成分中的氮、 磷、 钾对作物生长有重要的作用， 适量的成分对作物的光合作用有一定的促进作用， 快速准确地估算土壤的氮、 磷、 钾含量有助于促进处方施肥。 土壤成分中不同营养元素在不同光谱波段上的敏感度不同， 如近红外波段对土壤中的氮、 磷元素比较敏感， 而可见光能更好的检测土壤中的钾元素[35]。

Thielebruhn等[36]利用可见-近红外光谱估算土壤总氮， 使用竞争性自适应重加权采样和迭代保留信息变量算法， 结合偏最小二乘回归， 提高预测精度。 Zhang等[37]测量土壤样品的实时近红外吸收光谱， 使用连续体去除土壤水分干扰， 利用小波分析确定6个敏感波段， 建立总氮含量回归模型实时预测土壤中的总氮含量。 Morellos等[38]利用移动式纤维型可见光和近红外光收集土壤光谱， 建立非线性数据模型。 结果表明Cubist方法提供了土壤总氮的最佳模型， RMSEP=0.071, RPD=1.96。 万余庆等[39]通过分析0.4～2.5 μm之间光谱反射率与土壤氮磷钾含量的相关性， 得到最密切反射率波长为500， 800和760 nm， 为大范围处方施肥提供依据。 Xu等[51]利用随机森林方法建立了土壤总氮预测模型， 探讨了不同超分辨率贝叶斯法(PAN)光谱指数对土壤总氮预测模型的影响， 其中基于多光谱和GS PAN光谱指数的土壤总氮模型获得较高的预测精度。

光谱技术能够快速、 方便的实时检测土壤的氮磷钾含量， 但是土壤氮的检测仍然受许多因素影响， 其中土壤水分含量的多少、 多光谱采集图像光照强度的影响以及土壤颗粒的大小等因素干扰土壤总氮含量的检测， 影响模型的校准精度。

2.4 其他成分检测

除上述成分含量检测之外， 土壤中的盐分和重金属也进行了相关的测定。

土壤盐分含量过多会形成土壤盐渍化， 引起土壤板结、 肥力下降。 土壤盐渍化是全世界发生的主要土壤退化威胁之一[40-41]， 土壤盐分含量监测同样受到许多专家学者的关注。 王明宽等[42]以黄河三角洲为采样点， 实地采集土壤样本， 采用Landsat8影像建立相应的模型得出光谱反射率与土壤盐分含量并不仅仅是线性关系。 张雅莉等[43]利用多光谱影像提高土壤盐分反演精度， 建立多元逐步回归土壤盐分估算模型。 研究表明： 重采样光谱数据的估算模型的决定系数达到0.779。 Zhang等[44]使用多光谱相机获取黄河三角洲土壤盐度数据， 建立土壤盐度估算模型， 最佳模型为SAVI线性模型， 相关系数达到0.747， 为土壤盐分估算提供了一种快速有效的方法。 土壤盐碱化跟干旱一样影响农作物生产， 因此土壤酸碱度也是一个重要指标。 李云等[45]利用多光谱图像的颜色特征建立土壤酸碱度偏最小二乘预测模型， 结果显示决定系数达到0.91。

土壤中含有多种重金属成分[46-47]， 随着土壤重金属含量的增加， 造成土壤污染， 导致微生物量下降， 影响作物生长， 甚至间接影响人类健康， 人们逐渐意识到土壤重金属含量的检测与治理的重要性。 郭云开等[48]根据多光谱传感器的光谱响应函数， 采用实测ISI921VF反射光谱数据模拟Landsat卫星ETM+传感器多光谱数据， 在模拟光谱的基础上， 通过光谱特征提取、 构建土壤指数对土壤重金属进行预测分析。 Chen等[49]基于重金属与各种土壤性质和反射光谱的关系， 表明土壤光谱更适合作为内插重金属浓度的协变量。 Khumaeni等[50]利用发射光谱线信噪比的方法检测土壤中Cu和Hg含量。

土壤中的重金属元素在其存在的环境中并不能稳定下来， 由于植物对重金属的累积、 自净等原因， 重金属会随着条件的变化进行迁移或富集。 多光谱技术的快速性能够实时检测土壤中的重金属。 但是， 土壤重金属检测不仅要注意它们的总量， 还必须分析各种形态的含量， 因此， 土壤重金属的检测还有有待进一步研究。

3 发展趋势

多光谱技术在土壤成分含量检测中的应用前景十分广泛， 利用多光谱技术检测土壤成分含量的研究很多， 取得了一定的成果， 为土壤成分含量的检测提供了良好的技术支持。

(1)机器学习算法建模可增加检测的准确性、 快速性、 鲁棒性。 在多光谱土壤成分检测中， 结合机器学习算法中的无监督和监督模型分析各种数据驱动模型与不同实际测量环境下使用的输入数据相结合的潜力， 机器学习提供了一种数据驱动的解决方案， 减少了土壤成分分布不均的光谱数据对建模结果的影响。

(2)多光谱图像融合在数字土壤制图中是一个新的研究方向， 正处于起步阶段， 多光谱图像与全色图像相结合获取多光谱全色波段， 多光谱图像融合也称全色锐化(PanSharpening)。 PanSharpening技术在土壤成分含量检测中， 能够提高预测模型的分析精度和准确度， 需要更多的研究来探索多光谱数字图像融合对土壤预测模型的影响， 为加快实现实时在线土壤成分含量检测提供技术参考。

(3)存在多种光谱数据算法， 目前将多种算法结合处理数据占据主导， 主要以2种或2种以上的算法的结合的数据处理方法为主， 在今后较长一段时间内， 这种方法仍是处理数据的主流方法。 将多种算法相结合， 对传统数据处理的算法进行改进， 能够更好地发挥这些算法的优势， 获取一种更好的结合算法能够精确、 快速的解决更多的土壤成分含量检测的问题。