高光谱土壤元素估测

1、高光谱遥感的在估测土壤有机质含量上的应用摘要：高光谱遥感技术以其光谱分辨率高、波段连续性强、数据丰富的特 点，因而在土壤养分研究中得到广泛应用。通过土壤有机质的高光谱遥感 分析，可以充分了解土壤养分的状况及动态变化，为指导农业生产及保护 农业生态环境提供科学依据。本文分析了土壤有机质的光谱特征，介绍了 一种基于高光谱的土壤有机质估测方法。关键词：高光谱遥感；土壤有机质；光谱分析；含量估测Hyperspectral remote sensing in the estimation ofsoil organic matter contentAbstract： Hyperspectral remot

2、e sensing technology with its high spectral resolution band continuity, rich data, which has been widely applied in the soil nutrient research .With the hyperspecttal remote sensing analysis of soil organic metter, we can fully understand the condition and dynamic changes of soil nutrients, and also

3、 can provide a scientific basis for the guidanee of agricuItural produetion and the protection of agricuItural ecological environment. This article analyzes the spectral characteristics of the soil organic matter, and an estimation method based on hyperspectral.Key words: remote sensing soil organic

4、 matter; spectral analysis content prediction1高光谱简介im高光谱遥感是高光谱分辨率遥感(Hyperspectral Remote Sensing) 的简称，是上世纪80年代发展以来的遥感科学的前沿领域,她的出现是遥 感界的-场革命。它是指利用很多窄的电磁波波段获取物体有关数据的技 术，它可在电磁波的紫外、可见光、近红外、中红外以至热红外区域，获 取许多非常窄且光谱连续的图像数据。这样，在传统的二维遥感的基础上 增加了光谱维，形成了一种独特的三维遥感。鉴于高光谱遥感的诸多优点， 人们逐渐尝试在植被、生态、大气、环境、考古、农业等方面进行应用， 取得

5、了一系列的进展和很多成就，使高光谱遥感成为科研必不可少的工具 和技术。高光谱遥感的基础是测谱学，主要在电磁波谱的紫外、可见光、近红外 和屮红外区域，获取许多非常窄口光谱连续的图像数据的技术。成像光谱仪 为每个像元提供数十个至数百个窄波段(通常波段宽度10nm)的光谱信息, 能产生一条完整而连续的光谱曲线。国际遥感界公认:光谱分辨率在1(T入 次数量级范围的定为多光谱(Multlspectral),这样的遥感器在可见光和近 红外光谱区域通常只有几个波段,如美国LandsatMSS. TM,法国的SPOT等; 而光谱分辨率在10_2A认的遥感信息称之为高光谱(Hyperspectral);随着 遥

6、感光谱分辨率的进一步提高,当达到10)以上时，遥感即进入超高光谱 (Ul-traspectral)阶段，许多物质的特征往往表现在一些狭窄的光谱范围 内，如何捕获物体的这种特征而又不失其整体形态及与周围地物的关系，这 就是高光谱分辨率遥感产生和发展的基础和要回答的问题。高光谱遥感的 特点是光谱分辨率高，波段连续性强，能获得多光谱传感器无法获得的精细 的光谱信息，由于其光谱分辨率高达纳米数量级，因此传感器在0.4- 2. 5p m范围内可细分成几十个，甚至几百个波段,光谱分辨率为5-10nmo2高光谱遥感的土壤有机质含量估测土壤是万物生存的物质基础，是人类财富的主要来源，是我们祖先最早 利用的生产

7、资料。土壤是地壳表层岩石矿物的风化产物(母质)，在气候、生 物、地形等环境条件和时间综合因索的作用下形成的一种特殊自然体。土 壤是一种复杂的混合体,它的物理与化学性质是由其组分所决定，这些有着 不同理化特征的丄壤成分及其质地和结构决定着各种上壤类型及其独特的 光谱特性。土壤有机质是存在于土壤中所有含碳的有机物质，包括土壤中各种动 植物残体、微生物体及其分解和合成的各种有机物质，是土壤中各种营养 元素的重要来源(如:氮、磷等)。它即含有刺激植物生长的胡敏酸等物质， 又具有胶体特性能吸附较多的阳离了，使土壤具有保肥力和缓冲性。同时 它还能使土壤疏松并形成结构，从而改善土壤理化特性。它是微生物必不

8、可少的碳源和能源。因此，土壤有机质含量的多少，是土壤肥力的一个重 要指标。通过高光谱估测上壤有机质含量，建立有机质高光谱预测模型, 旨在快速评价有机质含量，可以为提高上壤有机质含量水平，改善土壤条 件服务；同时对提高作物产量有着重要意义。2.1 土壤有机质光谱特征分析早在1965年Bowers等指出，影响土壤反射率的因素主要有土壤 质地，表面粗糙度、氧化铁含量、土壤水分及有机质含量，其中有机质是 影响土壤色调进而影响光谱反射特性的重要因索.Gunasulis等研究表 明，土壤有机质含量与红光波段反射率有较强的相关性.Chang等从有 机质结构组成的角度解释了近红外区的光谱特征，认为有机质增加一

9、般会 降低谱线的反射率值.徐彬彬等7J通过比较去除有机质前后的土壤光谱反 射率曲线，发现有机质的影响主要在可见光和近红外波段.贺军亮等利 用统计单相关分析方法确定511nm为诊断土壤有机质含量的敏感波段，并 对507516nm敏感范围反射率做均值标准化处理后获得有机质诊断指 数.刘焕军等通过研究东北地区黑土土壤有机质反射光谱特征，发现黑 土土壤有机质敏感波段主要在350llOOnm,其屮680810nm最为 显著，并利用多元逐步回归分析建立了丄壤有机质预测模型.张娟娟等“ 构建了近红外区反射率差值指数预测上壤有机质含量沈润平等应用多 元线性逐步回归法和人工神经网络建立了有机质反演模型，结果表明

10、人工 神经网络模型普遍优于回归模型。综上所述，土壤有机质敏感波段主要位于可见光及近红外区，多变量 回归分析方法常用于土壤有机质含量分析，包括线性和非线性分析方法. 其中非线性分析方法，例如人工神经网络，可能得到的精度会较高，但模 型复杂且不具备稳定性和可靠性.Vasques等切用高光谱数据预测土壤含 碳量时比较了各种线性及非线性的多变量分析方法，发现采用多元线性逐 步回归法和偏最小二乘法得到的预测效果最好且最稳定.高光谱预测上壤 有机质含量的研究虽然很多，但针对不同地区建立的模型往往差别很人， 主要原因是成上母质对上壤光谱影响很人;同时建模波段及建模方法的选 择也会影响研究结果。2. 2 土壤

11、有机质含量估测方法2. 2.1 土壤样品采集采集020cm深的土层土壤。土壤样本经风干、研磨后过2mm孔筛。 每个样品分2份，一份用于光谱测量；一份用于有机质含量分析。2. 2. 2 土壤有机质分析土壤有机质含量用高温外热重珞酸钾氧化-容量法测定。2. 2. 3光谱测量运用便携式光谱仪ASD FieldSpec Pro测量空样品Q勺反射光谱，波 长范围为3502500nm（350lOOOnm之间光谱采样间隔为1. 4nm,光 谱分辨率为3nm； 10002500nm范围内光谱采样间隔为2nm,光谱分辨 率为10nm；光谱仪最后将数据重采样成lnm）。将每个土样测得的10条反 射光谱曲线经算术平

12、均后则为该土样的实际反射光谱数据。2. 2. 4光谱数据预处理（1）光谱数据的数学变换光谱一阶微分和倒数的对数是两种常用的光谱处理方法間，其中，对 光谱曲线作一阶微分变换，可以对重叠混合光谱进行分解以便识別，扩大 样品之间的光谱特征差异；光谱反射率经倒数的对数变换后，不仅趋向于 增强可见光区的光谱差异，而且趋向于减少因光照条件变化引起的乘性因 素影响剛。（2）特征吸收带的提取剔除水吸收峰影响严重的波段区域。进行相应去包络线操作，获得光 谱波段深度曲线每一样点土壤反射率归一到对应的光谱背景上，有利于光 谱曲线Z间特征波段的比较，并且可以有效地突出光谱曲线的吸收和反射 特征。2. 2. 5反演模型

13、多元线性逐步回归模型多元线性逐步回归分析(stepwise multiple linear regression, SMLR) 的基本思想是:对全部的自变量，按它们对因变量贡献的人小进行比较，并 通过F检验法，选择偏回归平方和显著的变量进入回归方程，每一步只引 入1个变量，同时建立一个偏回归方程。当1个变量被引入后，对原已 引入回归方程的变量，逐个检验它们的偏回归平方和。如果由于引入新的 变量而使得已进入方程的变量变为不显著时，则及时从偏回归方程中剔除。 在引入了 2个自变量以后，便开始考虑是否有需要剔除的变量。只有当回 归方程中的所有自变量对因变量都有显著影响而不需要剔除时，再考虑从 未选入

14、方程的自变量中，挑选有显著影响的新的变量进入方程。不论引入 还是剔除一个变量都称为一步。断重复这一过程，直至无法剔除已引入的 变量，也无法再引入新的自变量时，逐步回归过程结束。偏最小二乘回归模型偏最小二乘回归(partial least squares regression, PLSR)考虑 p 个因变量与m个自变量的建模问题，基本做法是首先在自变量集中提取第 1成分t”同时在因变量集中提取第1成分5,要求仁和5相关程度达 到最大.然后建立因变量与仁的回归，如果回归方程已达到满意的精度， 则算法终止.否则继续对第2对成分的提取，直到达到满意的精度为止。 最终对白变量集提取r个成分t】，tr,偏

15、最小二乘回归将建立因 变量与ti, t2,，t的回归式，然后再表示为因变量与原自变量的回归方 程式，即偏最小二乘回归方程式17。为避免过拟合，检验所建模型的有 效性，PLSR所选择的成分数量h可通过交叉有效性检验确定例如对于 某个应变量，具体做法为:将m个样本中m-1个用作训练样本，的1个样 本作检验样本.第1次先将第1个样本作检验样本，用其余m-1个样本 建模，用检验样本代入模型，可以求得1个估计值，记为刃，”第2次 再将笫2个样本留下作检验样本，用其余样本建模，将2个样本代入模 型求得估计值y2“如果循环m次，将这m个残差值平方求和，用press 表示该值，表达式为：PRESS(h) =(

16、yiy严i=l另外，得到拟合h个成分的回归方程后，汁算所有样本点预测值y i,并 求其误差平方和SS(h) :SS(h)二器丄-“)2 o PRESS(h)/ss(h-l)值越小越 好，一般设定限制值为0. 05,当PRESS(h)/ss(h-l)(1-0. 05)2时，增加新 的成分，对减小方程预测误差无明显改善作用。如此，可确定偏最小二乘 法所提取的因了数。偏最小二乘回归分析在建模的过程中集屮了主成分分 析，典型相关分析和线性回归分析方法的特点，因此在分析结果中，除了 可以提供一个更为合理的回归模型外，还可以同时完成一些类似主成分分 析和典型相关分析的研究内容，提供更丰富、深入的一些信息。

17、3问题与展望土壤光谱是其各种成分的光谱响应的综合反映,土壤是由固、液、气三 相物质共同组成的多相分散体系,光有机质的组成就相当复杂，铁在土壤屮 存在也形式多样，每种成分的光谱响应都是不一样的，如有机质屮的胡敏酸 和富里酸就有不同光谱反射曲线。在以后的研究中，应该进行更细致的研究, 了解丄壤各种成分光谱响应机理，区分出哪些特征是其本质的特征表现。这 样便于光谱分析研究。尽管国内外对土壤有机质含量的光谱效应进行了大量的研究工作，但 在认识上还存在一定的分歧，主要是由于以上研究大多基于多种土壤类型， 不同的土壤，由于受气候、母质、地形、生物等成土因素不同和成土年龄 的差异，以及人类活动的影响，其理化

18、特性具有明显差别。土壤组成物质 的多样性以及每种组成物质独特光谱特性，使得各类上壤光谱具有自己的 特点。因此，将其他研究结果应用到特定上壤时，都可能产生较人误差。 对于丄壤有机质估测，模型的适用性往往较差，建立-种实用而且通用的 模型是现阶段探索的目标之一。要建立准确的土壤反射率特征参量的反射 率反演模型需要进行大量而精确的实验室与田间的光谱测量，若将其推广 到图像上则需进行严格地校正和大量反复的实验。参考文献：黄昌勇.土壤学J.北京，中国农业出版社,2000鲍士旦土壤农化分析,2000,中国农业出版社武彦清，张柏，宋开山，刘焕军等.松嫩平原土壤有机质含量高光谱反演研究J,中 国科学院研究生院

19、学报，2011, 5, 28 (2), 187-19Bowers S A Hanks R J. Reflection of radiant energy from soils J . Soil Science, 1965 100(2) : 130-138.Gunsaulis F R, Kocher M F, Griffis C L. Surface structure effects on close-range reflectance as a function of soil organic matter content J . Soil Chemistry and Physics, 1991 34(2):641-649.Chang C W, Lai rd D A. Near-infrared reflectance spec tgscopic anal y si s of soil C and N J . Soil Science, 2002 167(2):110-116.徐彬彬，戴昌达.南疆土壤光谱反射特性与有机质含量的相关分析J.科学通报， 1980, 25(6) : 282-284.贺军亮，蒋建军，周生路，等.土壤有机质含量的高