朱昕宇， 孫明偉， 王樹(shù)根， 余 磊， 郭 清

(武漢大學(xué) 遙感信息工程學(xué)院， 武漢 430079)

?

衛(wèi)星影像反立體改正算法比較探究

朱昕宇*， 孫明偉， 王樹(shù)根， 余 磊， 郭 清

(武漢大學(xué) 遙感信息工程學(xué)院， 武漢 430079)

太陽(yáng)同步衛(wèi)星成像時(shí)的光照方向與人類視覺(jué)習(xí)慣相反，導(dǎo)致影像中出現(xiàn)地形反轉(zhuǎn)的感知現(xiàn)象，即反立體現(xiàn)象.現(xiàn)有的一些方法能夠改正這種嚴(yán)重誤導(dǎo)影像解譯的現(xiàn)象，但針對(duì)算法的系統(tǒng)性比較研究較為缺乏.以SPOT-5與WorldView-2衛(wèi)星影像為例，從反立體改正程度、空間細(xì)節(jié)的保持、光譜特性的保持方面選取幾種代表性算法進(jìn)行比較，進(jìn)而歸納出各個(gè)算法的特點(diǎn)及適用情況，為實(shí)際應(yīng)用打下基礎(chǔ).

反立體； 衛(wèi)星影像； 空間細(xì)節(jié)； 光譜特性

由于可見(jiàn)光遙感衛(wèi)星影像能夠直觀地呈現(xiàn)地球表面概況，并且具有時(shí)相多、成像范圍廣的特點(diǎn)，在土地規(guī)劃、能源管理等領(lǐng)域有廣泛的應(yīng)用.多數(shù)可見(jiàn)光遙感衛(wèi)星與太陽(yáng)同步，于當(dāng)?shù)貢r(shí)間的上午成像，此時(shí)光線來(lái)自影像右下方、陰影落在左上方[1-2]；人類大腦認(rèn)為自然界存在唯一光源，且光源來(lái)自上方、陰影落在下方[3].衛(wèi)星成像時(shí)的光照方向與人類視覺(jué)感知習(xí)慣相反，導(dǎo)致人在深度感知上出現(xiàn)錯(cuò)誤，將山谷解譯為山脊，反之亦然.這種地形反轉(zhuǎn)的現(xiàn)象即為可見(jiàn)光遙感影像的反立體現(xiàn)象，也稱為虛假地形感知現(xiàn)象[4].以谷歌地球發(fā)布的藏南北部地區(qū)喜馬拉雅山部分影像(圖1)為例，影像中南北向河流在山脊流淌，明顯與人類認(rèn)知習(xí)慣不符.將影像與相應(yīng)數(shù)字高程模型(Digital Elevation Model， DEM)或地形圖比較時(shí)，即可發(fā)現(xiàn)反立體現(xiàn)象.

雖然通過(guò)文字標(biāo)注或參考地形信息的方式可以提示反立體現(xiàn)象的存在，但直接對(duì)影像進(jìn)行改正能夠提供更直觀的目視效果，進(jìn)而提升遙感影像地圖的產(chǎn)品質(zhì)量.現(xiàn)有改正方法大致可分為以下幾類：基于觀察者視角的改正方法，如旋轉(zhuǎn)法[5]；基于地表特征的改正方法，主要包括灰度反轉(zhuǎn)法[5]、亮度分量旋轉(zhuǎn)法[1]、SRM(Shaded Relief Model)替代法[2]、小波變換融合法[6]、SRM疊加法[7]及SRM乘法融合法[8]；基于光照方向的改正方法，以純地形糾正法[8]為代表；人工改正方法，如使用圖像處理軟件的人工改正[9]等.

圖1 谷歌地球喜馬拉雅山部分影像Fig.1 Image of part of the Himalayas from Google Earth

現(xiàn)階段存在一些反立體改正方法，但目前缺乏針對(duì)已有算法的系統(tǒng)性比較研究，為實(shí)際應(yīng)用造成了障礙.因此，全面了解各個(gè)算法的特點(diǎn)與適用情況具有重要意義.本文以SPOT-5與WorldView-2影像為數(shù)據(jù)源，對(duì)幾種代表性的反立體改正方法進(jìn)行測(cè)試，從反立體改正程度、原始影像空間細(xì)節(jié)及光譜特性的保持3個(gè)方面進(jìn)行比較，進(jìn)而歸納出各個(gè)算法的特點(diǎn)及適用情況，為影像解譯及進(jìn)一步分析打下基礎(chǔ).

1 反立體改正方法

鑒于旋轉(zhuǎn)法、灰度反轉(zhuǎn)法與人工處理法的局限性，本文選擇以下5種算法進(jìn)行對(duì)比與分析：SRM替代法、亮度分量反轉(zhuǎn)法、SRM乘法融合法、純地形糾正法、SRM疊加法.

1.1 SRM替代法

Saraf等人[2]最早將SRM引入反立體改正過(guò)程，提出一種基于SRM的改正算法.SRM表示地形在特定太陽(yáng)高度及方位下的陰影分布情況，在一定程度上可以為反立體改正提供參考.該方法首先將原始影像進(jìn)行HIS變換，得到原始影像的亮度分量I0；然后將西北向光照下的SRM代替I0，其中SRM的太陽(yáng)高度角與原始影像一致，太陽(yáng)方位角在原始高度角的基礎(chǔ)上加180°獲得；最后將SRM與原始影像的飽和度及色調(diào)分量進(jìn)行HIS逆變換，獲得反立體改正后影像.

1.2 亮度分量反轉(zhuǎn)法

亮度分量反轉(zhuǎn)法[1]通過(guò)反轉(zhuǎn)影像HIS變換后的亮度分量來(lái)獲得正確立體視覺(jué)效果的影像.算法首先將原始影像的亮度分量I0做灰度反轉(zhuǎn)處理，得到反轉(zhuǎn)后亮度I1.由于直接反轉(zhuǎn)亮度往往會(huì)造成亮度增加，將I1與I0進(jìn)行直方圖匹配，得到新的亮度分量I2，再將其與原始影像的飽和度及色調(diào)分量進(jìn)行HIS逆變換.處理后影像往往會(huì)產(chǎn)生較為嚴(yán)重的色彩偏差，需要進(jìn)行色彩增強(qiáng)處理，本文中采用Photoshop軟件的顏色匹配功能進(jìn)行處理.

1.3 SRM乘法融合法

SRM乘法融合法[8]假設(shè)地表為理想的朗伯體，即地物表面在各個(gè)方向反射等量的光照.算法的核心為將西北光照方向下的SRM與原始影像相乘.由于SRM本質(zhì)上為特定光照方向下的太陽(yáng)入射角余弦，SRM乘法融合法相當(dāng)于將西北光照下的太陽(yáng)入射角與原始影像DN值相乘.通過(guò)引入代表陰影分布的SRM，改變影像中光照與陰影的分布情況，進(jìn)而改正影像中的反立體現(xiàn)象.

1.4 純地形糾正法

地形糾正，也稱為地形歸一化或地形光照改正，是指通過(guò)建立一系列模型來(lái)改正由地形起伏引起的地表光照不均的現(xiàn)象.純地形糾正法[8]引入地形糾正中的cosine模型，將影像的光照方向由原始東南向改為垂直入射，再由垂直入射改為西北向入射.算法主要涉及的公式如下：

(1)

其中，LT1為原始影像輻射值；LT2為反立體改正后影像輻射值；i1為東南光照下的太陽(yáng)入射角；i2為西北向光照下的太陽(yáng)入射角.

1.5 SRM疊加法

Gil等人[7]提出，可以通過(guò)在原始影像上疊加西北向光照下SRM來(lái)改正反立體現(xiàn)象.通過(guò)對(duì)西北光照下的SRM賦以一定的透明度并疊加至原始影像，即可獲得正確立體解譯的影像.其中，SRM的太陽(yáng)方位角為315°；對(duì)于平坦地面，太陽(yáng)高度角為45°，否則為60°；對(duì)于全色、多光譜、融合后影像而言，SRM的透明度應(yīng)分別設(shè)置為50%、70%與55%.

2 實(shí)驗(yàn)結(jié)果

采用SPOT-5與WorldView-2衛(wèi)星影像對(duì)算法進(jìn)行測(cè)試與評(píng)估.影像均選自廣東省山區(qū)，部分參數(shù)如表1所示.從表中可知，影像獲取時(shí)太陽(yáng)方位角處于135°與180°之間，即太陽(yáng)位于東南偏南方向，在光照方向上符合反立體形成的條件.實(shí)驗(yàn)所用DEM為廣東省1∶10000基礎(chǔ)測(cè)繪成果.將經(jīng)過(guò)外業(yè)控制點(diǎn)精確定向后的影像利用DEM進(jìn)行正射糾正處理，得到與DEM匹配的正射影像.將影像分別利用上述5種方法進(jìn)行反立體改正，結(jié)果參見(jiàn)圖2、圖3.

表1 實(shí)驗(yàn)數(shù)據(jù)參數(shù)

3 比較與分析

在上述實(shí)驗(yàn)結(jié)果的基礎(chǔ)上，從反立體改正程度、空間細(xì)節(jié)的保持、光譜特性的保持3個(gè)方面比較與分析每種算法的改正結(jié)果；然后從原理出發(fā)，進(jìn)一步對(duì)算法的性能進(jìn)行衡量.

3.1 算法性能比較

3.1.1 反立體改正程度 由圖2、圖3可知，幾種方法對(duì)于反立體現(xiàn)象均有不同程度的改正.就反立體改正程度而言，SRM替代法、SRM乘法融合法與純地形糾正法性能較優(yōu)，改正前后影像中出現(xiàn)明顯的立體反轉(zhuǎn)，且改正后影像的立體效果有所提升；其次為亮度分量反轉(zhuǎn)法與SRM疊加法，改正后影像中存在原始光照與陰影位置關(guān)系的變化，但影像立體感較弱.

圖2 SPOT-5影像改正結(jié)果Fig.2 Experimental results of SPOT-5 image

圖3 WorldView-2改正結(jié)果Fig.3 Experimental results of WorldView-2 image

對(duì)反立體改正程度的定量衡量采用影像與西北向光照下SRM的相關(guān)系數(shù)指標(biāo)[8，10].若相關(guān)系數(shù)為正，表明影像中不存在反立體現(xiàn)象，反之則相反.圖4為不同方法的計(jì)算結(jié)果，在一定程度上驗(yàn)證了上述目視判斷：原始影像與西北向SRM的相關(guān)系數(shù)為負(fù)，影像中存在反立體效應(yīng)；SRM替代法影像相關(guān)系數(shù)最高，立體效果最為明顯；SRM乘法融合法、純地形糾正法與SRM疊加法其次；亮度分量反轉(zhuǎn)法對(duì)反立體的改正程度較弱.

3.1.2 空間細(xì)節(jié)的保持 由于缺乏衡量影像間空間細(xì)節(jié)相關(guān)程度的定量指標(biāo)，本文采用目視判讀的方法比較算法性能.為了便于觀察及說(shuō)明，以SPOT-5影像為例，選擇四種地物類型進(jìn)行細(xì)節(jié)對(duì)比，分別為植被、裸地、水體及村落，如圖5所示.由圖可知，就空間細(xì)節(jié)信息的保持而言，亮度分量反轉(zhuǎn)法、SRM乘法融合法、純地形糾正法性能較優(yōu)，能夠在改正反立體的同時(shí)最大程度保留原始影像的細(xì)節(jié)信息；SRM疊加法性能其次，損失了影像中部分細(xì)節(jié)；SRM替代法效果較差，4種地物的空間細(xì)節(jié)信息均遭到損失，其中水體及村落尤為嚴(yán)重，甚至達(dá)到難以辨認(rèn)的程度.

圖4 影像與西北向SRM相關(guān)系數(shù)Fig.4 Correlation coefficient between image and SRM in the NW direction

圖5 不同算法細(xì)節(jié)對(duì)比圖(以SPOT-5數(shù)據(jù)為例)Fig.5 Details of the experimental results (taking SPOT-5 data as an example)

3.1.3 光譜特性的保持 對(duì)算法光譜特性保持的比較從整體及局部?jī)煞矫孢M(jìn)行.就整體而言，與原始影像相比，幾種算法改正后影像均存在一定程度的偏色，如圖2、圖3所示.其中，SRM乘法融合法性能最優(yōu)，改正前后影像整體色調(diào)基本一致，但亮度略有下降；其次為SRM疊加法，改正后影像的整體色調(diào)略有偏差，飽和度和對(duì)比度有所減少；純地形糾正法影像在向陽(yáng)面偏色程度較大；SRM替代法及亮度分量反轉(zhuǎn)法的效果較差，地物光譜信息發(fā)生了嚴(yán)重扭曲.

對(duì)于局部地物光譜特性保持能力的比較仍依據(jù)圖5進(jìn)行說(shuō)明.針對(duì)不同地物的比較進(jìn)一步驗(yàn)證了上述判斷：SRM乘法融合法對(duì)地物光譜信息的保持能力較強(qiáng)，改正后各地物的整體偏色較小，但亮度略有降低；SRM疊加法改正影像中的地物色調(diào)基本保持不變，但飽和度有所降低；純地形糾正法對(duì)于裸地、水體及村落的光譜特性有一定保持能力，但對(duì)原始影像背陽(yáng)面植被的改正效果較差；四種地物在SRM替代法與亮度分量反轉(zhuǎn)法改正后影像中均發(fā)生了嚴(yán)重偏色.

3.2 分析與討論

SRM替代法將西北向光照下的地形陰影模型代替原始影像HIS變換后的亮度分量進(jìn)行反立體改正.其中，SRM代表特定地形在特定光照條件下的光照及陰影分布，僅考慮反射光，代表部分地形信息；而原始影像的亮度分量不僅包含由地形引起的亮度信息，還包括由地物引起的亮度信息，如地物的紋理、結(jié)構(gòu)等.將僅代表部分地形的亮度代替原始影像的亮度，會(huì)不可避免地造成大部分空間信息及一部分光譜信息的損失.

亮度分量反轉(zhuǎn)法直接反轉(zhuǎn)原始影像的亮度分量，從原理上保持了影像的空間細(xì)節(jié)信息.但是，亮度反轉(zhuǎn)操作會(huì)造成地物光譜特性的嚴(yán)重扭曲，即使引入直方圖匹配進(jìn)行顏色調(diào)整，也無(wú)法彌補(bǔ)色彩信息的損失.原始影像中亮度較高的裸地、村落等地物，經(jīng)亮度反轉(zhuǎn)后失去了原始光譜結(jié)構(gòu)，在改正后影像中偏暗；原始影像中亮度較低的植被、水體等地物，經(jīng)反轉(zhuǎn)后亮度增加，導(dǎo)致反立體改正后影像偏色.

SRM乘法融合法將西北向光照下SRM與原始影像灰度值相乘.由于SRM為西北向光照下太陽(yáng)入射角的余弦值，其范圍在0～1之間，與原始影像灰度值相乘會(huì)造成改正后影像整體變暗.雖然改正后影像能夠保留大部分空間信息及光譜特性，但算法在改變光照方向的同時(shí)減少了光照強(qiáng)度，改正后影像能夠滿足目視需求卻不能應(yīng)用于后續(xù)定量處理.

純地形糾正法以地形糾正后影像為中介，將光照方向由東南向改為垂直入射，再改為西北向入射.由于地形糾正可以較大程度地保留影像的空間細(xì)節(jié)和光譜特性，此方法能夠從原理上保留原始影像的大部分信息.但是，地形糾正中采用的cosine模型在影像亮度較低時(shí)往往會(huì)造成地形過(guò)度改正[11-12]，導(dǎo)致原始影像背陽(yáng)面地物的光譜信息發(fā)生嚴(yán)重扭曲.此外，cosine模型對(duì)于由附近高大地物遮擋而產(chǎn)生的陰影無(wú)效，進(jìn)而限制了純地形糾正法的適用范圍.

SRM疊加法將半透明的西北向SRM與原始影像疊加，是一種在改正反立體現(xiàn)象的同時(shí)盡量保留原始影像其它信息的折衷方法.影像半透明化與疊加的過(guò)程中，反立體現(xiàn)象能夠在一定程度上得到改正，但影像的細(xì)節(jié)及光譜信息會(huì)遭到損失.此外，SRM透明度的選取至關(guān)重要，較高的透明度對(duì)反立體的改正效果不明顯，較低的透明度則會(huì)犧牲原始影像的空間信息與光譜特性.

4 結(jié)論

本文通過(guò)對(duì)比SRM替代法、亮度分量反轉(zhuǎn)法、SRM乘法融合法、純地形糾正法、SRM疊加法，總結(jié)每種方法的特點(diǎn)及適用范圍如下．

1)SRM替代法原理簡(jiǎn)單，改正后影像立體效果明顯.但是，原始影像的空間細(xì)節(jié)與光譜特性均會(huì)產(chǎn)生嚴(yán)重?fù)p失，且算法需要與影像同分辨率的DEM，對(duì)數(shù)據(jù)要求較高.因此，算法適用于僅對(duì)反立體改正程度有較高要求、且相應(yīng)數(shù)據(jù)充足的情況；

2) 亮度分量反轉(zhuǎn)法簡(jiǎn)單易行，能夠從原理上保持原始影像的空間細(xì)節(jié)，對(duì)數(shù)據(jù)要求低.但是，算法會(huì)造成地物光譜特性的扭曲，且改正后影像的立體效果較弱.因此，算法在僅對(duì)空間細(xì)節(jié)要求較高，且相應(yīng)數(shù)據(jù)不足的情況下適用.

3)SRM乘法融合法整體效果最優(yōu)，能夠在改正反立體的同時(shí)兼顧原始影像的空間細(xì)節(jié)及光譜特性，且算法易于實(shí)現(xiàn).但是，改正后影像亮度降低，且算法對(duì)數(shù)據(jù)要求較高.因此，當(dāng)相應(yīng)數(shù)據(jù)充足，且對(duì)改正后影像的目視效果要求較高時(shí)，可選擇此算法進(jìn)行反立體改正.

4) 純地形糾正法能夠從原理上保持影像的空間細(xì)節(jié)及光譜特性，但算法中引入的cosine校正模型易引起過(guò)度地形改正，導(dǎo)致改正后影像偏色嚴(yán)重，且算法對(duì)于由遮擋而產(chǎn)生的陰影無(wú)效.此外，較高的數(shù)據(jù)要求問(wèn)題仍無(wú)法避免.因此，算法適用于相應(yīng)數(shù)據(jù)充足，對(duì)空間細(xì)節(jié)有一定要求、色彩信息要求不高，且地形較平坦的情況.

5)SRM疊加法能夠保留影像的部分空間細(xì)節(jié)及光譜特性，但影像疊加的原理使得算法無(wú)法在改正反立體的同時(shí)兼顧其它信息，且對(duì)數(shù)據(jù)要求較高.因此，算法適用于對(duì)反立體改正程度、空間細(xì)節(jié)及光譜特性有一定要求，且相應(yīng)數(shù)據(jù)充足的情況.

此外，由于反立體改正后影像主要用于輔助目視解譯，且?guī)追N改正算法均會(huì)造成原始影像空間細(xì)節(jié)與光譜特性的損失， 反立體改正后影像不適用于后續(xù)定量處理，如土地覆蓋分類、數(shù)字高程模型提取等.

[1]SarafAK，SinhaST，GhoshP，etal.Anewtechniquetoremovefalsetopographicperceptionphenomenonanditsimpactsinimageinterpretation[J].InternationalJournalofRemoteSensing， 2007， 28(5): 811-821．

[2]SarafAK，GhoshP，SarmaB，etal.Cover:Developmentofanewimagecorrectiontechniquetoremovefalsetopographicperceptionphenomena[J].InternationalJournalofRemoteSensing， 2005， 26(8): 1523-1529．

[3]RamachandranVS.Perceptionofshapefromshading[J].Nature， 1988， 331(6152): 163-166．

[4]Bernabé-PovedaMA， ??ltekinA.Prevalenceoftheterrainreversaleffectinsatelliteimagery[J].InternationalJournalofDigitalEarth， 2015, 8(8): 640-655.

[5]SarafAK，DasJD，AgarwalB，etal.Falsetopographyperceptionphenomenaanditscorrection[J].InternationalJournalofRemoteSensing， 1996， 17(18): 3725-3733．

[6]WuB，LiH，GaoY.InvestigationandremediationoffalsetopographicperceptionphenomenaobservedonChang’E-1lunarimagery[J].PlanetaryandSpaceScience， 2013， 75: 158-166．

[7]GilML，ArzaM，OrtizJ，etal.DEMshadingmethodforthecorrectionofpseudoscopiceffectonmulti-platformsatelliteimagery[J].GIScience&RemoteSensing， 2014， 51(6): 630-643．

[8] 章皖秋， 袁 華， 徐天蜀. 遙感影像正立體化研究[J]. 北京林業(yè)大學(xué)學(xué)報(bào)， 2010 (3): 41-46．

[9]PattersonT，KelsoNV.HalSheltonrevisited:Designingandproducingnatural-colormapswithsatellitelandcoverdata[J].CartographicPerspectives， 2004 (47): 28-55．

[10] 陳建珍， 何 超， 岳彩榮. 基于反DEM和HSV變換的遙感圖像反立體校正[J]. 東北林業(yè)大學(xué)學(xué)報(bào)， 2011， 39(3): 125-126．

[11]MeyerP，IttenKI，KellenbergerT，etal.RadiometriccorrectionsoftopographicallyinducedeffectsonLandsatTMdatainanalpineenvironment[J].ISPRSJournalofPhotogrammetryandRemoteSensing， 1993， 48(4): 17-28．

[12]SmithJA，LinTL，RansonKL.TheLambertianassumptionandLandsatdata[J].PhotogrammetricEngineeringandRemoteSensing， 1980， 46(9): 1183-1189．

Comparison of terrain reversal effect correction methods in satellite imagery

ZHU Xinyu， SUN Mingwei， WANG Shugen， YU Lei， GUO Qing

(School of Remote Sensing and Information Engineering, Wuhan University， Wuhan 430079)

Most sun-synchronous satellites acquire images with a SE illumination source， which goes against human visual habits. Hence， when observing such images， the terrain reversal effect occurs particularly in large-relief regions. Many methods have been developed to solve the problem which is seriously misleading imagery interpretation. However， no systematic comparisons have been made. Experiments were conducted on SPOT-5 and WorldView-2 satellite images， and comparisons of several typical algorithms were made through three aspects: correction extent， preservation of spatial details， and preservation of spectral characteristics. The feature and application situation of each approach were further concluded for practical application．

terrain reversal effect; satellite imagery; spatial details; spectral characteristics

2015-04-19.

國(guó)家973計(jì)劃項(xiàng)目(2012CB719904);國(guó)家自然科學(xué)基金項(xiàng)目(41301519).

1000-1190(2015)05-0805-06

TP751.1

A

*E-mail: zhuxinyu@whu.edu.cn.