張 衡 潘 潔 巨云為 廖振峰

(南京林業(yè)大學(xué)，南京，210037)

責(zé)任編輯:程 紅。

松萎蔫病是由松材線蟲(chóng)引起的一種對(duì)松屬樹(shù)種的毀滅性災(zāi)害，該病最早于1905年在日本發(fā)現(xiàn)［1］。我國(guó)自1982年在南京中山陵風(fēng)景區(qū)首次發(fā)現(xiàn)該病以來(lái)迅速蔓延至全國(guó)各地，對(duì)森林資源和生態(tài)環(huán)境造成嚴(yán)重破壞。據(jù)了解，在我國(guó)每年由松材線蟲(chóng)病造成的直接和間接經(jīng)濟(jì)損失高達(dá)上億元［2］，因此，如何實(shí)現(xiàn)對(duì)松萎蔫病實(shí)時(shí)有效的監(jiān)測(cè)預(yù)報(bào)尤其是早期預(yù)測(cè)是迫切需要解決的問(wèn)題。

近年來(lái)，隨著信息技術(shù)的快速發(fā)展，高光譜遙感作為一種新的監(jiān)測(cè)手段被廣泛運(yùn)用到森林病蟲(chóng)害的監(jiān)測(cè)中［3］。高光譜遙感是高光譜分辨率遙感的簡(jiǎn)稱，它是指利用很多很窄(通常波段寬度＜10 nm)的電磁波波段從感興趣的物體上獲取相關(guān)數(shù)據(jù)，形成一條完整且連續(xù)的光譜曲線。它不僅具有波段窄、通道多、圖像與光譜能夠合二為一等優(yōu)點(diǎn)，而且它獲取的地物光譜信息連續(xù)且精細(xì)，因此高光譜遙感應(yīng)用于地物精細(xì)識(shí)別和分類的研究越來(lái)越引起人們的廣泛關(guān)注［4－7］。

本研究以成年馬尾松(Pinus massoniana)為接種材料，利用全波段(350～2 500 nm)高光譜儀于野外實(shí)測(cè)獲取冠層光譜數(shù)據(jù)，分析感病植株與對(duì)照植株在各個(gè)時(shí)期的光譜特征，尋找感病植株與健康植株在感病早期差異顯著的光譜特征參數(shù)，同時(shí)跟蹤測(cè)量各時(shí)期針葉的葉綠素質(zhì)量分?jǐn)?shù)，通過(guò)光譜特征參數(shù)結(jié)合感病植株的葉綠素質(zhì)量分?jǐn)?shù)定量反演，以期達(dá)到松材線蟲(chóng)病各病程定量監(jiān)測(cè)，特別是實(shí)現(xiàn)植株在出現(xiàn)肉眼可見(jiàn)感病特征前的早期監(jiān)測(cè)。

1 材料與方法

1.1 試驗(yàn)設(shè)計(jì)

試驗(yàn)地點(diǎn)位于江蘇省南京市棲霞山馬尾松林。采用貝曼氏漏斗法從南京地區(qū)自然死亡的黑松木段上分離出松材線蟲(chóng)，然后經(jīng)人工挑選，培養(yǎng)，供接種。根據(jù)當(dāng)?shù)氐臍夂蛱卣?，選擇在2012年6月份接種松材線蟲(chóng)，在7—9月份進(jìn)行光譜測(cè)量。在樣地中隨機(jī)選擇7 棵樹(shù)高相近年齡相同的健康成年馬尾松作為試驗(yàn)對(duì)象，其中5 棵接種松材線蟲(chóng)，2 棵為對(duì)照。

1.2 數(shù)據(jù)采集

光譜數(shù)據(jù)的采集過(guò)程包括測(cè)量健康與感病植株的冠層光譜反射率。數(shù)據(jù)采集分為野外冠層光譜數(shù)據(jù)采集和采樣針葉室內(nèi)葉綠素質(zhì)量分?jǐn)?shù)測(cè)定兩部分。采用美國(guó)ASD 公司生產(chǎn)的全波段光譜儀，在350～2 500 nm 波段范圍內(nèi)對(duì)選定的健康與感病的馬尾松植株冠層進(jìn)行連續(xù)光譜測(cè)量，自接種之日起，每隔10d 左右(依天氣而定)進(jìn)行定株光譜測(cè)量。測(cè)量時(shí)選擇晴朗無(wú)風(fēng)無(wú)云的天氣，同時(shí)為了減少太陽(yáng)高度角的變化對(duì)測(cè)量結(jié)果的影響，測(cè)量時(shí)間選擇在10:00—14:00。光譜數(shù)據(jù)采集過(guò)程中，每株選擇5 個(gè)冠層測(cè)量點(diǎn)，每次取10 組數(shù)據(jù)，每次測(cè)量前后均用硫酸鋇標(biāo)準(zhǔn)白板進(jìn)行校正，最后取平均值作為該馬尾松植株的冠層光譜反射率。光譜測(cè)量結(jié)束后，從光譜測(cè)量時(shí)對(duì)應(yīng)的相應(yīng)位置采集小段植株樣品裝入冰袋，帶回室內(nèi)即時(shí)進(jìn)行葉綠素質(zhì)量分?jǐn)?shù)的測(cè)定。

葉綠素質(zhì)量分?jǐn)?shù)的測(cè)定采用80%丙酮提取比色法［8］。將針葉剪碎、稱質(zhì)量并磨碎，用80%丙酮浸泡后過(guò)濾，利用UV－2102C 型紫外可見(jiàn)分光光度計(jì)(上海龍尼柯儀器有限公司)分別測(cè)量葉綠素a、葉綠素b 的吸光度，再依據(jù)公式Ca=12.21A663－2.81A646;Cb=20.13A646－5.03A663分別計(jì)算出其相應(yīng)的質(zhì)量分?jǐn)?shù)，然后將兩者相加即為葉綠素總量，為減小誤差，每組做3 次重復(fù)，最后取其平均值作為測(cè)試植株的葉綠素質(zhì)量分?jǐn)?shù)。

1.3 數(shù)據(jù)處理

在高光譜數(shù)據(jù)處理過(guò)程中進(jìn)行了預(yù)處理和光譜參數(shù)提取。首先，對(duì)采集的各時(shí)期高光譜數(shù)據(jù)進(jìn)行逐一核查，剔除其中的異常數(shù)據(jù)，然后求取50 組測(cè)量原數(shù)據(jù)的平均值作為每株馬尾松冠層反射率，然后對(duì)接種的5 株馬尾松植株光譜數(shù)據(jù)進(jìn)行差異分析。剔除與其他植株感病差異明顯的植株后，取其余感病株光譜數(shù)據(jù)均值作為感病馬尾松植株的冠層反射率。對(duì)照健康馬尾松的冠層反射率數(shù)據(jù)同理獲取。第二步，對(duì)健康及感病馬尾松植株的冠層反射率數(shù)據(jù)借助ASD 公司提供的View Specpro 軟件和EXCEL 數(shù)據(jù)分析工具進(jìn)行處理，提取各種光譜參數(shù)，分析其變化特征;然后在MATLAB 軟件中編程實(shí)現(xiàn)各波段原始反射率、導(dǎo)數(shù)光譜、不同光譜指數(shù)與葉綠素質(zhì)量分?jǐn)?shù)間的相關(guān)關(guān)系，選擇相關(guān)系數(shù)最大的光譜指數(shù)建立葉綠素質(zhì)量分?jǐn)?shù)的定量反演模型并進(jìn)行精度檢驗(yàn)。

具體光譜參數(shù)及處理方法如下。①基于原始反射率的參數(shù):510～570 nm 最大反射率值及對(duì)應(yīng)的波長(zhǎng)(綠峰);640～680 nm 最小反射率及對(duì)應(yīng)波長(zhǎng)(紅谷);反射特征區(qū)域的反射率，如反射峰與吸收谷中心波長(zhǎng)的位置，分別為455、550、680、980 nm。②植被指數(shù):兩波段組合的比值、差值和歸一化植被指數(shù)的計(jì)算方法參考表1，計(jì)算各植被指數(shù)的入選波段見(jiàn)表2。③導(dǎo)數(shù)光譜:采用一階導(dǎo)數(shù)光譜指數(shù)作為篩選參數(shù)［16］。

表1 多種植被指數(shù)計(jì)算方法及出處

表2 植被指數(shù)的波段組合

2 結(jié)果與分析

2.1 不同感病時(shí)間馬尾松原始反射光譜特征

由圖1可見(jiàn)，健康馬尾松各時(shí)期的光譜變化曲線遵循植被的反射光譜曲線特征，即在綠光區(qū)有一個(gè)較低的反射峰值，在近紅外區(qū)存在一個(gè)較高的反射峰值。在光譜測(cè)量的6 個(gè)時(shí)期內(nèi)反射光譜曲線變化規(guī)律相似，無(wú)明顯差異。其中，只有8月18日光譜曲線的反射率明顯高于其它各時(shí)期，尤其在近紅外區(qū)，其它各時(shí)期的反射率都逐漸趨于平穩(wěn)，而它的反射率卻呈現(xiàn)出線性增加的趨勢(shì)。究其原因，主要是因?yàn)樵谶@個(gè)時(shí)期對(duì)照植株長(zhǎng)出了較多的新葉，葉綠素整體含量明顯升高，使得這個(gè)時(shí)期的光譜反射特征出現(xiàn)了這樣的變化趨勢(shì)。

圖1 健康植株(左)與感病植株(右)各時(shí)期反射光譜

感病馬尾松光譜曲線與健康植株相比有著明顯的差異，尤其表現(xiàn)在感病中后期，感病植株在綠光區(qū)的反射特征逐漸減弱，反射峰值逐漸減小，直至病株完全枯萎，光譜曲線上的反射峰或吸收谷特征明顯消失，趨向于一條平直曲線。其中，8月18日測(cè)量的感病植株冠層反射率和對(duì)照類似，也表現(xiàn)出反射率的絕對(duì)值顯著高于其它測(cè)量時(shí)間的反射率值，原因與對(duì)照相同，此時(shí)的感病植株在枯萎變黃的同時(shí)，也有部分新葉長(zhǎng)出，使得感病植株的反射率值明顯要高于其它時(shí)期。

由此可見(jiàn)，健康馬尾松與感病馬尾松植株冠層原始反射率光譜的變化差異比較顯著，特別是在感病的中后期，感病植株的反射峰與吸收谷特征都逐漸減弱直至消失。

2.2 不同感病時(shí)間馬尾松特征光譜指數(shù)差異

通過(guò)對(duì)健康植株與感病植株原始反射光譜的對(duì)比分析，發(fā)現(xiàn)感病植株與健康植株的光譜曲線在感病的中后期差異較為顯著，但在感病早期，冠層葉片顏色變化不明顯，兩者原始反射光譜曲線間的差異也不顯著，單從光譜反射率來(lái)看并不能判別馬尾松是否感染了松材線蟲(chóng)病，因此，需要進(jìn)一步地尋找能在植株出現(xiàn)肉眼可見(jiàn)感病特征前有明顯差異的光譜指數(shù)來(lái)進(jìn)行早期的區(qū)分與判別。通過(guò)查閱資料，首先選擇了表1中所列舉的7 個(gè)植被指數(shù)作為敏感參數(shù)來(lái)進(jìn)行松萎蔫病的早期判別，這幾個(gè)指數(shù)是目前光譜分析較為常用的植被指數(shù)，在分析植被的反射率時(shí)能夠較好地反映出植被的光譜特征，具有一定的代表性［9－15］。

表3 健康植株與感病植株不同時(shí)間的各植被指數(shù)

根據(jù)表1和表2中所給的公式分別計(jì)算感病植株與健康植株在不同感病階段對(duì)應(yīng)的植被指數(shù)。因?yàn)楦鞑ǘ谓M合的對(duì)比結(jié)果呈類似的變化規(guī)律，因此選取一種作為代表在表3中列出。通過(guò)對(duì)比分析(表3)可以看出，感病植株與健康植株的比值植被指數(shù)和差值植被指數(shù)的差異主要表現(xiàn)在感病的中后期，結(jié)構(gòu)敏感色素指數(shù)和抗大氣植被指數(shù)2 直到感病末期才有較為明顯的差異，而感病植株的歸一化植被指數(shù)、加強(qiáng)植被指數(shù)和土壤調(diào)整植被指數(shù)在整個(gè)感病階段都與健康植株沒(méi)有太大差異。由此可見(jiàn)，感病植株的這些植被指數(shù)在植株出現(xiàn)肉眼可見(jiàn)感病特征前與健康植株相比無(wú)明顯差異，因此不能用來(lái)作為馬尾松松萎蔫病早期判別的光譜指數(shù)。

2.3 不同感病時(shí)間馬尾松一階導(dǎo)數(shù)光譜

一階微分光譜在辨別細(xì)微光譜差異上具有一定的優(yōu)勢(shì)，能夠較好地反映出植物的本質(zhì)特征［17］，因此將原始反射光譜進(jìn)行一階微分處理，進(jìn)一步探尋在感病早期能夠表現(xiàn)出明顯差異的光譜參數(shù)。經(jīng)過(guò)仔細(xì)的篩查與對(duì)比分析，發(fā)現(xiàn)感病植株與健康植株在波長(zhǎng)593 nm 處一階導(dǎo)數(shù)光譜(FD593)有較明顯的差異(表4)，健康植株在波長(zhǎng)593 nm 處的一階導(dǎo)數(shù)光譜有波動(dòng)性變化但不顯著，其平均值為(0.025±0.005)，但感病植株的值呈現(xiàn)出一種先減小后逐漸增大的趨勢(shì)，這種變化規(guī)律與不同感病階段植株葉綠素含量的變化特征吻合(表4)，同時(shí)在植株出現(xiàn)肉眼可見(jiàn)感病特征前的發(fā)病早期階段，F(xiàn)D593值會(huì)減小到0.01 以下，因此，可以利用在波長(zhǎng)593 nm 處的一階導(dǎo)數(shù)光譜來(lái)對(duì)馬尾松是否感染了松萎蔫病進(jìn)行早期定性判別。

表4 健康與感病馬尾松在波長(zhǎng)593 nm 處的發(fā)射率一階微分與葉綠素質(zhì)量分?jǐn)?shù)

2.4 馬尾松不同病程葉綠素質(zhì)量分?jǐn)?shù)定量反演

研究發(fā)現(xiàn)，在測(cè)量的感病植株各種生理生化指標(biāo)中，葉綠素質(zhì)量分?jǐn)?shù)的變化動(dòng)態(tài)與光譜變化特征有明顯的相關(guān)關(guān)系，并且感病植株在感病過(guò)程中葉綠素質(zhì)量分?jǐn)?shù)的變化規(guī)律明顯，更能從外觀上直接反應(yīng)感病植株不同感病階段的變化［18］。因此，本研究選取植株的葉綠素質(zhì)量分?jǐn)?shù)進(jìn)行定量反演以實(shí)現(xiàn)松材線蟲(chóng)病的定量監(jiān)測(cè)。

由圖2可以看出，隨著馬尾松的生長(zhǎng)，對(duì)照組葉綠素質(zhì)量分?jǐn)?shù)有波動(dòng)性變化但不顯著，感病植株的葉綠素質(zhì)量分?jǐn)?shù)在發(fā)病初期會(huì)略有升高，但隨著病情的加重，葉綠素質(zhì)量分?jǐn)?shù)會(huì)顯著降低，在發(fā)病中后期，隨著針葉的不斷失水、褪綠、褐變，葉綠素質(zhì)量分?jǐn)?shù)會(huì)呈明顯下降趨勢(shì)，這一結(jié)果與很多研究中顯示葉綠素質(zhì)量分?jǐn)?shù)變化的情況［18－19］相吻合。

利用單波段光譜反射率與葉綠素質(zhì)量分?jǐn)?shù)進(jìn)行相關(guān)性分析，同時(shí)也將波段進(jìn)行組合后獲取的多種光譜變換參數(shù)與葉綠素質(zhì)量分?jǐn)?shù)進(jìn)行相關(guān)性分析，計(jì)算二者的相關(guān)系數(shù)，將最終結(jié)果進(jìn)行篩選，發(fā)現(xiàn)在波長(zhǎng)1 405 nm 處的反射率與葉綠素質(zhì)量分?jǐn)?shù)的相關(guān)系數(shù)最大，達(dá)到了0.884，呈現(xiàn)出顯著的正相關(guān)關(guān)系(圖2)，因此用波長(zhǎng)1 405 nm 處的反射率來(lái)建立葉綠素質(zhì)量分?jǐn)?shù)的定量反演模型。模型為:CHL=－1.74(X1405)2+4.72X1405－0.76。式中:CHL為感病植株的葉綠素質(zhì)量分?jǐn)?shù)。

圖2 馬尾松葉綠素質(zhì)量分?jǐn)?shù)隨波長(zhǎng)1 405 nm 反射率的變化

通過(guò)在波長(zhǎng)1 405 nm 處的反射率獲取，可以直接計(jì)算出感病馬尾松植株的葉綠素質(zhì)量分?jǐn)?shù)。而研究發(fā)現(xiàn)，導(dǎo)數(shù)光譜指數(shù)對(duì)傳感器的靈敏度和信噪比的要求很高，用微分因子構(gòu)建的模型很難在影像反演中獲得令人滿意的結(jié)果［19］，因此，不考慮用導(dǎo)數(shù)光譜指數(shù)作為敏感波段來(lái)構(gòu)建反演模型。

綜上所述，可以利用在波長(zhǎng)593 nm 處的導(dǎo)數(shù)光譜并結(jié)合相應(yīng)的葉綠素質(zhì)量分?jǐn)?shù)的定量演算來(lái)早期判別馬尾松是否感染了松萎蔫病，在植株出現(xiàn)肉眼可見(jiàn)感病特征前，利用高光譜遙感獲取其光譜反射率進(jìn)行分析，當(dāng)FD593的值為0～0.01，相應(yīng)的葉綠素質(zhì)量分?jǐn)?shù)超過(guò)(1.50 ±0.10)mg·g－1時(shí)，就可以確定馬尾松感染了松萎蔫病。

3 結(jié)論與討論

結(jié)合馬尾松松萎蔫病發(fā)生時(shí)的生態(tài)學(xué)特征與相應(yīng)的高光譜特征，同時(shí)結(jié)合植株不同感病時(shí)期的光譜特征與葉綠素質(zhì)量分?jǐn)?shù)的相關(guān)性分析，從定性與定量?jī)蓚€(gè)方面判斷馬尾松是否感染了松萎蔫病，以此實(shí)現(xiàn)植株在出現(xiàn)肉眼可見(jiàn)感病特征前的早期監(jiān)測(cè)。研究表明:①感病與健康馬尾松的原始光譜反射率差異明顯。主要表現(xiàn)在，感病中后期感病植株在綠光區(qū)的反射峰逐漸減弱，而在紅光區(qū)的吸收谷會(huì)逐漸增強(qiáng)，隨著病害程度的進(jìn)一步加深，反射率、吸收谷都會(huì)慢慢消失，反射光譜曲線也逐漸變?yōu)橐粭l平滑的直線。但在感病早期二者之間的差異不顯著。②馬尾松在感染松萎蔫病后，葉綠素質(zhì)量分?jǐn)?shù)呈現(xiàn)先略微升高后下降的趨勢(shì)，并且植株的葉綠素質(zhì)量分?jǐn)?shù)與波長(zhǎng)1 405 nm 處的光譜反射率有顯著的相關(guān)性，可以據(jù)此來(lái)實(shí)現(xiàn)感病植株葉綠素質(zhì)量分?jǐn)?shù)的定量預(yù)測(cè)。模型為:CHL= －1.74(X1405)2+4.72X1405－0.76，式中，CHL為感病植株的葉綠素質(zhì)量分?jǐn)?shù)。③通過(guò)對(duì)多種光譜特征參數(shù)的篩選，選擇波長(zhǎng)593 nm 處反射率的一階微分作為敏感指數(shù)，并結(jié)合相應(yīng)的葉綠素質(zhì)量分?jǐn)?shù)來(lái)判斷馬尾松在出現(xiàn)肉眼可見(jiàn)感病特征前是否感染了松萎蔫病。當(dāng)FD593的值為0～0.01，相應(yīng)的葉綠素質(zhì)量分?jǐn)?shù)超過(guò)(1.50±0.10)mg·g－1時(shí)，就可以確定馬尾松感染了松萎蔫病。

我國(guó)對(duì)松萎蔫病的致病機(jī)理及相應(yīng)的防治措施和控制方法的研究已日漸成熟，但對(duì)該病的早期監(jiān)測(cè)預(yù)報(bào)的研究還比較少見(jiàn)，而利用高光譜遙感技術(shù)對(duì)該病進(jìn)行早期的監(jiān)測(cè)預(yù)報(bào)仍處于探索階段。由于松材線蟲(chóng)病在發(fā)病早期病害木自身并無(wú)明顯癥狀，而等到植株顯現(xiàn)出肉眼可見(jiàn)的感病特征則會(huì)在短時(shí)期內(nèi)迅速萎蔫死亡，因此，在未知感病與否的情況下，如何判斷感病并實(shí)現(xiàn)感病植株的早期監(jiān)測(cè)對(duì)松材線蟲(chóng)病的預(yù)防與治理具有重要的實(shí)踐意義與應(yīng)用價(jià)值。

本研究的結(jié)論表明，通過(guò)對(duì)健康與感病馬尾松連續(xù)的光譜監(jiān)測(cè)，結(jié)合相應(yīng)時(shí)期植株的生理生化參數(shù)(如葉綠素質(zhì)量分?jǐn)?shù))的同步分析，在發(fā)病早期可以較為準(zhǔn)確地判別出被感染的馬尾松，這為深入研究利用高光譜遙感技術(shù)來(lái)實(shí)現(xiàn)松萎蔫病的早期監(jiān)測(cè)預(yù)報(bào)提供了依據(jù)。但是，由于時(shí)間和試驗(yàn)條件的限制，本研究在建模時(shí)所用的數(shù)據(jù)有限，需要在將來(lái)的研究中對(duì)定量模型進(jìn)行進(jìn)一步的檢驗(yàn)和驗(yàn)證，同時(shí)，不同樹(shù)種間、不同生境、不同病蟲(chóng)害的判別是否會(huì)有差異都有待深入的研究與探索。

［1］ Kwon T S，Shin J H，Lim J H，et al.Management of pine wilt disease in Korea through preventative silvicultural control［J］.Forest Ecology and Management，2011，261(3):562 －569.

［2］ 趙博光.松萎蔫病致病新機(jī)理論文選編［M］.香港:中國(guó)科學(xué)文化出版社，2008:279 －289.

［3］ Sankaran S，Mishra A，Ehsani R，et al.A review of advanced techniques for detecting plant diseases［J］.Computers and Electronics in Agriculture，2010，72(2):1 －13.

［4］ 陳鵬程，張建華，雷勇輝，等.高光譜遙感監(jiān)測(cè)農(nóng)作物病蟲(chóng)害研究進(jìn)展［J］.中國(guó)農(nóng)學(xué)通報(bào)，2006，22(2):388 －391.

［5］ 何誠(chéng)，馮仲科，袁進(jìn)軍，等.高光譜遙感技術(shù)在生物多樣性保護(hù)中的應(yīng)用研究進(jìn)展［J］.光譜學(xué)與光譜分析，2012，32(6):1628－1632.

［6］ Vane G，Goetz A F H.Terrestrial imaging spectrometry:current status，future trends［J］.Remote Sensing of Environment，1993，44(2/3):117 －126.

［7］ Green R O，Eastwood M L，Sarture C M，et al.Imaging spectroscopy and the airborne visible/infrared imaging spectrometer(AVIRIS)［J］.Remote Sensing Environment，1998，65(3):227－248.

［8］ 郝再彬.植物生理實(shí)驗(yàn)［M］.哈爾濱:哈爾濱工程大學(xué)出版社，2004.

［9］ Pearson R L，Miller L D.Remote mapping of standing crop biomass for estimation of the productivity of the shortgrass prairie［R］.Michigan:Willow Run Laboratories，1972:1357 －1381.

［10］ Richardson A J，Wiegand C L.Distinguishing vegetation from soil background information［J］.Photogrammetric Engineering and Remote Sensing，1977，43(12):1541 －1552.

［11］ Rouse J W Jr.Monitoring the vernal advancement of retrogradation of natural vegetation［R/OL］.［2014 －02 －11］.http://ntrs.nasa.gov/search.jsp?R=19730009608.

［12］ Justice C O，Vermote E，Townshend J R G，et al.The moderate resolution imaging spectroradiometer (MODIS):land remote sensing for global change research［J］.IEEE Transactions on Geoscience and Remote Sensing，1998，36(4):1228 －1250.

［14］ Gitelson A A，Kaufman Y J，Stark R，et al.Novel algorithms for remote estimation of vegetation fraction［J］.Remote Sensing of Environment，2002，80(1):76 －87.

［15］ Huete A R.A soil-adjusted vegetation index (SAVI)［J］.Remote Sensing of Environment，1988，25:295 －309.

［16］ 喬紅波，程登發(fā)，孫京瑞，等.麥蚜對(duì)小麥冠層光譜特性的影響研究［J］.植物保護(hù)，2005，31(2):21 －26.

［17］ 許勇.江蘇沿海無(wú)機(jī)氮磷營(yíng)養(yǎng)鹽遙感監(jiān)測(cè)研究［D］.南京:南京師范大學(xué)，2009.

［18］ 王曉堂.基于高光譜的松萎蔫病動(dòng)態(tài)變化的研究［D］.南京:南京林業(yè)大學(xué)，2011.

［19］ 陳玉惠，葉建仁，魏初獎(jiǎng)，等.松材線蟲(chóng)侵染對(duì)馬尾松、黑松水分及其相關(guān)代謝的影響［J］.植物病理學(xué)報(bào)，2005，35(3):201－207.