摘要：地表溫度是描述地表過(guò)程變化的重要參數(shù)，在全球氣候變化、地表輻射能量平衡、城市熱環(huán)境、資源環(huán)境監(jiān)測(cè)等研究領(lǐng)域中扮演著重要的角色。隨著遙感技術(shù)的發(fā)展，利用遙感技術(shù)估算陸地表面溫度的方法不斷成熟，熱紅外遙感技術(shù)是核心方法，眾多熱紅外數(shù)據(jù)中Landsat系列憑借連續(xù)性好、易獲取的優(yōu)勢(shì)被廣泛應(yīng)用于地表溫度反演研究。本文梳理了目前利用Landsat 8熱紅外數(shù)據(jù)反演地表溫度的基礎(chǔ)理論、反演算法以及反演參數(shù)計(jì)算等，并指出基于Landsat 8數(shù)據(jù)的地表溫度反演有待解決的問(wèn)題，最后對(duì)利用Landsat 8數(shù)據(jù)反演地表溫度研究進(jìn)行了展望。

關(guān)鍵詞：地表溫度;熱紅外遙感技術(shù);Landsat 8影像;反演

中圖分類號(hào)：X830.2 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：2095-672X（2020）05-0-02

DOI：10.16647/j.cnki.cn15-1369/X.2020.05.149

Discussion on land surface temperature inversion based on Landsat 8 imagery

Li Yanan1，2，3

（1.Shaanxi Provincial Land Engineering Construction Group Co.，Ltd.，Xian Shaanxi 710075，China;2.Institute of Land Engineering and Technology， Shaanxi Provincial Land Engineering Construction Group Co.，Ltd.，Xian Shaanxi 710075，China;3.Key Laboratory of Degraded and Unused Land Consolidation Engineering，the Ministry of Natural Resources，Xian Shaanxi 710075，China）

Abstract：Surface temperature is an important parameter describing the change of surface processes， and it plays an important role in research fields such as global climate change， surface radiant energy balance，urban thermal environment， and resource environment monitoring.With the development of remote sensing technology，the method of estimating land surface temperature using remote sensing technology has continued to mature.Thermal infrared remote sensing technology is the core method.Among the many thermal infrared data，the Landsat series has been widely used for land surface temperature reflection due to its good continuity and easy access.Act Research. In this paper， the basic theory，inversion algorithm and inversion parameter calculation of Landsat 8 thermal infrared data inversion are summarized.The problems of land surface temperature inversion based on Landsat 8 data are pointed out.Finally，Landsat 8 data The inversion of surface temperature research is prospected.

Key words：Land surface temperature;Thermal infrared remote sensing technology;Landsat 8 image;Inversion

陸地表面溫度LST（Land Surface Temperature）即地表溫度，在地表與大氣能量交換中起著關(guān)鍵作用，是分析和模擬區(qū)域尺度地表過(guò)程的關(guān)鍵參數(shù)[1]。充分了解區(qū)域乃至全球水-熱循環(huán)的內(nèi)在機(jī)理與聯(lián)系，并定量化模擬這些循環(huán)過(guò)程是地理學(xué)研究的核心和熱點(diǎn)。地表溫度是連接地表能量平衡與水分循環(huán)的關(guān)鍵參數(shù)之一。在全球氣候變化、地表輻射能量平衡、城市熱環(huán)境、資源環(huán)境監(jiān)測(cè)等研究領(lǐng)域具有極為重要的意義。被反射的太陽(yáng)輻射和被地面吸收的太陽(yáng)輻射這兩部分熱能影響著地表溫度，使地面增熱，對(duì)地面的溫度進(jìn)行測(cè)量后得到的溫度就是地表溫度。

地表溫度的準(zhǔn)確估算是國(guó)內(nèi)外研究的熱點(diǎn)。傳統(tǒng)的地表溫度觀測(cè)方法大多是基于單點(diǎn)的局地測(cè)量，主要包括玻璃液體溫度計(jì)、鉑電阻測(cè)溫法和紅外測(cè)溫法。氣象站主要通過(guò)前兩種接觸式儀器進(jìn)行測(cè)量地表溫度。傳統(tǒng)的觀測(cè)手段獲取的地表溫度信息適用于單點(diǎn)及田間尺度的應(yīng)用，對(duì)于實(shí)際的區(qū)域尺度研究具有一定的局限性，難以使用常規(guī)的空間插值直接應(yīng)用于區(qū)域尺度的研究，不僅會(huì)引起較大的精度誤差，而且不能滿足實(shí)際的需求。伴隨著新型對(duì)地觀測(cè)衛(wèi)星的星載傳感器的種類、探測(cè)精度、實(shí)用性能的不斷提升，多時(shí)相、大范圍、多波譜、高分辨率的遙感數(shù)據(jù)的快速獲取，使得遙感技術(shù)成為高時(shí)效、高精度獲取地球表面區(qū)域尺度溫度分布的最有效手段。利用遙感技術(shù)可以獲得區(qū)域尺度下多種水熱情況參數(shù)，有利于更精確的描述區(qū)域乃至全球地氣之間的水熱交換過(guò)程，使得區(qū)域和全球尺度下的地表溫度反演精度適用于實(shí)際的科研生產(chǎn)需要。

1地表溫度反演理論

普朗克黑體輻射定律、斯特藩-玻爾茲曼定律和維恩位移定律這三大定律是熱輻射基本定律，地表溫度反演算法的提出均是基于此。熱紅外輻射傳輸方程是熱紅外遙感反演地表溫度的基礎(chǔ)，熱紅外輻射傳輸方程認(rèn)為衛(wèi)星傳感器接收到總輻射包括地表熱輻射經(jīng)大氣衰減后到達(dá)傳感器的輻射、大氣直接熱輻射即大氣上行熱輻射和經(jīng)地表發(fā)射后通過(guò)大氣傳輸?shù)絺鞲衅鞯妮椛浼创髿庀滦袩彷椛洹?/p>

2 地表溫度反演算法

地表溫度反演算法主要有單通道算法、多通道算法、多角度算法、多時(shí)相算法、高光譜反演算法等，以上算法均可用于反演區(qū)域地表溫度。載有熱紅外波段的衛(wèi)星遙感數(shù)據(jù)是獲取區(qū)域地表溫度的主要數(shù)據(jù)來(lái)源，目前常用的載有熱紅外波段的衛(wèi)星傳感器包括：NOAA搭載的AVHRR3;TERRA搭載的MODIS;CBERS搭載的IRMSS;Landsat 5/7搭載的TM/ETM+;Landsat 8搭載的TIRS;FY-2和FY-3搭載的VIRR、HJ-1B搭載的IRS等。

眾多熱紅外數(shù)據(jù)中，Landsat系列具有連續(xù)性好、空間分辨率高以及數(shù)據(jù)易獲取的特點(diǎn)，被廣泛應(yīng)用于地表溫度研究中。Landsat 8衛(wèi)星成功發(fā)射以來(lái)，為地表溫度的研究帶來(lái)了新的機(jī)遇，傳感器TIRS包含第10和第11兩個(gè)波段，但由于第11波段正處于不穩(wěn)定階段，所以很多學(xué)者將第10波段作為單熱紅外波段進(jìn)行地表溫度的研究，同時(shí)這使得劈窗算法不適用以Landsat 8作為基礎(chǔ)數(shù)據(jù)進(jìn)行地表溫度的反演，本文將主要介紹覃志豪單窗算法和Juan C. Jiménez-Mu?oz 單通道算法。

3 地表溫度反演參數(shù)計(jì)算

3.1 亮溫值

亮溫值Ti是指觀測(cè)物體和黑體輻射出輻射能量相等時(shí)，黑體的溫度就被稱為觀測(cè)物體的亮溫值，可以作為一個(gè)指標(biāo)來(lái)衡量物體溫度，但不是物體的真實(shí)溫度，亮溫值的計(jì)算公式如下[4]：

式中：Ti是為波段i對(duì)應(yīng)的亮度溫度（K）;Lλ為熱輻射強(qiáng)度值（W·m-2·sr-1·μm-1）;h為普朗克常數(shù)、c為光速、k為玻爾茲曼常數(shù)，Landsat 8 TIRS 10波段的中心波長(zhǎng)λ為10.9，根據(jù)公式可以得出常數(shù)K1為774.89 W·m-2·sr-1·μm-1，K2為1321.08K。

3.2 地表比輻射率

地表比輻射率指相同溫度條件下，地表發(fā)射的輻射量與黑體發(fā)射的輻射量的比值?；诜诸惖谋容椛渎史椒?、基于NDVI的比輻射率方法以及獨(dú)立溫度光譜指數(shù)法是目前地表比輻射率計(jì)算的主流方法。其中基于NDVI的比輻射率方法具有易獲取、限制少的優(yōu)勢(shì)被廣泛使用于地表比輻射率的計(jì)算中。

混合像元，包含不同程度的植被和非植被，像元組成復(fù)雜。 （3.4）

（3.5）

純植被像元，植被覆蓋度為1。

其中dε表征地表幾何形狀，地表相對(duì)平坦地區(qū)取0，對(duì)于粗糙地表計(jì)算方法如下：

（3.6）

3.3 大氣透過(guò)率

大氣透過(guò)率是指電磁波經(jīng)大氣衰減后的電磁輻射通量與入射時(shí)電磁輻射通量的比值。利用遙感技術(shù)估算地表溫度時(shí)，大氣水汽含量是計(jì)算大氣透過(guò)率的關(guān)鍵參數(shù)。通常利用大氣模型軟件模擬大氣透過(guò)率與水汽含量的關(guān)系。使用現(xiàn)有大氣水汽含量產(chǎn)品或氣象站氣象數(shù)據(jù)，通過(guò)大氣模擬來(lái)進(jìn)一步確定大氣透過(guò)率。Landsat 8 ITRS 10波段 0.5～3 g/cm2水汽范圍內(nèi)大氣透過(guò)率與水汽含量的關(guān)系如下：

對(duì)于Landsat 8數(shù)據(jù)來(lái)說(shuō)，大氣作用指數(shù)算法[7]如下：

其中：ω為大氣水汽含量。

3.4 大氣平均作用溫度

大氣剖面氣溫分布和大氣狀態(tài)是決定大氣平均作用溫度的主要因素，由于衛(wèi)星經(jīng)過(guò)研究區(qū)上空的時(shí)間很短，很難實(shí)時(shí)觀測(cè)到大氣剖面數(shù)據(jù)和大氣狀態(tài)。可以通過(guò)近地表近似推算得到大氣平均作用溫度，方便后續(xù)推算地表溫度。

4結(jié)論與展望

利用Landsat系列數(shù)據(jù)進(jìn)行地表溫度反演研究已經(jīng)取得了顯著成效，目前已成為反演地表溫度的主流。針對(duì)不同的熱紅外數(shù)據(jù)，已有很多反演地表溫度的算法，這些算法精度和適用性也都具有差異，今后在地表溫度反演研究中，提高算法反演精度，減少算法參數(shù)，使得模型簡(jiǎn)單易行，是今后研究的方向。

參考文獻(xiàn)

[1]宋挺，嚴(yán)飛，黃君，等.基于環(huán)境小衛(wèi)星反演地表溫度以及與實(shí)測(cè)數(shù)據(jù)的比較與分析[J].環(huán)境科技，2015，28（01）：39-44.

[2]覃志豪，Zhang Minghua，Arnon Karnieli，等.用陸地衛(wèi)星TM 6數(shù)據(jù)演算地表溫度的單窗算法[J].地理學(xué)報(bào)，2001，56（4）：456-466.

[3]Jiménez-Mu?oz JC，Cristóbal J，Sobrino J A，Sòria G，Ninyerola M and Pons X.2009.Revision of the single-channel algorithm for land surface temperature retrieval from Landsat thermal-infrared data.IEEE Transactions on Geoscience and Remote Sensing，2009，47（1）：339-349.

[4]趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京：科學(xué)出版社，2013.

[5]宋挺，段崢，劉軍志，等.Landsat 8數(shù)據(jù)地表溫度反演算法對(duì)比[J].遙感學(xué)報(bào)，2015，19（03）：451-464.

[6]Rozenstein O.，Qin Z. H.，Derimian Y.，et al. Derivation of land surface temperature for Landsat-8 TIRS using a split window algorithm [J].Sensors，2014，14（4）：5768-5780

[7]Jiménez-Mu?oz J C，Sobrino J A，Skokovic'D，Mattar C and Cristóbal J.2014. Land Surface temperature retrieval methods from landsat-8 thermal infrared sensor data[J].IEEE Geoscience and Remote Sensing Letters，2014，11（10）.

[8]Wan Z.，Li Z.L.A physics-based algorithm for retrieving land-surface emissivity and temperature from EOS/MODIS data. IEEE Trans [J].Geosci.Remote Sens，1997（35）：980–996.

收稿日期：2020-03-09

作者簡(jiǎn)介：黎雅楠（1993-），女，碩士，助理工程師，研究方向?yàn)榄h(huán)境遙感及3S 應(yīng)用。