基于Sentinel-2A的孙吴地区土壤有机质反演研究

doi:10.11720/wtyht.2022.0038

基于Sentinel-2A的孙吴地区土壤有机质反演研究

陈超群^,¹^,²^,³, 戴慧敏¹^,²^,³, 冯雨林¹, 杨泽¹^,²^,³, 杨佳佳^,¹

1.中国地质调查局沈阳地质调查中心,辽宁沈阳 110034

2.自然资源部黑土地演化与生态效应重点实验室,辽宁沈阳 110034

3.辽宁省黑土地演化与生态效应重点实验室,辽宁沈阳 110034

Sentinel-2A based inversion of the organic matter content of soil in the Sunwu area

CHEN Chao-Qun^,¹^,²^,³, DAI Hui-Min¹^,²^,³, FENG Yu-Lin¹, YANG Ze¹^,²^,³, YANG Jia-Jia^,¹

1. Shenyang Center of China Geological Survey, Shenyang 110034, China

2. Key Laboratory of Black Soil Evolution and Ecological Effect, Ministry of Natural Resources, Shenyang 110034, China

3. Key Laboratory of Black Soil Evolution and Ecological Effect, Liaoning Province, Shenyang 110034, China

通讯作者: 杨佳佳(1984-),男,博士,高级工程师,主要从事遥感应用研究工作。Email:haixianxiaomei@163.com

责任编辑: 蒋实

收稿日期: 2022-01-25 修回日期: 2022-06-27

基金资助:

中国地质调查局项目“东北黑土地1:25万土地质量地球化学调查”(121201007000161312)
“兴凯湖平原及松辽平原西部土地质量地球化学调查”(DD20190520)

Received: 2022-01-25 Revised: 2022-06-27

作者简介 About authors

陈超群(1996-),女,硕士研究生,主要研究方向为生态环境遥感与地理信息系统。Email:522110156@qq.com

摘要

利用Sentinel-2A多光谱遥感影像,结合实测土壤信息,对黑龙江省孙吴县黑土区土壤有机质含量进行反演研究。对影像进行预处理后,通过相关分析和随机森林(RF)选取特征波段,采用偏最小二乘法和BP神经网络构建土壤有机质含量多光谱模型反演红旗林场土壤有机质含量。研究表明:相关性选取的倒数对数一阶微分反射率波段和RF选择的组合波段能够有效提高土壤反演精度,组合波段的RF-BP神经网络模型反演效果最佳,R²=0.724 5,RMSE=1.312 7%。本次研究可为实现土壤有机质动态监测提供技术支持和参考。

关键词： 黑土; 有机质; Sentinel-2A; 随机森林; BP神经网络

Abstract

This study conducted the inversion of the organic matter content in the soil of the black soil area in Sunwu County, Heilongjiang Province using the Sentinel-2A multispectral remote sensing images and the surveyed soil data. After preprocessing the images, the characteristic bands were selected through correlation analysis and using the random forest (RF) method. Subsequently, a multispectral inversion model for the organic matter content of the soil was built using the partial least square method and the BP neural network, and the inversion of the organic matter content of the soil in the Hongqi Forest Farm was conducted. According to the obtained results, the bands selected based on the reciprocal of the logarithm of the first-order differential of reflectance through the correlation analysis and the combined bands selected using the RF method can effectively improve the inversion precision of the organic matter content in the soil, and the RF-BP neural network model for the combined bands yielded the optimal inversion performance (R²=0.7245 and RMSE=1.3127%). The results of this study will provide technical support and reference for the dynamic monitoring of the organic matter content in soils.

Keywords： black soil; organic matter; Sentinel-2A; random forest; BP neural network

PDF (1863KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

陈超群, 戴慧敏, 冯雨林, 杨泽, 杨佳佳. 基于Sentinel-2A的孙吴地区土壤有机质反演研究[J]. 物探与化探, 2022, 46(5): 1141-1148 doi:10.11720/wtyht.2022.0038

CHEN Chao-Qun, DAI Hui-Min, FENG Yu-Lin, YANG Ze, YANG Jia-Jia. Sentinel-2A based inversion of the organic matter content of soil in the Sunwu area[J]. Geophysical and Geochemical Exploration, 2022, 46(5): 1141-1148 doi:10.11720/wtyht.2022.0038

0 引言

土壤有机质是土壤质量的重要参数,可为农作物提供各类养分,同时对元素表生地球化学特征有重要影响。黑土作为珍贵的土壤资源,其有机质含量是反映土壤质量的重要指标参数^[1-2]。近年来,随着黑土逐渐退化,土壤中有机质成分明显减少,估算黑土有机质含量,扭转含量下降趋势,是黑土地保护的重要举措^[3]。传统的土壤有机质监测主要是通过对监控区进行大量野外土壤样品采集和室内化学实验分析进行反演,这种方法周期较长,费时费力,精度受样品密度控制,难以满足现代农业快速发展的需求^[4]。随着遥感技术日益成熟,通过有机质含量的光谱差异来测定土壤有机质的含量已成为一种有效手段。

土壤有机质含量遥感反演主要包含两个研究方向:光谱信息的处理与选择和反演模型的构建。常采用的光谱处理方法有倒数、对数、去包络线变换等,但选取的有机质特征波段因影像数据源不同而有所差异。屈冉等^[5]选取Landsat TM 影像反演广西壮族自治区富川县的有机质含量,认为土壤有机质含量与Landsat TM 波段5和波段7 的DN值相关性最高。陈德宝等^[6]借助Landsat 8遥感影像对农安县黑土区有机质进行建模反演,表明短波红外B6波段反射率所建模型拟合效果最好。陈思明等^[7]对Landsat 7 土壤光谱进行线性波谱分离重建,认为重建光谱能显著增强与土壤有机质含量的相关性,提高土壤有机质反演精度。在以往研究中,线性回归和偏最小二乘回归模型(PLSR)常被用于土壤有机质含量反演。Dhawale 等^[8]结合土壤样品有机质含量和相应的土壤反射率,选用PLSR建模,均方根误差不超过2.24%。马驰^[9]对比Sentinel-2A遥感影像不同波段组合的多元回归模型,R²均大于0.7。目前针对有机质敏感波段的选择主要采用Pearson相关分析法,反演模型也多选择线性拟合。本次研究借助Sentinel-2A遥感影像,结合黑河市孙吴县实测土壤有机质含量,通过Peason 相关分析和随机森林(RF)选择不同特征波段作为模型输入量,采用PLSR和BP 神经网络建模,以期研究土壤表层有机质含量与遥感影像关系,并实现地面黑土区红旗林场的土壤有机质高精度快速反演。

1 研究区概况

研究区孙吴县地处黑龙江省黑河市中部,位于东经126°39'35″~128°1'6″,北纬48°59'00″~49°41'55″(图1)。东部紧靠逊克县,西边为嫩江县,南侧与五大连池市相挨,北方为黑河市爱辉区,总面积4 318.9 km²。孙吴县海拔110~755 m,属于低山丘陵区,地势总体呈西南高东北低趋势。地貌分界清晰,从西到东分别为低山沟谷区、丘陵河谷地区和沿江平原。土壤类型以暗棕壤和草甸黑土为主^[10]。气候属于寒温带大陆性季风气候,年均气温-0.6 ℃,年均降雨约550 mm,冻结期较长,无霜期短^[11]。本文选择孙吴县红旗林场地区进行遥感反演,红旗林场位于孙吴县西北方向,范围为东经126°41'25″~127°14'34″,北纬49°16'32″~49°30'58″,界内发育孙吴县最高山峰松木山。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 孙吴县遥感影像(a)及红旗林场位置(b)

Fig.1 Remote sensing image of Sunwu County(a)and the location of Hongqi Forest Farm(b)

2 数据采集与处理

2.1 土壤采集与有机质测定

按照《土地质量地球化学评价规范》(DZ/T 0295—2016)采样要求,在孙吴县采集土壤时去除表面枯枝落叶等杂物,用刻槽法垂直采集地表至20 cm深土样,保证上下均匀采集,并去除动、植物残留体、砾石、肥料团块等。土壤有机质含量采用硫酸—重铬酸钾法测定。共计采集806个样品,其中564个土样作为建模集,242个样品为测试集,统计信息如表1所示。

表1 土壤样品中有机质含量统计信息

Table 1 Statistical information of organic matter content in soil samples

	个数	最小值/%	最大值/%	均值/%	标准差/%
建模集	564	0.8620	11.8266	5.7226	2.0316
测试集	242	1.1896	11.9128	5.8393	2.1074

新窗口打开| 下载CSV

2.2 遥感数据获取与处理

选取研究区内2018年11月7日裸土无雪时期的Sentinel-2A影像,云覆盖0%。影像的预处理包括几何校正、大气校正、图像镶嵌及图像剪裁等操作。所选影像为Level-1C 上层大气反射产品,已经过系统几何精校正处理,其精度在一个像元内,满足研究需求。借助SNAP软件中Sen2cor280工具箱实现大气校正,校正后丢失卷云波段B10。为提高土壤有机质与光谱反射率(R)相关性,对遥感影像进行倒数(1/R)、对数(lgR)、幂函数(R^a)、一阶微分(FDR)、二阶微分(SDR)及倒数对数一阶微分(FDLR)处理。

3 算法原理

3.1 Pearson相关分析

为获取土壤有机质光谱响应波段,研究中采用Pearson相关判断反射率与有机质含量之间的线性相关性。公式为:

(1)

r = \frac{\sum_{i = 1}^{n} (X_{i} - \bar{X}) (Y_{i} - \bar{Y})}{\sqrt{{(X_{i} - \bar{X})}^{2}} \sqrt{{(Y_{i} - \bar{Y})}^{2}}}

式中:r表示相关系数; $X_{i}$ 为样品 $i$ 的波段的反射率; $Y_{i}$ 为样品 $i 的$ 有机质含量; $\bar{X} 为反射率均值$ ; $\bar{Y}$ 为有机质含量均值。为防止出现虚假相关,引入显著水平sig进行判断。将通过显著性检验的波段及其变换作为土壤有机质的特征波段。

3.2 随机森林

2001年Breiman提出随机森林(random forest,RF)算法,主要优势体现在处理多维数据集时无需降维,最大程度的保留数据集原始信息^[12-13]。本次研究中以X.IncMSE指标为重要性选择依据,X.IncMSE值越大,表明该波段信息越能反应有机质含量。

3.3 偏最小二乘回归模型

偏最小二乘回归分析(partial least square regression,PLSR)常用于遥感光谱反演建模,优势在于建模过程中集中了主成分分析、典型相关分析和线性回归分析方法的优点^[14]。其建模思路为:设自变量 $X = {[x_{1}, x_{2}, \dots, x_{p}]}_{n \times p}$ ,因变量为 $Y = {[y]}_{n \times 1}$ ,其中X为波谱曲线反射率及其变换形式,Y为土壤成分含量,n为采集样本数,p为特征波段数。从自变量信息中提取最大变异信息成分 $μ_{1}$ ,且与因变量呈最大相关性。在提取第一主成分 $μ_{1}$ 后,建立Y与 $μ_{1}$ 的回归方程,若模型精度满意,则算法停止;否则继续利用X和Y被成分解释后的残余信息进行多次迭代提取,直到回归方程达到满意精度。

为了检验反演模型的精度及稳定性,借助决定系数 $R^{2}$ 和均方根误差RMSE作为模型评价指标,公式如下:

(2)

R^{2} = \sum_{i = 1}^{n} {(\overset{︿}{y_{i}} - \bar{y})}^{2} / \sum_{i = 1}^{n} {(y_{i} - \bar{y})}^{2}

(3)

R M S E = \sqrt{\frac{\overset{n}{\sum_{i = 1}} (\overset{︿}{y_{t}} - y_{i})^{2}}{n}}

式中:n为样本个数; $y_{i}$ 为土壤有机质含量真实值; $\bar{y}$ 为有机质含量平均值; ${\overset{︿}{y}}_{i}$ 为有机质含量预测值。R²越大,表明模型越稳定,拟合效果越好。RMSE值越小,表明模型的估测能力越稳定^[15]。

3.4 BP神经网络

BP神经网络(误差反向传播)是人工神经网络的一种,由输入层—隐含层—输出层组成 ^[16]。学习过程包含信号的正向传播与误差的反向传播,正向传播时,信号从输入层传入,经隐含层逐层处理后,最后达到输出层。若输出层的实际输出与期望输出不符,则转向误差的反向传播阶段^[17]。BP神经网络具有较强的非线性处理能力和自适应特点,能够较好地拟合光谱反射率与土壤有机质含量的关系。通过经验公式(4)确定隐含层的节点数取值范围,结合训练结果的精度选择最佳隐含层的节点数:

(4)

N = \sqrt{n + m} + k

式中:n为输入层节点数;m为输出层节点数;k为1~10之间的常数;N为隐藏层节点数。

4 黑土区土壤有机质反演

4.1 特征波段选择

4.1.1 相关性波段选择

在SPSS 26平台下计算土壤有机质含量与Sentinel-2A遥感影像反射率及其变换间的相关性。如图2所示,有机质含量与光谱反射率呈现负相关,但相关性不高,各波段不同数学变换的相关系数最高值绝大多数出现在FDLR变换,表明该预处理方法能有效提高Sentinel-2A反射率与土壤有机质的相关性。将不同数学变换中通过显著性检验的波段作为相应反演的特征波段,同时组合各波段中最高相关系数的变换形式作为一种响应波段参考,其中B8和B11无通过显著性检验的变换,故不予讨论(表2)。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 波段反射率及其变换与土壤有机质含量相关性

Fig.2 Correlation between band reflectivity and transformations and soil organic matter content

表2 相关性选取特征波段

Table 2 Feature bands selected by correlation analysis

数学变换	波段
R	B1、B2、B3、B4、B5、B6、B7、B8、B9、B12
1/R	B1、B2
lgR	B1、B2、B3、B4、B5、B12
R^a	B4、B5、B6、B7、B8、 B9、B12
FDR	B1、B2、B3、B4、B5、B6
SDR	B1、B4、B5
FDLR	B1、B2、B3、B4、B5、B6、B8、B12
组合波段	1/B1、FDL(B2)、lg (B3)、B4、FDL(B5)、FDL(B6)、B7、FDL(B8)、(B9)^a、FDL(B12)

新窗口打开| 下载CSV

4.1.2 RF重要性选择

采用R语言的randomForest包实现土壤有机质特征波段选取,其中默认生成500棵决策树,并进行5次重复十折交叉验证,结合最简原则选择不同光谱变换下的特征波段。以对数变换为例,图3交叉验证曲线展示了模型误差与用于拟合的自变量数量的关系,当波段数为6时,误差下降幅度基本保持不变,结合简约性原则,选择重要程度值从大到小排序前6的波段作为有机质反演建模的输入参数,实验中R、1/R、R^a、FDR、SDR以及FDLR均需要6个重要变量表示土壤有机质含量。为提高变量表达精度,将所有波段的变换作为RF的因变量,并需要X.IncMSE值前26个波段变换作为建模输入参数集,以精准表达有机质含量信息(图4)。RF重要波段选取结果如表3所示。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 对数变化交叉验证曲线

Fig.3 Cross validation curve of lgR

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 所有波段交叉验证曲线

Fig.4 Cross validation curve of all bands

表3 RF重要波段

Table 3 Important bands of RF

数学变换	波段
R	B1、B2、B3、B5、B7、B8A
1/R	B1、B2、B3、B5、B7、B8A
lgR	B1、B2、B3、B5、B6、B7
R^a	B2、B3、B5、B6、B7、 B8A
FDR	B1、B2、B3、B4、B6、B12
SDR	B1、B3、B4、B5、B11、B12
FDLR	B1、B2、B3、B4、B8、B12
组合波段	B1、B7、B11、1/B1、1/B2、1/B3、1/B11、1/B1、lg(B2)、FDL (B1)、FDL (B2)、FDL (B3)、FDL (B5)、FDL (B6)、FDL (B8)、FD (B1)、FD (B5)、FD (B6)、FD (B12)、SD(B1)、SD(B2)、SD(B3)、SD(B5)、SD(B8)、SD(B11)、SD(B12)

新窗口打开| 下载CSV

4.2 土壤有机质反演

4.2.1 PLSR模型反演

将Pearson相关分析(表2)和RF(表3)提取的特征波段作为自变量,土壤有机质含量作为因变量,建立有机质含量PLSR反演模型,如表4所示,结果显示相关-PLSR模型和RF-PLSR模型反演精度结果相近,FDLR变换和组合波段都能有效提高模型反演精度,其中针对传统的相关-PLSR模型,FDLR光谱变换的拟合程度最好,RF-PLSR模型中组合波段的效果更为显著。但PLSR模型下建模集和测试集的决定系数R²均未超过0.1。

表4 基于PLSR模型的土壤有机质反演

Table 4 Inversion of soil organic matter by PLSR

	数学变换	拟合模型	建模集		测试集
	数学变换	拟合模型	R²	RMSE/%	R²	RMSE/%
相关性	R	$y = 6.3252 - 0.1679 x_{1} 1 - 0.3836 x_{2} - 0.4512 x_{3} - 0.6285 x_{4} - 0.6827 x_{5} - 0.6340 x_{6} + 0.5943 x_{7} - 0.5943 x_{8} - 0.5199 x_{9} - 0.5565 x_{10}$	0.0275	2.0263	0.0155	2.0925
	1/R	$y = 5.2818 + 0.0001 x_{1} + 0.000017 x_{2}$	0.0439	1.9847	0.0657	2.0375
	lgR	$y = 5.4429 - 0.5202 x_{1} - 0.2827 x_{2} + 1.1558 x_{3} - 1.5976 x_{4} - 0.5254 x_{5} + 3.3470 x_{6}$	0.0525	1.9757	0.0609	2.0406
	R^a	$y = 6.0894 + 48.0295 x_{1} - 110.9565 x_{2} + 21.8818 x_{3} + 78.5713 x_{4} - 33.5780 x_{5} - 10.4766 x_{6} - 5.2279 x_{7}$	0.0159	2.0136	0.0121	2.0917
	FDR	y=5.8376-23.9683x₁+67.1411x₂-63.9701x₃+44.2371x₄	0.0223	2.0071	0.0409	2.0632
	SDR	$y = 6.1351 + 4.19 x_{1} + 14.4844 x_{2} + 59.2323 x_{3}$	0.0200	2.0093	0.0238	2.0800
	FDLR	$y = 5.0094 + 1.7352 x_{1} - 6.2089 x_{2} + 6.6447 x_{3} - 5.4056 x_{4} + 3.6822 x_{5} - 11.4086 x_{6} - 3.5433 x_{7} - 3.0895 x_{8}$	0.0524	1.9779	0.071	2.0366
	组合	$y = 8.8137 - 0.0162 x_{1} + 1.7147 x_{2} + 0.4484 x_{3} - 0.1764 x_{4} - 1.1268 x_{5} - 0.4075 x_{6} - 0.5233 x_{7} - 1.5618 x_{8} - 0.8610 x_{9} - 1.2728 x_{10}$	0.0434	1.9852	0.0400	1.9933
随机森林	R	$y = 5.77391 - 8.6669 x_{1} - 0.6491 x_{2} + 9.5163 x_{3} - 26.4628 x_{4} + 8.9020 x_{5} + 6.6144 x_{6}$	0.0203	2.0318	0.0035	2.0539
	1/R	$y = 5.80955 + 0.0002 x_{1} + 0.0001 x_{2} + 0.0009 x_{3} - 0.1832 x_{4} - 0.0226 x_{5} + 0.2048 x_{6}$	0.0234	2.2798	0.09	2.0390
	lgR	$y = 4.40343 - 0.4786 x_{1} - 0.0567 x_{2} + 0.0415 x_{3} + 0.0548 x_{4} + 0.0692 x_{5} + 0.0786 x_{6}$	0.0480	2.0029	0.0408	2.0124
	R^a	$y = 5.95664 - 0.4164 x_{1} - 0.8538 x_{2} - 2.1607 x_{3} - 2.2302 x_{4} - 2.3068 x_{5} - 2.6665 x_{6}$	0.0051	2.0475	0.0023	2.0525
	FDR	$y = 5.65109 - 33.6899 x_{1} + 85.8885 x_{2} - 75.5978 x_{3} - 13.8874 x_{4} + 28.8994 x_{5} - 6.2782 x_{6}$	0.0365	2.0149	0.0207	2.0356
	SDR	$y = 5.55202 + 13.2022 x_{1} - 29.3378 x_{2} + 42.5214 x_{3} + 53.8592 x_{4} - 3.5595 x_{5} - 2.0586 x_{6}$	0.0331	2.0496	0.0200	2.0520
	FDLR	$y = 4.69522 - 0.0386 x_{1} - 0.7984 x_{2} - 0.6553 x_{3} - 0.1570 x_{4} - 0.0148 x_{5} + 0.0701 x_{6}$	0.0463	2.0046	0.0418	2.0118
	组合	$\begin{array}{l} y = 1.2967 x_{1} + 0.2947 x_{2} + 2.6967 x_{3} + 0.5528 x_{4} + 0.8238 x_{5} + 0.7684 x_{6} - 0.2272 x_{7} - 0.6389 x_{8} + 1.0292 x_{9} + 1.3686 x_{10} - 1.8264 \\ x_{11} + 0.7614 x_{12} - 0.1926 x_{13} + 1.2489 x_{14} - 0.1014 x_{15} - 0.3382 x_{16} + 0.9709 x_{17} - 0.3017 x_{18} + 1.5403 x_{19} + 2.7079 x_{20} + 0.0221 x_{21} - 3.4156 x_{22} + 3.8703 x_{23} + 4.6421 x_{24} - 0.5339 x_{25} - 1.746 x_{26} - 1.7177 \end{array}$	0.0760	1.9518	0.0363	2.0183

新窗口打开| 下载CSV

4.2.2 BP神经网络模型反演

结合式(4),根据相关性和RF重要程度选取获得的特征波段数,确定BP神经网络模型中隐含层层数。网络的训练函数为Trainlm,输入层和输出层传递函数分别为Tansig和Purelin,表5为相关分析和RF两种特征响应波段模拟结果,对比PLSR拟合结果(表4),非线性拟合的多光谱遥感影像反射率与土壤有机质含量模型精度能得到显著提高。由于多光谱遥感的光谱分辨率较低,光谱包含的土壤信息较为复杂,因此无法类比高光谱土壤有机质遥感反演,线性回归拟合模型不能有效提取影像上土壤有机质含量信息^[18⇓-20]。相关性提取波段与RF提取的重要波段在进行BP神经网络建模时,建模集和测试集的R²主要集中在0.2~0.5,RMSE集中在1.3%~1.4%。相关-BP神经网络模型中FDLR建模拟合程度最高,建模集R²为0.623 7,RMSE为1.354 8%,测试集 R²为0.444 6,RMSE为1.266 4%。RF-BP神经网络模型中组合波段建模拟合程度最高,建模集R²为0.724 5,RMSE为1.312 7%,测试集 R²为0.541 8,RMSE为1.372 2%。

表5 基于BP神经网络模型的土壤有机质反演

Table 5 Inversion of soil organic matter by BP neural network

	数学变换	隐藏层个数	建模集		测试集
	数学变换	隐藏层个数	R²	RMSE/%	R²	RMSE/%
相关性	R	5	0.3635	1.3998	0.2711	1.4472
	1/R	11	0.2816	1.4057	0.2291	1.4283
	lgR	6	0.3392	1.3967	0.2818	1.4207
	R^a	11	0.2726	1.4041	0.2388	1.4555
	FDR	10	0.2074	1.4149	0.1697	1.4368
	SDR	9	0.2005	1.4106	0.1977	1.4388
	FDLR	6	0.6237	1.3548	0.4446	1.2664
	组合	16	0.5637	1.3548	0.4305	1.3659
随机森林	R	7	0.2906	1.4000	0.2601	1.4362
	1/R	7	0.2603	1.4068	0.2241	1.4262
	lgR	14	0.2883	1.4068	0.2664	1.4231
	R^a	14	0.1663	1.4152	0.2499	1.4200
	FDR	7	0.3751	1.3980	0.2860	1.4090
	SDR	11	0.2783	1.4107	0.2057	1.4520
	FDLR	13	0.4544	1.3750	0.3653	1.2420
	组合	16	0.7245	1.3127	0.5418	1.3722

新窗口打开| 下载CSV

结合表4中对比R、1/R等波段变换的不同提取方法,Pearson相关分析中选择的FDLR反射率变换在线性回归和非线性回归中都展现较高的拟合效果。因为Pearson相关性分析获取的是有机质含量与反射率间简单直线性相关的方向和密切程度,因此光谱处理的程度直接决定了与有机质含量的相关性,进而影响了模型反演精度。FDLR变换中光谱倒数对数计算可以有效放大相似光谱间的差别,再经过一阶微分处理后消除部分线性的背景,同时降低噪声光谱对目标光谱的影响程度。但相关性选择的组合波段包含了反射率变量处理的不同级别,当加入相关性低的波段变换,很有可能引入了土壤其他成分的特征信息,导致有机质反演精度降低。而在RF重要程度选择中存在误差验证,结合每次选择的特征集计算袋外误差率,最后选择袋外误差率最低的特征集作为回归模型的输入集。筛选出的波段可能与有机质含量相关性低,但叠加其他波段光谱特征反而提高了有机质估测精度。因此RF波段选取时,组合所有光谱反射率及变换信息后筛选出的特征响应波段更能充分反映有机质含量信息,建模精度显著提升。

4.3 红旗林场有机质空间分布

对比8种影像变换的相关-PLSR、相关-BP神经网络、RF-PLSR、RF-BP神经网络建模情况,选择模型拟合度最高、稳定性最好的RF-BP神经网络模型作为多光谱遥感数据的土壤有机质反演模型,并预测孙吴县红旗林场的有机质含量分布(图5)。红旗林场的有机质范围大致在0.1%~18.8%,主要集中在3.742 6%~12.455 2%,平均含量为7.939 9%,呈现出中间高、四周低的分布趋势。该地区土壤有机质地球化学分析结果显示,红旗林场有机质平均含量为8.51%,分布趋势为中西部含量偏高,向北侧逐渐降低,该结果与遥感反演程度几乎吻合。但由于在林场附近采集土壤样品较少,因此地球化学实测结果分布较为粗糙,而遥感反演获取的土壤有机质含量更为细致。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 红旗林场土壤有机质遥感反演和地球化学对比

Fig.5 Distribution of soil organic matter in Hongqi Forest Farm

5 结论

对Sentinel-2A多光谱遥感影像反射率进行1/R、lgR、R^a、FDR、SDR及FDLR变换,结合不同模型实现了土壤有机质的反演,取得很好应用成效。主要得出以下结论:

1)通过相关性分析法建模时反射率的FDLR变换模型拟合程度最好,而采用RF算法筛选的组合波段在反演时能有效提高土壤有机质含量建模精度。

2)多光谱遥感影像光谱分辨率较低,因此线性拟合模型无法准确估测土壤有机质含量,需要非线性模型以实现光谱信息与有机质含量的有效拟合。

3)对比不同遥感影像预处理操作下的相关-PLSR、相关-BP神经网络、RF-PLSR、RF-BP神经网络建模情况,RF-BP神经网络模型反演土壤有机质含量拟合程度最高,建模集R²为0.724 5,RMSE为1.312 7%,测试集 R²为0.541 8,RMSE为1.372 2%。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

孔牧, 杨少平.

森林沼泽景观区有机质对元素表生地球化学特征的影响机制

[J]. 物探与化探, 2008, 32(1):31-32,74.