《利用高光谱数据估测不同温度胁迫下的水稻籽粒中粗蛋白和直(共6页).doc》由会员分享,可在线阅读,更多相关《利用高光谱数据估测不同温度胁迫下的水稻籽粒中粗蛋白和直(共6页).doc(6页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上利用高光谱数据估测不同温度胁迫下的水稻籽粒中粗蛋白和直链淀粉含量基金项目:国家公益性行业科研专项项目(编号:GYHY);江苏省农业气象重点实验室开放基金(编号:JKLAM)。作者简介:谢晓金(1979-),女,安徽安庆人,讲师,博士,主要从事农业遥感、植物生理生态等研究。收稿日期:2012-03-15; 修回日期:谢晓金1,2,李秉柏1,朱红霞1(1.江苏省农业科学院农业经济与信息研究所,江苏 南京;2.南京信息工程大学江苏省农业气象重点实验室,江苏 南京 )摘 要:水稻籽粒粗蛋白和直链淀粉含量是评价稻米品质的两个主要指标。本文以不同温度胁迫下的水稻田间试验为基础,分
2、析水稻成熟籽粒粗蛋白和直链淀粉含量与不同生育期冠层反射光谱的相关性。结果显示,水稻籽粒的粗蛋白和直链淀粉含量与不同生育期冠层原始光谱和一阶导数在某些波段达到极显著或显著相关水平。进一步分析比值植被指数、差值植被指数、归一化植被指数和红边参数等光谱参数与成熟籽粒粗蛋白和直链淀粉含量的相关性,用回归分析方法对相关拟合较好的参数进行筛选,建立了水稻成熟籽粒粗蛋白含量( GCPC)和直链淀粉含量( GAC)监测模型,运用独立观测数据对模型进行检验,其预测值与实测值的精确度为0.3930.683,准确度为0.7080.923,RMSE值为8.706%11.296%。表明模型预测值与实测值之间具有较高的符
3、合度,对水稻籽粒粗蛋白质和直链淀粉含量具有较好的预测性。关键词: 水稻籽粒;粗蛋白质含量;直链淀粉含量;反射光谱;光谱参数冠层反射光谱能够反映作物群体面源信息,通过对冠层光谱信息的解析,可以有效提取和突出作物目标的实时信息。多年来,反射光谱数据已成为估测叶面积指数、地上生物量、叶绿素含量和氮素营养等农作物长势监测和遥感估产的重要手段1-2。粗蛋白和直链淀粉含量是谷类作物营养品质的两个重要指标,近年来,国内外许多研究表明,利用反射光谱来预测谷类籽粒的粗蛋白和直链淀粉含量也是可行的。Hansen等报道了利用冠层光谱反射率和偏最小二乘法预测小麦籽粒蛋白质含量3。李映雪等指出, 灌浆中期冠层光谱指数
4、RVI (1 220, 710)能预测不同小麦品种成熟期籽粒蛋白质含量的变化4。刘芸等利用原始光谱反射率和一阶导数光谱来构建植被指数来预测籽粒粗蛋白、粗淀粉和直链淀粉含量,其决定系数达0.7以上5。但以上研究主要偏重于氮素胁迫下的谷类籽粒品质的光谱反演,而定量反演高温胁迫条件下水稻籽粒品质的研究未见报道。目前,随着全球气候变化的日益加剧,夏季极端高温和持续高温频繁出现,我国长江流域水稻发生热害的频率与受害程度随之加大,对水稻品质损害的程度也越来越严重,而利用高光谱遥感技术进行高温热害长势和品质监测的研究必是未来农业遥感重心之一6,7。本研究以不同年份、不同品种类型在不同温度胁迫下的水稻试验为基
5、础,测定不同生育期水稻冠层高光谱反射率以及成熟籽粒的粗蛋白和直链淀粉含量,讨论水稻籽粒的粗蛋白和直链淀粉含量与冠层反射光谱参数的定量关系,为快速、低成本地检测稻米品质提供依据,同时也为用高光谱遥感方法在水稻高温热害中的品质监测应用提供重要的理论依据。1材料与方法1.1试验设计试验1:试验于2007、2008年在江苏省农业科学院试验站进行,供试水稻品种为扬稻6号(常规籼稻,全生育期约138d)和南粳43(常规粳稻,全生育期约160d)。5月15日育秧,6月18日移栽,扬稻6号于8月18日开始抽穗,南粳43于9月1日开始抽穗。采用浙江余姚生产的远红外加热灯管(长1.5m,额定功率1000W)和自动
6、控温系统,自行设计高1.7m,长、宽各1.5m的柱型钢管支架,安置在水稻试验田中,远红外加热管离地1.5m高。在温控和通风系统的调控下,使加热区(离地面0.81.2m高的空间)的温度变幅控制在0.2左右。在两个水稻品种处于孕穗期时,随机选取10箱水稻进行高温处理,另外10箱为常温对照(CK),加热设施提前30min通电预热,以确保09:00-14:00时段内穗层温度为35、38与41,处理时间为3d,每天高温处理5h后停止增温。在处理时段内,常温对照的气温为3031,相对湿度为65%,高温处理的相对湿度为68%。试验2:在以上两个水稻品种处于孕穗期时,随机选取5箱水稻移入RXZ型人工气候光照培
7、养箱中,进行35、38与41 3个温度梯度处理,处理时间为3d,每天处理5h(09:00-14:00),每天高温处理结束之后放置于网室内。以自然环境条件下温度为对照。试验资料用于检验预测模型。1.2试验测定1.2.1光谱测定于水稻开花期、灌浆期及蜡熟期采用美国Analytical Spectral Device(ASD)公司生产的FieldSpec Pro FR2500型背挂式野外高光谱辐射仪测定冠层光谱反射率,光谱仪的波段范围为3502 500 nm,其中,3501 000 nm光谱采样间隔为1.4 nm,光谱分辨率为3nm; 1 0002 500 nm光谱采样间隔为2 nm,光谱分辨率为1
8、0 nm。冠层光谱测定选择在天气晴朗、无风或风速很小时进行,时间范围为10: 0014: 00,传感器探头垂直向下,距冠层顶垂直高度约1. 0m,光谱仪视场角为25,地面视场范围直径为0.44m。每个小区记录10个采样光谱,取平均值作为该处理的高光谱反射值。测量过程中及时进行标准白板校正。用ASD ViewSpec Programs 软件读取反射光谱原始数据。1.2.2 粗蛋白和直链淀粉测定将光谱测量区域的水稻收获后碾磨去壳,得大米。大米直链淀粉的提取参照文献8方法。用凯氏定氮法测全氮含量,根据一般大米蛋白质含约16.80%氮的原理,将所得的全氮含量乘以换算因子K=5.95,推算为粗蛋白含量。
9、1.3数据分析参照已有的植物光谱参数计算方法9-12,基于所有波段的反射率、一阶导数构建两波段组合的比值、差值和归一化植被指数以及红边幅值等特征光谱参数,以试验1为基础,对水稻冠层各项光谱参数与成熟籽粒粗蛋白质和直链淀粉含量进行相关分析,选择与籽粒粗蛋白质和直链淀粉含量显著相关的波段及植被指数进行回归拟合,并通过回归分析确定监测水稻籽粒蛋白质和直链淀粉含量的最佳光谱参数,建立预测模型。利用试验2的观测数据对以上模型进行测试和检验。采用国际上常用的根均方差法(RMSE , root mean square error ) 对模型进行检验。RMSE以模拟值与观测值的一致性来反映模型的预测能力,RM
10、SE值小于10 %,表明预测性强;RMSE值大于30 %,则表明预测性差。并引用Massart 等13提出的精确度(Precision)和准确度(Accuracy)概念对模型的预测精度进行评价,精确度即实测值与预测值之间的决定系数,准确度即实测值与预测值拟合方程的斜率。2 结果与分析2.1高温胁迫下籽粒粗蛋白和直链淀粉含量的变化图1 高温胁迫下水稻籽粒粗蛋白和直链淀粉含量的变化图1展示出不同温度胁迫条件下扬稻6号和南粳43的粗蛋白和直链淀粉含量的变化特点。可以看出,与对照相比,两个水稻品种的籽粒粗蛋白和直链淀粉含量均达到显著或极显著水平(p0.05或P0.01),并且随胁迫温度的增加,两个品种
11、的籽粒粗蛋白含量呈下降趋势,而直链淀粉含量呈上升趋势。2.2冠层反射光谱与水稻籽粒粗蛋白和直链淀粉含量的关系2.2.1冠层单波段反射率与籽粒粗蛋白和直链淀粉含量的相关性图2不同生育期冠层光谱反射率与水稻籽粒粗蛋白和直链淀粉含量的相关系数A:开花期B:灌浆期C:蜡熟期从图2可以看出,水稻籽粒的粗蛋白和直链淀粉含量与3个生育期冠层光谱反射率的相关性具有相反的趋势,其中粗蛋白与开花期和灌浆期冠层光谱反射率在整个光谱区域都达从图2可以看出,水稻籽粒的粗蛋白和直链淀粉含量与3个生育期冠层光谱反射率的相关性具有相反的趋势,其中粗蛋白与开花期和灌浆期冠层光谱反射率在整个光谱区域都达到正相关,直链淀粉含量与开
12、花期冠层光谱反射率仅在红光(642-691nm)范围内呈正相关,在其他光谱范围内均呈负相关,而与灌浆期冠层光谱反射率在整个光谱区域都达到正相关。对于蜡熟期,粗蛋白在350-532 nm、575-695 nm以及1 760 nm以上光谱区域呈负相关,在其他光谱范围内均呈正相关。直链淀粉含量在350-706 nm光谱区域呈正相关,而706 nm以上可见光和红外区域光谱为负相关。另外,水稻籽粒的粗蛋白和直链淀粉含量与不同生育期叶片冠层光谱一阶导数在某些波段(如376、416、628和715 nm等)也达到了显著或极显著相关水平(p0.05或p0.01)。这表明,可以通过冠层原始高光谱反射率和一阶导数
13、来估测籽粒的粗蛋白和直链淀粉的相对含量。本试验还显示,籽粒的粗蛋白与直链淀粉含量之间呈显著负相关,相关系数达-0.356,表明水稻籽粒粗蛋白含量越高,其直链淀粉含量相对就越低。2.2.2各类光谱变量与籽粒粗蛋白和直链淀粉含量的相关性将水稻冠层原始光谱反射率、一阶导数以及组合成的光谱变量与粗蛋白和直链淀粉含量进行相关性分析,并将相关系数较大的光谱变量列于表1,部分估算方程列于表2。结果显示:3个生育期冠层光谱变量对粗蛋白和直链淀粉含量预测性都较好,均达到显著或极显著水平(p0.05或p0.01),但开花期的冠层光谱变量对粗蛋白和直链淀粉含量的预测性优于灌浆期,而灌浆期的预测性又优于蜡熟期,光谱变
14、量中的差值植被指数DVI(810,450)、DVI(810,680)和670-755nm面积可以同时预测籽粒的粗蛋白和直链淀粉含量。表1不同生育期冠层光谱参数与水稻籽粒粗蛋白和直链淀粉含量的相关性(n=40)开花期 灌浆期 蜡熟期光谱参数相关系数光谱参数相关系数光谱参数相关系数粗蛋白DVI(810,450)0.812*DVI(810,450)0.753*DVI(810,450)0.567*DVI(810,560)0.722*DVI(810,560)0.656*DVI(810,560)0.459*DVI(810,680)0.771*DVI(810,680)0.746*DVI(810,680)0.
15、398*PVI(810,680)0.848*DVI(560,450)0.383*PVI(810,680)0.380*Area 670-7550.850*DVI(560,680)0.466*PVI0.397*R7430.842*Area 670-755 0.757*Area 670-7550.371*Dr0.661*Dr0.644*r 8970.791*FD7220.722*FD7230.624*FD715-0.471*直链淀粉DVI(810,450)-0.853*DVI(810,450)-0.771*DVI(810,450)-0.595*DVI(810,560)-0.783*DVI(560,6
16、80)-0.485*DVI(810,560)-0.490*DVI(810,680)-0.832*DVI(810,560)-0.647*DVI(810,680)-0.457*PVI-0.815*PVI(810,680)-0.683*Dr-0.361*Dr0.732*PVI-0.686*Area 670-755-0.448*Area 670-755 -0.834*Area 550-5820.468*R1077-0.671*R800-0.854*Area 670-755-0.749*R1106-0.652*FD722-0.804*FD723-0.663*Fd715-0.521*和*分别表示0.05和
17、0.01显著水平。表2不同生育期冠层光谱参数与水稻籽粒粗蛋白和直链淀粉含量的回归方程(n=40)开花期 灌浆期 蜡熟期 光谱参数回归方程光谱参数回归方程光谱参数回归方程 粗蛋白DVI(810,450)y = 0.010x + 4.386PVIy =0.144x0.750DVI(810,450)y = 1.152e0.011xPVI(810,680)y = 1.021Ln(x) -1.963R901y =1.467Ln(x) +2.902R1076y = 1.077Ln(x) + 2.812R880y = 1.002Ln(x) +2.191FD723y =0.951x0.711FD997y =
18、1.55x0.694直链淀粉DVI(810,450)y =74.586x-0.334DVI(810,450)y = 31.171e-0.011xDVI(810,450)y = 24.170e-0.006xPVI(810,680)y = -7.022Ln(x) + 44.601DVI(810,680)y = 29.700e-0.010xFD715y = 21.882e-0.223xR750y =15.797x-0.334R929y = 26.9701e-0.415xR897y = 13.483 x-0.3522.3籽粒粗蛋白与直链淀粉含量的监测模型及检验基于单波段反射率、一阶导数以及表1所列的光
19、谱参数,对籽粒粗蛋白和直链淀粉含量进行回归分析,筛选DVI(810,450) 作为最佳光谱参数,并且开花期的DVI(810,450)对籽粒粗蛋白和直链淀粉含量的预测优于灌浆期和蜡熟期。图3展示了开花期DVI(810,450)与成熟籽粒粗蛋白和直链淀粉含量的关系,并得出基于开花期冠层反射光谱的籽粒粗蛋白含量(GCPC)和直链淀粉含量(GAC)的预测方程:GCPC=0.010 DVI(810,450) + 4.386 R2=0.662 n=40 (1)GAC=74.586 DVI(810,450)-0.334 R2=0.754 n=40 (2)图3 水稻籽粒粗蛋白和直链淀粉含量与开花期DVI(81
20、0,450)的关系利用试验2开花到蜡熟期的数据对水稻成熟籽粒粗蛋白和直链淀粉含量预测方程进行测试和检验(见表3),结果显示,预测值与实际值的精确度为0.3930.683,均达到极显著水平;准确度为0.7080.923,RMSE值为8.706%11.296%。说明粗蛋白和直链淀粉含量预测模型对不同时期水稻的估测效果具有一定的可靠性。表3籽粒粗蛋白和直链淀粉含量预测模型的预测表现生育时期开花期灌浆期蜡熟期粗蛋白精确度0.6780.6830.547准确度0.7080.7080.743均方根差9.23610.1079.576直链淀粉精确度0.3930.4780.640准确度0.7610.9230.76
21、3均方根差8.70611.46311.2963结论与讨论3.1结论与常温对照相比,两个水稻品种的籽粒粗蛋白和直链淀粉含量均达到显著或极显著水平,并且随胁迫温度的增加,两个水稻品种的籽粒粗蛋白含量呈上升趋势,而直链淀粉含量呈下降趋势。将水稻冠层原始光谱反射率、一阶导数以及组合成的光谱变量与粗蛋白与直链淀粉含量进行相关性分析,结果表明:三个生育期冠层光谱变量对粗蛋白与直链淀粉含量预测性都较好,但开花期的预测性优于灌浆期与蜡熟期,其中光谱变量中的DVI(810,450)、DVI(810,680)和670-755nm面积三个参数可以同时预测籽粒的粗蛋白与直链淀粉含量。3.2讨论前人研究表明,在氮素和干
22、旱胁迫下,谷类籽粒的粗蛋白和直链淀粉含量都有显著变化14,15,本研究发现,与常温对照相比,高温胁迫下水稻籽粒的粗蛋白和直链淀粉含量均达到显著或极显著水平。冠层反射光谱能够反映作物群体面源信息,在一定程度上能反映籽粒的品质信息,因此利用反射光谱预测籽粒的粗蛋白和直链淀粉含量具有一定的可行性。本研究发现,基于开花期的冠层光谱参数对籽粒粗蛋白和直链淀粉含量的预测均优于灌浆期和蜡熟期,原因可能由于水稻籽粒营养成分的形成大部分来自抽穗前茎、叶营养贮存的转移,即主要取决于灌浆前期,而蜡熟期水稻叶、穗大部分变黄,叶面积与叶片叶绿素含量明显降低,这与粗蛋白和淀粉正相关的叶片叶绿素对整个冠层光谱的贡献相对减弱
23、,因此,相比之下,灌浆期和蜡熟期冠层光谱预测水稻籽粒的粗蛋白和淀粉的精度降低。本研究基于开花期的冠层反射光谱参数DVI(810,450)构建了成熟期籽粒粗蛋白和直链淀粉含量的预测模型,通过检验表明,其估测效果较好,因此该模型具有一定得的准确性,此模型为高光谱遥感方法在水稻高温热害品质监测中的应用提供了方法和依据。另外,由于水稻籽粒粗蛋白质和直链淀粉含量的形成是一个动态过程,因此要利用遥感手段进行水稻籽粒粗蛋白质和直链淀粉含量的精确监测预报,还必须深入探明稻穗的特征光谱及不同叶位叶片和茎秆对冠层反射光谱的影响及其与籽粒蛋白质形成的机理关系。参考文献1冯伟,朱艳,田永超,等.基于高光谱遥感的小麦叶
24、片氮积累量J.生态学报,2008,28(1):23-322Huang Z, Turner B J, Dury S J, et al. Estimating foliage nitrogen concentration from HYMAP data usingcontinuum removal analysis. Remote Sens Environ,2004,93:18-29.3Hansen P M, Jorgensen J R , Thomsen A. Predicting grain yield and protein content in winter wheat and sprin
25、g barley using repeated canopy reflectance measurements and partial least squares regressionJ. J Agri Sci, 2002,139(3):307-318.4李映雪,朱艳,田永超,等.小麦冠层反射光谱与籽粒蛋白质含量及相关品质指标的定量关系J.中国农业科学,2005,38(7):1332-1338.5刘芸,唐延林,黄敬峰,等.利用高光谱法估测稻穗稻谷的粗蛋白质和粗淀粉含量J.中国农业科学,2008,41(9):2617-2623.6雷东阳,陈立云,李稳香,等.杂交水稻抽穗扬花期高温对结实率及相关生
26、理特陛的影响J杂交水稻,2006,21(3):68-77程高峰,张佳华,李秉柏,等.不同温度处理下水稻高光谱及红边特征分析J.江苏农业学报,2008,24(5):573-5808洪 雁,顾正彪,刘晓欣.直链淀粉和支链淀粉纯品的提取及其鉴定J.食品工业科技, 2004,25(4):86-88.9 Rouse J W, Haas R H, Schell J A. Monitoring the vernal advancement of retrogradation of natural vegetationJ. NASA/GSFC,Type III, Final Report, Greenbelt
27、, MD, USA, 1974.10Richardson A J, Wiegand C L. Distinguishing vegetation from soil background information. Photogrammetric Engineering and Remote Sensing, 1977,43(12): 1541-1552.11Horler D N H, Dockray M, Barber J. The red edge of plant reflectance J. International Journal of Remote Sensing.1983,4:2
28、73-288.12Jordan C F. Derivation of leaf area index from quality of light on the forest floor J.Ecology,1969,50:663-666. Gupta R K,Vijayan D, Prasad T S.Comparative anaylysis of red edge hyperspectral indices. Advance in Space Research, 2003,32:2217-2222.13Massart D L, Buydens L, Peeters A.Knowledge representation for the selection of methods in analytical chemistry J. Chemometrics and Intelligent Laboratory Systems,1988,5: 73-79.14冯 伟,姚 霞,田永超,等.J.作物学报,2007,33(12):1935-194.15赵 辉,荆 奇,戴廷波,等.J.作物学报,2007,33(12):2021-2027.专心-专注-专业