亢菊俠 趙家奇 吳雅茹
摘要:為構(gòu)建基于高光譜參數(shù)的設(shè)施黃瓜霜霉病發(fā)生級(jí)別預(yù)測(cè)模型,使用手持光譜儀ASD FieldSpec HandHeld檢測(cè)不同黃瓜霜霉病為害下黃瓜葉片高光譜反射率,并分析400~900 nm波段內(nèi)的光譜反射率參數(shù)與黃瓜霜霉病發(fā)生級(jí)別的相關(guān)關(guān)系和回歸關(guān)系。結(jié)果表明,在400~718 nm波段范圍內(nèi),隨著霜霉病級(jí)別的增加,黃瓜光譜反射率均有明顯上升,并在568~687 nm波段范圍內(nèi)達(dá)到極顯著相關(guān)水平;在825~900 nm波段范圍內(nèi),隨著霜霉病級(jí)別的增加,黃瓜光譜反射率均有明顯下降,在837~895 nm波段范圍內(nèi)達(dá)到極顯著相關(guān)水平。利用可見(jiàn)光區(qū)特征波長(zhǎng)點(diǎn)679 nm、近紅外區(qū)特征波長(zhǎng)點(diǎn)861 nm的光譜反射率以及歸一化植被指數(shù)(NDVI)、系統(tǒng)比值植被指數(shù)(RVI)構(gòu)建黃瓜霜霉病發(fā)生級(jí)別監(jiān)測(cè)預(yù)警方程,其中前3個(gè)預(yù)測(cè)模型均達(dá)到極顯著水平(P<0.01)。因此,可以利用高光譜信息數(shù)據(jù)來(lái)預(yù)測(cè)設(shè)施黃瓜霜霉病發(fā)生的級(jí)別。
關(guān)鍵詞:設(shè)施黃瓜;霜霉病;光譜測(cè)量;預(yù)測(cè)模型
中圖分類號(hào):S127;S436.421.1+1 ?文獻(xiàn)標(biāo)志碼: A ?文章編號(hào):1002-1302(2019)17-0117-03
設(shè)施黃瓜在生長(zhǎng)過(guò)程中,由于受到寡照、低溫、高濕等不良生態(tài)環(huán)境的影響,容易發(fā)生病害。黃瓜霜霉病是設(shè)施黃瓜生產(chǎn)中發(fā)生最廣泛、為害最嚴(yán)重的病害,病原為鞭毛菌亞門(mén)假霜霉屬古巴假霜霉菌(Pseudoperonospora cubensis)[1]。黃瓜霜霉病主要為害葉片。黃瓜苗期感染此病害后,子葉首先呈褪綠斑,慢慢發(fā)展成黃色斑,如果濕度大,子葉背面出現(xiàn)黑色霉層,隨后子葉逐漸全部變黃并干枯。黃瓜成株期感染此病害后,首先葉片上出現(xiàn)綠色水浸斑,慢慢變成黃綠色、淡褐色,并呈多角形,最后病斑合在一起,使得整個(gè)葉片干枯,并向上卷縮,如果濕度大,葉片背面呈灰黑色霉層,容易導(dǎo)致整個(gè)葉片枯死[2]。傳統(tǒng)黃瓜霜霉病診斷和預(yù)測(cè)預(yù)報(bào)方法是病理分析或經(jīng)驗(yàn)判斷,該方法不僅費(fèi)時(shí)、費(fèi)力,且效果差。隨著科技進(jìn)步,逐漸發(fā)展出高光譜技術(shù)來(lái)監(jiān)測(cè)農(nóng)業(yè)病蟲(chóng)害的危害程度,原理為作物受到為害時(shí),葉片細(xì)胞的活性、含水量、葉綠素含量等生理特性發(fā)生變化,表現(xiàn)為作物冠層光波反射特性發(fā)生變化,尤其是可見(jiàn)光波段和短波紅外波段的光波特征變化。因此,可以通過(guò)受害作物某些特征波長(zhǎng)的光波特性與正常作物光波特性的區(qū)別,來(lái)監(jiān)測(cè)作物病蟲(chóng)害危害級(jí)別[3-4]。因此,高光譜技術(shù)為快速、方便、準(zhǔn)確地進(jìn)行設(shè)施黃瓜霜霉病的監(jiān)測(cè)和預(yù)測(cè)預(yù)報(bào)提供了可能。
國(guó)內(nèi)外相關(guān)文獻(xiàn)已報(bào)道了高光譜遙感技術(shù)對(duì)重要病害發(fā)生程度的監(jiān)測(cè)。黃木易等測(cè)定了小麥感染條銹病后冠層光波反射率,結(jié)果發(fā)現(xiàn),光波反射率與條銹病的發(fā)生級(jí)別顯著相關(guān),據(jù)此用光譜反射率參數(shù)構(gòu)建條銹病病情指數(shù)的預(yù)測(cè)預(yù)報(bào)模型,從而為條銹病的適時(shí)、適量防治提供依據(jù)[5]。陳兵等在棉花各個(gè)生育期條件下,對(duì)感染黃萎病棉花葉片的光波反射率與葉綠素含量進(jìn)行相關(guān)分析,并同時(shí)分析光波反射率的一階微分方程與葉綠素含量的相關(guān)關(guān)系,在此基礎(chǔ)上構(gòu)建黃萎病侵染條件下,棉花葉片葉綠素含量的預(yù)測(cè)預(yù)報(bào)模型[6]。Franke等采用手持式光譜儀分析了冬小麥?zhǔn)艿桨追鄄?、條銹病危害后的光譜反射率,其結(jié)果表明,采用光譜技術(shù)能監(jiān)測(cè)、預(yù)警白粉病和條銹病的發(fā)生情況[7]。
當(dāng)光照射在植物上時(shí),近紅外波段的光大部分被植物反射回來(lái),紅光波段的光則大部分被植物吸收,通過(guò)近紅外波段和紅光波段反射率的線性或非線性組合,可以消除土地光譜的影響,得到的特征指數(shù)被稱為植被指數(shù)。植被指數(shù)的種類很多,其中比值植被指數(shù)(RVI)和歸一化植被指數(shù)(NDVI)[8-9]是目前應(yīng)用廣泛的2種植被指數(shù)。喬紅波對(duì)田間不同程度麥蚜為害初期的冬小麥品種鄭州891冠層反射光譜特征進(jìn)行分析,結(jié)果表明,麥蚜為害后小麥冠層的光譜反射率在近紅外區(qū)產(chǎn)生顯著變化,且隨著危害程度的加重其值逐步下降,并構(gòu)建了NDVI與小麥百株蚜量回歸預(yù)測(cè)方程[10]。楊可明等采用高光譜影像技術(shù)分析了條銹病侵染條件下的小麥冠層光譜,構(gòu)建了條銹病的預(yù)測(cè)預(yù)報(bào)模型,并提出新的植被指數(shù) MI-NDVI,該指數(shù)可消除作物不同生育階段的影響[11]。
目前,利用高光譜遙感數(shù)據(jù)構(gòu)建適用于設(shè)施黃瓜霜霉病為害監(jiān)測(cè)模型的研究還鮮有報(bào)道。為此,本研究采用高光譜技術(shù)探討黃瓜霜霉病為害后黃瓜反射光譜的變化規(guī)律,并選擇特征波長(zhǎng)點(diǎn),同時(shí)計(jì)算植被指數(shù),構(gòu)建設(shè)施黃瓜霜霉病發(fā)生級(jí)別的預(yù)測(cè)方程,以期為設(shè)施黃瓜霜霉病為害的快速評(píng)估提供理論依據(jù),進(jìn)而為生產(chǎn)上及時(shí)有效地防治設(shè)施黃瓜霜霉病打下基礎(chǔ)。
1 材料與方法
1.1 試驗(yàn)時(shí)間與地點(diǎn)
試驗(yàn)于2017年在陜西省西安市楊凌區(qū)設(shè)施黃瓜基地(34.31°N,108.07°E)進(jìn)行。溫室大棚長(zhǎng)為50 m,寬為6 m。
1.2 試驗(yàn)材料
供試黃瓜品種為目前陜西省關(guān)中地區(qū)主栽品種津春4號(hào)。
1.3 試驗(yàn)方法
1.3.1 試驗(yàn)設(shè)計(jì) 采用隨機(jī)區(qū)組設(shè)計(jì),將溫室大棚劃分為 6 m×3 m的15個(gè)小區(qū),在黃瓜霜霉病發(fā)生的始盛期通過(guò)人工接種和噴灑殺菌劑的方法,形成5個(gè)不同發(fā)病級(jí)別的小區(qū):0級(jí),全株無(wú)病;1級(jí),全株1/4及以下的葉片有少量病斑;2級(jí),全株1/2及以下、1/4以上的葉片有少量病斑或1/4及以下的葉片有較多的病斑;3級(jí),全株3/4及以下、1/2以上的葉片發(fā)病或全株1/4及以下的葉片全部枯黃;4級(jí),全株3/4以上的葉片發(fā)病或全株1/2及以上的葉片枯黃到整株枯黃。每個(gè)處理重復(fù)3個(gè)小區(qū)。整個(gè)黃瓜生長(zhǎng)期間無(wú)其他病蟲(chóng)害發(fā)生。
1.3.2 光譜測(cè)量 黃瓜光譜測(cè)量使用便攜式野外光譜儀ASD FieldSpec HandHeld(325~1 075 nm),該光譜儀采樣間隔為1.0 nm。在黃瓜霜霉病發(fā)病盛期測(cè)量光波反射率,測(cè)量時(shí)光譜儀傳感器探頭豎直朝下距黃瓜葉片20 mm,測(cè)量葉片前,要先對(duì)參考板和太陽(yáng)輻射光譜進(jìn)行標(biāo)定,測(cè)定中間要用硫酸鋇(BaSO4)板校正。每個(gè)樣方測(cè)定10次。測(cè)量后立即調(diào)查全株病害發(fā)生程度。
1.3.3 數(shù)據(jù)處理 采用ASD Viewspec-pro軟件讀取和分析源光譜數(shù)據(jù),將每個(gè)樣方測(cè)定的10條光譜數(shù)據(jù)的平均值作為該樣方光譜值,并將原始光譜反射強(qiáng)度數(shù)據(jù)轉(zhuǎn)化為反射率數(shù)據(jù)。從測(cè)量數(shù)據(jù)來(lái)看,325~399 nm和901~1 075 nm范圍內(nèi)的數(shù)據(jù)噪聲較大,因此只處理了400~900 nm范圍內(nèi)的光譜數(shù)據(jù)。利用400~900 nm波段光譜特征波長(zhǎng)點(diǎn)的光譜反射率計(jì)算RVI和NDVI[12],公式為
式中:RR、RNIR分別為620~700 nm紅光波段、701~900 nm近紅外波段中與黃瓜霜霉病發(fā)生級(jí)別相關(guān)性最為顯著的特征波長(zhǎng)點(diǎn)的光譜反射率。
利用上述植被指數(shù)及最佳波長(zhǎng)點(diǎn)的光波反射率數(shù)據(jù)作為預(yù)測(cè)預(yù)報(bào)模型參數(shù),同時(shí)采用SPSS-Statistics 17.0系統(tǒng)對(duì)霜霉病級(jí)別與400~900 nm范圍內(nèi)的光譜反射率進(jìn)行相關(guān)和回歸分析,用SigmaPlot 10.0系統(tǒng)進(jìn)行作圖。
2 結(jié)果與分析
2.1 不同級(jí)別病害光譜反射率分析
圖1為不同級(jí)別黃瓜霜霉病為害后溫室黃瓜光譜反射率曲線比較。健康植株與不同級(jí)別黃瓜霜霉病為害植株的光譜曲線在400~718 nm可見(jiàn)光波段的光譜反射率有明顯差別。這是由于在可見(jiàn)光譜段內(nèi),葉綠素的含量決定了作物的光譜特征,葉綠素能夠大量吸收可見(jiàn)光中藍(lán)光區(qū)(450 nm為中心)和紅光區(qū)(670 nm為中心)的輻射能,從而在這2個(gè)區(qū)域表現(xiàn)為吸收谷;在吸收谷中間波段吸收輻射能較少,形成綠色反射峰而使植物呈現(xiàn)綠色。當(dāng)黃瓜發(fā)生霜霉病時(shí),葉部出現(xiàn)黃褐色病斑,此病斑處葉綠素含量降低,造成對(duì)藍(lán)光區(qū)、紅光區(qū)輻射能的吸收降低,反射增強(qiáng),特別是紅光波段的反射率提高。因此在400~718 nm可見(jiàn)光譜段范圍內(nèi)黃瓜霜霉病光譜反射強(qiáng)度比正常部分有明顯提高,且隨著危害級(jí)別的升高,光譜反射率也逐漸提高。另外在825~900 nm的近紅外譜段內(nèi),健康植株與不同級(jí)別黃瓜霜霉病為害植株的光譜反射強(qiáng)度也有明顯的差異。這是由于在近紅外波段內(nèi),植物葉片細(xì)胞構(gòu)造是形成葉片光波特性的主要原因。在此波段內(nèi),葉片病斑處內(nèi)部結(jié)構(gòu)與正常部分有所不同,因此黃瓜霜霉病植株與健康植株在825~900 nm的近紅外波段內(nèi)反射強(qiáng)度有一定的差異,黃瓜霜霉病為害植株比健康植株的光譜反射強(qiáng)度有所降低,且隨著為害級(jí)別的升高,光譜反射率逐漸降低。
2.2 黃瓜霜霉病級(jí)別與光譜反射率相關(guān)性分析
對(duì)黃瓜霜霉病級(jí)別與光譜數(shù)據(jù)進(jìn)行相關(guān)分析。由圖2可知,400~900 nm波段范圍內(nèi),黃瓜霜霉病發(fā)病級(jí)別與不同波長(zhǎng)光譜反射率的相關(guān)性差異很大。在719~824 nm波段范圍內(nèi),不同級(jí)別黃瓜霜霉病的光譜反射率均不變,因此不作分析。在400~718 nm波段范圍內(nèi),黃瓜霜霉病級(jí)別與光譜反射率呈正相關(guān),其中568~687 nm波段范圍內(nèi),達(dá)到了極顯著水平。
2.3 基于高光譜數(shù)據(jù)的黃瓜霜霉病級(jí)別預(yù)測(cè)
以“2.2”節(jié)中確定的黃瓜霜霉病高光譜特征波長(zhǎng)點(diǎn)(可見(jiàn)光波段范圍內(nèi)的679 nm及近紅外波段范圍內(nèi)的861 nm)上的光譜反射率以及系統(tǒng)RVI和NDVI這2個(gè)植被指數(shù)為自變量,黃瓜霜霉病發(fā)生級(jí)別為因變量,建立黃瓜霜霉病發(fā)生情況的預(yù)測(cè)模型(表1)。
3 結(jié)論與討論
本研究分析了霜霉病為害設(shè)施黃瓜后,黃瓜光譜反射率對(duì)霜霉病級(jí)別的響應(yīng)規(guī)律,并明確了在400~718 nm波段范圍內(nèi),隨著霜霉病級(jí)別的增加,黃瓜光譜反射率均明顯上升,并在568~687 nm波段范圍內(nèi)達(dá)到極顯著相關(guān)水平;在 825~900 nm波段范圍內(nèi),隨著霜霉病級(jí)別的增加,黃瓜光譜反射率均明顯下降,在837~895 nm波段范圍內(nèi)達(dá)到極顯著相關(guān)水平。將可見(jiàn)光波段的679 nm、近紅外波段的861 nm作為設(shè)施黃瓜光譜特征波長(zhǎng)點(diǎn),計(jì)算系統(tǒng)RVI和NDVI,構(gòu)建黃瓜霜霉病監(jiān)測(cè)預(yù)警方程。通過(guò)對(duì)模型的F檢驗(yàn)表明,通過(guò)679、861 nm波長(zhǎng)點(diǎn)光譜反射率及NDVI構(gòu)建的模型均能很好地揭示霜霉病級(jí)別與設(shè)施黃瓜光譜反射率之間的關(guān)系,且通過(guò)861 nm波長(zhǎng)點(diǎn)光譜反射率構(gòu)建的模型優(yōu)于通過(guò)679 nm光譜反射率和NDVI構(gòu)建的模型。
農(nóng)作物的光波反射率具有很好的規(guī)律性:綠光波段(560 nm)光波吸收少,反射率較高;紅光波段(680 nm)光波吸收多,反射率較低;近紅外波段(720 nm)光波吸收顯著減少,反射率增加,隨后此區(qū)域反射率一直表現(xiàn)較高。病蟲(chóng)害為害農(nóng)作物時(shí),上述3個(gè)波段的光波反射率改變顯著,前期文獻(xiàn)表明,作物葉片葉綠素含量及葉肉細(xì)胞的構(gòu)造變化是其改變的根本原因[13-16]。本研究結(jié)果顯示,隨著黃瓜霜霉病發(fā)生級(jí)別的增加,綠光區(qū)和紅光區(qū)光波反射率明顯增加,而近紅外區(qū)光波反射率明顯下降,此結(jié)果與其他作物上的光譜反射率結(jié)果表現(xiàn)出一致性??捎命S瓜霜霉病病害生理學(xué)的現(xiàn)象來(lái)解釋:霜霉病菌為害黃瓜后,附著胞生成的侵入絲侵入葉片細(xì)胞間隙,而此時(shí)細(xì)胞不斷喪失水分,從而對(duì)近紅外波段反射率持續(xù)下降;隨著為害程度增加,侵入絲進(jìn)一步穿透葉片細(xì)胞,葉綠素含量進(jìn)一步減小,細(xì)胞水分持續(xù)喪失,使得黃瓜葉片綠光區(qū)光波反射率上升,而近紅外區(qū)光波反射率下降。本研究結(jié)果顯示,利用高光譜信息監(jiān)測(cè)黃瓜霜霉病是完全可行的。
本研究中的高光譜系統(tǒng)短時(shí)間采集到的黃瓜植株數(shù)據(jù)多,能很好地實(shí)現(xiàn)快速診斷的目的。通過(guò)該研究可以得到黃瓜霜霉病的最有效特征波長(zhǎng),然后根據(jù)該特征波長(zhǎng)可以進(jìn)行病害預(yù)測(cè)預(yù)報(bào),進(jìn)而提高診斷效率,真正達(dá)到快速、無(wú)損診斷農(nóng)作物病害的目的。
參考文獻(xiàn):
[1]熊 艷,王鶴冰,向華豐,等. 黃瓜霜霉病研究進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào),2016,32(1):130-135.
[2]劉建鳳,吉春明,蘇建坤,等. 幾種殺菌劑對(duì)黃瓜霜霉病菌室內(nèi)毒力測(cè)定[J]. 江蘇農(nóng)業(yè)科學(xué),2015,43(12):152-154.
[3]陳鵬程,張建華,雷勇輝,等. 高光譜遙感監(jiān)測(cè)農(nóng)作物病蟲(chóng)害研究進(jìn)展[J]. 中國(guó)農(nóng)學(xué)通報(bào),2006,22(2):388-391.
[4]Mathus T J,Maderia A C. High resolution spectroradiometry:spectral reflectance of field bean leaves infected by Botrytis fabae [J]. Remote Sensing of Environment,1993,45(1):107-116.
[5]黃木易,王紀(jì)華,黃文江,等. 冬小麥條銹病的光譜特征及遙感監(jiān)測(cè)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2003,19(6):154-158.
[6]陳 兵,王克如,李少昆,等. 蚜蟲(chóng)脅迫下棉葉光譜特征及其遙感估測(cè)[J]. 光譜學(xué)與光譜分析,2010,30(11):3093-309.
[7]Franke J,Menz G. Multi-temporal wheat disease detection by multi-spectral remote sensing[J]. Precision Agriculture,2007,8(3):161-172.
[8]Luo J H,Huang W J,Zhao J L,et al. Detecting aphid density of winter wheat leaf using hyperspectral measurements[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2013,6(2):690-698.
[9]Merrill S C,Holtzer T O,Peairs F B,et al. Modeling spatial variation of russian wheat aphid overwintering population densities in colorado winter wheat[J]. Journal of Economic Entomology,2009,102(2):533-541.
[10]喬紅波. 麥蚜、白粉病危害后冬小麥冠層高光譜的測(cè)量與分析[D]. 北京:中國(guó)農(nóng)業(yè)科學(xué)院,2004.
[11]楊可明,陳云浩,郭達(dá)志,等. 基于高光譜影像的小麥條銹病光譜信息探測(cè)與提取[J]. 光子學(xué)報(bào),2008,37(1):145-151.
[12]聶克艷,岳延濱,孫長(zhǎng)青,等. 辣椒茶黃螨危害的冠層光譜特性與危害級(jí)別反演[J]. 西南農(nóng)業(yè)學(xué)報(bào),2012,25(5):1689-1693.
[13]Aggarwal P K,Kalra N,Chander S,et al. Infocrop a generic simulation model for assessment of crop yields,losses due to pests and environmental impact of agro - ecosystems in tropical environments 1 model description[J]. Agricultural Systems,2006,89:1-25.
[14]浦瑞良,宮 鵬. 高光譜遙感及其應(yīng)用[M]. 北京:高等教育出版社,2000:52-96.
[15]楊可明,郭達(dá)志. 植被高光譜特征分析及其病害信息提取研究[J]. 地理與地理信息科學(xué),2006,22(4):31-34.
[16]Tsai F,Philpot W . Derivative analysis of hyperspectral data[J]. Remote Sensing of Environment,1998,66(1):41-51.