馬云強趙寧王華張海燕

(1.西南林業(yè)大學(xué)生物多樣性保護學(xué)院/云南省森林災(zāi)害預(yù)警與控制重點實驗室，云南 昆明 650224；2.西南林業(yè)大學(xué)生命科學(xué)學(xué)院，云南 昆明 650224；3.云南這里信息技術(shù)有限公司，云南 昆明 650224)

松材線蟲病(Bursaphelenchusxylophilus)，亦稱松樹萎蔫病或“松樹癌癥”[1]，是國際上公認(rèn)的重要檢疫性有害生物，松樹感染后數(shù)月枯死，3～5a可蔓延傳播至整片松林，導(dǎo)致大面積的生態(tài)破壞，給國民經(jīng)濟造成重大損失[2]。根據(jù)國家林業(yè)局2020年松材線蟲病疫區(qū)公告，截至2020年已有20個省(市)，比2018年新增2個省(市)。松材線蟲病防治難度大，常規(guī)的監(jiān)測采用單一的人工地面調(diào)查手段進行[3]，這種方式耗時長、數(shù)據(jù)準(zhǔn)確度不高，難以有效控制松材線蟲病的蔓延[4]。隨著衛(wèi)星遙感影像分類技術(shù)的提高，黃慧萍[5]、齊樂[6]利用面向?qū)ο蠓诸惙椒ǚ謩e對分類尺度、支持向量機的參數(shù)進行研究，提高了遙感影像分類精度。張學(xué)敏[7]、馬宇龍[8]、梁守真[9]、呂利利[10]、董連英[11]利用CART決策樹分類法分別對病害松樹、龍祥島濕地、橡膠林、沙漠化土地、植被類型進行分類，提取效果較好。遙感影像可以實現(xiàn)大范圍的松林監(jiān)測，但受分辨率和時效性影響，難以做到實時監(jiān)測和單株精準(zhǔn)監(jiān)測。張紅梅[12]對近年來無人機遙感技術(shù)在松線蟲病監(jiān)測研究進行了總結(jié)，得出無人機遙感在松線蟲病的監(jiān)測上有巨大空間。近年來移動設(shè)備的日益更新，在監(jiān)測方面也有很大的作用。

此次研究以衛(wèi)星遙感監(jiān)測為基礎(chǔ)，利用無人機遙感監(jiān)測和便攜式移動設(shè)備進行驗證，有效提高松材線蟲病疫木提取的效率和精度?；跀?shù)據(jù)挖掘、數(shù)據(jù)融合、數(shù)據(jù)協(xié)同和數(shù)據(jù)同化等關(guān)鍵技術(shù)，獲得更加準(zhǔn)確數(shù)據(jù)支持的林業(yè)監(jiān)測感知體系。

1 研究區(qū)域概況

東錢湖鎮(zhèn)地處浙江鄞州東南，處于全國海岸線中段。區(qū)域內(nèi)自然資源豐富，植被種類300余種，山地森林覆蓋率92.4%。生態(tài)環(huán)境優(yōu)美，湖面開闊，岸線曲折，四周群山環(huán)抱，森林蒼郁；氣候良好，屬亞熱帶季風(fēng)氣候，全年溫和濕潤，雨量充沛，年平均氣溫15.4℃，面積145km2，規(guī)劃控制面積230km2，人口4.49萬，轄2個社區(qū)、4個居民區(qū)、36個行政村。

所選研究區(qū)域位于東錢湖正南方，面積為71km2，見圖1，該研究區(qū)林地集中，多為馬尾松林地，且該區(qū)域是東錢湖鎮(zhèn)中的松材線蟲病感染重災(zāi)區(qū)，具有代表性。因此，利用衛(wèi)星遙感影像，通過面向?qū)ο蟮腃ART決策樹方法進行遙感解譯，對東錢湖旅游度假區(qū)松材線蟲病發(fā)生區(qū)域的監(jiān)測提供了科學(xué)依據(jù)。

2 衛(wèi)星遙感影像分類

2.1 數(shù)據(jù)源

BJ-2號星座系統(tǒng)設(shè)計壽命7a，由3顆0.8m全色、3.2m多光譜的光學(xué)遙感衛(wèi)星組成(藍(lán)波段：0.44～0.51nm，綠波段：0.51～0.59nm，紅波段：0.60～0.67nm，近紅外波段：0.76～0.91nm)，可提供覆蓋全球、空間和時間分辨率俱佳的遙感衛(wèi)星數(shù)據(jù)。此外，松材線蟲病的的病發(fā)時間多發(fā)生在5—10月，尤其在10月發(fā)病最為明顯，發(fā)病樹木呈紅色或紅褐色，健康樹木呈綠色，為了能把健康樹木和發(fā)病樹木分開，故以東錢湖鎮(zhèn)2018年9月28日的BJ-2號影像為信息源進行研究，如圖2所示的原始全色影像和高光譜影像。另外，將野外實際取樣數(shù)據(jù)作為輔助數(shù)據(jù)，用于典型地物的樣本選擇和精度驗證。

2.2 技術(shù)流程

對BJ-2號遙感影像進行預(yù)處理，處理過程包括輻射定標(biāo)、正射校正、影像融合。接著對預(yù)處理的結(jié)果影像進行面向?qū)ο蠓指?、?gòu)建特征指數(shù)、提取分類閾值、進行CART分類，提取病疫木。圖3為完成預(yù)處理的遙感影像，圖4為技術(shù)流程。

2.3 研究方法

2.3.1 面向?qū)ο蠓指?/p>

面向?qū)ο蠓椒ㄒ酝|(zhì)性像元組成的“對象”為分類基礎(chǔ)，其核心算法是影像分割[13]。在合理的分割尺度下得到的“對象”破碎化程度不一，分割尺度直接決定了分類結(jié)果的精度。本文使用eCognition軟件進行遙感影像的分割。eCognition中有3種分割模式，分別是棋盤式分割、四叉樹分割與多尺度分割，其中多尺度分割的效果最優(yōu)[14]。多尺度分割可充分利用目標(biāo)地物的光譜和形狀、紋理和尺寸等空間信息，提高了分類的準(zhǔn)確性和穩(wěn)定性，使分類效率和分類精度得以兼顧，因此確定多尺度分割的參數(shù)尤為重要。本文以分割對象異質(zhì)性最小、分割后對象與目標(biāo)地物吻合度最高為依據(jù)，根據(jù)目標(biāo)地物的影像特征，經(jīng)多次實驗確定：影像分割尺度選取為80，形狀因子、緊致度因子分別選取為0.5，各波段權(quán)重為1。

2.3.2 特征指標(biāo)選擇

CART決策樹是使用監(jiān)督學(xué)習(xí)的方法，從1組無規(guī)律分布的復(fù)雜數(shù)據(jù)中尋找最佳分類體系，構(gòu)建1個二分支模型，實現(xiàn)對數(shù)據(jù)的分類和預(yù)測。模型構(gòu)建時，每1分裂節(jié)點的異質(zhì)性下降程度作為該節(jié)點是否生成的判別準(zhǔn)則，節(jié)點異質(zhì)性越低，雜質(zhì)含量就越少，該節(jié)點的分類精度越高[6]；每個節(jié)點的特征參數(shù)選擇尤為重要，特征參數(shù)作為分類屬性用來整個影像的反演。本文研究區(qū)域遙感影像包含建筑物、道路、水系、耕地、植被等，特征參數(shù)的選擇要充分考慮遙感影像的地物環(huán)境。本文最終目的是提取松木病疫木，即分類結(jié)果精確到林木植被即可。

基于以上考慮，在常見的遙感特征中，選擇歸一化水體指數(shù)(MNDWI)和歸一化植被指數(shù)(NDVI)分別對水體、植被進行分類，選擇30m高程DEM數(shù)據(jù)對耕地進行分類，選擇BJ-2號影像的紅綠波段比率對病疫木進行分類。特征指標(biāo)選擇充分考慮高分辨率遙感影像的光譜信息，結(jié)合各波段均值及標(biāo)準(zhǔn)差等多種特征參數(shù)進行決策樹的構(gòu)建。

2.3.3 決策樹模型構(gòu)建

基于對以上特征指標(biāo)的分析結(jié)果，利用上述特征指標(biāo)構(gòu)建5層結(jié)構(gòu)的決策樹，采用基于決策樹規(guī)則的面向?qū)ο蟮姆椒?，逐層分類，提取松材線蟲病病疫木信息，如圖5所示。節(jié)點層0為總節(jié)點；節(jié)點層1，應(yīng)用MNDWI通過設(shè)定閾值分離水體和非水體；節(jié)點層2，應(yīng)用NDVI分離植被和非植被；節(jié)點層3，設(shè)定DEM閾值區(qū)分植被中的林地和其它林地；節(jié)點層4，基于BJ-2號影像的紅綠波段比率、各波段均值及標(biāo)準(zhǔn)差提取松材線蟲病病疫木對象。

2.4 分類結(jié)果與精度評價

基于BJ-2號影像數(shù)據(jù)，采用面向?qū)ο蟮腃ART決策樹分類方法，選擇MNDWI、NDVI、DEM地形數(shù)據(jù)以及影像的紅綠波段比率作為特征參數(shù)對試驗區(qū)地物進行逐級區(qū)分，提取松材線蟲病病疫木分布信息，結(jié)果如圖6所示。從空間分布上看，病疫木提取結(jié)果和遙感影像的表現(xiàn)一致。根據(jù)野外調(diào)查點及使用ArcGIS生成的隨機樣本點作為檢驗樣點(共450個)，對分類結(jié)果進行精度評價，結(jié)果如表1所示。試驗區(qū)總體分類精度為85.70%，Kappa系數(shù)為0.8213。病疫木驗證樣本共179個，正確分類的有164個；健康木驗證樣本176個，正確分類的有166個。這2種目標(biāo)地物的生產(chǎn)精度都在91.6%以上，用戶精度都在90.7%以上，分類精度符合要求。

考慮到檢驗樣點代表的范圍太小有局域性，為了更精確地評價提取結(jié)果，在結(jié)果圖上隨機選擇100個病疫木樣區(qū)作為檢驗樣本，采用目視判讀的方式對選擇樣本的病疫木提取精度進行評價。考慮到區(qū)域地物的破碎性，將每個樣區(qū)大小設(shè)置為125×125像元(即100m×100m)。相比于點數(shù)據(jù)的精度評價結(jié)果，基于樣區(qū)數(shù)據(jù)得到的病疫木信息提取精度有所降低，為80.2%。

表1 分類精度評價表

分析錯分地類的地理位置分布，可看出錯分的地類主要分布在地類的交界處和地物破碎處。在山坡地帶，病疫木和健康木交錯生長，無明顯空間界限，且病疫木感染無規(guī)律，分布凌亂，不易分類；在平緩地帶，多種人類活動導(dǎo)致土地的種類更加多樣，地類混合的情況更為復(fù)雜，地類的提取精度也受到一定影響。但總體來說，面向?qū)ο蟮腃ART決策樹分類提取松材線蟲病疫木具有較高可行性。

3 結(jié)論與討論

3.1 結(jié)論

在衛(wèi)星遙感加無人機航空遙感加地面移動調(diào)查中，調(diào)查人員發(fā)現(xiàn)，松材線蟲病災(zāi)害的區(qū)域大多數(shù)集中于山頂部分以及山林深處，無人機在調(diào)查驗證過程中起到了關(guān)鍵作用。調(diào)查人員實地前往感染區(qū)域利用無人機對山頂樹木進行飛行作業(yè)，根據(jù)實地調(diào)查驗證結(jié)果以及無人機飛行得到的驗證資料，對影像區(qū)域重新進行面向?qū)ο蟮腃ART決策樹分類，分類結(jié)果如圖7所示。

表2 精度對比表

從表2可以得出，利用無人機遙感影像和地面移動調(diào)查對衛(wèi)星遙感監(jiān)測的補充，能有效提高松材線蟲病疫木提取的精度。遙感衛(wèi)星、無人機和人工輔助3種技術(shù)相結(jié)合，實現(xiàn)了區(qū)域范圍內(nèi)松材線蟲病病疫木的快速識別，準(zhǔn)確率提升到90.65%，有利于防控指揮部門及時掌握各疫點疫情發(fā)生情況和風(fēng)險程度，有效提高了松材線蟲病防控工作和疫情評估的科學(xué)性，為進一步防止松材線蟲病的擴散和蔓延提供科學(xué)的技術(shù)支持。

3.2 討論

本研究利用0.8m空間分辨率的BJ-2號遙感影像，通過面向?qū)ο蟮腃ART決策樹分類提取松材線蟲病病疫木區(qū)域，后利用無人機和地面移動設(shè)備進行調(diào)查驗證，將三者結(jié)合得到的數(shù)據(jù)進行綜合分析，普查松材線蟲病的時間縮短，所需人力物力等資源減少，預(yù)警時效性增強，并且可以實時監(jiān)測預(yù)警松材線蟲病蟲害的發(fā)生地點、時間，在此基礎(chǔ)上得出松材線蟲病實際發(fā)生區(qū)域圖，再將所有結(jié)果進行統(tǒng)計分析，得出松材線蟲病普查預(yù)警相關(guān)材料和報告。