張麗娟 ZHANG Li-juan；溫興平 WEN Xing-ping；周楊 ZHOU Yang；王軍 WANG Jun

（①昆明理工大學(xué)國(guó)土資源工程學(xué)院，昆明 650093；②云南省礦產(chǎn)資源預(yù)測(cè)評(píng)價(jià)工程實(shí)驗(yàn)室，昆明 650093）

（①Faculty of Land Resource Engineering，Kunming University of Science and Technology，Kunming 650093，China；②Mineral Resources Prediction and Evaluation Engineering Laboratory of Yunan Province，Kunming 650093，China）

0 引言

地表溫度[1-2]是研究區(qū)域地表能量平衡和資源環(huán)境變化的重要參數(shù)之一，它在眾多領(lǐng)域有著廣泛的需求，因此準(zhǔn)確獲取其分布的時(shí)空狀況具有非常重要的意義。要想獲得大區(qū)域的地表溫度，就需要通過(guò)遙感手段，因?yàn)闇囟鹊倪b感反演是遙感技術(shù)的關(guān)鍵環(huán)節(jié)之一，具有非常重要的現(xiàn)實(shí)意義。

MODIS（中分辨率成像光譜儀）[3]是新一代地球觀測(cè)傳感器，搭載于NASA的地球觀測(cè)系統(tǒng)（EOS）TERRA和AQUA兩顆衛(wèi)星之上，其數(shù)據(jù)在全世界范圍都可免費(fèi)接收。MODIS數(shù)據(jù)采用HDF格式進(jìn)行存儲(chǔ)、發(fā)布。HDF的數(shù)據(jù)結(jié)構(gòu)是一種分層式數(shù)據(jù)管理結(jié)構(gòu)，存儲(chǔ)了MODIS數(shù)據(jù)以及相應(yīng)元數(shù)據(jù)，提供了MODIS數(shù)據(jù)輻射校正所需的相應(yīng)衛(wèi)星參數(shù)。

1 研究區(qū)概況

會(huì)澤縣位于云南省東北部、金沙江東岸、曲靖市西北部，全縣占地面積6077平方公里，地處東經(jīng)103°03′～103°55′、北緯 25°48′～27°04′之間。而縣城所在地金鐘鎮(zhèn)，海拔達(dá)到2120米，距離曲靖市區(qū)有187公里，距省會(huì)昆明205公里。

2 數(shù)據(jù)來(lái)源

本文采用2010年1月8日12時(shí)25分 (本地時(shí)間)Terra(AM)衛(wèi)星MODIS數(shù)據(jù)，影像覆蓋范圍為北緯26°～28°N，東經(jīng)102°～103°E。選取云南省會(huì)澤縣為研究區(qū)域，沒(méi)有云層覆蓋，無(wú)雪覆蓋，質(zhì)量良好。本次試驗(yàn)選取MODIS的第31波段（空間分辨率為1km）為試驗(yàn)對(duì)象。

3 數(shù)據(jù)處理

3.1 單窗算法反演溫度 單窗算法是覃志豪[4-6]針對(duì)Landsat TM影像在只有一個(gè)熱紅外波段TM6的時(shí)候提出的，該算法最大的優(yōu)點(diǎn)就是大氣參數(shù)少，但是精度高。因此，在對(duì)MODIS數(shù)據(jù)利用ENVI自帶的GEOREFERENCE MODIS模塊作幾何校正，去除蝴蝶結(jié)效應(yīng)的基礎(chǔ)上，運(yùn)用單窗算法對(duì)MODIS的31通道影像反演了會(huì)澤地區(qū)的溫度分布。

如果大氣透射率τ6，大氣平均作用溫度Ta和地表輻射率ε6已知，則可用如下單窗算法從像元的亮點(diǎn)溫度值T6中推算該像元的實(shí)際地表溫度TS：

式中，Ts的單位為K；a6和b6為常量，在一般情況下（即當(dāng)?shù)乇頊囟仍?～70℃范圍內(nèi)時(shí)），取值a6=-67.355351，b6=0. 458606；C6和D6是中間變量分別用下式表示：

因此，只要知道了參數(shù) Ta、ε6、τ6，就可用上述單窗算法推算任何像元的實(shí)際地表溫度。

MODIS的第31通道覆蓋的波段范圍與Landsat的波段范圍基本是一致的，同樣可以用這種思想來(lái)推導(dǎo)適用于MODIS通道的溫度反演方程。大氣輻射傳輸能量平衡方程如下：

為了簡(jiǎn)化計(jì)算，覃志豪將Planck方程進(jìn)行泰勒級(jí)數(shù)展開求近似表達(dá)式。根據(jù)Planck方程，輻射率與溫度的關(guān)系為：

3.2 亮度溫度的計(jì)算

MODIS影像是用DN值表示的，然后再用Planck方程求解星上亮溫。MODIS第31波段的輻射強(qiáng)度值計(jì)算公式如下：

式中，rad31是MODIS第31波段的熱輻射強(qiáng)度；band31是MODIS第31波段的DN值；scale31和offset31是MODIS第31波段的輻射定標(biāo)常量，可從MODIS數(shù)據(jù)集的屬性數(shù)據(jù)中查出。

我們可以使用Planck方程求解出星上亮度溫度。對(duì)31波段應(yīng)用Planck方程，并化簡(jiǎn)可得31波段的亮溫計(jì)算公式：

式中：K31,1=1304.413871，K31,2=729.541636。

反演后的研究區(qū)亮溫如圖1。

圖1 2010年1月8日12時(shí)25分亮溫反演影像

本文用到的MODIS影像是攝取于1月上旬，對(duì)會(huì)澤地區(qū)而言，通過(guò)利用上述方法對(duì)31波段進(jìn)行地面溫度反演發(fā)現(xiàn)，31波段的星上亮度溫度非常接近。通過(guò)用單通道方法反演，第31波段的最高溫度為22.85℃，平均溫度為8.15℃。符合影像攝取時(shí)的當(dāng)?shù)貙?shí)際溫度范圍，這種大面積溫度分布規(guī)律在TM影像上呈現(xiàn)不出來(lái)，但是反映大區(qū)域的溫度分布差異是MODIS影像的一個(gè)優(yōu)勢(shì)。如圖2。

圖2 溫度/K

圖3 輻射強(qiáng)度與溫度變化的關(guān)系圖

如圖3，散點(diǎn)圖中可以清晰的看出研究區(qū)的溫度區(qū)間為-8℃-22℃之間，為我們研究該區(qū)域給出了直觀的分析。

由圖4會(huì)澤區(qū)空間溫度分布圖看出，溫度的地域分布很有特點(diǎn)，總體的特征是隨著維度的地表溫度增大的，而溫度大致還可分為3個(gè)區(qū)間：＜0℃、0℃～11℃、＞11℃。這種地溫產(chǎn)生的的原因主要是由一些云塊引起的異常，是因?yàn)檩^濃厚的云層遮蔽，使得傳感器接受到地面極低的熱輻射量，從而反演出的溫度值極低。

圖4 會(huì)澤區(qū)空間溫度分布圖

4 結(jié)論

本文利用基于單窗算法的地表溫度反演算法，通過(guò)采用MODIS的第31波段提取地表比輻射率，結(jié)合熱紅外波段的亮溫信息，反演了會(huì)澤地區(qū)的地表溫度。結(jié)果表明：反演精度比較高。同時(shí)對(duì)MODIS的第31通道的參數(shù)進(jìn)行了計(jì)算，對(duì)會(huì)澤地區(qū)的整個(gè)MODIS用該方法進(jìn)行了反演，并對(duì)會(huì)澤地區(qū)反演的地表溫度進(jìn)行了分析。結(jié)果表明反演的地表溫度比較合理。該方法適宜于該地區(qū)。

[1]Becker,F,Li,Z,L..Towards a local split window method over land surface[J].Intermati-onalJourmal of Remote Sensing,1990(11):369-393.

[2]LijuanZhang,XingpingWen,JunWang,etal.Applied Materiala and Manufacturing[C].Trans Tech Publications,Vols.2014，889-890：1634-1637.

[3]覃志豪.Li Wenjuan.張名華等.單窗算法的大氣參數(shù)估計(jì)方法[J].國(guó)土資源遙感，2003（56）.