劉華堯，李舟，潘真，文思

(1.桂林理工大學(xué)測繪地理信息學(xué)院，廣西桂林541004； 2.廣州南方測繪儀器有限公司，廣東廣州 510665；3.桂林理工大學(xué)土木與建筑工程學(xué)院，廣西桂林 541004)

基于TM影像的FLAASH模塊大氣糾正及評價(jià)

劉華堯1?，李舟2，潘真3，文思1

(1.桂林理工大學(xué)測繪地理信息學(xué)院，廣西桂林541004； 2.廣州南方測繪儀器有限公司，廣東廣州 510665；3.桂林理工大學(xué)土木與建筑工程學(xué)院，廣西桂林 541004)

在高光譜遙感數(shù)據(jù)處理中，由于受到大氣的影響，傳感器接收到的輻射信息不能真實(shí)地反映地表反射光譜信息，因此對遙感影像進(jìn)行大氣糾正去除其影響，是高光譜遙感數(shù)據(jù)處理中極為重要的環(huán)節(jié)。文章介紹了FLAASH模塊的基本原理和TM數(shù)據(jù)在進(jìn)行大氣糾正時(shí)的參數(shù)設(shè)置，并用FLAASH模塊對豫西北區(qū)TM遙感影像進(jìn)行大氣校正以及對處理結(jié)果進(jìn)行了薄云糾正程度，光譜糾正情況和歸一化植被指數(shù)等的評價(jià)，并與ATCOR模型進(jìn)行了對比。證明FLAASH模塊大氣糾正效果良好。

遙感影像；大氣糾正；FLAASH；TM數(shù)據(jù)

1 引 言

太陽輻射在真空中傳播時(shí)幾乎不受大氣的影響。然而當(dāng)太陽輻射到達(dá)地面在進(jìn)人傳感器以前，首先經(jīng)過大氣層，在這個(gè)過程中會受到大氣中分子、氣溶膠和云粒子等大氣吸收與散射，這對于定量研究和分類精度造成很大的影響[1]。目前國內(nèi)外使用的大氣糾正模型主要有CRORN、ACTOR、MODTRAN和FLAASH等[2]。FLAASH大氣糾正模型在國內(nèi)應(yīng)用較少，F(xiàn)LAASH是目前精度較高的大氣輻射糾正模型，能有效消除大氣和光照等因素對地物反射的影響，獲得地物較為準(zhǔn)確的信息。本文對豫西北地區(qū)TM高光譜數(shù)據(jù)通過FLAASH模塊進(jìn)行大氣糾正。

2 FLAASH模塊大氣糾正原理及算法

FLAASH是基于MODTRAN4模型的大氣糾正模塊，主要用于從高光譜遙感影像中還原出無大氣影響或受大氣影響較小的地物地表反射率，它糾正的波長范圍為0.4 um～3 um。與其他大氣糾正模型的計(jì)算方法不同，F(xiàn)LAASH直接與MODTRAN4的大氣輻射傳輸編碼相結(jié)合，而不采用模型數(shù)據(jù)庫中加入的輻射參數(shù)對大氣進(jìn)行糾正，因此任何標(biāo)準(zhǔn)MODTRAN4大氣模型和氣溶膠類型都可以被直接選用，并計(jì)算處無大氣影響的地表反射率。FLAASH能夠生成薄云和卷云的分類影像，對光譜進(jìn)行平滑，消除噪聲FLAASH大氣糾正基于太陽波譜范圍內(nèi)(不包括熱輻射)標(biāo)準(zhǔn)的平面朗伯體在傳感器處接收到的單個(gè)像元光譜輻射亮度按式(1)進(jìn)行計(jì)算[3]:

式中L為傳感器接收到的某個(gè)像元的輻射強(qiáng)度；ρ表示該像元的地表反射率；ρe為該像元及周邊像元的混合平均地表反射率；S為大氣的球面反射率；A，B是由大氣條件及地表下墊面幾何條件所決定的系數(shù)。La為太陽輻射經(jīng)大氣散射后再由地表向上反射通過大氣進(jìn)入接收傳感器單元的一部分輻照度。大氣散射會引起“鄰近像元效應(yīng)”，而大多數(shù)的大氣輻射糾正模塊中一般假設(shè)ρ=ρe，這樣的糾正方式忽略了“鄰近像元效應(yīng)”，這些模型在地物類型單一而且具有較高能見度的情況下是可行的，但在有薄云或地表反射對比強(qiáng)烈的條件下會導(dǎo)致在短波范圍內(nèi)大氣糾正結(jié)果誤差較大。FLAASH是利用大氣點(diǎn)擴(kuò)散函數(shù)進(jìn)行空間均衡化處理對鄰近像元效應(yīng)進(jìn)行了糾正。在La、S、A、B、ρe已知的前提下，可以使用式(2)計(jì)算出像元的空間平均輻射率[4]。e

可以通過原始影像的基本參數(shù)計(jì)算得出。從而反演出影像的真實(shí)輻射率。在計(jì)算鄰近像元反射率時(shí)，傳感器接收到直接來自目標(biāo)地物的光子的信號被分為經(jīng)大氣散射后進(jìn)入傳感器和經(jīng)鄰近像元散射后進(jìn)入傳感器兩部分即:

式中ρRA為大氣散射后直接進(jìn)入傳感器的輻射率，ρ為地物真實(shí)反射率，TRA(us)為傳感器接收到的大氣向上散射透過率，eτ/uv為大氣向上散射透過率，td(uv)為大氣向上漫射透過率，〈ρe〉為鄰近像元反射率，由大氣點(diǎn)擴(kuò)散函數(shù)與像元空間平均反射率的乘積在徑向距離上的二重積分運(yùn)算獲取，表示為:

其中x，y表示鄰近像元到中心像元的幾何距離，在FLAASH中大氣點(diǎn)擴(kuò)散函數(shù)fr(x，y)用一個(gè)徑向距離的近似指數(shù)函數(shù)代替。ρe(x，y)為像元的空間平均反射率。鄰近像元反射率獲取之后，即可推得式(5)從而得到地物真實(shí)反射率ρ。

3 TM數(shù)據(jù)影像處理和參數(shù)設(shè)置

3.1 TM影像數(shù)據(jù)預(yù)處理

本文使用的TM影像來自于1984年3月發(fā)射的Landsat5衛(wèi)星，它的衛(wèi)星高度為7.5 km，半主軸為7285.438 km，傾角98.2°，覆蓋周期為16天掃描寬度185 km，有7個(gè)波段。研究區(qū)位于河南省西北部，地理坐標(biāo)為東經(jīng)112°33′40″～113°38′42″，北緯34°48′55″～35°29′59″。東西長120 km，南北寬77 km。地勢北高南低，地形變化較大地貌類型多樣，總體上可以分為山地、丘陵、平原三部分。平原地區(qū)的平均海拔為80 m，丘陵海拔為200 m，山區(qū)平均海拔為1 400 m[5]。區(qū)域內(nèi)植被種類分布多樣，使用TM影像時(shí)各種地物區(qū)分度高，符合實(shí)驗(yàn)條件。

首先對影像數(shù)據(jù)進(jìn)行波段疊加，將TM影像的多個(gè)波段數(shù)據(jù)集中為一個(gè)包含多光譜的文件。但是由于波段疊加，影像文件所包含的波段信息被清空，為了便于計(jì)算機(jī)分析，要將相關(guān)的信息(衛(wèi)星類型、波段)進(jìn)行重新輸入，依據(jù)TM衛(wèi)星各波段波長分布重新輸入波長，為確保研究普遍適用性本文使用的TM影像各波段范圍的中間值。

完成波段參數(shù)修改后，影像顯示的是傳感器的感應(yīng)值，而非輻射的絕對值，需要對影像進(jìn)行反射輻射值的計(jì)算，本文根據(jù)影像數(shù)據(jù)的采集時(shí)間將參數(shù)設(shè)為2004年8月30日，太陽高度角:55.233 228進(jìn)而進(jìn)行輻射糾正。由于模塊中要求輻射能量的量綱是μW/(cm2·nm·sr)，而經(jīng)輻射定標(biāo)TM數(shù)據(jù)的量綱為W/(cm2·um·sr)相差10倍的關(guān)系，所以我們得到的輻射圖像還要進(jìn)行尺度因子轉(zhuǎn)換。最后將得到的輻射文件轉(zhuǎn)換成FLAASH要求的輸入BIL或者BIP格式。

3.2 糾正參數(shù)設(shè)置

FLAASH模塊大氣糾正首先要對影像中的大氣參數(shù)，包括氣溶膠光學(xué)厚度、氣溶膠類型等參數(shù)進(jìn)行設(shè)置。因?yàn)檠芯繀^(qū)域大致為于北緯34°地區(qū)，時(shí)間為8月，所以選擇MLS模型；因?yàn)閭鞲衅鳑]有適當(dāng)?shù)牟ǘ蝸硌a(bǔ)償水汽的影響，所以水汽去除參數(shù)保持默認(rèn)；氣溶膠模型選擇Urban；氣溶膠獲取選擇2-Band(K-T)；初始可見度值設(shè)置為30 km；光譜平滑選擇Yes；使用鄰域糾正；平滑操作的幅寬為9；Aerosol scale height定義為2 km；CO2mixing ratio定義390；Reuse MODTRAN Calculations選擇NO；Modtran Resolution選擇5 cm；Modtran Multiscatter Model選擇Scaled DISORT這個(gè)模型對波長小于1 000 nm的波段處理經(jīng)度較高，其Streams值設(shè)為8；Output scale factor定義為10 000。由于本文影像中零背景比較單一所以tile size值設(shè)置為100 M。

4 大氣糾正應(yīng)用及結(jié)果分析

本文選用采集時(shí)間為2004年8月30日的Landsat衛(wèi)星TM影像為例進(jìn)行FLAASH大氣糾正試驗(yàn)，糾正結(jié)果如圖1、2所示。由圖可見糾正前后的圖像視覺清晰度有了明顯提高，經(jīng)FLAASH大氣糾正后的TM影像對比度增強(qiáng)，圖像較糾正前更加清晰。根據(jù)影像不同位置的比較，發(fā)現(xiàn)經(jīng)過糾正的山區(qū)和城鎮(zhèn)區(qū)域清晰度都有所提高。說明FLAASH有效去除了大氣中氣溶膠、水汽等大氣因素的影響。

圖1 糾正前后山區(qū)影像對比

圖2 糾正前后山區(qū)影像對比

為進(jìn)一步對Landsat衛(wèi)星TM影像的FLAASH大氣糾正結(jié)果進(jìn)行驗(yàn)證與評價(jià)，利用典型地物波譜曲線和歸一化植被指數(shù)(NDVI)對原始圖像與大氣糾正后的圖像進(jìn)行對比分析。

(1)典型地物波譜曲線分析，在原始圖像和大氣糾正后的圖像中分別隨機(jī)選取同名點(diǎn)地物水體、植被和建筑3種典型地物進(jìn)行光譜曲線分析(圖3)。由圖3可見，隨機(jī)選取的3種典型地物中在可見光波段原始圖像的反射率總體比經(jīng)大氣糾正后圖像的反射率要高，這是因?yàn)樵紙D像的反射率不是地物的真實(shí)反射率，它包含大氣中煙、塵埃、小水滴和氣溶膠等引起的瑞利散射、米氏散射和多次散射等作用，從而導(dǎo)致地物在可見光波段的反射率總體偏高，所以經(jīng)大氣糾正后的各典型地物在該波段的反射率均比原始圖像反射率降低；而經(jīng)大氣糾正后植被反射率在近紅外波段上升明顯，這是因?yàn)樵诮t外波段對其的散射較小。因此從地物光譜分析的角度來看，證明FLAASH大氣糾正可以有效地去除Landsat衛(wèi)星TM數(shù)據(jù)中大氣因素的影響。

圖3 糾正前后地物波譜曲線對比

(2)歸一化植被指數(shù)NDVI分析，NDVI是一種由遙感傳感器所接收的地物光譜信息推算而得的反映地表植被狀況的定量值。所以就可以利用歸一化植被指數(shù)分別分析原始圖像與大氣糾正后的圖像。NDVI的計(jì)算公式為[6]:

NDVI=(NIR-R)/(NIR+R)(6)

式中R和NIR分別為紅光和近紅外波段，分別是Landsat衛(wèi)星TM數(shù)據(jù)的第3和第4波段。在糾正前后的兩幅影像中隨機(jī)選取水體、植被、建筑等同名點(diǎn)的歸一化植被指數(shù)(NDVI)值進(jìn)行對比分析如表1所示。

糾正前后歸一化植被指數(shù)比較 表1

由表1可以得到3種地物經(jīng)FLAASH大氣糾正后的NDVI值總體上比原始影像的NDVI值高。特別是植被地區(qū)，其糾正后的歸一化植被指數(shù)提高幅度較大，說明FLAASH大氣糾正模塊對于植被覆蓋較密的地區(qū)的TM影像的糾正效果更為明顯。

(3)與ATCOR糾正結(jié)果的對比分析，為了對校正結(jié)果進(jìn)行評價(jià)以植被的地物反射曲線為例，從圖像中分別按FLAASH和ATCOR模型糾正后提取植被的反射率曲線進(jìn)行對比。

圖4 FLAASH和ATCOR校正后植被的地物反射率曲線對比

圖4分別是采用FLAASH和ATCOR校正后的植被地物反射率曲線。由圖4分析FLAASH和ATCOR糾正都能夠較好地糾正植被的反射率光譜特征，可以去除絕大多數(shù)的大氣影響，但是由于FLAASH綜合考慮了大氣中的分子散射和氣溶膠散射的影響，在綠光波段FLAASH糾正后的反射率比ATCOR校正結(jié)果更接近真實(shí)反射率；紅外波段FLAASH較ATCOR也更為精確。將兩種方法的校正結(jié)果與真實(shí)結(jié)果進(jìn)行誤差分析，在整個(gè)波譜范圍內(nèi)FLAASH與真實(shí)反射率曲線的絕對誤差平均值在0.007 3～0.015之間，而ATCOR與真實(shí)反射率曲線的絕對誤差均值在0.003 8～0.063之間，這在一定程度上表明FLAASH糾正的結(jié)果更為穩(wěn)定。

5 結(jié) 論

本文通過對豫西北地區(qū)的高光譜數(shù)據(jù)進(jìn)行FLAASH大氣糾正及評價(jià)，結(jié)果表明，F(xiàn)LAASH能夠有效去除TM影像的大氣影響和糾正漫反射引起的領(lǐng)域效應(yīng)，經(jīng)糾正后的影像地表地物的影像清晰度有明顯提升。FLAASH能夠有效改善典型地物的光譜特征以及NDVI值，經(jīng)糾正后的典型地物光譜特征和NDVI值都有明顯的提升。經(jīng)與ATCOR模型進(jìn)行對比，F(xiàn)LAASH也體現(xiàn)出了一定的優(yōu)勢。由此可見，應(yīng)用FLAASH模塊能較好地消除TM影像的大氣影響，從而更好地從高光譜遙感影像中獲取地物的地表反射率信息。

[1] 宋曉宇，王紀(jì)華，劉良云等.基于高光譜遙感影像的大氣糾正:用AVIRIS數(shù)據(jù)評價(jià)大氣糾正模塊FLAASH[J].遙感技術(shù)與應(yīng)用，2005，20(4):393～398.

[2] Matthew MW，Adler-Golden SM，Berk A，etal.Atmospheric correction of spectral imagery:Evaluation of the FLAASH algorithm with AVIRIS data[J].Applied Imagery Patten Recognition Workshop，2002，31st:157～163.

[3] 何海艦.基于輻射傳輸模型的遙感圖像大氣校正方法研究[D].長春:東北師范大學(xué)，2006.

[4] Fallah-Adl H，JaJ，Liang S.Efficientalgorithms for estimating atmospheric parameters for surface reflectance retrieval[C]. 1996 International Conference on Parallel Processing，1996: 132～140.

[5] 中科學(xué)院數(shù)據(jù)庫國際科學(xué)數(shù)據(jù)服務(wù)中心LANDSAT5TM衛(wèi)星數(shù)字產(chǎn)品[EB/OL].2012.

[6] 李苗苗，吳炳方，顏長珍.密云水庫上游植被覆蓋度的遙感估算[J].資源科學(xué)，2004，26(4):135～158.

FLAASH M odule Atmospheric Correction and Evaluation Based on TM Image

Liu Huayao1，Li Zhou2，Pan Zhen3，Wen Si1
(1.College of Geomatics and Geoinformation，Guilin University of Technology，Guilin 541004，China；2.Guanzhou South Surveying＆Mapping Instrument Co.，Ltd.Guangzhou 510665，China；3.College of Geomatics and Geoinformation，College of Civil Engineening and Architecture，Guilin 541004，China)

In hyper spectral remote sensing data processing，due to the influence of the atmosphere，the sensor can not receive the truly reflection spectrum information reflect form the surface.Therefore for remote sensing image for atmospheric correction to remove its influence is very important link in hyper spectral remote sensing data processing.This paper introduces the basic principle of FLAASH module and TM data parameter settings for atmospheric correction.This paper use FLAASH module to correct the TM remote sensing image ofwestern north of Henan，It dealwith the correct degree of cloud，spectrum correction and Normalized Difference Vegetation Index.The results show that FLAASH module atmospheric correction effect is good.

remote sensing image；atmospheric correction；FLAASH；TM data

1672-8262(2013)05-63-04

TP75

A

2012—12—13

劉華堯(1986—)，男，碩士研究生，研究方向?yàn)镚IS基礎(chǔ)理論與應(yīng)用方法。

廣西自然科學(xué)基金重點(diǎn)項(xiàng)目(2011GXNSFD018003)；廣西空間信息重點(diǎn)實(shí)驗(yàn)室基金(1103108-03)。