白粉病胁迫下小麦冠层叶绿素密度的高光谱估测.pdf
《白粉病胁迫下小麦冠层叶绿素密度的高光谱估测.pdf》由会员分享,可在线阅读,更多相关《白粉病胁迫下小麦冠层叶绿素密度的高光谱估测.pdf(10页珍藏版)》请在咨信网上搜索。
1、第 29 卷 第 13 期 农 业 工 程 学 报 Vol.29 No.13 114 2013 年 7 月 Transactions of the Chinese Society of Agricultural Engineering Jul.2013 白粉病胁迫下小麦冠层叶绿素密度的高光谱估测 冯 伟,王晓宇,宋 晓,贺 利,王晨阳,郭天财(河南农业大学国家小麦工程技术研究中心,郑州 450002)摘 要:为了明确病害胁迫下作物生长特征及其危害程度,基于大田小区和盆栽小麦白粉病接种试验,在病害胁迫下不同生育时期测定群体光谱及叶绿素密度。综合分析群体光谱反射率、一阶微分及传统光谱特征参数与冠层
2、叶绿素密度间关系,建立病害叶绿素密度估算模型并检验。结果表明,随病情指数增加,叶绿素含量下降,不同感性品种均如此,对白粉病易感品种的危害较重。病害冠层叶绿素密度与红光 600630 nm 和红边 690718 nm的反射率及红边长波段(718756 nm)的一阶微分间相关性最显著。在传统植被指数中,以 SDr/SDb 和 VOG3为变量的估测模型拟合精度较高,决定系数 R2分别为 0.752 和 0.723,模型检验相对误差(RE)最小,RE 分别为18.0%和 18.6%。利用红边区域(680760 nm)波段差异特性,选取 680、718 和 756 nm 波段新建红边角度指数(REAI)
3、,较传统植被指数的模型拟合精度更高,归一化角度指数 NDAI(,)和比值角度指数 RAI(,)的 R2分别为 0.783 和 0.776,模型检验误差更小,RE 分别为 16.8%和 17.5%。因此,NDAI(,)是估测病害小麦冠层叶绿素密度的可靠指标,对利用该模型监测小麦光合潜力和病害影响评价具有积极意义。关键词:光谱分析,叶绿素,模型,估测,白粉病胁迫,冬小麦 doi:10.3969/j.issn.1002-6819.2013.13.016 中图分类号:Q657.3 文献标志码:A 文章编号:1002-6819(2013)-13-0114-10 冯 伟,王晓宇,宋 晓,等.白粉病胁迫下小
4、麦冠层叶绿素密度的高光谱估测J.农业工程学报,2013,29(13):114123.Feng Wei,Wang Xiaoyu,Song Xiao,et al.Hyperspectral estimation of canopy chlorophyll density in winter wheat under stress of powdery mildewJ.Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE),2013,29(13):114123.(in Chine
5、se with English abstract)0 引 言 光合作用中吸收光能的植物色素有叶绿素和类胡萝卜素,其中叶绿素是主要的吸收光能物质,其含量高低与光合能力、发育时期和营养水平关系密切,可敏感指示环境胁迫、光合潜力和阶段发育1-2。小麦白粉病是世界性病害,中国在 20 世纪 70 年代以前主要在西南及山东沿海局部发病严重,随后发生范围不断扩大,危害程度加重,导致叶片退绿发黄,叶绿素含量降低,严重影响小麦光合作用3-4。因此,通过快速光谱测试实时无损估算胁迫下叶绿素含量,对作物病害及群体光合能力进行有效评价,在作物精确诊断和高效管理中具有重要意义。绿色植物典型的光谱特征,在可见光范围内主
6、 收稿日期:2012-12-15 修订日期:2013-05-31 基金项目:国家自然科学基金(30900867)和现代农业(小麦)产业技术体系(MATS)专项资助。作者简介:冯 伟(1976),男,博士,副研究员,主要从事小麦生长监测与遥感应用研究。郑州 国家小麦工程技术研究中心,450002。Email: 通信作者:郭天财(1953),男,教授,博导,主要从事小麦生态与高产栽培技术研究郑州。郑州 国家小麦工程技术研究中心,450002。Email: 要受色素影响,而在近红外区域主要受叶片内部结构、生物量和含水量等影响。因此,可以利用植被光谱信息特征来估算叶绿素含量5-6。Horler 等7研
7、究了植被光谱与叶绿素浓度的关系,提出反射光谱红边位置可指示植被叶绿素浓度变化。由于红边一阶微分数据的离散性,传统算法的红边位置对叶绿素含量反演产生较大偏差,Miller 等和 Cho 等8-9分别利用倒高斯光学模型和线性外推法修正了红边位置算法,对叶绿素状况的估算取得了较好效果。Blackburn10-11发现,冠层和叶片尺度上色素浓度与反射光谱的一阶和二阶导数极显著相关,并利用叶绿素特征吸收波段 470、635 和 680 nm 构造了适宜估算叶绿素状况的光谱指数,但 Sims 等12进一步研究表明,550、700 nm 较 680 nm 反射率对叶绿素更敏感。孟庆野等13在 PROSPEC
8、T+SAIL 辐射传输物理模型基础上对红边特征进行相关分析,引入红边、红谷和绿峰 3 个敏感波段,提出了改进转换叶绿素吸收反射指数(modified transformed chlorophyll absorption ratio index,MTCARI)。以上研究均基于植被健康生长条件下利用光谱分析技术估算叶绿素含量,但当作物遭受病害胁迫后,植株生物量下降,叶绿素和水分减少,叶片结构破农业信息与电气技术第 13 期 冯 伟等:白粉病胁迫下小麦冠层叶绿素密度的高光谱估测 115 坏,在光谱特征上常表现为可见光和短波红外波段光谱反射率增加,近红外波段减小,这与绿色植被典型光谱特征存在很大差异1
9、4-15。前人研究确立的健康植被条件下其叶绿素遥感估测指标及模型很难适用于病虫危害后植被叶绿素状况的遥感监测与评价。近年来作物病虫害光谱分析与监测评价相结合逐渐受到众多学者的广泛关注,目前有关病虫危害后植被理化参数遥感监测的研究报道较少,而白粉病胁迫下小麦冠层叶绿素密度的光谱估测尚未见报道。本文开展白粉病胁迫下小麦冠层叶绿素密度高光谱估测研究,建立小麦冠层叶绿素密度的反演模型,以期用简便、无损的遥感技术估测病害小麦冠层叶绿素密度,为大面积遥感监测小麦病害、光合能力及估产提供理论依据和技术支持。1 材料与方法 1.1 试验设计 试验 1(EXP.1)于 2009-2010 年在河南农业大学科教示
10、范园区大田进行。供试材料为对白粉病表现易感的品种偃展 4110 和豫麦 34,中感品种周麦18 和郑麦 366。前茬为玉米,土壤为砂壤土,播前030 cm 土壤有机质质量分数 8.3 g/kg,全氮0.87 g/kg,速效磷24.44 mg/kg,速效钾 124.32 mg/kg,施肥量为 270 kg/hm2纯氮、120 kg/hm2 P2O5和90 kg/hm2 K2O。试验设计为 3 次重复,小区面积为20.3 m2。10 月 15 日播种,三叶期定苗,各小区栽培管理均一致,其他同一般高产麦田。同时在病圃田培养偃展 4110 等易感材料,作为试验接种的病菌来源。依据接种频次设定重(1 d
11、 接种一次,I3)、中(2 d 接种一次,I2)和轻(3 d 接种一次,I1)3个接种等级处理,另设不接种的对照区(I0)。小麦拔节期开始接种白粉病菌,一直持续至最低级别处理开始表现明显发病症状时停止人工接种。小麦抽穗期开始染病,冠层光谱测试和采样时期为花后3 d、花后 15 d 和花后 25 d。试验 2(EXP.2)于 20092010 年在河南农业大学科教示范园区温棚进行。试验材料为易感品种偃展 4110 和豫麦 49,中感品种矮抗 58 和郑麦 366。土壤为过筛后砂壤土,每盆装土 15 kg,盆面直径25 cm,花盆埋于陆地,盆面高出地面 2 cm,所用肥料同大田,施用量为 2 倍大
12、田。10 月 18 日播种,三叶期定苗。在拔节后期进行白粉病接种,依据接种量大小和频次设定重、中和轻 3 个接种等级,另设不接种的对照区。小麦孕穗期开始染病,冠层光谱测试和采样时期为抽穗期和花后 10 d,测定时分单盆、2 盆并和 4 盆并进行,3 种方式均测定 2 组数据。试验 3(EXP.3)于 20102011 年在河南农业大学国家小麦工程技术研究中心病圃田进行。供试材料为易感品种偃展 4110 和西农 979,中感品种矮抗 58。土壤为砂壤土,播前 030 cm 土壤有机质质量分数 16.3 g/kg,全氮 1.22 g/kg,速效磷35.65 mg/kg,速效钾 114.68 mg/
13、kg,施肥量为270 kg/hm2纯氮、120 kg/hm2 P2O5和90 kg/hm2 K2O。依据接种量大小和频次设置 4 个接种梯度处理,2次重复,小区面积为 8.1 m2。10 月 15 日播种,三叶期定苗,各小区栽培管理均一致,其他同一般高产麦田。小麦孕穗期开始染病,冠层光谱测试和采样时期为抽穗期和花后 15 d。1.2 测定项目与方法 1.2.1 光谱数据获取 小麦冠层光谱测定采用 FieldSpec HandHeld 手持便携式光谱分析仪(美国 ASD 公司生产),波段值为 3251 075 nm,光谱采样间隔为 1.6 nm,光谱分辨率为 3.5 nm,探头视场角 25。冠层
14、光谱测定选择在天气晴朗、无风或风速很小时进行,北京时间 10:0014:00,传感器探头垂直向下,距冠层顶垂直高度约 1.0 m,地面视场范围直径为0.44 m。测量过程中及时进行标准白板校正,以 10 条光谱为一采样光谱,每个采样点记录 10 个光谱,以其平均值作为该点的光谱反射率。1.2.2 病情指数调查 在各小区随机调查 30 株小麦发病程度,严重度具体分 9 个等级:即 0、1%、10%、20%、30%、45%、60%、80%和 100%,分别记录各严重度的小麦叶片数,然后利用式(1)算出病情指数16:()(%)100 xfDInf=(1)式中,x 为各梯度的级值,n 为最高梯度值为
15、9,f为各梯度的叶片数。1.2.3 冠层叶绿素密度测定 与光谱测量同步,每小区取有代表性单茎 12个,分 2 份分别测定叶绿素含量和生物量。叶绿素含量测定时,取叶片中部 0.1 g 剪成细丝,用无水乙醇和丙酮(v/v=1/1)混合液暗处提取,用紫外可见分光光度计测定 470、645 和 663 nm 处的吸光度,然后参考 Lichtenthaler 方法17计算叶绿素含量。将叶片质量乘以对应叶片叶绿素含量,即得单位土地面积绿叶片的叶绿素总量(干物质量),作为冠层叶片叶绿素密度(kg/hm2)。1.3 数据分析与利用 由于光谱仪采集的是离散数据,因此采用式(2)计算光谱反射率的一阶微分值18。农
16、业工程学报 2013 年 116 11()()()()2iiiidRRRRd+=(2)式中,R(i)为i nm处的光谱反射率,R(i)是i处的一阶微分,i是第i个通道的波长。在一阶微分中提取了红边特征,其在680760 nm之间反射率增长很快,为健康植被的显著光谱特征。参考相关文献14,18选用了不同类型光谱参数,在Matlab语言环境下编程实现了光谱指数与叶绿素密度间相关系数运算,选择8个与叶绿素含量最显著相关的光谱特征参数列于表。以试验1和试验2资料为基础,对常见特征光谱指数与冠层叶绿素密度间进行相关回归分析,优选和构建敏感光谱参数,建立小麦冠层叶绿素密度监测模型。利用试验3独立资料对所筛
17、选建立的模型进行测试和检验,采用预测精度(R2)、均方根误差(RMSE,式(3)和相对误差(RE,式(4)等统计指标进行评价。211()niiiRMSEPOn=(3)()211100%niiiiOPREnO=(4)式中,Pi和Oi分别为估测值和观察值,n为样本数。表 1 本研究采用的高光谱参数算法及出处 Table 1 Summary of different hyperspectral parameters used in this study 光谱参数 Spectral parameter 定义与算法 Definition and algorithm 参考文献 Reference RVI(
18、810,560)810 nm 与 560 nm 处反射率比值 Ratio of reflectance at 810 nm to 560 nm R810/R560 Rouse et al.(1974)19 GM1 750 nm 与 550 nm 处反射率比值 Ratio of reflectance at 750 nm to 550 nm R750/R550 Gitelson and Merzlyak(1994)20 VOG3 红边光谱比值植被指数 Red edge reflectance ratio indices (R734-R747)/(R715+R720)Zarco-Tejada et
19、 al.(2001)21 mSR705 修正 SR705 简单比值指数 Simple ratio indices modified on SR705(R750-R445)/(R705-R445)Sims&Gamon(2002)12 PSSRb 特定色素简单比值法 Pigment specific simple ratio R800/R635 Blackburn(1998)10 ZM 750 nm 与 710 nm 处反射率比值 Ratio of reflectance at 750 nm to 710 nm R750/R710 Gupta et al.(2003)22 SDr/SDb 红蓝边面
20、积比 Radio of the red-edge integral areas and the blue-edge integral areas 7555301111680490()/2()/2iiiiiiRRRR+=Feng et al.(2008)23 REPIG 基于倒高斯模型的红边位置 Red edge position obtained from inverted gaussian fitting technique 2002()()()exp2ssRRRR=REPIG=i+Miller et al.(1990)8 注:R 和 分别表示原始光谱反射率和波长,Rs为近红外区域肩反射率,
21、R0为近红外区域低谷反射率,为高斯模型变异参数。Notes:R and denote spectral reflectance and wavelength respectively.Rs is the“shoulder”spectral reflectance within the red edge(680760 nm),R0 is the minimum spectral reflectance at wavelength 0corresponding to the chlorophyll absorption well,and is the Gaussian function devia
22、tion parameter.2 结果与分析 2.1 不同接种等级下小麦叶绿素含量和病情指数的的变化 以试验1的易感品种偃展4110和中感品种周麦18为例,说明不同接种处理对叶绿素含量和病情指数的影响。从表2可以看出,在病菌感染的灌浆期随灌浆进程小麦叶绿素含量呈下降趋势,2品种一致。随接种等级增加,植株发病严重,病情指数升高,叶绿素含量下降;且随发病日期延长,病情指数增加,病害对叶绿素的影响加重。品种间比较,易感品种偃展4110发病速度更快,叶绿素含量损失更多。就单个品种单个时期而言,病情指数与叶绿素含量呈显著负相关(r0.65*)。在红边的近红外端(730750 nm)相关性很弱(|r|0.
23、2)。在750950 nm近红外波段内则呈显著正相关,但相关系数相对较小(r0.4*),这是由于近红外波段主要受冠层及叶片组织结构影响,而对叶绿素密度的变化不敏感。由此可知,可见红光及红边波段对小麦叶绿素密度的变化较为敏感,通过对此区域光谱信息的挖掘将有助于进一步解析白粉病侵染后小麦叶绿素密度与对应光谱之间的定量关系。图 2 小麦叶片叶绿素密度与原始光谱反射率的相关性 Fig.2 Correlation patterns of canopy chlorophyll density to spectral reflectance of wheat 2.4 病害小麦叶绿素密度与常见光谱参数的定量关
24、系 以前人提出的不同类型光谱参数(比值、差值、归一化、红边特征及微分等)为基础,综合试验1和试验2的不同生育时期小麦冠层叶绿素密度与光农业工程学报 2013 年 118 谱参数进行指数及线性模型估算(n=183),将其中表现较好的8个光谱参数、方程拟合的决定系数及标准误差列于表3,并选用2个光谱参数为代表,将回归方程作散点图直观展示方程拟合效果(图3)。结果显示,光谱参数RVI(810,560)、GM1和PSSRb的线性回归方程的拟合精度(决定系数)优于指数模型,其他光谱参数的指数模型拟合精度均高于线性模型;从方程拟合标准误差来看,所有优选参数的方程拟合RMSE均为指数模型低于线性模型。综合而
25、言,选用误差较小的指数模型可以较好地对叶绿素密度进行建模估算,指数模型中特征光谱参数VOG3和SDr/SDb对叶绿素密度拟合效果较好,决定系数分别为0.723(P0.001)和0.752(P0.001),拟合误差分别为1.311和1.240 kg/hm2。表 3 小麦叶片叶绿素密度与光谱参数的定量关系 Table 3 Quantitative relationships of canopy chlorophyll density(y)to different spectral indices(x)in wheat 指数模型 Index model 线性模型 Linear model 光谱参数
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 白粉病 胁迫 小麦 叶绿素 密度 光谱 估测
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【二***】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【二***】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。