馬 術(shù) 蒲 智

（新疆農(nóng)業(yè)大學(xué)計(jì)算機(jī)與信息工程學(xué)院 烏魯木齊 830052）

1 引言

紅外遙感是指遙感傳感器的工作波段限制在紅外范圍內(nèi)。熱紅外信息可以被傳感器捕獲并保存，熱紅外信息被用于區(qū)分地面物體并反演常用數(shù)據(jù)，如氣溫以及熱量慣性。地表溫度（LST）即通過(guò)熱紅外遙感技術(shù)得到的地溫［2］。其是水量和能量平衡的組成部分［2］。準(zhǔn)確的測(cè)量和估算有利于保護(hù)生態(tài)環(huán)境，保護(hù)水資源和改善生活的質(zhì)量［3］。

遙感應(yīng)用在我們生活中，通過(guò)熱紅外信息反演地表溫度的算法已經(jīng)成熟［3～4］。自20世紀(jì)中后期至今，通過(guò)遙感數(shù)據(jù)來(lái)反演地表溫度，許多學(xué)者研究出了陸地表面溫度反演算法［4～6］。2001 年，覃志豪等得出了一種基于表面熱輻射傳導(dǎo)方程的簡(jiǎn)單而準(zhǔn)確的單窗算法。并且單窗口算法直接包括大氣和表面的影響在計(jì)算公式中［7］；2003 年，Sobrino 和Jimenez-Munoz［7］為任何傳感器找到了通用的SC 方法，并使用TM6數(shù)據(jù)對(duì)。陸地表面溫度的反演產(chǎn)生了更高的精度結(jié)果。隨后的研究基于單通道算法的比較。孟憲紅采用大氣校正法，單通道算法和單窗算法，利用遙感數(shù)據(jù)反演甘肅省金塔地區(qū)的地表溫度［8］。單窗算法因其準(zhǔn)確和簡(jiǎn)單的優(yōu)勢(shì)，是目前進(jìn)行地表溫度反演的主要算法。

本文以烏魯木齊市為研究區(qū)，利用Landsat8 數(shù)據(jù)，反演研究區(qū)的陸地表面溫度，結(jié)合大氣校正方法和單窗算法，分析并比較結(jié)果。

2 研究區(qū)域和數(shù)據(jù)來(lái)源

2.1 研究區(qū)概況

烏魯木齊是新疆的首府，位于天山北部和準(zhǔn)格爾盆地南緣［10］，位于86°37′E～88°58′E，42°45′N(xiāo)～44°08′N(xiāo)，屬于中溫帶干旱氣候，年平均降水量194mm［10］，春秋和秋季短，冬季和夏季較長(zhǎng)，季節(jié)分布不均，熱量和濕度變化劇烈，晝夜溫差較大，逆溫層出現(xiàn)在冬季。最溫暖的七八月這兩個(gè)月的平均氣溫為25.7℃，1 月份最冷，平均氣溫為-15.2℃，年平均氣溫7.2℃，年平均日照時(shí)間為4.6h。由于烏魯木齊位于42°45′N(xiāo)～44°08′N(xiāo) 之間，夏天太陽(yáng)較高，地處內(nèi)陸，遠(yuǎn)離海洋，再加上城市熱島效應(yīng)的影響，天氣炎熱，極端溫度達(dá)到40.8℃。研究區(qū)域圖如圖1所示。

圖1 研究區(qū)位置及Landsat8影像

2.2 研究數(shù)據(jù)和處理

本文選取2016 年7 月28 日的Landsat8 遙感 數(shù)據(jù)，條帶號(hào)為142，行號(hào)為30。Landsat8數(shù)據(jù)從地理空間數(shù)據(jù)云網(wǎng)站下載，其云量為0.63，質(zhì)量較好。該衛(wèi)星于2013年2月11日發(fā)射，有兩個(gè)傳感器，操作陸地成像儀（OLI）和熱紅外傳感器（TIRS）［9］。除Landsat5 的所有光譜波段之外，OLI 還增加了四個(gè)波段［35］。Landsat8的每個(gè)波段的輻射分辨率從8位提高至12 位，這增加了圖像的灰度級(jí)并改善了信噪比［10］。TIRS 有band10 和band11，分辨率為100米，這使得大氣校正更容易，并且可以用劈窗算法［11］或ENVI 的熱傳感器校正工具進(jìn)行熱紅外校正。三個(gè)衛(wèi)星參數(shù)在表1中進(jìn)行了比較。

研究區(qū)內(nèi)的最高氣溫、最低氣溫、平均高溫和平均低溫等數(shù)據(jù)主要通過(guò)歷史天氣查詢網(wǎng)（mip.li?shi.tianqi.com）查詢。Landsat8 數(shù)據(jù)使用WGS84 坐標(biāo)系。文中研究區(qū)域行政邊界矢量圖是從相關(guān)網(wǎng)站下載的。使用ENVI Radiometric Calibration 進(jìn)行輻射定標(biāo)，并選擇熱紅外數(shù)據(jù)。由于band11 精度不高，因此選擇中心波長(zhǎng)為10.90μm 的波段。由于要排除大氣輻射的影響，因此選擇輻射定標(biāo)類(lèi)型的亮度值。公式如下：

ML和AL分別是增益參數(shù)和偏移參數(shù)，Qcal是灰度值［12］。兩個(gè)熱紅外波段輻射亮度公式為

3 地表溫度反演算法

3.1 大氣校正法

基本原理：去除大氣對(duì)地表熱輻射產(chǎn)生的誤差，把熱輻射強(qiáng)度轉(zhuǎn)換為相應(yīng)的地表溫度。熱紅外輻射亮度值Lλ由三部分組成：通過(guò)大氣的地面的真實(shí)輻射的能量，輻射向上的能量和向下輻射到達(dá)地面后反射的能量。公式如下：

在等式（3）中，τ 是熱紅外波段中的大氣透過(guò)率，ε 是地表比輻射率，Ts是陸地表面真實(shí)溫度（K），B(Ts)是黑體輻射亮度。溫度為T(mén)的黑體的輻射亮度B(Ts)：

在等式（4）中，透過(guò)率τ，大氣向上輻射亮度L↑和大氣向下輻射亮度L↓，這三個(gè)參數(shù)可以在NASA官方網(wǎng)站中信息獲得［16］。烏魯木齊市圖像的成像時(shí)間為2016 年7 月28 日04 時(shí)50 分、圖像的中心緯度為86.6833E 和43.7667N，氣壓為800 百帕，相對(duì)濕度為39%，得到的大氣參數(shù)透過(guò)率τ 為0.7，大氣向上輻射亮度L↑為1.29，大氣向下輻射亮度L↓為0.75。Ts用普朗克公式獲得：

對(duì)于band10，K1=774.89 W/（㎡*μm*sr），K2=1321.08K。

表1 Landsat5，Landsat7和Landsat8數(shù)據(jù)參數(shù)對(duì)比

3.2 單窗算法

大氣校正法受大氣剖面數(shù)據(jù)和探測(cè)數(shù)據(jù)的限制，覃志豪提出以TM 數(shù)據(jù)作為數(shù)據(jù)源的地表溫度反演算法，即單窗算法［13］。公式如下：

等式（6）中，Ta是大氣平均作用溫度（K），Tb是亮度溫度（K），Ts是反演溫度（K），a、b 值是-67.35535、0.45861；C 和D 是中間量，由式（7）和（8）計(jì)算。Ta是從式（9）獲得：

在公式中，T0為近地面溫度，單位為K，τ 為大氣透射率，ε 是地表比輻射率。從當(dāng)日的氣象報(bào)告中獲得烏魯木齊市當(dāng)時(shí)近地面氣溫是24℃，算出Ta是291.23K。大氣透射率τ 需通過(guò)大氣含水量W 與表2一起估算。

表2 大氣透射率估算方程

計(jì)算絕對(duì)水氣壓e：

在式（10）中，T0是溫度（K），RH 是相對(duì)濕度，e是以千帕為單位。利用地表水壓力和絕對(duì)水氣壓e和楊景梅確定的大氣含水量之間的關(guān)系計(jì)算大氣水汽含量。公式如下：

在式（11）中，e 是絕對(duì)水壓，a1和a0是0.1978和0.1788。在該式中，e 單位是100Pa，通過(guò)公式參數(shù)得出W=3.06g/cm2。

使用楊槐Landsat8 數(shù)據(jù)的透過(guò)率和水汽計(jì)算大氣透過(guò)率τ：

為估計(jì)亮度溫度，首先通過(guò)式（2）將像元的灰度值轉(zhuǎn)化為熱輻射強(qiáng)度值。其次將熱輻射強(qiáng)度值轉(zhuǎn)化為亮度溫度，公式如下：

式（13）中，K1，K2值與式（5）中相同。

3.3 地表比輻射率的計(jì)算

計(jì)算地表比輻射率對(duì)于陸地表面溫度的反演是必要步驟，其他方法也要有此步驟。TIRS 的熱紅外波段與TM/ETM+6 熱紅外波段具有相近的波普范圍［14］，本文使用與TM/ETM+6 一樣的地表比輻射率計(jì)算方法，公式如下：

其中，Pv是植被覆蓋度，公式如下［9］：

NDVI 是歸一化植被指數(shù)，NDVIveg是植被完全覆蓋的NDVI 值，NDVIsoil是裸土的NDVI 值或沒(méi)有植被覆蓋區(qū)域的NDVI 值，NDVIveg=0.70，NDVIsoil=0.05［7］。

4 結(jié)果與分析

通過(guò)研究區(qū)域的Landsat8 數(shù)據(jù)（波段是B10 波段），利用單窗算法和大氣校正法得到烏魯木齊地表溫度反演圖像，并利用ArcMap軟件繪圖［15］。圖2是研究區(qū)亮溫圖；圖3 是輻射傳導(dǎo)方程法地表溫度圖；圖4是單窗算法地表溫度圖。

圖2 烏魯木齊市亮溫圖（℃）

圖3 烏魯木齊市大氣校正法地表溫度反演圖（℃）

從圖2、圖3 和圖4 可以看出，由亮度溫度得到的陸地表面溫度與輻射傳導(dǎo)方程法和單窗算法的陸地表面溫度反演的空間分布基本相同，但變化程度略有不同。由于城市擴(kuò)張，工業(yè)產(chǎn)熱等因素，導(dǎo)致城市的溫度較高；裸地與綠地相比，溫差由比熱容和含水量的差異引起［13］。白天日照時(shí)水體溫度較低［13］；晚上，由于水的熱容量較大，溫度較高。

圖4 烏魯木齊市單窗算法地表溫度反演圖（℃）

根據(jù)反演溫度的影像統(tǒng)計(jì)，圖2 是亮溫圖，最高溫度是15.69℃，最低溫度是0℃，平均溫度是3.71℃；圖3 是大氣校正法地表溫度反演圖（LST-RTE），最高溫度是36.21℃，最低溫度是-34.02℃，平均溫度是3.31℃；圖4是單窗算法地表溫度反演圖（LST-MW），最高溫度是42.32℃，最低溫度是-39.16℃，均值溫度是3.66℃。LST-RTE影像平均溫度和亮度溫度的平均值相差0.4℃，LST-MW 影像的平均溫度與亮度溫度的平均溫度相差0.05℃，LST-RTE 影像的平均溫度與LST-MW影像的平均溫度相差0.35℃。

上述分析表明，大氣校正法和單窗算法得到的平均溫度高于亮度溫度。將大氣校正法和單窗算法得到的平均溫度與亮度溫度比較，單窗算法反演精度優(yōu)于大氣校正法，因?yàn)檩椛鋫鬏斔惴ǜ蕾?lài)于探空數(shù)據(jù)和大氣剖面數(shù)據(jù)，很少數(shù)據(jù)符合條件，本文由美國(guó)宇航局網(wǎng)站獲取精確度較低的大氣剖面數(shù)據(jù)；單窗算法的參數(shù)主要取決于大氣含水量，具有較高的精度。

5 結(jié)語(yǔ)

本文基于Landsat8 數(shù)據(jù)，采用單通道算法和大氣校正法兩種算法，反演烏魯木齊市陸地表面溫度，并對(duì)研究結(jié)果和研究區(qū)的亮溫進(jìn)行了比較。得到兩點(diǎn)結(jié)論。

1）大氣校正法和單窗算法的陸地表面溫度和亮溫反演在地理上的分布很接近，但也有差別；

2）反演的平均溫度與亮度溫度相比，大氣校正法低出0.4℃，單窗算法低出0.05℃，單窗算法和大氣校正法平均溫度相差0.35℃。