一种基于光学模型的活体叶片色素重叠吸收系数分离方法与流程

1.本发明属于遥感监测技术领域，具体涉及一种基于光学模型的活体叶片色素重叠吸收系数分离方法。


背景技术：

2.植物叶片色素信息特征是植物生理生态状况的重要表征之一，在植物叶片内，叶片光合色素和光保护色素的体内吸收特征与叶片400-800nm区域光谱密切相关。
3.在400-800nm区域的叶片光谱包含了多种光合色素的信息，包括叶绿素a(chla)、叶绿素b(chlb)和类胡萝卜素(cars)，以及光保护色素，如花青素(ants)。这些色素与植被的生理生态功能密切相关。其中，chla和chlb分子进行叶片光合电子的收获、运输和吸收，chla分子进行叶片荧光发射，cars进行叶黄素循环的叶片热耗散，ants进行叶片过剩光能的淬灭。利用遥感数据反演叶片色素与体内叶片色素的光学性质密切相关。因此，通过更好地测量和了解植物叶片中存在的色素，可以更好地在体内测定叶片中色素的光学性质，从而更好地了解植物的生长状况。
4.随着遥感技术的发展，高光谱遥感技术为植被色素信息特征检测或监测提供了一种有效的途径，因此，使用高光谱遥感技术在色素信息特征波段400-800nm的高能光谱区间检测或监测叶片色素变化，可间接地提供植被生理生态的信息特征。由于叶片色素在可见光谱中的选择性光吸收，单个叶片的颜色变化遵循其色素的比例，这提供了一个利用叶片反射率/透射率来检索叶片色素浓度的机会。但位于400-800nm的特定波段色素基团的吸收重叠特性掩盖了或影响了单个色素对叶片光谱的贡献，限制了单一色素含量的光谱反演。因此，精确反演各种单一色素含量中，分离叶片色素基团吸收之间的重叠特征是提取叶片色素含量的关键因素。
5.基于光学辐射传输模型在叶片色素的波段特征区域定量描述影响叶片光学属性因子是准确获取植物叶片色素信息特征的有效手段。通过光学辐射传输模型与色素波段重叠特征分离的耦合为活体叶片中提取各种色素吸收系数提供了一个机遇。


技术实现要素：

6.针对现有技术的不足，本发明提出了一种基于光学模型的活体叶片色素重叠吸收系数分离方法，改进叶片光学辐射传输模型prospect，从活体叶片光谱中提取各种单一色素的吸收系数，为利用光学特性监测植被生理生态特性提供更多的信息。
7.一种基于光学模型的活体叶片色素重叠吸收系数分离方法，包括以下步骤：
8.步骤一、数据采集
9.采集各种纯色素在有机溶剂中的吸收峰位和吸收峰个数，作为数据集i。采集多种不同状态下的植物叶片，记录植物叶片的方向半球反射光谱与透射光谱，以及各种色素的含量，作为数据集ii。
10.作为优选，将乙腈、甲醇、二氯甲烷按照体积比60∶20∶20混合后加入少量37％的盐
酸溶液，配置成ph值为5的有机溶剂。
11.作为优选，数据集ii中还包括植物叶片的含水量以及各种色素含量的信息。
12.步骤二、数据预处理
13.根据数据集ii中叶片的反射光谱与透射光谱，计算得到每个叶片的叶片结构参数n。选择数据集ii中的白化叶片，根据这些白化叶片的光谱数据、生化组分与叶片结构参数，计算得到叶片的基线吸收系数k0(λ)。
14.作为优选，选用反射率最大波段、透射率最大波段和吸收率最小波段的光谱特征进行叶片结构参数的计算。
15.步骤三、模型建立与系数分离
16.对传统proapect模型中的叶片吸收系数进行扩展，同时考虑基线吸收系数k0(λ)，将叶片吸收系数k(λ)分解为：
[0017][0018]
其中，λ表示波长；k1(λ)、k2(λ)、k3(λ)和k4(λ)分别表示叶片中chla、chlb、cars和ants的特定吸收系数；c
chla
，c
chlb
，c
cars
和c
ants
分别表示chla，chlb，cars和ants在活体叶片中的浓度。
[0019]
由于色素的吸收特征来源于各峰的吸收特征之和，因此色素的特定吸收系数ki(λ)由各对应色素的所有吸收峰的总和得出：
[0020][0021]
其中j表示第i种色素的吸收峰个数，k
i，j
(λ)表示第i种色素的的第j个吸收峰系数。使用gauss-lorentz函数对每种色素的每个吸收峰k
i，j
(λ)进行统一表征：
[0022][0023]
其中，k
i,j,v
，k
i,j,h
和k
i,j,w
为活体叶片中第i种色素的第j个吸收峰的高斯比、峰高和半峰全宽；a
i,j,p
是第i种色素在有机溶液中的第j个吸收峰的峰值位置，k
i,j,δλ
是叶片中第i种色素的第j个吸收峰的光谱位移。使用数据集i和数据集ii，通过最小二乘法优化迭代获取chla、chlb、cars和ants在活体叶片中吸收系数，完成活体叶片中单一色素吸收系数的分离。
[0024]
本发明具有以下有益效果：
[0025]
通过一系列吸收特征参数，可以清晰地描绘出活体叶片中分离色素之间的重叠特征，并为未来的细化描述描述叶光学性质的框架，同时，首次分离获得叶绿素a和叶绿素b在活体叶片中的吸收系数。
附图说明
[0026]
图1为实施例中测得的纯色素在有机溶剂中的吸收光谱；
[0027]
图2为实施例中测得的纯色素在有机溶剂中400-800nm范围的吸收峰的数目和位置；
[0028]
图3为实施例收集的数据集ii中叶片的色素含量；
[0029]
图4为实施例收集的数据集ii中叶片的反射和透射特性；
[0030]
图5为实施例中各种色素在活体叶片中吸收光谱与有机溶液中吸收光谱的光学性质比较；
[0031]
图6为实施例中基于光学辐射传输模型率定的活体叶片色素吸收系数参数特征；
[0032]
图7为实施例中得到的活体叶片色素特异性吸收特性。
具体实施方式
[0033]
以下结合附图对本发明作进一步的解释说明；
[0034]
一种基于光学模型的活体叶片色素重叠吸收系数分离方法，包括以下步骤：
[0035]
步骤一、数据采集
[0036]
s1.1、数据集i的采集
[0037]
将乙腈、甲醇、二氯甲烷按照体积比60:20:20混合后加入少量37％的盐酸溶液，配置成ph值为5的有机溶剂。使用高效液相色谱系统中的岛津紫外-可见检测器测量纯色素chla、chlb、β-car、vi、an、ze、ne、lu和ants在该有机溶液中的吸收光谱，如图1所示，采用一阶导数和二阶导数联合方法对吸收光谱进行求导，获取对应色素在400-800nm波长范围下特定吸收系数明显的吸收峰个数与峰位，如图2所示。记录各种纯色素在有机溶剂中的吸收光谱、吸收峰位和吸收峰个数，作为数据集i。
[0038]
s1.2、数据集ii的采集
[0039]
采集多种不同状态下的植物叶片，这些叶片来自不同的常绿和落叶的乔木、灌木、亚灌木以及草本阔叶植物，包括其健康和患有白化病的幼叶、成熟叶和衰老叶。选择agilent 1100高效液相色谱仪进行叶片色素含量与含水量的测定。由于叶片中花青素的含量需要获取其总含量，因此使用分光光度计进行测定。测定结果如图3所示，其中cars的浓度表示为叶片中lu、an、ze、vi、ne和β-car浓度的总和。使用uv-3600分光光度计和积分球测定叶片的半球反射光谱与透射光谱，并记录如图4所示，400-800nm波长范围的反射率与透射率。多种不同状态下的植物叶片信息作为数据集ii。
[0040]
步骤二、数据预处理
[0041]
s2.1、计算叶片结构参数n
[0042]
叶片解剖结构的变化能够产生对400-2500nm区间整个光谱特征的影响，其中在800-1350nm区间存在一个反射率最大的波段、一个透射率最大的波段和一个吸收率最小的波段。在这个光谱区域中低的吸收率，也是叶片叶肉组织细胞在叶片光谱特征曲线中响应最敏感的区域。因此，选用这些近红外波段的光谱特征获取模型叶片结构参数，根据数据集ii中叶片的反射光谱与透射光谱，计算得到每个叶片的叶片结构参数n。
[0043][0044]
其中，r
mea
(λi)和t
mea
(λi)分别表示在波段λi叶片实测的反射率和透射率；r
mod
(λi)和t
mod
(λi)分别表示相应叶片在波段λi模拟的反射率和透射率；λ1、λ2、λ3分别是反射率最大的波段、透射率最大的波段和吸收率最小的波段。基于光谱最小拟合法，通过迭代使不同生化组含量的叶片实测光谱与相应的模拟光谱达到最小差异，得到估算参数，然后使用最小二乘法令χ(n)函数达到最小值，获取每个叶片的叶片结构参数n。
[0045]
s2.2、计算基线吸收系数k0(λ)
[0046]
基线吸收系数为非色素光敏物质在叶片体内吸收特性。白化叶片中含有痕量植物
色素，在400-800nm范围内叶片的吸收特性主要受叶片中蛋白质、酚类物质的影响。选择数据集ii中枫树的白化叶片，计算白化叶片的叶片吸收系数k
albino
(λ)：
[0047][0048]
其中，c
i，0
表示枫树白化叶片中第i种色素的浓度，n0表示结构参数；n0表示枫树白化叶片色素种类的数量。根据这些白化叶片的光谱数据、生化组分与叶片结构参数，基于公式(3)所述的最小迭代函数，得到叶片的基线吸收系数k0(λ)：
[0049][0050]
其中r
mea，0
、t
mea，0
和r
mod，0
、t
mod，0
分别表示白化叶片的实测和模拟的反射和透射光谱。
[0051]
步骤三、模型建立与系数分离
[0052]
s3.1、模型建立
[0053]
prospect模型能够通过色素的吸收特性很好地定量光在也叶片中的辐射传输特性。在prospect模型中，叶片吸收系数k(λ)函数是直接表征色素信息。由于ants、chlb、cars和chla的吸收特征存在明显的重叠带，为了同时获取chla、chlb、cars和ants的吸收系数，需要对传统proapect模型中的叶片吸收系数进行扩展，同时考虑基线吸收系数k0(λ)，将叶片吸收系数k(λ)分解为：
[0054][0055]
其中，λ表示波长；k1(λ)、k2(λ)、k3(λ)和k4(λ)分别表示叶片中chla、chlb、cars和ants的特定吸收系数；c
chla
，c
chlb
，c
cars
和c
ants
分别表示chla，chlb，cars和ants在活体叶片中的浓度。
[0056]
由于色素的吸收特征来源于各峰的吸收特征之和，因此色素的特定吸收系数ki(λ)由各对应色素的所有吸收峰的总和得出：
[0057][0058]
其中j表示第i种色素的吸收峰个数，k
i，j
(λ)表示第i种色素的的第j个吸收峰系数。
[0059]
s3.2、系数分离
[0060]
高斯-洛伦茨(gauss-lorentz，g-l)函数能够很好地分离具有波段重叠特征的吸收光谱，因此使用g-l函数化各种色素的特定吸收系数，使k(λ)函数中各种色素特定吸收系数的同性参数特征转化为相应色素特定吸收系数的g-l参数的非同性参数关系，消除由于目标色素之间在叶片吸收特征的重叠在分离过程中的遮挡现象。
[0061]
使用gauss-lorentz函数对每种色素的每个吸收峰k
i，j
(λ)进行统一表征：
[0062][0063]
其中，k
i，j，v
，k
i，j，h
和k
i，j，w
为活体叶片中第i种色素的第j个吸收峰的高斯比、峰高和半峰全宽；a
i，j，p
是第i种色素在有机溶液中的第j个吸收峰的峰值位置，k
i，j，δλ
是叶片中第i种色素的第j个吸收峰的光谱位移。
[0064]
由统一表征后的公式(3)可知，g-l函数使用需要已知每个色素的吸收峰个数和吸收峰位，这些数据来源于步骤一收集的数据集i和数据集ii。由于有机溶液中的极性小于水相的极性，而且，极性高的溶剂中吸收峰个数小于等于极性低的溶剂，因此，本实施例设定活体叶片各种色素吸收系数中吸收峰个数小于等于有机溶液。另外，考虑溶剂溶液与活体叶片极性的差异对吸收峰位置的红移或蓝移，本实施例使用各色素吸收峰在有机溶液中吸收峰位置 红移/蓝移的变量。
[0065]
将公式(4)代入公式(5)中，结合多种光合作用色素和光保护色素的活体叶片内吸收系数，通过最小二乘优化最小化价值函数χ(ki(λ))，确定了chla、chlb、cars和ants在体内叶片中的吸收系数(ki(λ))。价值函数χ(ki(λ))为：
[0066][0067]
式中m为数据集ii中的叶片样本数；r
mea
和t
mea
，r
mod
和t
mod
分别为所选叶片的实测半球反射率和半球透射率，模型反射率和透射率。
[0068]
然后通过输入有机溶液中色素各吸收峰的峰位、叶片结构参数、基线吸收系数及其每个叶片的反射和透射光谱，利用最小二乘法的优化方法，进行活体叶片单一色素各个吸收的参数迭代获取，最终获取活体叶片中chla，chlb，cars和ants的吸收系数。
[0069]
步骤四、实验验证
[0070]
s4.1、与纯色素在有机溶剂中的吸收光谱对比
[0071]
通过公式(7)确定色素在活体叶片中的吸收系数，并与纯色素在有机溶液中的吸收光谱进行了光学性质比较，如图5所示，k
chla
和k
chlb
有两个明显的吸收峰，并有各自的吸收峰位置。k
chla，1
的位置为419nm，k
chla，4
位置为679nm，k
chlb，1
的位置为468nm，k
chlb，3
的位置为661nm。k
cars
和k
ants
有一个明显的吸收峰，其吸收峰位分别位于482nm和544nm。模型测定的色素吸收系数的形状与相应纯色素在有机溶液中的吸收光谱具有相似的起伏形状。
[0072]
本方法得到的色素特定吸收系数中吸收峰位置与色素在有机溶液中吸收峰位置出现相关的不同位移，k
chla，1
位置移动到较短的波长方向，即蓝移；k
chla，4
位置移动到更长的波长，即红移。使用光谱位移参数δλ
i，j
来描述峰位置的变化，结果如图6所示，k
chla，1
和k
chla，4
的位移分别为-13nm和15nm，k
chlb，1
和k
chlb，3
分别为4nm和11nm，k
cars
和k
ants
分别为39nm和14nm。在数据集ii使用的叶片的平均含水量通常超过50％。此外，鲜叶可视为被水介质包围。活体叶片和有机溶液的极性不同，因此，在体内叶片中测定的叶色素的环境极性大于在混合有机溶液。同时，由于活体叶片和有机溶液中发色团的电子跃迁模式的极性差异的影响，可以引起单个色素吸收峰的蓝移和红移所涉及的不同位移。因此，在活体叶片和有机溶液之间的叶片色素的吸收光谱/系数的相应吸收峰位置有蓝移或红移。
[0073]
g-l函数可以准确地表征色素的吸收光谱，因为该函数的参数具有吸收光谱的明确物理意义，其中k
i，j，w
和k
i，j，p
可以描述吸收峰的主要吸收特征。采用由k
i，j，w
和k
i，j，p
驱动的raf参数来描述所测定的色素吸收系数的主要特征。如图6所示，k
chla
的两个明显的吸收峰位于400-434纳米和659-699纳米的区域。k
chlb
在442-495nm和639-683nm处有明显的吸收峰。k
cars
和k
ants
分别在447-517nm和494-594nm区域。
[0074]
s4.2、与活体叶片对比
[0075]
此处选用两组叶片的反射率来评估活体叶片吸收峰位置的光谱位移的有效性。从数据集ii中提取了一组具有一定范围ants浓度的叶子和一组缺乏chlb的叶子。在图7的(a)
表格中，ants的吸收在500-580nm区域占主导地位，通过计算这些反射光谱的一阶和/或二阶导数获得活体叶片内ants的吸收峰位置(545nm)。在图7的(b)中，显示了叶片内中chla和cars的吸收峰位置：a
cars
＝480nm，a
chla，2
＝590nm，a
chla，3
＝628nm，a
chla，4
＝676nm。而本方法测定的吸收峰位置：a
ants
＝544nm，a
cars
＝482nm，a
chla，2
＝591nm，a
chla，3
＝627nm，a
chla，4
＝679nm。这与活体叶片内吸收峰相对应，证明本方法成功实现了色素特定吸收系数的分离。