基于適宜尺度選擇的冬小麥赤霉病遙感監(jiān)測

李衛(wèi)國，陳 華，金正婷，張琤琤，葛廣秀，嵇福建

(江蘇省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)信息研究所，江蘇南京 210014)

由于作物病蟲害的發(fā)生與氣候環(huán)境條件變化關(guān)系明顯，近年來，有學(xué)者開始考慮將遙感光譜信息與氣候、環(huán)境因素結(jié)合進行作物病蟲害的監(jiān)測預(yù)報研究。如Bhattacharya等[7]將氣象數(shù)據(jù)及多個遙感植被指數(shù)結(jié)合進行芥菜腐爛病多階段遙感追蹤監(jiān)測。Zhang等[8]綜合遙感信息及氣象數(shù)據(jù)預(yù)報北京麥區(qū)區(qū)域尺度內(nèi)冬小麥白粉病。Otuka等[9]結(jié)合區(qū)域氣候環(huán)境條件研究了越南湄公河三角洲褐飛虱遙感監(jiān)測的可行性。Silva等[10]將MSG 衛(wèi)星數(shù)據(jù)和地表溫度數(shù)據(jù)結(jié)合進行番茄蟲害風(fēng)險區(qū)劃研究。當(dāng)冬小麥發(fā)生病蟲害時，不但植株長勢狀況(如郁閉程度、葉色、形態(tài)等)會發(fā)生改變，而且其生長的農(nóng)田氣候、環(huán)境也與正常年份或正常區(qū)域有所不同[11-14]，因此，利用遙感信息和生長狀況進行作物病害識別與估測，從研究機理和方法上也就成為可能[15-16]。

本研究在對高空間分辨率GF-1全色影像進行重采樣基礎(chǔ)上，分別將2 m×2 m、8 m×8 m和16 m×16 m的 GF-1全色波段影像與HJ-1衛(wèi)星30 m×30 m的多光譜影像融合，通過融合質(zhì)量評價及光譜特征比較，篩選適宜區(qū)域冬小麥田塊分布特征的融合遙感影像。進一步分析冬小麥長勢指標(biāo)、光譜信息與赤霉病病情指數(shù)之間的互作關(guān)系，篩選與赤霉病病發(fā)生關(guān)系較為緊密的主要長勢因子(或農(nóng)學(xué)參數(shù)指標(biāo))。基于主要長勢因子建立冬小麥赤霉病病情指數(shù)估測模型，結(jié)合適宜尺度遙感影像對研究區(qū)域冬小麥赤霉病進行空間變化監(jiān)測，以期為區(qū)域冬小麥生產(chǎn)上病害防治中的準(zhǔn)確信息獲取提供方法參考。

1 材料與方法

1.1 研究區(qū)域

研究區(qū)1(東經(jīng)120.417°～120.79°，北緯32.757°～33.156°)為江蘇省東臺市東北部及大豐市東南部交界處，位于江蘇省中部沿海，屬北亞熱帶季風(fēng)氣候，四季分明，雨量充沛，全年無霜期220 d左右，作物一年兩熟，是典型的稻麥連作區(qū)。該區(qū)域雖地勢平坦，但河網(wǎng)密布，田塊破碎，對冬小麥面積的遙感提取造成一定的困難，這種地理狀況在江蘇省中部具有普遍性。研究區(qū)2為江蘇省泰州興化市區(qū)(東經(jīng)119°48′～120°18′，北緯32°20′～32°42′)，該區(qū)域地理屬性與研究區(qū)域1基本相近。

1.2 遙感影像數(shù)據(jù)與預(yù)處理

選用2014年3月16日空間分辨率為2 m×2 m的GF-1全色波段影像和8 m×8 m多光譜影像，HJ-1選用同年同月21日空間分辨率為30 m×30 m的多光譜影像。由于數(shù)據(jù)獲取期間天氣晴好，成像質(zhì)量佳，有利于冬小麥生長監(jiān)測。

在ERDAS IMAGINE軟件中，運用多項式幾何校正模型，依據(jù)采集地面控制點對各遙感影像進行幾何校正，校正誤差控制在0.5個像元內(nèi)。在ENVI軟件中，利用FLAASH大氣校正模型對所需影像進行大氣校正。最后，利用研究區(qū)域的AOI文件，分別裁剪經(jīng)過大氣校正HJ-1多光譜影像和GF-1多光譜影像，得到研究區(qū)域的遙感(反射率)影像。

1.3 遙感影像融合

在ERDAS IMAGINE軟件中利用正射校正以及雙線性插值(bilinear interpolation)對GF-1全色波段影像重采樣，形成空間分辨率為8 m×8 m 和16 m×16 m的全色影像。通過主成分變換法，將3景全色波段影像分別與HJ-1多光譜影像進行融合，得到2 m×2 m、8 m×8 m、16 m×16 m的多光譜融合影像。另外，用相同融合方法將GF-1空間分辨率為2 m×2 m的全色波段影像與8 m×8 m多光譜影像進行融合得到2 m×2 m 的多光譜影像，用于研究區(qū)域冬小麥種植面積提取精度的驗證參考。

1.4 融合影像定量評價

選用均值、標(biāo)準(zhǔn)差及平均梯度值對影像進行定量(客觀)評價，其運算是通過Matlab編程計算實現(xiàn)。

在民宿這樣廣泛的概念基礎(chǔ)上，通過對“攜程”“去哪兒”等旅游網(wǎng)站的觀察，不難看出，民宿可分為城市民宿和鄉(xiāng)村民宿兩種，城市民宿主要存在于大型城市中，通過展現(xiàn)外來文化形成主題房間，體現(xiàn)著新潮的流行文化；鄉(xiāng)村民宿多依托景區(qū)或者地域特色資源而發(fā)展，核心內(nèi)涵體現(xiàn)為鄉(xiāng)村文化。

用影像中所有像元的灰度平均值(μ)，度量光譜信息的轉(zhuǎn)變，反映影像亮度大小，如果均值適中，則表示影像視覺效果良好。其表達式為：

(1)

標(biāo)準(zhǔn)差(Std)描述影像中像元值與均值的離散程度，可以評價影像的反差程度。一般情況下，標(biāo)準(zhǔn)差大，則影像反差大，所含信息多；反之，則說明影像反差小，色調(diào)單一，所含信息量小。其表達式為：

(2)

式(1)和(2)中，M、N分別為影像的行列數(shù)，F(xiàn)是融合影像灰度值，i、j為同一波段影像中各像元的行列號。

平均梯度(g)是影像在多維度方向上密度變化的速率，能夠敏感地反映影像中的微小細節(jié)差異，平均梯度越大，影像層次越多，視覺越清晰。其表達式為：

(3)

1.5 試驗區(qū)數(shù)據(jù)調(diào)查

在研究區(qū)1利用GPS定位建立試驗監(jiān)測樣點30個，每個監(jiān)測樣點均選擇冬小麥種植面積大，能代表附近冬小麥長勢的田塊。在冬小麥抽穗-揚花期調(diào)查各監(jiān)測樣點冬小麥赤霉病病情、生物量，采用對角線測定方法，分別用Green Seeker 冠層光譜儀、Sun Scan葉面積指數(shù)儀和SPAD儀分別測量試驗樣點冬小麥的地物光譜信息(包括紅光反射率和近紅外反射率)、葉面積指數(shù)(LAI)以及葉片葉綠素含量。取植株地上部于取樣袋中，置室內(nèi)烘箱105 ℃殺青20 min，75 ℃烘干并稱取其作為生物量數(shù)據(jù)。該數(shù)據(jù)用于分析建模。研究區(qū)2建立試驗監(jiān)測樣點15個，數(shù)據(jù)調(diào)查與處理同研究區(qū)1，數(shù)據(jù)用于模型驗證。

2 結(jié)果與分析

2.1 不同空間尺度融合影像質(zhì)量評價

三種空間尺度(2 m×2 m、8 m×8 m和16 m×16 m)影像融合的均值分別為81.2、79.26和79.44，三者之間相差不大，說明這三幅融合影像的平均亮度相近(表1)。三種融合尺度的平均梯度存在明顯差異，以16 m×16 m融合影像清晰度最好，細節(jié)表現(xiàn)力最強。標(biāo)準(zhǔn)差作為衡量信息豐富程度的重要指標(biāo)，其值越大，反差就越大，視覺信息就越明顯。三種空間尺度融合影像的平均標(biāo)準(zhǔn)差也存在差異， 且以16 m×16 m融合影像的信息量最多。總體上，16 m×16 m融合影像較適合研究區(qū)域冬小麥種植的田塊分布特征。

表1 不同空間尺度融合影像定量評價Table 1 Quantitative evaluation of fusion image on differet spatial scales

2.2 不同空間尺度融合影像植被指數(shù)特征

利用ENVI遙感軟件在三景融合影像中通過建立ROI感興趣點隨機提取88個冬小麥種植樣點，形成冬小麥樣點的光譜信息基礎(chǔ)數(shù)據(jù)，然后計算歸一化差值植被指數(shù)(NDVI)和比值植被指數(shù)(RVI)[11]。

以88個冬小麥種植樣點順序作為橫坐標(biāo)，兩種植被指數(shù)分別作為縱坐標(biāo)，依次生成散點圖(圖1)。圖1(a)和(b)分別是88個樣點所對應(yīng)的NDVI和RVI值結(jié)果，每幅圖都包含2 m×2 m、8 m×8 m、16 m×16 m三種尺度融合影像的樣點光譜信息。為了更直觀地展現(xiàn)三種不同尺度融合影像的植被指數(shù)差異，在圖1中添加線性趨勢線，輔以說明融合影像的植被指數(shù)的變化特征。從圖1(a)可以看出，冬小麥拔節(jié)期的NDVI處于0.19～0.54范圍，多數(shù)樣點的NDVI集中在0.30～0.43，樣點間NDVI存在明顯差異，說明不同田塊冬小麥之間的長勢不盡相同。結(jié)合冬小麥樣點趨勢線看，16 m×16 m融合影像的NDVI值明顯高于2 m×2 m和8 m×8 m融合影像，說明16 m×16 m融合影像光譜信息量較豐富。在圖1(b)中，RVI介于1.49～3.45范圍，密集分布于2.3左右，趨勢線也表明16 m×16 m融合影像所含光譜信息量多于2 m×2 m和8 m×8 m融合影像，融合影像光譜信息增加，有利于冬小麥的識別。

趨勢線由上到下分別為16 m×16 m、2 m×2 m和8 m×8 m。The trend lines belong to 16 m×16 m,2 m×2 m and 8 m×8 m from the top.

2.3 冬小麥長勢指標(biāo)與赤霉病發(fā)生間的關(guān)系

葉面積指數(shù)(LAI)、葉綠素含量(SPAD值)和地上部生物量是反映冬小麥長勢好壞的重要指標(biāo)。相關(guān)分析表明，在冬小麥抽穗-揚花期，這三個指標(biāo)與赤霉病病情指數(shù)均呈顯著正相關(guān)(P<0.05，n=30)(圖2)，線性擬合的決定系數(shù)R2分別為0.688、0.709和0.669，說明在冬小麥抽穗-揚花期， 這三個生長指標(biāo)越大，赤霉病的發(fā)病幾率越高，病情越嚴重，這可能因為試驗區(qū)4月份的氣溫較高，小麥生長較快，如果長勢過旺，田間郁閉程度會較高，冠層下面光照不足，容易引起赤霉病發(fā)生，同時雖然葉綠素含量較高，植株干物質(zhì)累積較多，但營養(yǎng)生長和生殖生長不協(xié)調(diào)，導(dǎo)致植株抗病性下降。因此，可將LAI、SPAD值和地上部生物量確定為容易導(dǎo)致冬小麥赤霉病發(fā)生的生長指標(biāo)。

圖2 冬小麥LAI(a)、SPAD(b)和地上部生物量(c)與赤霉病病情指數(shù)間的關(guān)系Fig.2 Relationship between LAI(a), SPAD(b) and aboveground biomass(c) and scab disease index of winter wheat

2.4 縣域冬小麥赤霉病發(fā)生的遙感估測

基于上述分析的冬小麥長勢指標(biāo)與赤霉病病情指數(shù)之間的特征關(guān)系，采用多元回歸方法，綜合LAI、SPAD值和地上部生物量，構(gòu)建冬小麥赤霉病發(fā)生的農(nóng)學(xué)參數(shù)估測模型(WSEMgi，winter wheat scab estimation model based on growth index)：

DIws = 2.01×LAI+1.25×SPAD+0.000 68×ABG-59.07

(4)

式中，DIws(disease index of winter wheat scab)是冬小麥赤霉病病情指數(shù)，LAI是冬小麥葉面積指數(shù)，SPAD是葉片葉綠素含量，AGB是地上部生物量(kg·hm-2)。

為實現(xiàn)對縣域空間尺度冬小麥赤霉病發(fā)生的遙感估測，選用經(jīng)過影像融合質(zhì)量評價較好的16 m×16 m重采樣影像作為影像數(shù)據(jù)源。引用金正婷[12]構(gòu)建的冬小麥長勢指標(biāo)遙感監(jiān)測模型，對WSEMgi模型進行修訂，即分別用不同的植被指數(shù)算法替代LAI、SPAD值和地上部生物量，以實現(xiàn)對進行冬小麥赤霉病病情指數(shù)的遙感估測。

LAI = 3.48×NDVI+2.79

(5)

SPAD =18.04×NDVI+35.29

(6)

AGB = -42.62×RVI2+1 119.4×RVI+5 365.42

(7)

在遙感影像處理軟件ERDAS IMAGINE的Modeler模塊中分別調(diào)入16 m×16 m重采樣影像、WSEMgi模型和生長指標(biāo)監(jiān)測模型，編程運行后即可生成冬小麥赤霉病病情指數(shù)遙感監(jiān)測圖。冬小麥病等級按照赤霉病病情指數(shù)大小初步分成4個等級：0級用綠色顯示，為長勢正常(無病害)；Ⅰ級為輕度病害(病情指數(shù)0～15%)，用黃色顯示；Ⅱ為較重病害(15%～25%)，用桔紅色表示；Ⅲ為嚴重病害(25%～35%)，用紅色表示。再經(jīng)過上述分級運算，添加圖例、比例尺和指北針信息，最終生成縣域冬小麥赤霉病病情遙感監(jiān)測信息專題圖(圖3)。利用研究區(qū)域2內(nèi)建立的15個試驗監(jiān)測樣點對冬小麥赤霉病病情遙感監(jiān)測信息進行精度驗證，REMS(平方根誤差)為10.5%，相對誤差為14.6%，說明冬小麥赤霉病病情遙感監(jiān)測算法的精度較為可靠。

圖3較為直觀地顯示出研究區(qū)(縣域)冬小麥抽穗-揚花期冬小麥赤霉病發(fā)生狀況，其中，縣域南部的周莊、茅山、陳堡、沈倫、臨城以及大朵鎮(zhèn)和東北部大營鎮(zhèn)的顏色顯示較紅黃，大多數(shù)田塊赤霉病發(fā)病較重，其東南部的一些鄉(xiāng)鎮(zhèn)如戴南、張郭、狄垛、陶莊、戴窯和合陳等鄉(xiāng)鎮(zhèn)的田塊均有不同程度發(fā)生，其他的鄉(xiāng)鎮(zhèn)也有零星發(fā)生。冬小麥揚花期是進行赤霉病防治的理想時期，針對田塊冬小麥的不同發(fā)病等級，可以分鄉(xiāng)鎮(zhèn)或分片區(qū)采用相應(yīng)的植保防治措施，以實現(xiàn)高效、環(huán)保的防治目的。

3 討 論

傳統(tǒng)的人工實地抽樣調(diào)查一直是冬小麥赤霉病監(jiān)測的主要手段，雖已在冬小麥病害防治方面發(fā)揮了重要作用，但對于大面積冬小麥病害發(fā)生而言，傳統(tǒng)監(jiān)測方法不僅需花費大量人力物力成本，且取樣范圍和樣本量較為有限，很難及時獲得大范圍的病情數(shù)據(jù)信息，在很大程度上影響到防治措施的有效實施。因此，需要采用及時、大范圍、低成本的遙感監(jiān)測方法來解決大面積農(nóng)田冬小麥赤霉病信息獲取問題。遙感方法是基于單個或多個植被光譜信息指標(biāo)進行冬小麥病情估算[2,4,6]，當(dāng)植被光譜信息指標(biāo)與冬小麥病情指標(biāo)關(guān)聯(lián)程度高時，估測精度較好，當(dāng)植被光譜信息指標(biāo)與冬小麥病情指標(biāo)關(guān)聯(lián)性弱時，幾乎無法估測。麥田偏施氮肥、密度大以及田間郁閉極發(fā)會誘發(fā)赤霉病的發(fā)生。本研究結(jié)合赤霉病發(fā)生的物候與生態(tài)特點，綜合冬小麥生長指標(biāo)(葉面積指數(shù)、葉片葉綠素含量和地上部生物量)和遙感光譜信息來監(jiān)測冬小麥的病發(fā)情況，在增強遙感監(jiān)測模型解釋性的同時，明顯提高了病情監(jiān)測的準(zhǔn)確性。

圖3 縣域冬小麥赤霉病發(fā)生狀況遙感監(jiān)測Fig.3 Remote sensing monitoring of winter wheat scab in county area

位居江淮區(qū)域的江蘇省興化市內(nèi)部水網(wǎng)密集，田塊較為破碎，種植結(jié)構(gòu)多樣。利用中低空間分辨率衛(wèi)星影像，由于其混合像元較多，常常會造成“異物同譜”、“同物異譜”現(xiàn)象，較難滿足該類縣區(qū)冬小麥生長監(jiān)測精度要求[11-12, 14]。因此，適宜空間分辨率(空間尺度)的遙感影像數(shù)據(jù)選取成為進行冬小麥長勢與病害有監(jiān)測的前提與必要條件。本研究通過對空間分辨率2 m×2 m的 GF-1衛(wèi)星全色影像進行重采樣，分別與空間分辨率30 m×30 m 的HJ-1衛(wèi)星多光譜影像進行融合，生成2 m×2 m、8 m×8 m 和16 m×16 m 3 種尺度多光譜影像。通過對多尺度融合影像質(zhì)量評價比較，發(fā)現(xiàn)16 m×16 m融合影像質(zhì)量較原始多光譜影像質(zhì)量發(fā)生明顯改善，突出豐富的植被光譜信息，比較適合江蘇省冬小麥田塊分布特征，有利于冬小麥遙感識別與信息提取。以16 m×16 m 影像作為影像數(shù)據(jù)源，引用WSEMgi模型和長勢指標(biāo)監(jiān)測模型，對縣域冬小麥抽穗-揚花期赤霉病病情進行遙感監(jiān)測，有效獲取到赤霉病發(fā)生狀況與空間分布信息，該方法對于縣域冬小麥病害防治措施制定有很好參考作用。

影響冬小麥赤霉病發(fā)生的因素很多[4，15,17]，本文僅對與冬小麥赤霉病發(fā)生較緊密的主要長勢指標(biāo)(葉面積指數(shù)、葉片葉綠素含量和生物量)進行了分析研究，其他影響因素的存在及其不確定性也是影響本研究所建模精度的主要原因。今后仍有必要綜合氣候、環(huán)境以及物候等影響因素，完善或改進冬小麥赤霉病遙感估測模型與方法，以更好實現(xiàn)對縣域冬小麥病害的精準(zhǔn)化監(jiān)測。