閆利，李青山，王嫣然，葉志云

(武漢大學(xué) 測(cè)繪學(xué)院，武漢 430079)

0 引言

中波紅外譜段(3～5 μm)與熱紅外譜段(8～12 μm)一直是紅外技術(shù)發(fā)展應(yīng)用的重點(diǎn)。面陣型紅外焦平面探測(cè)器技術(shù)出現(xiàn)后，第二代紅外成像技術(shù)已將工作波段由長波擴(kuò)展至中波和短波紅外。目前，紅外成像技術(shù)已全部覆蓋短波紅外至長波紅外譜段[1]。相較于可見光譜段，紅外波段對(duì)煙霧穿透能力更強(qiáng)，可全天候成像，其中中波紅外在制導(dǎo)頭、自動(dòng)跟蹤系統(tǒng)以及對(duì)空中飛行目標(biāo)探測(cè)領(lǐng)域作用距離更遠(yuǎn)[2]。

中波紅外譜段在對(duì)地觀測(cè)中的優(yōu)勢(shì)可廣泛應(yīng)用于目標(biāo)提取[3]、溫度反演[4]、圖像融合[5]、工業(yè)排放氣體監(jiān)測(cè)[6]、火災(zāi)監(jiān)測(cè)[7]、武器制導(dǎo)[8]等應(yīng)用領(lǐng)域。在遙感觀測(cè)領(lǐng)域，國外已有包含中波紅外譜段的機(jī)載高光譜成像系統(tǒng)[9]?；谥胁t外譜段，已有海岸線分割[10]以及植物類別提取[11]等地物分類相關(guān)研究。由圖1可知，我國在中波紅外波段的對(duì)地觀測(cè)上空間分辨率較低，與國外存在較大差距。同時(shí)由于政策原因，國內(nèi)可采購中波紅外相機(jī)最大面陣大小為640像素×512像素。故目前國內(nèi)尚未展開中波紅外在遙感影像分類中的研究。開展中波紅外遙感相關(guān)研究可以有效利用頻譜資源，填補(bǔ)相應(yīng)領(lǐng)域的空白。

圖1 紅外波段衛(wèi)星對(duì)地觀測(cè)現(xiàn)狀示意圖

本文的項(xiàng)目背景是國內(nèi)首套機(jī)載中波紅外成像系統(tǒng)的研制，該系統(tǒng)同時(shí)裝配有中波紅外和可見光光學(xué)相機(jī)。本文旨在探索地物在中波紅外波段的輻射特性，進(jìn)而融合地物的中波紅外與可見光的特征，提高遙感地物分類的精度，探索中波紅外譜段在對(duì)地觀測(cè)領(lǐng)域的潛在價(jià)值。

1 研究區(qū)及數(shù)據(jù)源概況

1.1 研究區(qū)概況

研究區(qū)位于湖北省荊門市漳河機(jī)場(chǎng)。荊門為湖北省地級(jí)市，位于湖北省中部，漢江中下游，介于111°51′E～113°29′E，30°32′N～31°36′N之間。荊門市屬于亞熱帶季風(fēng)氣候，其植被主要以常綠闊葉、落葉闊葉混交林為主，松樹、杉木、櫟樹分布普遍。建筑物主要由零散村鎮(zhèn)中的低層民居以及城鎮(zhèn)居民區(qū)和低層廠房構(gòu)成。影像中所包含道路類型很豐富，有柏油馬路、水泥路面以及飛機(jī)跑道。同時(shí)，影像中還包含有林地、草地、耕地等植被類型和湖泊、池塘等水體類型。

1.2 數(shù)據(jù)源及預(yù)處理

研究區(qū)航拍影像為2017年2月26日早9時(shí)至10時(shí)于荊門漳河機(jī)場(chǎng)通過運(yùn)-5飛機(jī)搭載自主研發(fā)航拍系統(tǒng)得到。其可見光影像由哈蘇H3D相機(jī)拍攝，影像大小為7 216像素×5 412像素，分辨率為0.136 6 m，中紅外影像由Onca-MWIR-MCT640相機(jī)拍攝，其接收光譜范圍為3.7～4.8 μm，影像大小為640像素×512像素，分辨率為0.625 m。

航拍得到的中紅外影像為小像幅影像，與光學(xué)相機(jī)的像幅相差很大。因此，中波紅外相機(jī)采用面陣擺掃成像的方式，以實(shí)現(xiàn)寬幅機(jī)載中波紅外成像。利用定位定姿系統(tǒng)(positioning and orientation system，POS)以及地面實(shí)測(cè)控制點(diǎn)實(shí)施POS輔助空中三角測(cè)量，制作可見光及中紅外正射影像(digital orthophoto map，DOM)。然后將可見光與中紅外影像生成的正射影像進(jìn)行配準(zhǔn)，裁剪二者的重疊區(qū)域得到最終使用的研究數(shù)據(jù)(圖2)，其中圖2(a)為可見光(visible，VIS)三波段彩色顯示，圖2(b)為中波紅外譜段(mid-wave infrared，MID)灰度顯示影像。

圖2 研究數(shù)據(jù)

2 中紅外波段特性分析

2.1 中波紅外光譜輻射特性

中波紅外譜段(3～5 μm)與熱紅外波段(8～12 μm)同屬于紅外譜段，但其光譜輻射特性不同于可見光與熱紅外譜段，既包含地面反射輻射，也包含地面物體的發(fā)射輻射能量，可實(shí)現(xiàn)全天候成像。圖3表明[12]，不同溫度目標(biāo)在中波紅外波段具有更高的熱對(duì)比度，更容易區(qū)分不同溫度目標(biāo)，同時(shí)在600～1 000 K物體輻射功率峰值對(duì)應(yīng)波長為中波，具有優(yōu)于熱紅外的高溫物體監(jiān)測(cè)能力。此外，中波紅外波長較長，不易受到霧霾的散射，能有效突破低能見度氣候?qū)b感影像獲取時(shí)效性的限制。

圖3 不同溫度目標(biāo)對(duì)比示意圖

2.2 研究影像中波紅外特性

光學(xué)遙感傳感器的照明源主要是太陽輻射，因此高分辨率遙感影像中存在很多由建筑、橋梁、樹木等地物的遮擋形成的陰影區(qū)域。這些陰影區(qū)域雖然對(duì)于影像的信息提取及密集匹配等處理有很大的干擾[13]，但研究表明[14]，同一地物類型非陰影區(qū)域與陰影區(qū)域輻射值之比有與波長成正比的物理特性。

結(jié)合遙感影像非陰影區(qū)域與陰影區(qū)域接收到的光照輻射能量比值的特性與中波紅外譜段包含地物發(fā)射輻射能量的特性，分析研究影像不同地物非陰影區(qū)域與陰影區(qū)域的光譜特性曲線。如圖 4所示，中波紅外波段顯示出不同于可見光譜段的輻射特性。波譜特性曲線(1)表明，裸露土地與水泥道路的波譜特性曲線在可見光譜段呈上升趨勢(shì)而后在中波紅外波段下降，瀝青道路則呈現(xiàn)相反的特征。由波譜特性曲線(2)可知，中波紅外波段不同地物的亮度比值(非陰影區(qū)/陰影區(qū))不符合隨波長增大而增大的趨勢(shì)。其原因?yàn)橹胁t外波段傳感器接收到的輻射能量不只包括地物反射的太陽直射輻射與大氣漫反射，還包含有地物的發(fā)射輻射能量(圖 5)，傳感器在中紅外譜段接收到的太陽直射輻射能量在總輻射能量中占比低于可見光譜段，故地物的亮度比值(非陰影區(qū)/陰影區(qū))不遵循可見光譜段的規(guī)律。

圖4 研究區(qū)影像中波紅外特性

圖5 可見光譜段與中波紅外譜段輻射能量構(gòu)成

3 研究方法

3.1 多尺度分割

影像分割是面向?qū)ο蠓诸惖幕A(chǔ)。通過圖像分割可以得到同質(zhì)像元組成的多邊形對(duì)象。多尺度分割[15](multiresolution segmentation)算法連續(xù)地合并像元或現(xiàn)有的影像對(duì)象。該過程首先基于像素層，根據(jù)同質(zhì)性標(biāo)準(zhǔn)合并像元，然后基于合成的對(duì)象結(jié)合設(shè)置的尺度參數(shù)繼續(xù)合并影像對(duì)象，循環(huán)此步驟直至同質(zhì)性指標(biāo)大于設(shè)定的尺度參數(shù)，分割即停止。多尺度分割過程識(shí)別一個(gè)像元或現(xiàn)有的影像對(duì)象并根據(jù)相關(guān)的同質(zhì)性標(biāo)準(zhǔn)與它們相鄰的對(duì)象合并，是一種自下而上的分割。

3.2 隨機(jī)森林分類器及特征評(píng)估

隨機(jī)森林(random forests)算法是Leo Breiman繼Bagging算法之后，提出的另一種利用多個(gè)樹學(xué)習(xí)器進(jìn)行分類和回歸的組合算法[16-17]。該算法利用bootstrap重抽樣方法從原始樣本中抽取多個(gè)重抽樣樣本集，對(duì)每個(gè)bootstrap樣本進(jìn)行決策樹建模得到?jīng)Q策樹分類模型集：

{h(X，Θk)，k=1，2，…}

(1)

式中：X為給定自變量，本文中具體為所選訓(xùn)練樣本對(duì)應(yīng)的特征向量集；參數(shù)集{Θk}是獨(dú)立同分布的隨機(jī)向量；k為隨機(jī)森林中決策樹棵樹，本文實(shí)驗(yàn)中k為100。

隨機(jī)森林分類器可通過OOB平均下降精度、Gini指數(shù)來度量特征重要性。

OOB平均下降精度衡量特征重要性FIM(feature importance measures)：

FIM=∑(errOOB2-errOOB1)/Ntree

(2)

式中：errOOB1為每棵決策樹的袋外數(shù)據(jù)誤差；errOOB2為對(duì)袋外數(shù)據(jù)OOB所有樣本的特征X加入噪聲干擾，再次計(jì)算得到的袋外數(shù)據(jù)誤差。

(3)

式中：K表示有K個(gè)類別；|K|為類別個(gè)數(shù)；k′表示除k之外的類別；Pmk′表示節(jié)點(diǎn)m中類別k′所占的比例；Pmk表示節(jié)點(diǎn)m中類別k所占的比例。

特征Xj在節(jié)點(diǎn)m的重要性，即節(jié)點(diǎn)m分枝前后的Gini指數(shù)變化量為：

(4)

式中：GIl和GIr分別表示分枝后2個(gè)新節(jié)點(diǎn)的Gini指數(shù)。如果，特征Xj在決策樹i中出現(xiàn)的節(jié)點(diǎn)在集合M中，那么Xj在第i棵樹的重要性為：

(5)

假設(shè)RF中共有n棵樹，那么

(6)

最后，所有求得的重要性評(píng)分歸一化：

(7)

3.3 特征空間分析與優(yōu)化

通過分割形成的影像對(duì)象在面向?qū)ο蟮姆诸惙椒ㄖ袨樾碌男畔⑤d體，其包含有光譜、形狀、紋理、層次等特性，這些特性稱為特征(features)。充分利用影像對(duì)象的這些特征信息能夠產(chǎn)生更加精確的分類結(jié)果，對(duì)初始特征空間進(jìn)行合理的特征選擇，將經(jīng)過優(yōu)選之后的特征空間作為分類特征，有助于提高分類器的分類效率及精度[18]

表1描述了本文所選擇的初始特征空間，涵蓋波段灰度值、亮度、標(biāo)準(zhǔn)差等光譜特征，長寬比、形狀指數(shù)等幾何特征以及GLCM紋理特征。

表1 初始特征空間

實(shí)驗(yàn)中基于分割形成的對(duì)象選取分析樣本，提取可見光影像與可見光+中紅外融合四波段影像特征向量集。其中HIS顏色模型特征選取可見光與中紅外四波段中任意三波段計(jì)算各HIS分量，其余特征單獨(dú)提取中紅外波段影像特征，構(gòu)建可見光+中紅外融合影像特征向量集。分別對(duì)可見光影像特征集與融合影像特征集通過OOB下降精度與Gini指數(shù)排序特征重要性，通過驗(yàn)證樣本集精度確定各自的最優(yōu)特征向量集，實(shí)施隨機(jī)森林分類。

4 實(shí)驗(yàn)與分析

4.1 實(shí)驗(yàn)流程

本文實(shí)施例中，基于多尺度分割與隨機(jī)森林分類器實(shí)現(xiàn)面向?qū)ο蟮匚锓诸?，融合了地物的可見光與中波紅外特征，并且評(píng)估了光譜、形狀、紋理等特征在分類中的重要程度，定量分析了融合中波紅外波段后的特征空間。如圖6所示，包括以下步驟：

步驟1，數(shù)據(jù)準(zhǔn)備。包括利用POS系統(tǒng)及地面實(shí)測(cè)控制點(diǎn)制作可見光與中紅外正射影像，將生成的正射影像分為可見光三波段影像與可見光與中紅外融合四波段影像，2組數(shù)據(jù)進(jìn)行對(duì)照實(shí)驗(yàn)。

步驟2，影像分割。基于多尺度分割算法分割影像得到對(duì)象圖層，實(shí)驗(yàn)中2組數(shù)據(jù)均以可見光三波段影像等權(quán)重進(jìn)行多尺度分割，中紅外波段不參與分割，以保證特征分析選取樣本一致。形狀異質(zhì)性指標(biāo)權(quán)重為0.1，平滑與緊密2個(gè)子異質(zhì)性指標(biāo)權(quán)重均為0.5，分割參數(shù)為45。

步驟3，初始特征提取。基于分割形成的對(duì)象圖層選取分析樣本及初始特征向量集。將選取的樣本分為訓(xùn)練樣本與驗(yàn)證樣本，分別用于模型訓(xùn)練與精度分析。分類體系根據(jù)研究區(qū)影像內(nèi)容，確定為裸露土地、建筑物、林地、道路交通、草地和水體6種地物類型，各地物均選取1 000個(gè)訓(xùn)練樣本，500個(gè)驗(yàn)證樣本。根據(jù)初始特征空間分別提取可見光與融合影像的訓(xùn)練樣本與測(cè)試樣本的初始特征向量集?？梢姽庥跋衽c融合影像初始特征向量維度分別為88、120。

圖6 實(shí)驗(yàn)流程圖

步驟4，優(yōu)化特征空間。分析步驟3提取的訓(xùn)練樣本原始特征向量集，利用隨機(jī)森林OOB平均下降精度與Gini指數(shù)衡量特征向量重要性，達(dá)到優(yōu)化特征空間的目的。

步驟5，隨機(jī)森林分類。基于步驟2得到的分割結(jié)果與步驟3確定的優(yōu)化特征空間利用隨機(jī)森林分類器進(jìn)行分類并分析分類精度，對(duì)比可見光影像分類結(jié)果與融合影像分類結(jié)果，分析中波紅外特征在分類中的作用。

4.2 優(yōu)化特征空間

利用訓(xùn)練樣本初始特征向量集，訓(xùn)練隨機(jī)森林分類器，并根據(jù)隨機(jī)森林OOB平均下降精度與Gini指數(shù)衡量特征向量重要性，可見光+中紅外融合影像初始特征向量集重要性衡量如圖7所示，分別為OOB精度下降示意圖與Gini指數(shù)示意圖。根據(jù)兩指數(shù)特征重要性排序如表2所示，涉及中波紅外波段的特征由黑體顯示，針對(duì)中紅外特征，最有效特征為中紅外與可見光其中兩波段組合HIS空間各分量特征，其次為GLCM紋理信息。

結(jié)合2種指數(shù)，通過排序特征變量重要性指標(biāo)優(yōu)化選擇特征空間，訓(xùn)練隨機(jī)森林分類器，以驗(yàn)證樣本分類總精度選擇優(yōu)化特征空間維度，選取總精度最高的維度，如表3所示，實(shí)驗(yàn)中可見光影像優(yōu)化特征空間維度為35，融合影像優(yōu)化特征空間維度為40。

圖7 融合影像特征重要性示意圖

表2 融合影像特征重要性排序表

表3 不同維度優(yōu)化特征空間分類結(jié)果總精度表 %

4.3 分類結(jié)果分析與精度評(píng)價(jià)

基于步驟2分割結(jié)果與步驟4優(yōu)化特征空間，分別對(duì)可見光影像與融合影像實(shí)現(xiàn)基于隨機(jī)森林的面向?qū)ο蠓诸?，分類結(jié)果細(xì)節(jié)如圖8所示。以步驟3選取的驗(yàn)證樣本評(píng)估可見光影像單獨(dú)分類與融合影像分類精度差異，分類混淆矩陣見表4、表5。分析分類結(jié)果細(xì)節(jié)圖8及易錯(cuò)分地物波譜特性曲線圖9可知：

圖9 易錯(cuò)分地物波譜特性曲線

區(qū)域(1)(2)中植被陰影區(qū)、瀝青道路在可見光影像中易被錯(cuò)分為水體，引入中波紅外波段可以改善此類錯(cuò)分現(xiàn)象。結(jié)合圖8中VIS、MID影像與圖 9(a)分析：三類地物在可見光譜段亮度值均較低，色彩暗淡。但瀝青道路與植被陰影區(qū)在中波紅外譜段相較于紅光波段明顯上升，水體則沒有此趨勢(shì)。因此引入中波紅外譜段可避免植被陰影區(qū)和瀝青道路被錯(cuò)分為水體的現(xiàn)象，提高道路交通類別的制圖精度。

區(qū)域(3)中磚瓦屋頂在可見光影像中被錯(cuò)分為裸露土地。結(jié)合圖8中VIS、MID影像與圖9(b)分析：磚瓦屋頂在可見光譜段與裸露土地顏色近似，但在中波紅外譜段，磚瓦屋頂向陽面亮度值明顯高于裸露土地以及背陰面磚瓦屋頂，故引入中波紅外譜段有助于此類房屋的提取。同時(shí)，對(duì)比區(qū)域(1)VIS+MID分類結(jié)果與VIS分類結(jié)果，水泥材質(zhì)屋頂在引入中波紅外后分類結(jié)果更加完整。因此，引入中波紅外譜段有助于提取磚瓦材質(zhì)屋頂以及水泥材質(zhì)屋頂?shù)耐暾裕岣呓ㄖ镱悇e的制圖精度。

由表2及分類結(jié)果混淆矩陣(表4、表5)可知：

1)可見光影像地物分類總精度最高為87.87%，融合影像地物分類總精度最高為89.63%，引入中波紅外波段后地物分類總精度可提高1.76%；

2)人工地物(建筑物、道路交通)制圖精度提升了2.8%，非人工地物制圖精度僅提升了1.25%，故中波紅外波段對(duì)于人工地物的分類效果優(yōu)于非人工地物類型。由表6可知，建筑物的精度提升最多，其制圖精度與用戶精度分別提高3.00%、3.44%；其次為道路交通，其制圖精度與用戶精度分別提升了2.60%、2.70%。

表4 可見光影像分類結(jié)果混淆矩陣 %

表5 可見光&中紅外融合影像分類結(jié)果混淆矩陣 %

表6 融合影像分類地物精度變化表 %

5 結(jié)束語

本文使用自主研發(fā)的機(jī)載中波紅外成像系統(tǒng)，航拍得到荊門市漳河機(jī)場(chǎng)區(qū)域影像資料，生成可見光及中波紅外正射影像后，基于隨機(jī)森林分類器研究融合影像特征空間及中波紅外影像對(duì)地物分類精度的提升效果。根據(jù)目前的實(shí)驗(yàn)結(jié)果，中波紅外波段更建議應(yīng)用于城市區(qū)域遙感相關(guān)應(yīng)用領(lǐng)域。

研究結(jié)果表明：

①針對(duì)中紅外特征，最有效特征為其與可見光其中兩波段組合HIS空間各分量特征，其次為GLCM紋理信息。

②中波紅外可以穩(wěn)定得提高地物分類的總精度，可見光影像與可見光+中紅外影像的最高分類精度分別為87.87%、89.63%，引入中波紅外波段提高了1.76%。

③中波紅外波段對(duì)于人工地物的分類效果優(yōu)于非人工地物類型。引入中波紅外波段后，人工地物制圖精度提升了2.8%，非人工地物制圖精度僅提升了1.25%。其中建筑物的精度提升最多，其制圖精度與用戶精度分別提高3.00%、3.44%；其次為道路交通，其制圖精度與用戶精度分別提升了2.60%、2.70%。