徐永強(qiáng)，周　俊，趙海富，仉樹軍，吳金亮

(1. 總參陸航研究所，北京 101121； 2. 中國電子科技集團(tuán)公司第五十四研究所，河北 石家莊 050081)

A Semi-inverse Approach to Detect Cloud

XU Yongqiang，ZHOU Jun，ZHAO Haifu，ZHANG Shujun，WU Jinliang

?

半逆法自動(dòng)云檢測(cè)

徐永強(qiáng)1，周俊1，趙海富1，仉樹軍1，吳金亮2

(1. 總參陸航研究所，北京 101121； 2. 中國電子科技集團(tuán)公司第五十四研究所，河北 石家莊 050081)

A Semi-inverse Approach to Detect Cloud

XU Yongqiang，ZHOU Jun，ZHAO Haifu，ZHANG Shujun，WU Jinliang

摘要：針對(duì)高分影像質(zhì)量判定的需求，提出色度過濾法和半逆法相結(jié)合的云檢測(cè)算法。根據(jù)色度過濾，識(shí)別可疑云區(qū)，生成搜索掩碼；利用半逆法，生成原圖的半逆圖像，度量色度差異，識(shí)別高亮云區(qū)；由高亮云區(qū)在掩碼區(qū)迭代搜索，進(jìn)一步達(dá)到對(duì)薄云的檢測(cè)，迭代策略的使用避免了城區(qū)對(duì)薄云區(qū)檢測(cè)的干擾。該方法在大量高分影像數(shù)據(jù)上進(jìn)行了試驗(yàn)驗(yàn)證， 均得到較好結(jié)果。

關(guān)鍵詞：色度過濾；半逆法；圖像分割；顏色空間

一、引言

高分一號(hào)自投入使用以來，在土地利用變更調(diào)查、地質(zhì)災(zāi)害監(jiān)測(cè)、生態(tài)環(huán)境監(jiān)測(cè)、耕地?cái)?shù)量調(diào)查、農(nóng)業(yè)設(shè)施分布調(diào)查等方面發(fā)揮著越來越重要的作用。首批高分影像數(shù)據(jù)包括2 m全色、8 m多光譜、2 m全色8 m多光譜融合、16 m多光譜4類。

衛(wèi)星在獲取影像時(shí)，很容易受到拍攝地區(qū)天氣狀況的影響，其中云區(qū)對(duì)地物的遮擋直接影像拍攝影像的可用性。云區(qū)遮擋地物信息不僅造成采集信息的缺失，而且給影像后續(xù)的使用帶來極大干擾，直接影響地物的分類、目標(biāo)的檢測(cè)、地物的識(shí)別等。影像數(shù)據(jù)接收后，對(duì)其進(jìn)行質(zhì)量檢測(cè)，識(shí)別其被云影響的區(qū)域范圍，有很大實(shí)用價(jià)值。高分影像數(shù)據(jù)量巨大，如果質(zhì)量檢測(cè)完全依賴人工操作，將消耗大量人力和時(shí)間。云自動(dòng)檢測(cè)算法自動(dòng)對(duì)影像中的云區(qū)進(jìn)行檢測(cè)，篩選排除被云高度污染的影像，將極大地減少人工操作，提高效率。

二、研究背景介紹

云主要由空氣中的水汽凝結(jié)而成，云中粒子直徑遠(yuǎn)大于太陽入射波長(zhǎng)，對(duì)可見光無選擇性反射及散射，即紅、綠和藍(lán)光的反射量相同，散射量也相同。因此在可見光波段中，云色調(diào)較亮呈現(xiàn)白色或灰白色，反射率接近1。

目前有大量云檢測(cè)算法被提出。Cao等[1]利用積分法對(duì)高分辨率影像進(jìn)行自動(dòng)云檢測(cè)；Smith等[2]利用天頂角和云透射率結(jié)合閾值分割云區(qū)；Souza-echer等[3]在HIS顏色空間中利用平行六面體法識(shí)別云區(qū)；楊俊等[4]利用插值法獲取局部閾值來識(shí)別云區(qū)； 曹瓊等[5]分析圖像的紋理特征進(jìn)行云區(qū)的提取。目前的云檢測(cè)算法主要是針對(duì)國外衛(wèi)星(Landsat、MODIS、NOAA等)影像進(jìn)行處理。國內(nèi)高分影像數(shù)據(jù)近期才發(fā)布，針對(duì)其影像數(shù)據(jù)的云檢測(cè)算法研究較少，并且空間分辨率、波譜參數(shù)不同，算法處理也不相同。本文算法主要針對(duì)國內(nèi)高分多光譜影像數(shù)據(jù)進(jìn)行自動(dòng)云檢測(cè)。GF-1衛(wèi)星有效光譜參數(shù)見表 1。

表1　GF-1衛(wèi)星有效光譜參數(shù)　μm

三、云區(qū)識(shí)別及提取

算法的整體處理流程如圖1所示。對(duì)于待檢測(cè)的高分影像數(shù)據(jù)，首先利用meanshift算法[7]對(duì)其進(jìn)行分塊預(yù)處理，以每個(gè)分割塊作為基本處理單元， 利用色度過濾法去除彩色區(qū)域，利用半逆法識(shí)別可信云區(qū)，再擴(kuò)延識(shí)別出薄云區(qū)域。

圖1　半逆法自動(dòng)云檢測(cè)流程

1. 分塊預(yù)處理

圖像分割作為云檢測(cè)預(yù)處理，合并相鄰的相似像元，分離相異像元。相比于影像中單個(gè)像素，分割出的區(qū)域具有相同或相似的光譜及紋理特性。本文采用meanshift算法對(duì)圖像進(jìn)行分塊處理。meanshift算法通過不斷迭代偏移均值，將相似像素進(jìn)行聚類合并。分割算法中最小分割區(qū)域的設(shè)置可以去除掉孤立的高亮像素點(diǎn)，提高算法的穩(wěn)定性及魯棒性。圖2(a)為輸入的高分影像；圖2(b)為其對(duì)應(yīng)的分割影像，每個(gè)分割塊的顏色值是區(qū)塊內(nèi)像素點(diǎn)的顏色值平均值。

圖2　CIE Lab色度過濾

2. CIE Lab色度過濾

將分塊的高分影像(圖2(b))由RGB顏色空間映射到CIE Lab 顏色空間。CIE Lab顏色空間為均勻顏色空間，視覺感知距離與顏色空間中點(diǎn)的距離成比例。其中，L代表亮度;a正值代表紅色，負(fù)值代表綠色；b正值代表黃色，負(fù)值代表藍(lán)色。CIE Lab顏色空間中的3個(gè)坐標(biāo)是相互獨(dú)立的。色度度量a,b絕對(duì)值的大小，度量了顏色偏移白色或灰白色的距離。從RGB顏色空間到CIE Lab顏色空間的變換公式如下

(1)

(2)

(3)

圖2(c)顯示了映射的強(qiáng)度分量，圖2(d)、(e)對(duì)a、b取絕對(duì)值，并映射到[0, 255]進(jìn)行顯示。圖2(d)、(e)顯示每個(gè)分割塊的色度跟強(qiáng)度大小無關(guān)，真實(shí)地反映了色彩強(qiáng)度。 圖2(f)為色度過濾生成的可疑云區(qū)掩碼圖，在算法處理中，掩碼圖生成的判斷標(biāo)準(zhǔn)為Dct<10。生成的掩碼圖覆蓋了圖像中表現(xiàn)接近灰色或白色的區(qū)域，過濾了山區(qū)植被區(qū)， 縮小后續(xù)云區(qū)識(shí)別的范圍。

3. 半逆法識(shí)別

筆者受He等工作[8]的啟發(fā)，并分析了大量含云高分多光譜影像數(shù)據(jù)，發(fā)現(xiàn)在云區(qū)中，每個(gè)像素點(diǎn)的每一個(gè)顏色通道都有較高值，而非云區(qū)中，像素分布滿足暗點(diǎn)原理, 即每個(gè)像素點(diǎn)的周圍都會(huì)包含暗點(diǎn)像素，即其至少在某一個(gè)顏色通道上有很小值(r、g或b)，接近于0, 具體原理細(xì)節(jié)可參考文獻(xiàn)[8]。Ancuti等[9]進(jìn)一步延伸暗點(diǎn)原理，提出了像素半逆法，用于霧效去除。 經(jīng)過大量高分多光譜影像的測(cè)試發(fā)現(xiàn)，半逆法在檢測(cè)云區(qū)方面也有較好效果，云區(qū)內(nèi)每個(gè)像素點(diǎn)各通道值均較大，進(jìn)行半逆處理，仍得到原始值，而對(duì)于在某一個(gè)通道上有較低值的像素點(diǎn)，進(jìn)行半逆處理后，色度有很大變化。

(4)

圖3　半逆法提取云區(qū)

4. 薄云區(qū)檢測(cè)

半逆法檢測(cè)了高亮顯示的云區(qū)，色度過濾排除了色彩明顯的地物的干擾。從高亮云區(qū)向周圍擴(kuò)延，進(jìn)一步對(duì)薄云區(qū)進(jìn)行檢測(cè)。如圖4(a)所示, 從半逆法檢測(cè)的高亮云區(qū)開始，基于色度過濾產(chǎn)生的掩碼(圖2(f)), 向周邊迭代遍歷檢測(cè)薄云區(qū)。 遍歷范圍為圖2(f)對(duì)應(yīng)的掩碼區(qū)域，在每個(gè)分割塊周圍臨域?qū)ふ冶≡茀^(qū)域，并迭代進(jìn)行。采用閾值法對(duì)薄云區(qū)進(jìn)行判別，RGB顏色通道的最小值大于90。 該判斷基于暗點(diǎn)原理，即非云區(qū)在某一個(gè)通道上值較小，接近于0。在擴(kuò)延過程中考慮近紅外波段的影響(近紅外波段為0.77～0.89 μm)，植被葉片結(jié)構(gòu)產(chǎn)生的多重反射，導(dǎo)致其反射率較高，通常用于識(shí)別植被覆蓋。擴(kuò)延中如果近紅外的值大于RGB通道的最小值，即停止外擴(kuò)。檢測(cè)結(jié)果如圖4(b)所示。

圖4　薄云區(qū)檢測(cè)

四、結(jié)果及比較

檢測(cè)結(jié)果圖4(b)所示， 不僅檢測(cè)了影像中的高亮云區(qū)，對(duì)其周圍的薄云區(qū)也有較好檢測(cè)結(jié)果。半逆法及薄云區(qū)檢測(cè)策略的使用使算法排除了城區(qū)對(duì)薄云區(qū)檢測(cè)的干擾，實(shí)現(xiàn)了魯棒的云檢測(cè)效果。圖5列出了其他一些檢測(cè)結(jié)果。圖5(a)為原始高分影像， 圖5(b)為檢測(cè)結(jié)果圖。為更好地顯示檢測(cè)結(jié)果，在檢測(cè)結(jié)果圖中添加了一個(gè)半透明的白色圖層，用于標(biāo)識(shí)檢測(cè)到的云區(qū)。在圖像分割中，使用同一組參數(shù)對(duì)原始影像進(jìn)行分塊預(yù)處理。在色度過濾中，閾值統(tǒng)一設(shè)置為10，該閾值的選擇考慮了云區(qū)顏色可能出現(xiàn)的色偏現(xiàn)象，提高了算法的適用性。半逆法中，閾值也統(tǒng)一設(shè)置為10， 對(duì)云區(qū)的顏色偏移，也有一定的適應(yīng)性。而薄云區(qū)閾值的選擇，是作了大量測(cè)試得出的經(jīng)驗(yàn)值。但大量試驗(yàn)證明，檢測(cè)效果很好。對(duì)于大小為4500像素×4548像素的多光譜影像，算法處理時(shí)間在10 s以內(nèi)。該時(shí)間主要消耗在對(duì)影像的分割預(yù)處理上，有待改近。算法的處理時(shí)間相對(duì)于人工檢測(cè)云覆蓋，已經(jīng)大大減少了。人工檢測(cè)涉及打開影像、查看、量算云區(qū)面積，整個(gè)工作流程繁瑣，耗時(shí)，還易出現(xiàn)漏查的情況。相對(duì)人工檢測(cè)，自動(dòng)云檢測(cè)提供了效率更高的選擇。

圖5　檢測(cè)結(jié)果

五、結(jié)束語

色度過濾法及半逆法的結(jié)合，實(shí)現(xiàn)了對(duì)高分影像數(shù)據(jù)的自動(dòng)云區(qū)檢測(cè)。檢測(cè)算法可以較好地排除城區(qū)對(duì)薄云區(qū)檢測(cè)的干擾。算法目前還達(dá)不到實(shí)時(shí)檢測(cè)的程度，瓶頸在于分割預(yù)處理，筆者正在尋找更快速的分割算法對(duì)其進(jìn)行替代，以實(shí)現(xiàn)對(duì)高分影像的實(shí)時(shí)自動(dòng)云檢測(cè)處理。

參考文獻(xiàn)：

[1]CAOQ，ZHENGH.AMethodforDetectingCloudinSatelliteRemoteSensingImageBasedonTexture[J].ActaAeronauticaetAstronauticaSinica, 2007, 28(3):661-665.

[2]SMITHS，TOUMIR．MeasuringCloudCoverandBrightnessTemperaturewithaGroundBasedThermalInfraredCamera[J]．JournalofAppliedMeteorolClimital，2008，47(2):683-693．

[3]SOUZA-ECHERMP，PEREIRAEB，BINSLS，etal．ASimpleMethodfortheAssessmentoftheCloudCoverStateinHighLatitudeRegionsbyaGroundBasedDigitalCamera[J]．JournalofAtmosphericandOceanicTechnology，2006，23(3):437-447.

[4]楊俊，呂偉濤，馬穎，等． 基于局部閾值插值的地基云自動(dòng)檢測(cè)方法[J]． 氣象學(xué)報(bào)，2010，68(6): 1007-1017．

[5]曹瓊,鄭紅,李行善.一種基于紋理特征的衛(wèi)星遙感圖像云探測(cè)方法[J].航空學(xué)報(bào)，2007, 28(3): 661-666.

[6]楊寧. 高分辨率影像面向?qū)ο蠓诸愄卣鬟x擇方法研究[D]. 西安: 西安科技大學(xué), 2012.

[7]COMMANICIUD，MEERP.MeanShift:ARobustApproachtowardFeatureSpaceAnalysis[J].IEEETransactionsonPatternAnalysisandMachineIntelligence, 2002, 24(5)：630-619.

[8]HEKM,SUNJ,TANGXO.SingleImageHazeRemovalUsingDarkChannelPrior[C]∥IEEEConferenceonComputerVisionandPatternRecognition.Miami,FL:IEEE, 2009: 1956-1963.

[9]ANCUTICO,ANCUTIC,HERMANSC,etal.AFastSemi-InverseApproachtoDetectandRemovetheHazefromaSingleImage[C]∥10thAsianConferenceonComputerVision.BerlinHeidelberg:Springer,2011: 501-514.

引文格式： 徐永強(qiáng)，周俊，趙海富，等. 半逆法自動(dòng)云檢測(cè)[J].測(cè)繪通報(bào)，2015(1)：63-65.DOI:10.13474/j.cnki.11-2246.2015.0012

作者簡(jiǎn)介：徐永強(qiáng)(1974—)，男，高級(jí)工程師，研究方向?yàn)榭仗煨畔?yīng)用。E-mail： xyq77@126.com

收稿日期：2014-07-10

中圖分類號(hào)：P237

文獻(xiàn)標(biāo)識(shí)碼：B

文章編號(hào)：0494-0911(2015)01-0063-03