基于最优数学和小波变换的绿洲耕层土壤盐分RF高光谱反演

    Hyperspectral inversion of RF model of soil salinity in oasis tillage layer based on optimal mathematics and wavelet transform

    • 摘要: 土壤盐渍化会导致土壤结构恶化、作物减产,影响农业生态系统的可持续性,而其监测和管理已成为盐渍化地区亟待解决的问题。为精准识别盐渍化区域以优化农田管理,以新疆渭干河-库车河绿洲耕层土壤为研究对象,测定其光谱反射率并进行数学变换、连续小波变换(continuous wavelet transformation,CWT)、离散小波变换(discrete wavelet transformation,DWT),以及数学变换结合CWT和DWT处理,并基于Pearson相关分析提取与土壤盐分含量(soil salt content,SSC)密切相关的特征波段,从而构建随机森林(random forest,RF)模型定量估算SSC。结果表明:1)在350~2450 nm范围内,不同盐渍化等级的土壤光谱反射率随SSC增加而升高,但达到极重度盐渍化时反射率降低。2)一阶微分变换、CWT、DWT以及一阶微分结合CWT和DWT处理均能显著提升光谱数据与SSC的相关性,其中DWT效果最优,相关系数最高达−0.621(P<0.01)。3)一阶微分结合CWT处理显著提升了模型的估测能力,以(1/R)′_CWT_28的全波段构建的估测模型最优,其验证集的决定系数为0.715,均方根误差为17.467 g/kg,相对分析误差为1.56。研究表明数学变换结合小波变换能提升模型的估测效果,可为SSC定量监测提供参考。

       

      Abstract: Soil salinization can often result in the continuous accumulation of soluble salt in the surface layer of soil, due mainly to natural and human activities. Saline soil has then posed a serious threat to the soil structure and crop yield in sustainable agriculture. It is very urgent to effectively monitor and manage the soil salinization in saline-alkali areas. Fortunately, hyperspectral imaging can be expected to estimate the soil salinity in recent years, due to its rich spectral information and high resolution. Among them, soil salt content (SSC) is one of the most important indicators to evaluate the degree of soil salinization. This study aims to estimate the soil salinity in the tillage layer of the oasis during hyperspectral imaging. The salinized areas were also precisely identified to optimize farmland management. A total of 193 soil samples were collected from the field in the Weigan-Kuqa River Oasis in Xinjiang, China. The salt content was then measured in the laboratory. Additionally, the spectrum of each soil sample was measured ten times and then averaged to obtain the spectral reflectance. The noisy spectral bands were also removed at the ends of the curves and two water absorption bands, followed by smoothing preprocessing using the Savitzky-Golay method. The spectral data was obtained to further transform using mathematical transformations, continuous wavelet transformation (CWT), and discrete wavelet transformation (DWT). According to Pearson correlation analysis, the spectral bands passing the significance test (P<0.01) were extracted as the SSC characteristic bands. A random forest model was then constructed to quantitatively estimate the SSC. The results indicate that: 1) There was a general consistency in the spectral curve shapes of soils with different salinization levels. There was a great variation in the patterns of reflectance with the wavelength and SSC. As the wavelength increased from the short (350 nm) to the long (2 450 nm), the reflectance increased in the range of 350~800 nm and then remained relatively constant in the range of 800~2 130 nm. While an absorption peak appeared around 2 130 nm corresponding to Na2SO4. After that, the reflectance started to decrease. Furthermore, the reflectance increased with the SSC in the entire range of wavelengths. There was also a positive correlation between them. But there was a decreasing trend at extreme levels of high salinization. 2) First-order differential transformation, CWT, DWT, as well as first-order differential combined with CWT and first-order differential combined with DWT, all significantly enhanced the sensitivity of the spectra to the SSC. Specifically, the DWT was performed the best among them. The highest correlation coefficient between the processed spectra and SSC reached −0.621 (P<0.01), indicating a negative correlation. The performance of the rest was ranked in the descending order of the first-order differential combined with CWT, first-order differential combined with DWT, CWT, and first-order differential transformation. 3) The first-order differential transformation, CWT, DWT, first-order differential combined with CWT, and first-order differential combined with DWT all improved the accuracy and stability of the model. The combination of first-order differential with CWT significantly enhanced the accuracy of SSC estimation. The optimal model used the full-band spectral of (1/R)'_CWT_28, with the coefficients of determination (R2) of 0.821 and 0.715 for the training and validation sets, respectively, the root mean square errors (RMSE) of 16.049 and 17.467 g/kg, respectively, and the relative predictive determinant (RPD) values of 1.65 and 1.48, respectively. As such, the effective SSC estimation was achieved in the study area. Therefore, the mathematical transformations with the wavelet transformation were better utilized to process the smoothed original spectral data. The sensitivity of spectra to SSC was significantly enhanced after processing. Compared with the mathematical or wavelet transformations alone, the combined approach was more effective in quantitatively monitoring the SSC. The finding can also provide a valuable reference to alleviate the soil structure deterioration for agricultural production and environmental protection

       

    /

    返回文章
    返回