張昊,焦瑞莉,喬聰聰,霍娟,宗雪梅

(1.北京信息科技大學(xué) 信息與通信工程學(xué)院,北京 100101;2.中國科學(xué)院大氣物理研究所,北京 100029)

0 引言

利用云和云陰影區(qū)域與其他區(qū)域的差別實現(xiàn)氣溶膠參數(shù)的反演是一種可行的方法[1],并且由于云陰影檢測通常由云檢測實現(xiàn),因此在遙感影像中準(zhǔn)確的云檢測是后續(xù)氣溶膠研究的重要前提。

近年來,國內(nèi)外在遙感影像的云檢測方面進行了大量的研究,并提出多種不同的方法。云檢測的方法可分為兩大類:光譜閾值方法和深度學(xué)習(xí)方法。

光譜閾值方法利用衛(wèi)星遙感影像的若干個光學(xué)波段的反射率進行閾值判別,從而進行遙感影像的云檢測。光譜閾值方法可分為單時相方法和多時相方法。單時相方法通過對單幅遙感影像進行光譜反射率的閾值判別得到云檢測結(jié)果。Luo等[2]提出使用閾值決策矩陣來識別MODIS衛(wèi)星影像中的云。Zhu等[3]提出Fmask方法,對Landsat衛(wèi)星影像采用光譜閾值和概率閾值相結(jié)合的兩步閾值法進行云檢測,并結(jié)合卷云波段和多種輔助數(shù)據(jù)提升檢測準(zhǔn)確率[4-6]。多時相方法是利用多幅不同時間、相同地點的衛(wèi)星影像組成時間序列,并結(jié)合光譜反射率閾值進行判別得到云檢測結(jié)果。文獻(xiàn)[7-8]分別在Fmask方法的云檢測結(jié)果基礎(chǔ)上,利用時間序列進一步檢測。然而,單時相光譜閾值方法主要通過光譜特征進行檢測,因此針對某些和云具有相似光譜特征的區(qū)域效果較差;多時相光譜閾值方法雖然通過時間序列在一定程度上提高了云檢測的準(zhǔn)確率,但是其對數(shù)據(jù)量需求較大,運算效率較低。

深度學(xué)習(xí)方法利用龐大的數(shù)據(jù)集對云檢測模型進行訓(xùn)練,然后通過訓(xùn)練完成的模型對被測影像的每一個像素進行預(yù)測分類。Wieland等[9]使用U-Net模型進行遙感影像中的云檢測。Li等[10]提出了一種以多個殘差單元為編碼器,以相同結(jié)構(gòu)的多個殘差單元和反卷積為解碼器的多尺度特征融合云檢測方法。Zhang等[11]提出了一種以多個殘差單元和小波變換模塊為編碼器,以ASPP金字塔和空間注意力機制為解碼器的多尺度特征融合模型。Mehdi等[12]提出了一種以不同尺寸的卷積核組成的并聯(lián)結(jié)構(gòu)為編碼器,以反卷積為解碼器的多尺度特征融合模型。Qu等[13]提出了一種以條狀殘差網(wǎng)絡(luò)為編碼器,以雙線性插值和雙注意力機制為解碼器的多尺度特征融合模型。然而,深度學(xué)習(xí)方法在具有良好檢測結(jié)果的同時,需要龐大的數(shù)據(jù)量和較高的硬件需求,且對訓(xùn)練參數(shù)的初始化較為敏感。

多特征聯(lián)合(multi-feature combined,MFC)是一種針對單幅遙感影像,并結(jié)合了光譜特征、幾何特征和紋理特征,應(yīng)用在高分一號的云檢測方法[14]。由于MFC方法在高分一號數(shù)據(jù)上表現(xiàn)出良好的性能,因此本文采用波段范圍更加豐富的Landsat 8數(shù)據(jù)進行MFC方法的性能評價,并提出兩點改進。一是針對幾何特征與云相似的非云區(qū)域,利用灰度共生矩陣提取紋理特征,提高云識別正確率。二是針對冰雪檢測的局限性,引入一種光譜特征對冰雪區(qū)域進行檢測和濾除。

1 理論基礎(chǔ)

1.1 光譜特征和引導(dǎo)濾波生成初步云檢測結(jié)果

MFC方法首先通過3個光譜特征的閾值判別生成一個包含大部分厚云的云檢測結(jié)果,在最大限度上減小非云像素引起的錯分誤差。

第1個光譜特征是霧霾最優(yōu)變換(haze optimized transformation,HOT)。其原理是在藍(lán)色波段和紅色波段中,云和大部分地表以及薄霧有著不同的光譜特性,前者的HOT值通常大于后者的HOT值。

第2個光譜特征是可見光波段比值(visible band ratio,VBR)。VBR能夠排除具有顯著顏色特征的非云像素。云在可見光波段通常為灰色或白色,其VBR值接近1。

第3個光譜特征是紅色波段反射率,其原理是云在紅色波段的反射率較高。

MFC方法分別將3個光譜特征的閾值設(shè)置為0.13、0.7和0.07。如果像素同時滿足上述3個閾值條件,則將其檢測為云。

在得到包含大部分厚云的云檢測結(jié)果之后,MFC方法認(rèn)為薄云通常分布在厚云周圍,即在云邊界周圍存在一個從厚云到薄云的過渡區(qū)域。因此在光譜特征檢測的基礎(chǔ)上,采用引導(dǎo)濾波器檢測每個云邊界遺漏的薄云。

引導(dǎo)濾波器包括引導(dǎo)圖像、輸入圖像和輸出圖像。在MFC方法中,引導(dǎo)圖像是藍(lán)色波段、綠色波段和紅色波段的合成影像,輸入圖像是光譜特征閾值判別得到的二值云檢測圖像,輸出圖像是引導(dǎo)濾波后的二值云檢測圖像。

引導(dǎo)濾波器描述了引導(dǎo)圖像和輸出圖像之間的局部線性關(guān)系,即輸出圖像是引導(dǎo)圖像以某個像素為中心的窗口的線性變換。對于以某個像素為中心的窗口,線性變換的系數(shù)包括ak和bk。其中ak與引導(dǎo)圖像在此窗口內(nèi)的所有像素的總數(shù)、均值、方差,以及正則化參數(shù)有關(guān);bk與引導(dǎo)圖像和輸入圖像在此窗口內(nèi)的所有像素的均值,以及ak有關(guān)。

對于任意像素i,其包含在多個不同的窗口,因此像素i的輸出應(yīng)考慮到包含其本身的所有窗口,即系數(shù)a和b分別為包含其本身的所有窗口的ak和bk的平均值。

MFC方法分別針對水區(qū)域和陸地區(qū)域設(shè)置了不同的云邊界檢測方法。首先采用歸一化植被指數(shù)(normalized difference vegetation index,NDVI)和近紅外波段反射率將水區(qū)域和陸地區(qū)域進行劃分,原理在于后者的NDVI值通常高于前者的NDVI值,且前者在近紅外波段的反射率較低。MFC方法使用了兩組不同的閾值判別進行水區(qū)域的檢測。第1組將NDVI閾值設(shè)置為0.15,近紅外波段閾值設(shè)置為0.2;第2組將NDVI閾值設(shè)置為0.2,近紅外波段閾值設(shè)置為0.15。如果滿足其中任意一組閾值條件,則檢測為水區(qū)域,否則為陸地區(qū)域。

對于水區(qū)域,MFC方法采用大津閾值法對引導(dǎo)濾波器輸出的灰度圖像進行閾值分割生成二值云檢測圖像,閾值設(shè)置為0.12;對于陸地區(qū)域,在對引導(dǎo)濾波器輸出圖像進行閾值分割的基礎(chǔ)上,再次使用HOT進行閾值分割,以防止云周圍非云像素的干擾,閾值設(shè)置為0.08。

然而,初步云檢測結(jié)果可能包含一些非云區(qū)域,如冰雪、建筑物等,它們與云區(qū)域具有相似的光譜特征,易與云區(qū)域混淆。因此在初步云檢測結(jié)果中對這些非云區(qū)域進行濾除是十分必要的。

1.2 濾除非云區(qū)域

在濾除非云區(qū)域之前,MFC方法首先在初步云檢測二值圖像中進行8鄰域連通,形成若干個連通域,以連通域為單位提取多種特征并逐一進行判別。

MFC方法通過利用連通域提取分形維數(shù)指數(shù)(fractal dimension index,FRAC)和長寬比(length width ratio,LWR)這兩種幾何特征對海岸線、道路、建筑物等具有幾何形狀的區(qū)域進行濾除。FRAC由連通域的周長和面積確定,反映了連通域形狀的復(fù)雜性。LWR由連通域最小外接矩形的長和寬確定,反映了連通域最小外接矩形的長與寬之間的關(guān)系。相比于云區(qū)域,這些非云區(qū)域的FRAC值和LWR值較大。

然而,其中一些非云區(qū)域的FRAC和LWR和云區(qū)域較為相似,僅通過幾何特征無法將它們與云進行區(qū)分。因此,本文在幾何特征的基礎(chǔ)上采用灰度共生矩陣進一步提取連通域的紋理特征進行判別。

灰度共生矩陣反映了圖像灰度關(guān)于方向、距離、變化幅度的綜合信息,它是分析圖像局部模式結(jié)構(gòu)及其排列規(guī)則的基礎(chǔ),其統(tǒng)計了圖像在某一方向上相距一定距離的兩個灰度值同時出現(xiàn)的概率。

作為紋理分析的特征量,往往不是直接應(yīng)用灰度共生矩陣,而是在灰度共生矩陣的基礎(chǔ)上再提取紋理特征量。均值和[15]是紋理特征量之一,是區(qū)域內(nèi)像素平均灰度值的度量,反映了區(qū)域的明暗深淺。對于具有較高灰度值的區(qū)域,均值和較大;對于具有較低灰度值的區(qū)域,均值和較小。

在可見光波段,云具有較高反射率,即較高灰度值,均值和較大;大部分地表具有較低反射率,即較低灰度值,均值和較小。因此選取藍(lán)色波段用于計算灰度共生矩陣的均值和,提取紋理特征。

考慮到方法的效率以及灰度共生矩陣的計算量,本文從4個方向中選取45°方向,相對距離設(shè)置為1,灰度級數(shù)設(shè)置為16,窗口大小設(shè)置為11。為了避免個別異常值的影響,采用以連通域為單位的判別方式,即統(tǒng)計每個連通域中灰度共生矩陣均值和大于閾值T的像素總數(shù),并計算其與所在連通域的面積的比值,如果比值小于0.5,則判別為非云區(qū)域并進行濾除,否則判別為云區(qū)域。對于閾值T的選擇,本文選取了239個誤判為云的非云連通域和276個云連通域,并計算在不同閾值T的選取下分類的正確率。首先在區(qū)間[2,3]之間以0.1為間隔對T進行取值,計算分類正確率,結(jié)果如圖1所示?？梢钥闯?閾值在區(qū)間[2.5,2.6]時,云分類正確率所受影響較小,且非云分類正確率較高。為了得到更加精確的閾值,在區(qū)間[2.5,2.6]內(nèi)以0.01為間隔對閾值T進行取值,計算分類正確率,結(jié)果如圖2所示。可以看出,若要同時兼顧二者的分類正確率,則將閾值T設(shè)置為2.56時最佳。

圖1 閾值區(qū)間為[2,3]的分類正確率曲線

圖2 閾值區(qū)間為[2.5,2.6]的分類正確率曲線

對于冰雪區(qū)域,MFC方法選取了圖像尺寸為100像素×100像素的84個樣本,計算它們的局部二值模式(local binary patterns,LBP)紋理直方圖,并依據(jù)卡方距離進行判斷。分析認(rèn)為此方法具有一定的局限性。一是由于MFC方法針對面積小于40 000個像素點的連通域進行判別,而遙感影像中有可能存在大面積的冰雪區(qū)域,因此可能會造成漏判。二是檢測結(jié)果取決于樣本的選擇,如果樣本選擇不佳,則可能影響最終的檢測結(jié)果。因此,本文以像素點為單位,引入一種光譜特征對冰雪區(qū)域進行判別并濾除。

圖3是云和雪從可見光波段到短波紅外波段的反射率折線圖[16],其中藍(lán)色波段的波長范圍為0.45～0.52 μm,短波紅外波段的波長范圍為1.56～1.66 μm。由圖3可知,在藍(lán)色波段,云和雪的反射率較高,通常大于0.6;在短波紅外波段,云和雪的反射率相對于可見光波段有所降低。不同的是,云的反射率在0.6浮動,下降量較小,而雪的反射率小于0.2,下降量較大。因此,通過藍(lán)色波段和短波紅外波段進行光譜特征提取,特征表達(dá)式見式(1)。

圖3 云和冰雪的反射率曲線

SI=2×B6-B2

(1)

式中:B6表示短波紅外波段反射率;B2表示藍(lán)色波段反射率。如果像素的SI值小于閾值0,則將其檢測為冰雪,否則將其檢測為云。

2 實驗與分析

2.1 實驗方案

實驗所采用的數(shù)據(jù)是Landsat 8 OLI衛(wèi)星遙感數(shù)據(jù)。為了驗證改進MFC方法在不同地表類型下的性能,本文選取了涵蓋森林、冰雪、水、城市、瘠地等不同地表類型的遙感影像進行性能評價。

用于計算性能參數(shù)的參照結(jié)果采用LandsatLook 8-bit質(zhì)量評估圖像。LandsatLook 8-bit質(zhì)量評估圖像采用不同顏色標(biāo)注了遙感影像中每個像素的類別。因此,通過對LandsatLook 8-bit質(zhì)量評估圖像中的云像素進行標(biāo)記得到云參照結(jié)果,將其作為真實標(biāo)記,并分別與原始MFC方法和改進MFC方法的云檢測結(jié)果進行性能參數(shù)的計算。

實驗通過計算總體精度(overall accuracy,OA)、用戶精度(user accuracy,UA)、生產(chǎn)精度(producer accuracy,PA)這3個性能參數(shù)對改進MFC方法進行性能評價,并與原始MFC方法進行對比。將所有像素劃分為云類或非云類,通過衡量云檢測結(jié)果和云參照結(jié)果中像素類別的一致性和差異性進行精度計算,計算方法如式(2)、式(3)、式(4)所示。

(2)

(3)

(4)

式中:MP和MN分別表示參照結(jié)果中云的像素數(shù)量和參照結(jié)果中非云的像素數(shù)量;MTP表示參照結(jié)果中標(biāo)記為云且檢測結(jié)果標(biāo)記為云的像素數(shù)量;MTN表示參照結(jié)果標(biāo)記為非云且檢測結(jié)果標(biāo)記為非云的像素數(shù)量;MFP表示參照結(jié)果標(biāo)記為非云且檢測結(jié)果標(biāo)記為云的像素數(shù)量;MFN表示參照結(jié)果標(biāo)記為云且檢測結(jié)果標(biāo)記為非云的像素數(shù)量。

總體精度OA衡量了總體分類的正確性;生產(chǎn)精度PA衡量了遺漏分類的程度;用戶精度UA衡量了錯誤分類的程度。

圖4是實驗流程圖,總共分為4個步驟。第1步是數(shù)據(jù)預(yù)處理,將遙感影像的灰度值轉(zhuǎn)換為反射率;第2步是結(jié)合光譜閾值分割和引導(dǎo)濾波器生成初始云檢測結(jié)果;第3步是利用幾何特征、灰度共生矩陣的紋理特征,以及本文引入的光譜特征濾除初始云檢測結(jié)果中的建筑物、冰雪等非云區(qū)域;第4步是進行實驗結(jié)果分析與對比。

圖4 實驗流程圖

2.2 實驗結(jié)果與分析

圖5、圖6、圖7列出了3組同時含有云和建筑物、道路等非云區(qū)域的不同地點的結(jié)果,從左至右依次為遙感影像、原始MFC方法的云檢測結(jié)果以及改進MFC方法的云檢測結(jié)果,其中白色表示云區(qū)域,黑色表示非云區(qū)域。每組圖像中,圖(a)是區(qū)域尺寸為1 024像素×1 024像素的影像,其中紅色矩形標(biāo)記了區(qū)域尺寸為256像素×256像素的部分區(qū)域。圖(b)是圖(a)中紅色矩形區(qū)域的放大圖像,其中黃色矩形標(biāo)記了其中部分非云區(qū)域。

圖5 含有建筑物等非云區(qū)域的遙感影像云檢測結(jié)果一

圖6 含有建筑物等非云區(qū)域的遙感影像云檢測結(jié)果二

圖7 含有建筑物等非云區(qū)域的遙感影像云檢測結(jié)果三

從上述圖中可以看出,原始MFC方法的云檢測結(jié)果中仍然存在著一些非云區(qū)域,說明僅通過幾何特征并不能有效地從云檢測結(jié)果中濾除這些非云區(qū)域。而改進MFC方法的云檢測結(jié)果中存在著更少的非云區(qū)域,說明進一步通過灰度共生矩陣提取紋理特征能夠?qū)@些非云區(qū)域進行更加有效的濾除。

本文針對冰雪區(qū)域的檢測和濾除所引入的光譜特征也表現(xiàn)出良好的效果。圖8列出了3組不同地點的遙感影像,每組從左至右依次為尺寸1 024像素×1 024像素的遙感影像及對應(yīng)的改進MFC方法的云檢測結(jié)果,其中白色表示云區(qū)域,黑色表示非云區(qū)域。圖中紅色矩形標(biāo)記了部分冰雪區(qū)域。圖8(a)是一幅含有大面積冰雪區(qū)域的遙感影像,可以看出,改進MFC方法濾除了大面積冰雪區(qū)域,防止其與云區(qū)域混淆。圖8(b)和圖8(c)是兩幅同時含有云和冰雪的遙感影像,可以看出,改進MFC方法有效地區(qū)分了云和冰雪,并將冰雪區(qū)域進行了濾除。

圖8 含有冰雪的遙感影像云檢測結(jié)果

表1和表2分別針對云量少于50%和云量高于50%的遙感影像,列出了原始MFC方法和改進MFC方法的總體精度、用戶精度和生產(chǎn)精度。表3列出了不同云量條件下的改進MFC方法性能參數(shù)的增長百分比。

表1 云量>50%的遙感影像云檢測結(jié)果的性能參數(shù) %

表2 云量<50%的遙感影像云檢測結(jié)果的性能參數(shù) %

表3 改進MFC方法云檢測結(jié)果性能參數(shù)的增長百分比 %

根據(jù)表1和表2可得出,相比于原始方法,改進方法在總體精度、用戶精度和生產(chǎn)精度上都有所提升;且相比于總體精度和生產(chǎn)精度,用戶精度的增長百分比較大,原因在于改進方法能夠濾除更多的建筑物、道路等非云區(qū)域,使得云檢測結(jié)果中包含更少的非云像素。

根據(jù)表3可得出,對于云量大于50%的遙感影像,改進MFC方法的總體精度、用戶精度和生產(chǎn)精度的增長百分比分別為0.46%、0.59%和0.27%,說明在云量較多的情況下改進方法的提升效果不顯著;對于云量小于50%的遙感影像,總體精度、用戶精度和生產(chǎn)精度的增長百分比分別為1.21%、4.37%和1.29%,說明在云量較少的情況下改進方法的提升效果較為顯著。

分析獲得上述結(jié)論的原因在于,通常相比于云區(qū)域,遙感影像中建筑物、道路等非云區(qū)域的面積較小。對于云量較多的遙感影像,這些非云區(qū)域所占的像素數(shù)量比例較小,其誤檢對云檢測結(jié)果的影響較小,即對精度的影響較小,從而改進方法的提升效果不明顯;反之,對于云量較少的遙感影像,這些非云區(qū)域所占的像素數(shù)量比例較大,其誤檢對云檢測結(jié)果的影響較大,即對精度的影響較大,從而改進方法的提升效果較為明顯。

3 結(jié)束語

本文驗證了MFC方法在Landsat 8遙感影像數(shù)據(jù)上的性能,并針對該方法進行了兩點改進。一是采用灰度共生矩陣提取紋理特征,并結(jié)合原始方法采用的幾何特征對初始云檢測結(jié)果中檢測為云的非云區(qū)域進行檢測和濾除。二是引入一種光譜特征,對初始云檢測結(jié)果中可能存在的冰雪區(qū)域進行檢測和濾除。通過實驗數(shù)據(jù)的收集和分析,驗證了MFC方法能夠有效檢測Landsat 8遙感影像中的云,并且驗證了改進方法的有效性,在一定程度上彌補了原始方法的缺點。但是改進方法仍存在局限性。一是雖然使用灰度共生矩陣提取紋理特征濾除了更多的非云區(qū)域,但是同時也會濾除一些極薄的薄云,原因在于一些薄云區(qū)域在藍(lán)色波段的反射率較低,與大部分非云區(qū)域相似。二是雖然引入的光譜特征能夠有效區(qū)分云區(qū)域和冰雪區(qū)域,但是其無法使用在一些沒有短波紅外波段的衛(wèi)星遙感影像上。三是對于云量較多或者含有較少易與云區(qū)域混淆的非云區(qū)域的遙感影像,改進方法對云檢測結(jié)果的提升效果有限。在后續(xù)的研究中,可以結(jié)合時間序列或其他方法,以進一步提高云檢測精度。