专题一:土壤有机质含量高光谱预测模型及其差异性研究

合集下载

基于高光谱遥感的土壤有机质反演模型研究进展

基于高光谱遥感的土壤有机质反演模型研究进展

设施农业2023-0727农业工程技术(设施农业)基于高光谱遥感的土壤有机质反演模型研究进展*摘要:利用高光谱遥感反演土壤有机质与传统的化学分析方式有很大的优势,促进现代农业的发展。

国内外利用高光谱反演土壤有机质含量模型进行了大量的研究,从最初的线性回归模型发展到如今较为复杂的非线性模型,该文总结了常见的高光谱对有机质含量建模方法与模型特点。

通过分析总结前人的研究发现,目前土壤有机质含量模型研究进展主要有运用的模型复杂化、反演模型的评价方法多样化、多种建模方法综合应用及深度学习大数据分析新技术更加普遍。

关键词: 高光谱遥感;土壤有机质含量;反演模型土壤有机质是衡量土壤肥力的重要指标,是维持农业健康发展、土地质量评价的重要考虑因素[1]。

传统的土壤有机质化学测定方法试剂消耗量大、具有一定危险性、对操作人员有专业性要求等原因,难以做到迅速、大规模和周期性地土壤有机质含量测定[2]。

而高光谱遥感可以快速地进行大范围、周期性、快速的特点受到人们的关注。

利用高光谱反演土壤有机质含量最重要的是建立土壤高光谱信息与有机质含量之间的联系。

目前利用高光谱遥感反演土壤有机质含量的模型中,大体可以分为两大类。

一类是线性模型。

例如:偏最小二乘法,多元线性回归模型;另一类是非线性模型。

例如:反向神经网络,随机森林,支持向量机等。

根据目前的研究进展,总结了现阶段利用高光谱遥感反演土壤有机质含量方法和模型的特点。

1 土壤有机质的高光谱反演模型发展历程1965年,Bowers 等[3]开始研究土壤与光谱反射之间的关系时就发现土壤有机质含量与光谱反射之间存在负相关的关系。

21世纪前,主要是研究土壤有机质与高光谱反射率之间的定性关系上的研究,处于探索阶段,研究样本数量较少,使用的模型主要以多元线性回归,偏最小二乘回归等线性方法为主[4]。

近年来随着研究的深入与研究工具的进步,更多的模型方法,研究方式的韦森超1,梁立恒1,董炜华1**,李晓强1,邓守奇1,2(1.长春师范大学地理科学学院,吉林 长春 130123;2.长春文理高中,吉林 长春 130062)改进,比如机器学习与深度学习[5]。

基于不同模型的土壤有机质含量高光谱反演比较分析

基于不同模型的土壤有机质含量高光谱反演比较分析
第3 3 卷, 第1 期
2 0 1 3 年 1月








Vo 1 . 3 3 , No . 1 , p p 1 9 6 — 2 0 0
S p e c t r o s c o p y a n d S p e c t r a l An a l y s i s
4 .新疆大学资源与环境科学学院 , 新疆 乌鲁木齐
8 3 0 0 4 6


以新疆奇 台县 为研究 区域 ,选取 该县 4 O个土壤样本 , 采用多元线性 逐步 回归 法和人工神经 网络法
两种方法分别建立了土壤 有机质含 量的反演 模型 , 并 对模 型进行 了检验 。 结 果发 现 : 不 同模型 的精度值 各 异, 其拟合效果从高到低依次为人工神经网络( ANNs ) 集 成模 型>单个人 工神经 网络 ( ANNs ) 模型 >多元逐 步 回归 ( ML S R) 模 型。人工神经网络的线性和非线性 逼近能力 较强 ,而其集成模型作为提高反演模 型精度 的 重要手段 , 相关 系数高达 0 . 9 3 8 , 均方根误差和总均方根 误差 最小 , 分别仅 为 2 . 1 3和 1 . 4 0 4 , 对土 壤有机质 含量 的预测能力与实测光谱非常接近 , 分析结果达到 了较实用的预测精度 , 为最优拟合模型 。 关键词 奇 台;土壤有机质 ;高光谱 反演 分析 ; 多元逐步 回归 ; 人工
基 于不 同模 型 的 土壤 有机 质含 量 高 光谱 反演 比较分 析
栾福 明 , 张小雷 , 熊黑钢。 一, 张 芳 , 王 芳
8 3 0 0 1 1
1 .中国科学院新疆生态与地理研究所 , 新疆 乌鲁木齐 2 .中国科学 院大学 , 北京 1 0 0 0 4 9 1 0 0 0 8 3 3 .北 京联合大学应用文理学院 , 北京

利用光谱数据进行土壤有机质评价

利用光谱数据进行土壤有机质评价

利用光谱数据进行土壤有机质评价随着世界人口的增长和农业生产方式的改变,土壤有机质的质量越来越受到人们的关注。

土壤有机质是一种关键的土壤物质,它对土壤肥力和作物生长有着非常重要的影响。

其中,土壤的颜色是反映土壤有机质含量的重要指标之一。

然而,传统的土壤有机质检测方法需要取样,样品处理等繁琐的操作,需要耗费大量的时间和人力,而光谱数据则提供了一种快速、准确的土壤有机质评价方法。

本文将从光谱数据的基本原理、光谱仪使用原理以及光谱数据处理方法等方面来阐述如何利用光谱数据进行土壤有机质评价。

一、光谱数据的基本原理光谱数据是指一段特定波长内的光传输信息,包括红外线、可见光和紫外线等,其中红外线光谱数据在土壤有机质评价中具有非常重要的地位。

土壤有机质含量对土壤红外线光谱信号的吸收有很大的影响,因此,利用光谱仪测量土壤的红外光谱信号,可以获取土壤有机质的含量信息。

二、光谱仪的使用原理光谱仪是一种专门用于测量光谱数据的仪器设备,主要由光源、样品室、检测器等几个主要组成部分构成。

在进行土壤有机质评价时,需要将测量所需的光源辐射到样品上,样品会吸收部分光线,并将其所具有的能量转换为热能或散发出去。

光谱仪接收样品经过后发出的光线,通过检测器将光谱信号进行记录和处理,最终形成光谱数据以进行分析。

三、光谱数据处理方法在利用光谱数据进行土壤有机质评价时,首先需要进行数据处理。

光谱数据处理方法主要有以下几种:1.预处理对于光谱信号进行预处理以去除干扰因素。

预处理方法包括基线平滑、二次导数变换、离散傅里叶变换等。

2.特征提取将处理后的光谱数据提取真正反映土壤有机质含量指示的特征谱线,通过这些谱线可以对土壤有机质的含量进行定量描述。

3.模型建立将提取的光谱特征作为输入,以土壤有机质含量作为输出建立起对应的数学模型,通过模型来计算土壤有机质含量。

四、结论和建议通过光谱数据进行土壤有机质评价可以快速、准确地获取土壤有机质含量信息,从而为农业生产提供有力的支撑。

土壤有机质不同测定方法的对比6篇

土壤有机质不同测定方法的对比6篇

土壤有机质不同测定方法的对比6篇第1篇示例:土壤中的有机质含量是评价土壤肥力和质量的重要指标之一。

有机质对土壤的保肥补肥、促进土壤肥力的提高起着重要的作用。

准确测定土壤有机质含量对于科学施肥、提高耕地利用率,保护生态环境等方面具有重要意义。

目前,测定土壤有机质含量的方法繁多,各具特点,我们可以根据实际需要选择适合的方法进行测定。

常见的测定土壤有机质含量的方法主要有以下几种:全硫酸铵法、加热法、碱解-蒒灰法、湿法酸解法、直接抽提法等。

不同的方法在原理、操作步骤、准确度和适用范围等方面有所区别,下面我们来对比分析一下。

1. 全硫酸铵法:全硫酸铵法是一种常用的土壤有机质测定方法,其原理是土壤样品在加热时有机质会被分解成氨气和水,测定生成氨气的量来反映土壤中的有机质含量。

此方法操作简单,结果准确,但需要使用昂贵的仪器设备,且操作过程中需注意安全。

2. 加热法:加热法是一种常见的土壤有机质测定方法,其原理是通过将土壤样品在高温下加热,使有机质分解成气体脱除,通过称量前后的重量差来计算有机质含量。

该方法操作简单,成本低,但相对来说准确度略有不足。

3. 碱解-蒒灰法:碱解-蒒灰法是通过将土壤样品与氢氧化钠碱解后,经加热-石蕊灰处理得出的有机质量。

该方法操作相对复杂,但准确度高,适用范围广。

4. 湿法酸解法:湿法酸解法是以稀硝酸和硫酸为酸液,在高温下酸解土壤有机质,通过测定生成的氨气量来计算土壤有机质含量。

该方法操作简单,准确度高,但酸性强,需注意安全。

5. 直接抽提法:直接抽提法是通过将土壤样品与有机溶剂反应,使有机质在溶剂中溶解,通过溶液的浓度来计算有机质含量。

该方法适用范围广,操作简单,准确度高。

不同的测定土壤有机质含量的方法各有特点,选择合适的方法取决于具体的实际情况。

在进行测定时,应根据试验目的、仪器设备的条件和测定精度的要求等因素进行选择,以保证结果的准确性和可靠性。

希望通过以上对比分析,能够为大家在测定土壤有机质含量时提供一些参考。

利用高光谱遥感预测土壤有机碳

利用高光谱遥感预测土壤有机碳

利用高光谱遥感与可见光-近红外段光谱预测土壤有机碳(Soil organic carbon prediction by hyperspectral remote sensing and field vis-NIR spectroscopy)摘要本文利用可见光和近红外反射(vis-NIR)高光谱图像的最接近遥感数据。

预测土壤有机碳(SOC)的结果进行比较,土样采集在Narrabri地区,在澳大利亚新南威尔士(NSW)的北方占据优势的是Vertisols(变性土),这个地区Vis-NIR的光谱采集用AgriSpec便携式光谱仪(350 - 2500nm)和远距离机载高光谱传感器卫星(400 - 2500nm)。

,利用偏最小二乘法回归法(PLSR),偏最小二乘回归法(PLSR:partial least squares regression):是一种新型的多元统计数据分析方法,它主要研究的是多因变量对多自变量的回归建模,特别当各变量内部高度线性相关时,用偏最小二乘回归法更有效。

另外,偏最小二乘回归较好地解决了样本个数少于变量个数等问题。

偏最小二乘法是集主成分分析、典型相关分析和多元线性回归分析3种分析方法的优点于一身。

它与主成分分析法都试图提取出反映数据变异的最大信息,但主成分分析法只考虑一个自变量矩阵,而偏最小二乘法还有一个“响应”矩阵,因此具有预测功能。

)使用的近距离与星载遥感光谱资料预测土壤有机碳含量。

近距离和远距离遥感数据的光谱分辨率是不会影响预测精度的。

然而,在相同光谱分辨率情况下,利用高光谱预测土壤有机碳含量精度没有用Agrispec 便携式光谱仪(resampled )的精度高。

结果显示,用土卫七(Hyperion就是土卫七,围绕着土星运转,1848年9月16日发现Hyperion 是世界上第一个成功发射的星载民用成像光谱仪,也是是目前少数几个仍在轨运行的星载高光谱成像仪。

其搭载卫星EO-1(EarthObserving-1)是美国NASA为接替Landsat7而研制的新型地球观测卫星,于2000年11月21日发射升空。

基于高光谱的土壤有机碳含量预测研究

基于高光谱的土壤有机碳含量预测研究

基于高光谱的土壤有机碳含量预测研究张静【摘要】[Objective] To predict soil organic carbon content.[Method] Surface soil was detected by high spectrometer spectrometric and spectral data was treated,through stepwise multiple linear regression (SMLR) and partial least-squares regression (PLSR) method,soil organic carbon content was predicted,and the accuracy of the two models was compared.[Result] The accuracy of PLSR model was higher than SMLR model.[Conclusion] PLSR method is better than SMLR method in forecasting organic carbon.%[目的]对土壤有机碳含量进行预测研究.[方法]利用高光谱仪对表层土壤进行光谱测定并且进行光谱数据的预处理,通过多元线性逐步回归(SMLR)和偏最小二乘回归(PLSR)方法对土壤有机碳含量进行预测,并对2种模型的精度进行比较.[结果]LSR模型的精度高于SMLR模型.[结论]偏最小二乘回归法优于多元逐步回归法,对有机碳的预测具有更好的效果.【期刊名称】《安徽农业科学》【年(卷),期】2018(046)002【总页数】4页(P1-3,7)【关键词】土壤;有机碳;高光谱;多元线性逐步回归;偏最小二乘回归【作者】张静【作者单位】安徽师范大学,安徽芜湖241000【正文语种】中文【中图分类】S126土壤是人类生活和农业生产的支撑和基础,对土壤类型的区分可以帮助人类因地制宜地利用土壤,同时也给土地利用规划、土地评价、农业技术推广提供重要的依据。

土壤有机质高光谱估算研究进展

称, 包 括 土 壤 中的 动 物 、 植 物 及 微生 物 残 体 的处 于不 同分 解 、 合 成 阶 段 的各种 产物 … . 土 壤有 机质 不但 是 土 壤 的重 要 组成 部 分 , 而 且 是土 壤发 育 过 程 的 重要 标 志 , 对 土壤 的物 理 、 化 学 性 质有 重 要 影 响. 对 土
加 强 土壤信 息 获取新 技术 的研 究 , 实现 土 壤 信息 快 速 获取 . 土壤 学
家 们 对实现 利 用遥感 方法 进行 土壤研 究 向往 已久 . 自从 1 9 7 2年 L a n d -
s a t . 1 发 射 升空 以来 , 人 们 就 非 常希 望 能 够 利 用 土壤 遥 感 进 行 土 壤 环 境 评价 与制 图 _ 4 J . 为实 现这 一 目标 , 土 壤学 家 主要 进行 了 2个方 面 的 研究 : 一 是 室 内光谱 理 论 与 应 用 研究 , 估 算 土壤 成 分 及 含 量 , 并 试 图 分 析及 解释 其机 理 ; 二 是 实现 真 正 的 高光 谱 遥感 或 高 光谱 制 图 , 利 用
郑光辉 , 等. 土壤有机质高光谱估算研究进展
4 8 2
Z H EN G Gu a n g h u i . e t a 1 . Ch a r a c t e r i s t i c s 0 f e v a p o r a t i o n c h a n g e a n d i t s i n l f u e n c e f a c t o r s d u r i n g 1 9 6 9 —2 0 l 1 i n C h a o z h o u
1 南京信息工程大学 遥感学院, 南京 , 2 1 0 0 4 4 2 南京 市浦 口区气象局 , 南京 , 2 1 1 8 0 0

遥感技术在农业方面的应用进展汇总

遥感技术在农业方面的应用进展汇总推动智慧农业、设施农业的发展,全国各地的科研院所都投入不少人力物力,跨行业跨学科展开了研究应用,利用物联网、遥感遥测、人工智能、机器视觉、深度学习、影像采集等技术,为实现农业生产管理的数字化、智能化、自动化而努力。

无数关于农业的科研论文也陆续发表,其中,遥感技术在农业方面的应用已然有了一定进展。

一、机载遥感系统应用精准农业的有人机载成像系统,由安装在农用飞机上的消费级相机组成的系统,详细描述了多光谱相机、高光谱相机和热成像相机等部分定制,和商用机载成像系统。

并举例应用实例,说明如何将不同类型的遥感图像用于精准农业应用中的作物生长评估和作物病虫害管理。

二、大尺度区域水田空间格局及生态服务基于1990—2015年土地利用遥感监测数据,利用GIS的空间分析功能,探究长江经济带水田空间格局动态变化特征。

结果表明水田规模持续缩减,与经济建设及水产养殖的发展、其他生态系统转化、及生态系统服务,有助于揭示长江流域水田的时空变化过程,及其对各项生态系统服务的影响,可为区域土地利用规划、农业政策与生态可持续发展提供理论支持。

三、水稻含水量无人机遥感监测利用多旋翼无人机低空遥感平台,获取不同生育期水稻冠层的RGB图像和多光谱图像,通过提取植被指数和纹理特征,分析水稻的动态生长变化,并构建了基于随机森林回归方法的含水量预测模型。

试验结果表明,基于无人机遥感技术监测水稻含水量是可行的,可为农田精准灌溉、田间管理决策提供新思路。

四、植被分类中的对比分析利用一景AVIRIS高光谱植被影像,从分类精度的角度,提取方法在高光谱影像植被分类中的性能。

实验结果为后续改进空-谱特征方法及其两者有效结合,进一步提高植被分类正确率提供了参考。

五、寒地水稻叶片叶绿素含量遥感反演研究通过分析寒地水稻关键生育期叶片高光谱反射率信息,同时结合PROSPECT模型叶绿素含量吸收系数,参考借鉴现有高光谱植被指数的构造方法和形式,利用相关性分析、连续投影法、遗传算法优化的粗糙集属性简约法,进行高光谱特征选择,结果表明;ORVI能够作为快速反演水稻叶绿素含量的高光谱植被指数,为寒地水稻叶绿素含量高光谱遥感诊断,及管理决策提供了的客观数据支撑和模型参考。

高光谱遥感在土壤质地识别中的应用研究

高光谱遥感在土壤质地识别中的应用研究概述:土壤是地球上非常重要的资源之一,对于农业生产、生态环境和水资源管理具有重要的意义。

而土壤的质地是土壤的基本特性之一,对于农业生产和土地利用规划起着至关重要的作用。

传统的土壤质地识别方法需要进行大量的野外调查和实验室测试,耗时耗力且成本高昂。

而高光谱遥感技术能够获取土壤光谱信息,可通过对光谱数据的分析,实现对土壤质地的自动、快速、准确识别。

一、高光谱遥感技术简介高光谱遥感是指利用遥感卫星或航空遥感平台获取地物连续光谱的技术。

与传统遥感技术相比,高光谱遥感能够提供地物更丰富、更准确的光谱信息,有效地提高了对地物的识别能力。

高光谱遥感数据通常包含了成百上千个窄带的光谱波段,覆盖了可见光、近红外和短波红外等范围。

二、高光谱遥感在土壤质地识别中的应用1. 土壤光谱特征分析高光谱遥感技术通过获取土壤的光谱特征,能够发现土壤质地与光谱之间的关系。

对土壤光谱数据进行分析可以获得不同波段对土壤质地的敏感性信息。

例如,通过对可见光和近红外波段的反射率进行分析,可以发现不同质地土壤在光谱曲线形状上的差异,从而实现对土壤质地的初步区分。

2. 光谱指数方法光谱指数是一种通过计算光谱波段间的比值或差值来表征地物特征的方法。

在土壤质地识别中,常用的光谱指数包括归一化差异植被指数(NDVI)、土壤调查等级指数(SCI)等。

这些指数通过光谱数据进行计算,能够从土壤光谱中提取出与土壤质地相关的特征,达到快速准确识别土壤质地的目的。

3. 光谱分类方法光谱分类方法是利用计算机算法对高光谱遥感数据进行处理和分析,将地物进行分类。

其中最常用的方法是基于支持向量机(SVM)算法的分类。

该方法通过训练样本对土壤质地进行分类器的训练,然后使用训练好的分类器对未知样本进行分类,即可以对土壤质地进行自动识别。

光谱分类方法能够在较高的精度下对土壤质地进行分析,实现了对大规模土地的快速识别。

三、高光谱遥感在土壤质地识别中的优势与挑战1. 优势(1)高光谱遥感技术可以快速获取大面积土地的光谱信息,提供了便捷的数据来源。

基于连续小波变换的土壤有机质含量高光谱反演

中国农业科技导报ꎬ2021ꎬ23(5):132-142JournalofAgriculturalScienceandTechnology㊀收稿日期:2020 ̄08 ̄24ꎻ接受日期:2020 ̄10 ̄24㊀基金项目:内蒙古自治区科技重大专项(2019ZD003)ꎮ㊀联系方式:陈昊宇E ̄mail:chenhaoyu0807@163.comꎻ∗通信作者杨光E ̄mail:yg331@126.com基于连续小波变换的土壤有机质含量高光谱反演陈昊宇ꎬ㊀杨光∗ꎬ㊀韩雪莹ꎬ㊀刘昕ꎬ㊀刘峰ꎬ㊀王宁(内蒙古农业大学沙漠治理学院ꎬ内蒙古自治区风沙物理与防沙治沙工程重点实验室ꎬ呼和浩特010010)摘㊀要:以托克托县境内120个土壤有机质含量以及对应光谱数据为数据源ꎬ探究了不同土壤类型与土地利用类型下土壤有机质高光谱反演研究的可行性ꎬ采用连续小波变换对原始光谱(R)㊁光谱倒数(1/R)㊁光谱对数(LnR)㊁光谱一阶微分(Rᶄ)进行分解生成小波系数并与土壤有机质进行相关系分析ꎬ提取特征波段建立BP神经网络与支持向量机模型(SVM)ꎮ结果表明:①R㊁1/R㊁LnR㊁Rᶄ与土壤有机质相关系数经过连续小波变换后ꎬ较之前增加了0.204㊁0.090㊁0.199㊁0.252ꎬ表明连续小波变换可深度挖掘光谱潜在信息ꎬ提升与有机质含量之间的相关系数ꎮ②未经过连续小波处理前ꎬSVM无法实现对当地土壤有机质含量的预测ꎬ经过处理后ꎬ模型SVM ̄CWT ̄R与SVM ̄CWT ̄Rᶄ的精度决定系数分别达到了050㊁0.56ꎬ均方根误差为0.17㊁0.15ꎬ相对分析误差为1.62㊁1.53ꎬ实现了对土壤有机质的有效估算ꎮ③经过连续小波变换后BP神经网络预测模型结果得到提升ꎬ其中BP ̄CWT ̄LnR预测模型效果最佳ꎬ精度决定系数达到0.76ꎬ较之前BP ̄LnR提升了0.2ꎻ均方根误差达到0 15ꎬ降低0.04ꎻ相对分析误差为2.12ꎬ增加了0.87ꎮ因此利用BP ̄CWT ̄LnR高光谱反演模型进行区域土壤有机质遥感监测ꎬ可为当今精准农业提供理论参考与技术支持ꎮ关键词:连续小波变换ꎻBP神经网络ꎻ支持向量机ꎻ精准农业doi:10.13304/j.nykjdb.2020.0742中图分类号:S127㊀㊀㊀文献标识码:A㊀㊀㊀文章编号:1008 ̄0864(2021)05 ̄0132 ̄11HyperspectralInversionofSoilOrganicMatterContentBasedonContinuousWaveletTransformCHENHaoyuꎬYANGGuang∗ꎬHANXueyingꎬLIUXinꎬLIUFengꎬWANGNing(KeyLaboratoryofAeolianPhysicsandDesertificationControlEngineeringfromInnerMongoliaAutonomousRegionꎬCollegeofDesertControlScienceandEngineeringꎬInnerMongoliaAgriculturalUniversityꎬHohhot010010ꎬChina)Abstract:Takingorganicmattercontentsof120soilsamplesandcorrespondingspectraldatainTuoketuoCountyasdatasourcesꎬthefeasibilitiesofhyperspectralinversionofsoilorganicmattersunderdifferenttypeofsoilsandlandsofdifferentusewereexplored.Theoriginalspectrum(R)ꎬspectralreciprocal(1/R)ꎬspectrallogarithm(LnR)andspectralfirst ̄orderdifferential(Rᶄ)weredecomposedbycontinuouswavelettransformtogeneratewaveletcoefficientsꎬandthecorrelationbetweensoilorganicmatterandwaveletcoefficientswasanalyzedꎬandBPneuralnetworkandsupportvectormachine(SVM)wereestablishedbyextractingthecharacteristicbands.Theresultswerefollowed.①ThecorrelationcoefficientsbetweenRꎬ1/RꎬLnRꎬRᶄandsoilorganicmatterwereincreasedby0.204ꎬ0.090ꎬ0.199and0.252aftercontinuouswavelettransformꎬrespectivelyꎬwhichshowedthatcontinuouswavelettransformcoulddeeplyminethepotentiallyspectralinformationandenhancethecorrelationwithorganicmattercontent.②BeforecontinuouswaveletprocessingꎬSVMcouldnotpredictthecontentofsoilorganicmatterꎬwhileafterprocessingꎬtheaccuracies(R2)ofSVM ̄CWT ̄RandSVM ̄CWT ̄Rᶄwere0.50and0.56ꎬRootmeansquareerrors(RMSE)were0.17and0.15ꎬresidualpredictivedeviations(RPD)were1.62and1.53ꎬrespectivelyꎬwhichrealizedtheeffectiveestimationofSOM.③AftercontinuouswavelettransformꎬtheresultsofBPneuralnetworkpredictionmodelwereimproved.AmongthemꎬBP ̄CWT ̄LnRpredictionmodelhadthebesteffectꎬR2was0.76ꎬwhichwashigherthanBP ̄lnRꎬRMSEwas0.15reducedby0.04ꎬRPDwas2.12increasedby0.87.ThereforeꎬtheBP ̄CWT ̄LnRhyperspectralinversionmodelcouldprovidetheoreticalreferenceandtechnicalsupportforprecisionagriculture.Keywords:continuouswavelettransformꎻBPneuralnetworkꎻsupportvectormachineꎻprecisionagriculture㊀㊀精准农业作为目前农业发展的主要方向ꎬ是一种基于信息和知识管理的现代化生产系统ꎬ主要是通过3S(GPS㊁GIS和RS)技术与现代农业相结合ꎬ最大限度地提高农业生产力ꎮ所以快速㊁无损㊁精确地获取土壤中水分㊁养分的空间分布成为了实现精准农业的关键环节ꎬ近年来ꎬ光谱分析在土壤化学分析领域得到了迅猛发展ꎬ为实现土壤养分的快速诊断提供了新思路[1]ꎮ有机质是土壤养分供应能力和肥力的重要指标之一ꎬ在全球碳循环中发挥着重要作用ꎮ因此ꎬ快速准确地估测土壤有机质含量对于发展精准农业具有重要意义[2]ꎮ传统的土壤有机质测定方法虽然精度比较高ꎬ但周期较长㊁成本较高ꎬ只能达到瞬测量ꎬ很难进行长时间大面积测量ꎮ高光谱遥感具有波段多㊁波段窄㊁信息丰富和实时高效等特点ꎬ为快速测量土壤有机质含量提供了一种新的方法和手段[3]ꎮ国内外已经有大量研究表明ꎬ通过对光谱数据进行不同的数学变换(主要通过对光谱进行倒数㊁对数㊁微分㊁平方根㊁吸收峰深度㊁包络线去除等方法)可以有效提高光谱数据与土壤有机质含量之间的相关系数ꎬ有效筛选出光谱信息中的敏感波段[4]ꎮ现在各学者主要将研究重心放到了模型建立上[5]ꎬ普遍运用的线性模型有多元逐步回归与偏最小二乘回归[6]ꎻ常见的非线性模型包括BP神经网络[7]㊁支持向量机[8]㊁决策树[9]等ꎬ而且随着非线性模型算法的逐步改良与完善ꎬ在土壤有机质含量估算中已经成为不可取代的一部分ꎮ随着小波算法的改进与发展ꎬ最初仅运用于植物叶绿素㊁冠层成分含量预测中[10 ̄11]ꎬ目前已成为土壤养分预测的热点问题[12 ̄13]ꎬ连续小波变换是目前被广泛应用的一种方法ꎮ王祥浩[14]选择土地裸露地区为样区ꎬ利用神经网络算法对光谱连续小波变换㊁一阶导数㊁对光谱的平均值处理㊁光谱背景及深度4种方法建模ꎬ模型结果表明ꎬ小波变换方法得到的神经网络模型精度最高ꎻ包青岭等[15]选择渭干河-库车河三角洲具有代表性的干旱区绿洲为研究区ꎬ对光谱进行8层分解ꎬ结果表明小波变换不同分解层ꎬ从低频到高频范围内与土壤有机质含量的相关性呈现先减后增的趋势ꎬ结合随机森岭模型可以对干旱区土壤有机质含量进行有效的估算ꎻ王延仓等[16]以北京东部区潮土为例ꎬ对不同梯度重采样的光谱进行连续小波变换后ꎬ利用偏最小二乘法建立模型ꎬ结果表明连续小波分析算法可深入挖掘土壤光谱内的有益信息ꎬ提升对有机质含量的估测能力ꎬ与土壤高光谱反射率相比ꎬ经连续小波技术处理后ꎬ模型精度得到了有效的提升ꎻ叶红云等[17]同样针对干旱区土壤ꎬ通过对两种常用光谱变换Rᶄ㊁Ln(1/R)进行连续小波变换建立偏最小二乘模型ꎬ结果表明连续小波变换不会因人类干扰程度的提高而使模型精度大幅度降低ꎬ更加适用于干旱区有机质含量的预测ꎻ林鹏达等[18]通过解决黑土有机质高光谱野外反演的困难ꎬ同样证明了连续小波变换可有效提升模型精度ꎮ小波技术在土壤有机质高光谱反演研究中逐渐趋于成熟ꎬ但目前学者的研究多数都在同一土壤类型下或同一区域内ꎬ对于不同土壤类型及土地利用下土壤有机质高光谱反演是否存在影响的研究目前并不多ꎮ本文研究区内土壤类型主要包括3类:沙壤土㊁栗钙土㊁盐碱土ꎬ且部分区域土壤盐渍化程度严重ꎬ导致土壤养分空间分布上存在较大差异ꎬ取样表层土地利用类型主要包括:耕地㊁林地㊁草地㊁盐渍地㊁荒地ꎮ通过对原始光谱(R)㊁原始光谱倒数(1/R)㊁原始光谱对数(LnR)以及原始光谱一阶微分(Rᶄ)4种不同情况进行连续小波变换ꎬ利用BP神经网络以及支持向量机2种模型ꎬ探究了不同土壤类型与不同土地利用类型下是否会对土壤有机质高光谱反演模型产生影响ꎬ小波变换前后土壤有机质反演模型的精度ꎬ旨为区域土壤有机质含量监测及实现精准农业提供理论与技术支持ꎮ1㊀材料与方法1.1㊀研究区概况托克托县隶属于内蒙古自治区呼和浩特市ꎬ位于自治区中部㊁大青山南麓㊁黄河上中游分界处北岸的土默川平原上(图1)ꎮ地理坐标东经3315期陈昊宇等:基于连续小波变换的土壤有机质含量高光谱反演111ʎ2ᶄ30ᵡ 111ʎ32ᶄ21ᵡ㊁北纬40ʎ5ᶄ55ᵡ 40ʎ35ᶄ15ᵡꎬ总面积1409.67km2ꎬ平均海拔1117mꎬ属于温带大陆性干旱气候ꎬ年均气温7.3ħꎬ年均降雨362mmꎮ托克托县耕地总面积达400km2ꎬ其中古城镇㊁新营子镇和五申镇的耕地较多ꎬ占全县耕地面积的60%以上[19]ꎬ主要作物包括小麦㊁玉米㊁莜麦ꎮ工农业及生产生活用水主要来源于大黑河和黄河水资源ꎬ整个地形以大黑河为轴ꎬ呈现由丘陵向平原过渡的趋势ꎬ地势为东南高㊁西北和西南低ꎮ东南向西北土壤类型依次为栗钙土㊁砂壤石灰性冲积土㊁盐渍化石灰性冲积土[20]ꎬ土壤类型的不同导致土壤养分存在差异性分布ꎮ植被类型从西向东依次为草甸草原㊁干草原和退化灌丛草原分布ꎮ以Landsat8OLI影像为基础数据源ꎬ运用人工目视解译与BP神经网络分类法得到托克托县2019年7月份土地利用数据ꎬ其中耕地面积最大为730.12km2ꎬ占51.79%ꎻ林草地338.7km2ꎬ占24.02%ꎻ盐碱地141.1km2ꎬ占10.00%ꎮ详细土地利用空间分布见图1ꎮ图1㊀土样采集点及土地利用空间分布Fig.1㊀Collectionpointsofsoilsamplesandspatialdistributionoflanduse1.2㊀研究方法1.2.1㊀土样采集与处理㊀土壤样本点均匀地分布在托克托县境内ꎬ采集方法为五点采样法ꎬ采集深度为0 20cmꎬ共采集120个点ꎮ采集的土样置于通风干燥室内进行自然风干㊁研磨ꎬ过10目筛ꎬ进行土壤光谱测定ꎻ过100目筛ꎬ采用重铬酸钾外加热法进行土壤有机质含量测定ꎮ1.2.2㊀光谱测量及光谱处理㊀土壤光谱于暗室内测量ꎬ采用SVCHR ̄1024(北京东方佳气科技有限公司)便携式光谱仪ꎬ光谱范围在350~2500nmꎮ在350~1000nm波段之间光谱分辨率ɤ3.5nmꎻ在1000~1850nm波段之间ꎬ光谱分辨率ɤ9.5nmꎻ在1850~2500nm波段之间ꎬ光谱分辨率ɤ6.5nmꎮ光源采用与太阳光接近的50W卤素灯ꎬ将土壤样品放入深2cm㊁宽10cm的黑色器皿内ꎬ用直尺将土壤表面刮平ꎬ探头距离土样10cmꎬ光源距离土壤表面30cmꎬ天顶角为15ʎꎮ测量前用白板进行标定ꎬ每个土样采集5条光谱作为该土样的光谱数据ꎮ由于受噪音与仪器暗电流的的影响ꎬ导致光谱数据混入噪音等信息ꎬ因此删除350~399nm和2400~2500nm的波段ꎬ采用五点平滑法对光谱进行平滑处理ꎬ并将光谱重采样至5nmꎬ同时对原始光谱(R)进行一阶微分(Rᶄ)㊁倒数(1/R)㊁对数(LnR)等传统数学变换ꎮ1.2.3㊀连续小波变换㊀采用连续小波变换ꎬ并用Mexh小波母函数对原始光谱㊁原始光谱的倒数㊁对数㊁一阶微分进行10层小波变换ꎬ生成一系列小波系数ꎮΨaꎬb=1㊀aΨλ-baæèçöø÷(1)式中ꎬa为伸缩因子ꎬb为平移因子ꎬλ为土壤高光谱数据的波段数ꎮWfaꎬb()=fꎬΨaꎬb()=ʏ+ɕ-ɕfλ()Ψaꎬbλ()dy(2)式中ꎬfλ()为土壤光谱反射率ꎬ小波系数Wfaꎬb()包含二维ꎬ分别为波长(350~2500)与分解尺度(1ꎬ2ꎬ3 10)ꎬ故小波系数行为尺度数ꎬ列为波长数的矩阵[16]ꎮ1.2.4㊀模型及精度验证㊀采用BP神经网络与支持向量机模型(supportvectormachineꎬSVM)建立土壤有机质预测模型ꎬ支持向量机采用线性核函数ꎬ相对于径向基函数(radialbasisfunctionꎬ431中国农业科技导报23卷RBF)来说计算高效ꎬ不易过拟合ꎮBP神经网络的迭代次数设置为1000ꎬ学习率0.01ꎬ训练的均方根误差(rootmeansquareerrorꎬRMSE)小于0.001ꎮ依据相关系数筛选的特征波段以及小波系数作为自变量ꎬ土壤有机质含量为因变量ꎬ分别建立模型ꎬ模型精度采用决定系数(R2)㊁均方根误差(RMSE)㊁相对分析误差(relativepercentdeviationꎬRPD)以及1ʒ1线共同评价ꎮR2表征模型的稳定性ꎬ越接近于1模型越稳定ꎬ拟合程度越好ꎮ均方根误差(RMSE)用来检验模型的预报能力ꎬRMSE越小则表明模型的估测能力越好ꎮRPD是样本的标准差与RMSE的比值ꎬRPD<1.4时ꎬ模型无法对样品进行预测ꎻ1.4ɤRPD<2时ꎬ模型效果一般ꎬ可以用来对样品进行粗略评估ꎻRPDȡ2时ꎬ模型具有极好的预测能力ꎮ1ʒ1线表示实测值与预测值构成的点偏离y=x线的程度[21]ꎮ2㊀结果与分析2.1㊀土壤有机质含量统计分析建模样品集㊁不同土地利用方式㊁不同土壤类型下土壤有机质含量描述性统计见表1ꎮ本研采样点内土地利用方式主要包括林地㊁草地㊁耕地㊁盐渍地ꎬ土壤有机质在草地内均值含量最大(0 80%)ꎬ其次为林地(0.72%)㊁耕地(0.67%)㊁盐渍地有机质含量最低(0.63%)ꎻ土壤有机质含量最大值位于耕地(1.28%)ꎬ最小值位于林地(0 19%)ꎮ采样点内主要土壤类型为栗钙土㊁沙壤土㊁盐碱土ꎬ沙壤土有机质含量最高(0.77%)ꎬ其次为盐碱土(0.68%)和栗钙土(0.67%)ꎬ土壤有机质含量最大值位于沙壤土内(1.28%)ꎬ最小值位于盐碱土内(0.19%)ꎮ表1㊀土壤有机质含量描述性统计结果Table1㊀Descriptivestatisticsresultsoforganicmattercontentinsoilsamples项目Item样品集及类型Samplesetandtype土样数Numberofsamples最大值Maximum/%最小值Minimumvalue/%均值Meanvalue/%标准差Standarddeviation/%模型样品集Modelsampleset样品全集Wholeset1201.200.1940.710.276建模集Modelingset901.200.190.720.227验证集Validationset301.280.250.70.257土地利用方式Landusepattern耕地Cultivatedland551.280.200.670.20林地Woodland201.150.190.720.23草地Grassland251.180.250.800.21盐渍地Salinesoil201.030.250.630.22土壤类型Soiltype栗钙土Chestnutsoil541.20.20.670.22沙壤土Sandyloam511.280.250.770.24盐碱土Saline ̄alkalisoil151.150.190.680.222.2㊀土壤反射光谱特征对R㊁1/R㊁LnR㊁Rᶄ进行小波变换ꎬ变换结果如图2所示ꎬR㊁1/R㊁LnR光谱曲线较为平滑ꎬ分解曲线随波峰波谷变化.Rᶄ其光谱曲线并不规则存在较多波峰波谷ꎬ分解小波系数与前三者不同ꎮR㊁1/R㊁LnR㊁Rᶄ分解后ꎬ小波系数均随分解尺度的增加而增加ꎬ同时可以看出ꎬ由Mexh小波母函数进行的连续小波变换ꎬ对于光谱波峰与波谷有较高的敏感性ꎬ对于放大㊁挖掘光谱信息有着显著的作用ꎮ2.3㊀相关性分析2.3.1㊀不同导数变换光谱与土壤有机质含量相关性㊀土壤有机质含量与光谱相关性曲线及敏感波段见图3ꎮR与土壤有机质含量呈负相关关系(相关系数r=-0.463)ꎬ主要集中于735~780nm处波段ꎻ1/R与土壤有机质的相关性则与R相反ꎬ呈正相关关系(r=0.462)ꎬ集中于600~800nm与1800~2200nm处波段ꎻLnR的相关性曲线图与R相关性曲线类似ꎬ总体呈现负相关关系ꎬ相关系数(r=-0.465)ꎬ主要集中于745~7955315期陈昊宇等:基于连续小波变换的土壤有机质含量高光谱反演nm处的波段ꎻRᶄ相关性在500nm(r=-0.589)与1400nm(r=-0.411)处为负相关ꎬ在800nm(r=0.408)与1380nm(r=0.412)处为正相关ꎬ相关系数曲线变换趋势与前三者不同ꎬ呈无规律变化ꎮ2.3.2㊀不同分解尺度小波系数与土壤有机质含量的相关性㊀图4为不同光谱变换方式经过连续小波变换后与土壤有机质含量的相关系数矩阵图ꎬ其中红色代表相关性高的区域ꎬ蓝色代表相关性低的区域ꎮR在800~1000㊁1400~1600nm处相关性明显增加ꎬ在500㊁800㊁2200nm波段处相关系数达到最大值(r=0.667)ꎻ1/R在800~1200nm处相关系数达到最大值(r=0.552)ꎬ在2400~2500nm处相系数达到0.4ꎬ受噪音和仪器本身的影响ꎬ此波段的相关系数不进行相关性参考ꎻLnR在分解尺度1下相关性较低ꎬ在2~10尺度下ꎬ相关性出现最大值(r=0.664)ꎻRᶄ相关性主要集中在500~900㊁1200~1600㊁2100~2300nm处ꎮ筛选的敏感波段与尺度如表2所示ꎮ有效的光谱信息主要存在于低分解尺度ꎬ随分解尺度的增加呈递减趋势ꎬ相关性最大值较未处理前分别增加了0.204㊁0.09㊁0.199㊁0.252ꎬ对于挖掘潜在光谱信息有着重要意义ꎮ2.4㊀土壤有机质高光谱模型建立2.4.1㊀BP神经网络预测模型㊀采用BP神经网络构建反演模型ꎬ结果如表3所示ꎮ未进行连续小波变换处理的模型中ꎬBP ̄R与BP ̄Rᶄ效果较好ꎬR2分别为0.69和0.73ꎬRPD为1.45与1.53ꎬ模型能粗略估算土壤有机含量ꎬBP ̄LnR与BP ̄1/R样本外预测能力较差ꎬ同时RPD未达到1.4以上ꎬ不能对土壤有机质未能进行有效预测ꎻ连续小图2㊀连续小波变换光谱特性Fig.2㊀Spectralcharacteristicsofcontinuouswavelettransform631中国农业科技导报23卷图3㊀土壤光谱相关性曲线及敏感波段Fig.3㊀Correlationcurveandsensitivebandofsoilspectrum图4㊀土壤有机质与小波系数相关性Fig.4㊀Correlationbetweensoilorganicmatterandwaveletcoefficients7315期陈昊宇等:基于连续小波变换的土壤有机质含量高光谱反演表2㊀筛选的敏感波段Table2㊀Sensitivebandforscreening处理方式Treatmentmethod相关系数Correlationcoefficient分解尺度Decompositionscale敏感波段Sensitiveband/nmCWT ̄R0.6671~10214㊁212㊁91㊁91㊁109㊁109㊁109㊁110㊁111㊁112CWT ̄1/R0.5521~7400㊁212㊁108㊁108㊁108㊁108㊁109CWT ̄LnR0.6642~8213㊁212㊁45㊁45㊁45㊁109㊁110CWT ̄Rᶄ0.6641~10212㊁215㊁215㊁216㊁99㊁99㊁99㊁98㊁98㊁98波变换处理之后的模型ꎬ仅BP ̄CWT ̄1/R模型RPD未达到预测水平ꎬ其余3种模型R2与RPD较未处理前均有所增加ꎬRMSE均减少ꎬ其中BP ̄CWT ̄LnR模型预测效果较好ꎬRPD达到2.12可以有效地对土壤有机质进行预测ꎮ将BP ̄CWT处理的4个模型的实测值与预测值进行1ʒ1线分析ꎮ由图5可知ꎬ除BP ̄CWT ̄1/R模型外ꎬ其余模型的实测值与预测值样点基本分布在1ʒ1线附近ꎬBP ̄CWT ̄LnR效果较为明显ꎬ且估算精度高ꎬ可较好地进行土壤有机质含量的估算ꎮ2.4.2㊀支持向量机预测模型㊀SVM构建反演模型ꎬ结果如表4所示ꎮ未经过连续小波处理的光谱特征波段未能较好地对土壤有机质进行预测反演ꎬ经过CWT后模型SVM ̄CWT ̄R与SVM ̄CWT ̄Rᶄ预测结果较之前有较大的提升ꎬR2分别达到了0.50与0.56ꎬ二者RPD均达到1.4以上ꎬ可以粗表3㊀土壤有机质BP神经网络估测模型结果Table3㊀ResultsofBPneuralnetworkestimationmodelforsoilorganicmatter模型Model建模集ModelingsetR2RMSE验证集ValidationsetR2RMSERPDBP ̄R0.690.170.520.191.45BP ̄1/R0.680.170.330.221.25BP ̄LnR0.560.190.240.221.25BP ̄Rᶄ0.730.160.520.181.53BP ̄CWT ̄R0.800.140.540.171.62BP ̄CWT ̄1/R0.640.180.210.280.98BP ̄CWT ̄LnR0.760.150.740.132.12BP ̄CWT ̄Rᶄ0.770.140.660.161.72表4㊀土壤有机质支持向量机估测模型结果Table4㊀Supportvectormachineestimationmodelresultsofsoilorganicmatter模型Model建模集ModelingsetR2RMSE验证集ValidationsetR2RMSERPDSVM ̄R0.210.200.190.201.38SVM ̄1/R0.210.200.200.211.31SVM ̄LnR0.210.200.200.211.31SVM ̄Rᶄ0.430.170.270.201.38SVM ̄CWT ̄R0.500.160.480.171.62SVM ̄CWT ̄1/R0.290.190.160.211.31SVM ̄CWT ̄LnR0.490.160.270.201.38SVM ̄CWT ̄Rᶄ0.560.150.410.181.53831中国农业科技导报23卷图5㊀BP ̄CWT模型土壤实测值与预测值对比Fig.5㊀ComparisonofmeasuredvalueandpredictedvalueofBP ̄CWTmodel略地对土壤有机质进行预测ꎮ同时根据图6ꎬSVM ̄CWT模型进行1:1线分析ꎬ二者实测值与预测值分布情况在4种模型下较好ꎬ虽然模型SVM ̄CWTLnR分布同样较为集中ꎬ但其样本外预测情况较差(RPD=1.38)ꎬ综合考虑不对其进行土壤有机质预测ꎮ结合表3和表4的结果分析ꎬ连续小波变换能够有效地提升模型精度与模型泛化能力ꎬ对于光谱信息挖掘有着重要意义ꎬBP神经网络与支持向量机对CWT ̄R与CWT ̄Rᶄ都能够提升R2减少RMSEꎬ可对土壤有机质做出较好的预测ꎮ虽然BP神经网络与支持向量机在处理非线性回归问题中有较强的能力ꎬ但本身模型中存在不稳定性ꎬ对模型的环境设置同样要求较高ꎬ所以未能对所有数据集进行良好的预测ꎮ3㊀讨论本研究采用连续小波变换对光谱进行处理ꎬ用BP神经网络与支持向量机(SVM)两种模型对土壤有机质含量进行反演预测ꎮ未经过连续小波变换前ꎬR㊁1/R㊁LnR㊁Rᶄ与土壤有机质的相关系系数最大值分别为-0.463㊁0.462㊁-0.465㊁0.589ꎬ可以看出ꎬRᶄ与土壤有机质的相关系数最高ꎬ与吴倩等[22]㊁张新乐等[23]的研究结果相同ꎻ经过连续小波变换后ꎬCWT ̄R㊁CWT ̄1/R㊁CWT ̄LnR㊁CWT ̄Rᶄ相关系数最大值分别为0 667㊁0.552㊁0 664㊁0.662ꎬ较之前分别增加了0 20㊁0.09㊁0.19㊁0.07ꎮ王延仓等[1]㊁于雷等[4]㊁叶红云等[17]等同样证明连续小波变换可有效提高与土壤有机质含量的相关系数ꎮ不同分解尺度对于光谱数据的深度挖掘有着重要意义ꎬ本研究只利用Mexh小波母函数进行处理ꎬ未对其他函数进行考虑ꎬ分解层数同样是根据前人经验所得[4ꎬ10]ꎬ小波技术的研究与发展仍然有很大的探索空间ꎮ相对于两种模型来看ꎬ未进行连续小波处理9315期陈昊宇等:基于连续小波变换的土壤有机质含量高光谱反演图6㊀利用SVM ̄CWT模型土壤实测值与预测值的对比Fig.6㊀ComparisonofsoilmeasuredvalueandpredictedvalueofSVM ̄CWTmodel的支持向量机模型中ꎬ只有SVM ̄Rᶄ模型R2最高达到0.43ꎬ其余三者均未到达0.4ꎮ综合多种模型评价方法ꎬ由于其RPD未达到1.4以上ꎬ无法对土壤有机质含量进行预测ꎮ经过连续小波处理后ꎬ各模型的R2有明显提高ꎬ其中SVM ̄CWT ̄R与SVM ̄CWT ̄Rᶄ模型效果较好ꎬR2分别提高了0.29㊁0.13ꎬRPD达到1.62与1.53实现了对土壤有机质有效的预测ꎬ但预测结果较BP神经网络较低ꎮ在BP神经网络预测模型中ꎬ未进行连续小波变换前ꎬBP ̄R与BP ̄Rᶄ预测效果较好ꎬR2达到0.69与0.73ꎬRPD为1.45与1.53ꎻ进行连续小波处理后ꎬ除SVM ̄CWT ̄1/R模型未到达预测效果ꎬ其余3种模型预测结果较之前均有明显改善ꎬ可实现对土壤有机质较好的预测ꎬ其中BP ̄CWT ̄LnR预测模型效果最佳R2达到0.76ꎬRPD达到2.12ꎮ根据1:1线分析图也可看出ꎬ其实测值与预测值分布较为集中ꎬ于雷等[4]㊁叶红云等[17]㊁林鹏达等[18]同样通过连续小波变换有效提升了模型的精度与泛化能力ꎮ针对土壤有机质高光谱反演研究中ꎬ姚聪[24]对耕层土壤通过BP神经网络与支持向量机模型ꎬ反演精度R2分别为0.42与0.67ꎻ叶红云等[17]采用连续小波变换对干旱区土壤有机质反演ꎬ模型精度R2=0.75㊁EMSE=0.71ꎻ谢文[25]在森林土壤有机质反演研究中ꎬBP神经网络模型R2=0 78㊁EMSE=0.77ꎬ支持向量机模型R2=0.87㊁EMSE=0.76ꎮ本研究对耕地㊁林草地㊁盐碱地㊁栗钙土㊁沙壤土㊁盐渍土等不同土地利用类型与土壤类型进行综合反演ꎬ最佳反演模型为BP ̄CWTLnRꎬR2=0.76㊁EMSE=0.15㊁RPD=2.12ꎬ与前人研究的结果基本相符ꎬ证明通过连续小波变换处理ꎬ不同土壤类型与土地利用类型未对土壤反演模型精度产生影响ꎮ所以采用连续小波变换进行光谱数据挖掘ꎬ采用BP ̄CWT ̄LnR神经网络建041中国农业科技导报23卷立反演模型ꎬ可对不同土地利用于土壤类型条件下土壤有机质高光谱反演提供一定的理论支持与应用价值ꎮ参㊀考㊀文㊀献[1]㊀王延仓ꎬ杨秀峰ꎬ赵起超ꎬ等.二进制小波技术定量反演北方潮土土壤有机质含量[J].光谱学与光谱分析ꎬ2019ꎬ39(9):2855-2861.WANGYCꎬYANGXFꎬZHAOQCꎬetal..Quantitativeinversionofsoilbasedonbinarywavelettransform[J].SpectroscopySpectralAnal.ꎬ2019ꎬ39(9):2855-2861. [2]㊀于雷ꎬ洪永胜ꎬ耿雷ꎬ等.基于偏最小二乘回归的土壤有机质含量高光谱估算[J].农业工程学ꎬ2015ꎬ31(14):103-109.YULꎬHONGYSꎬGENGLꎬetal..Hyperspectralestimationofsoilorganicmattercontentbasedonpartialleastsquaresregression[J].Trans.CSAEꎬ2015ꎬ31(14):103-109. [3]㊀钟浩ꎬ李西灿ꎬ翟浩然ꎬ等.耕层土壤有机质高光谱间接估测模型[J].测绘科学技术学报ꎬ2019ꎬ36(1):74-78ꎬ85.ZHONGHꎬLIXCꎬZhAIHRꎬetal..Hyperspectralindirectestimationmodelofsoilorganicmattercontentinploughlayer[J].J.GeomaticsTechnol.ꎬ2019ꎬ36(1):74-78ꎬ85. [4]㊀于雷ꎬ洪永胜ꎬ周勇ꎬ等.连续小波变换高光谱数据的土壤有机质含量反演模型构建[J].光谱学与光谱分析ꎬ2016ꎬ36(5):1428-1433.YULꎬHONGYSꎬZHOUYꎬetal..Inversionofsoilorganicmattercontentusinghyperspectraldatabasedoncontinuouswavelettransformation[J].SpectroscopySpectralAnal.ꎬ2016ꎬ36(5):1428-1433.[5]㊀STENBERGBꎬVISCARRARRAꎬMOUAZENAMꎬeta1..Visibleandnearinfraredspectroscopyinsoilscience[J].Adv.Agron.ꎬ2010ꎬ107:163-215.[6]㊀聂哲ꎬ李秀芬ꎬ吕家欣ꎬ等.东北典型黑土区表层土壤有机质含量高光谱反演研究[J].土壤通报ꎬ2019ꎬ50(6):1285-1293.NIEZꎬLIUXFꎬLYUJXꎬetal..HyperspectralretrievalofsurfacesoilorganicmattercontentinatypicalblacksoilregionofnortheastChina[J].Chin.J.SoilSci.ꎬ2019ꎬ50(6):1285-1293.[7]㊀沈润平ꎬ丁国香ꎬ魏国栓ꎬ等.基于人工神经网络的土壤有机质含量高光谱反演[J].土壤学报ꎬ2009ꎬ46(3):391-397.SHENRPꎬDINGGXꎬWEIGSꎬetal..Retrievalofsoilorganicmattercontentfromhyper ̄spectrumbaseonAnn[J].ActaPedol.Sin.ꎬ2009ꎬ46(3):391-397.[8]㊀沈强ꎬ张世文ꎬ夏沙沙ꎬ等.基于支持向量机的土壤有机质高光谱反演[J].安徽理工大学学报(自然科学版)ꎬ2019ꎬ39(4):39-45.SHENQꎬZHANGSWꎬXIASSꎬetal..Hyperspectralinversionofsoilorganicmatterbasedonsupportvectormachine[J].J.AnhuiUniv.Sci.Technol.(Nat.Sci.)ꎬ2019ꎬ39(4):39-45.[9]㊀BREIMANLꎬFRIEDMANJHꎬOLSHENRAꎬetal..ClassificationandRegressionTrees[M].Belmont:WadsworthInternationalGroupꎬ1984.[10]㊀方圣辉ꎬ乐源ꎬ梁琦.基于连续小波分析的混合植被叶绿素反演[J].武汉大学学报ꎬ2015ꎬ40(3):296-302.FANGSHꎬLEYꎬLIANGQ.Retrievalofchlorophyllcontentusingcontinuouswaveletanalysisacrossarangeofvegetationspecies[J].GeomaticsInform.Sci.WuHanUniv.ꎬ2015ꎬ40(3):296-302.[11]㊀何汝艳ꎬ乔小军ꎬ蒋金豹ꎬ等.小波法反演条锈病胁迫下冬小麦冠层叶片全氮含量[J].农业工程学报ꎬ2015ꎬ31(2):141-146.HERYꎬQIAOXJꎬJIANGJBꎬetal..Retrievingcanopyleaftotalnitrogencontentofwinterwheatbycontinuouswavelettransform[J].Trans.CSAEꎬ2015ꎬ31(2):141-146. [12]㊀陈红艳ꎬ赵庚星ꎬ李希灿ꎬ等.小波分析用于土壤速效钾含量高光谱估测研究[J].中国农业科学ꎬ2012ꎬ45(7):1425-1431.CHENHYꎬZHAOGXꎬLIXCꎬetal..Applicationofwaveletanalysisforestimationofsoilavailablepotassiumcontentwithhyperspectralreflectance[J].Sci.Agric.Sin.ꎬ2012ꎬ45(7):1425-1431.[13]㊀高洪智ꎬ卢启鹏.土壤主要养分近红外光谱分析及其测量系统[J].光谱学与光谱分析ꎬ2011ꎬ31(5):1245-1249.GAOHZꎬLUQP.Nearinfraredspectralanalysisandmeasuringsystemforprimarynutrientofsoil[J].SpectroscopySpectralAnal.ꎬ2011ꎬ31(5):1245-1249.[14]㊀王祥浩.土壤有机质高光谱反演模型研究[J].黑龙江工程学院学报ꎬ2019ꎬ33(5):34-39.WANGXH.Researchonhighspectralinversionmodelofsoilorganicmatter[J].J.HeilongjiangInstituteTechnol.ꎬ2019ꎬ33(5):34-39.[15]㊀包青岭ꎬ丁建丽ꎬ王敬哲ꎬ等.基于随机森林算法的土壤有机质含量高光谱检测[J].干旱区地理ꎬ2019ꎬ42(6):1404-1414.BAOQLꎬDINGJLꎬWANGJZꎬetal..Hyperspectraldetectionoforganicmattercontentbasedonrandomforestalgorithm[J].AridLandGeographyꎬ2019ꎬ42(6):1404-1414.[16]㊀王延仓ꎬ张兰ꎬ王欢ꎬ等.连续小波变换定量反演土壤有机质含量[J].光谱学与光谱分析ꎬ2018ꎬ38(11):3521-3527.WANGYCꎬZHANGLꎬWANGHꎬetal..Quantitativeinversionofsoilorganicmattercontentbasedoncontinuouswavelettransform[J].SpectroscopySpectralAnal.ꎬ2018ꎬ38(11):3521-3527.[17]㊀叶红云ꎬ熊黑钢ꎬ张芳ꎬ等.基于CWT的人类不同程度干扰下干旱区土壤有机质含量估算研究[J].激光与光电子学进展ꎬ2019ꎬ56(5):115-124.YEHYꎬXIONGHGꎬZHANGFꎬetal..CWT ̄Basedestimationofsoilorganicmattercontentinaridareaunderdifferenthumandisturbancedegrees[J].LaserOptoelectronicsProgr.ꎬ2019ꎬ56(5):115-124.[18]㊀林鹏达ꎬ佟志军ꎬ张继权ꎬ等.基于CWT的黑土有机质含量野外高光谱反演模型[J].水土保持研究ꎬ2018ꎬ25(2):46-52ꎬ57.LINPDꎬTONGZJꎬZHANGJQꎬetal..Inversionofblacksoilorganicmattercontentwithfieldhyperspectralreflectancebasedoncontinuouswavelettransformation[J].Res.SoilWaterConnserv.ꎬ2018ꎬ25(2):46-52ꎬ57.1415期陈昊宇等:基于连续小波变换的土壤有机质含量高光谱反演[19]㊀范晓冰ꎬ鲁丽波ꎬ范海娇.托克托县耕地动态变化及驱动因子分析[J].内蒙古师范大学学报ꎬ2017ꎬ46(1):151-155.FANGXBꎬLULBꎬFANGHJꎬetal..AnanalysisofdynamicchangesanddrivingfactorsofarablelandinTogtohcounty[J].J.InnerMongoliaNorm.Univ.ꎬ2017ꎬ46(1):151-155. [20]㊀李政葵.内蒙古托克托县潜水与土壤中氟化物的分布规律及其相关性研究[D].呼和浩特:内蒙古大学ꎬ硕士学位论文ꎬ2015.LIKZ.DistributionandcorrelationoffluoridebetweenunconfinedwaterandsoilofInnerMongoliaTuoketuocounty[D].Hohhot:InnerMongoliaUniversityꎬMasterDissertationꎬ2015.[21]㊀VISCARRARAꎬMCGLYNNRNꎬMCBRATNEYAB.Determiningthecompositionofmineral ̄organicmixesusingUV ̄vis ̄NIRdiffusereflectancespectroscopy[J].Geodermaꎬ2007ꎬ137(1/2):70-82.[22]㊀吴倩ꎬ姜琦刚ꎬ史鹏飞ꎬ等.基于高光谱的土壤碳酸钙含量估算模型研究[J/OL].国土资源遥感ꎬ2020:[2020-09-16].https://kns.cnki.net/kcms/detail/11.2514.P.20200810.1518.004.html.WUQꎬJIANGQGꎬSHIPFꎬetal..Estimationofsoilcalciumcarbonatecontentbasedonhyperspectraldata[J/OL].RemoteSensingLandResourcesꎬ2020:[2020-09-16].https://kns.cnki.net/kcms/detail/11.2514.P.20200810.1518.004.html.[23]㊀张新乐ꎬ于滋洋ꎬ李厚萱ꎬ等.东北水稻叶片SPAD遥感光谱估算模型[J].中国农业大学学报ꎬ2020ꎬ25(1):66-75.ZHANGXLꎬYUZYꎬLIHXꎬetal..RemotesensingestimationmodelofSPADforriceleavesinNortheastChina[J].J.Chin.Agric.Univ.ꎬ2020ꎬ25(1):66-75.[24]㊀姚聪.基于卷积神经网络的耕层土壤有机质含量估测模型研究[D].山东泰安:山东农业大学ꎬ硕士学位论文ꎬ2020.YAOC.Studyonestmationmodelofsoilorganicmattercontentinplowedlayerbasedonconvolutionalneuralnetwork[D].ShandongTaian:ShandongAgriculturalUniversityꎬMasterDissertationꎬ2020.[25]㊀谢文.基于高光谱技术的森林土壤不同养分含量光谱特征及估测模型研究[D].南昌:江西农业大学ꎬ博士学文论文ꎬ2017.XIEW.Studyonspectralcharacteristicsandestimationmodelsofdifferentnutrientcontentsinforestsoilsbasedonhyperspectraltechnology[D].Nanchang:JiangsuAgriculturalUniversityꎬDoctorDissertationꎬ2017.(责任编辑:陈凌云)241中国农业科技导报23卷。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

由于土壤是一个由多种物质组成的复杂综合体,而土壤光谱 又是各种组份的综合表现,所以要准确定量描述土壤中某种 物质与光谱反射系数的关系有很大的难度,因此目前虽然光 谱在植被等方面的应用研究取得诸多重要进展,但光谱在土 壤方面的研究仍然停留在宏观及区域性研究的水平,一些定 量研究的实用及推广价值也相当有限。
式中,Y为有机质含量 (%), h为弓曲差。
弓 曲 差
图5 样品的散点图及回归曲线 (徐彬彬,1980) 湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
1.2 研究目的和意义
基于土壤高光谱变化规律,建立各类试验土壤相应SOM含 量高光谱预测模型,并探求所建模型之间的差异性。
人工模拟SOM梯度试验 探索在较为严格地控制土壤其他理化特性不变的条件,不同 土种SOM含量与土壤光谱特性之间的关系。 野外实地土壤样品SOM数据 研究建立自然条件下以土种为建模单元的SOM含量高光谱 预测模型。
数据处理 光谱微分 方差分析 多元线性回归分析
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
3、研究结果及分析
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
3.1 室内高光谱数据的不确定性及其影响因素分析
室内几何测试条件对土壤高光谱质量的影响
光源为50W的卤素灯、探头视场角为8°的条件下,以光源入 射角度15°、光源距离30 cm、探头距离15 cm作为室内几何测试 条件较为理想。
国内外已有的研究表明,土壤理化性质中,土壤水分、土壤 有机质含量、铁的氧化物(矿物组份)、土壤质地等诸多因 素对土壤光谱特性有明显影响。
土壤中一些理化特性只 有在 高光谱下才能得到反 映
图3、光谱分辨率对水铝石光谱的影响 (Cloutis,1996) 湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
已有的土壤有机质光谱研究
土壤反射光谱与土壤有机质含量呈显著负相关( Al-Abbas, A.H., 1972;Krishnan P.等 ,1980); 有机质中不同的成份光谱特性不一样(徐彬彬,1991; Demattê 等,1999;李震宇等, 1999; 牛灵安等, 2001 ) ; 有机质含量可以从土壤反射光谱中得到一定程度上的反应 ( Krishnan等, 1980;徐彬彬,戴昌达,1980;Henderson等, 1992;Sudduth等, 1993;Bon-Dor等,1995;Chang, C. W.等, 2002;McCarty等,2002 );
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
0.6
1 2
富里酸 胡敏酸
Humic acid Fulvic acid
反射系数 Reflectance
0.4
0.2
0.0 400
600
800
1000
波长 Wavelength(nm)
图4 胡敏酸和富里酸反射光谱曲线(徐彬彬等,1991)
图6
技术路线图 湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
2.3 试验设计
(1)室内光谱测试过程的不确定性研究
室内几何测试条件对高光谱质量的影响
表1 试验Leabharlann 素与水平设置土壤表面处理及粒径对室内土壤高光谱数据的可重复性影响 表面处理:表面刮平处理,表面稍压平处理 粒 径:5mm,2mm,1mm, 0.5mm, 0.28mm, 0.145mm
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
2、技术路线及试验设计
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
2.1 技术路线
应用研究得到的室内土壤高光谱数据获取的统一方法,结合 传统土壤理化特性分析技术,研究基于土种的SOM含量高光 谱预测模型及其差异性(技术路线如图6所示)。
土壤有机质含量高光谱预测模型 及其差异性研究
湖南农业大学资源环境学院
土壤有机质含量高光谱预测模型及其差异性研究
主要内容 1、研究背景、目的和意义 2、技术路线及试验设计 3、研究结果及分析 4、结论
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
1、研究背景、目的和意义
湖南农业大学资源环境
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
室内土壤高光谱测试 采用 ASD公司生产的ASDFieldSpec-FRTM型光谱仪测试, 该光谱仪光谱范围350~2500 nm,其中350~1000 nm 范围内光谱 采样间隔(波段宽)为1.4 nm,光谱分辨率为3 nm,在 1000~2500 nm范围内光谱采样间隔(波段宽)2.0 nm,光谱分辨 率为10 nm。测试具体操按研究中提出的规范进行。光谱反射系 数经专用平面白板标准化。
用过量双氧水氧化法弃除样品中的SOM,有机质被完全氧化后,过 量的过氧化氢用加热法排除(中国科学院南京土壤研究所, 1978)。 制备SOM梯度样品 用未弃SOM的样品和氧化后的样品按比例混合, 可按需要获取SOM含量范围在氧化后样品SOM含量和原土壤样品SOM含量 之间的任一SOM梯度样品。 SOM含量高于原土壤SOM含量的梯度样品, 用在原土壤样品中加入腐殖酸的方法获得。
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
有机质含量=
d log(1 / r0.6 2 3 6) / d0.6 2 3 6 K d log(1 / r0.5 6 4 4) / d 0.5644
式中,K为回归常数,r 为反射值,λ 为波长 ( Krishnan 等,1980) 。
土壤有机质含量高光谱预测模型及其差异性研究
1.1 研究背景
土壤在发育过程中,因为母质、气候、地形、生物等成土环境不同 和成土时间长短的差异,再加上人类活动的影响,其理化特性,如 机械组成、有机质含量、粘土矿物类型等有明显差别,光谱特性也 因此各有特点(如图1、图2)。
反射系数 Reflectance(%)
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
(2)基于土种的SOM含量高光谱预测模型及其差异性研究
人工模拟SOM梯度: 12种土壤(表3) 野 外 土 样: 两种土壤(青紫泥和红黄泥)
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
表2 被试土壤名称及相应发育母质
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
室内土壤高光谱数据采集规范
田间土壤样品采集 保证田间土壤样品采集时的剖面位置一致, 建议表层土样都采集0~15 cm深度范围内的土壤,且所采集的 样品应能尽量具备研究对象的代表性。 土壤样品制备 在土壤样品制备过程中,样品前处理条件(主 要包括土样粒径尺度、土样水分控制、土样表面处理等)应 一致,建议以1 mm作为土壤光谱测试样品的粒径尺度,风干 土壤样品测试前在40℃恒温条件下烘24小时,然后密闭冷却 24小时,测试时土壤表面进行稍稍压平处理。 测试条件的控制 测试时除遵循光谱测试指南的要求外,还 必须保证测试光照条件的一致性,建议在50W卤素灯、8°探 头角条件下,采用15°的光源照角度、30 cm的光源距离和15 cm探头距离作为室内土壤高光谱测试几何条件,同时还须保 证参考板表面与测试土壤样品的表面在同一平面上。
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
土壤有机质含量高光谱预测模型研究
土壤样品的采集、风干与制备
室内光谱测试过程的不确定性研究
SOM 含量高光谱预测模型研究
不同几何条件测试 条件对土壤光谱质 量的影响
样品不同表面 处理过程的光 谱差异性
不同粒径土样 品重复测试
利用较佳粒径样品 制作各个土类(或 土属)SOM 梯度
土壤编号 01 02 03 04 05 06 07 08 09 10 11 12 土 种(相当于系统分类中的土系或所属类别) 麻沙泥(铁渗水耕人为土) 麻沙菜园土(肥熟旱耕人为土) 花岗岩红壤(富铝湿润富铁土) 黄泥田(铁聚水耕人为土) 板岩红壤(粘化湿润富铁土) 黄菜园土(土垫旱耕人为土) 红黄泥(铁聚水耕人为土) 红菜园土(土垫旱耕人为土) 河潮泥(简育水耕人为土) 潮菜园土(肥熟旱耕人为土) 青紫泥(吴山系) 紫泥田(简育水耕人为土) 母质名称 花岗岩风化物 花岗岩风化物 花岗岩风化物 板岩风化物 板岩风化物 板岩风化物 第四纪红色粘土 第四纪红色粘土 河流冲积物 河流冲积物 海湖沉积物 紫色砂页岩风化物 采样地点 湖南长沙 湖南长沙 湖南长沙 湖南长沙 湖南长沙 湖南长沙 湖南长沙 湖南长沙 湖南长沙 湖南长沙 浙江湖洲 湖南长沙
湖南农业大学资源环境
土壤有机质含量高光谱预测模型及其差异性研究
(3)研究方法
土样的采集和制备 人工模拟SOM梯度土壤原始样品按常规方法采集(鲁如坤, 2000),野外实地SOM梯度土壤样品采样范围在所选采样点的 1m×1m面积内,采样深度0~15 cm。 土壤样品理化性质的测定 对所有土样的机械组成、SOM含量及速效N、P、K等理化性质 采用常规方法进行测试(鲁如坤,2000)。 SOM梯度样品的制备
野外实地土样
方差分析
光谱数据质量及可重复性分析
土样室内 光谱测试
土壤属 性测试
较佳的几何测 试条件
统一的表面处理 方法
较佳的土壤 样品粒径
多 元 统 计 分 析 方 法
统一室内土壤光谱测试方法及流程
基于人工 SOM 梯 度的 SOM 含量高 光谱参数模型
基于野外实地土样 的 SOM 含量高光 谱参数模型
由于受诸多因素的影响,土壤光谱研究具有风险大、难度高 的特点,为此,研究选择土壤众多理化特性中对土壤光谱反 射影响较大,同时又是肥力重要指标之一的土壤有机质 (SOM)作为研究对象,应用高光谱数据获取手段,对以土 种为单元的SOM含量高光谱预测模型作探索性的研究。
相关文档
最新文档