91制片厂果app,麻豆麻豆必出精品入口,718吃瓜网

Landsat 8 OLI多光譜與全色影像融合算法的比較*

姜曉晨,，鄧正棟，武國瑛,，王東豪

（陸軍工程大學(xué) 國防工程學(xué)院,，江蘇南京 210007）

摘要： 以Landsat8 OLI衛(wèi)星遙感影像為數(shù)據(jù)源，以大伙房水庫為研究區(qū),，采用Brovery融合法,、NNDiffuse Pan Sharpening融合法和Gram-Schmidt融合法對比研究了多光譜和全色數(shù)據(jù)的融合問題。經(jīng)過定性和定量分析發(fā)現(xiàn),，Gram-Schmidt融合法和NN融合法擁有更豐富的光譜信息,，且光譜保真度最優(yōu)的是Gram-Schmidt融合法。針對Gram-Schmidt融合法,，選取了水體邊界線處的影像像元進行光譜分析,，研究結(jié)果發(fā)現(xiàn),，融合后的水體邊界線呈現(xiàn)得更為清晰,，有助于細(xì)化水體提取邊界,。

關(guān)鍵詞： 融合算法 Landsat8 多光譜與全色光譜分析

中圖分類號：TP301
文獻標(biāo)識碼：A
DOI： 10.19358/j.issn.2096-5133.2018.08.008
中文引用格式：姜曉晨,，鄧正棟,，武國瑛,，等.Landsat 8 OLI多光譜與全色影像融合算法的比較［J］.信息技術(shù)與網(wǎng)絡(luò)安全,，2018,37（8）：31-35.

Comparisons on fusion algorithms of Landsat 8 OLI multi-spectral and panchromatic images

Jiang Xiaochen, Deng Zhengdong, Wu Guoying, Wang Donghao

(College of National Defense Engineering, The Army Engineering University of PLA, Nanjing 210007, China)

Abstract： Three image fusion algorithms of Brovery, NNDiffuse Pan Sharpening(NND) and Gram-Schmidt(GS) were compared when the multi-spectral and panchromatic images of Landsat 8 OLI were merged，taking the Dahuofang Reservoir as the study area. From the qualitative and quantitative analyses, it was concluded that GS and NND had richer spectral information and GS had the best spectrum fidelity. Based on the GS method, the image elements at the borderline of the water body were selected for spectral analysis. The results showed that the boundary line of the water after fusion appeard more clearly, which helped to refine the water extraction boundary.

Key words : fusion algorithm; Landsat 8; multi-spectral and panchromatic; spectral analysis

0 引言

圖像數(shù)據(jù)融合是一個對多傳感器的圖像數(shù)據(jù)和其他信息的處理過程,，重點是按照一定的規(guī)則或算法對那些在時空上冗余或互補的多元數(shù)據(jù)作運算處理,，獲取較單一數(shù)據(jù)更精確,、豐富的信息,，生成具有新的空間,、波譜時間特征的合成圖像。按照融合水平從低到高,，圖像數(shù)據(jù)融合可在像素,、特征,、決策層三個層次上進行^［1］,。目前基于像素的融合應(yīng)用較為廣泛，雖然涉及的數(shù)據(jù)處理量大,，但它基于最原始的圖像數(shù)據(jù),，能更多地保留圖像信息,，具有最高的精度^［2］。常用的像素級圖像融合算法有HIS融合法,、主成分變換融合法,、Brovery融合法,、Gram-schidt融合法,、小波變換融合法、NNDiffuse融合法等,。針對這些融合算法，一些學(xué)者利用Quickbird^［3］,、Landsat 7 ETM+^［4］,、SPOT6^［5］,、高分一號^［6］,、高分二號^［7］等多分辨率傳感器的多光譜與全色影像融合進行了研究,，取得了顯著的成果,，但對 Landsat 8 衛(wèi)星影像開展的研究較少,。

Landsat 8是由NASA發(fā)射于2013年的新型遙感衛(wèi)星,，集成了兩個傳感器，即陸地成像儀（OLI）和熱紅外傳感器（TIRS）,。OLI傳感器共有9個波段,，為了與Landsat 7 ETM+波段保持?jǐn)?shù)據(jù)一致性，Landsat 8的幾個反射波段和ETM+的反射波段是相同的或非常接近的,。然而，也有若干OLI波段的寬度經(jīng)過了改善,，以減輕或避免在ETM+波段中出現(xiàn)的各種大氣吸收特征的影響,。例如OLI全色波段8被縮小，以提高植被和非植被區(qū)域的對比度^［8］,。本文以Landsat 8數(shù)據(jù)為例,，采用3種不同的算法對其多光譜和全色數(shù)據(jù)進行融合研究，并從定性和定量兩個方面分析比較這3種融合算法,。對最佳融合算法,，選取了水體邊界線處的影像像元進行光譜分析，進一步研究融合效果,。

1 影像融合方法和原理

本文主要利用了傳統(tǒng)的Brovery融合法,、NNDiffuse融合法和Gram-Schmidt融合法，融合原理簡要介紹如下,。

1.1 Brovery變換^［1］

Brovery變換是通過對多波段數(shù)據(jù)進行歸一化處理，再乘以高分辨率數(shù)據(jù),，如式（1）所示,。在本次研究中,，選取Landsat 8 數(shù)據(jù)中的中紅外,、近紅外、綠波段為R,、G、B通道,，選取同一時相的全色波段為高分辨率波段進行數(shù)據(jù)融合。

微信截圖_20180911163546.png

式中：DN_1~DN₃為多波段圖像數(shù)據(jù)值,；DNh為高分辨率圖像數(shù)據(jù)值,；DN_f1~DN_f3為結(jié)果值,。

1.2 NNDiffuse變換^［9］

NNDiffuse Pan Sharpening（Nearst Neighbor Diffusion Pan Sharpening）圖像融合算法由美國羅徹斯特理工學(xué)院（RIT）SUN W等人最先提出,。該方法首先建立低分辨率多波段數(shù)據(jù)與重采樣后全色波段間（重采樣后分辨率與多波段相同）的線性響應(yīng)向量T,，建立9個興趣像元與超像素區(qū)分布計算全色波段的像元差異系數(shù)N,，結(jié)合差異系數(shù)N與多波段數(shù)據(jù)建立高分辨率多波段數(shù)據(jù),。ENVI5.2及其以后版本支持該融合算法,。它支持眾多傳感器類型，如Landsat 8,、SPOT、WorldView-2/3及國產(chǎn)衛(wèi)星等，融合結(jié)果能很好地保留色彩,、紋理和光譜信息,。

1.3 Gram-Schmidt 變換^［10］

Gram-Schmidt變換是統(tǒng)計學(xué)中常用的一種多維線性正交變換,，采用GS變換對遙感影像多維數(shù)據(jù)進行正交化處理,，可以有效地去除相鄰波段間較強的相關(guān)性,，最大程度地消除影像的信息冗余,。該變換首先使用多光譜低空間分辨率影像模擬出一個全色高分辨率影像,；隨后，以模擬的全色高分辨率影像作為GS變換的第一個分量對低分辨率影像進行GS變換；再次,，將高分辨率影像與GS變換后的第一分量進行匹配,，以產(chǎn)生經(jīng)過修改的高分辨率影像,，用經(jīng)過修改的高分辨率影像替換GS變換后的第一個分量,，產(chǎn)生一個新的數(shù)據(jù)集；最后,，將新的數(shù)據(jù)集進行GS逆變換即可產(chǎn)生空間分辨率增強的多光譜影像,。該算法可以對任意波段數(shù)的圖像進行融合，且產(chǎn)生的融合影像保持了低空間分辨率波段的光譜特性,，信息失真小,。

2 影像融合質(zhì)量評價

2.1 主觀評價

主觀評價主要觀察圖像的清晰度、色調(diào)以及地物的文理信息等方面的因素,，屬于定性評價,。

2.2 定量評價

客觀評價的指標(biāo)主要從信息量、清晰度和逼真度三個角度進行評價,。圖像融合的目的之一在于增加圖像的有效信息量,，熵和標(biāo)準(zhǔn)差可以反映圖像的信息量；影像的清晰度是指其邊界或影線兩側(cè)附近灰度有明顯區(qū)別,，這種區(qū)別可以用梯度表示,；影像的逼真度是指融合圖像與原始圖像的偏離程度，本文選用偏差表示,。假設(shè)影像大小為m×n,，影像灰度變化范圍為(0,255)，以下為各指標(biāo)的計算和評價方法,。

2.2.1 均值

均值在目視上表現(xiàn)為影像的像素平均值，通俗地說就是平均亮度,。融合后影像像素平均值變化越小,，則融合后影像的平均亮度與原始影像越接近，該融合方法的光譜保真性就越好^［11］,。計算公式為：

微信截圖_20180911164730.png

式中,，微信截圖_20180911164204.png 為影像的灰度平均值；m為影像的行數(shù)；n為影像的列數(shù),；M(x,y)為影像中原影像坐標(biāo)為(x,y)的像元灰度值,。

2.2.2 標(biāo)準(zhǔn)差S

標(biāo)準(zhǔn)差指的是影像各個像元的灰度值相對于全景像元平均值的離散程度。標(biāo)準(zhǔn)差越大,，影像反差就越大,，也越容易區(qū)分地物，有利于提取影像的信息,。計算公式為：

微信截圖_20180911164833.png

式中,，微信截圖_20180911164204.png 為影像的灰度平均值；m為影像的行數(shù),；n為影像的列數(shù),；M(x,y)為影像中原影像坐標(biāo)為(x,y)的像元灰度值。

2.2.3 信息熵Hx

信息熵由香農(nóng)于1948年提出,，可用來衡量影像信息豐富程度,。熵越大，說明影像包含的信息越多,，融合效果越好^［12］,。根據(jù)香農(nóng)的信息論，圖像的信息熵可按下式進行計算：

微信截圖_20180911165333.png

式中,，Pi為灰度值為i的概率,，x為像元編號。Pi通過Quick Stats直接導(dǎo)出,。

2.2.4 偏差 D

偏差反映的是影像處理前后的變化程度,，該值越大，圖像畸變越大,，圖像信息損失越多^［13］,。偏差即融合前后每個波段影像對應(yīng)像元的差值的平均值，按式(5)進行計算：

微信截圖_20180911165632.png

式中,，M(x,y)為原影像坐標(biāo)為(x,y)的像元灰度值,；F(x,y)為融合影像坐標(biāo)為(x,y)的像元灰度值。

2.2.5 平均梯度 g

平均梯度是通過計算融合影像相鄰像元間的插值來反映影像對微小細(xì)節(jié)的反差,，對于兩地類的交界處相鄰像元差值較大,。研究融合后影像橫向與縱向兩方向的相鄰像元差值，進而評價其清晰度,。平均梯度按式(6)進行計算：

微信截圖_20180911170018.png

式中,，ΔFx(x,y)、ΔFy(x,y)分別為x軸和y軸方向的相鄰像元灰度值差值,。

3 融合效果評價

以遼寧省撫順縣的大伙房水庫為實驗區(qū),，選取成像于2015年9月27日的Landsat 8 影像,，對經(jīng)過預(yù)處理的多光譜波段和全色波段采用上述3種方法進行融合。經(jīng)大氣處理后的多光譜影像值為反射率值的10 000倍,，定標(biāo)后的全色波段數(shù)據(jù)為輻射亮度值,。因不同的融合前后數(shù)據(jù)單位及數(shù)據(jù)類型不同，導(dǎo)致均值,、標(biāo)準(zhǔn)差等指標(biāo)不具有可比性,。為了使3種融合結(jié)果具有可比性，首先按式(7)將融合影像像元值拉伸到0~255之間,。同時利用三次采樣法將多光譜影像的空間分辨率采樣至15 m,，該過程通過ENVI中的Layer stack模塊實現(xiàn)。

微信截圖_20180911170202.png

式中,，Bi為融合后波段像元值,；min、max分別為該波段影像的最小與最大值,；B′i為將融合影像拉伸后的值,。

3.1 主觀評價

主觀評價主要觀察圖像的清晰度、色調(diào)以及地物的文理信息等方面的因素,，圖1所示為原始影像和不同融合方法處理后的影像對比,。

微信截圖_20180911170319.png

圖 1 原始影像與3種不同的融合算法的融合影像

其中圖1（a）、（c）,、（d）采用R,、G、B通道選擇紅綠藍(lán)波段,，圖1（b）采用SWIR,、Red、Green顯示,。從清晰度角度分析,，融合后影像紋理更加清晰，水體邊界更容易辨別,；從色調(diào)方面觀察,，Brovery融合和NNDiffuse融合使得植被、裸地,、水體比融合前色彩更為明亮,，而Gram-Schmidt融合與原影像的色差最小。

3.2 客觀評價

為了客觀,、定量地評價融合結(jié)果,，計算了3種融合算法獲得的融合影像在各個波段的均值、標(biāo)準(zhǔn)差,、信息嫡,、偏差、平均梯度,，詳細(xì)的評價指標(biāo)如表1,、表2、表3所示,。

表 1 Brovery 融合法評價結(jié)果

微信截圖_20180911170501.png

表 2 NNDiffuse 融合法評價結(jié)果

微信截圖_20180911170511.png

表 3 Gram-Schmidt 融合法評價結(jié)果

微信截圖_20180911170522.png

從信息熵看,，Gram-Schmidt融合法和NNDiffuse融合法的值大于Brovery融合法，說明前兩種融合方法包含有更豐富的光譜信息,；對比標(biāo)準(zhǔn)差,，任一種融合方法相比于其他兩種融合方法均沒有明顯優(yōu)勢；Gram-Schmidt融合法與原圖像的偏差小于其他兩種方法,，說明其光譜畸變較小,，符合主觀評價結(jié)果；Brovery融合法的平均梯度優(yōu)于Gram-Schmidt融合法和NNDiffuse融合法,；相比于其他兩種融合方法,，Brovery融合法必須也只能同時對三個波段進行融合處理，會增加處理難度,。

綜上所述,，Gram-Schmidt融合法擁有更加豐富的光譜信息、較高的光譜保真和清晰度,。為進一步研究融合算法的融合效果,，本文對水體邊界線處的影像像元進行光譜分析，檢驗其對水體邊界提取的精細(xì)化程度,。Gram-Schmidt融合前后像元的位置如圖2所示,。圖中，邊界細(xì)線為利用MNDWI指數(shù)將閾值設(shè)為1所提取的水體邊界線,。

微信截圖_20180911171253.png

圖 2 融合前后影像像元空間分布

對該部分像元進行編號分析如圖3所示,。圖中，灰色像元為融合后15 m空間分辨率像元,，虛線表示的為融合前30 m空間分辨率像元,。

融合前后對應(yīng)影像位置的像元反射率光譜特征曲線如圖4和圖5所示。根據(jù)融合前像元的水體和植被的反射率波譜特征,，可將a,、b像元判定為水體，c,、d像元判定為植被,。

微信截圖_20180911171301.png

圖 3 融合前后像元空間分布分析圖

微信截圖_20180911171309.png

圖 4 融合前像元反射率光譜特征

微信截圖_20180911171317.png

圖 5 融合后各像元反射率光譜特征

對融合后的像元分析發(fā)現(xiàn)，2,、3像元具有明顯的水體波譜特征,；4,、7像元符合植被光譜特征；5,、8像元與植物光譜特征類似,，但是其數(shù)值與植被相差較大；其余像元均未表現(xiàn)出明顯符合某種光譜特征的屬性,，這些像元融合全色影像的光譜特征,，近紅外和中紅外的反射率高于水體，但又遠(yuǎn)低于植被,。利用MNDWI指數(shù)法對影像進行提取發(fā)現(xiàn),，1、2,、3,、5、6,、9像元被歸為水體,。通過光譜分析，相較于融合前的邊界線,，融合后的水體邊界線也更為清晰地呈現(xiàn)出來,。

4 結(jié)論

通過主觀評價和客觀評價兩個角度選取融合算法，使經(jīng)過融合的影像具有全色圖像數(shù)據(jù)的高空間分辨率和多光譜圖像數(shù)據(jù)的多光譜特征,。對NNDiffuse,、Gram-Schmidt、Brovery 3種融合算法進行比較發(fā)現(xiàn),，Gram-Schmidt融合算法具有更高的光譜保真度,，在標(biāo)準(zhǔn)差、信息熵,、偏差等客觀指標(biāo)上優(yōu)于另外兩種算法,。對水體邊界處的像元進行光譜分析顯示，融合后的影像可以提高水體邊界提取的準(zhǔn)確程度,。

參考文獻

［1］趙英時.遙感應(yīng)用分析原理與方法［M］.北京: 科學(xué)出版社,2013.

［2］花利忠,鄧琪琪,張敏莉,，等.基于Landsat 8 OLI遙感影像的融合算法比較［J］.廈門理工學(xué)院學(xué)報,2017,25(1): 5-11.

［3］譚永生,沈掌泉,賈春燕，等. QuickBird全色與多光譜影像融合方法比較研究［J］. 科技通報,2008,24(4):498-503.

［4］徐涵秋.Landsat 7 ETM+影像的融合和自動分類研究［J］.遙感學(xué)報,2005,9(2):186-194.

［5］郭蕾,楊冀紅,史良樹,，等. SPOT6遙感圖像融合方法比較研究［J］. 國土資源遙感, 2014, 26(4): 71-77.

［6］劉錕,付晶瑩,李飛. 高分一號衛(wèi)星4種融合方法評價［J］. 遙感技術(shù)與應(yīng)用, 2015,30(5): 160-166.

［7］薛晶,官云蘭,李宜龍,，等.“高分二號”衛(wèi)星影像融合方法對比和評價［J］. 航天返回與遙感,2017,38(2):91-99.

［8］ THOMAS L, RALPH W. K, JONATHAN C. 遙感與圖像解譯（原書第7版）［M］.北京: 電子工業(yè)出版社,2016.

［9］ SUN W,MESSINGER D. Nearest-neighbor diffusion-based pan-sharpening algorithm for spectral images［J］.Optical Engineering, 2014, 53(1):013107.

［10］楊麗萍,夏敦勝,陳發(fā)虎. Landsat 7 ETM+全色與多光譜數(shù)據(jù)融合算法的比較［J］. 蘭州大學(xué)學(xué)報(自然科學(xué)版), 2007,43(4):7-11, 17.

［11］李瑩瑩,洪金益,張建國,等. IKONOS影像融合方法對比研究及應(yīng)用效果評價［J］. 礦產(chǎn)勘查,2014,5(3):499-504.