不同降噪方式下基于高分五号影像的土壤有机质反演

通过遥感技术预测土壤有机质（Soil Organic Matter，SOM）的空间分布是精准农业和土地资源管理研究的重要内容，与粮食安全及环境监测密切相关。该研究主要研究采用高分五号（GF-5）反演土壤有机质的最佳降噪方式。以黑龙江省典型黑土区明水县为研究对象，获取GF-5高光谱遥感影像，对影像进行不同降噪处理，包括奇异值分解（Singular Value Decomposition，SVD），离散小波变换（Discrete Wavelet Transform，DWT）及中值滤波（Median Filtering，MF）降噪。然而，分别结合二维光谱指数，应用随机森林（Random Forest，RF）方法预测不同降噪方式的SOM含量。结果表明：1）所选择的不同降噪方法中，与SOM含量的相关性由高到低依次为DWT、SVD、MF，其中，基于MF降噪后的光谱反射率与SOM含量相关性低于原始反射率与SOM含量的相关性。2）基于降噪方式下的光谱曲线更加平滑，且不同有机质含量对应的光谱曲线形状相似。3）采用DWT降噪方式，基于影像波段和光谱指数，以RF为预测模型的SOM最优反演模型精度R2为0.69，均方根误差为2.26%。研究成果可为利用高光谱卫星数据实现大尺度范围内SOM的数字土壤制图提供参考，为实时定量监测土壤肥力变化提供依据。

Inversion of soil organic matter based on GF-5 images under different noise reduction methods

Predicting the spatial distribution of Soil Organic Matter(SOM)through remote sensing images is an important part of precision agriculture and land resource management research.The emergence of Gaofen5(GF-5)has facilitated quantitative soil monitoring based on hyperspectral remote sensing.As the first full spectrum hyperspectral satellite for comprehensive observations of the atmosphere and land,GF-5 has a high spatial resolution of 30 m,including a total of 330 bands,and the coverage is consistent with the range obtained by spectro-radiometerFieldSpec3(Analytical Spectral Devices,Longmont,USA).However,there is no clear and appropriate noise reduction method in the relevant soil research.Therefore,this study was to explore the best noise reduction method for GF-5 application in soil organic matter inversion.The experimental area was Mingshui County(124°18′-125°21′E,46°44′-47°29′N),a typical black soil region of Heilongjiang Province.A total of 38 sampling points were selected at the depth of 0-20 cm in the study area and the spatial coordinates of sampling points were recorded by GPS.The non-cloud GF-5 hyperspectral remote sensing image was obtained from the website of China Resources Satellite Application Center.According to our prior knowledge,radiometric calibration and atmospheric correction pre-processing of hyperspectral remote sensing images in ENVI5.1 were performed.Singular Value Decomposition(SVD),Discrete Wavelet Transformation(DWT)and Median Filtering(MF)noise reduction method were used.The two-dimensional spectral indexes were calculated under different noise reduction methods.The Principal Component Analysis(PCA)method was used to decrease the dimensionality of the inputs,and the Random Forest(RF)method was used to predict the SOM content.A total of 26 data were used for model establishment and the rest 12 data were used for model validation.The SOM content of these samples were 3.45%-5.53%with mean of 4.44%,standard deviation of 4.28%and coefficient of variance of 9.6.The results showed that the best decomposition layer of DWT was 1 since the correlation coefficient between reflectance and SOM was the highest and the wavelength locations with high correlation were consistent with those in the original reflectance curves.Among different noise reduction methods,the correlation between the spectral curve and the SOM content was the highest for DWT,followed by SVD and MF.The small"burr"in the soil reflection spectrum curve was removed by using different noise reduction methods and the shapes of the spectral curves under different SOM contents were consistent.The principal component analysis results showed that the cumulative contribution rates for data experiencing noise reduction by SVD,DWT and MF methods were all about 95%.The DWT method combined with RF prediction model achieved the highest accuracy of SOM estimation for validation dataset(R^2:0.69,RMSE:2.26%,RPD:1.80)and the noise reduction effect was the most significant.The SOM spatial distribution under the DWT noise reduction model was mapped and the result was consistent with the real condition.The SOM content was relatively high and its distribution was uniform in the central areas of Mingshui County.However,it was low in the southeast region because of severe soil erosion.The research results can provide an effective way of the realization of remote sensing prediction of SOM at a large scale by using GF-5 hyperspectral satellite data,expanding the application scope of satellite images.