宋麗瑤，張 浩，林 峰

(1. 中國石油大學(華東)，山東 青島 266555； 2. 中國科學院遙感與數(shù)字地球研究所，北京 100094； 3. 青島海洋地質(zhì)研究所，山東 青島 266555)

陰影區(qū)提取對高分辨率影像地形輻射校正的影響與分析

宋麗瑤1，張 浩2，林 峰3

(1. 中國石油大學(華東)，山東 青島 266555； 2. 中國科學院遙感與數(shù)字地球研究所，北京 100094； 3. 青島海洋地質(zhì)研究所，山東 青島 266555)

利用輻射傳輸機理對復雜山區(qū)進行地形輻射校正，可以全面地考慮太陽直接輻射亮度、天空散射輻射亮度及鄰近地表反射輻射亮度3部分因素，從而取得較好的校正結(jié)果。遮蔽因子的準確提取是決定地形校正結(jié)果的關(guān)鍵因素。本文使用同一個輻射傳輸模型，針對同一片地表裸露山區(qū)的SPOT- 5全色影像，分別使用對應(yīng)的DEM和直方圖閾值分割法計算遮蔽因子，并將得到的校正結(jié)果進行對比發(fā)現(xiàn):在DEM精度不夠高的情況下，使用直方圖閾值分割法計算遮蔽因子，可以得到更好的地形校正效果，充分證明了直方圖閾值分割法計算遮蔽因子的優(yōu)勢及使用價值。

輻射傳輸；遮蔽因子；地形校正

在定量遙感中，影像中的地物光譜信息受大氣、地形效應(yīng)等多種因素影響[1]，特別是復雜的山區(qū)影像，由于地形的起伏，各個像元接收到的有效光照差別較大[2]，并且影像空間分辨率越高，這種效應(yīng)就越明顯，為了消除這種地形效應(yīng)的影響，需要進行山區(qū)地形輻射校正。對于山區(qū)遙感影像，地表接收到的輻射主要有太陽直射輻射、天空的散射輻射、鄰近地表的反射和散射4個部分。由于受到坡度、坡向等因素影響，山地陰影區(qū)域不能接收到太陽直接輻射，其接收到的輻射主要為天空的散射及鄰近地表的反射與散射，因此在對高空間分辨率數(shù)據(jù)進行地形校正時，以上兩項必須加以考慮[3]。

地形校正算法主要分為3種類型：經(jīng)驗?zāi)Ｐ汀⑽锢砟Ｐ?、半?jīng)驗?zāi)Ｐ?。其中全色影像的地形校正主要是運用經(jīng)驗?zāi)Ｐ托Ｕā＝?jīng)驗?zāi)Ｐ蜑樘柸肷浣桥c衛(wèi)星傳感器接收的輻射之間的經(jīng)驗關(guān)系的總結(jié)。其優(yōu)點是所用到的地形因子少，簡單易懂且適用性很強，缺點是理論性不完善。波段比模型、Minnaer校正模型都屬于該模型。這些模型用于全色影像的地形校正，往往出現(xiàn)過校正或欠校正的效果。物理模型為通過研究光與地表的物理相互作用，基于輻射傳輸理論建立的模型。其優(yōu)點是理論基礎(chǔ)完善，模型參數(shù)物理意義明確，缺點是模型復雜，所需參數(shù)較多。這類模型有余弦校正法、SCS校正模型等。半經(jīng)驗?zāi)Ｐ蛣t是綜合了經(jīng)驗?zāi)Ｐ秃臀锢砟Ｐ偷膬?yōu)點[2,4]。相比之下，基于輻射傳輸理論的模型能克服理論性不完善和過校正的缺點，利用定量化信息對影像進行校正[5]。與此同時，基于輻射傳輸模型可以在更好地恢復地物像元之間相對關(guān)系的基礎(chǔ)上對全色遙感影像進行山區(qū)地形校正。

利用輻射傳輸模型校正時，遮蔽因子的精確計算對校正的效果至關(guān)重要。通常采用DEM和太陽有效入射角、坡度和坡向角度等地形參量結(jié)合的方法計算遮蔽因子，該方法受DEM精度影響較大，且陰影的邊緣位置判斷不準確，從而會導致地形校正不完整的現(xiàn)象。本文首先用影像頭文件的定標系數(shù)對全色影像進行定標，然后運用直方圖的特殊性計算遮蔽因子并進行地形校正，最后與DEM計算了遮蔽因子得到的地形校正結(jié)果作對比，分別從目視效果和統(tǒng)計指標等方面進行對比分析。

1 原理與方法

地形輻射校正的實質(zhì)是通過各種方法變換，將像元的輻亮度變換到某一相同基準面[6]。對復雜山區(qū)地形校正采用的輻射傳輸模型需要考慮天空散射輻射和鄰近地表反射輻射兩個重要影響因素。將地表視為朗伯體，且假定在所有高程下太陽輻射及大氣為均勻的情況下，地表反射率公式如下[7]

(1)

cosβ=cosθzcosS+sinSsinθzcos(φS-A)

(2)

式中，S為坡度；A為坡向,可通過DEM計算得到；φS為太陽方位角。以上為本文用來進行地形校正所使用的輻射傳輸模型。本文中b因子的提取采取兩種方法并對其進行分析。

DEM計算遮蔽因子主要是在ENVI的Topographic Modeling- Shaded Relief模塊中輸入太陽高度角、太陽方位角，通過對仰角的判斷遮蔽情況判別陰影區(qū)，即可通過波段運算得到值為0- 1的遮蔽因子。

直方圖計算遮蔽因子主要采用閾值分割法。閾值法是一類被廣泛使用的分割方法，它具有操作簡便、易于理解并能實現(xiàn)快速分割的優(yōu)勢，其中直方圖閾值分割是最為常用的方法?；叶戎狈綀D橫坐標代表圖像的像素灰度值，縱坐標表示具有該灰度級所包含的像素數(shù)或出現(xiàn)這個灰度級的概率，直方圖曲線可以體現(xiàn)出圖像對比度、明亮程度及灰度分布等方面的差異[9]。當圖像中目標物與背景出現(xiàn)明顯的谷峰形狀時，利用直方圖分割可以取得較好的效果[10]，本文基于該特征通過直方圖提取遮蔽因子。

總體來說，圖像直方圖可以分為單峰直方圖、雙峰直方圖及多峰直方圖3種類型。單峰直方圖比較常見，通常是由于目標區(qū)域與背景區(qū)域面積相差太大，從而像素數(shù)量相差懸殊造成的。雙峰直方圖是選取分割閾值的最理想情況。此類圖像明顯可以區(qū)分兩種地物。目標對象與背景區(qū)域各自對應(yīng)一個波峰，兩峰之間存在一個波谷，波谷區(qū)域的像素量是最少的，同時也代表著兩類區(qū)域的邊界處，可認為谷值為最佳的分割點。多峰直方圖則表示該類圖像包含的信息非常巨大且結(jié)構(gòu)復雜[10]。在進行山區(qū)遮蔽因子的提取時，對陰影區(qū)域的識別尤為重要。直方圖閾值法可以運用到山區(qū)地形校正中的陰影區(qū)域判別與提取。對于裸露的復雜山區(qū)遙感影像，地物信息較為單一，山體(陽坡、陰坡)、陰影邊界及陰影是影像的重要組成部分，三者的影像灰度值逐級遞減，而且在直方圖上表現(xiàn)為明顯的雙峰直方圖特征。兩峰之間的波谷對應(yīng)影像中類別為陰影邊界的像元，此類像元在頻數(shù)上明顯少于其余兩種地物類別。此時選取雙峰之間的谷值作為閾值，通過波段運算將檢測出的陰影區(qū)域設(shè)為0，非陰影區(qū)值設(shè)為1，便可得到更為精確的遮蔽因子。

2 試驗與分析

2.1 試驗過程

試驗所用數(shù)據(jù)為SPOT- 5 HRG2傳感器的全色影像數(shù)據(jù)，獲取日期為2013年6月13日。研究區(qū)域位于甘肅省張掖市(39°38′06″N，100°35′53″E)，海拔1630 m，地形以戈壁及山區(qū)為主。SPOT- 5的全色數(shù)據(jù)頭文件中獲取定標系數(shù)C0為1.320 660(1/W·m2·sr·μm)。試驗當天在地面同步測量了大氣數(shù)據(jù)。另外，試驗所用DEM數(shù)據(jù)(SRTM DEM，分辨率30 m)下載自地理空間數(shù)據(jù)云(http:∥www.gscloud.cn/)，如圖1所示。

圖1 對應(yīng)區(qū)域DEM數(shù)據(jù)

根據(jù)上述輻射傳輸模型，本文選用美國RSI(Research System Inc)的IDL(interactive data language)交互式數(shù)據(jù)語言作為開發(fā)平臺，同時結(jié)合6SV模型計算的輻射分量進行試驗。表1為利用6SV輻射傳輸模型計算時采用的輸入?yún)?shù)。

表1 SPOT- 5傳感器6S輻射傳輸模型輸入?yún)?shù)

在計算遮蔽因子時，分別采用DEM和直方圖兩種方法，并對比其獲取遮蔽因子的差異性。利用ENVI的Topographic- Topographic Modeling選取DEM數(shù)據(jù)，選取shaded relief生成陰影地貌圖像。進行波段運算，得到的二值化遮蔽因子如圖2所示。

隨后通過直方圖的方法計算遮蔽因子，打開全色影像，統(tǒng)計出直方圖信息，如圖3所示。

圖2 DEM計算的遮蔽因子

圖3 原始影像的灰度直方圖

從曲線的形狀可以看出該直方圖為雙峰直方圖且有明顯的波谷。對波谷進行量取可知影像灰度值小于或大于50的像元均符合正態(tài)分布，故選取雙峰之間的谷值50作為閾值，通過波段運算計算出遮蔽因子，結(jié)果如圖4所示。

圖4 直方圖統(tǒng)計圖計算遮蔽因子

2.2 結(jié)果分析

分別采用兩種方法計算遮蔽因子并利用式(1)進行地形校正，得到結(jié)果如圖5、圖6所示。

圖5 DEM計算遮蔽因子地形校正結(jié)果

圖6 直方圖計算遮蔽因子地形校正結(jié)果

對上述結(jié)果分別從目視判讀和指標統(tǒng)計兩個方面進行分析評價。

2.2.1 目視判讀評價

從校正結(jié)果可以看出，使用不同方法進行地形校正后，原始影像的陰影均得到了消除，且影像的凹凸感也得到了平滑[11]。從圖5中原始影像和校正結(jié)果影像的對比可以看出，大部分陰影區(qū)域并沒有得到很好的恢復，且有部分陰影區(qū)出現(xiàn)過校正的現(xiàn)象。從圖6中校正前后的對比可以看出，陰影區(qū)域均得到了很好的恢復，且陰影區(qū)域的山坡、山脊線輪廓清晰。

2.2.2 統(tǒng)計指標評價

本文采用離散分析、回歸分析和直方圖3個指標對校正結(jié)果進行指標統(tǒng)計評價。

2.2.2.1 離散分析

為了檢驗影像的整體效果，采用離散度(方差/均值)指標對研究區(qū)影像進行統(tǒng)計分析。離散度由下式計算[3]得出

式中，DI表示離散度；M與SD分別表示SPOT- 5影像校正前后的輻亮度均值及標準差。

表2 地形校正前后影像離散度統(tǒng)計值比較 (%)

地形校正前，影像離散度相對較大，通過不同方法校正后，離散度均得到了不同程度的減小，說明地形校正減弱了地形對地物信息遮擋的效應(yīng)，從整體上改善了圖像質(zhì)量[3,7]。直方圖法得到的遮蔽因子進行地形校正的結(jié)果離散度小于DEM法，這充分說明當DEM數(shù)據(jù)精度不高時，直方圖在提取遮蔽因子時的準確程度，可以大大提高地形校正的效果。

2.2.2.2 回歸分析

影像與太陽有效入射角之間的回歸關(guān)系反映了地形效應(yīng)的程度[12]。一般情況下，地形校正前后擬合的線性回歸方程斜率k越大，地形效應(yīng)越明顯[13- 14]。圖7、圖8分別為原始影像和不同校正方法所得結(jié)果與太陽入射角的回歸曲線示意圖。

圖7 原始影像與太陽入射角回歸曲線

圖8 地形校正結(jié)果與太陽入射角回歸曲線

從圖7、圖8可以看出，經(jīng)過校正后的回歸曲線斜率小于原始影像與太陽入射角的回歸曲線斜率，其中直方圖法校正結(jié)果(圖8(b))的斜率最小，說明經(jīng)過校正后，地形效應(yīng)得到明顯減小，且直方圖法的校正效果優(yōu)于DEM法。

2.2.2.3 統(tǒng)計直方圖

圖9 校正結(jié)果統(tǒng)計直方圖

圖9是兩種方法校正結(jié)果影像的灰度直方圖。對比圖3和圖9(a)可以看出，通過DEM計算出的遮蔽因子對地形校正有一定的效果，但是曲線的波谷說明該方法沒有很好地去除陰影與非陰影的邊界過渡問題。從圖3和圖9(b)的對比可以看出，直方圖呈正態(tài)分布，說明校正后影像的灰度值過渡平滑自然，陰影區(qū)域得到了很好地消除。與DEM計算遮蔽因子不同，直方圖計算遮蔽因子的方法是通過判斷在全色影像中，不同灰度值出現(xiàn)的頻率來分離待校正的陰影區(qū)域的值。之所以可以使用直方圖分割的方法，是因為在裸露山體的山區(qū)地物信息單一，并且在高分辨率的全色影像上能明顯區(qū)分陰影區(qū)域與非陰影區(qū)域。將陰影區(qū)域看作目標地物，非陰影區(qū)域作為背景信息。試驗中所選取的影像全部為山區(qū)且沒有植被等其他地物信息覆蓋。陰影區(qū)域與非陰影區(qū)域的比例均衡，不存在大面積都為陰影區(qū)域或只有極小一部分陰影區(qū)域的現(xiàn)象。當影像中出現(xiàn)部分城市或部分區(qū)域為植被所覆蓋，直方圖曲線都不出現(xiàn)雙峰，也不呈正態(tài)分布形狀；若影像中絕大多數(shù)面積都為陰影區(qū)域或非陰影區(qū)域，那么直方圖曲線只呈正態(tài)分布，同樣不出現(xiàn)雙峰。這樣均無法得到分割值，則無法通過該方法計算遮蔽因子。

使用直方圖方法計算遮蔽因子進行地形校正，一方面由于計算遮蔽因子的過程中不受DEM的影響；另一方面由于山區(qū)這種特殊的地形，能明顯區(qū)分陰影區(qū)域與非陰影區(qū)域。特別是在地表裸露地物信息單一的情況下，統(tǒng)計直方圖為雙峰狀，此時可以快速準確地找到分割值，進而計算遮蔽因子，得到較好的校正效果。

3 結(jié) 語

DEM計算遮蔽因子普遍適用于輻射傳輸方法的地形校正，但這種方法存在以下問題：

(1) DEM精度將影響地形陰影判斷的精度，如果能獲取高精度DEM數(shù)據(jù)，該方法可以很好地提取陰影區(qū)。由于大多數(shù)情況無法獲取高精度DEM，在選擇使用30 m DEM時，DEM精度對陰影判斷的精度有很大的影響，導致無法準確地提取陰影區(qū)域，造成校正結(jié)果不佳。

(2) 在判斷陰影區(qū)域的過程中一般均使用高程采樣點內(nèi)插方法。由于山區(qū)特殊的復雜地形，在高程插值時容易造成高程值計算不準確，并且無法準確確定地形遮蔽角的搜索半徑。這些都會造成無法相對準確地提取遮蔽因子，從而導致地形校正效果不理想。雖然延長地形遮蔽角的搜索半徑可以使對陰影區(qū)域的判斷更加精確，但是會增大計算量，降低校正時的計算效率[13,15]。

本文采用直方圖計算遮蔽因子的方法進行地形校正，取得了較好的效果。通過與DEM法校正結(jié)果的比較分析可以得出，對于特殊的山區(qū)地形影像，使用直方圖方法進行地形校正的結(jié)果要優(yōu)于DEM法。

[1] 趙英時. 遙感應(yīng)用分析原理與方法[M]. 北京：科學出版社, 2003.

[2] 段四波, 閻廣建. 山區(qū)遙感圖像地形校正模型研究綜述[J]. 北京師范大學學報(自然科學版)，2007, 43(3): 362- 366.

[3] 高永年. 遙感影像地形校正理論基礎(chǔ)與方法應(yīng)用[M]. 北京：科學出版社, 2013.

[4] 張兆明, 何國金, 劉定生, 等. 一種改進的遙感影像地形校正物理模型[J]. 光譜學與光譜分析, 2010, 30(7): 1839- 1842.

[5] 孫源, 顧行發(fā), 余濤, 等. 山地森林地區(qū)遙感影像地形輻射校正研究[J]. 遙感信息,2009(2): 7- 12.

[6] 陳志明, 李家國, 余濤, 等. TM遙感影像的地形輻射校正研究 [J]. 遙感信息, 2009(2): 29- 33.

[7] RICHTER R. Correction of Satellite Imagery over Mountainous Terrain [J]. Applied Optics, 1998, 37(18): 4004- 4015.

[8] 梁順林. 定量遙感 [M]. 北京：科學出版社, 2009.

[9] 侯越. 基于灰度直方圖的閾值分割算法 [J]. 硅谷, 2010(23): 165- 168.

[10] 盧建華. 基于直方圖閾值法的遙感圖像分割算法研究[D].福州：福建農(nóng)林大學, 2013.

[11] 李翠翠, 樊基倉, 付瀟華, 等. 復雜地形山區(qū)Landsat TM影像C校正策略與實驗 [J]. 地球信息科學學報, 2014(1): 134- 141.

[12] RICHTER R, MüLLER A. De- shadowing of Satellite/Airborne Imagery [J]. International Journal of Remote Sensing, 2005, 26(15): 3137- 3148.

[13] 高永年, 張萬昌. 遙感影像地形校正研究進展及其比較實驗 [J]. 地理研究, 2008, 27(2): 467- 477.

[14] LI A, WANG Q, BIAN J, et al. An Improved Physics- based Model for Topographic Correction of Landsat TM Images [J]. Remote Sensing, 2015, 7(5): 6296- 6319.

[15] 羅慶洲, 劉順喜, 曾齊紅, 等. 基于DEM的本影與落影判斷研究 [J]. 國土資源遙感, 2009(2):29- 31.

Impact and Analysis of Shadow Area Extracting to the Terrain RadiationCorrection for High Resolution Image

SONG Liyao1,ZHANG Hao2,LIN Feng3

(1. China University of Petroleum(Huadong),Qingdao 266555,China; 2. Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences,Beijing 100094,China; 3. Qingdao Institute of Marine Geology, Qingdao 266555,China)

By using transporting mechanism to correct topographic radiation, we can fully consider direct solar radiance, sky diffuse radiance and near surface reflectance radiance, and get a better result. The key factor that determines topographic correction result is accurate extraction of shadow factor. The article use DEM method and histogram threshold segmentation method to compute shadow factor of the same area of SPOT- 5 panchromatic imagery, and then use the same radiation transmission model to correct topographic radiation. By contrast, the result that using the histogram threshold method is better than DEM method when DEM accuracy is not high enough, and these have fully proved its value and advantage.

radiative transfer； shadow factor； topographic correction

2016- 06- 26；

2017- 01- 23

國家自然科學基金(41325004)；國家重點研發(fā)計劃(2016YFB0501501) 作者簡介： 宋麗瑤(1990—)，女，碩士生，研究方向為地形輻射校正。E- mail：631212606@qq.com 通信作者： 張 浩。E- mail：haozhang@ceode.ac.cn

宋麗瑤，張浩，林峰.陰影區(qū)提取對高分辨率影像地形輻射校正的影響與分析[J].測繪通報，2017(3)：29- 33.

10.13474/j.cnki.11- 2246.2017.0078.

P237

A

0494- 0911(2017)03- 0029- 05