孟 鵬 燕，徐 元 進(jìn)，陳 曦

(1.中國地質(zhì)大學(xué)數(shù)學(xué)地質(zhì)遙感地質(zhì)研究所，湖北 武漢 430074;2.湖北省地質(zhì)調(diào)查院，湖北 武漢430034)

?

地形對熱紅外數(shù)據(jù)反演林地表面溫度的影響

孟 鵬 燕1，徐 元 進(jìn)1，陳 曦2

(1.中國地質(zhì)大學(xué)數(shù)學(xué)地質(zhì)遙感地質(zhì)研究所，湖北 武漢 430074;2.湖北省地質(zhì)調(diào)查院，湖北 武漢430034)

基于同日過境的Landsat ETM+和TERRA ASTER 影像，評估地形對林地表面溫度的影響，并針對復(fù)雜崎嶇地形建立一個林地表面溫度歸一化模型，集中研究復(fù)雜地形(數(shù)字高程模型DEM和山體陰影HS)對地表溫度的影響，因此只選擇同種地物(林地)用于建模；此外，利用逐步掩膜的方法分析了林地邊緣混合像元帶(過渡帶)對模型反演精度的影響。研究結(jié)果表明，地表溫度和DEM呈負(fù)相關(guān)關(guān)系，與山體陰影呈正相關(guān)關(guān)系；兩種衛(wèi)星數(shù)據(jù)都證實了地形對地表溫度的影響，且林地邊緣過渡帶寬度為30 m時模型的精度達(dá)到最高。

地形；熱紅外遙感；地表溫度；林地；DEM

0 引言

熱紅外(TIR)遙感數(shù)據(jù)是模擬地表能量平衡和反演地表溫度的重要參數(shù)[1]，對于理解地表熱輻射過程不可或缺?，F(xiàn)有許多星載和機(jī)載傳感器可以獲取不同空間或輻射分辨率的熱紅外輻射，如ETM+和ASTER空間分辨率分別是60 m和90 m，輻射分辨率分別是8 bit(量化級為255)和12 bit(量化級為4 095)。利用熱紅外遙感數(shù)據(jù)反演的溫度通常被稱作地表溫度[2]，在本研究中地表溫度是指森林冠層表面溫度。

許多學(xué)者研究了地表溫度反演的相關(guān)影響因素。劉向培等[3]研究了植被覆蓋對中國東部地表溫度變化的影響，認(rèn)為植被覆蓋與地表溫度變化之間呈負(fù)相關(guān)關(guān)系；李彩霞等[4]評價了不同點(diǎn)位的太陽輻射對地表溫度的影響；彭菁菁等[5]研究了熱輻射方向性與地表溫度的關(guān)系；曹麗琴等[6]模擬了下墊面覆蓋類型對熱導(dǎo)效應(yīng)的影響；魏合理等[7]分析了風(fēng)場條件下的地表熱紅外輻射統(tǒng)計特征；Price[8]提出了影響地表溫度反演的各種因素：植被的蒸散(蒸發(fā))、大氣、地面風(fēng)的速度、表面粗糙度、熱儲存能力、反照率、發(fā)射率和地表含水量。

地形也是一個控制地表溫度的重要影響因素[9,10]，地形因素一般是指高程、緯度、坡向和坡度。高程和緯度影響一般的氣候特征，如年均氣溫、年降水量等，但這種差別與區(qū)域有關(guān)；坡向和坡度主要影響潛在的輻射和熱負(fù)荷[11]，山體陰影對能量平衡有重要的影響[12]，特別是在復(fù)雜的地形環(huán)境下，這種影響表現(xiàn)得更為明顯。

本研究主要探討地形對地表溫度建模的影響(由于研究區(qū)面積較小，故忽略緯度的影響)，獨(dú)特之處在于使用遙感數(shù)據(jù)對地形相對復(fù)雜崎嶇的林地表面溫度進(jìn)行建模，這是利用實地測量手段很難做到的，這種建模有助于理解地形和太陽的位置(天頂角和方位角)對森林冠層熱特性的影響。

1 研究區(qū)概況及數(shù)據(jù)源

研究區(qū)位于云南紅河自治州個舊市西部山地，地理范圍為102.93°～103.1°E，23.35°～23.49°N，面積為225 km2；該地地形復(fù)雜，最低海拔為856 m，最高海拔為2 518 m，平均海拔1 890.7 m。研究區(qū)林地覆蓋較為連續(xù)，植被覆蓋率大于35%，植被類型為亞熱帶常綠闊葉林。

本研究使用的數(shù)據(jù)是Landsat ETM+影像和搭載于Terra之上的ASTER影像，將兩傳感器的可見光/近紅外反射率波段用于林地分類，將熱紅外波段用于地表溫度模型反演；ETM+熱紅外傳感器譜段范圍為10.4～12.5 μm，ASTER所攜帶的熱紅外系統(tǒng)有5個波段，分別為：8.125～8.475 μm、8.475～8.825 μm、8.925～9.275 μm、10.25～10.95 μm、10.95～11.65 μm(表1)。

研究中使用兩種遙感數(shù)據(jù)，其衛(wèi)星在研究區(qū)的過境日期相同(2000年11月2日)，過境時間接近(北京時間11∶25、12∶02)，這就保證了太陽高度角、太陽方位角及大氣和光照條件一致，兩種衛(wèi)星都是近極太陽同步軌道，軌道高度相同(赤道處為705 km)且其重訪周期都為16 d。研究中使用的地形數(shù)據(jù)是從美國地質(zhì)調(diào)查局(USGS)對地觀測和科學(xué)中心獲取的分辨率為30 m的高精度DEM數(shù)據(jù)，并且將它和兩景衛(wèi)星影像轉(zhuǎn)換成相同的投影和坐標(biāo)系統(tǒng)(基準(zhǔn)面為WGS-84，投影為UTM)，以確保研究結(jié)果的精度。

表1 波段空間分辨率和波段間像元大小比值

Table 1 Spatial resolution and pixel size ratio between bands

衛(wèi)星傳感器空間分辨率(TIR)空間分辨率(V/NIR：用于林地分類)ratioETM+60m30m4ASTER90m15m36

注：ratio表示一個溫度像元在空間大小上對應(yīng)的林地像元數(shù)。

2 數(shù)據(jù)和研究方法

2.1 林地分類提取

使用ENVI4.8遙感處理軟件的監(jiān)督分類方法(光譜角制圖)對研究區(qū)的兩種影像(使用相同位置的感興趣區(qū)域)進(jìn)行分類，獲得研究區(qū)的兩個林地覆蓋專題圖。在分類過程中需要考慮林地邊緣的混合像元帶，即林地到非林地的過渡帶，過渡帶包含林地?zé)彷椛涮卣?，與林地區(qū)域產(chǎn)生熱通量偏差，影響林地表面溫度模型的精度，因此需要去除這部分影響。利用逐步掩膜方法去除林地邊緣過渡帶(對于ASTER影像，15 m為一個單位寬度；對于ETM+，30 m為一個單位寬度)，并逐步評估去除后的影響，選擇最佳去除寬度，使其對地表溫度模型的影響最小。圖1為由ASTER數(shù)據(jù)所得林地邊緣過渡帶的一個示例(過渡帶寬度為180 m)。

圖1 林地邊緣過渡帶示意

Fig.1 The schematic diagram of woodland edge buffer

2.2 地表溫度反演

利用ETM+熱紅外通道第6波段和ASTER熱紅外通道第13波段計算地表溫度。采用輻射傳輸方程法[13]，通過大氣輻射傳輸模型(6S、MODTRAN等)根據(jù)標(biāo)準(zhǔn)大氣輪廓線數(shù)據(jù)或?qū)崟r的大氣探空數(shù)據(jù)估計大氣對地表熱輻射的影響；然后把這部分大氣影響從衛(wèi)星傳感器觀測的熱輻射總量中去除，獲得地表熱輻射亮度；最后根據(jù)地表比輻射率轉(zhuǎn)化為地表溫度。

衛(wèi)星傳感器接收到的熱輻射亮度可通過DN值轉(zhuǎn)化而來，公式如下：

Lλ=gain·DN+offset

(1)

式中：DN 是像元的灰度值，gain 是該波段所對應(yīng)的增益量，offset 是該波段的偏移量。

衛(wèi)星傳感器接收到的熱輻射亮度Lλ包括：大氣上行輻射亮度 L↑，地表熱輻射亮度經(jīng)大氣吸收后到達(dá)衛(wèi)星傳感器的能量，大氣下行輻射亮度L↓。因此，衛(wèi)星傳感器接收到的熱輻射亮度可表示為：

Lλ=[ε L(λ,T)+(1-ε)L↓]τ+L↑

(2)

式中：ε表示地表比輻射率(常溫下，闊葉林比輻射率為0.98)，T表示地表溫度，L(λ,T)是溫度為T的黑體輻射亮度，τ為熱紅外波段的大氣透過率。

大氣上行輻射亮度、大氣下行輻射亮度和大氣透過率可利用NASA官網(wǎng)(http://atmcorr.gsfc.nasa.gov/)的大氣校正參數(shù)計算器計算得出。則溫度為T時的黑體輻射亮度L(λ,T)表達(dá)式為：

L(λ,T)=[Lλ-L↑-τ(1-ε)L↓]/(τ·ε)

(3)

由于黑體是郎伯體，其輻射亮度L(λ,T)與輻射功率W(λ,T)的關(guān)系為：

L(λ,T)=W(λ,T)/ π

(4)

根據(jù)普朗克黑體輻射規(guī)律，W(λ,T)表達(dá)式為：

W(λ,T)=2πhc2/[λ5(ehc /λkT-1)]

(5)

式中：W(λ,T)為單位面積黑體單位波長間隔向2π空間輻射的功率，h為普朗克常數(shù)(h=6.62606896×10-34J·s)，c為真空中光速(c=2.99792458×108m/s)，λ為熱紅外波段波長，k為波爾茲曼常數(shù)(k=1.3806505×10-23J/K)，T為地表溫度，T表達(dá)式為：

T=C1/[λ×ln(1+C2/(λ5·L(λ,T)))]

(6)

式中：C1、C2為常數(shù)，C1=hc/k=1.43877505×10-2m·k，C2=2hc2=1.19104276×10-16W·m2。

通過以上參數(shù)計算得出地表溫度(K)。由于無法獲取整個研究區(qū)的實測地表溫度，無法對研究區(qū)每個位置的地表溫度進(jìn)行驗證，但可從NASA官網(wǎng)獲取當(dāng)天的探空數(shù)據(jù)，并利用探空數(shù)據(jù)對ETM+和ASTER反演的地表溫度進(jìn)行驗證。驗證結(jié)果可得：LandsatETM+反演研究區(qū)地表溫度的精度為0.42℃，ASTER反演研究區(qū)地表溫度精度為0.31℃。由于熱紅外波段的空間分辨率是可見光/近紅外波段的空間分辨率的整數(shù)倍(表1)，而利用逐步掩膜方法去除林地邊緣要求地表溫度的像元尺寸和可見光/近紅外波段的尺寸一致，因此，需對地表溫度進(jìn)行重采樣，使其與分類結(jié)果像元一一對應(yīng)。

2.3 地形數(shù)據(jù)獲取

為了精確估算地形的影響，需要對高程數(shù)據(jù)進(jìn)行處理。首先將DEM數(shù)據(jù)重采樣，使其分別對應(yīng)ETM+和ASTER的空間分辨率(分別是30m和15m)；然后使用ArcGIS的空間分析模塊分別獲取對應(yīng)兩種影像空間分辨率的坡度和坡向數(shù)據(jù)；最后在DEM的基礎(chǔ)上計算出山體陰影[11]HS，公式為：

HS=255[cos(90-Z)sin(S)cos(α-A)+

sin(90-Z)cos(S)]

(7)

式中：Z表示太陽天頂角，S為坡度，α為坡向，A為太陽方位角。

山體陰影將假想的日照賦給柵格圖像中的每一個像元，因此，山體陰影對應(yīng)于不同地形區(qū)域的能量平衡差異。太陽天頂角和方位角表示太陽在一年和一天中特定的位置，由于衛(wèi)星過境具體時間不同，對應(yīng)的太陽位置也不同，因此，對應(yīng)的山體陰影值也略有差異。最后，分別分析兩種衛(wèi)星影像對應(yīng)的地表溫度、山體陰影和高程三者之間的關(guān)系。

2.4 定量轉(zhuǎn)換關(guān)系分析

使用線性回歸的方式計算了研究區(qū)內(nèi)地形對地表溫度的影響，回歸方程中主要有高程(DEM)和山體陰影(HS)兩個獨(dú)立變量，地表溫度(T)與二者之間的線性回歸方程為：

T=B0+B1×HS+B2×DEM+err

(8)

T由衛(wèi)星熱紅外波段記錄輻射值計算所得；HS和DEM在上文中已給出；誤差項(err)主要是指系統(tǒng)誤差和隨機(jī)誤差，系統(tǒng)誤差主要依賴地表溫度像元的齊性特征，一般齊性像元對于線性方程有最好的擬合優(yōu)度，因此考慮去除林地邊緣的混合像元以減少方程中系統(tǒng)誤差。

要建立地表溫度的歸一化模型，首先計算回歸方程，然后將計算線性方程所得回歸系數(shù)值(B0,B1,B2)代入式(8)中，得到地表溫度的歸一化模型。

3 模型回歸結(jié)果

地表溫度歸一化模型的擬合優(yōu)度(R2)很大程度上取決于回歸方程計算過程中地表溫度(T)的選擇。由于林地邊緣過渡帶影響邊緣地表溫度像元的齊性特征，因此，應(yīng)考慮林地邊緣過渡帶的寬度對地表溫度所產(chǎn)生的影響。通過線性回歸，計算出不同去除寬度對應(yīng)的R2值，圖2展示了林地邊緣過渡帶的去除寬度和R2值之間的依賴關(guān)系。從圖2中可以看出，起初，在ETM+和ASTER兩個傳感器對應(yīng)的方程中，R2值隨著林地邊緣過渡帶去除寬度的增加而增加；當(dāng)去除寬度達(dá)到一個峰值時，R2開始減小，直至去除寬度達(dá)到最大值180m；R2達(dá)到峰值時對應(yīng)的林地邊緣去除寬度為150m，即過渡帶寬度為30m。因此，當(dāng)溫度歸一化模型的擬合優(yōu)度達(dá)到最好時，林地邊緣過渡帶的最佳寬度為30m。

圖2 林地邊緣去除寬度和模型的擬合優(yōu)度之間關(guān)系

Fig.2TherelationshipbetweenthewidthofremovingwoodlandedgesandR2

表2展示了LandsatETM+和ASTER地表溫度模型中最高R2值及其對應(yīng)的線性回歸方程系數(shù)值；將最高R2值對應(yīng)線性回歸系數(shù)(B0,B1,B2)代入兩傳感器對應(yīng)的線性回歸方程：

ASTER：

T=18.472+0.014×HS-0.003×DEM

(9)

LandsatETM+：

T=18.503+0.017×HS-0.002×DEM

(10)

結(jié)果顯而易見，兩個方程表現(xiàn)出一定的相似性。首先，地表溫度(T)和山體陰影(HS)之間呈正相關(guān)關(guān)系，即光照條件好會導(dǎo)致地表溫度升高；其次，地表溫度(T)和高程(DEM)呈負(fù)相關(guān)關(guān)系，驗證了溫度隨高程增加而遞減的理論[14]。

表2 不同傳感器對應(yīng)的線性回歸方程參

Table2Thelinearregressionparameterscorrespondingtodifferentsensors

傳感器類型B0B1B2RR2LandsatETM+18．5030．017-0．0020．4680．219ASTER18．4720．014-0．0030．5050．255

注：B0、B1、B2為線性回歸方程系數(shù)，R2為擬合優(yōu)度。

4 討論

在本研究中，由于地表溫度模型只針對林地，忽略了地表其他地物的發(fā)射率問題，因此，在計算地表溫度時，掩膜處理了其他地物覆蓋。此外，不同樹的參數(shù)差異會引起林地比輻射率的差異(常溫下林地的比輻射率接近1[15])，故研究假設(shè)林地是一個相對均勻的土地利用/覆蓋單元，做出該假設(shè)主要考慮以下兩個因素的影響：第一，所選取感興趣區(qū)域內(nèi)樹梢的密度和高度不應(yīng)有很明顯的差異，差異較大，會明顯影響可見光和近紅外反射率值，產(chǎn)生較多的陰影，這會額外增加林地表面溫度反演的影響因素；第二，冠層覆蓋不應(yīng)該有太明顯的變化，因為冠層和背景值的改變也會導(dǎo)致反射率的變化[16]。因此，在進(jìn)行林地分類提取時，應(yīng)嚴(yán)格控制閾值的設(shè)置，但這會增強(qiáng)林地像元的齊性特征，而高的齊性特征又會增加式(8)中的系統(tǒng)誤差。

為了提高回歸模型的精度，減少因像元的齊性特征引起的系統(tǒng)誤差，研究中增加林地邊緣過渡帶，并將逐步掩膜過渡帶單位寬度后的結(jié)果用于回歸模型；通過R2評價回歸模型的精度，Landsat ETM+和ASTER熱紅外數(shù)據(jù)對應(yīng)回歸模型的R2值最高時的過渡帶寬度均為30 m，并在此條件下建立最佳回歸模型。從模型結(jié)果看，Landsat ETM+和ASTER熱紅外數(shù)據(jù)對應(yīng)回歸模型的R2值有差異，這種差異可能是兩種傳感器的空間分辨率和輻射分辨率的不同所引起的，因為像元空間尺度大小和熱輻射量化級的差異會影響像元的齊性特征，具體其間的定量關(guān)系還有待深入研究。研究最后利用ASTER遙感數(shù)據(jù)所建立地表溫度歸一化模型計算了研究區(qū)的地表溫度(圖3)，圖3展示了地表溫度與高程之間的變化趨勢，海拔較高的位置溫度偏低，海拔較低的位置溫度較高。

圖3 根據(jù)地表溫度歸一化模型(式(9))推演的研究區(qū)地表溫度

Fig.3 The estimated surface temperature of the whole area according to Eq.(9)

另外，在模型中除高程對地表溫度的影響外，坡度和坡向也是重要的潛在影響因素。Deng等[17]指出地表面溫度是緯度、坡度、坡向、地形遮擋作用和時間的函數(shù)，Pierce等[12]將這些地形參數(shù)的影響歸結(jié)為山體陰影指數(shù)，從反演的地表溫度圖(圖3)可以看出，山體陰影覆蓋較多的地方溫度較低，證實了地表溫度和山體陰影之間的關(guān)系。

5 結(jié)論

本研究使用Landsat ETM+和ASTER熱紅外數(shù)據(jù)驗證了地形對林地表面溫度反演的影響，得出的結(jié)論都是地表溫度隨著高程數(shù)值的增加而減小，山體陰影越強(qiáng)地表溫度越低。通過在林地邊緣設(shè)置一定寬度的過渡帶可以有效減弱林地邊緣的混合像元所導(dǎo)致的熱通量偏差，有利于增加地表溫度模型的精度。本研究涉及的地形和太陽位置等因素是許多學(xué)者在反演地表溫度時未曾考慮到的，引入這些因素建立地表溫度模型對復(fù)雜地形地表溫度的反演有一定的輔助作用。本研究沒有深入研究不同衛(wèi)星傳感器對地表溫度的具體影響，在后續(xù)的實驗中可進(jìn)一步研究地表溫度與傳感器的空間分辨率和輻射分辨率之間的定量變化關(guān)系，使復(fù)雜地形地表溫度的反演更加精確。

[1] ANDERSON M C,KUSTAS W P,NORMAN J M,et al.Mapping daily evapotranspiration at field to continental scales using geostationary and polar orbiting satellite imagery[J].Hydrology and Earth System Sciences,2011,15(1):223-239.

[2] DASH P,G?TTSCHE F M,OLESEN F S,et al.Land surface temperature and emissivity estimation from passive sensor data:Theory and practice-current trends[J].International Journal of Remote Sensing,2002,23(13):2563-2594.

[3] 劉向培，王漢杰，史小康，等.植被覆蓋對中國東部地面溫度變化的影響[J].解放軍理工大學(xué)學(xué)報(自然科學(xué)版)，2013，14(5):585-590.

[4] 李彩霞，周新國，孫景生，等.不均勻下墊面的地面光溫模擬研究[J].農(nóng)業(yè)機(jī)械學(xué)報，2013，44(12):140-147.

[5] 彭菁菁，劉強(qiáng)，柳欽火，等.多角度熱紅外亮溫值的模型擬合與應(yīng)用[J].紅外與毫米波學(xué)報，2011，30(4):361-371.

[6] 曹麗琴，張良培，李平湘，等.城市下墊面覆蓋類型變化對熱島效應(yīng)影響的模擬研究[J].武漢大學(xué)學(xué)報(信息科學(xué)版)，2008，33(12):1229-1232.

[7] 魏合理，宋正方.隨機(jī)風(fēng)場對自然地表紅外輻射統(tǒng)計分布的影響[J].紅外與毫米波學(xué)報，1995，14(6):424-428.

[8] PRICE J C.Quantitative aspects of remote sensing in the thermal infrared[A].Theory and Applications of Optical Remote Sensing[M].1989.578-603.

[9] FLORINSKY I V,KULAGINA T B,MESHALKINA J L.Influence of topography on landscape radiation temperature[J].International Journal of Remote Sensing,1994,15(16):3147-3153.

[10] GEIGER R,ARON R H,TODHUNTER P.The Climate Near the Ground[M].Lanham MD:Rowman and Littlefield Publishers,2003.

[11] 趙登忠，張萬昌，劉三超.基于 DEM 的地理要素 PRISM 空間內(nèi)插研究[J].地理科學(xué)，2004，24(2):205-211.

[12] PIERCE K B,LOOKINGBILL T,URBAN D.A simple method for estimating potential relative radiation(PRR)for landscape-scale vegetation analysis[J].Landscape Ecology,2005,20(2):137-147.

[13] QIN Z,KARNIELI A,BERLINER P.A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region[J].International Journal of Remote Sensing,2001,22(18):3719-3746.

[14] 傅抱璞.山地氣候[M].北京:北京出版社，1983.94-137.

[15] 梅安新，彭望琭，秦其明，等.遙感導(dǎo)論[M].北京：高等教育出版社，2001.19-23.

[16] 楊貴軍，柳欽火，劉強(qiáng)，等.植被冠層 3D 輻射傳輸模型及熱輻射方向性模擬[J].紅外與毫米波學(xué)報，2010，29(1):38-44.

[17] DENG Y,WILSON J P,GALLANT J C.Terrain analysis[A].The Handbook of Geographic Information Science[C].2008.417-435.

Influence of Topography on the Woodland Surface Temperature Retrieved from Thermal Infrared Remote Sensing Data

MENG Peng-yan1,XU Yuan-jin1,CHEN Xi2

(1.InstituteofMathematicalandRemoteSensingGeology,ChinaUniversityofGeosciences,Wuhan430074; 2.HubeiInstituteofGeologicalSurvey,Wuhan430034,China)

The purpose of this study was to assess the influence of topography on the woodland surface temperature,and to build a normalized model for surface temperature for complex and rugged woodland terrain.In this study,TERRA ASTER and Landsat ETM+images which transited in the same day were carried out.These models were mainly focused on the effects of complex topography(Digital Elevation Model(DEM) and illumination(HS))on surface temperature.So only homogeneous objects(woodland)were used for surface temperature modeling.In addition,this study analyzed the effects of the woodland edge mixed pixels on inversion accuracy model by using progressive masking.The results show that there is a negative correlation between surface temperature and DEM,and a positive correlation between surface temperature and HS.The influence of topography on surface temperature was confirmed in two types of satellite data used,and the precision of the model reached the highest when transition width of the woodland edge was 30 m.The research results have certain reference value for the complex terrain of land surface temperature inversion.

topography;thermal infrared remote sensing;surface temperature;woodland;DEM

2014-05-08；

2014-06-24

國家自然科學(xué)基金項目(41072246)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項資金項目(CUG120116)

孟鵬燕(1988-)，男，碩士研究生，主要研究方向是國土資源遙感、地質(zhì)遙感。E-mail：meng_pengyan@sina.com

10.3969/j.issn.1672-0504.2015.01.007

TP722.5;TP732.2;P404

A

1672-0504(2015)01-0032-05