羅冬　羅紅霞　劉光鵬　雷茜　馮華梅

摘要：為了提高多云霧地區(qū)光學(xué)遙感影像利用率、探索不同云量的光學(xué)影像和SAR影像的最佳融合尺度，選取重慶市渝中地區(qū)為研究區(qū)，基于Sentinel-1A極化合成孔徑雷達(dá)（SAR）影像與不同云量的Sentinel-2A（云量分別為0、10%、20%、30%）光學(xué)影像，進(jìn)行小波變換融合、乘積變換融合、高通濾波融合，再利用圖像評(píng)價(jià)方法評(píng)價(jià)影像融合以后的效果，最后利用面向?qū)ο蟮姆椒▽?duì)融合前后的影像進(jìn)行分類，利用混淆矩陣比較最終的分類精度。結(jié)果顯示，無(wú)云情況下，小波變換融合效果最好，最大限度地保留了原多光譜影像的亮度、反差，有效防止了影像信息的丟失，對(duì)植被的解譯能力有了明顯改善，而融合后的影像保真度較差，其余兩種融合效果相對(duì)次之，因此，在多云霧地區(qū)進(jìn)行多源遙感數(shù)據(jù)融合時(shí)，尤其是異質(zhì)數(shù)據(jù)融合時(shí)，優(yōu)先推薦小波融合算法;當(dāng)云量為10%以上時(shí)，3種融合算法雖然信息熵略有增加，但平均梯度、標(biāo)準(zhǔn)差減少了，造成解譯困難，導(dǎo)致最終分類精度均略低于Sentinel-2A影像，遠(yuǎn)低于Sentinel-1A影像，基本不能滿足使用要求，因此，在進(jìn)行地表覆被解譯時(shí)，推薦用SAR影像替代光學(xué)影像。

關(guān)鍵詞：影像融合;合成孔徑雷達(dá)（SAR）;地表覆被分類;面向?qū)ο蠓诸惙椒?多云霧地區(qū)

中圖分類號(hào)：P237 ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2020）05-0028-09

Abstract： This research aimed to improve the utilization ratio of optical remote sensing image in cloudy areas and explore the best fusion scales of different cloud cover optical remote sensing image and SAR image. Taking the Yuzhong area of Chongqing city as the research area， wavelet transform fusion， multiplicative transform fusion and high-pass filter fusion was carried out based on sentinel-1A polarized synthetic aperture radar （SAR） image and different cloud cover of sentinel-2A （cloud volume of 0， 10%， 20%， 30%） multi-spectral images， then the effect of image fusion was evaluated by image evaluation method， finally all images were classified by object-oriented methods，and the final classification accuracy was compared by the confusion matrix. The results show that， in the cloudless case， the wavelet transform has the best fusion effect， and the brightness and contrast of the original multi-spectral image are preserved to the maximum extent， which effectively prevents the loss of image information， and the interpretation ability of the vegetation is obviously improved， but after the fusion the image fidelity is poor， and the other two fusion effects are relatively inferior. Therefore， when multi-source remote sensing data fusion is performed in a cloudy fog region， especially when heterogeneous data is fused， the wavelet fusion algorithm is preferred. When the cloud cover is more than 10%， although the information entropy is slightly increased in the four fusion algorithms， the average gradient and standard deviation are reduced， which makes the interpretation difficult， resulting in the final classification accuracy is slightly lower than the sentinel-2A image， much lower than sentinel-1A images basically cannot meet the requirements for use. Therefore， when performing surface overlay interpretation， it is recommended to replace the optical images with SAR images.

Key words： image fusion; polarimetric synthetic aperture radar（sar）; land cover classification; object-oriented classification; cloudy areas

多云霧地區(qū)一直是遙感研究的瓶頸之一，這種地區(qū)云霧較多，其形成過(guò)程復(fù)雜，涉及物理、化學(xué)、湍流、輻射、動(dòng)力學(xué)以及下墊面狀況等諸多因素[1]，典型區(qū)域有四川盆地及周邊山區(qū)、東南沿海、西南季風(fēng)迎風(fēng)坡的云貴高原部分地區(qū)等[2]。由于云霧較多，這些區(qū)域可利用的光學(xué)遙感影像數(shù)量較少。由重慶市2017—2018年（行：39-40，列：127-129）的8景Landsat-8 OLI存檔數(shù)據(jù)（USGS）以及12景Sentinel-2 MSI1C多光譜存檔數(shù)據(jù)發(fā)現(xiàn)，Landsat-8影像云量低于5%的影像數(shù)量不到全年總量的4%，云量低于10%的影像數(shù)量只有全年總量的5%左右;而Sentinel-2影像云量低于5%的影像數(shù)量低于全年總量的5%，云量低于10%的影像數(shù)量不到全年總量的7%;結(jié)果顯示該地區(qū)可供利用的優(yōu)質(zhì)光學(xué)遙感影像數(shù)量較少。利用光學(xué)影像進(jìn)行資源、環(huán)境、災(zāi)害調(diào)查與監(jiān)測(cè)一直都是遙感研究的重要課題[3]，而由于氣候和地表環(huán)境等因素的影響，使得多云霧地區(qū)光學(xué)影像的利用效率受到限制。

合成孔徑雷達(dá)（Synthetic aperture radar，SAR）影像由于其高分辨率、多模式、多極化方向、數(shù)據(jù)獲取過(guò)程與自然光源無(wú)關(guān)，具有穿透云霧、全天候和全天時(shí)對(duì)地觀測(cè)、對(duì)地物形狀敏感的特點(diǎn)，彌補(bǔ)了光學(xué)遙感的不足[4]。融合光學(xué)遙感影像和SAR影像已成為提高光學(xué)遙感影像利用效率的重要手段之一，且大量研究證明融合光學(xué)影像和SAR影像能提高光學(xué)影像的利用率，尤其是在地表解譯時(shí)的分類精度[5，6]。

目前，光學(xué)遙感影像和SAR影像的融合算法主要還是借鑒了光學(xué)遙感影像這類同質(zhì)影像的融合算法，主要分為像元級(jí)、特征級(jí)、決策級(jí)3個(gè)級(jí)別[7-10]，大部分的融合算法還是基于像元級(jí)的融合，特征級(jí)和決策級(jí)的融合算法相對(duì)較少，像元級(jí)典型的融合算法有IHS（Intensity hue saturation，IHS）變換、乘積變換、Brovey變換、高通濾波、主成分分析、小波變換、GS（Gram-schmidt Pan，GS）以及支持向量機(jī)等。IHS融合方法簡(jiǎn)單，容易實(shí)現(xiàn)。賈永紅等[11]最早進(jìn)行了4種IHS變換方法用于多光譜影像與SAR影像融合的試驗(yàn)比較。Brovey變換對(duì)融合后的影像光譜信息保持較好，但受波段限制;乘積變換對(duì)大的地貌類型效果好，同時(shí)可用于多光譜與高光譜的融合。賈永紅等[12，13]進(jìn)行了分辨率變換融合（Brovey變換、乘積變換等）、小波分析、證據(jù)推理的遙感融合分析，發(fā)現(xiàn)融合算法能明顯提高分類精度;向海燕[14]采用乘積變換、Brovery變換將ALOS-2全極化影像和Landsat OLI影像進(jìn)行融合并對(duì)融合后的影像進(jìn)行解譯發(fā)現(xiàn)，融合后的影像分類總精度和Kappa系數(shù)均有明顯提高。融合后的影像既有光學(xué)影像的光譜信息，又有全極化SAR影像的極化分解特征，具有很強(qiáng)的信息互補(bǔ)性。小波變換融合對(duì)圖像的視覺特征、明暗色調(diào)、紋理特征反映較好，既提高了光學(xué)影像的空間分辨率，又保留了其光譜特征[15，16]。王智均等[17]研究了不同小波變換的遙感影像數(shù)據(jù)融合算法，證明小波變換在融合時(shí)比別的算法更具優(yōu)勢(shì);以上這些研究都取得了較好的效果。Garzelli[18]、Wang等[19]以及王宇航[20]的研究表明，與單一采用多光譜影像進(jìn)行地表覆被分類相比，SAR和多光譜影像融合以后的影像更有優(yōu)勢(shì)，兩種數(shù)據(jù)具有互補(bǔ)性和合作性;在地形測(cè)繪、土地利用分類、農(nóng)作物分類、森林和濕地監(jiān)測(cè)、城市監(jiān)測(cè)、洪澇災(zāi)害監(jiān)測(cè)等實(shí)際運(yùn)用方面具有重要意義[21-25]。而對(duì)于多云霧地區(qū)，可利用的光學(xué)影像數(shù)量稀少，多源遙感數(shù)據(jù)融合是解決這一問題的重要途徑之一。因此，本研究選擇了2018年不同云量（0、10%、20%、30%）的4景Sentinel-2A數(shù)據(jù)和相對(duì)應(yīng)的4景Sentinel-1A，進(jìn)行多云霧地區(qū)不同云量光學(xué)影像和SAR影像融合的適應(yīng)性研究。原則上云量越少融合的效果越好，但由于多云霧地區(qū)光學(xué)影像數(shù)量稀少，有時(shí)不得不使用部分的有云霧數(shù)據(jù)，所以考慮到現(xiàn)實(shí)條件以及適用性，本研究選擇小波變換（多分辨率）、乘積變換（波段運(yùn)算）、高通濾波（濾波模型）3種典型的融合算法融合4景不同云量的Sentinel-2A多光譜影像和相對(duì)應(yīng)的4景Sentinel-1A SAR影像，以求提高多云霧地區(qū)光學(xué)遙感影像的利用率、彌補(bǔ)光學(xué)遙感影像的缺點(diǎn)，探索不同云量光學(xué)遙感影像和SAR影像最佳融合方法，發(fā)揮各自的長(zhǎng)處，實(shí)現(xiàn)多源遙感影像的綜合運(yùn)用。

1 ?研究區(qū)與數(shù)據(jù)源

1.1 ?研究區(qū)

研究區(qū)位于重慶市的部分主城區(qū)，主要包括渝中、渝北、江北、沙坪壩、南岸、北碚、九龍坡，介于106°23′9.58″—106°36′36.36″E，29°28′16.50″—29°39′48.95″N;東西被中梁山和銅鑼山包圍，中間為嘉陵江和長(zhǎng)江的交匯地帶;屬于典型的山地城市，地勢(shì)起伏大，最低海拔108 m，最高海拔663 m;地形以山地、丘陵為主;屬亞熱帶濕潤(rùn)性季風(fēng)氣候，降水豐富。

1.2 ?數(shù)據(jù)源

Sentinel-1A是歐空局ERS-1/2和Envisat的后續(xù)衛(wèi)星，Sentinel-1A搭載了C波段SAR傳感器，具有干涉寬（IW，Interferometric wide swath）、條帶成像（SM，Stripmap）、波浪（Wave mode）、超寬幅（EW，Extra wide swath）4種成像模式;可實(shí)現(xiàn)單極化、雙極化，它和Sentinel-1B組網(wǎng)后重訪周期可由12 d縮短成6 d。本研究選取了4景2018年的Sentinel-1A地距影像（GRD，Ground range detected）產(chǎn)品，具體參數(shù)見表1。

Sentinel-2A衛(wèi)星是全球環(huán)境與安全監(jiān)測(cè)計(jì)劃的第二顆衛(wèi)星，于2015年6月23日發(fā)射。Sentinel-2A攜帶一枚多光譜成像儀，可覆蓋13個(gè)光譜波段，最高空間分辨率10 m、單顆重訪周期10 d，2顆衛(wèi)星組網(wǎng)重訪周期可達(dá)5 d。從可見光和近紅外到短波紅外具有不同的空間分辨率，本研究選取了4景2018年的Sentinel-2A MSL1C產(chǎn)品，具體參數(shù)見表1。

2 ?研究方法

2.1 ?數(shù)據(jù)預(yù)處理

2.1.1 ?Sentinel-1A數(shù)據(jù)預(yù)處理 ?利用歐空局的ESA_SNAP6.0版本軟件對(duì)Sentinel-1A數(shù)據(jù)進(jìn)行軌道校正、熱噪聲去除、輻射定標(biāo)、多視處理、濾波、地形校正以及后向散射歸一化等預(yù)處理，濾波選用的是GAMMA-MAP濾波，濾波窗口大小設(shè)置為5×5，然后將圖像重采樣為10 m×10 m分辨率的影像，用RGB 3個(gè)通道顯示預(yù)處理后的VV極化、VH極化、VV極化，結(jié)果如圖1所示。

2.1.2 ?Sentinel-2A數(shù)據(jù)預(yù)處理 ?利用歐空局的Sen2cor插件對(duì)Sentinel-2A影像進(jìn)行輻射定標(biāo)、大氣校正等預(yù)處理，再利用ESA_SNAP、ENVI將Sentinel-1A和Sentinel-2A影像進(jìn)行配準(zhǔn)，重投影，裁剪出研究區(qū)范圍[本研究盡量裁剪出云量平均分布的區(qū)域，但受數(shù)據(jù)影響，只能裁剪出大概的符合梯度（0、10%、20%、30%）標(biāo)準(zhǔn)的研究區(qū)]，裁剪后的Sentinel-2A RGB合成圖如圖2所示。

2.2 ?圖像融合方法

2.2.1 ?小波變換 ?小波變換融合與傅里葉變換類似。首先，對(duì)光學(xué)影像和合成孔徑雷達(dá)影像進(jìn)行匹配，使這兩個(gè)影像具有相似的概率分布，相應(yīng)的小波系數(shù)具有相同的等級(jí);其次，對(duì)合成孔徑雷達(dá)影像和光學(xué)影像分別進(jìn)行小波分解，對(duì)兩幅圖像的小波系數(shù)進(jìn)行了比較，選取較大的部分組成新的小波分解系數(shù);最后通過(guò)新系數(shù)重建圖像得到融合圖像。小波變換融合算法常用于多光譜影像（如TM）與合成孔徑雷達(dá)影像的融合，它具有提高影像空間分辨率的同時(shí)又能保持影像飽和度、色調(diào)等不變的優(yōu)點(diǎn)[26]。

2.2.2 ?乘積變換 ?乘積變換融合是一種典型的代數(shù)運(yùn)算融合算法，相似的還有比值變換等，其原理是通過(guò)對(duì)兩幅影像的像元值進(jìn)行相乘得到融合圖像，缺點(diǎn)是融合圖像的亮度值高速增長(zhǎng)，其計(jì)算公式如下：

Rnew=I×RSentinel-2A+RSentinel-1A

Gnew=I×GSentinel-2A+GSentinel-1A

Bnew=I×BSentinel-2A+BSentinel-1A

式中，Rnew、Gnew、Bnew為融合圖像紅、綠、藍(lán)波段的像元值，I為經(jīng)驗(yàn)常數(shù)，RSentinel-2A、GSentinel-2A、BSentinel-2A為Sentinel-2A圖像紅、綠、藍(lán)波段的像元值，RSentinel-1A、GSentinel-1A、BSentinel-1A為Sentinel-1A合成圖像中的紅、綠、藍(lán)波段的像元值。

2.2.3 ?高通濾波 ?高通濾波融合算法在融合SAR影像和多光譜影像時(shí)可以簡(jiǎn)單地理解為兩個(gè)影像信息的加權(quán)求和，對(duì)高空間分辨率影像進(jìn)行高通濾波以提取線性特征和邊緣等空間結(jié)構(gòu)信息，對(duì)低分辨率影像進(jìn)行低通濾波以提取其光譜、后向散射強(qiáng)度等信息，然后將高通濾波和低通濾波的結(jié)果進(jìn)行加權(quán)求和，從而得到融合影像。融合后的影像中高分辨率影像的高頻信息顯而易見，低頻部分卻減少了[27]。

2.2.4 ?圖像融合質(zhì)量評(píng)價(jià) ?對(duì)融合后影像的質(zhì)量進(jìn)行定量評(píng)價(jià)更能直接反映融合效果，衡量融合后的影像在光譜信息、空間信息和信息量等方面的改善情況。本研究主要統(tǒng)計(jì)了融合前后影像的平均梯度、標(biāo)準(zhǔn)差、信息熵3個(gè)參數(shù)以對(duì)融合前后影像的質(zhì)量進(jìn)行定量評(píng)價(jià)[28]。

3 ?結(jié)果與分析

3.1 ?融合結(jié)果分析與評(píng)價(jià)

Sentinel-1A和Sentinel-2A影像融合結(jié)果如圖3所示。從色彩上看，3種融合算法中只有乘積變換、高通濾波融合較好地保留了原始多光譜影像的色彩;小波變換則相反，其對(duì)不同地表的覆被賦予了不同的色彩，且不同云量也被賦予了不同的色彩，側(cè)面印證了其保真度較差。從云霧改善狀況來(lái)看，小波變換大面積改善了云霧狀況，但部分厚云區(qū)域仍未完全去除。乘積變換和高通濾波對(duì)云霧區(qū)域的改善僅限部分面積較小的薄云區(qū)域，厚云區(qū)域基本未得到改善，但后者對(duì)薄云區(qū)域的融合效果明顯優(yōu)于前者。因此，綜合來(lái)看，小波變換融合后的影像效果較好，其實(shí)際效果需要進(jìn)一步通過(guò)質(zhì)量定量評(píng)價(jià)參數(shù)進(jìn)行評(píng)價(jià)。

融合后影像的質(zhì)量定量評(píng)價(jià)結(jié)果如表2所示。從表2可以看出，在無(wú)云條件下，3種融合算法的平均梯度都明顯增加，說(shuō)明融合后的影像更清晰，層次感更豐富;信息熵除高通濾波外，其余2種融合算法均明顯增加，說(shuō)明這3種融合算法能更多地保留原影像的紋理等細(xì)節(jié)信息，有利于目標(biāo)分割;標(biāo)準(zhǔn)差僅小波變換融合有明顯增加，乘積變換和高通濾波大幅下降，說(shuō)明只有小波變換較好地保留了原影像的亮度和光譜信息。從以上3個(gè)指標(biāo)來(lái)看，云量為0時(shí)，3種融合算法中小波變換融合最優(yōu)。當(dāng)云量增加到10%和20%以后，3種融合算法的平均梯度、標(biāo)準(zhǔn)差均明顯下降，說(shuō)明3種融合算法部分保留了原影像的光譜、亮度等信息;信息熵只有小波變換有明顯增加，表明這2種算法繼承了原影像的紋理等細(xì)節(jié)信息?？偟膩?lái)說(shuō)，云量為10%～20%時(shí)，小波變換融合效果較好（圖4）。當(dāng)云量為30%時(shí)，小波變換融合的平均梯度、標(biāo)準(zhǔn)差均下降，信息熵有所增加，但乘積變換和高通濾波卻相反，平均梯度有明顯增加，而標(biāo)準(zhǔn)差和信息熵均明顯下降。由此可知，云量為30%時(shí)，3種融合算法的融合效果均較差。

總之，無(wú)論是目視效果還是定量評(píng)價(jià)，結(jié)果都顯示小波變換對(duì)融合以后影像的質(zhì)量改善最大，其保留了原多光譜影像的亮度與反差，融合圖像能更好地反映圖像的細(xì)節(jié)特征，對(duì)植被、河流的解譯能力有了很大的提高;乘積變換和高通濾波在云量大幅增加后對(duì)原圖像的亮度和反差保留較少，紋理細(xì)節(jié)丟失比較嚴(yán)重，地物間可分辨性降低，地物解譯難度變大。

3.2 ?分類結(jié)果分析

3.2.1 ?面向?qū)ο蠓诸惙椒▍?shù)設(shè)置 ?為了評(píng)價(jià)融合后影像的實(shí)際運(yùn)用效果，采用面向?qū)ο笞钹徑诸惙▽?duì)融合前后的影像進(jìn)行解譯[31]，解譯軟件采用eCognition9.0。因?yàn)檠芯繀^(qū)土地類型比較單一，加上研究區(qū)具體的區(qū)域狀況，經(jīng)過(guò)多次試驗(yàn)，影像分類參數(shù)設(shè)置如表3所示。

3.2.2 ?分類結(jié)果精度評(píng)價(jià) ?對(duì)影像進(jìn)行分類評(píng)價(jià)的目的在于再次驗(yàn)證融合后影像的效果，驗(yàn)證樣本主要來(lái)源于研究區(qū)野外實(shí)地調(diào)查點(diǎn)、土地利用現(xiàn)狀數(shù)據(jù)和Google Earth衛(wèi)星影像等共230個(gè)驗(yàn)證樣本點(diǎn)（圖4），采用混淆矩陣進(jìn)行分類精度評(píng)價(jià)，同時(shí)還分析了云量為0時(shí)的影像融合前后的具體地類分類精度及其分類結(jié)果。

1）分類精度總體評(píng)價(jià)。在質(zhì)量評(píng)價(jià)結(jié)果的基礎(chǔ)上，進(jìn)一步利用影像解譯后的分類精度對(duì)融合后的影像進(jìn)行總體分析評(píng)價(jià)，結(jié)果見表4。無(wú)云時(shí)，3種融合算法解譯結(jié)果與單一采用多光譜影像或SAR影像相比，總精度提高了近1個(gè)百分點(diǎn)或以上。當(dāng)云量增加到10%和20%時(shí)，3種融合算法融合后影像的總精度與光學(xué)影像相比，下降了最多5個(gè)百分點(diǎn)，與SAR影像相比最多下降了10個(gè)百分點(diǎn)，結(jié)合圖像質(zhì)量評(píng)價(jià)結(jié)果來(lái)看，其原因是融合后的影像雖然信息量有增加，但圖像清晰度和反差變小了，地物可分離性降低，圖像層次模糊，不利于地表覆被解譯。

當(dāng)云量為30%時(shí)，小波變換融合后影像的總精度與光學(xué)影像相比下降了不到1個(gè)百分點(diǎn)，與SAR影像相比總精度下降了14個(gè)百分點(diǎn)，Kappa系數(shù)減少了0.18，而乘積變換、高通濾波兩種算法的精度則大幅度減少，與光學(xué)影像相比下降了約2個(gè)百分點(diǎn)，與SAR影像相比下降了16個(gè)百分點(diǎn)。結(jié)合前面表2分析其原因在于，小波變換的平均梯度、標(biāo)準(zhǔn)差明顯減少了，信息熵略有增加，但仍不利于解譯，而乘積變換和高通濾波的平均梯度雖然有明顯增加，但標(biāo)準(zhǔn)差和信息熵均減少，圖像清晰度變差，影像信息量減少，更不利于地表覆被解譯。

小波變換在3種云量情況下的總精度和Kappa系數(shù)都表現(xiàn)較好，而圖像質(zhì)量評(píng)價(jià)結(jié)果也印證了這一點(diǎn)，從而表明小波變換方法對(duì)融合后影像的視覺特征、明暗色調(diào)、紋理特征有較好改善。而乘積變換和高通濾波較之差一些，尤其是當(dāng)云量增加后，融合后的影像結(jié)構(gòu)比較破碎，色調(diào)層次感差，增加了地表覆被解譯的難度。

2）云量為0時(shí)的影像分類精度分析。云量為0時(shí)，從總精度和Kappa系數(shù)來(lái)看，融合后影像對(duì)解譯的精度有提升但不是特別明顯（圖5），因此進(jìn)一步分析了云量為0%時(shí)4種融合算法對(duì)具體地類的分類效果（圖6）。

從圖5可看出，云量為0時(shí)，融合后的影像對(duì)林地和灌草地的解譯精度都有明顯提高。對(duì)人工建筑和水域的解譯，融合后的影像沒有明顯優(yōu)勢(shì)，多光譜影像對(duì)人工建筑解譯精度最高，SAR影像對(duì)水域解譯的效果最好。

綜上所述，無(wú)云時(shí)，3種融合算法的解譯精度較好，能基本滿足現(xiàn)實(shí)需要;當(dāng)云量增加到10%以上時(shí)，光學(xué)數(shù)據(jù)的分類精度急劇下降，3種融合算法不但沒有明顯改善分類效果且還有小幅度下降;當(dāng)云量增加到20%以上時(shí)，融合效果和分類精度大幅下降。因此，當(dāng)光學(xué)影像云量超過(guò)10%進(jìn)行地表覆被解譯時(shí)，建議用SAR影像替代光學(xué)影像，基本可以獲得80%左右的分類精度，或者尋找別的替代數(shù)據(jù)。

通過(guò)對(duì)融合后影像的質(zhì)量評(píng)價(jià)和分類結(jié)果及云量為0時(shí)影像分類后的制圖精度分析可知，小波變換在最大限度地保留了原多光譜影像的亮度和光譜信息的同時(shí)，還提高了總體分類精度，對(duì)植被解譯精度較大，融合效果明顯優(yōu)于其余2種融合算法。

影像質(zhì)量評(píng)價(jià)結(jié)果顯示，高通濾波、乘積變換在融合后影像信息量增加了，但是卻不能顯著提高融合后的圖像分類精度，這可能與光學(xué)影像上厚云和薄云的占比和分布情況有關(guān)，同時(shí)還可能受分類方法、分類參數(shù)選擇等因素影響。

SAR影像雖然基本不受云層影響，但是從融合影像分類結(jié)果來(lái)看，微波影像對(duì)薄云區(qū)域的融合效果相較于厚云區(qū)域來(lái)說(shuō)有一定提升，這可能是因?yàn)楸≡茀^(qū)域地表覆被信息頻率較多，而厚云區(qū)域地表覆被信息頻率幾乎為零的原因。

無(wú)云時(shí)，融合后影像的分類總精度與多光譜影像和SAR影像相比雖有一定提高，但沒有顯著提高，原因可能是受分類方法和驗(yàn)證樣本點(diǎn)的影響。

4 ?小結(jié)與討論

多源遙感影像數(shù)據(jù)在地形測(cè)繪、土地利用分類、地圖更新、農(nóng)作物監(jiān)測(cè)分類、森林監(jiān)測(cè)分類、洪澇災(zāi)害監(jiān)測(cè)等領(lǐng)域發(fā)揮了重要作用。本研究選擇了2018年不同云量（0、10%、20%、30%）的4景Sentinel-2A影像，首先對(duì)雙極化C波段Sentinel-1A SAR影像進(jìn)行軌道校正、熱噪聲去除、輻射校正、多視處理及濾波、地形校正和輻射歸一化等預(yù)處理，再與輻射定標(biāo)和大氣校正后不同云量的Sentinel-2A影像進(jìn)行配準(zhǔn)并進(jìn)行小波融合、乘積融合、高通濾波融合，然后評(píng)價(jià)融合前后影像的質(zhì)量。最后基于面向?qū)ο蟮淖钆R近分類法對(duì)融合前后的影像進(jìn)行分類，并利用混淆矩陣進(jìn)行精度分析與評(píng)價(jià)，比較融合前后影像的質(zhì)量和分類精度，得出結(jié)論如下。

1）無(wú)云時(shí)，Sentinel-1A SAR影像的數(shù)據(jù)質(zhì)量、分類精度和Kappa系數(shù)一般都低于Sentinel-2A多光譜影像。3種融合算法中，小波變換融合效果有明顯提升，最大限度地保留了原多光譜影像的亮度和光譜信息， 防止了影像信息的丟失，更好地反映了圖像的細(xì)節(jié)特征， 對(duì)植被的解譯效果有明顯改善， 缺點(diǎn)是圖像的保真度較差，但3種融合算法相比，建議在進(jìn)行光學(xué)影像和SAR影像融合時(shí)，優(yōu)先采用小波變換融合。

2）與單一遙感數(shù)據(jù)相比，多源遙感數(shù)據(jù)融合后所提供的信息具有互補(bǔ)性和合作性，但在多云霧地區(qū)，云量低于10%時(shí)，多源遙感數(shù)據(jù)融合能在一定程度上提高數(shù)據(jù)的質(zhì)量和利用價(jià)值。當(dāng)云量大于10%后，數(shù)據(jù)融合的質(zhì)量和運(yùn)用價(jià)值逐漸減小，因此，在多云霧地區(qū)進(jìn)行地表覆被解譯時(shí)，當(dāng)光學(xué)影像云量大于10%以上時(shí)，不建議進(jìn)行光學(xué)影像和SAR影像融合，可直接用SAR影像代替光學(xué)影像或者尋找別的替代數(shù)據(jù)進(jìn)行地表覆被解譯。

3）融合后的影像對(duì)植被的解譯精度有較大提升，對(duì)人工建筑和水域的解譯則沒有明顯優(yōu)勢(shì)。

本研究存在的不足：①本研究區(qū)屬典型的山地地形，部分地區(qū)地形起伏較大，坡度達(dá)到了60°，雖進(jìn)行了地形校正等預(yù)處理，但并未完全消除SAR影像上頂?shù)椎怪煤屠走_(dá)陰影等現(xiàn)象，對(duì)SAR影像和多光譜影像的融合效果以及分類結(jié)果都可能造成干擾，而2種影像的分辨率也略有差異，雖重采樣為相同的分辨率，但均可能會(huì)對(duì)分類結(jié)果和分類精度產(chǎn)生一定影響。②受數(shù)據(jù)限制，本研究不同梯度的Sentinel-2A影像云量沒有完全在圖像上均勻分布，薄云和厚云均有，這在一定程度上會(huì)影響影像融合效果和具體地物的分類精度。本研究對(duì)云量為10%以上時(shí)的建議針對(duì)的是厚云，且云霧的具體分布也得視具體情況而定，如在具體運(yùn)用時(shí)，影像以薄云為主，那么這個(gè)尺度就不能作為參考，還需要進(jìn)一步研究。

參考文獻(xiàn)：

[1] 林 ?琿，邵 ?蕓，楊立民.多云多雨地區(qū)遙感初探[J].遙感信息，2004（4）：71-74.

[2] VIHMA T，SCREEN J，TJERNSTR？魻M M，et al. The atmospheric role in the Arctic water cycle：A review on processes，past and future changes，and their impacts[J].Journal of geophysical research：Biogeosciences，2016，121（3）：586-620.

[3] 賈永紅，李德仁，孫家柄.多源遙感影像數(shù)據(jù)融合[J].遙感技術(shù)與應(yīng)用，2000，15（1）：41-44.

[4] 張良培，沈煥鋒.遙感數(shù)據(jù)融合的進(jìn)展與前瞻[J].遙感學(xué)報(bào)，2016， 20（5）：1050-1061.

[5] GAETANO R，COZZOLINO D，AMIANO L D，et al. Fusion of sar-optical data for land cover monitoring[A].2017 IEEE International geoscience and remote sensing symposium （IGARSS）[C].Fort Worth，TX，USA：IEEE，2017.5470-5473.

[6] ZHU L，TATEISHI R. Fusion of multisensor multitemporal satellite data for land cover mapping[J].International journal of remote sensing，2007，27（5）：903-918.

[7] GHASSEMIAN H. A review of remote sensing image fusion methods[J].Information fusion，2016，32：75-89.

[8] LIU Z，BLASCH E，JOHN V. Statistical comparison of image fusion algorithms：Recommendations[J].Information fusion，2017， 36：251-260.

[9] POHL C，VAN GENDEREN J L. Review article multisensor image fusion in remote sensing：Concepts，methods and applications[J].International journal of remote sensing，2010，19（5）：823-854.

[10] ZHANG J. Multi-source remote sensing data fusion：Status and trends[J].International journal of image and data fusion，2010， 1（1）：5-24.

[11] 賈永紅，李德仁.SAR與TM影像的IHS變換復(fù)合及其質(zhì)量定量評(píng)價(jià)[J].國(guó)土資源遙感，1997（3）：36-41.

[12] 賈永紅，李德仁.多源遙感影像像素級(jí)融合分類與決策級(jí)分類融合法的研究[J].武漢大學(xué)學(xué)報(bào)（信息科學(xué)版），2001（5）：430-434.

[13] 賈永紅，李德仁，孫家柄.一種基于小波變換的圖像融合算法[J].電子學(xué)報(bào)，2004（5）：750-753.

[14] 向海燕.基于光學(xué)和SAR數(shù)據(jù)的多云霧山區(qū)土地覆被分類[D].重慶：西南大學(xué)，2018.

[15] TANG Z W，WANG J G，HUANG S J. Wavelet transform application for image fusion[J/OL].Wavelet applications VII，2000，4056：https：//doi.org/10.1117/12.381706.

[16] SANLI F B，ABDIKAN S，ESETLILI M T，et al. Evaluation of image fusion methods using PALSAR，RADARSAT-1 and SPOT images for land use/land cover classification[J].Journal of the Indian society of remote sensing，2017，45（4）：591-601.