資源三號(hào)測(cè)繪衛(wèi)星自動(dòng)云檢測(cè)

陳振煒，張 過(guò)，寧津生，唐新明

1.武漢大學(xué)測(cè)繪學(xué)院，湖北 武漢430079；2.武漢大學(xué)測(cè)繪遙感信息工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，湖北 武漢430079；3.國(guó)家測(cè)繪地理信息局衛(wèi)星測(cè)繪應(yīng)用中心，北京100830

1 引 言

由于云層覆蓋了地表上空50%以上的面積，資源三號(hào)測(cè)繪衛(wèi)星的數(shù)據(jù)在進(jìn)行傳感器校正產(chǎn)品［1］的生產(chǎn)之前，需要剔除云量覆蓋過(guò)大的影像，以減輕生產(chǎn)負(fù)擔(dān)，提高生產(chǎn)效率；另一方面，估計(jì)出影像的云量值，有助于后續(xù)的應(yīng)用。

目前遙感影像的云量檢測(cè)方法主要有兩大類［2］：一類是基于遙感影像的光譜特性［3－6］，這種方法主要是利用紅外波段，通過(guò)分析云和地物在同一個(gè)波長(zhǎng)上的反射率或亮溫值不同，以及云和地物的反射率或亮溫值隨著波長(zhǎng)的變化而變化情況不同，從而進(jìn)行云地分離。這種方法的云檢測(cè)效果雖然較好，但是要求影像有足夠的波段范圍，適用性有限。另一類方法是基于圖像處理的手段［7－8］，通過(guò)分析影像上云和地物的紋理特征，提取合適的特征或特征組合，區(qū)分云和地物。這種方法的問(wèn)題主要在于，遙感影像上云的種類相當(dāng)繁雜，形態(tài)和紋理上沒(méi)有統(tǒng)一的特征，在各個(gè)特征空間內(nèi)的分布都不集中，而地物的種類也是如此，因此使用紋理特征進(jìn)行云地區(qū)分有較大的難度［10］。文獻(xiàn)［11］采用小波SCM 提取紋理特征，并利用基于ASM和熵的雙判別方式辨別云層，取得了一定的成果，但受到傳感器類型不同、影像分辨率不同等因素的限制，用于資源三號(hào)衛(wèi)星影像時(shí)仍存在較高的誤判率。

本文研究對(duì)象為資源三號(hào)測(cè)繪衛(wèi)星原始數(shù)據(jù)編目處理時(shí)所生成的瀏覽圖，尺寸均為1024像素×1024像素。針對(duì)云地種類繁雜的問(wèn)題，首先對(duì)待檢測(cè)的影像進(jìn)行了增強(qiáng)處理，然后在多尺度空間下分析圖像的紋理特征，實(shí)現(xiàn)云地的分離。最后選取了993幅瀏覽圖作為待測(cè)樣本進(jìn)行試驗(yàn)。

2 資源三號(hào)測(cè)繪衛(wèi)星自動(dòng)云檢測(cè)

2.1 樹(shù)判別結(jié)構(gòu)云檢測(cè)流程

對(duì)資源三號(hào)測(cè)繪衛(wèi)星編目生成的瀏覽圖采用樹(shù)狀判別結(jié)構(gòu)［12］的方式進(jìn)行高效率的云量檢測(cè)，如圖1所示。

算法采用了灰度均值、方差、一階差分、二階矩等圖像特征，經(jīng)過(guò)大量樣本的訓(xùn)練獲取特征最佳閾值。首先將云圖進(jìn)行分割成64×64的子圖塊，對(duì)于每個(gè)子圖，根據(jù)灰度均值，分為3種情況：灰度均值高于235的子圖，圖像近乎全白，無(wú)明顯的紋理和灰度變化，云的概率極高，此時(shí)將該子圖直接歸為云類；對(duì)于灰度均值低于80的子圖，圖像整體表現(xiàn)較暗，地物概率較高，但也可能包含了所占子圖面積比例較小的云塊，這種情況相當(dāng)于暗地表上空存在著亮度明顯高于地表的云，此時(shí)圖像會(huì)呈現(xiàn)較高的方差值，因此對(duì)于灰度均值低于80的子圖，若其方差小于所設(shè)閾值，則認(rèn)為是下墊面，反之則說(shuō)明有云的存在，可根據(jù)實(shí)際工程進(jìn)行歸類，對(duì)于資源三號(hào)的實(shí)際工程應(yīng)用，一般將這種情況下的子圖歸為云，或者可根據(jù)灰度值的高低測(cè)定子圖中云的比例，參與最后的統(tǒng)計(jì)；對(duì)于灰度均值介于80至235之間的子圖，首先對(duì)其進(jìn)行直方圖均衡化，以強(qiáng)調(diào)其隱含的紋理細(xì)節(jié)，經(jīng)過(guò)該過(guò)程，地物的紋理細(xì)節(jié)得到強(qiáng)化，而云類由于本身不含較細(xì)的紋理，因此變化較小，此時(shí)可以結(jié)合圖像的灰度值和二階矩進(jìn)行初步分類，二階矩的臨界取值f0為該子圖未經(jīng)直方圖均化時(shí)的灰度均值gm的函數(shù)，通過(guò)大量樣本可以粗略確定

圖1 樹(shù)狀判別結(jié)構(gòu)云檢測(cè)流程Fig.1 Cloud detection workflow with tree discriminate structure

即意味著在較高灰度均值的條件下，云的概率高于地物，則二階矩的閾值降低。當(dāng)子圖二階矩高于相應(yīng)的閾值時(shí)，歸為云類，否則還有待確定。經(jīng)過(guò)大量樣本的試驗(yàn)可以發(fā)現(xiàn)，待確定的子圖比例約為52%，其中主要為地物，還包括一定比例的云。對(duì)于這部分在特征上云地混雜的子圖，進(jìn)行一階圖像尺度分解，并求取分解前后圖像的一階差分退化率，將退化率與子圖未經(jīng)直方圖均化的原始灰度均值相結(jié)合，通過(guò)大量樣本訓(xùn)練出退化率臨界值d0與灰度均值gm的關(guān)系，最終區(qū)分云地，此時(shí)退化率的臨界值仍可以表示為灰度均值的函數(shù)

該流程中，灰度均值雖然不能直接作為判別云地的特征，但是貫穿了整個(gè)算法流程，起到了重要的參考作用。

2.2 子圖切分

在實(shí)際工程運(yùn)用中，云量值作為評(píng)價(jià)數(shù)據(jù)質(zhì)量的一個(gè)指標(biāo)，并不需要很高的精度，因此本文先對(duì)待檢測(cè)的影像進(jìn)行子圖分割，對(duì)分割后的小尺寸子圖，進(jìn)行云地分類，即將云目標(biāo)識(shí)別的問(wèn)題轉(zhuǎn)換為子圖的云地分類問(wèn)題，降低了云檢測(cè)的復(fù)雜程度。然后統(tǒng)計(jì)全圖中被分為云類的子圖所占的比例，即為原圖的云量值。同時(shí)根據(jù)云類子圖的分布情況，也能得到影像上云層的大致覆蓋區(qū)域。

不同的子圖尺寸會(huì)對(duì)云檢測(cè)的準(zhǔn)確性帶來(lái)一定的影響，若子圖的尺寸過(guò)大，則包含的信息較多，影像中可能同時(shí)含有云類和地物類，或者含有不同類型的云，由于圖像的特征是全圖的綜合反應(yīng)，此時(shí)的子圖不能作為一個(gè)云類或地物類的單元，且子圖尺寸越大時(shí)，全圖的取樣數(shù)就越少，會(huì)導(dǎo)致準(zhǔn)確率受到影響；若子圖的尺寸過(guò)小，則所含信息量小，不足以體現(xiàn)云和地物的紋理區(qū)別，甚至出現(xiàn)部分云和地物的子圖表現(xiàn)幾乎一致，人眼也難以區(qū)分的情況。針對(duì)資源三號(hào)衛(wèi)星影像的編目瀏覽圖，本文通過(guò)試驗(yàn)后選擇尺寸為64×64的大小進(jìn)行子圖分割。

2.3 特征增強(qiáng)和提取

經(jīng)典的基于紋理分析的云檢測(cè)算法，主要思路是［8，12－15］：首先對(duì)影像進(jìn)行子圖切分，然后提取圖像的紋理特征，并采用合適的分類法，對(duì)子圖進(jìn)行云和地的區(qū)分。傳統(tǒng)紋理分析法云檢測(cè)依賴于提取的特征，在這方面前人已作了大量的工作，如使用分?jǐn)?shù)維描述紋理復(fù)雜度和粗糙度［13］、使用灰度梯度共生矩陣描述紋理統(tǒng)計(jì)特征［14］等。另外，熵、方差、一階差分、分?jǐn)?shù)維等反應(yīng)圖像紋理粗細(xì)或細(xì)節(jié)復(fù)雜度的特征［8，15］也是經(jīng)常用于云檢測(cè)的圖像特征。因?yàn)橐话愣栽圃趫D像上的灰度分布均勻、跳變程度小，紋理較粗且模糊。反之，地物的紋理復(fù)雜度和灰度跳變程度較大，但是對(duì)于本文研究的對(duì)象資源三號(hào)衛(wèi)星編目生成的瀏覽圖，其采樣率較低，在這個(gè)分辨率下，建筑物、道路、植被等地物不可見(jiàn)，僅山脈、河流、沙漠、田野等大面積的地物類型能夠呈現(xiàn)在影像上。而這其中有不少的地物類型包括了以下特征：梯度較小、紋理較粗、灰度值較高。這與云類的圖像特征相似，容易產(chǎn)生混淆。如沙漠地區(qū)的影像往往就有這個(gè)現(xiàn)象（圖2）。另外，在該分辨率下的云類，尤其是層云以及體積較小的塊狀云，會(huì)呈現(xiàn)出較大的梯度和較細(xì)的紋理，甚至是較低的灰度值（圖3），加大了直接通過(guò)紋理特征來(lái)區(qū)分云地的難度。

圖2 典型的沙地影像Fig.2 Typical sand images

圖3 紋理較細(xì)的云影像Fig.3 Cloud with thin textures

基于上述考慮，首先對(duì)分割后得到的每個(gè)子圖進(jìn)行了直方圖均衡化的增強(qiáng)處理，以有目的性地改變圖像的紋理特征，均衡化通過(guò)控制圖像灰度級(jí)概率密度函數(shù)，改變圖像的灰度層次，可表示如下［16］

式中，sk為原圖灰度值為k的像素經(jīng)過(guò)變換后的新灰度值；T（rk）表示變換函數(shù)；pr（ri）是灰度值為i的像素頻率；L即圖像的灰度等級(jí)。

均衡化的目的主要在于突出隱含有紋理細(xì)節(jié)的圖像。由于地物包含著較豐富的紋理信息，其細(xì)節(jié)清晰度比均衡化之前會(huì)有較大提高，如圖4和圖5所示，沙漠和城市地區(qū)的原始影像灰度變化程度小，其中沙地的紋理顯得十分平坦，特征類似于云，但經(jīng)過(guò)直方圖均化之后，呈現(xiàn)出明顯的顆粒感，灰度跳變幅度大，與云類有顯著的差異。城市地區(qū)經(jīng)過(guò)直方圖均化則突出了建筑物等細(xì)節(jié)。而云類本身不隱含較復(fù)雜的細(xì)節(jié)，因此均化前后的紋理特性變化較小，如圖6和圖7所示，云類在經(jīng)過(guò)直方圖均化后，只有對(duì)比度增強(qiáng)，紋理卻沒(méi)有變細(xì)，灰度的空間分布仍保持較好連續(xù)性。

圖4 沙地影像直方圖均衡化前后對(duì)比Fig.4 Comparison results before－and－after histogram equalization for sand images

圖5 城市影像直方圖均衡化前后對(duì)比Fig.5 Comparison results before－and－after histogram equalization for city images

圖6 積云影像直方圖均衡化前后對(duì)比Fig.6 Comparison results before－and－after histogram equalization for cumulus cloud images

圖7 薄云影像直方圖均衡化前后對(duì)比Fig.7 Comparison results before－and－after histogram equalization for thin cloud images

2.4 多尺度特征分析

對(duì)同一事物進(jìn)行不同尺度的觀測(cè)，可以得到不同的結(jié)果。當(dāng)遙感影像尺度較大時(shí)，主要看到的是地塊和海域，在尺度縮小、分辨率隨之升高的過(guò)程中，山脈、湖泊等地物則顯現(xiàn)出來(lái)，而當(dāng)尺度繼續(xù)縮小之后，則可以看到街道、樓房、樹(shù)木等細(xì)節(jié)。相應(yīng)的，在不同尺度不同分辨率下，圖像的特征在會(huì)產(chǎn)生變化。經(jīng)過(guò)大量的試驗(yàn)觀測(cè)可發(fā)現(xiàn)：對(duì)于某些類型的目標(biāo)，其圖像特征會(huì)隨圖像的尺度變化而發(fā)生明顯的變化；而對(duì)于另外一些類型的目標(biāo)，則對(duì)圖像的尺度變化具有較強(qiáng)的魯棒性，即其特征受圖像尺度變化的影響很小。因此，可以考慮通過(guò)改變圖像分辨率，研究不同目標(biāo)的各類特征隨圖像尺度變化而變化的情況，來(lái)實(shí)現(xiàn)分類的目的。由于單一尺度下的特征空間中的云類和地物存在一定的混雜區(qū)域，為了減少誤判和漏判，本文將特征的多尺度延拓引入光學(xué)遙感影像云量檢測(cè)的研究，在多尺度下分析圖像特征的變化率，通過(guò)提取云類和地物在尺度變化情況下圖像特征變化的差異，從而進(jìn)行進(jìn)一步的云檢測(cè)算法優(yōu)化。

首先，通過(guò)低通抽樣濾波對(duì)圖像進(jìn)行多尺度分解，盡管每一級(jí)影像都是對(duì)上一級(jí)的向下采樣，但是由于采樣之前進(jìn)行了濾波，因此采樣間隔的點(diǎn)對(duì)于采樣后的圖像也有貢獻(xiàn)，在一定程度上保留了圖像的信息。本文對(duì)圖像先進(jìn)行5×5的卷積濾波再進(jìn)行隔點(diǎn)采樣［19］。設(shè)第k階的濾波圖像為gk，尺寸為Rk×Ck，則各階圖像有遞推公式

式中，需滿足歸一性、對(duì)稱性和均勻分布性，在這些條件約束下，此時(shí)圖像的多尺度分解類似于高斯金字塔分解。

由于云固有的物理特性，在圖像尺度不斷降低的過(guò)程中，云的紋理特性變化不大，而地物由于含有較豐富的細(xì)節(jié)，隨著圖像尺度降低，其紋理細(xì)節(jié)必然會(huì)不斷損失，且灰度的跳變程度也由于低通濾波的緣故不斷變小。這從直觀的角度解釋了云和地物在圖像分解過(guò)程中的變化情況是不同的。

從圖像特征的層面分析，對(duì)于每一級(jí)經(jīng)過(guò)分解的影像，都可以提取出前文所述的各種特征，并且這些特征會(huì)隨著尺度的變化而變化，這個(gè)過(guò)程稱為特征退化，并稱某個(gè)特征在第k＋n級(jí)分解圖與第k級(jí)分解圖中的取值之比為該特征的k級(jí)n階退化率

本文利用特征退化率，在多尺度下進(jìn)行進(jìn)一步的圖像分析，提高了云檢測(cè)的準(zhǔn)確率。

3 試驗(yàn)分析

3.1 特征增強(qiáng)比較試驗(yàn)

以下針對(duì)前文所示的圖4—圖7中4張影像，分析它們?cè)谥狈綀D均衡化前后的圖像特征變化情況，如表1所示。

表1 幾張影像經(jīng)過(guò)直方圖均衡化前后的特征變化情況Tab.1 Changes of image feature parameters before－andafter histogram equalization

圖8—圖11對(duì)表1中的數(shù)據(jù)進(jìn)行直觀展示，形象地體現(xiàn)了圖像特征經(jīng)過(guò)直方圖均衡化之后的變化情況。

結(jié)合表1和圖8—圖11可以看出，經(jīng)過(guò)直方圖均衡化處理之后，一些圖像特征發(fā)生了明顯的變化。而在均衡化處理之前，云和地物的這些特征都沒(méi)有分離，但經(jīng)過(guò)直方圖均衡化之后，在一些圖像特征上，圖6、圖7的云和圖4、圖5的地物有了極其明顯的分離，例如，均衡化之前沙地的一階差分值最小，為0.755，但經(jīng)過(guò)均衡化之后，增加到41.93；同樣城市的一階差分值經(jīng)過(guò)均衡化也從3.178急劇上升到39.00；相反云類的一階差分值，雖然也有所增加，但分別只增加到6.34和10.87，尤其是圖6的厚云，均衡化前的一階差分是4組數(shù)據(jù)中最大的，但均衡化后增量很小，結(jié)果遠(yuǎn)小于沙地和城市兩類地物；再比如二階矩，該特征越小表示紋理越細(xì)，經(jīng)過(guò)均衡化之后，兩類地物的二階矩均遠(yuǎn)小于兩類云，還原了地物紋理細(xì)的特點(diǎn)；同樣，經(jīng)過(guò)均衡化之后，兩類地物的慣性矩增量明顯，達(dá)到8以上，而云類的慣性矩變化則相對(duì)較小，均衡化后的慣性矩值均小于1，與地物產(chǎn)生明顯差異，由于慣性矩體現(xiàn)了紋理的溝紋深淺，由此可見(jiàn)，經(jīng)過(guò)均化的地物很好地還原了紋理較深的特點(diǎn)。

圖8 云地的方差經(jīng)均衡化前后的變化Fig.8 Variance changes of clouds and grounds before－and－after histogram equalization

圖9 云地的一階差分經(jīng)均衡化前后的變化Fig.9 First difference changes of clouds and grounds before－and－after histogram equalization

圖10 云地的二階矩經(jīng)均衡化前后的變化Fig.10 Secondary moment changes of clouds and grounds before－and－after histogram equalization

圖11 云地的慣性矩經(jīng)均衡化前后的變化Fig.11 Inertia moment changes of clouds and grounds before－and－after histogram equalization

通過(guò)試驗(yàn)分析可知，直方圖均衡化后的子圖有以下優(yōu)點(diǎn)：

（1）增強(qiáng)了部分紋理特征，拉大了云地差距，使兩者的混雜區(qū)域明顯變小。圖12和圖13比較了195個(gè)云類樣本和181個(gè)地物類樣本經(jīng)過(guò)均衡化前后，云和地物在一階差分特征空間中分布的變化情況。可以看出，經(jīng)過(guò)直方圖均衡化之后，云類和地物類在一階差分這個(gè)特征上的分布區(qū)域被拉開(kāi)，地物的一階差分值整體高于云類。

圖12 均衡化前云類和地物類子圖的一階差分值分布情況Fig.12 Distribution of first difference of clouds and grounds before histogram equalization

圖13 均衡化后云類和地物類子圖的一階差分值分布情況Fig.13 Distribution of first difference of clouds and grounds after histogram equalization

（2）提高了云檢測(cè)算法的適用性。傳統(tǒng)紋理分析法受以下幾點(diǎn)因素限制：①受傳感器限制，不同傳感器的影像可能存在分辨率不同、輻射情況不同等區(qū)別，使某些云檢測(cè)方法只適用于個(gè)別傳感器，對(duì)其他傳感器的衛(wèi)星影像則效果較差；②受產(chǎn)品生產(chǎn)和瀏覽圖生成方式的限制，0級(jí)數(shù)據(jù)或?yàn)g覽圖的重采樣方式、色階壓縮變換方式等處理方法不同，使用于云判的影像的圖像特征有區(qū)別，這也是導(dǎo)致云判算法適用性有限的原因之一。但經(jīng)過(guò)色調(diào)均化后，在一定程度上統(tǒng)一了圖像，使云檢測(cè)算法的適用性得到提高。

（3）有利于薄云的檢測(cè)。薄云由于具有一定的透明性，一直是云檢測(cè)的難點(diǎn)。首先其透明性導(dǎo)致圖像特征值往往介于云類和地物類之間，容易誤判；其次，對(duì)于透明度較高的薄云，可以看見(jiàn)其下墊面，往往歸為地物，而透明度較低的薄云，下墊面模糊不清，則一般歸為云類，這為薄云的檢測(cè)增加了復(fù)雜性。通過(guò)直方圖均衡化的圖像，薄云的存在被強(qiáng)化，且對(duì)于透明度高的薄云，其下墊面的紋理會(huì)得到一定程度的恢復(fù)，而對(duì)于透明度低的薄云，則呈現(xiàn)出云的形態(tài)。如圖14所示，經(jīng)過(guò)直方圖均衡化后，薄云透明度高的區(qū)域，地物的紋理得到了較好的還原，而薄云透明度低的區(qū)域，則呈現(xiàn)出云的形態(tài)，這就便于薄云的檢測(cè)。

圖14 薄云影像均衡化前（a）后（b）對(duì)比Fig.14 Comparison of thin cloud images before（a）－and－after（b）histogram equalization

對(duì)子圖進(jìn)行直方圖均化處理是資源三號(hào)測(cè)繪衛(wèi)星編目瀏覽圖云檢測(cè)中最關(guān)鍵的步驟之一。通過(guò)這一處理，可以使錯(cuò)綜復(fù)雜的云類和地物類的區(qū)分工作得到極大的簡(jiǎn)化，不僅提高了計(jì)算效率，也提高了分類結(jié)果的準(zhǔn)確率。

本文選取了195個(gè)云類子圖樣本和181個(gè)地物類子圖樣本，經(jīng)過(guò)直方圖均衡化預(yù)處理之后，其中一階差分值分布如前文圖13所示?？梢钥闯?，云類樣本在經(jīng)過(guò)直方圖均衡化之后的一階差分值主要分布在0～5之間，而地物樣本則分布在5～30之間，說(shuō)明經(jīng)過(guò)直方圖均衡化后，圖像的一階差分值能夠較好地體現(xiàn)云地特性，在一定程度上起到云地區(qū)分的作用，但在一階差分值為5的上下區(qū)域，仍有不少云類和地物混合在一起，因此僅靠一階差分并不能準(zhǔn)確地區(qū)分云地。

另外，圖15和圖16分別反映了經(jīng)過(guò)直方圖均衡化處理后樣本中，云地的二階矩和慣性矩分布情況。由圖15可以看出，二階矩具有良好的云地區(qū)分能力，經(jīng)過(guò)直方圖均衡化之后的大多數(shù)云的二階矩值均小于地物。但仍有約10%的云類和地物混雜在一起。而圖16中反映出慣性矩也具有一定的區(qū)分云地的能力，但效果遠(yuǎn)不如二階矩，因?yàn)榈匚锏膽T性矩值分布很不集中，廣泛地分布在0～4之間，而云類的慣性矩值普遍較低，但也有一定比例的云具有較高的慣性矩值。從圖像的整體角度上看，說(shuō)明地物的圖像既有較深的溝紋，也有較淺的溝紋，圖像從清晰到模糊均有；而云類的圖像溝紋普遍較淺，灰度變化較緩且小，效果模糊，但也存在一部分云圖具有較深的溝紋和較清晰的效果。

圖15 云和地物的二階矩值分布圖Fig.15 Distribution of secondary moment of clouds and grounds

圖16 云和地物的慣性矩值分布圖Fig.16 Distribution of inertia moment of clouds and grounds

3.2 多尺度特征提取試驗(yàn)

灰度均值對(duì)尺度的變化具有很強(qiáng)的魯棒性，圖像分解的過(guò)程中灰度均值變化很小；而一階差分在圖像分解中呈現(xiàn)出較大的變異性，且該特征的退化率與云地分類具有較高的契合度。本文選取了在單一尺度下誤判的云地樣本各50張，圖17為這些樣本的一階差分特征的0級(jí)1階退化率分布情況?？梢钥闯?，通過(guò)多尺度的分析，可以進(jìn)一步改善云檢測(cè)效果。

3.3 云檢測(cè)綜合試驗(yàn)

為了驗(yàn)證該云檢測(cè)算法，本文選取了多幅不同資源三號(hào)測(cè)繪衛(wèi)星編目瀏覽圖進(jìn)行試驗(yàn)。

圖17 云和地物的一階差分0級(jí)1階退化率分布情況Fig.17 First－order degeneration distribution of first difference of clouds and grounds of level 0images

圖18（a）為一幅云含量約為60%的遙感影像，影像上的云灰度大多在145～210之間，地物的灰度值大多在70～185之間，圖上包含了層云和薄云，云與地物沒(méi)有明顯的分界。根據(jù)本文的云檢測(cè)算法進(jìn)行處理，得到檢測(cè)結(jié)果如圖18（b）所示。其中黑色部分為灰度均值低于80的子圖，由于在這個(gè)灰度水平下，云的概率極低，因此可以直接認(rèn)為是地物。白色部分代表了檢測(cè)出來(lái)的云，其余部分為地物，計(jì)算得到云量比例為57.13%，符合預(yù)期結(jié)果。

圖18 原始影像（a）和檢測(cè)結(jié)果示意圖（b）Fig.18 Original image（a）and result image（b）

為了分析特征多尺度延拓對(duì)云量檢測(cè)效果的影響，試驗(yàn)選取了一幅包含了層云、卷云和積云的遙感影像，如圖19（a）所示。該圖的云含量為100%，在實(shí)際生產(chǎn)中，一般視為廢片。對(duì)該圖在單一尺度下進(jìn)行特征分析，檢測(cè)出屬于云類的子圖，如圖19（b）所示，灰色區(qū)域?yàn)樵谶@一步尚未檢測(cè)出云的子塊。這部分子塊內(nèi)的影像，溝紋較細(xì)，紋理復(fù)雜度較高，和地物存在一定的相似性，這體現(xiàn)了單一尺度下特征分析的局限性。對(duì)該圖進(jìn)行特征的多尺度分析，即求取一階差分在各個(gè)子圖經(jīng)過(guò)分解之后的0級(jí)1階退化率，灰色區(qū)域?yàn)橥ㄟ^(guò)征多尺度延拓分析后所檢測(cè)出來(lái)的云子圖，經(jīng)過(guò)圖像的多尺度分析，云檢測(cè)的效果得到了很大的提升。

圖19 原始影像（a）和多尺度檢測(cè)結(jié)果（b）Fig.19 Original image（a）and result image in multi－scale（b）

為了進(jìn)一步驗(yàn)證算法對(duì)薄云的檢測(cè)能力，試驗(yàn)選取了兩張含有大量薄云的影像，其中圖20（a）中的薄云絕大部分不具備透明性，其下墊面已模糊不清，只有陸地和水域的交界還依稀可見(jiàn)；而圖20（b）的薄云透明程度不均勻，其透明的部分可以看到地面信息，尤其是對(duì)于傳感器校正產(chǎn)品［1］，具有遠(yuǎn)高于相應(yīng)縮略圖的分辨率，能夠通過(guò)去云的方式還原下墊面的影像，基于這點(diǎn)考慮，當(dāng)云的透明度高到一定程度時(shí)，可歸為地物類，以減少不必要的數(shù)據(jù)浪費(fèi)。

圖20 均勻薄云（a）和不均勻薄云（b）Fig.20 Uniformly thin clouds and non－uniformly thin clouds

對(duì)圖20使用本文的檢測(cè)算法進(jìn)行云檢測(cè)，得到結(jié)果如圖21所示，同樣白色部分代表云子圖。圖21（a）所計(jì)算出的云量值為94.34%，而圖21（b）所計(jì)算得到的云量值為62.11%，均基本與預(yù)期相符。

另外，本文選取了993幅資源三號(hào)測(cè)繪衛(wèi)星編目生成的瀏覽圖作為試驗(yàn)待測(cè)樣本，通過(guò)人工目視給定云量值，然后根據(jù)所提出的云量檢測(cè)算法計(jì)算所有圖中所含云量。待測(cè)樣本的總數(shù)據(jù)量為160MB，總計(jì)算時(shí)間為74s，平均每副影像檢測(cè)時(shí)間為0.075s。對(duì)于實(shí)際工程應(yīng)用而言，通常評(píng)價(jià)云量自動(dòng)檢測(cè)效果的指標(biāo)并非檢測(cè)精度，而是檢測(cè)結(jié)果的準(zhǔn)確率。若設(shè)定與人工結(jié)果相差±5%以內(nèi)為準(zhǔn)確，否則為不準(zhǔn)確，則符合準(zhǔn)確條件的影像個(gè)數(shù)為852幅，準(zhǔn)確率為85.80%；若設(shè)定與人工結(jié)果相差±10%以內(nèi)為準(zhǔn)確，則符合準(zhǔn)確條件的影像個(gè)數(shù)為931幅，準(zhǔn)確率為93.76%。

圖21 均勻（a）和不均勻（b）薄云檢測(cè)結(jié)果Fig.21 Result image of uniformly thin clouds（a），result image of non－uniformly thin clouds（b）

4 結(jié) 論

本文對(duì)資源三號(hào)測(cè)繪衛(wèi)星編目瀏覽圖的云量檢測(cè)進(jìn)行了分析，對(duì)待檢測(cè)的影像進(jìn)行子圖分割，將子圖作為云量檢測(cè)的基本單元，把云目標(biāo)識(shí)別問(wèn)題轉(zhuǎn)化為子圖的云地分類問(wèn)題。在使用圖像特征進(jìn)行分類之前，采用直方圖均衡化對(duì)原影像進(jìn)行預(yù)處理，以突出地物的紋理細(xì)節(jié)，解決了傳統(tǒng)的紋理分析方法進(jìn)行云量檢測(cè)時(shí)所面臨的云地在各個(gè)特征空間中分布不集中且混雜程度高的問(wèn)題，分析并提取了灰度均值、方差、一階差分、二階矩等圖像特征作為云地分類的有效特征，并對(duì)圖像特征進(jìn)行多尺度延拓分析，將特征退化率與云地分類相結(jié)合，采用樹(shù)狀判別結(jié)構(gòu)的方式進(jìn)行云量檢測(cè)。試驗(yàn)結(jié)果表明，云檢測(cè)算法達(dá)到了較高的檢測(cè)準(zhǔn)確率。

［1］ PAN Hongbo，ZHANG Guo，TANG Xinming，et al.The Geometrical Model of Sensor Corrected Products for ZY－3Satellite［J］.Acta Geodaetica et Cartographica Sinica，2013，42（4）：516－522.（潘紅播，張過(guò)，唐新明，等.資源三號(hào)測(cè)繪衛(wèi)星傳感器校正產(chǎn)品幾何模型［J］.測(cè)繪學(xué)報(bào)，2013，42（4）：516－522.）

［2］ HARUMA I，TAKASHI Y N.Development of an Unbiased Cloud Detection Algorithm for a Spaceborne Multispectral Imager［J］.Journal of Geophysical Research－Atmospheres，2009.DOI：10.1029／2008JD010710.

［3］ YANG Tieli，HE Quanjun.Cloud Detection in MODIS Data［J］.Journal of Anshan University of Science and Technology，2006（2）：162－166.（楊鐵利，何全軍.MODIS數(shù)據(jù)的云檢測(cè)處理［J］.鞍山科技大學(xué)學(xué)報(bào)，2006（2）：162－166.）

［4］ PASQUALE L S，MICHELLE A B，JOHN L D.Devel－opment of the Landsat Data Continuity Mission Cloud－Cover Assessment Algorithms［J］.IEEE Transactions on Geoscience and Remote Sensing，2012，50（4）：1140－1154.

［5］ PAN Cong，XIA Bin，CHEN Yu，et al.Cloud Detection of MODIS Data Based on FCM［J］.Microcomputer Information，2009（4）：124－125，131.（潘聰，夏斌，陳彧，等.基于模糊聚類的MODIS云檢測(cè)算法研究［J］.微計(jì)算機(jī)信息，2009（4）：124－125，131.）

［6］ GUO Hongtao，WANG Yi，LIU Xiangpei，et al.Integrated Optimal Method of Cloud Detection with Meteorological Satellite Data［J］.Journal of PLA University of Science and Technology（Natural Science Edition），2010，11（2）：221－227.（郭洪濤，王毅，劉向培，等.衛(wèi)星云圖云檢測(cè)的一種綜合優(yōu)化方法［J］.解放軍理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，11（2）：221－227.）

［7］ YAN Yusong，LONG Teng.Real－Time Cloud Detection in Optical Remote Sensing Image［J］.Transactions of Beijing Institute of Technology，2010，30（7）：817－821.（閆宇松，龍騰.遙感圖像的實(shí)時(shí)云判技術(shù)［J］.北京理工大學(xué)學(xué)報(bào)，2010，30（7）：817－821.）

［8］ ZHOU Quan.The Classification of Remote Sensing Cloud Images Based on Multiscale Features［D］.Hefei：University of Science and Technology of China，2010.（周全.基于多尺度特征的遙感云圖分類［D］.合肥：中國(guó)科學(xué)技術(shù)大學(xué)，2010.）

［9］ ZHOU Lijuan.Cloud Detection Algorithm Research For Optical Satellite Image［D］.Xi’an：XiDian University，2012.（周麗娟.可見(jiàn)光衛(wèi)星圖像的云檢測(cè)算法研究［D］.西安：西安電子科技大學(xué)，2012.）

［10］ CHEN P Y，SRINIVASAN R，F(xiàn)EDOSEJEVS G，et al.An Automated Cloud Detection Method for Daily NOAA－14AVHRR Data for Texas，USA［J］.Remote Sensing，2002，23（15）：2939－2950.

［11］ TAO Shuping，JIN Guang，ZHANG Guixiang，et al.A Wavelet SCM Algorithm Used to Detect Cloud in Remote Sensing Cameras［J］.Acta Geodaetica et Cartographica Sinica，2011，40（5）：598－603.（陶淑蘋，金光，張貴祥，等.實(shí)現(xiàn)遙感相機(jī)自主辨云的小波SCM算法［J］.測(cè)繪學(xué)報(bào)，2011，40（5）：598－603.）

［12］ SHEN Na，ZHENG Tianyao，WANG Zhensong.High－speed and High－accuracy Algorithm for Cloud Detection and Its Application［J］.Journal of Remote Sensing，2009，13（6）：1138－1155.

［13］ XIA Deshen，JIN Sheng，WANG Jian.Fractal Dimension and GGCM Based Meteorology Cloud Pictures Recognition－Description of Texture Complexity and Roughness by Fractal Dimension（Ⅰ）［J］.Journal of Nanjing University of Science and Technology，1999，23（3）：90－93.（夏德深，金盛，王健.基于分?jǐn)?shù)維與灰度梯度共生矩陣的氣象云圖識(shí)別（Ⅰ）——分?jǐn)?shù)維對(duì)紋理復(fù)雜度和粗糙度的描述［J］.南京理工大學(xué)學(xué)報(bào)，1999，23（3）：90－93.）

［14］ XIA Deshen，JIN Sheng，WANG Jian.Fractal Dimension and GGCM Based Meteorology Cloud Pictures Recognition－Description of Texture Statistical Characteristics by Gray Level－gradient Co－occurrence Matrix（Ⅱ）［J］.Journal of Nanjing University of Science and Technology，1999，23（4）：289－292.（夏德深，金盛，王健.基于分?jǐn)?shù)維與灰度梯度共生矩陣的氣象云圖識(shí)別（Ⅱ）——灰度梯度共生矩陣對(duì)紋理統(tǒng)計(jì)特征的描述［J］.南京理工大學(xué)學(xué)報(bào)，1999，23（4）：289－292.）

［15］ CHEN Peng，ZHANG Rong，LIU Zhengkai.Feature Detection for Cloud Classification in Remote Sensing Images［J］.Journal of University of Science and Technology of China，2009（5）：484－488.（陳鵬，張榮，劉政凱.遙感圖像云圖識(shí)別中的特征提?。跩］.中國(guó)科學(xué)技術(shù)大學(xué)學(xué)報(bào)，2009（5）：484－488.）

［16］ GONZALEZ C R.Digital Image Processing［M］.2nd ed.Beijing：Electronic Industry Press，2008：72－74.（岡薩雷斯.數(shù)字圖像處理［M］.第2版.北京：電子工業(yè)出版社，2008：72－74.）

［17］ WALKER R F，JACKWAY P T，LONGSTAFF I D.Recent Developments in the Use of the Co－occurrence Matrix for Texture Recognition［C］∥Proceedings of 13th International Conference on Digital Image Processing.［S.l.］：IEEE，1997：63－65.

［18］ FENG Jianhui，YANG Yujing.Study of Texture Images Extraction Based on Gray Level Co－occurence Matrix［J］.Beijing Surveying and Mapping，2007（3）：19－22.（馮建輝，楊玉靜.基于灰度共生矩陣提取紋理特征圖像的研究［J］.北京測(cè)繪，2007（3）：19－22.）

［19］ ZHOU Quan，ZHANG Rong，YIN Dong.The Extracting of Multiscale Features of Remote Sensing Cloud Images Based on Gaussian Pyramid［J］.Remote Sensing Technology and Application，2010（5）：604－608.（周全，張榮，尹東.基于高斯金字塔的遙感云圖多尺度特征提?。跩］.遙感技術(shù)與應(yīng)用，2010（5）：604－608.）

［20］ JIA Yonghong.Digital Image Processing［M］.Wuhan：Wuhan University Press，2003：182－184.（賈永紅.數(shù)字圖像處理［M］.武漢：武漢大學(xué)出版社，2003：182－184.）