中国农业科学 ›› 2016, Vol. 49 ›› Issue (11): 2126-2135.doi: 10.3864/j.issn.0578-1752.2016.11.009

• 土壤肥料·节水灌溉·农业生态环境 • 上一篇    下一篇

黄土高原煤矿区复垦农田土壤有机质含量的高光谱预测

南 锋,朱洪芬,毕如田   

  1. 山西农业大学资源环境学院,山西太谷030801
  • 收稿日期:2016-01-28 出版日期:2016-06-01 发布日期:2016-06-01
  • 通讯作者: 毕如田,Tel:0354-6288912;E-mail:birutian@163.com
  • 作者简介:南锋,Tel:0354-6286586;E-mail:nanfeng24@126.com
  • 基金资助:
    国土资源部公益性行业科研专项(201411007)、山西农业大学科技创新基金项目(201307)

Hyperspectral Prediction of Soil Organic Matter Content in the Reclamation Cropland of Coal Mining Areas in the Loess Plateau

NAN Feng, ZHU Hong-fen, BI Ru-tian   

  1. College of Resources and Environment, Shanxi Agricultural University, Taigu 030801, Shanxi
  • Received:2016-01-28 Online:2016-06-01 Published:2016-06-01

摘要: 【目的】针对黄土高原丘陵地多、地形复杂、有机质含量低、采样困难以及因采煤活动引起大面积土地损毁等问题,在土地复垦与综合整治过程中,为快速定量监测与评估复垦农田土壤质量提供一种新的方法。【方法】以山西省襄垣县复垦农田土壤为研究对象,选取由北向南土地损毁中间条带状区域采集样品152个,进行室内土壤农化分析、光谱测定,运用ParLes 3.1软件对光谱曲线进行多元散射校正(multipication scatter correction,MSC)、基线偏移(baseline offset correction,BOC)和Savitzky-Golay filter平滑去噪预处理。对土壤原始光谱反射率(raw spectral reflectance,R)作一阶微分(first order differential reflectance,D(R))和倒数的对数变换(inverse-lg reflectance ,lg(1/R)),分析3种不同变换形式的光谱数据与土壤有机质含量的相关性,相关系数通过P=0.01水平显著性检验来确定显著性波段的范围。基于全波段(400—2400 nm)和显著性波段利用偏最小二乘回归(partial least squares regression,PLSR)分析方法建立该区域土壤有机质含量高光谱预测模型,通过模型精度评价指标:决定系数(coefficient of determination,R2)、均方根误差(root mean square error,RMSE)和相对预测偏差(residual prediction deviation,PRD)确定最优模型。【结果】通过P=0.01水平显著性检验的波段范围为:R的400—1 800、1880—2 400 nm;D(R)的420—790、1 020—1 040、2 150—2 200 nm;lg(1/R)的400—1 830、1 860—2 400 nm。光谱与有机质含量的相关系数绝对值最大的波段是R的800 nm;D(R)的600 nm;lg(1/R)的760 nm。进行D(R)变换,光谱曲线的吸收特征更加明显,相关系数在可见光(400—800 nm)波段范围内有所增加,其最大值由0.72提高到了0.82;基于显著性波段的PLSR建模效果优于全波段,其中lg(1/R)变换的预测精度为最佳,具有很好的预测能力,其校正模型的R2和RMSE分别为0.95、7.64,预测模型的R2、RMSE和RPD分别为0.85、3.00、2.56;基于全波段的R-PLSR和lg(1/R)-PLSR模型具有较好的预测能力,其预测模型的R2、RMSE和RPD分别为0.79、3.64、2.10和0.79、3.53、2.17,而D(R)-PLSR模型只能进行粗略估测,其预测模型的R2、RMSE和RPD分别为0.61、5.43、1.41。综合分析全波段和显著性波段3种光谱数据的预测精度,发现基于显著性波段的R-PLSR、D(R)-PLSR、lg(1/R)-PLSR模型均取得了显著的预测效果。【结论】研究区土壤光谱反射率与土壤有机质含量具有高度的相关性,应用偏最小二乘回归分析方法可以很好地建立土壤有机质含量反演模型。

关键词: 煤矿区, 复垦农田, 土壤有机质, 高光谱, 偏最小二乘回归

Abstract: 【Objective】 In terms of the problems in the Loess Plateau, such as many hills, complex topography, low soil organic matter content (SOMC), sampling difficulties, large areas of land damage caused by mining activities and so on, the object of this study is to provide an alternative method for the rapidly quantitative monitoring and evaluation of the SOMC in the process of land reclamation and comprehensive renovation. 【Method】 Taking the cropland soil in the coal mining areas in Xiangyuan County, Shanxi Province was picked as research object, 152 soil samples were collected from the intermediate strip area of land destruction region in a north to south direction. The physical and chemical properties of the soil samples were analyzed. At the same time, the raw hyperspectral reflectance (R) of the soil samples was measured by the standard procedure with an ASD FieldSpec 3 instrument equipped with a high intensity contact probe under the laboratory conditions. The raw spectral reflectance (R) were pretreated by the smoothing or denoising methods of multiplication scatter correction (MSC), baseline offset correction (BOC) and Savitzky-Golay filter in the ParLes 3.1 software. And the raw spectral reflectance (R) was transformed into two types of spectra, which were first order differential reflectance (D (R)) and inverse-log reflectance (lg (1/R)), to analyze the correlation coefficients between the three spectra and their SOMC. Then the significant bands were extracted by the significant correlation coefficients (P=0.01) of the three spectra with the SOMC. Finally, based on the full bands (400-2 400 nm) and significant bands of the three spectra, the hyperspectral predicting models of the SOMC were established by the method of partial least squares regression (PLSR). The optimal models were determined by the assessing indices of predicting accuracies, including coefficient of determination (R2), root mean square error (RMSE), and residual prediction deviation (RPD). 【Result】 The spectra in the bands of 400-1 800 and 1 880-2 400 nm for the raw spectral reflectance (R), 420-790, 1 020-1 040, and 2 150-2 200 nm for D (R), and 400-1 830 and 1 860-2 400 nm for lg (1/R), were significantly correlated with SOMC (P=0.01). And the maximum correction coefficients between the three spectra and their SOMC were 800 nm of the raw spectral reflectance (R), 600 nm of D (R), and 760 nm of lg (1/R). After the transformation of D (R), there were prominent differences among the absorption peaks of the spectral curves in different soil samples, and their correlation coefficients were improved from the value of 0.72 to that of 0.82 in the range of visible bands (400-800 nm). The models of significant bands could obtain better predicting accuracies compared with that of full bands by the method of PLSR. Among the three spectra, the predicting accuracy of lg (1/R) was the best, and R2, RMSE of the calibration dataset were 0.95 and 7.64, while R2, RMSE, and RPD of the validation dataset were 0.85, 3.00, and 2.56, respectively. For the models of R-PLSR and lg (1/R)-PLSR of full bands, the predicting abilities were good. The R2, RMSE, and RPD of R-PLSR were 0.79, 3.64, and 2.10, respectively. And the coefficient of R2, RMSE and RPD of lg (1/R)-PLSR were 0.79, 3.53, and 2.17, respectively. However, for the model of D (R)-PLSR, the predicting SOMC were only roughly estimated, and the indices of the predicting accuracies were not satisfying. R2, RMSE and PRD of the D (R)-PLSR were 0.61, 5.43, and 1.41, respectively. Finally, by analyzing the predicting accuracies of the three spectra in both full bands and significant bands, it was found that the models of R-PLSR, D (R)-PLSR and lg (1/R)-PLSR in significant bands achieved desirable predicting effect. 【Conclusion】 In the study area, soil spectral reflectance has a high correlation with SOMC, and PLSR is a good method to establish the predicting model of SOMC.

Key words: coal mining areas, reclamation cropland, soil organic matter content (SOMC), hyperspectral, partial least squares regression (PLSR)