免费文献传递   相关文献

Research on the hyperspectral remote sensing estimation models for the fresh yield of alfalfa grassland

苜蓿人工草地高光谱遥感估产模型的研究



全 文 :书苜蓿人工草地高光谱遥感估产模型的研究
吕小东1,王建光1,孙启忠2,姚贵平3,高凤芹2
(1.内蒙古农业大学生态环境学院,内蒙古 呼和浩特010019;2.中国农业科学院草原研究所,内蒙古 呼和浩特010010;
3.内蒙古农业大学理学院,内蒙古 呼和浩特010019)
摘要:实地测量了10个苜蓿品种在不同生长时期的冠层高光谱数据,用以消除不同品种间光谱的差异性,并以多
个波段的反射率、一阶导数和多种光谱吸收特征参数为光谱参数,运用多种单变量回归模型,对人工苜蓿草地的鲜
草产量进行了估算。结果表明,复合、指数等非线性模型要优于线性模型,在同类型的光谱参数中,线性模型决定
系数高的参数,其二次型、三次型多项式和复合、乘幂、指数等非线性回归模型的决定系数通常也比较高,且通常高
于线性模型;诸多估算模型中,以747nm处一阶导数为自变量的复合、指数2种形式的估算模型,其相关系数最高
为狉=0.852,均方根误差为0.466kg/m2,相对误差为21.14%,其估算精度最高,可作为多个苜蓿品种统一使用的
鲜草产量高光谱估算模型。
关键词:苜蓿;品种;鲜草产量;高光谱遥感;回归模型
中图分类号:S812.3;S551+.7;S127  文献标识码:A  文章编号:10045759(2014)01008408
犇犗犐:10.11686/cyxb20140111  
  目前,高光谱分辨率遥感(hyperspectralremotesensing)技术在植被上的应用与发展,已经使植被监测和研
究进入到了精准化、定量化和机理化的新阶段。植被光谱与植物品种、植株密度、冠层结构、叶片形状、叶组织结
构、植物生化组分及比例、光谱测量条件(如气象条件、光谱仪分辨率、测量日期、背景)等因素有关[1]。所以,可以
应用高光谱遥感数据对重要的植物生长信息(如覆盖度、叶面积、生物量、叶绿素含量等)进行反演,拟合出植物生
长参数的定量模型,从而实现对植被冠层快速、有效、非接触、无破坏的野外信息采集与处理,加强对植被长势的
实时诊断研究,乃是国内外农业遥感学者们的研究热点之一[2]。
在草地资源领域,高光谱研究以其快速、便捷、准确等优势,受到世界各国学者的普遍关注。目前,国外已将
高光谱遥感技术广泛应用于草地植被生化参数的无损伤估算牧草品质参数测定、含水率估测、退化草地指示、草
种识别等领域[3]。我国应用高光谱,还主要集中在小麦(犜狉犻狋犻犮狌犿犪犲狊狋犻狏狌犿)、水稻(犗狉狔狕犪狊犪狋犻狏犪)、玉米(犣犲犪
犿犪狔狊)、棉花(犌狅狊狊狔狆犻狌犿犫犪狉犫犪犱犲狀狊犲)、大豆(犌犾狔犮犻狀犲犿犪狓)等大宗农作物长势的监测和比较,较广泛地应用于主
要产粮县旱地和水田作物的物候期监测,在草地资源上应用还不是很多,人工牧草地植被监管方面更少。
在国外,Kawamura等[4]通过高光谱来预测禾草与白三叶(犜狉犻犳狅犾犻狌犿狉犲狆犲狀狊)混播草地中豆科植物的组分含
量。Zhao等[5]发现波段1145和1205nm 的反射率可以用于狗牙根(犆狔狀狅犱狅狀犱犪犮狋狔犾狅狀)草地生物量的估算。
Mutanga等[6]利用草地的高光谱吸收特征来估测牧草品质。Post等[7]以多种禾草与苜蓿(犕犲犱犻犮犪犵狅狊犪狋犻狏犪)混
播地为对象,发现红光区域内、红光和近红外区域之间的导数光谱,与牧草营养品质和生物量显示出的相关性为
最高。国内的研究也认为,光谱一阶导数在一定程度上能较好的消除土壤环境的影响,并结合线性、对数、指数、
乘幂和二次、三次多项式等多种单变量函数回归模型,来估算植被的生物量、叶面积等[8]。张凯等[9]以723nm
处一阶导数的对数模型来反演甘南草地地上生物量;杨红丽等[10]应用高光谱探讨了多花黑麦草(犔狅犾犻狌犿犿狌犾狋犻
犳犾狅狉狌犿)植株含氮量的估测模型;纳钦[11]对苜蓿和缘毛雀麦(犅狉狅犿狌狊犮犻犾犻狋狌狊)混播地中营养成分的高光谱估算模
型进行了研究。周宇庭等[12]发现,在估算藏北典型高寒草甸地上生物量的多种回归函数中,以乘幂函数的预测
84-91
2014年2月
   草 业 学 报   
   ACTAPRATACULTURAESINICA   
第23卷 第1期
Vol.23,No.1
收稿日期:20130307;改回日期:20130422
基金项目:现代农业产业体系建设专项资金(CARS35)和科技部“十二五”国家科技支撑计划项目(2012BAD12B09)资助。
作者简介:吕小东(1984),男,河北涞源人,在读硕士。Email:lxd635@163.com
通讯作者。Email:wangjg8580@163.com,sunqz@126.com
效果最好。张艳楠等[13]认为,随着草地生物量范围的增大,估算模型逐渐由线性趋近于指数模型,当生物量较大
时以指数函数为佳。
一直以来,对农学参数与光谱特性的研究,大多以单一品种为材料,而未考虑到不同品种间的光谱特性差异
对估算模型的影响,从而在生产应用上造成一定的局限性。因此,本试验以10个苜蓿品种在不同时期的冠层高
光谱和鲜草产量数据为材料,统一分析苜蓿鲜草产量的高光谱特性,并尝试建立关于对多个苜蓿品种统一使用的
鲜草产量高光谱估算模型。
1 材料与方法
1.1 试验地概况
试验地位于内蒙古包头市土默特右旗萨拉齐镇(东经110.5°、北纬40.5°、海拔1067.2m)的内蒙古农业大学
职业技术学院科技园区。全年少雨雪,降水多集中在7-8月,年均降水量339.8mm;日温差较大,年均气温
7.1℃,无霜期年均132d;土壤为栗钙壤土,pH值7.80,有机质15.97g/kg,有效氮347.95mg/kg,有效磷54.33
mg/kg,有效钾327.67mg/kg;属典型大陆性半干旱季风气候。
1.2 田间试验设计
以苜蓿王(Alfaking)、大富豪(Milionaire)、农宝(FarmersTreasure)、费纳尔(Vernal)、阿尔冈金(Algon
quin)、金皇后(GoldenEmpress)、WL232、WL232HQ、WL323等9个国外引进品种及1个国产品种敖汉
(Aohan)为供试材料,采用随机区组设计,重复3次,小区数共30个,每个小区面积为160m2(4m×40m)。于
2009年7月15日播种,南北向条播种植,行距45cm,播深2cm,播量12kg/hm2。
1.3 数据采集方法
分别于2012年6月18日、8月4日、8月21日、9月6日,采用美国光谱分析仪器公司(AnalyticalSpectral
Device)生产的ASDFieldSpecHandHeld手持式高光谱仪,在光谱325~1075nm波长范围内对生长第4年的各
苜蓿品种处理小区进行连续测量,采样间隔为1.5nm,光谱分辨率3.5nm,光纤探头视场角(FOV)10°。所有的
测试均选择在晴朗无风天气进行,每次测定时间为北京时间10:00-14:00时。观测时传感器探头垂直向下,以
苜蓿冠层为基准,与冠层顶相距约1m。在每个处理小区里,选定生长均匀且长为1m、宽为3行的3个样方测量
其冠层光谱值,每个样方采集3组光谱值,以其平均值作为该样方的光谱反射值。测量过程中及时进行标准白板
校正。
每次光谱数据采集后,立即将对应样方内苜蓿鲜草按留茬5cm进行刈割,并立即装入密闭自封袋中,带回实
验室称重并计算出单位面积内的鲜草产量(kg/m2)。
1.4 数据分析方法
以第1、第3重复区数据为建模数据集(狀=240)用以建立估算模型,以第2重复区数据为验证数据集(狀=
120)用以检验估算模型。采用冠层光谱原始反射率、一阶导数及多种光谱吸收特征参数作为冠层光谱特征参数,
研究苜蓿草地高光谱数据与鲜草产量数据之间的相关性,然后运用单变量线性和非线性回归方法建立估算模型,
并通过检验最终确立苜蓿鲜草产量的估算模型。
采用 MicrosoftExcel2003、MicrosoftExcel2007和SPSS18.0等软件,对试验所得数据进行图表绘制和相
关分析、回归分析等数理统计。
2 结果与分析
2.1 苜蓿鲜草产量性状的高光谱特征分析
2.1.1 苜蓿鲜草产量与反射率的相关性分析 采用相关分析方法,得到冠层反射率与苜蓿鲜草产量的相关系数
曲线(图1)。
从图1看出,在338~757nm波段的冠层反射光谱与鲜草产量呈负相关,其中在445~684nm波段处形成1
个相关系数的高平台区,并在615nm波段达到最高,为-0.791,在615~642nm形成一个小波谷,其中以633
58第23卷第1期 草业学报2014年
nm波段相关系数最低,为-0.730,都达到极显著水平(犘<0.01)。在520~600nm的绿光波段中,以600nm
处相关系数最高,为-0.778;在630~690nm的红光波段中,以645nm处相关系数最高,为-0.774,均达到极
显著水平(犘<0.01)。在758~1075nm波段区间,冠层光谱反射率与鲜草产量呈正相关,其中在近红外758~
770nm波段,相关系数迅速上升,在近红外波段760nm处达到最大(狉=0.582),后急速下降,在770nm达到最
低(狉=0.149)之后开始缓慢上升,并在813~890nm波段处形成一个相关系数的高平台区,均在0.550以上,之
后开始缓慢下降。
2.1.2 苜蓿鲜草产量与一阶导数光谱变量的相关性分析  一般认为,对反射率采用微分处理提取一阶导数,可
以部分地去除线性或接近线性的背景和噪声光谱对植被光谱(必须为非线性)的影响,能增强光谱曲线在坡度上
的细微变化[14]。图2反映了苜蓿鲜草产量与一阶导数光谱变量间的相关性。
图1 苜蓿鲜草产量与反射率的相关性分析
犉犻犵.1 犆狅狉狉犲犾犪狋犻狅狀犫犲狋狑犲犲狀狋犺犲犳狉犲狊犺狔犻犲犾犱狅犳
犪犾犳犪犾犳犪犪狀犱狉犲犳犾犲犮狋犪狀犮犲 
图2 苜蓿鲜草产量与一阶导数的相关性
犉犻犵.2 犆狅狉狉犲犾犪狋犻狅狀犫犲狋狑犲犲狀狋犺犲犳狉犲狊犺狔犻犲犾犱狅犳
犪犾犳犪犾犳犪犪狀犱狋犺犲犳犻狉狊狋犱犲狉犻狏犪狋犻狏犲 
由图2可知,在400~453nm、487~553nm和637~650nm、688~710nm、917~929nm、947~959nm波
段为负相关,以波长552nm 处相关系数最大(狉=-0.806),其次为443nm(狉=-0.803),均达到了极显著水平
(犘<0.01);在461~470nm、554~636nm、735~761nm、771~817nm波段为正相关,最大相关系数波长位于
746nm处(狉=0.847),达到了极显著水平(犘<0.01)。
2.1.3 苜蓿鲜草产量与不同吸收特征参数的相关性分析 常见的高光谱遥感吸收特征变量包括从原始光谱、一
阶导数光谱提取的基于光谱位置的特征变量、基于高光谱面积的特征变量和基于高光谱植被指数的吸收特征参
数,见表1。
如表1所示苜蓿鲜草产量与光谱吸收特征参数的相关系数,都达到了0.01检验水平。基于光谱位置的几种
特征参数,其与苜蓿鲜草产量的相关系数普遍较高,都在0.7以上,以犇y 最高(狉=0.776);而基于高光谱面积的
特征参数中,仅有2个参数的相关系数较高,其中犛犚g的相关系数最高,其次为犛犇b,都大于0.7;基于高光谱植
被指数的特征参数中,除了犚g/犚r、犛犇r/犛犚r、(犚g-犚r)/(犚g+犚r)、(犛犇r-犛犚r)/(犛犇r+犛犚r)等参数的相关系数
很低外,其余特征参数均较高,其中尤以(犛犇r-犛犚g)/(犛犇r+犛犚g)最高,超过0.8,(犛犇r-犛犇y)/(犛犇r+犛犇y)
和犛犇r/犛犚g的相关系数也都接近0.8。
2.2 模型建立
在同类型光谱参数中,选取相关系数较大的光谱参数,作为建模参数,与苜蓿鲜草产量间进行单变量回归分
析;其选用的回归方程形式包括线性、对数、逆模型、二次型、三次型多项式和复合、乘幂、S型、指数等9种。
2.2.1 应用反射率建模 从反射率中,挑选相关系数大于0.8的波段参数,进行回归分析,其决定系数见表2。
由表2看出,不同波段反射率参数中,线性回归模型决定系数高的参数,其非线性模型的决定系数也比较高,
68 ACTAPRATACULTURAESINICA(2014) Vol.23,No.1
以615nm反射率为自变量的回归方程的决定系数较高,并以其复合、指数2种形式为最高,达到0.788。在众多
回归形式中,除逆模型和对数、S型的回归形式外,其余非线性回归模型都比线性模型的决定系数要高;其中复
合、指数最高,其次为乘幂。
2.2.2 应用一阶导数建模 从一阶导数中,挑选相关系数大于0.8的波段参数,进行回归分析,其决定系数见表3。
表1 苜蓿鲜草产量与光谱吸收特征参数间的相关性
犜犪犫犾犲1 犆狅狉狉犲犾犪狋犻狅狀犫犲狋狑犲犲狀狋犺犲犳狉犲狊犺狔犻犲犾犱狅犳犪犾犳犪犾犳犪犪狀犱狋犺犲狊狆犲犮狋狉犪犾犪犫狊狅狉狆狋犻狅狀犳犲犪狋狌狉犲狆犪狉犪犿犲狋犲狉狊
类型Variety 参数Parameter 定义与算法 Definitionandalgorithm 相关系数狉
基于光谱位置变量
Basedonspectral
position
犇b 蓝边内一阶导数光谱中的最大值,蓝边覆盖490~530nm。Thefirstorder
maximumderivativein490-530nm.
-0.706
犇y 黄边内一阶导数光谱中的最大值,黄边覆盖550~582nm。Thefirstorder
maximumderivativein550-582nm.
-0.776
犇r 红边内一阶导数光谱中的最大值,红边覆盖680~780nm。Thefirstorder
maximumderivativein680-780nm.
-0.737
犚g 绿峰反射率,即波长510~560nm 范围内最大的波段反射率。Themaxi
mumreflectancein510-560nm.
-0.750
犚r 红谷反射率,即波长640~680nm范围内最小的波段反射率。Theminimum
reflectancein640-680nm.
-0.770
基于光谱面积变量
Basedonspectral
area
犛犇b 蓝边面积,即蓝边内一阶导数的总和,蓝边覆盖490~530nm。Summation
ofthefirstorderderivativesin490-530nm.
-0.724
犛犇y 黄边面积,即黄边内一阶导数的总和,黄边覆盖550~582nm。Summationof
thefirstorderderivativesin550-582nm.
0.514
犛犇r 红边面积,即红边内一阶导数的总和,红边覆盖680~780nm。Summation
ofthefirstorderderivativesin680-780nm.
0.499
犛犚g 绿峰内(波长510~560nm范围内)一阶导数波段值总和。Summationofthe
firstorderderivativesin510-560nm.
-0.744
犛犚r 红谷内(波长640~680nm范围内)一阶导数波段值总和。Summationofthe
firstorderderivativesin640-680nm.
-0.257
基于植被指数
Basedonvegetation
index
犚g/犚r 绿峰反射率与红谷反射率的比值。Radioof犚gand犚r. 0.354
犛犇r/犛犇b 红边内一阶导数的总和与蓝边内一阶导数的总和的比值。Radioof犛犇rand犛犇b. 0.764
犛犇r/犛犇y 红边内一阶导数的总和与黄边内一阶导数的总和的比值。Radioof犛犇rand犛犇y.-0.724
犛犇r/犛犚g 红边内一阶导数总和与绿峰内一阶导数总和的比值。Radioof犛犇rand犛犚g. 0.791
犛犇r/犛犚r 红边内一阶导数总和与红谷内一阶导数总和的比值。Radioof犛犇rand犛犚r. -0.279
(犚g-犚r)/(犚g+犚r) 绿峰反射率与红谷反射率的归一化值。Normalizeddifference犚gand犚r. 0.388
(犛犇r-犛犇b)/(犛犇r+犛犇b) 红边内一阶导数的总和与蓝边内一阶导数的总和的归一化值。Normalized
difference(犛犇r-犛犇b)and(犛犇r+犛犇b).
0.724
(犛犇r-犛犇y)/(犛犇r+犛犇y) 红边内一阶导数的总和与黄边内一阶导数的总和的归一化值。Normalized
difference(犛犇r-犛犇y)and(犛犇r+犛犇y).
-0.795
(犛犇r-犛犚g)/(犛犇r+犛犚g) 红边内一阶导数总和与绿峰内一阶导数总和的归一化值。Normalizeddiffer
ence(犛犇r-犛犚g)and(犛犇r+犛犚g).
0.821
(犛犇r-犛犚r)/(犛犇r+犛犚r) 红边内一阶导数总和与红谷内一阶导数总和的归一化值。Normalizeddiffer
ence(犛犇r-犛犚r)and(犛犇r+犛犚r).
0.228
 注:和分别表示通过0.05和0.01水平的显著性检验。下同。
 Note:andindicate0.05and0.01significantlevel,respectively.Thesamebelow.
78第23卷第1期 草业学报2014年
表2 苜蓿鲜草产量与反射率间回归分析的决定系数(狀=240)
犜犪犫犾犲2 犜犺犲犮狅犲犳犳犻犮犻犲狀狋狅犳犱犲狋犲狉犿犻狀犪狋犻狅狀狅犳狉犲犵狉犲狊狊犻狅狀犪狀犪犾狔狊犻狊犫犲狋狑犲犲狀狋犺犲犳狉犲狊犺狔犻犲犾犱狅犳犪犾犳犪犾犳犪犪狀犱狋犺犲狉犲犳犾犲犮狋犪狀犮犲(狀=240)
反射率
Reflectance
线性
Linear
对数
Logarithmic
逆模型
Inverse
二次型
Quadratic
三次型
Cubic
复合
Compound
乘幂
Power
S型

指数
Exponential
R614nm 0.626 0.586 0.463 0.627 0.644 0.787 0.697 0.507 0.787
R615nm 0.626 0.588 0.468 0.627 0.644 0.788 0.700 0.514 0.788
R616nm 0.625 0.586 0.463 0.626 0.644 0.788 0.698 0.510 0.788
表3 苜蓿鲜草产量与一阶导数间回归分析的决定系数(狀=240)
犜犪犫犾犲3 犜犺犲犮狅犲犳犳犻犮犻犲狀狋狅犳犱犲狋犲狉犿犻狀犪狋犻狅狀狅犳狉犲犵狉犲狊狊犻狅狀犪狀犪犾狔狊犻狊犫犲狋狑犲犲狀狋犺犲犳狉犲狊犺狔犻犲犾犱狅犳
犪犾犳犪犾犳犪犪狀犱狋犺犲犳犻狉狊狋犱犲狉犻狏犪狋犻狏犲狊狆犲犮狋狉狌犿(狀=240)
一阶导数
Derivation
线性
Linear
对数
Logarithmic
逆模型
Inverse
二次型
Quadratic
三次型
Cubic
复合
Compound
乘幂
Power
S型

指数
Exponential
D433nm 0.644 - - 0.656 0.658 0.786 - - 0.786
D451nm 0.640 - - 0.640 0.640 0.775 - - 0.775
D543nm 0.643 - - 0.645 0.645 0.699 - - 0.699
D552nm 0.649 - - 0.658 0.658 0.801 - - 0.801
D746nm 0.717 - - 0.717 0.728 0.848 - - 0.848
D747nm 0.716 - - 0.716 0.721 0.831 - - 0.831
D748nm 0.706 - 0.561 0.708 0.716 0.820 - 0.701 0.820
D749nm 0.702 - 0.529 0.714 0.722 0.811 - 0.677 0.811
D750nm 0.713 - - 0.720 0.720 0.823 - - 0.823
D751nm 0.697 - 0.542 0.699 0.708 0.799 - 0.695 0.799
D752nm 0.716 - - 0.717 0.729 0.810 - - 0.810
D753nm 0.714 - - 0.715 0.726 0.806 - - 0.806
不同波段的一阶导数参数中,线性回归模型决定系数高的参数,其二次型、三次型多项式和复合、指数等非线
性模型的决定系数也比较高,以746nm处一阶导数为自变量的回归方程的决定系数较高,并以其复合、指数2
种形式为最高,达到0.848。在不同的回归形式中,二次型、三次型多项式和复合、指数等非线性回归模型的决定
系数,都比线性模型要高;其中复合、指数的决定系数最高。
2.2.3 应用高光谱吸收特征参数建模 基于光谱位置的特征变量、基于高光谱面积的特征变量和基于高光谱植
被指数的3种吸收特征参数,作为自变量,分别与苜蓿鲜草产量进行回归分析,其决定系数见表4。
由表4可知,基于光谱位置的几种特征参数,其回归模型的决定系数普遍较高,以犇y 最高;而基于高光谱面
积的特征参数中,犛犚g的决定系数最高,其次为犛犇b;基于高光谱植被指数的特征参数中,犚g/犚r、犛犇r/犛犚r、(犚g-
犚r)/(犚g+犚r)、(犛犇r-犛犚r)/(犛犇r+犛犚r)等参数的决定系数很低外,其余特征参数均较高,以(犛犇r-犛犚g)/(犛犇r
+犛犚g)最高,其次为(犛犇r-犛犇y)/(犛犇r+犛犇y)和犛犇r/犛犚g。
从不同的回归方程中,基于光谱位置的几种特征参数,二次型、三次型多项式和复合、乘幂、指数等非线性回
归模型的决定系数,都比线性模型要高;而基于高光谱面积的特征参数,二次型、三次型多项式和复合、乘幂、指数
等非线性回归模型的决定系数,都比线性模型要高;基于高光谱植被指数的特征参数,二次型、三次型多项式和复
合、乘幂、指数等非线性回归模型的决定系数,都比线性模型要高,此外,多个参数的对数、S型回归模型的决定系
数也比线性模型要高。综合来看,在同类型的参数中,线性模型决定系数高的参数,往往其二次型、三次型、复合、
乘幂、指数等非线性回归模型的决定系数也比较高,且比其线性模型要高;其中复合、指数的决定系数最高,其次
为三次型、二次型多项式。
88 ACTAPRATACULTURAESINICA(2014) Vol.23,No.1
表4 苜蓿鲜草产量与光谱吸收特征参数间回归分析的决定系数(狀=240)
犜犪犫犾犲4 犜犺犲犮狅犲犳犳犻犮犻犲狀狋狅犳犱犲狋犲狉犿犻狀犪狋犻狅狀狅犳狉犲犵狉犲狊狊犻狅狀犪狀犪犾狔狊犻狊犫犲狋狑犲犲狀狋犺犲犳狉犲狊犺狔犻犲犾犱狅犳犪犾犳犪犾犳犪犪狀犱
狋犺犲狊狆犲犮狋狉犪犾犪犫狊狅狉狆狋犻狅狀犳犲犪狋狌狉犲狆犪狉犪犿犲狋犲狉狊(狀=240)
参数
Parameter
线性
Linear
对数
Logarithmic
逆模型
Inverse
二次型
Quadratic
三次型
Cubic
复合
Compound
乘幂
Power
S型

指数
Exponential
犇b 0.498 0.439 0.332 0.504 0.525 0.542 0.459 0.311 0.542
犇y 0.602 - - 0.674 0.674 0.666 - - 0.666
犇r 0.543 0.538 0.443 0.572 0.621 0.625 0.627 0.517 0.625
犚g 0.562 0.482 0.355 0.580 0.599 0.649 0.541 0.375 0.649
犚r 0.592 0.545 0.394 0.597 0.625 0.783 0.686 0.473 0.783
犛犇b 0.523 0.457 0.344 0.533 0.552 0.587 0.496 0.340 0.587
犛犇y 0.265 - 0.214 0.265 0.272 0.227 - 0.191 0.227
犛犇r 0.249 0.276 0.271 0.285 0.286 0.390 0.445 0.453 0.390
犛犚g 0.554 0.486 0.359 0.558 0.578 0.605 0.504 0.334 0.605
犛犚r 0.066 - - 0.079 0.137 0.178 - - 0.178
犚g/犚r 0.125 0.144 0.153 0.147 0.164 0.231 0.282 0.313 0.231
犛犇r/犛犇b 0.584 0.652 0.620 0.618 0.687 0.630 0.790 0.830 0.630
犛犇r/犛犇y 0.524 - 0.631 0.607 0.646 0.570 - 0.741 0.570
犛犇r/犛犚g 0.626 0.689 0.639 0.662 0.708 0.638 0.810 0.843 0.638
犛犇r/犛犚r 0.078 - - 0.080 0.081 0.082 - - 0.082
(犚g-犚r)/(犚g+犚r) 0.150 0.153 0.150 0.154 0.155 0.302 0.316 0.315 0.302
(犛犇r-犛犇b)/(犛犇r+犛犇b) 0.523 - 0.344 0.533 0.552 0.587 - 0.340 0.587
(犛犇r-犛犇y)/(犛犇r+犛犇y) 0.632 0.631 0.628 0.632 0.632 0.747 0.742 0.734 0.747
(犛犇r-犛犚g)/(犛犇r+犛犚g) 0.674 0.550 0.182 0.687 0.687 0.849 0.762 0.301 0.849
(犛犇r-犛犚r)/(犛犇r+犛犚r) 0.052 0.055 0.059 0.117 0.117 0.173 0.179 0.184 0.173
  结合吸收特征参数相关系数来看,就会发现,在同类型的吸收特征参数中,相关系数高的参数,往往其回归模
型的决定系数也较高。综合来看,在3种吸收特征参数的多种回归方程中,以(犛犇r-犛犚g)/(犛犇r+犛犚g)的复合、
指数形式的估算方程的决定系数最高,为0.849,其次为犛犇r/犛犚g的S型回归形式的决定系数、犛犇r/犛犇b的S型
回归形式和犛犇r/犛犚g的乘幂形式,其决定系数都在0.8以上。
2.3 苜蓿鲜草产量高光谱估算模型的检验
从上述表中挑选决定系数在0.82以上的回归方程,拟作为估算模型进行检验。一个最优的估算模型,不仅
应该使估算值与实际值具有较高的相关系数(狉),同时还应具有较小的均方根误差(rootmeansquareerror,简称
RMSE)和平均相对误差(relativeerror,简称RE)[8]。因此,综合采用狉、RMSE和 RE这3个统计量对苜蓿鲜草
产量的高光谱估算模型进行检验,如表5所示。
总体来看,与以犛犇r/犛犇b、犛犇r/犛犚g、(犛犇r-犛犚g)/(犛犇r+犛犚g)为自变量的估算模型相比,分别以746,747,
748nm处一阶导数为自变量的复合、指数2种形式的估算模型,其相关系数普遍较高,且RMSE和RE也相对较
低。也就是说,经过一阶导数处理后的光谱变量,要比犛犇r/犛犇b、犛犇r/犛犚g、(犛犇r-犛犚g)/(犛犇r+犛犚g)等光谱吸
收参数等的估算精度要高。其中,以747nm处一阶导数为自变量的复合、指数2种形式的估算模型,其相关系
数最高(狉=0.852),且RMSE和RE也最低,分别为0.466kg/m2 和21.14%,由于这几个模型的估算精度最高,
且最能较准确地反演苜蓿鲜草产量,故以747nm处一阶导数为自变量的模型作为苜蓿鲜草产量估算模型。其
公式表示如下:
犉(狓=D747nm)=0.602×(3.018×10172)狓=0.602×e293.533狓
98第23卷第1期 草业学报2014年
表5 苜蓿鲜草产量高光谱估算模型检验结果(狀=120)
犜犪犫犾犲5 犜犺犲狏犪犾犻犱犪狋犻狅狀狉犲狊狌犾狋狅犳犺狔狆犲狉狊狆犲犮狋狉犪犾犲狊狋犻犿犪狋犻狅狀犿狅犱犲犾狊犳狅狉犳狉犲狊犺狔犻犲犾犱狅犳犪犾犳犪犾犳犪犻狀狋犲狊狋犱犪狋犪(狀=120)
参数(狓)Parameter(狓) 回归方程Regressionequation 相关系数狉 均方根误差RMSE(kg/m2) 相对误差RE(%)
D746nm 犢=0.413×(1.616×10155)狓 0.789 0.549 23.16
犢=0.413×e353.381狓 0.789 0.549 23.16
D747nm 犢=0.602×(3.018×10172)狓 0.852 0.466 21.14
犢=0.602×e293.533狓 0.852 0.466 21.14
D748nm 犢=0.792×(0.012×10106)狓 0.833 0.490 22.97
犢=0.792×e244.086狓 0.833 0.490 22.97
D750nm 犢=0.807×(2.557×10122)狓 0.824 0.504 22.94
犢=0.807×e281.854狓 0.824 0.504 22.94
犛犇r/犛犇b 犢=e(1.214-3.565/狓) 0.714 0.677 45.15
犛犇r/犛犚g 犢=e(1.239-3.252/狓) 0.757 0.637 42.74
(犛犇r-犛犚g)/(犛犇r+犛犚g) 犢=0.242×18.092狓 0.756 0.644 40.17
犢=0.242×e2.895狓 0.756 0.644 40.17
3 讨论与结论
本试验以10个苜蓿品种在不同时期的冠层高光谱数据为材料,并以多种高光谱参数作为变量进行分析,统
一分析苜蓿鲜草产量的高光谱特性,并建立了可供多个苜蓿品种统一使用的鲜草产量高光谱估算模型。
对光谱进行微分处理,是在实际应用中非常有效的技术。在一定程度上,微分处理可以去除土壤和噪声等多
种背景对植被光谱的部分影响,明显消除光谱基线平移及其旋转,从而提高了光谱对植被的敏感性。马勤建
等[15]发现,棉花冠层的一阶导数与叶面积指数及地上干物质量的相关系数要比反射光谱高得多。张凯等[9]发
现,在457,557和723nm处的一阶导数光谱,与甘南草地的地上鲜生物量的相关系数较高,通过了0.01极显著
性检验水平;其中以723nm处一阶导数为自变量的对数回归模型的估算精度较高。本研究中,在敏感波段的一
阶导数与苜蓿鲜草产量的相关系数及回归方程的决定系数,也比光谱反射率要高一些,说明导数光谱能较好地反
映植被的生长信息。
正确选取回归曲线形式能提高模型预测精度。许多学者以线性、对数、指数、乘幂和二次型、三次型多项式等
多种回归模型对植被生物量进行估算,结果显示非线性回归模型比线性回归模型更为精准,尤以乘幂模型突出。
黄春燕等[16]发现,在对棉花地上鲜生物量的高光谱估算模型中,以指数、乘幂和双曲线函数构建的估测模型,其
精度要比线性、对数要高,其中以指数模型的精度最高。柏军华等[17]应用5种回归模型对棉花地上鲜生物量进
行高光谱估算的研究中,发现指数、乘幂模型比线性、二次型、对数等函数模型的决定系数要高。谭昌伟等[18]采
用738nm处光谱一阶导数建立的指数回归模型,对夏玉米叶绿素密度估算精度最好。本研究以线性、对数、逆
模型、二次型、三次型多项式和复合、乘幂、S型、指数等9种回归形式,建立估算模型并加以验证,发现复合、指数
等非线性模型要优于线性模型。同时,在同类型的光谱参数中,线性模型决定系数高的参数,其二次型、三次型、
复合、乘幂、指数等非线性回归模型的决定系数通常也比较高,且通常高于线性模型,其中复合、指数最高。该结
果与前人研究结果大体相同。
在本研究的供试参数和回归模型中,以747nm处一阶导数为自变量的复合、指数2种形式的苜蓿产量估算
模型,其相关系数最高为狉=0.852,且均方根误差为0.466kg/m2,相对误差为21.14%,估算精度最高,因此该
模型可用于苜蓿鲜草产量的高光谱估算。
参考文献:
[1] 祁亚琴.基于高光谱数据提取棉花冠层特征信息的研究[D].石河子:石河子大学,2006.
09 ACTAPRATACULTURAESINICA(2014) Vol.23,No.1
[2] 马勤建.基于高光谱植被指数的棉花冠层结构参数的估算研究[D].石河子:石河子大学,2008.
[3] 钱育蓉,于炯,贾振红,等.新疆典型荒漠草地的高光谱特征提取和分析研究[J].草业学报,2013,22(1):157166.
[4] KawamuraK,WatanabeN,SakanoueS,犲狋犪犾.Wavebandselectionusingaphasedregressionwithabootstrapprocedurefor
estimatinglegumecontentinamixedsownpasture[J].GrasslandScience,2011,57(2):8193.
[5] ZhaoD,StarksPJ,BrownMA,犲狋犪犾.Assessmentofforagebiomassandqualityparametersofbermudagrassusingproximal
sensingofpasturecanopyreflectance[J].GrasslandScience,2007,53(1):3949.
[6] MutangaO,SkidmoreAK,PrinsHHT.PredictinginsitupasturequalityintheKrugerNationalPark,SouthAfricausing
continuumremovedabsorptionfeatures[J].RemoteSensingofEnvironment,2004,89(3):393408.
[7] PostCJ,DegloriaSD,CherneyJH,犲狋犪犾.Spectralmeasurementsofalfalfa/grassfieldsrelatedtoforagepropertiesandspe
ciescomposition[J].PlantNutrition,2007,30(11):17791789.
[8] 张凯,王润元,王小平,等.黄土高原春小麦地上鲜生物量高光谱遥感估算模型[J].生态学杂志,2009,28(6):11551161.
[9] 张凯,郭铌,王润元,等.甘南草地地上生物量的高光谱遥感估算研究[J].草业科学,2009,26(11):4450.
[10] 杨红丽,陈功,吴建付.施氮水平对多花黑麦草植株氮含量及反射光谱特征的影响[J].草业学报,2011,20(3):239244.
[11] 纳钦.苜蓿和缘毛雀麦高光谱与营养成分的相关性研究[D].呼和浩特:内蒙古农业大学,2010.
[12] 周宇庭,付刚,沈振西,等.藏北典型高寒草甸地上生物量的遥感估算模型[J].草业学报,2013,22(1):120129.
[13] 张艳楠,牛建明,张庆,等.植被指数在典型草原生物量遥感估测应用中的问题探讨[J].草业学报,2012,21(1):229
238.
[14] 浦瑞良,宫鹏.高光谱及其应用[M].北京:高等教育出版社,2000:5253.
[15] 马勤建,王登伟,黄春燕,等.棉花叶面积指数和地上干物质积累量的高光谱估算模型研究[J].棉花学报,2008,20(3):
217222.
[16] 黄春燕,王登伟,曹连莆,等.棉花地上鲜生物量的高光谱估算模型研究[J].农业工程学报,2007,23(3):131135.
[17] 柏军华,李少昆,王克如,等.基于近地高光谱棉花生物量遥感估算模型[J].作物学报,2007,33(2):311316.
[18] 谭昌伟,王纪华,黄文江,等.夏玉米叶片全氮、叶绿素及叶面积指数的光谱响应研究[J].西北植物学报,2004,24(6):
10411046.
犚犲狊犲犪狉犮犺狅狀狋犺犲犺狔狆犲狉狊狆犲犮狋狉犪犾狉犲犿狅狋犲狊犲狀狊犻狀犵犲狊狋犻犿犪狋犻狅狀犿狅犱犲犾狊
犳狅狉狋犺犲犳狉犲狊犺狔犻犲犾犱狅犳犪犾犳犪犾犳犪犵狉犪狊狊犾犪狀犱
LVXiaodong1,WANGJianguang1,SUNQizhong2,YAOGuiping3,GAOFengqin2
(1.ColegeofEcologyandEnvironment,InnerMongoliaAgriculturalUniversity,Hohhot010019,China;
2.GrasslandResearchInstituteofChineseAcademyofAgriculturalSciences,Hohhot010010,
China;3.ColegeofScience,InnerMongoliaAgricultural
University,Hohhot010019,China)
犃犫狊狋狉犪犮狋:Inthefield,hyperspectraldataandfreshyieldsindifferentgrowthperiodsof10alfalfa(犕犲犱犻犮犪犵狅
狊犪狋犻狏犪)varietieswerecolectedandthenthereflectance,firstderivativespectrumandspectralabsorptionfea
tureparameterswereusedasindependentvariables,tobuildunivariateregressionmodelsforestimatingthe
freshyieldofalfalfa.Thenonlinearregressionequationmodels,suchasquadratic,cubic,compound,power
andexponentialwerebetterthanlinearmodels.Amongthespectralparametersofthesametype,theparame
terswhichhadhighcoefficientofdeterminationsofthelinearmodel,alwayshadhighercoefficientsofdetermi
nationofquadratic,cubic,compound,powerandexponentialmodels.Amongtheseestimationmodels,the
compoundandexponentialestimationmodelsbasedonthefirstderivativevalueatwavelength747nmwasthe
mostaccurateforestimatingthefreshyieldofalfalfa;Thecorrelationcoefficient(狉)was0.852,therootmean
squareerror(RMSE)was0.466kg/m2,andtherelativeerror(RE)was21.14%.Themodelscanbeusedto
estimatethefreshyieldofthe10alfalfavarieties.
犓犲狔狑狅狉犱狊:alfalfa;variety;freshyield;hyperspectralremotesensing;regressionmodel
19第23卷第1期 草业学报2014年