aes 应用和环境土壤学 1687 - 7675 1687 - 7667 Hindawi出版公司 971252年 10.1155 / 2012/971252 971252年 研究文章 基于特征的比较高,请回归预测技术的三个南非退化生态系统中土壤成分 拜耳 安妮塔 1 巴赫曼 马丁 1 穆勒 安德烈亚斯 1 考夫曼 赫尔曼 2 Ben-Dor Eyal 1 部门的地表 德国遥感数据中心(过程) 德国航空航天中心(DLR) Oberpfaffenhofen 82234年Weßling 德国 corine.dfd.dlr.de 2 遥感部分 大地测量学和遥感 德国地球科学研究中心(德国) Telegrafenberg 14473年波茨坦 德国 gfz-potsdam.de 2012年 17 9 2012年 2012年 16 02 2012年 20. 04 2012年 21 05年 2012年 2012年 版权©2012 Anita拜耳等。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

选择土壤成分的准确评估可以提供有价值的指标来识别和监控土地变化加上退化在半干旱地区频繁的现象。两种方法的量化土壤有机碳、铁氧化物和粘土含量根据现场和实验室自然表面的光谱测试。(1)物理方法基于光谱吸收特性分析。对于每个土壤组成,一组诊断光谱特征选择和与化学参考数据由多个线性回归(高)技术。(2)偏最小二乘回归(PLS)作为一个专门统计多元方法适用于比较。回归模型开发基于大量的地面参考数据收集的163个采样站点的丛林生物群系,南非,土地变化观察到由于密集的过度放牧。评估其预测性能的方法和意义关于未来的量化土壤成分在大面积使用成像光谱。

 1。介绍

土壤上层的地球表面是最重要的能源和营养流层植被的发展,因此所需的关键重要性进行景观分析。描述生态系统的土壤条件及其时空变化土壤健康的重要指标,尤其是在农业生态系统与作物生产直接相关。在半干旱地区,加上退化和水土流失土地覆盖变化频繁的现象可能是长期管理实践的结果或可能与气候变化有关。特别是在土壤碳的消耗库存以土壤退化和侵蚀为主( 1)和直接原因不仅环保而且经济问题。半干旱的亚热带丛林生物群系的南非东开普省土地变化观察到由于几十年的过度放牧山羊。这造成了独特的生态系统改变从茂密的灌木地植被的固碳率高,周边土壤系统开放的研制。土壤化学和物理属性可以作为示踪剂,评估和监测等现象。然而,映射的空间分布和时间发展有限使用传统的土壤分析,因为这需要密集采样和分析工作。

尽管如此,字段和成像光谱提供了一个——高效的工具,用于映射的选择对大面积土壤化学和物理属性。土壤成分可以根据诊断确定固有的土壤的反射光谱吸收特征。一般来说,量化诊断光谱特性的直接应用目的是有限的许多特性受到(1)的重叠光谱特性和(2)附加其他的影响,主要是物理,土壤性质(如表面粗糙度)影响光谱反射率无意义的。

介绍了两种方法的量化土壤有机碳、铁氧化物和粘土含量基于实地和实验室光谱。首先是物理方法基于光谱吸收特性分析。对于每个土壤组成,一组独特的光谱特性在光谱的选择根据他们的存在和他们各自的文档文献。他们是与化学有关的地面通过多元线性回归技术参考信息。这种方法建立了利用几个光谱特性和特点以及各种属性描述这些收集最物理特征所提供的信息内容。测试如果结合几个光谱特性的优势,虽然每个光谱特性可能影响,将导致预测模型及统计技术的性能。第二种方法,常用的偏最小二乘回归(PLS)作为一个专门统计多元方法适用于比较。请技术是否适合映射的几个土壤成分是众所周知的(例如, 2])。研究了两种方法的预测性能和意义三个光谱数据集以现场和实验室和两个光谱分辨率(全ASD HyMap传感器的分辨率和光谱分辨率)。本研究的目的是允许一个健壮的量化碳库存半干旱土壤和另外提供指标与土地退化的生态系统功能。

在下一步中,这些方法将被应用为未来大规模量化土壤成分南非研究网站的灌木丛生态区。高光谱图像(大约320公里2)HyMap机载传感器( 3)在2009年收购了这个区域的,加上一个广泛的光谱和化学地面参考数据的集合(163参考图)。在灌木丛生态区,准确的空间信息可以用来量化的后果不可持续的土地管理技术和检测侵蚀或退化地区的初级阶段。

 2。状态的艺术土壤光谱

许多研究存在的上部土层描述基于光谱数据从实验室或现场环境。例子包括土壤有机碳(例如, 4- - - - - - 7)、土壤铁氧化物(如[ 8- - - - - - 10]),土壤质地与粘土含量(例如, 7, 11, 12])。

与越来越多的最后几年里,机载光谱仪数据用于预测土壤成分。直接转移的方法在实验室或现场数据开发的机载扫描仪是有限的技术方面根据传感器系统等空间和光谱分辨率较低,信噪比增加,大气影响,双向反射效果,在像素光谱混合。最近的总结主要研究利用成像光谱研究土壤特性可以在Ben-Dor et al。 13]。

Viscarra Rossel et al。 2)大纲的潜力不同的谱域结合不同方法的推导定量土壤利用土壤光谱信息。偏最小二乘回归(PLS)为多元校正过程是一个行之有效的方法对这些目的,通常是用于以前的研究(例如, 14- - - - - - 16])。请统计导出适用于光谱特征之间的关系来描述土壤化学和光谱信息。然而,请模型的意义和可转移性建立在本地观测区域范围内往往是有限的,因为他们的高度的适应(见,例如, 15])。允许一个映射为生态系统的土壤成分的区域范围,需要更健壮的算法。在先前的研究中,方法直接应用光谱吸收特性或乐队指数被认为是更合适的,因为他们是基于物理(例如, 17])。这些方法看起来更健壮和允许转移到地区相似的环境条件。在以下论文中,一个特殊的焦点被设置为研究应用土壤光谱特征参数量化,特别是基于机载传感器。

土壤有机碳的稳定的重要性和生育能力及其对全球碳循环的贡献导致全球倡议量化生态系统碳储量及其变化。因此,有非常强烈的努力为这些目的使用光谱导致各种研究的量化处理土壤有机碳(例如, 14, 15, 18, 19])。旁边的多元技术,特别是估计土壤有机碳是基于全面降低反射率在可见光及近红外光谱区域(见[ 4])也使用光谱指数(例如, 4, 6])。

以类似的方式,几项研究已经应用吸收特征参数的量化铁氧化物含量。里希特et al。 20.]所以乐队使用900 nm的深度特性和评价光谱变量土壤质地对预测精度的影响,实验室的设置。铁氧化物的浓缩Israelian沙丘(rubification)被映射Ben-Dor et al。 21使用属于接近机载传感器),红色指数的应用程序作为一个iron-related乐队指数。

粘土矿物的大规模映射,Chabrillat et al。 22)应用的不同形状2200纳米蒙脱石的特性,伊利石、高岭石和检测和地图的空间分布从空中AVIRIS HyMap数据科罗拉多州,美国,而戈麦斯et al。 23和也Lagacherie et al。 24)成功应用的连续删除带深度量化粘土含量2200纳米特性。

 3所示。研究区和数据收集

亚热带丛林是一个独特的植被类型位于半干旱的山谷,覆盖大约17%的南非的东开普省( 25]。其中,70%以上的植被单位视为“适度”“严重”退化 26]。丛林是一个不寻常的植被类型为半干旱环境从一个几乎完全覆盖茂密的植被存在于原始的条件。丛林植被组成组件的多汁的灌木 Portulacaria身兼( 27),占大比例的碳生物量和周边的土壤。密集的山羊养殖自1900年代初以来导致的损失 身兼p .,宽的变换部分丛林植被开放研制系统,这是伴随着严重的生物多样性的丧失和生态系统碳储量 25]。栅栏农场之间不同的利用显示原始和改变了植被在直接附近(图 1)。

高对比度土地覆盖特征发现栅栏牧场和农场游戏之间的界线在东开普省,南非。牧场农场左边右边是高度退化而土地只有轻微影响,几乎是原始丛林灌木植被。持续土地利用的后果可以检测和监控使用。

研究区位于城市附近的伊丽莎白港(33.0°S / 25.3°E;参见图 2)。这是选择在灌木丛生态区的程度 75年 × 4 公里,包括最高的方差各种错综复杂的植被类。它将两部分分开的一个国家公园。地形高程变化从300年到930年的北部,130年到400年南方部分。研究区与月度温暖半干旱气候冬季平均气温从13°C到25°C(在夏季(数据 28])。平均年降水量200 - 400毫米。底层泥岩和砂岩的发展导致了肥沃的土壤和桑迪在表面(主要是始成土和淋溶土),(自己的调查和 29日])。

HyMap图像的详细视图的两个横断面萨默塞特郡东部和西部的柯克伍德在东开普省附近,南非,显示采样站点的分布。每个样部分由三个HyMap飞行线路。

 3.1。现场取样

163年两个抽样活动自然、非农网站采样为地面真理,96年9月2009年6月和67年/ 2009年10月。故事情节被随机分布在地形单元和研究区域的两个部分,如图 2。大部分的采集标本在裸露的土壤是占主导地位的土地覆盖的地区分数。网站的数量和可变性采样会导致完全覆盖土壤的方差。每个点被选为代表的典型条件调查周围的区域。对于每一个情节,一般的地址信息(GPS坐标、地形信息、土地利用等),现场接触土壤和土壤的光谱样本收集。额外的信息关于植被覆盖密度和类型,以及表面条件评估中心的点 1 × 1 米)及其周边( 10 × 10 米)。土壤样本的代表点包括顶层最大深度约1厘米,仅限于地壳如果物理表面土壤结皮出席。38块,野外光谱测量不可能是由于阴条件;因此土壤野外光谱采样站点的只有125。

机载HyMap成像光谱仪数据获得2009年10月在研究区地面分辨率为3.3 m。这些数据用于适当的光谱重采样点的光谱数据集而不是本研究的范围,因此这里描述不详细。

 3.2。实验室化学分析

土壤样品化学分析有机碳和亚硫酸氢可榨出的铁使用Walkley-Black分析[ 30.),分别citrate-bicarbonate-dithionite提取方法(CBD, [ 31日])。粒径在5个分数决定吸管法( 32]。分析显示桑迪和肥沃的土壤有机碳含量普遍偏低。没有一个直方图分布显示了正态分布(图 3)。所有参数,低浓度表示。对有机碳,只有两个样品超过3.5%,显示高5.9%和5.7的内容。这两个样本在不同寻常的网站,因为样本是在一个地区受到放牧动物的粪便,在另一个示例网站部分草层覆盖裸露的土壤表面,可能会导致增加土壤碳的输入。铁氧化物,它们主要是底层的风化过程的结果砂-泥岩范围高达10.6%。实地调查表明,粘土和砂的分布上覆土壤底层的材料直接相关。粘土和砂浓度的增加,泥岩、砂岩,分别发生接近水面。所有样本用于进一步分析,尽管一些样品的化学内容不同于其他人群的分布(见图 3)。

直方图分布统计参数164表土样品收集在灌木丛生态区,南非,分析土壤有机碳、铁氧化物和粘土含量。

 3.3。光谱参考数据

测量协议导致三个光谱数据集字段和实验室环境,(1)insitu公司现场原状土表面的光谱,(2)裸露的土壤光谱,在最终呈现小石块的封面,和(3)实验室光谱。

spectroradiometric测量的便携式FieldSpec Pro 2光谱仪(分析光谱设备,Inc .)使用。进行了测量与裸露的纤维(FOV 25°),直接转换为使用Spectralon反射面板作为参考。测量身高在1.15米保持不变,导致足迹直径大约50厘米。场所的安静的裸露的土壤部分非常小,另一个身高0.5米的应用。如果小石子(毫米厘米范围)覆盖表面,改变表面的测量进行了第一只(insitu公司),然后从裸露的土壤表面的石头盖小心地删除。石头覆盖大约一半的网站很低(< 3%)。尽管如此,特别是在地区土壤浅和发展在泥岩、它和范围的风险高50%左右。

样本干后,溶解聚合物,和土壤渗分数小于2毫米了光谱测量在实验室。恒照明条件由两个石英卤素灯(300 W) 30°的天顶角和谱仪探针和样品表面之间的距离15厘米。5测量后被裸露的纤维,再由90°和测量每个样本被降低土壤表面的光照效果。使用Savitzky-Golay平均场和实验室光谱平滑滤波器( 33]。此外,光谱测量领域,范围的大气水乐队插值从1345年到1445 nm和1800 - 1960 nm)。

自量化方法开发应用到高光谱图像,每个数据集用于完整的ASD分辨率为2087光谱实验室数据和1829年乐队的乐队现场数据集,分别为(382 - 2468 nm, 1海里重采样),另外重新取样116校准HyMap图像的光谱特征选择乐队(456 - 2455纳米,光谱分辨率13 - 17海里)。

不同的测量设置反映在统计(图的谱库 4)。的平均反射率光谱测量在实验室高于两个区域的数据集,因为事先预处理应用样品(均质化和筛分)防止阴影效应可能是由于土壤聚合或外来物质(图 4(一))。控制的实验室条件下导致非常低的标准差。尽管如此,实地测量的光谱差异影响拉伸,不特定的土壤成分的特征,如表面状况和照明效果(数字 4 (b)和 4 (c))。insitu公司实地测量,平均土壤反射率降低小石块上覆土壤表面。详细的统计数据显示,石头的存在还介绍了进一步变化的测量,因此,标准差是高于相应的裸露的土壤实地测量。裸露的土壤野外光谱特别是现在大极端最大和最小光谱(图 4 (c))。土壤物理外壳存在,光滑的土壤表面结合减少了土壤表面的粒子的大小增加土壤反射率。因此单个样品的反射率测量现场环境中超过了相应的实验室光谱反射率。

比较平均光谱insitu公司领域的三个光谱库(a)(固体),字段裸露的土壤(虚线)和实验室测量(虚线)。统计概述125年insitu公司现场土壤测量(b), 125场裸露的土壤(c),与163例实验室光谱(d)(固体),积极的和消极的标准差(虚线),最小和最大(虚线)。所有光谱给出了完整的自闭症光谱分辨率。

 4所示。方法

应用两种方法建立定量关系的化学内容调查土壤成分和光谱信号。是一个基于物理模型的方法,光谱特性分析是加上多元线性回归技术,而方法B应用偏最小二乘回归作为一个专门多元统计方法。75%的数据(94个样本数据集,123个样本进行实验室数据集)是由随机分层抽样选择基于每个土壤成分的化学参考培训和用于模型校准,而剩下的25%(31岁,分别地。,40个样品作为测试集验证模型。训练集和测试集的分布从而代表化学引用值发生在南非研究区域。即使测试样品不用于模型校准,他们不能被视为完全独立于训练集由于空间距离的测量。

4.1。方法一:多元线性回归的光谱特征参数

这种方法适用于一系列的光谱特性,发现在之前的文献特征对土壤有机碳、铁氧化物和粘土含量对随后的多元线性回归分析。因为特定的光谱属性相关的物理吸收过程,这种方法被认为是一个物理模型。这是利用开发的几个光谱特性和特征的组合。通过参数化的光谱特性,利用现有的基于物理信息。

以下4.4.1。选择的光谱特性

诊断光谱特性的选择是基于他们的存在在光谱和著名的文学作品,例如,( 7, 8, 34- - - - - - 40]。从现有的信息只选择最重要和强大的功能来提高算法的鲁棒性。功能独特的土壤成分是首选,但并不是所有的光谱特性可以满足这个需求。表 1列出了光谱特性,选择用于本研究的总结和引用先前的研究,这些特性被描述为丰富土壤或nonsynthetic矿物粉末。

光谱特性,用于界定土壤有机碳、铁氧化物和粘土含量方法a .引用先前的研究描述这些光谱特性。

特点和类型 波长 赋值 参考
C org (我)1730海里(AF)CR1:1600 - 1815海里 λ 马克斯 2 :没有定义 1650/1669 nm 2 υ3芳香的CH拉伸 ( 5, 7]
≈1700海里 2 υCH延伸 ( 34]
1706/1754 nm 2 υCH,烷基紧身上衣 ( 7]
1706海里 4 υ脂肪族的CH拉伸 ( 7]
1726海里 2 υ脂肪族的CH拉伸 ( 5]
1761/1769 nm 2 υ脂肪族的CH拉伸 ( 5]
1730 - 1852纳米 4 υ甲基CH ( 7]
(2)2330海里(AF)克雷格:2240 - 2410 nm λ 马克斯 :2330海里 2275/2279 nm 3 υCH的2,CH3 ( 5, 7]
≈2300海里 υ+ υ4CH延伸 ( 34]
2307 - 2469纳米 3 υ甲基CH拉伸 ( 7]
2309海里 3 υ脂肪族的CH拉伸 ( 5]
2331海里 3 υCH的2,首席运营官 ( 5]
2337/2386 nm 3 υ首席运营官,CH3 ( 5]
2347海里 3 υ脂肪族的CH拉伸 ( 5]
2381海里 3 υ的公司的碳水化合物 ( 7]
(3)450 - 740 nm(高频) 降低反射率的可见范围 (例如,( 4- - - - - - 6, 35])
(四)1460 - 1750海里(高频) 减少反射的短波红外范围附近 (例如,( 4, 6, 35])
铁氧化物 (我)550海里(AF)克雷格:450 - 680 nm λ 马克斯 :550海里 ≈490海里 等5乐队的铁3 + ( 7, 36]
503海里 针铁矿 ( 37]
≈510海里 ET的铁2 + ( 36, 38]
529海里 ET的赤铁矿 ( 7]
535海里 赤铁矿 ( 37]
≈550海里 ET的铁2 +、赤铁矿 (例如,( 36, 38, 39])
(2)700海里(AF)克雷格:580 - 800 nm λ 马克斯 :700海里 650海里 等群赤铁矿和针铁矿 ( 7, 39]
665海里 等的针铁矿 ( 37]
700海里 4 υ的地 ( 7]
≈700海里 ET的铁3 + (例如,( 35, 36, 38])
(3)900海里(AF)克雷格:750 - 1300 nm λ 马克斯 :870海里 ≈850 - 870纳米 ET的赤铁矿 ( 39, 40]
860海里 ET的铁2 +和菲3 + ( 34]
868海里 赤铁矿 ( 37]
≈870海里 ET的铁3 + ( 35, 36, 38]
884海里 ET的赤铁矿 ( 7]
≈900海里 过渡的铁2 +和菲3 + (例如,( 34, 38])
900 - 930纳米 等的针铁矿 ( 39]
920海里 等的针铁矿 ( 7]
≈930海里 3 υo - h键发出的,针铁矿 ( 7, 37, 40]
940海里 υ+ υo - h键发出的 ( 7]
1000 - 1100纳米 ET的铁2 + (例如,( 35, 36])
1025海里 ET的铁3 + ( 8]
1075海里 ET的铁3 + ( 8]
(四)550 - 590海里(CF) 减少反射蓝色波长范围对紫外线 (例如,( 34, 38, 39])
(v) 450 - 750 nm(高频) 降低反射率的可见范围 (例如,( 34, 38, 39])
粘土 (我)2200海里(AF)克雷格:2100 - 2290 nm λ 马克斯 :2206海里 2160/216 + 2208/2209 nm δ6高岭石AlOH弯曲的紧身上衣 υ+ υ哦的延伸 ( 7, 37, 50]
≈2200海里 δAlOH弯曲和 υ+ υ哦的蒙脱石和伊利石 (例如,( 35, 40, 50, 51])
2200/2204 nm 蒙脱石 ( 37]
2206海里 υ哦的蒙脱石和伊利石 ( 7]
2208海里 δ高岭石AlOH弯曲的紧身上衣 ( 7]
2216海里 伊利石 ( 37]
2230海里 δ蒙脱石的AlOH弯曲 ( 7]
(2)2340海里(AF)克雷格:2270 - 2410 nm λ 马克斯 :2340海里 2308/2312 nm 高岭石 ( 37]
2336海里 伊利石 ( 37]
2340海里 υ哦的伊利石 (例如,( 7, 22, 50])
2372/2376 nm 高岭石 ( 37]
(3)450 - 700 nm(高频) 增加反射率在可见范围内 (例如,( 34, 42])
(四)1460 - 1750海里(高频) 增加反射的短波红外范围附近 (例如,( 34, 42])

1克雷格:范围连续执行删除;2 λ 马克斯 :最大吸收波长主要发现在文学;3 υ:泛音吸收带;4 υ+ υ:乐队的基本组合,泛音吸收;5外星人:电子跃迁的乐队;6 δ:基本的吸收带。

选中的光谱特征分为三种特性:吸收特性(AF)的特征光谱曲线(CF)和特性的连续光谱的凸包(高频)。与这三种类型假设中描述的所有特定的光谱特征覆盖以前的研究作为重要的土壤有机碳、铁氧化物和粘土含量。

光谱检测和评估的土壤有机碳具有挑战性,因为许多有机材料所固有光谱特性通常是弱匹配的复杂化学有机物。在吸收以往的研究发现有机物分子和组合的特点,选择波长范围在1730和2330海里,因为他们被认为是最独特的先前的研究。尤其是吸收约1730纳米有机碳是独一无二的,因为在这个波长区域没有其他共同矿物质的吸收存在于土壤和岩石发生。在这两个波长区域,几个第二和第三的泛音吸收的功能组分配到纤维素、木质素、淀粉、果胶、葡聚糖、蛋白质、蜡和腐殖酸结合土壤有机质的组成部分重叠和导致明显的吸收(例如, 5, 34])。organic-carbon-rich土壤的一般减少反射,例如,可见由于腐殖酸的黑暗,包括光谱连续的一个特征。几项研究采用不同的波长区域为了描述该属性(例如, 5, 6])或测试各种光谱指数的性能检测土壤有机碳(例如, 4])。在这项研究中,这两个范围可见地区450至740海里,在不久的短波红外介于1460和1750 nm之间用来描述特征减少反射产生的土壤有机碳的存在(高频)。

两种有机物的选择特征模糊,但仍可以应用自当地条件研究区域排除任何影响。首先,这适用于附近的吸收2330海里,这也是固有的碳酸盐矿物,但由于没有碳酸盐岩基岩在研究区沉积物和碳酸盐沉淀只出现在本地和主要在更深的土壤层。第二,除了土壤有机碳的影响,总体反射率也受土壤含水量和粒度分布的影响。然而,土壤水分的影响可以忽略,因为所有土壤干燥的日子和时间前的调查。这也是由当地气象数据显示( 41]。晶粒尺寸的影响对整体反射和特性接近碳特性还包括在2330 nm作为粘土含量的测定的光谱特性,但它是不可能解决每个成分的比例。

测定土壤铁氧化物,我们应用强大的吸收发生在可见范围大约在550年,700 nm,接近900海里。他们由于电子跃迁的bi -和三价铁离子化合物最常见的铁氧化物针铁矿和赤铁矿(例如,[ 34, 35, 39])。电子跃迁乐队通常表现为广泛的吸收。减少反射蓝色对紫外线的地区的广泛吸收带紫外线范围是由不同的作者描述(例如, 38, 39]),包括光谱曲线的两个特性:首先,反射率的变化在蓝色波长范围光谱曲线的斜率之间的550和590海里(CF)第二,连续体的形状包含整个VIS范围(450到740纳米,高频)接曲线的形状与减少550纳米吸收的影响。

粘土含量的光谱测定使用知名AlOH吸收特性约2200海里,是内在的粘土矿物,和一个小羟基功能约2340海里。因为特性源于基本的官能团,泛音吸收固有的粘土矿物,它们的存在是明显的,窄的吸收。随着伊利石和蒙脱石作为主要粘土矿物在研究区,2206年被选为中心的波长2200 nm的吸收。此外,粒度分布的光谱效应作为土壤粘粒含量的间接指标。越来越内容的细粒度大小的影响导致反射率的增加(例如, 34, 42])分析了形状和意味着反射光谱曲线的活力和近红外光谱/短波红外成像范围(高频)。

 4.1.2。特征参数化

在特性参数,分析了光谱数据集选择的三种类型的光谱特征进行了介绍和转移到数值参数描述光谱的形状特性(见图 5)。这些光谱特征变量用于随后的回归分析。表 1应用光谱特性的列表参数用于参数化。

参数化变量描述的三种类型的光谱特性(AF、CF和高频)用于测定土壤有机碳、铁氧化物和粘土。实线:反射率曲线,虚线:连续的反射率曲线。

光谱吸收特性(AF)从连续删除参数化反射光谱分离从整体反射和允许相互比较的趋势(见图 5(一个))。连续去除(见[ 43])计算单独定义区间的特性映射(CR在表 1)。房颤是通过以下六个变量,描述相似的光谱吸收特性分析由不同的作者描述(例如,[ 37, 44, 45):深度( d 马克斯 )和波长( λ d 马克斯 )的最大吸收,吸收深度应根据各自的文学特征波长( d λ 点燃 )、特征宽度( w ),两者之间的距离确定功能的肩膀( 年代 左 / 正确的 之间的区域),归一化连续和光谱曲线( 一个 )和不对称( 作为 = 一个 左 / 一个 正确的 )。注意,对于附近的房颤1730海里,只有五个变量计算因为没有中心波长( λ 点燃 )自定义此功能由几种弱吸收波长重叠在这个地区(见表 1和 2)。描述重要的光谱曲线的形状在一个波长范围定义,介绍了两个特点。(1)曲线的特性(CF)描述反射率的变化发生在一个特定的波长范围(见图 5 (b)),例如,引发的一个强大的铁氧化物吸收紫外线。他们只有通过特征光谱曲线的斜率( 年代 从一条线),计算在给定的波长范围(见表 1)。(2)凸包的特点(高频)描述土壤成分的影响广泛的光谱范围,不能产生一个明显的吸收(见图 5 (c))。VIS地区之一(450到740海里)和第二个短波红外成像范围(1460 - 1750海里)。这些波长区域的两个高频用于多个土壤成分的测定。高频特性参数化意味着斜率( 年代 )和平均反射率( r )光谱的凸包的波长范围,这样他们只描述定义不同的光谱形状没有影响当地的吸收。这两个变量的高频特性计算通过使用一条线。

数量的光谱特性和总数计算变量用于关键土壤成分的测定。

吸收特性(6 var。) 曲线特性(1 var。) 船体的特性(2 var。) 变量的总数
C org 21 0 2 15
铁 3 1 1 21
粘土 2 0 2 16

1AF附近1730海里,只有五个变量计算(见文本)。

的参数化 n 光谱和 米 导致光谱变量 米 × n 大小的矩阵。表 2显示数量的光谱特性进行分析测定的土壤有机碳、铁氧化物、粘土、以及随之而来的数字变量描述光谱形状和随后用于回归分析模型参数的内容。

 4.1.3。MLR模型校准,即分析

多元线性回归分析是用于建立光谱变量之间的函数关系和化学参考数据。图 6总结关键步骤工作流程应用的基于功能的回归方法。高钙进行单独为每个土壤成分使用地面参考数据收集领域的活动。产生的群光谱变量参数化的最初选择的光谱特性对三个参数进行检查,以确保他们适合高分析:(1)每个变量的值的归一化标准差(标准差与均值)必须高于0.001排除非常小的变量只会导致适应不良的统计回归分析。如此低的标准差,例如,将达到的深度吸收功能,只有弱建立或只发生在一些光谱。(2)光谱变量导致冗余值可能是由于降低了光谱分辨率(例如, d 马克斯 和 d λ 点燃 ),(3)变量的方差分析确定变量与一个高度不变的部分,可能是由上级的控制因素(例如,传感器带位置)。这些类型的光谱变量不适合回归分析并排除在分析之外。

基于功能的高钙和PLS回归分析工作流程。

随后计算特征变量标准化允许回归系数的比较和排名。与化学相关的参考价值考虑土壤成分,进行多元回归分析导致最初的关系没有进一步干预。分析交叉验证进行分析的光谱变量集合意义和为了检查样本人口的离群值。光谱变量回归分析被认为是无关紧要的,如果一个回归系数的绝对值的平均值小于两倍标准偏差(见[ 46)确定分析交叉验证。光谱变量发现微不足道的被排除在进一步分析,因为他们为模型开发和提供无附加价值可能提供统计适应。一个示例 我 被确定为离群值的绝对偏差 R 我 2 从所有人的意思 R 我 2 高于两倍标准差 R 我 2 确定分析交叉验证。样品确认为异常值可以通过手动交互从人口在这一步如果一个合适的理由是著名的(例如,抽样问题或分析)。最后建立多元线性回归模型的基础上,显著的特征变量,包括所有样本。培训土壤参数预测模型校准光谱数据的光谱分辨率。模型在测试集验证评估来评估每个模型的可预测性。

 4.2。方法B:偏最小二乘回归

相比之下,校准模型都是建立在相同的光谱数据集使用偏最小二乘回归,因为它是一种行之有效的最优化技术,经常应用于土壤光谱。请基于投影的预测( x )和反应( y )变量为一组潜在的变量(或请因素)和相应的分数,减少数据的维数,同时最大化之间的协方差 x 和 y 变量。的详细描述请技术是由( 47]。请建模软件应用(就算 48]。

光谱数据的预处理与一个或两个总共11操作方法的结合(转换反射( R ) 日志 ( 1 / R ) 5光散射和基线修正如乘法散射校正(MSC)、标准正态变量变换(SNV),小波消除趋势,等等,第一和二阶导数的计算,意味着定心,方差数据的规模,后两种的组合)。每个预处理设置的性能评价分析交叉验证(简历)。对于每一个设置,最优数量的潜在变量用于建模的变化了均方根误差(均方根)和Akaike信息准则(AIC)的潜在变量的函数。请确定因素的最优数量内的局部最小值的均方根和AIC稳定这两个因素的趋势。一个预处理设置结合其相应的最优数量的请因素表现最好的交叉验证被选中,随后用于模型校准基于训练集的观察。校准模型进一步应用于测试集验证目的(见图 6)。

如果几个预处理设置提供类似的简历精度,然后为每个设置模型进行了标定和验证。精度最高的PLS模型校准和验证,以及它们之间的最小差异被选中。这种方式最重要和健壮请预测模型的检索和过度拟合模型是预防。

 5。结果与讨论

对于每个土壤组成、模型使用两种建模方法,建立三个光谱数据集,和两个光谱分辨率,导致12模型组成。所有样品,虽然分成训练集和测试集,是申请建模。为每个方法和模型性能评估数据集和比较基于模型的相关系数( R 2 )预测和测量成分、均方根误差(均方根),和性能比偏差(RPD)。RPD被定义为参考样本的标准差的比值除以RMS。这样做是对校准和验证使用相应的数据集的训练数据样本(94 123场和实验室测量)和测试数据(31/40)样品(见表 3和 5)。

校准和验证精度方法多重线性回归的光谱参数。模型进一步应用大规模预测土壤成分基于HyMap意象突出显示。

光谱数据集 不。光谱变量 校准(941/ 1232样品) 验证(311/ 402样品)
R 卡尔 2 RMS卡尔 RPD卡尔 R 瓦尔 2 RMS瓦尔 RPD瓦尔
自闭症光谱分辨率 (1)原位 12日16 0.77 0.47 2.09 0.49 0.54 1.26
C org (2)裸露的土壤 14日16 0.79 0.44 2.22 0.56 0.48 1.42
(3)实验室 13日16 0.75 0.47 2.01 0.74 0.36 1.93
(1)原位 19日21 0.63 0.80 1.64 0.24 0.90 1.07
铁氧化物 (2)裸露的土壤 15日21 0.64 0.78 1.68 0.21 0.93 1.04
(3)实验室 19日21 0.75 0.76 2.01 0.23 1.63 0.69
(1)原位 14日16 0.31 4.12 1.20 0.01 5.13 0.88
粘土 (2)裸露的土壤 11日16 0.21 4.39 1.13 0.03 4.80 0.94
(3)实验室 16日16 0.23 4.60 1.14 0.05 4.52 1.00
HyMap光谱分辨率 (1)原位 13日16 0.79 0.44 2.19 0.51 0.54 1.27
C org (2)裸露的土壤 12日16 0.81 0.43 2.28 0.62 0.43 1.57
(3)实验室 15日16 0.77 0.46 2.08 0.77 0.35 2.00
(1)原位 18日21 0.62 0.80 1.64 0.17 1.04 0.93
铁氧化物 (2)裸露的土壤 19日21 0.66 0.76 1.73 0.26 0.93 1.03
(3)实验室 20日,21日 0.73 0.79 1.94 0.25 1.58 0.71
(1)原位 14日16 0.28 4.20 1.18 0.00 5.05 0.89
粘土 (2)裸露的土壤 14日16 0.17 4.49 1.11 0.01 4.91 0.92
(3)实验室 12日16 0.25 4.53 1.16 0.05 4.62 0.98

1训练集和测试集的样本数量字段数据集。2训练集和测试集的样本进行实验室数据集。

善良的预测评估的定性分类Chang et al。 9]。他们建议 R 2 大于0.80和RPD值大于2.0作为优秀的指标预测模型。 R 2 在0.50和0.80之间和RPD值在2.0和1.4之间,被认为是中等品质的模型为定量预测在大多数应用程序是有用的。模型与 R 2 低于0.50和RPD低于1.4排名不是可用的。

5.1。方法一:多元线性回归的光谱特征参数

建立校正模型的细节给出了表使用方法 3。光谱变量发现微不足道的数量和排除在最终在交叉验证模型发展取决于它的变化。因此,重要的变量的数量变化对所有模型建立了一个特定的参数。表 4给五个光谱变量,显示每一个基于特征的回归方程系数最高回归模型,因此是最重要的发展模式。他们的影响力决定基于回归系数和给定的%总结所有回归系数的绝对值。签署的负面影响表示消极的回归系数。回归系数的绝对值,从而影响由此决定的,主要取决于参与的数量和特征光谱变量。因此,不能指望他们的绝对值为不同的模式是相同的。尽管如此,他们的范围是重要的和用于识别最重要的光谱变量。

影响( 我 )的五个最重要的光谱变量的回归方程,基于回归模型的方法。

光谱数据集 1级 等级2 等级3 等级4 等级5
我 (%) 变量 我 (%) 变量 我 (%) 变量 我 (%) 变量 我 (%) 变量
自闭症光谱分辨率 (1)原位 36.6 r 高频 SWIR1 −24.0 年代 高频 力 −10.9 r 高频 力 −6.0 年代 高频 SWIR1 −5.9 一个 AF2330
C org (2)裸露的土壤 30.5 r 高频 SWIR1 −21.1 年代 高频 力 −8.4 r 高频 力 8.0 d 马克斯 AF1730 −7.7 d 马克斯 AF2330
(3)实验室 27.6 一个 AF1730 16.9 r 高频 SWIR1 −14.1 d 马克斯 AF1730 −11.8 r 高频 力 −9.0 年代 高频 力
(1)原位 −23.7 一个 AF550 21.8 d 马克斯 AF550 −12.0 一个 AF900 4.9 d λ 点燃 AF900 4.9 d 马克斯 AF700
铁氧化物 (2)裸露的土壤 −17.6 一个 AF900 10.6 d λ 点燃 AF900 −10.4 d λ 点燃 AF550 9.9 d 马克斯 AF700 9.8 d 马克斯 AF900
(3)实验室 29.2 d 马克斯 AF900 −14.7 d λ 点燃 AF900 −10.3 一个 AF900 6.7 年代 高频 力 −5.6 d 马克斯 AF550
(1)原位 29.3 一个 AF2200 −29.1 d λ 点燃 AF2200 7.9 r 高频 SWIR1 5。6 r 高频 力 4.9 λ 距离 AF2200
粘土 (2)裸露的土壤 −26.4 d λ 点燃 AF2200 23.2 一个 AF2200 7.8 d 马克斯 AF2200 7.7 r 高频 SWIR1 7.5 r 高频 力
(3)实验室 −28.3 d 马克斯 AF2200 18.8 d λ 点燃 AF2200 −12.4 d λ 点燃 AF2340 12.0 一个 AF2200 7.5 d 马克斯 AF2340
HyMap光谱分辨率 (1)原位 34.5 r 高频 SWIR1 −20.8 年代 高频 力 −13.1 r 高频 力 8.1 一个 AF1730 5。7 作为AF2330
C org (2)裸露的土壤 30.4 r H F 年代 W 我 R 1 −18.6 年代 H F V 我 年代 −12.0 r H F V 我 年代 9.3 d 米 一个 x 一个 F 1730年 7.4 一个 年代 一个 F 2330年
(3)实验室 21.3 一个 AF1730 −14.2 d 马克斯 AF1730 10.1 r 高频 SWIR1 −7.6 d 马克斯 AF2330 −7.4 r 高频 力
(1)原位 −14.0 一个 AF900 12.3 d 马克斯 AF700 −10.4 d λ 点燃 AF700 8.8 d 马克斯 AF900 7.9 一个 AF550
铁氧化物 (2)裸露的土壤 22.0 一个 一个 F 550年 −16.7 d 米 一个 x 一个 F 550年 9.5 d 米 一个 x 一个 F 700年 −7.8 一个 一个 F 900年 −7.3 一个 一个 F 700年
(3)实验室 18.9 d 马克斯 AF900 −18.0 一个 AF550 11.7 d λ 点燃 AF550 −8.3 d λ 点燃 AF900 −7.1 一个 AF900
(1)原位 24.4 一个 AF2200 −22.2 d λ 点燃 AF2200 11.2 r 高频 SWIR1 −6.6 年代 高频 力 −6.5 年代 高频 SWIR1
粘土 (2)裸露的土壤 24.8 一个 一个 F 2200年 −21.7 d λ l 我 t 一个 F 2200年 11.8 λ d 米 一个 x 一个 F 2200年 7.7 r H F 年代 W 我 R 1 −5.3 一个 年代 一个 F 2200年
(3)实验室 26.3 一个 AF2200 −21.2 d 马克斯 AF2200 −14.3 d λ 点燃 AF2340 11.6 一个 AF2340 −6.2 λ d 马克斯 AF2340

符号:吸收功能(AFs): 一个 房颤 :区域, d 马克斯 房颤 :最大深度, λ d 马克斯 房颤 :波长的 d 马克斯 房颤 , d 点燃 房颤 :在文学,在波长位置给定深度房颤:不对称因素。船体特征(HFs): r 高频 :在时间间隔平均反射率, 年代 高频 :在间隔斜率。

校准和验证精度B-partial最小二乘回归的方法。模型进一步应用大规模预测土壤成分基于HyMap意象突出显示。

光谱数据集 预处理 不。的因素 校准(941/ 1232样品) 验证(311/ 402样品)
R 卡尔 2 RMS卡尔 RPD卡尔 R 瓦尔 2 RMS瓦尔 RPD瓦尔
自闭症光谱分辨率 (1)原位 一阶导 6 0.79 0.45 2.16 0.53 0.52 1.30
C org (2)裸露的土壤 MSC +一阶导 6 0.82 0.41 2.38 0.65 0.42 1.61
(3)实验室 意味着中心+一阶导 6 0.82 0.40 2.36 0.69 0.45 1.53
(1)原位 意思是中心 7 0.66 0.76 1.74 0.39 0.84 1.14
铁氧化物 (2)裸露的土壤 意思是中心 8 0.69 0.73 1.80 0.43 0.80 1.20
(3)实验室 日志(1 / R)+意味着中心 9 0.82 0.64 2.39 0.43 0.98 1.14
(1)原位 意思是中心 7 0.23 4.34 1.15 0.02 4.63 0.97
粘土 (2)裸露的土壤 日志(1 / R) 5 0.11 4.67 1.06 0.06 4.42 1.02
(3)实验室 意思是中心 8 0.33 4.31 1.22 0.08 4.43 1.02
HyMap光谱分辨率 (1)原位 日志(1 / R)+意味着中心 5 0.77 0.46 2.11 0.59 0.47 1.43
C org (2)裸露的土壤 日志(1 / R)+一阶导 5 0.82 0.41 2.38 0.62 0.45 1.51
(3)实验室 日志(1 / R) 9 0.79 0.43 2.20 0.79 0.32 2.11
(1)原位 意思是中心 7 0.67 0.75 1.75 0.43 0.76 1.26
铁氧化物 (2)裸露的土壤 意思是中心 7 0.66 0.76 1.73 0.45 0.77 1.25
(3)实验室 日志(1 / R)+意味着中心 9 0.81 0.67 2.28 0.49 0.88 1.27
(1)原位 意思是中心 4 0.14 4.58 1.09 0.01 5.06 0.89
粘土 (2)裸露的土壤 意思是中心 5 0.15 4.56 1.09 0.03 4.47 1.01
(3)实验室 一阶导 8 0.37 4.16 1.27 0.10 4.43 1.02

1训练集和测试集的样本数量字段数据集。2训练集和测试集的样本进行实验室数据集。

 土壤有机碳

土壤有机碳的所有预测模型提供优秀的标定精度 R 卡尔 2 在0.75和0.81之间,RMS卡尔在0.45%左右,和所有RPD卡尔高于2.0。良好的相关性也代表散点图的测量与计算土壤有机碳浓度作为模范地显示为模型建立在裸露的土壤光谱HyMap的光谱分辨率(图 8)。测试集验证的准确性好尤其是对实验室数据和校准(RPD在相同的范围瓦尔接近2.0)。这些模型特别显示良好的均方根瓦尔0.35%左右。然而验证精度略低,但合理的字段数据集( R 瓦尔 2 在0.49和0.62之间,RMS瓦尔在0.54%和0.43之间,RPD瓦尔在1.26和1.57之间)。

模型适用于12至15 16计算光谱变量。模型开发基于领域的光谱,专注于几个变量影响,光谱变量描述两个船体(见表特性是非常重要的 4)。在这些模型中, r 高频 SWIR1 , 年代 高频 力 , r 高频 力 主导因素的回归关系(见图 5光谱变量)的缩写。他们一起占超过60%的影响。两个最后的负面影响,因为它是期待,因为增加土壤有机碳含量减少反射明显,从而影响斜率和平均反射率。尽管预计相同的效果 r 高频SWIR1 ,其系数是正的。在实验室模型开发基于数据集,分布在多个光谱变量影响,特别是描述1730海里的吸收是主要的变量。 一个 AF1730 , r 高频SWIR1 , d maxAF1730 显示为重要因素。

校准和验证精度表明,所有预测模型对土壤有机碳的高质量根据应用的分类( 9]。所有模型达到优良的标准基于RPD预测模型卡尔,而他们在中等收入模型验证性能好。光谱变量的选择占最影响回归关系非常不变,似乎合理的尤其是对字段数据,变量的谱连续显示为最具影响力的人。然而,为实验室数据集的影响明显不同,更连接到特定的吸收约1730海里。

铁氧化物

校准铁氧化物的预测模型 R 卡尔 2 在0.62和0.75之间,RMS卡尔0.80%左右,RPD卡尔在1.64和2.01之间。再次与更好的性能实验室光谱。即使 R 卡尔 2 值为所有铁氧化物预测模型是合理的,RPD低卡尔值已经表明减少可预测性。与非常低的验证结果证实了这一点 R 瓦尔 2 (0.17 - 0.26),增加了RMS瓦尔和低RPD瓦尔。特别是模型建立在实验室土壤光谱提供显著降低验证精度。验证样本的散点图揭示了贫穷的预测一个样本的特征是异常高反射率的活力和光谱变化可能不代表。在基于功能的铁氧化物的总量预测模型不达到土壤有机碳的质量预测模型与验证精度更糟。

铁氧化物的预测模型的主要部分是由18到20的21个光谱变量描述光谱特征属性。一个模型只适用于15个变量。模型开发基于领域的光谱, 一个 AF900 , 一个 AF550 , d 马克斯AF550 , d 马克斯AF700 , d λ 点燃 AF900 具有突出的重要性(见表 4)。最大吸收深度通常积极签署,如预期。值得注意的是,900海里的面积吸收特性( 一个 AF900 )因此负面签署。基于实验室光谱模型的开发, d 马克斯AF900 , d λ 点燃AF900 , 一个 AF550 , 一个 AF900 似乎是最重要的光谱变量。还在模型建立在实验室光谱,负的系数 一个 AF900 也可以观察到的特征吸收900海里。这可能表明,随着铁氧化物的内容更加明显和陡峭的特性,尽管广泛的发生。

一起,铁氧化物的预测模型校准是中等质量。在预测精度验证显示缺陷。选择重要的光谱变量包括著名的铁吸收的属性,特别是一个约900海里。

粘土

基于粘土预测模型的校准会导致较低的相关性 R 卡尔 2 在0.17和0.31之间,RMS卡尔不断对光谱分辨率高于4.0%。结果验证使用独立的测试集是穷人和没有显著相关性(RMS瓦尔RPD多达5.13%卡尔低于1.00)。

土壤粘粒含量预测模型是建立基于11到16个派生的光谱变量。 一个 AF2200 , d 马克斯AF2200 , d λ 点燃AF2200 , r 高频 SWIR1 经常出现在最重要的光谱变量的多元线性回归方程中使用所有六个模型(见表 4)。不过,考虑所有的模型,观察影响光谱变量是高度变异。似乎没有特定变量的选择和领域之间的联系或实验室光谱基础数据的来源。一些变量(例如, 一个 AF2200 , r 高频 SWIR1 )出现因此积极签署,然而,例如, d λ 点燃AF2200 和 d 马克斯AF2200 发生积极的和消极的。

减少校准精度,验证精度差和内没有明显的模式中光谱变量回归方程及其具体影响表明了基于回归方法未能建立土壤粘粒含量的显著预测模型。发展表示高度的统计回归关系适应和最后一个低每个模型的重要性。

 5.2。方法B:偏最小二乘回归

表 5概述预测模型的选择土壤成分使用偏最小二乘回归技术。不同的预处理方法的应用允许一致请预测模型的推导过程,应用4和9之间请因素进行回归分析。在超过40%的模型中,意味着定心数据显示的是最好的预处理技术,而反射转换( R ) 日志 ( 1 / R ) 和一阶导数是适当的方法为每个模型的11%。在剩余的情况下,两种预处理方法表现最好的结合。请建立模型的回归系数提供信息的重要性一定波长的回归模型。尽管他们是高度变异根据每个模型的设置,他们的一些信息可以与土壤物理性质有关。

土壤有机碳的预测模型三个光谱数据集非常一致的预处理技术,只显示轻微的差异。校准精度最好的请为每个数据集的模型 R 卡尔 2 介于0.77和0.82之间,RPD卡尔介于2.11和2.38之间为应用光谱分辨率较低的RMS卡尔略高于0.4%。有机碳模型的测试集验证结果还质量好( R 瓦尔 2 在0.53和0.79之间,RMS瓦尔在0.52%和0.32之间,RPD瓦尔在1.30和2.11之间),最高精度为insitu公司现场光谱实验室光谱和低很多。尽管如此,有机碳由偏最小二乘方法建立预测模型表明良好的性能在校准和验证,不显示模型过度拟合。大多数土壤有机碳PLS模型的回归系数表示高度重视对波长的范围内,这可能与土壤颜色(见图 7(一)一个例子)。另外不同的特性发生在短波红外成像可以与吸收土壤有机碳官能团的(见表 1)。

请模型的回归系数的预测土壤有机碳(a)和铁氧化物(b)建立在裸露的土壤野外光谱HyMap 116波段的光谱分辨率。形状和范围的谱相关回归系数选择的数据预处理方法(答: 日志 ( 1 / R ) +一阶导数,B:意味着中心)。

相关图和直方图校正残留的预测土壤有机碳和铁氧化物使用两个提出的建模方法,开发了基于裸露的土壤野外光谱HyMap 116波段的光谱分辨率。方法一:多元线性回归谱参数,方法B:偏最小二乘回归。

铁氧化物预测模型的校准,因此请超过7个因素需要设置重要的模型。请模型显示相关性( R 卡尔 2 )在0.66和0.82之间的所有六个组合光谱数据和光谱分辨率。产生的均方根卡尔在0.64和0.76之间,RPD错误卡尔在1.73和2.39之间表明预测能力略低于有机碳。建立了铁氧化物的测试集验证预测模型显示精度降低,这是最明显的在减少 R 2 和RPD, RMS的大多数模型从校准验证(RMS仅略有增加瓦尔0.98%和0.76之间)。没有达到RPD铁氧化物模型瓦尔1.4指示中预测模型( 9]。请回归系数铁氧化物的模型(示例图 7 (b))在VIS显示增加的重要性波长范围,可与铁离子的电子跃迁乐队主要是一个大约900海里(表 1)。短波红外成像,2208纳米粘土吸收出现,代表一种间接关联与粘土有关的特性。事实上,铁氧化物和粘土含量的相关性可以确定基于化学参考( R 2 0.20)。

粘土的校准精度预测模型 R 卡尔 2 低于0.4和RPD卡尔分别低于1.3,低于中等预测模型的精度。均方根误差不断在4.0%以上。在这些模型中,最好的相关性都达到实验室数据集。相应的验证结果也同样糟糕。回归系数是高度变异。这些估计的结果表明,土壤中粘土含量没有显著预测模型提供足够的精度可以建立在使用偏最小二乘回归技术和基于地面南非研究网站的参考数据。

 5.3。比较建模方法对土壤有机碳和铁氧化物

所有功能关系建立的两种方法,光谱测量在实验室控制导致最好的模型(RMS最低瓦尔和最高RPD瓦尔),因为他们没有偏见的影响出现在现场环境(比较大的方差数据集在图 4)。这是预期,也观察到在先前的研究(例如, 24])。一般来说,insitu公司现场光谱调查通常执行最低。这很可能是由于小石头自然覆盖裸露的土壤表面和令人不安的光谱特征和化学浓度之间的相关性。的自闭症光谱分辨率重采样的现场和实验室数据集HyMap降低分辨率很低,对模型性能的影响。类似的损失模型精度减少光谱分辨率的各种报告的作者时(例如, 8, 24])。

有机碳预测模型,基于功能的高钙之间的性能差异(方法),请技术(方法B)很低。然而基于功能特性的方法的可预测性是略低于请技术。平均RPD卡尔基于功能的高模型的土壤有机碳是2.15和2.27的请模型,同时验证精度是两种方法(平均均方根可比瓦尔0.44%,RPD瓦尔1.58两种方法)。方法最重要的光谱变量的回归模型选择非常一致,意义重大。

预测模型对铁氧化物提供方法的平均校正相关系数约0.67和0.72与RPD方法B卡尔分别为1.77和1.95。这两种方法显示RMS卡尔0.75%左右。RPD低卡尔和高铁校准的均方根误差,这两种方法通常表明减少预测能力在验证结果证实使用独立的测试集( R 瓦尔 2 所有模型低于0.50,RPD瓦尔0.91的平均方法的模型和方法B 1.21模型)。也因此,在铁氧化物的预测模型,采用PLS回归建立的技术表现略好。与这些结果铁氧化物的预测模型是中等质量的两种方法,虽然他们不能达到有机碳模型的预测精度。

详细的结果表明,模型精度对土壤有机碳和铁氧化物的预测,由结合光谱特性分析和线性多元回归技术,在同一个范围统计请提供方法。不过,方法的性能略有降低对土壤有机碳和铁氧化物相比,请的方法。土壤有机碳模型的不同之处在于因此小但更明显的铁氧化物的预测。

图 8散点图显示测量和计算分布的预测土壤成分使用两种建模方法和基于光场光谱分辨率的HyMap图像。这些模型进一步应用到大规模的HyMap图像预测土壤参数。他们的数据在表中突出显示 3和 5。残留的直方图是正态分布的。这清楚地表明,建立的模型能够模型的化学差异选择参数,即使化学内容的参考样品不是正态分布(比较图 3)。

在化学参考数据,2样品对土壤有机碳和1样品铁氧化物引起注意,因为他们超过其他样本的总体布局(比较图 3)。这些样品出现在图的校准散点图 8对这两种方法接近1:1线,这表明它们可能有一个很大的影响在模型校准。但是当测试回归模型建立在光谱库中不包括这些样本,类似的模型和结果实现了。

 5.4。预测土壤的粘粒含量

回归模型建立和评价这项研究不充分预测土壤粘粒含量。这两种方法平均校准相关性很低( R 卡尔 2 0.23所有模型的平均),高RMS卡尔RPD(平均4.41%)和低卡尔(平均为1.15)。对这两种方法验证失败。

解释观测到的相关性差光谱和化学信息粘土含量可能导致地质变量的互层砂岩和泥岩。虽然在研究区开发在相同气候条件下,土壤的形成高度依赖土壤源材料,因此开发基于砂岩很可能有不同的化学和光谱特性的开发基于泥岩。在野外环境中存在的土壤物理外壳可能是一个额外的因素。在研究区土壤物理外壳在本地可以很好开发和达到超过1厘米的厚度。它们形成在雨滴的影响发现土壤表面,导致土壤总量和粒子运动的解体导致薄clay-rich表层低渗透的 49]。对光谱分析土壤外壳是有问题的,因为表层的化学测量的野外光谱仪并不对应于化学上1厘米的大块样品在实验室分析。

 6。摘要和结论

本研究提出了两种方法从光谱相关的土壤参数的量化数据。回归模型建立了土壤有机碳、铁氧化物和粘土含量基于163个样本测量在实验室环境和另外125个样本测量在两个字段设置。原始ASD的光谱数据集进行了光谱分辨率和降低分辨率匹配校准研究高光谱图像区域从空中HyMap传感器获得。第一种方法是一个物理模型基于光谱特征分析。的组合选择的光谱特性是由许多变量用于描述多元线性回归分析。第二种方法,传统的偏最小二乘回归分析被选作比较。最好的请每个光谱库模型被选中作为一个适当的预处理方法的组合和最优数量的潜在变量确定分析交叉验证。

结果表明,给出的两个方法提供类似的功能设置重要的预测模型尤其是对土壤有机碳和铁氧化物。好,中等收入预测模型可以建立(RPD卡尔的2.19 C org 1.83,铁氧化物平均模型的两种方法)。其中,有机碳模型两种方法显示最好的校准精度,RPD卡尔2.15基于功能的回归为2.27,请与相应的RMS回归卡尔(平均的0.44%左右 C org 每种方法的模型)。铁氧化物的模型提出了好的和介质标定精度尤其是显示减少验证性能。这两种方法未能建立重要的粘土含量预测模型(RPD卡尔的1.15和RMS卡尔4.41的平均模型),这可能是一个研究区域的地质变量的结果。

主要是模型建立的土壤有机碳和铁氧化物,基于功能特性的回归方法的预测性能是在相同的范围提供基于PLS回归的统计数据。同样也对这两种方法的模型,观察趋势如减少验证为铁氧化物相比精度校准模型和一般低可预测性的粘土模型。然而,基于功能特性的方法总体上略低于请执行的方法。土壤有机碳的预测差别很小(可比验证性能高和PLS模型),尽管对铁氧化物模型更明显。基于特征模型的光谱变量的回归关系非常一致,支持建立模型的意义,尽管不同的变量主要取决于实验室或现场数据为基础。特定波长的重要性请模型高度变异的结果统计建模过程。相关物理特性只出现在一些模型。

相比其他研究工作在农业环境( 15, 19预测模型的准确性,对这两种方法在本研究开发的略低。这可能是由于大尺寸和南非的高度变异特征研究区由于非农环境,地质和土壤类型的差异。大面积的困难来实现预测模型变化的条件(称为全球校准)在先前的研究(例如, 19])。根据各地- 160数据史蒂文斯et al。 19使用等请技术]预测土壤有机碳的农艺卢森堡的地区。一个RPD瓦尔1.47获得全球校准,在相同的范围作为土壤有机碳的预测模型建立与本研究中给出的两种方法。史蒂文斯等人很可能改善这种RPD瓦尔2.76使用本地校准基于agrogeological地区和土壤类型。在我们的例子中缺乏详细的空间连续信息匹配的变化在南非研究区域地质和土壤类型预防调查当地校准的性能。

进一步factor-lowering预测RPD有机碳的小变化测量内容在我们的研究区域。地面参考站点主要是选择有重大的裸露的土壤成分,这样也可以使用这些网站作为机载高光谱图像验证的目标。因此,没有在密集植被土壤样本地区的特点是输入的增加 C org 因此更高的浓度,增加了可变性在地面参考建模RPD直接相关。这符合报告的结果,例如,( 16, 19]。

这项研究的结果表明,有可能建立回归关系的物质基础上达到PLS-derived回归模型的可预测性。这可以通过一组光谱特性的应用为每个调查土壤成分和这些光谱特性的包含各种属性。统计适应在回归分析是减少到最低。因为一些每个土壤成分的光谱特性进一步结合在该方法建立一个重要的模型中,这种方法是建议类似的物理方法相比,更健壮的分析的基础上只有一个诊断光谱特性或乐队比率。此外,建立模型是相当简单和计算不成问题的有限数量的变量的回归关系。

基于物理的方法通常被认为是更健壮的专门请等统计方法相比,这种方法的可转让性,预计要高(例如,[ 17])。为新开发的方法,这将是测试基于高光谱图像的南非研究区域。为此,开发基于回归关系将被应用于大规模量化关键土壤参数高光谱图像的南非研究区域。要处理的图像数据获得了2009年的联合方法,允许提取土壤光谱信号从混合签名由小规模的土地覆盖的变化引起的。第一次大规模预测结果的表层土壤有机碳和铁氧化物使用校准在裸露的土壤光谱和开发使用了基于回归的方法非常有前途。然而,由于土壤粘粒含量的预测模型没有达到显著的精度,他们不适合申请的推导大规模的土壤信息。考试HyMap图像和在这项研究中提出的方法转移到图像数据超出本研究的范围。

导致土壤参数地图是有价值的信息量化南非研究区域内的土壤退化状态。生态区盘根错节,这些信息可以用来检测侵蚀和退化地区初级阶段关于直接选择性区域的恢复。方法一般是半干旱地区开发和不适应在研究区具体情况。转移到其他地区的类似环境条件将进一步调查。

 确认

这研究是作为博士项目资助的亥姆霍兹的一部分 EOS网络,协作的德国 亥姆霍兹研究中心。的 存在网络(一个 生活的土地倡议)支持强化现场取样。DLR的同事和德国承认他们的宝贵贡献的许多地方工作。

 拉尔 R。 土壤碳封存以减缓气候变化 Geoderma 2004年 123年 1 - 2 1 22 2 - s2.0 - 0842297310 10.1016 / j.geoderma.2004.01.032 Rossel r . a . Viscarra Walvoort d·J·J。 McBratney 答:B。 Janik l . J。 Skjemstad j . O。 可见、近红外、中红外或结合漫反射光谱同时评估各种土壤属性 Geoderma 2006年 131年 1 - 2 59 75年 2 - s2.0 - 31944445266 10.1016 / j.geoderma.2005.03.007 公鸡 T。 Jenssen R。 斯图尔特 一个。 威尔逊 我。 盾牌 T。 HyMap机载高光谱传感器:系统,校准和性能 学报第一EARSEL研讨会成像光谱 1998年 瑞士苏黎世 Bartholomeus h . M。 Schaepman m E。 Kooistra l 史蒂文斯 一个。 Hoogmoed w·B。 Spaargaren o . s . P。 基于光谱反射率的土壤有机碳的量化指标 Geoderma 2008年 145年 1 - 2 28 36 2 - s2.0 - 43649091529 10.1016 / j.geoderma.2008.01.010 Ben-Dor E。 Inbar Y。 陈 Y。 有机物的反射光谱中的可见光近红外、短波红外区期间(400 - 2500 nm)分解过程控制 环境遥感 1997年 61年 1 1 15 2 - s2.0 - 0031177876 10.1016 / s0034 - 4257 (96) 00120 - 4 山 J。 shutt B。 映射复杂的模式的侵蚀和干燥的地中海生态系统稳定 环境遥感 2000年 74年 3 557年 569年 2 - s2.0 - 0034476912 10.1016 / s0034 - 4257 (00) 00146 - 2 Rossel r . a . Viscarra behren T。 使用数据挖掘模型和解释土壤漫反射光谱 Geoderma 2010年 158年 1 - 2 46 54 2 - s2.0 - 77953621099 10.1016 / j.geoderma.2009.12.025 Ben-Dor E。 Banin 一个。 可见光和近红外(0.4 - -1.1 μ米)分析干旱和半干旱的土壤 环境遥感 1994年 48 3 261年 274年 2 - s2.0 - 0028667379 常 C.-W。 Laird d . A。 Mausbach m·J。 Hurburgh c·R。 近红外反射率spectroscopy-principal组件土壤属性的回归分析 美国土壤科学学会杂志》上 2001年 65年 2 480年 490年 2 - s2.0 - 0035042039 Palacios-Orueta 一个。 贾斯廷 s . L。 遥感的土壤属性圣塔莫尼卡山脉。光谱分析 环境遥感 1998年 65年 2 170年 183年 2 - s2.0 - 0032131988 10.1016 / s0034 - 4257 (98) 00024 - 8 Ben-Dor E。 Banin 一个。 近红外分析作为一个快速的方法来同时评估几个土属性 美国土壤科学学会杂志》上 1995年 59 2 364年 372年 2 - s2.0 - 0029500505 极冰原 t·H。 摩根 c·l·S。 布朗 d . J。 标志 c . T。 原位表征土壤粘粒含量的近红外漫反射光谱可见 美国土壤科学学会杂志》上 2007年 71年 2 389年 396年 2 - s2.0 - 33947411704 10.2136 / sssaj2006.0211 Ben-Dor E。 Chabrillat 年代。 Dematte j . a . M。 泰勒 g·R。 山 J。 怀廷 m . L。 大梁 年代。 利用成像光谱研究土壤特性 环境遥感 2009年 113年 1 S38 S55 2 - s2.0 - 67650245616 10.1016 / j.rse.2008.09.019 戈麦斯 C。 Viscarra Rossel r。 McBratney 答:B。 土壤有机碳由高光谱遥感和实地可见光谱预测:一个澳大利亚的案例研究 Geoderma 2008年 146年 3 - 4 403年 411年 2 - s2.0 - 49549111126 10.1016 / j.geoderma.2008.06.011 史蒂文斯 一个。 范Wesemael B。 Bartholomeus H。 Rosillon D。 Tychon B。 Ben-Dor E。 实验室、字段和机载光谱监测农业土壤有机碳含量 Geoderma 2008年 144年 1 - 2 395年 404年 2 - s2.0 - 39549105811 10.1016 / j.geoderma.2007.12.009 Udelhoven T。 艾默灵 C。 jarm T。 定量分析土壤化学性质与漫反射光谱法和偏最小二乘回归分析:可行性研究 植物和土壤 2003年 251年 2 319年 329年 2 - s2.0 - 0038027191 10.1023 /:1023008322682 穆德 诉L。 de熊先生 年代。 Schaepman m E。 娃 t·R。 使用遥感在土壤和地形映射审查 Geoderma 2011年 162年 1 - 2 1 19 2 - s2.0 - 79952959031 10.1016 / j.geoderma.2010.12.018 Selige T。 bohn J。 Schmidhalter U。 使用高光谱图像和高分辨率的表层土映射字段数据多元回归建模过程 Geoderma 2006年 136年 1 - 2 235年 244年 2 - s2.0 - 33751040555 10.1016 / j.geoderma.2006.03.050 史蒂文斯 一个。 Udelhoven T。 丹尼斯 一个。 Tychon B。 Lioy R。 霍夫曼 l 范Wesemael B。 测量在农田土壤有机碳在区域范围内使用机载成像光谱 Geoderma 2010年 158年 1 - 2 32 45 2 - s2.0 - 77953609554 10.1016 / j.geoderma.2009.11.032 里希特 N。 jarm T。 Chabrillat 年代。 Oyonarte C。 Hostert P。 考夫曼 H。 免费的地中海土壤中氧化铁的决心使用漫反射光谱 美国土壤科学学会杂志》上 2009年 73年 1 72年 81年 2 - s2.0 - 60549095930 10.2136 / sssaj2008.0025 Ben-Dor E。 莱文 N。 歌手 一个。 Karnieli 一个。 布劳恩 O。 汲沦谷 g . J。 在沙丘土壤rubification过程的定量映射使用机载高光谱传感器 Geoderma 2006年 131年 1 - 2 1 21 2 - s2.0 - 31944442340 10.1016 / j.geoderma.2005.02.011 Chabrillat 年代。 Goetz a . f . H。 Krosley l 奥尔森 h·W。 使用高光谱图像的识别和映射的膨胀性粘土土壤和空间分辨率的作用 环境遥感 2002年 82年 2 - 3 431年 445年 2 - s2.0 - 0036788020 10.1016 / s0034 - 4257 (02) 00060 - 3 戈麦斯 C。 Lagacherie P。 Coulouma G。 连续去除粘土和碳酸钙含量与PLSR方法估计从实验室和机载高光谱测量 Geoderma 2008年 148年 2 141年 148年 2 - s2.0 - 56249104807 10.1016 / j.geoderma.2008.09.016 Lagacherie P。 Baret F。 Feret j·B。 马德拉清爽的 J。 Robbez-Masson j . M。 使用实验室估计土壤的粘土和碳酸钙,领域和机载高光谱测量 环境遥感 2008年 112年 3 825年 835年 2 - s2.0 - 39749117390 10.1016 / j.rse.2007.06.014 Lechmere-Oertel r·G。 Kerley g . i . H。 整流罩 r·M。 在半干旱的多汁的灌木丛模式和转换的影响,南非 《干旱的环境 2005年 62年 3 459年 474年 2 - s2.0 - 19344363091 10.1016 / j.jaridenv.2004.11.016 劳埃德 j·W。 van den Berg e . C。 帕尔默 a。R。 模式转换和退化的丛林生物群系,南非 报告 2002年 39 南非伊丽莎白港 伊丽莎白港大学陆地生态研究中心 米尔斯 a·J。 异常兴奋的 m V。 转换的灌木丛在东开普省稀树大草原降低土壤质量,南非 植物和土壤 2004年 265年 1 - 2 153年 163年 2 - s2.0 - 14744271783 10.1007 / s11104 - 005 - 0534 - 2 ARC-ISCW 天气数据的三个选择在东开普省气候站 2011年 比勒陀利亚,南非 ARC-ISCW 米尔斯 a·J。 整流罩 r·M。 固碳速率在东开普省的灌木丛恢复两个地方,南非 恢复生态学 2006年 14 1 38 49 2 - s2.0 - 33646709385 10.1111 / j.1526 - 100 x.2006.00103.x Walkley 一个。 黑色的 我一个。 考试的Degtjareff方法确定土壤有机质和铬酸滴定法的修改建议 土壤科学 1934年 37 29日 37 10.1097 / 00010694-193401000-00003 用 O。 杰克逊 M。 氧化铁从土壤和粘土去除dithionite-citrate系统与碳酸氢钠缓冲 粘土和粘土矿物 1958年 7 317年 327年 10.1346 / CCMN.1958.0070122 哇 g·W。 波特 j·W。 悬疑类 一个。 粒度分析 土壤分析的方法。第1部分 1986年 2日 美国威斯康星州麦迪逊 美国社会的农学、土壤科学的社会 383年 411年 Savitzky 一个。 戈利 m·j·E。 平滑和分化的数据简化的最小二乘法程序 分析化学 1964年 36 8 1627年 1639年 2 - s2.0 - 0042553279 克拉克 R。 岩石和矿物的光谱,光谱的原则 遥感地球科学 1999年 3 纽约,纽约,美国 约翰威利& Sons 3 58 包姆加德纳 m F。 席尔瓦 l F。 贝赫尔 L . L。 斯通内尔 e·R。 土壤的反射特性 农学的进步。38卷。 1985年 1 44 2 - s2.0 - 0022240364 亨特 G。 索尔斯堡 J。 矿物和岩石的可见光和近红外光谱:1。硅酸盐矿物 现代地质学 1970年 1 283年 300年 格罗夫 C。 钩 年代。 Paylor二世 E。 160年实验室反射光谱矿物质,0.4到2.5微米 1992年 加利福尼亚州帕萨迪纳市,美国 喷气推进实验室,国家航空和宇宙航行局,喷气推进实验室发布92 - 2 亨特 G。 索尔兹伯里 J。 Lenhoff C。 矿物和岩石的可见光和近红外光谱:III。氧化物和氢氧化物 现代地质学 1971年 2 195年 205年 莫里斯 r . V。 劳尔 h . V。 劳森 c。 吉布森 e·K。 小 Nace g。 斯图尔特 C。 光谱和其他物理化学性质的亚微米赤铁矿粉( α菲2O3)、磁赤铁矿( γ菲2O3)、磁铁矿(Fe3O4)、针铁矿( α-FeOOH)和纤铁矿( γ-FeOOH) 地球物理研究杂志》 1985年 90年 4 3126年 3144年 2 - s2.0 - 0022268379 亨特 g·R。 卡迈克尔 r S。 岩石和矿物的光谱属性 岩石物理性质的手册 1982年 1 美国佛罗里达州波卡拉顿 CRC的新闻 295年 385年 CSIR 气象数据的气候站在灌木丛生态区,南非 2009年 比勒陀利亚,南非 CSIR 斯通内尔 e·R。 包姆加德纳 m F。 反射的表面土壤的变化特征 美国土壤科学学会杂志》上 1981年 45 6 1161年 1165年 2 - s2.0 - 0019728521 克拉克 r . N。 Roush t . L。 反射光谱:定量分析遥感应用的技术 地球物理研究杂志》 1984年 89年 7 6329年 6340年 2 - s2.0 - 0021644332 克拉克 r . N。 斯韦兹 g。 Livo k . E。 Kokaly r F。 萨特利 美国J。 道尔顿 j·B。 麦克道戈尔 R R。 绅士 c。 成像光谱:地球和行星遥感与美国地质调查局Tetracorder和专家系统 地球物理研究杂志》上 2003年 108年 12 5 1 2 - s2.0 - 1542318143 范德梅尔先生 F。 高光谱图像光谱吸收特性的分析 国际期刊《应用地球观测和地质信息 2004年 5 1 55 68年 2 - s2.0 - 1342331476 10.1016 / j.jag.2003.09.001 凯斯勒 W。 多元Datenanalyse 2007年 德国魏因海姆 Wiley-VCH 1 - 荒原 年代。 Sjostrom M。 埃里克森 l pls回归:化学计量学的一个基本工具 化学计量学和智能实验室系统 2001年 58 2 109年 130年 2 - s2.0 - 0035965476 10.1016 / s0169 - 7439 (01) 00155 - 1 Viscarra Rossel r·a·V。 说:软件进行最优化分析的光谱数据 化学计量学和智能实验室系统 2008年 90年 1 72年 83年 2 - s2.0 - 36749102560 10.1016 / j.chemolab.2007.06.006 麦金太尔 d S。 土壤渗透性测量外壳由雨滴形成的影响 土壤科学 1958年 85年 185年 189年 克拉克 r . N。 王 t . V V。 Klejwa M。 斯韦兹 g。 枝条 N。 高光谱分辨率反射光谱的矿物质 地球物理研究杂志》 1990年 95年 8 12 680年 2 - s2.0 - 0042894381 亨特 g·R。 颗粒矿物的光谱特征中的可见光和近红外 地球物理学 1977年 42 3 501年 513年 2 - s2.0 - 0017480895