免费文献传递   相关文献

Estimation of photosynthetic pigment of Festuca arundinacea using hyper-spectral data

利用高光谱数据快速估算高羊茅牧草光合色素的研究



全 文 :书利用高光谱数据快速估算高羊茅
牧草光合色素的研究
钱育蓉1,2,杨峰1,李建龙1,干晓宇1,杨齐1,王卫源3
(1.南京大学生命科学学院,江苏 南京210093;2.新疆大学软件学院,新疆 乌鲁木齐830046;3.新疆军区司令部机要处,新疆 乌鲁木齐830042)
摘要:光合色素是植物体进行初级生产的重要物质,能够间接反映植被的健康状况与光合能力,同时,高光谱遥感
为快速、大面积监测植被的色素变化提供了可能。本研究以高羊茅牧草为材料,实测了高羊茅冠层的高光谱反射
率与光合色素含量数据,对二者进行了相关分析,从5大类12个高光谱特征变量{犚犵、犚710、犇λ狉、犇λ700、犇[log(1/
λ犻)]、犛犵、犚犞犐犪、犚犞犐犫、犘犛犛犚、犘犛犖犇、犚犮犺、犆犃犚犐}中挑选了光合色素敏感参数,建立了植被指数光合色素估算模型。
结果表明,光合色素叶绿素a、叶绿素b、总叶绿素和类胡萝卜素与原始光谱进行相关性比较分析时,叶绿素b效果
最好;原始光谱、光谱一阶导和倒数对数一阶导3种光谱形式与光合色素进行比较分析时,光谱一阶导在700nm
附近与光合色素相关性最好,相关系数为-0.897;通过12个高光谱特征变量与光合色素相关性比较,选择达到极
显著水平的6个变量:犚犵、犇λ700、犇[log(1/犚730)]、犚犞犐犫、犚犮犺、犆犃犚犐进行光合色素含量的回归模型建立,这为利用物
理方法快速、无损伤探测高羊茅牧草的营养状况及牧草质量提供了理论支持。
关键词:高羊茅牧草;光合色素;高光谱数据;高光谱特征变量;回归模型
中图分类号:S688.4;S543+.9;Q945.11  文献标识码:A  文章编号:10045759(2009)04009409
  高羊茅(犉犲狊狋狌犮犪犪狉狌狀犱犻狀犪犮犲犪)是禾本科羊茅属的多年生草本植物,丛生型,须根发达,植株强健,具有耐寒、
耐热、耐潮湿、耐酸、耐干旱和产草量高等特点,是人工草原的优良牧草之一[1,2]。营养状况良好的高羊茅牧草产
草量高、草质优,对其提高载畜量和发挥生态功能均起着重要的作用,同时,由于植被光合色素(叶绿素a、叶绿素
b和类胡萝卜素)含量与其光合能力、发育阶段和氮素状况均有较好的相关性,因此,植被光合色素含量已经成为
评价植物长势和健康状况,检测光合作用、氮素状况和病害污染等胁迫的一种高效的植物健康指示剂[3,4]。
近年来,随着高光谱分辨率遥感的应用,使得利用反射光谱数据来估测植被叶绿素含量和氮素状况已经成为
植被遥感监测和遥感估产的一种重要手段,但其模型定量研究大多集中在玉米(犣犲犪犿犪狔狊)、水稻(犗狉狔狕犪狊犪狋犻
狏犪)、大豆(犌犾狔犮犻狀犲犿犪狓)等粮食作物上[5,6],且大多针对植被叶片或者单个植株进行,对草地群体冠层色素含量的
高光谱估算模型还鲜有报道。光合色素含量是反映植被群体光合作用强弱、植被长势的重要参量,它与叶面积指
数、生物量等参量共同反映了农作物的生长状况和生产能力。有研究表明,当叶片单位组成成分中干物质重量的
光合色素浓度相等,不同冠层色素浓度可能有很大的差异[7],而利用群体光谱数据估算植被叶绿素含量。一定程
度上克服了叶面积指数主要反映叶面积而忽略了植株状况的缺点;还弥补了叶片叶绿素含量只能反映单个植株
的长势而不能表达植被群体长势的不足。本研究结合高羊茅牧草的生长发育规律,从挑选高光谱特征变量和构
造光合色素拟合模型两方面着手深入研究了高羊茅牧草冠层的高光谱反射率和光合色素含量之间的相关关系和
估算模型,挑选能和光合色素有较好拟合的光谱指数,用于辐射光合效率和植被健康状态的估测,以此促进高光
谱遥感技术在牧草营养状况和草质监测中的应用。
1 材料与方法
1.1 材料
1.1.1 试验材料 供试材料是冷季型高羊茅牧草品种凌志(犉.犪狉狌狀犱犻狀犪犮犲犪cv.Barlexas),草种于2007年由
94-102
2009年8月
   草 业 学 报   
   ACTAPRATACULTURAESINICA   
第18卷 第4期
Vol.18,No.4
 收稿日期:20081028;改回日期:20081231
基金项目:国家高技术(863计划)专题项目(2007AA10Z231)资助。
作者简介:钱育蓉(1980),女,满族,山东武城人,讲师,博士。Email:QYR@xju.edu.cn
通讯作者。Email:jli2008@nju.edu.cn
克劳沃草业集团提供。
1.1.2 试验设计 试验在南京大学牧草试验地进行,采用盆栽施肥法进行不同施肥处理,盆径12cm,盆高10
cm,土壤为普通湿润砂质新成土,基础肥力为:有机质含量1.35%,水解氮71.25mg/kg,速效磷43.31mg/kg,
速效钾106.47mg/kg。选择健康的高羊茅草种于2008年3月15日播种,每盆均匀播种175~180粒种子。种
子出苗前隔天喷洒少量水,出苗后隔天浇1次透水,草种出齐苗后生长3周,进行施肥处理,肥料为常用牧草复合
肥(N∶P∶K=10∶9∶9),每周施肥1次,共施3次。试验设3个处理,1)对照处理:以施水为对照,记作CK;2)
肥1处理:根据花盆面积施5g/m2 的复合肥,记作F1;3)肥2处理:根据花盆面积施15g/m2 的复合肥,记作F2。
以上处理肥料均溶于150mL水中在早晨进行均匀施浇,每个处理4个重复。最后1次处理后1周内进行数据测
定。
1.2 测量方法
1.2.1 光谱测定 在晴朗无云无风或微风天气下(具体日期为4月28日、5月16日和5月29日),11:00-
13:00(太阳高度角为50°~55°),用地物光谱仪(ASD,FieldSpecRFR,384.7~1075.5nm,共512个波段,波段宽
度为1.3~1.5nm)分别测定不同施肥处理的高羊茅牧草冠层光谱反射率。测量时,高光谱仪视场角为15°,传感
器探头垂直向下,距冠层高度约20cm,探测面积直径约5.2cm。每个处理重复测9次,即在花盆的1个位置正
上方测3次,旋转花盆约120°后测3次,然后再旋转花盆约120°测3次,取平均值作为1个重复的光谱反射率值,
最后以3个重复的光谱反射率的平均值作为处理的光谱反射率值。
光谱测定是在高羊茅成坪后,盖度达到80%~95%时测定的。在室外太阳光下测量植被光谱反射率时,反
射率大小与太阳高度角有关,并且由于植被冠层不是绝对的漫反射体,而校正白板是标准漫反射体,测量时需要
及时进行标准白板校正,即每测完1个处理,用标准白板校正1次(标准白板反射率为1,因此所得目标物光谱为
无量纲的相对反射率)。
高羊茅群体冠层高光谱反射率谱线在不同的时期、温度、水肥条件下总体趋势大体保持一致,局部波段有差
异,谱线没有较大误差。由于叶绿素和类胡萝卜素特有的生化结构,其因电子跃迁造成的光谱吸收仅体现在可见
光区域,具体地说,叶绿素a(Chla)在380~450和675nm,叶绿素b(Chlb)在410~470和610nm,类胡萝卜素
(Cars)在470nm附近都有显著的吸收峰[7,8],赵祥等[9]的研究也表明,小麦(犜狉犻狋犻犮狌犿犪犲狊狋犻狏狌犿)冠层波谱与叶绿
素含量的回归反演计算模型在350~1060nm波段具有较高的反演精度,因此本研究只考虑植被色素与高光谱
反射率在400~1000nm波段的相关性。
1.2.2 光合色素测定 对应观测高羊茅冠层高光谱的叶片,随机剪切0.1g,按照丙酮∶无水乙醇∶蒸馏水=
4.5∶4.5∶1.0的混合液提取叶片色素,在波长663,646和470nm下测定光密度,然后计算叶绿素和类胡萝卜
素含量[10]。每个处理重复测3次,以其平均值作为色素含量。
1.2.3 高光谱数据特征参数及其提取 为了探究高羊茅牧草冠层反射光谱的特征,便于对其进行光谱匹配和分
解,以解释其光谱特征的物理学、生物化学、植物学和植物生理学的机理,需要提取目标物光谱的一些参数。常见
的高光谱参数有各类高光谱植被指数、各类高光谱数据变换形式(如对数变换、微分变换)构建的光谱参数、光谱
吸收指数、“蓝边”参数、“红边”参数、“绿峰”参数、红光吸收谷参数等。本研究所使用的高光谱变量主要有基于光
谱位置的变量、基于光谱面积的变量和基于高光谱植被指数的变量,其定义见表1。
2 结果与分析
2.1 光合色素差异分析
高羊茅经过施肥处理后(F1和F2)光合色素含量与未施肥处理(对照CK)相比显著提高(犘<0.05)(表2),
说明施肥处理后产生了明显的效果。中肥F1处理下,色素含量最高,随着施肥浓度的提高,光合色素Chla、
Chlb、Cars和叶绿素总量(Chls)的含量都有下降趋势。究其原因可能是施肥浓度过高,发生肥害,引起光合色素
浓度发生变化。
2.2 光合色素与高光谱特征数据的关系分析
2.2.1 光合色素与高光谱原始数据的相关性分析 Chla、Chlb、Chls(叶绿素总量)和Cars的含量与高羊茅冠层
59第18卷第4期 草业学报2009年
表1 高光谱特征变量的定义及描述
犜犪犫犾犲1 犜犺犲犱犲犳犻狀犻狋犻狅狀犪狀犱犱犲狊犮狉犻狆狋犻狅狀狅犳犺狔狆犲狉狊狆犲犮狋狉犪犾狋狔狆犻犮犪犾狏犪狉犻犪狀狋狊
光谱变量
Spectralvariant
名称
Name
定义和描述
Definitionanddescription
犚犻 光谱反射率Reflectance 原始光谱中波段犻处的反射率数值 Reflectanceatthespecifiedwavelength
犇λ犻 一阶导数光谱
Firstderivationof犚
波段犻处的一阶导数光谱数值,其导数光谱公式为:犚′=lim
Δλ→0
Δ犚
Δλ
,计算公式为犇λ犻=犚犻′=
犚犻+1-犚犻-1
λ犻+1-λ犻-1
Thefirstderivationof犚atthe犻wavelength
犇[log(1/λ犻)] 倒数对数一阶导
Thefirstderivativeofa
pseudoabsorbancespectral
signature
原始光谱中波段犻处的反射率数值倒数对数的一阶微分[7]Thefirstderivativeofapseudo
absorbancespectralsignatureatwavelength犻,log(1犚
)′=lim
Δλ→0
Δlog(1犚

Δλ
,为Itis:犇[log
(1
λ犻
)]=log(1犚犻
)′=
log(

犚犻+1
)-log(

犚犻-1

λ犻+1-λ犻-1
犇λ狉 红边幅值
Thefirstderivationofrededge
覆盖范围为680~760nm,犇λ狉是红边内一阶导数光谱的最大值Themaximumreflectance
offirstderivativespectrumbetween680and760nm
λ狉 红边位置 Rededgeposition 犇λ狉所对应的波长位置(nm)Wavelengthcorrespondingtothe犇λ狉
犚犵 绿峰反射率
Thereflectanceatgreenpeak
覆盖范围为510~560nm,犚犵 是绿光范围内最大的光谱反射率 Themaximumreflectance
between510and560nm
λ犵 绿峰位置 Greenpeakposition 犚犵 所对应的波长位置(nm)Wavelengthcorrespondingtothe犚犵
犚狉狑 红谷反射率
Thereflectanceatredvale
覆盖范围为650~690nm,犚狉狑是红光范围内最小的光谱反射率 Theminimumreflectance
between650and690nm
犛犵 绿峰面积 Theareaofthereflec
tancegreenpeak
510~560nm原始光谱曲线所包围的面积 Theareaofthereflectancebetween510~560nm
犚犞犐犻犼 比值植被指数
Ratiovegetationindex
犚犞犐犻犼=犚犻/犚犼,波段犻和犼的光谱反射率之比 Theratioofthereflectanceinwavelength犻
and犼nm
犘犛犛犚 色素简单比值指数
Pigmentspecificsimpleratio
犘犛犛犚=犚800犚680
[7]
犘犛犖犇 色素规一化差值指数 Pigment
specificnormalizeddifference
犘犛犖犇=犚800-犚680犚800+犚680
[7]
犚犮犺 叶绿素指数Chlorophylindex 犚犮犺=
犚640-犚673
犚673
[11]
犆犃犚犐 叶绿素吸收比值指数 Chlorophyl
absorptionratioindex
犆犃犚犐=犚700犚670
犆犃犚=犚700犚670
·670犪+犚670+犫
犪2槡 +1
,犪=犚700-犚550150
,犫=犚550-550犪 [12]
表2 不同施肥水平对高羊茅光合色素含量的影响
犜犪犫犾犲2 犈犳犳犲犮狋狊狅犳犱犻犳犳犲狉犲狀狋犳犲狉狋犻犾犻狕犲狉犾犲狏犲犾狅狀狆犺狅狋狅狊狔狀狋犺犲狋犻犮狆犻犵犿犲狀狋狊狅犳犉.犪狉狌狀犱犻狀犪犮犲犪狋狌狉犳 mg/g
色素Pigments
处理 Treatment
CK F1 F2
叶绿素a含量Chlacontent 1.549±0.090b 1.817±0.115a 1.783±0.114a
叶绿素b含量Chlbcontent 0.395±0.026b 0.495±0.032a 0.478±0.032a
类胡萝卜素含量Carscontent 0.418±0.016b 0.479±0.030a 0.475±0.033a
总叶绿素含量TotalChlscontent 1.944±0.116b 2.312±0.147a 2.262±0.146a
 注:不同小写字母表示0.05水平上差异显著。
 Note:Differentlowercaselettersmeansignificantdifferenceatthelevelof0.05.
69 ACTAPRATACULTURAESINICA(2009) Vol.18,No.4
光谱反射率的相关性在不同波段变化趋势一致,但相关性强度存在差异(图1)。光谱反射率与Chlb含量的相关
性最大,而与Cars含量相关性最差;原始光谱和光合色素在可见光波段500~710nm保持稳定的较高负相关性,
其相关性在波长450nm 处开始增加,并在波长540nm处达到持续稳定状态,在红光区域(710nm)的相关性最
大,而在波长710nm之后相关性下降,并在波长730nm左右达到最低点,然后开始呈正相关,在765nm处达到
稳定的正相关拟合度,并持续至1000nm;几种光合色素中,光谱反射率与Chlb的相关性最好,与Cars的相关
性最差,这一结果和宋开山等[13]利用小波分析所得出的结论相吻合。
高羊茅为多年生草本,丛生型,其冠层多层叶片间反射和多次散射现象在光谱测量过程中普遍存在,有研究
表明,随叶片叠加数增加,近红外光反射率增加,透射率降低;可见光反射率基本不变,透射率急剧降低[4]。基
于此有人指出红边与叶片叶绿素浓度之间有关系,还有人提出红边与牧草冠层叶绿素含量有关,进而提出红边是
测定叶绿素含量最好的遥感描述方法之一[14]。高羊茅冠层光谱曲线与叶绿素含量在波长680~760nm处的高
相关是因为此波段的光谱特征主要受叶绿素的吸收影响,而在730nm附近正处于叶绿素吸收和叶片反射以及
叶片的多次散射的交界处,因而相关性小,而在波长765nm之后,光谱数据和叶绿素含量都反映了叶片的反射
和多次散射的信息,所以二者的相关性较好。
2.2.2 光合色素与高光谱一阶导数的相关性分析 在植被高光谱的研究中,使用一阶导数已经成为一种减少低
频背景噪声和提高重叠光谱分辨率的常规技术,Blackburn[7]的研究表明光谱导数方法,特别是光谱反射率倒数
对数log(1/犚)的一阶导数和冠层叶片单位重量色素含量强相关。王秀珍等[14,15]的研究结果也证实了导数光谱
在消除背景信息的影响方面起到了很好的作用。根据上文Chlb与光谱反射率的相关性最佳,在此部分中选用
Chlb含量与原始光谱(犚)、光谱一阶导数[犇(犚)]、光谱倒数对数一阶导数{犇[log(1/犚)]}进行相关性比较分析。
高羊茅牧草冠层的原始光谱、光谱一阶导数、光谱倒数对数一阶导数与Chlb含量间的相关性变化比较复杂
(图2),但3种光谱形式在波长650~800nm具有较大的相关性,且其相关性较显著、稳定,可以用于构建色素含
量的模型。其中,冠层Chlb含量和光谱一阶导在700nm附近相关性最好,相关系数为-0.897。而与原始光谱
在710nm附近相关性较好,与光谱倒数对数一阶导在730~750nm相关性较好。光谱一阶导比原始光谱与
Chlb含量的相关性好,其结果可能与浦瑞良和宫鹏[16]、吴长山等[17]对森林冠层及作物群体导数光谱的分析结论
相类似。
图1 高羊茅冠层原始光谱与其光合色素的相关性分析
犉犻犵.1 犆狅狉狉犲犾犪狋犻狅狀犮狅犲犳犳犻犮犻犲狀狋犫犲狋狑犲犲狀狋犺犲狆犺狅狋狅狊狔狀狋犺犲狋犻犮狆犻犵犿犲狀狋犪狀犱狅狉犻犵犻狀犪犾
狊狆犲犮狋狉狌犿犻狀犉.犪狉狌狀犱犻狀犪犮犲犪犮犪狀狅狆狔
79第18卷第4期 草业学报2009年
图2 高羊茅冠层的犆犺犾犫含量与原始光谱(犚)、光谱一阶导数[犇(犚)]、光谱倒数对数一阶导
{犇[犾狅犵(1/犚)]}的相关性分析
犉犻犵.2 犆狅狉狉犲犾犪狋犻狅狀犮狅犲犳犳犻犮犻犲狀狋犫犲狋狑犲犲狀狋犺犲犆犺犾犫犮狅狀狋犲狀狋狊犪狀犱狅狉犻犵犻狀犪犾狊狆犲犮狋狉狌犿(犚),狋犺犲犳犻狉狊狋犱犲狉犻狏犪狋犻狏犲狊狆犲犮狋狉狌犿[犇(犚)]
犪狀犱狋犺犲犮狅狌狀狋犱狅狑狀狋狅狋犺犲狀狌犿犫犲狉狅犳犳犻狉狊狋犱犲狉犻狏犪狋犻狏犲{犇[犾狅犵(1/犚)]}犻狀犉.犪狉狌狀犱犻狀犪犮犲犪犮犪狀狅狆狔
2.3 光合色素与典型高光谱特征变量的相关性分析
2.3.1 典型高光谱特征变量的抽取 通常使用高光谱数据评价植被色素有2类常用方法:一是利用光谱数据、
植被指数、导数光谱评价植被色素含量;二是利用红边光学参数评价植被色素含量[18]。在高羊茅色素的估算中,
挑选典型的、与色素相关性较高的高光谱特征变量至关重要。吴长山等[5]和唐延林等[6]的研究揭示了玉米、水稻
等作物叶片和冠层中的色素浓度与红边位置、多种高光谱特征指数之间具有很好的相关性;Blackburn[7]的研究
表明,落叶树种冠层水平上单一色素含量和犘犛犛犚强相关,冠层色素含量和log(1/犚)的一阶导数也具有强相关
关系。
为了筛选典型的高光谱特征变量用于建模估算高羊茅群体冠层叶片色素含量,本研究选用以下5类高光谱
特征变量,①高光谱反射率:根据图1结果和已有试验数据,从犚狉狑、犚犵 和犚710中挑选了和色素含量相关系数较高
的犚犵 和犚710;②光谱一阶导数:根据图2结果和已有试验数据,从犇λ狉、犇λ700、犇[log(1/λ犻)]、犇λ犫(蓝边幅值)和
犇λ狔(黄边幅值)中挑选了与色素含量相关系数较高的犇λ狉、犇λ700和犇[log(1/λ犻)]三个典型光谱一阶导数变换;③
基于光谱面积的变量:根据实验数据从犛犵、犛狉(680~760nm的一阶导数光谱曲线所包围的面积)、犛狔(黄边面积,
560~640nm的一阶导数光谱曲线所包围的面积)和犛犫(蓝边面积,490~530nm的一阶导数光谱曲线所包围的
面积)中挑选了与色素含量相关系数最高的犛犵;④基于高光谱植被指数的变量:根据唐延林等[6]的研究结果,挑
选了2个光谱植被指数犚犞犐犪=犚800/犚550和犚犞犐犫=犚673/犚640;⑤各种色素指数:挑选了已知的多种色素指数,包括
色素简单比值指数(犘犛犛犚)、色素规一化差值指数(犘犛犖犇)、叶绿素指数犚犮犺、叶绿素吸收比值指数犆犃犚犐,其详细
定义见表1。
2.3.2 色素与典型高光谱特征变量的相关系数 通过相关系数筛选对高羊茅冠层光合色素含量敏感的典型高
光谱特征变量,结果见表3,所选用的植被指数与叶片叶绿素a、叶绿素b、总叶绿素及类胡萝卜素含量之间都具
有较好的相关性,抽取出6个相关性达到了0.01的极显著水平的变量:犚犵、犇λ700、犇[log(1/犚730)]、犚犞犐犫、犚犮犺、
犆犃犚犐,上述6个变量和光合色素的相关性强弱关系均表现为Chb>Chls>Cha>Cars。
89 ACTAPRATACULTURAESINICA(2009) Vol.18,No.4
表3 高羊茅群体冠层特征光谱植被指数与叶绿素、类胡萝卜素的相关系数
犜犪犫犾犲3 犆狅狉狉犲犾犪狋犻狅狀犮狅犲犳犳犻犮犻犲狀狋犫犲狋狑犲犲狀狋犺犲狋狔狆犻犮犪犾狏犲犵犲狋犪狋犻狅狀犻狀犱犻犮犲狊狅犳犮犪狀狅狆狔狊狆犲犮狋狉犪犪狀犱
犮犺犾狅狉狅狆犺狔犾犪狀犱犮犪狉狅狋犲狀狅犻犱犳狅狉犉.犪狉狌狀犱犻狀犪犮犲犪
色素
Pigment
犚犵 犚710 犇λ狉 犇λ700 犇[log(1/犚730)] 犛犵 犚犞犐犪 犚犞犐犫 犘犛犛犚 犘犛犖犇 犚犮犺 犆犃犚犐
Chla -0.739-0.744 0.599 -0.816-0.724-0.726 0.725 0.743 0.581 0.613 -0.747-0.773
Chlb -0.831-0.839 0.704 -0.897-0.833-0.818 0.829 0.831 0.681 0.717 -0.834-0.874
Cars -0.675 -0.681 0.591 -0.759-0.681 -0.663 0.685 0.708 0.542 0.569 -0.710-0.716
Chls -0.766-0.772 0.628 -0.840-0.754-0.753 0.755 0.768 0.609 0.642 -0.772-0.802
 注:和分别表示在0.05和0.01水平上达到显著;Chla、Chlb、Chls和Cars分别表示叶绿素a、叶绿素b、总叶绿素及类胡萝卜素含量。下
同。
 Note:and meanthesignificantatthelevelof0.05and0.01,respectively;Chla、Chlb、ChlsandCarsindicatethecontentsofchlorophyla,
chlorophylb,totalchlorophylandcarotenoid,respectively.Thesameasbelow.
上述结果中犇λ700和犚犞犐犫 中波段的选取是根据2008年4和5月的试验数据所得,但植被的不同物候期其
“红移”、“蓝移”现象普遍存在,因此,其典型高光谱变量的特征波段也会有所波动,犇λ700和犚犞犐犫 对植被色素含量
的指示性可能不具有普遍性。同时,表3结果也印证了犇λ狉(红边幅值)可能是单个叶片叶绿素含量较好的评价
器,而不是冠层叶绿素含量较好的评价器[18]。
2.4 色素的高光谱估算模型
2.4.1 Chb的单变量线性回归模型 在许多研究中使用高光谱遥感数据评价植被叶绿素含量的主要方法还是
利用回归方法(特别是逐步回归)找出与叶绿素含量相关的波段,然后建立这些波段(或者波段组合形式,包括植
被指数、微分光谱)与叶绿素含量之间的回归分析,建立各种回归模型。本研究采用筛选出的6个高光谱特征变
量,建立与它们相关性最好的Chlb的4种回归模型:
1.线性方程(Linear):犆犺=犪+犫狓
2.对数曲线(Logarithmic):犆犺=犪+犫ln狓
3.二次曲线(Auadratic):犆犺=犪+犫狓+犮狓2
4.S曲线:犆犺=e犪+犫/狓
如果方程采用犉检验,水准а=0.001,那么,犚犮犺和犚犵 的方程拟合不具有显著性。剩余4个变量中,与Chlb
拟合的犚2 关系表现为犇[log(1/犚730)]>犚犞犐犫>犆犃犚犐>犇λ700,其中拟合度较高的2个方程为Chlb=0.048-
0.269犇[log(1/犚730)]+0.26犇[log(1/犚730)]2 和Chlb=-0.976+7.579犚犞犐犫-7.757犚犞犐犫2,其犚2 均达到0.819
和0.896(表4)。用作逐步回归的4个模型中,二次曲线的拟合比其余3种模型的拟合度都高。
2.4.2 Chls、Chla和Cars的单变量线性回归模型 构造高羊茅冠层光合色素Chls、Chla和Cars的估测模型
时,利用和上面雷同的方法进行单变量线性拟合,结果如下:
①对Chls的拟合过程中,如果方程采用犉检验,水准а=0.001,那么,6个高光谱特征变量中拟合度较高的
2个为Chls=-1.855+2.3811犚犞犐犫-0.519犚犞犐犫2 和Chls=4.355-3.692犚犮犺+0.806犚犮犺2,其犚2 均达到0.846
和0.877。
②对Chla的拟合过程中,如果方程采用犉检验,水准а=0.001,那么,6个高光谱特征变量中拟合度较高的
2个为Chla=-2.185+3.391犚犞犐犫-0.939犚犞犐犫2 和Chla=4.866-5.259犚犮犺+1.459犚犮犺2,其犚2 均达到0.818
和0.851。
③对Cars的拟合过程中,如果方程采用犉检验,水准а=0.001,那么,拟合度较高的2个为Cars=-4.263
+21.697犚犞犐犫-22.861犚犞犐犫2 和Cars=8.014-33.317犚犮犺+35.145犚犮犺2,其犚2 均达到0.854和0.885。
④用作逐步回归的4个模型中,二次曲线的拟合比其余3种模型的拟合度都高。
99第18卷第4期 草业学报2009年
表4 高羊茅群体冠层叶片叶绿素犫含量与光谱特征植被指数的回归方程
犜犪犫犾犲4 犚犲犵狉犲狊狊犻狅狀犲狇狌犪狋犻狅狀狊狉犲犾犪狋犻狀犵狋犺犲狊狆犲犮狋狉犪狊狆犲犮犻犳犻犮狏犲犵犲狋犪狋犻狅狀犻狀犱犻犮犲狊狅犳犮犪狀狅狆狔犾犲犪狏犲狊狋狅
犆犺犾犫犮狅狀狋犲狀狋狊犳狅狉犉.犪狉狌狀犱犻狀犪犮犲犪
变量Variant 回归模型Regressionmodel 回归系数犚2 标准误差犛.犈. 犉检验犉test 显著性Sig.
犚犵
Chlb=0.716-0.06犚犵
=0.845-0.267ln犚犵
=0.565+0.008犚犵-0.007犚犵2
=e-1.414+2.65/犚犵
0.690 0.028 22.275 0.001
0.678 0.028 21.095 0.001
0.695 0.029 10.263 0.005
0.677 0.064 20.975 0.001
犇λ700
Chlb=17.731-20.138犇λ700
=1.366-9.084ln犇λ700
=50.031+163.431犇λ700-166.292犇λ7002
=e1.139+0.450/犇λ700
0.703 0.656 23.696 0.001
0.728 0.627 26.821 0.000
0.797 0.571 17.721 0.001
0.730 0.071 27.016 0.000
犚犞犐犫
Chlb=0.557+0.642犚犞犐犫
=1.081+0.292ln犚犞犐犫
=-0.976+7.579犚犞犐犫-7.757犚犞犐犫2
=e0.187-0.158/犚犞犐犫
0.691 0.022 22.346 0.001
0.727 0.020 26.963 0.000
0.896 0.013 38.820 0.000
0.768 0.022 30.080 0.000
犚犮犺
Chlb=0.582-0.728犚犮犺
=0.006-0.128ln犚犮犺
=0.538-0.244犚犮犺-1.216犚犮犺2
=e-1.071+0.045/犚犮犺
0.695 0.028 22.760 0.001
0.637 0.030 17.534 0.002
0.702 0.029 10.621 0.004
0.558 0.075 12.624 0.005
犆犃犚犐
Chlb=0.677-0.013犆犃犚犐
=1.116-0.232ln犆犃犚犐
=0.587-0.003犆犃犚犐+0.000犆犃犚犐2
=e-1.346+0.106/犆犃犚犐
0.763 0.024 32.218 0.000
0.750 0.025 30.009 0.000
0.766 0.026 14.731 0.001
0.746 0.057 29.439 0.000
犇[log(1/犚730)]
Chlb=-0.003-0.36犇[log(1/犚730)]
=-0.32-0.016ln{犇[log(1/犚730)]}
=0.048-0.269犇[log(1/犚730)]+0.26犇[log(1/犚730)]2

0.740 0.001 28.528 0.000
0.763 0.001 32.111 0.000
0.819 0.001 20.329 0.000
- - - -
 注:-表示无法拟合。
 Note:-indicatesthatnowaystofitting.
3 结论
本研究在高羊茅牧草关键生长期进行了3次试验,同步采集了高羊茅冠层光谱反射率与叶片色素含量数
据,分别对获取的数据进行了光合色素与原始光谱、导数光谱的相关分析,典型高光谱特征指数的挑选以及色素
的线性回归模拟,通过分析得出以下结论。
高羊茅牧草冠层光谱反射率与Chla、Chlb、Chls和Cars密切相关,其原始光谱和植被3种单一色素的相关
性在710nm处均达最高,且3种单一色素中,Chlb的相关性最好,Cars的相关性最差。
原始光谱、光谱一阶导数和光谱倒数对数一阶导在波长650~800nm与色素含量有较大的相关性,且其相
关性较显著、稳定。单就植被冠层Chlb含量而言,它和原始光谱在710nm附近相关性最好,和光谱一阶导在
700nm附近相关性最好,和光谱倒数对数一阶导在730~750nm最相关,其中和光谱一阶导的相关系数在700
nm附近达到最大,为-0.897,比原始光谱反射率与叶绿素含量的相关性更优。
从5大类高光谱变量中挑选了12个典型变量,与3种单一色素色素含量进行相关分析,得到相关性达到极
显著水平的6个变量,分别为犚犵、犇λ700、犇[log(1/犚730)]、犚犞犐犫、犚犮犺、犆犃犚犐,采用4种单变量线性回归模型进行植
001 ACTAPRATACULTURAESINICA(2009) Vol.18,No.4
被单一色素含量的模型拟合。
光合色素拟合结果表明,高光谱变量犇[log(1/犚730)]、犚犞犐犫 和犚犮犺与光合色素含量具有较好的拟合效果,且
4种拟合模型中二次曲线拟合模型相关性最高。利用高羊茅群体冠层光谱数据得到的光合色素含量的最优拟合
方程为Chlb=-0.976+7.579犚犞犐犫-7.757犚犞犐犫2,Chla=4.866-5.259犚犮犺+1.459犚犮犺2,Chls=4.355-3.692
犚犮犺+0.806犚犮犺2 和Cars=8.014-33.317犚犮犺+35.145犚犮犺2,其犚2 分别为0.896,0.851,0.877和0.885。
基于高光谱植被指数与作物生理参数,通过植被指数来建立的植被冠层色素的估算模型均有一定误差,误
差来源除了光谱仪自身的信噪比外,植被类型、生育阶段、生长条件以及测量环境等诸多因素都会影响统计模型
的建立。而基于辐射传输方程或几何光学原理建立的PROSPECT、SAIL等RT模型虽然有深厚的物理基础,
但是模型的参数过于复杂,有些参数无法实时获取,因而限制了其应用性。
参考文献:
[1] 王小利,刘正书,牟琼,等.高羊茅遗传多样性RAPD分析[J].草业学报,2007,16(4):8286.
[2] 杨志民,陈煜,韩烈保.不同光照强度对高羊茅形态和生理指标的影响[J].草业学报,2007,16(6):2329.
[3] FilelaD,PerruelasJ.Therededgepositionandshapeasindicatorsofplantchlorophylcontent,biomassandhydricstatus[J].Inter
nationalJournalofRemoteSensing,1994,15(7):14591470.
[4] MadeiraAC,MendoncaA,FerreiraME,犲狋犪犾.Relationshipbetweenspectroradiometricandchlorophyl measurementsin
greenBeans[J].CommunicationinSoilScienceandPlantAnalysis,2000,31(56):631643.
[5] 吴长山,童庆禧,郑兰芬,等.水稻、玉米的光谱数据与叶绿素的相关分析[J].应用基础与工程科学学报,2000,8(1):31
37.
[6] 唐延林,黄敬峰,王秀珍,等.玉米叶片高光谱特征及与叶绿素、类胡萝卜素相关性的研究[J].玉米科学,2008,16(2):7176.
[7] BlackburnGA.Quantifyingchlorophylsandcaroteniodsatleafandcanopyscales:Anevaluationofsomehyperspectralap
proaches[J].RemoteSensingofEnvironment,1998,66:273285.
[8] 程乾,黄敬峰,王人潮,等.MODIS植被指数与水稻叶面积指数及叶片叶绿素含量相关性研究[J].应用生态学报,2004,15
(8):13631367.
[9] 赵祥,刘素红,王培,等.基于高光谱数据的小麦叶绿素含量反演[J].地理与地理信息科学,2004,20(3):3639.
[10] 高俊风.植物生理学实验指导[M].北京:高等教育出版社,2006.101103.
[11] GuptaRK,VijayanD,PrasadTS.Newhyperspectralvegetationcharacterizationparameters[J].AdvancesinSpaceRe
search,2001,28(1):201206.
[12] BorgeNH,LebiancE.Comparingpredictionpowerandstabilityofbroadbandandhyperspectralvegetationindicesforesti
mationofgreenleafareaindexandcanopycholorophyldensity[J].RemoteSensingofEnvironment,2000,76:156172.
[13] 宋开山,张柏,王宗明,等.小波分析在大豆叶绿素含量高光谱反演中的应用[J].中国农学通报,2008,22(9):101108.
[14] 王秀珍,王人潮,黄敬峰.微分光谱遥感及其在水稻农学参数测定上的应用研究[J].农业工程学报,2002,18(1):913.
[15] 王秀珍,王人潮,李云梅,等.不同氮素营养水平的水稻冠层光谱红边参数及其应用研究[J].浙江大学学报(农业与生命科学
版),2001,27(3):301306.
[16] 浦瑞良,宫鹏.森林生物化学与CASI高光谱分辨率遥感数据的相关分析[J].遥感学报,1997,1(2):115123.
[17] 吴长山,项月琴,郑兰芬,等.利用高光谱数据对作物群体叶绿素密度估算的研究[J].遥感学报,2000,4(3):228232.
[18] 张金恒,王珂,王人潮.高光谱评价植被叶绿素含量的研究进展[J].上海交通大学学报(农业科学版),2003,21(1):74
80.
101第18卷第4期 草业学报2009年
犈狊狋犻犿犪狋犻狅狀狅犳狆犺狅狋狅狊狔狀狋犺犲狋犻犮狆犻犵犿犲狀狋狅犳犉犲狊狋狌犮犪犪狉狌狀犱犻狀犪犮犲犪狌狊犻狀犵犺狔狆犲狉狊狆犲犮狋狉犪犾犱犪狋犪
QIANYurong1,2,YANGFeng1,LIJianlong1,GANXiaoyu1,YANGQi1,WANGWeiyuan3
(1.SchoolofLifeScience,NanjingUniversity,Nanjing210093,China;2.SoftwareColege,Xinjiang
University,Urumqi830046,China;3.ConfidentialDepartment,Xinjiang
MilitaryAreaCommand,Urumqi830042,China)
犃犫狊狋狉犪犮狋:Photosyntheticpigmentsofvegetationareaprimaryproductandimportantmaterialsthatcanindi
rectlyreflecthealthstatusandphotosyntheticcapacityofvegetation.Hyperspectralremotesensingprovides
thepossibilityforrapid,largescalemonitoringofvegetationchlorophylchange.Inthisstudy,canopyreflec
tancespectraandchlorophylcontentofcoolseason犉犲狊狋狌犮犪犪狉狌狀犱犻狀犪犮犲犪weremeasuredandtheirrelationship
analyzed.Thetopfivechlorophylsensitiveparameterswereselectedfrom12hyperspectralcharacteristicvaria
bles,andthenusedtoestablishavegetationindexmodeltoestimatephotosyntheticpigments.1)Therelation
shipbetweenChlbandtheoriginalspectrumwasthebestamongthephotosyntheticpigments:Chla,Chlb,
Chls,andCars.2)Thefirstderivativespectragaveabetterrelationship(correlationcoefficient:-0.897)
withvegetationphotosyntheticpigmentsnear700nmwavelengththantheothertwohyperspectrumforms{R
and犇[log(1/犚)]}.Finaly,aftercomparingthecorrelationbetweenphotosyntheticpigmentconcentrations
andhyperspectraldata,themostsignificantvariables{犚犵,犇λ700,犇[log(1/犚730)],犚犞犐犫,犚犮犺,犆犃犚犐}were
selectedfromtwelvevariablestoestablishtheregressionmodelofphotosyntheticpigments.Alofthesepro
videatheoreticalbasisforrapid,noninvasivedetectionofnutritionalstatusandmeadowqualityof犉.犪狉狌狀犱犻
狀犪犮犲犪.
犓犲狔狑狅狉犱狊:犉犲狊狋狌犮犪犪狉狌狀犱犻狀犪犮犲犪;photosyntheticpigment;hyperspectraldata;typicalvariantofhyperspec
trum;regressionmodel
201 ACTAPRATACULTURAESINICA(2009) Vol.18,No.4