基于ETM+數(shù)據(jù)的和林格爾縣地表溫度反演

王云川

【摘 要】文章采用基于影像的反演方法（IB方法），以LandsatETM+為數(shù)據(jù)源，成功反演出研究區(qū)和林格爾縣地表溫度的分布狀況。結(jié)果顯示：（1）水體、農(nóng)田、林地地表溫度較低，裸地的溫度較高，與實(shí)際情況相符；（2）裸地的溫度分布呈現(xiàn)東高西低、南高北低的趨勢(shì)，與該區(qū)地形走勢(shì)一致。

【關(guān)鍵詞】地表溫度反演；ETM+數(shù)據(jù)；和林格爾縣

地表溫度是圈層系統(tǒng)中水、熱、碳等物質(zhì)交換的核心信息和主導(dǎo)因子，它在太陽(yáng)輻射-大氣-地面間物質(zhì)和能量的交換過(guò)程中扮演重要的角色，也是平衡地球圈層水熱平衡的關(guān)鍵因子。精確地反演地表溫度將促進(jìn)農(nóng)作物缺水和旱情預(yù)測(cè)、全球氣溫變化、數(shù)字天氣預(yù)報(bào)等研究。傳統(tǒng)的地表溫度監(jiān)測(cè)是通過(guò)地面氣象站來(lái)完成的，費(fèi)時(shí)費(fèi)力，且局限性較大。

隨著衛(wèi)星遙感技術(shù)的快速發(fā)展，地表溫度的大尺度、多時(shí)相監(jiān)測(cè)成為可能。廣泛使用熱紅外激光和雷達(dá)數(shù)據(jù)，紅外激光雷達(dá)使用主動(dòng)熱紅外輻射源，提高了空間分辨率，拓展了遙感信息源。推動(dòng)了熱紅外遙感的新發(fā)展，同時(shí)以此來(lái)反演的地表溫度數(shù)據(jù)精度也會(huì)大大提高。

一、研究區(qū)概況、數(shù)據(jù)來(lái)源

（一）研究區(qū)概況。和林格爾縣位于內(nèi)蒙古自治區(qū)中南部，是呼和浩特市所轄縣旗之一，是內(nèi)蒙古高原向黃土高原的過(guò)渡地帶，總體地形呈南高北低、東高西低的地形趨勢(shì)。該區(qū)屬于溫帶大陸性季風(fēng)氣候，夏季時(shí)短而溫?zé)岫韭L(zhǎng)而寒冷。年平均氣溫在6.2℃左右，一月平均氣溫為-12.8℃，七月平均氣溫為22.1℃，年平均降水量為392.8mm。

（二）數(shù)據(jù)來(lái)源與預(yù)處理。文章采用的數(shù)據(jù)源為L(zhǎng)andsat7 ETM+于 2002年9月16日獲取的影像，條帶號(hào)行列號(hào)為126/32。數(shù)據(jù)來(lái)自地理空間數(shù)據(jù)云，級(jí)別為L(zhǎng)evel 1T。在進(jìn)行反演算法前需要對(duì)影像進(jìn)行相關(guān)處理，將各像元的亮度值定標(biāo)到電磁輻射值，并進(jìn)行大氣校正，以消除水汽和地表等因素的干擾。

二、研究方法

利用遙感影像研究大范圍的地面溫度通常采用基于影像的反演算法（ Image-based Method）。算法的過(guò)程如下：

（一）計(jì)算地表比輻射率。地物的比輻射率用來(lái)表示地物向外輻射電磁波的能力。它受地表物體組成的制約，而且與物體的表面狀態(tài)（高度，走向，地面粗糙度等）及物理性質(zhì)有關(guān)?；谟跋竦姆囱菟惴ㄓ?jì)算地表比輻射率時(shí)通常根據(jù)與植被覆蓋度有關(guān)的經(jīng)驗(yàn)公式估算。植被覆蓋度與地物類型密切相關(guān)，文章對(duì)NDVI數(shù)據(jù)進(jìn)行密度分割，以獲取研究區(qū)遙感影像分為植被、水域、裸地3種類型。NDVI的計(jì)算公式如下：

NDVI=（TM4-TM3）/（TM4+TM3）

式中，TM3和TM4分別為大氣校正后的第三波段和第四波段。

植被覆蓋度在ENVI的波段運(yùn)算中完成，其代碼去下：

b2=（b1 gt NDVIv）×1+（b1 lt NDVIS）×0+（b1 ge NDVIS and b1 le NDVIv）×（（b1-NDVIS）/（ NDVIv-NDVIS））

式中，b1為NDVI圖像；NDVIv為植被的NDVI下限閾值，取0.58；NDVIS為水體的NDVI上限閾值，取0.14；b2為運(yùn)算得到的植被覆蓋度圖像。

水體的比輻射率賦值為0.995，自然表面和居民地像元的比輻射率則根據(jù)植被覆蓋度運(yùn)算得到。其代碼如下：

b3=（b1 le 0.14）×0.995+（b1 gt 0.14and b1 lt 0.58）×（0.9589 +0.086×b2 - 0.0671×b2^2）+（b1 ge 0.58）×（0.9625 + 0.0614×b2- 0.0461×b2^2）

式中，b1為NDVI值；b2為植被覆蓋度值。b3為運(yùn)算得到的地表比輻射率數(shù)據(jù)。

（二）計(jì)算黑體的輻射亮度值。衛(wèi)星傳感器接收到的熱紅外輻射亮度值Lλ由三個(gè)部分組成：大氣向上輻射亮度L↑，地面長(zhǎng)波熱輻射經(jīng)過(guò)大氣層之后到達(dá)衛(wèi)星傳感器的能量，大氣向下輻射輻射亮度L↓。因此，衛(wèi)星傳感器所接收到的輻射傳輸方程如下：

Lλ = [b3·b5 + （1-b3）L↓]·τ + L↑

則地面溫度T的黑體輻射亮度值b5的公式為

b5 = [Lλ - L↑- τ·（1-b3）L↓]/τ·b3

在NASA官網(wǎng)中通過(guò)申請(qǐng)得到該影像的輻射參數(shù)，大氣在熱紅外波段的透過(guò)率τ為0.83，大氣向上輻射亮度L↑為1.12W/（m2·sr·μm），大氣向下輻射亮輻射亮度L↓為1.89W/（m2·sr·μm）。因此，面溫度T的黑體輻射亮度值的波段運(yùn)算代碼如下

b5=（b4-1.12-0.83×（1-b3）×1.89）/（0.83×b3）

式中，b4為經(jīng)過(guò)ENVI輻射定標(biāo)并進(jìn)行FLAASH大氣校正后的熱紅外輻射亮度值；b5為黑體在熱紅外波段的輻射亮度值。

（三）反演地表溫度。在獲取溫度為T的黑體在熱紅外波段的輻射亮度b5后，可根據(jù)普朗克函數(shù)的反函數(shù)，求得地表真實(shí)溫度TS：

TS= K2/ln（K1/B（TS）+1）

在本次研究中，K1 =666.09W/（m2×sr×μm），K2= 1282.71K。利用波段運(yùn)算即可得到真實(shí)的地表溫度值b6，代碼如下：

b6=（1282.71）/alog（666.09/b5+1）-273

三、結(jié)果分析

地表溫度反演結(jié)果顯示，從地表覆蓋狀況來(lái)看，水體、農(nóng)田、林地的地表溫度值較低，裸地的地表溫度值較高；從空間分布狀況來(lái)看，東高西低、南高北低，與該區(qū)域的地形分布和海拔走勢(shì)一致；該區(qū)域的植被覆蓋較少，荒土和裸巖分布廣泛，農(nóng)田主要沿河谷分布（兩者同屬低值區(qū)且分布走向一致）。

參考文獻(xiàn)

[1] 鄭文武，曾永年.地表溫度的多源遙感數(shù)據(jù)反演算法對(duì)比分析[J].地球信息科學(xué)學(xué)報(bào)，2011，13（06）：840-847.

[2] Masuoka E，F(xiàn)leig A，Wolfe R E，et al.Key characteristics of MODIS data products[J].Remote Sensing， 1998，36（04）：1313-1323.

[3] 胡文星.基于不同算法的溫度反演比較[J].電子制作， 2014（21）：253-254.

[4] 廖麗.千島湖地形調(diào)查及湖面溫度場(chǎng)衛(wèi)星數(shù)據(jù)的反演[D].杭州：杭州電子科技大學(xué)，2014.