孫艷玲，郭 鵬，高曉燕

（天津師范大學(xué) 城市與環(huán)境科學(xué)學(xué)院，天津 300387）

基于Landsat TM／ETM＋ 的天津地區(qū)地表溫度時(shí)空分布特征研究

孫艷玲，郭 鵬，高曉燕

（天津師范大學(xué) 城市與環(huán)境科學(xué)學(xué)院，天津 300387）

以天津地區(qū)為研究對(duì)象，利用單窗算法對(duì)1993年、2001年和2009年3景Landsat TM／ETM＋遙感影像數(shù)據(jù)進(jìn)行地表溫度反演，將反演結(jié)果進(jìn)行標(biāo)準(zhǔn)差分類(lèi)后得到地表溫度等級(jí)圖.通過(guò)對(duì)不同時(shí)期地表溫度分布進(jìn)行時(shí)空分析，得到天津地表溫度高溫區(qū)在1993年到2009年期間主要集中在中心城區(qū)、濱海新區(qū)以及兩者之間道路沿線(xiàn)等區(qū)域，其中，中心城區(qū)高溫區(qū)不斷加劇，濱海新區(qū)高溫區(qū)從無(wú)到有、由弱到強(qiáng)，變化非常明顯.分析結(jié)果表明：天津地表溫度的空間分布與城市空間結(jié)構(gòu)基本吻合，變化趨勢(shì)與城市發(fā)展基本一致.

地表溫度；Landsat TM／ETM＋；天津；單窗算法

地表溫度（Land Surface Temperature，LST）是研究地表與大氣之間輻射和能量通量傳輸?shù)闹匾獏?shù).隨著城市規(guī)模不斷擴(kuò)大和城市化進(jìn)程不斷加快，通過(guò)城市地表溫度場(chǎng)了解城市熱環(huán)境及其熱島效應(yīng)的研究成為當(dāng)前城市氣候與環(huán)境研究中的熱點(diǎn)問(wèn)題之一［1－2］.研究城市地表溫度分布可以揭示城市的空間結(jié)構(gòu)和城市規(guī)模的發(fā)展變化，有助于引導(dǎo)城市可持續(xù)發(fā)展，提高人們的居住環(huán)境質(zhì)量［3］.

近年來(lái)，國(guó)內(nèi)外通過(guò)城市地表溫度了解城市熱環(huán)境的研究很多［4－5］，歸納起來(lái)，主要可以分為通過(guò)常規(guī)定點(diǎn)觀(guān)測(cè)和通過(guò)衛(wèi)星遙感反演2種實(shí)現(xiàn)途徑其中定點(diǎn)觀(guān)測(cè)難以滿(mǎn)足大空間范圍地表溫度獲取的需要，隨著遙感技術(shù)的發(fā)展，衛(wèi)星熱紅外傳感器已成為獲取地表溫度空間分布的主要途徑.目前，用于反演地表溫度的數(shù)據(jù)源主要有NOAA／AVHRR的熱通道數(shù)據(jù)、TM／ETM＋的第6波段和MODIS的熱紅外波段數(shù)據(jù)等.相比較而言，陸地衛(wèi)星的遙感圖像數(shù)據(jù)地面分辨率較高，其熱紅外波段可用來(lái)分析地表熱輻射和溫度區(qū)域差異.此外，Landsat陸地資源衛(wèi)星已有20年的歷史，積累了豐富的遙感圖像數(shù)據(jù)資料，為研究地表溫度的時(shí)間動(dòng)態(tài)變化提供了有利條件.在利用Landsat TM／ETM＋ 數(shù)據(jù)反演地表溫度的研究中，除了傳統(tǒng)上采用的大氣校正法之外，覃志豪等［6］提出了無(wú)需額外大氣參數(shù)的單窗算法，Jimёnez－Mu?oz等［7］提出了僅需單個(gè)熱紅外波段數(shù)據(jù)的普適性單通道算法.針對(duì)不同的具體實(shí)例，研究者們分別對(duì)大氣校正法、單窗算法和單通道算法進(jìn)行了比較分析，如黃妙芬等［8］利用北京地區(qū)5景Landsat／TM影像和實(shí)測(cè)探空數(shù)據(jù)，分別運(yùn)用大氣校正法、單窗算法和單通道算法反演了北京城區(qū)地表溫度，結(jié)果表明單窗算法與大氣校正法的結(jié)果具有較高的一致性，對(duì)只有一個(gè)熱紅外通道的Landsat／TM數(shù)據(jù)源利用單窗算法反演地表溫度，精度是可以接受的；白潔等［9］基于TM／ETM＋數(shù)據(jù)利用3種算法反演北京地區(qū)地表溫度，結(jié)果顯示，與地面實(shí)測(cè)溫度數(shù)據(jù)相比，單窗算法的結(jié)果與地面實(shí)測(cè)值一致性最好；杜嘉等［10］利用洪河濕地過(guò)境的Landsat／TM圖像、實(shí)測(cè)地面數(shù)據(jù)和MODIS地表發(fā)射率數(shù)據(jù)，運(yùn)用3種算法估算地表溫度，結(jié)果同樣表明單窗算法與實(shí)測(cè)地面數(shù)據(jù)估算結(jié)果一致.

隨著天津?yàn)I海新區(qū)的快速發(fā)展，天津市區(qū)與濱海新區(qū)的建設(shè)布局不斷完善，同時(shí)也帶動(dòng)了沿線(xiàn)地區(qū)的城市化進(jìn)程，這些均對(duì)天津市的熱環(huán)境造成了一定影響.目前，對(duì)天津地區(qū)城市熱環(huán)境的研究主要基于地面觀(guān)測(cè)數(shù)據(jù)，如郭軍等［11］利用天津市區(qū)、塘沽、寧河和薊縣4個(gè)臺(tái)站1964—2006年的氣象觀(guān)測(cè)資料，分析了天津地區(qū)近40年氣溫變化的特征及其與城市化的關(guān)系，得出天津市年平均氣溫呈顯著升高趨勢(shì)、市區(qū)年平均氣溫變化幅度明顯高于其他各區(qū)的結(jié)論.韓素芹等［12］和程晨等［13］利用 TM／ETM＋熱紅外波段數(shù)據(jù)分析了天津城市熱島效應(yīng)，但主要是直接利用遙感數(shù)據(jù)計(jì)算地表亮溫進(jìn)行分析，而亮溫（輻射溫度）僅是衡量物體溫度的一個(gè)指標(biāo)，不代表物體的真實(shí)溫度，兩者雖然密切相關(guān)，但是差異很大，所以存在一定局限性［14］.隨著熱紅外遙感基礎(chǔ)理論研究的深入，通過(guò)反演地表真實(shí)溫度來(lái)研究城市熱環(huán)境變化成為當(dāng)前研究的主流［15－16］.因此本研究利用Landsat TM／ETM＋遙感數(shù)據(jù)，運(yùn)用地理信息系統(tǒng)（Geographic Information Systems GIS）技術(shù)，針對(duì)天津市的實(shí)際情況，選用單窗算法對(duì)天津市的地表溫度進(jìn)行反演，分析了地表溫度的時(shí)間和空間分布特征，希望可以為天津城市建設(shè)與綠化規(guī)劃提供參考.

1 研究方法與數(shù)據(jù)

1.1 研究區(qū)域

本研究主要針對(duì)天津市快速城市化的地區(qū)進(jìn)行地表溫度反演，即選取天津市中心城區(qū)和濱海新區(qū)作為研究區(qū)域，如圖1所示.

圖1 天津市中心城區(qū)及濱海新區(qū)示意圖Fig.1 Location of central city and Binhai new area of Tianjin

研究區(qū)域中，中心城區(qū)包括市內(nèi)6區(qū)（河北區(qū)、紅橋區(qū)、南開(kāi)區(qū)、和平區(qū)、河?xùn)|區(qū)和河西區(qū)）和近郊4區(qū)（東麗區(qū)、津南區(qū)、西青區(qū)和北辰區(qū)）；濱海新區(qū)包括3個(gè)區(qū)，即塘沽區(qū)、漢沽區(qū)和大港區(qū).濱海新區(qū)位于天津市中心城區(qū)的東部沿海地區(qū)，距離市中心50 km，是繼深圳經(jīng)濟(jì)特區(qū)和上海浦東新區(qū)后，帶動(dòng)我國(guó)區(qū)域經(jīng)濟(jì)增長(zhǎng)的又一重大戰(zhàn)略部署區(qū)域，對(duì)環(huán)渤海經(jīng)濟(jì)圈的快速發(fā)展具有巨大的帶動(dòng)作用.

1.2 數(shù)據(jù)來(lái)源及預(yù)處理

本研究選取1993年6月15日、2001年9月1日和2009年8月30日獲取的軌道號(hào)為122／33的3景Landsat TM／ETM＋遙感第6波段數(shù)據(jù)（數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)）反演地表溫度，其中2001年9月1日的遙感數(shù)據(jù)為ETM＋影像，分辨率為60 m，其余兩景均為T(mén)M影像，分辨率為120 m.這3景遙感影像均包含了除北部薊縣外的天津市大部分地區(qū)，由于本研究的研究區(qū)域主要為天津市區(qū)和濱海新區(qū)，因此數(shù)據(jù)可以滿(mǎn)足研究要求.研究所用其他資料包括天津市1∶50 000地形圖和天津市行政區(qū)劃圖.

首先，在ArcGIS9.3支持下參照地形圖，對(duì)3景遙感影像進(jìn)行幾何校正，將其地理坐標(biāo)統(tǒng)一轉(zhuǎn)換為WGS 1984 UTM Zone 50 N，中央經(jīng)線(xiàn)為117°E，遙感數(shù)據(jù)格式由原始的＊.tiff格式統(tǒng)一轉(zhuǎn)換為＊.grid格式，以便進(jìn)行柵格運(yùn)算.然后，利用ArcGIS9.3中的Extract by mask工具，根據(jù)天津市行政區(qū)劃圖，對(duì)經(jīng)過(guò)校正的3景影像進(jìn)行裁切，提取出研究區(qū)域內(nèi)的遙感影像數(shù)據(jù)，后續(xù)工作將僅對(duì)研究區(qū)域進(jìn)行研究.

1.3 地表溫度反演算法

覃志豪等［6］通過(guò)引進(jìn)大氣平均溫度，提出根據(jù)TM6波段數(shù)據(jù)反演地表溫度的單窗算法：

式（1）中：Ts是地表溫度（K）；a和b為常量，a＝－67.355 351，b＝0.458 606；Ta是大氣平均作用溫度，中緯度地區(qū)夏季的Ta＝16.011 0＋0.926 21T0，T0為近地面氣溫（K）；C＝τε，τ是大氣透射率，ε是地表比輻射率；D＝ （1－τ）［1＋ （1－ε）τ］.

τ可以根據(jù)大氣總水汽含量w（g／cm2）來(lái)估算，夏季的w如果落在0.4～1.6 g／cm2之間，則

如果w落在1.6～3.0 g／cm2之間，則

ε可以通過(guò)植被指數(shù)來(lái)估計(jì)，采用覃志豪等［17提出的算法，首先對(duì)遙感影像進(jìn)行分類(lèi)，將地表分為不同的覆蓋類(lèi)型，然后計(jì)算各地表類(lèi)型的植被指數(shù)和植被覆蓋度，從而估算出地表比輻射率值，具體算法請(qǐng)?jiān)斠?jiàn)參考文獻(xiàn)［17］.

Tsensor是衛(wèi)星傳感器所測(cè)的像元亮度溫度（K）亮溫計(jì)算分為2步：首先將由灰度值表示的TM／ETM＋數(shù)據(jù)轉(zhuǎn)換成相應(yīng)的熱輻射強(qiáng)度，然后根據(jù)熱輻射強(qiáng)度推算出相應(yīng)的地表亮度溫度值.對(duì)于TM／ETM＋數(shù)據(jù)，所接收到的輻射強(qiáng)度與其灰度值存在如下關(guān)系：

式（2）中：L（λ）為傳感器所接收到的輻射強(qiáng)度，單位為 mW·cm－2·sr－1·μm－1；QDN是像元的灰度值Qmax＝255為最大灰度值；Lmin（λ）和Lmax（λ）分別為傳感器所接收到的最大和最小輻射強(qiáng)度，可從Landsat TM／ETM＋ 數(shù)據(jù)頭文件中獲取.

利用熱輻射強(qiáng)度求算亮度溫度值可用如下近似計(jì)算式：

2 結(jié)果與分析

2.1 地表溫度分布結(jié)果

利用單窗算法，在ArcGIS9.3支持下，反演出研究區(qū)域1993年6月15日、2001年9月1日和2009年8月30日的地表溫度，圖2為轉(zhuǎn)為攝氏溫度后的反演結(jié)果.針對(duì)反演得到的地表溫度分布圖（圖2），對(duì)結(jié)果進(jìn)行標(biāo)準(zhǔn)差法（Standard Deviation分類(lèi)，這種分類(lèi)方法首先計(jì)算柵格影像LST的平均值和標(biāo)準(zhǔn)差，然后根據(jù)這2個(gè)參數(shù)通過(guò)一系列計(jì)算將研究區(qū)的LST分為7個(gè)等級(jí)，按照溫度由低到高的順序，將分類(lèi)結(jié)果依次命名為極低溫區(qū)、低溫區(qū)較低溫區(qū)、中溫區(qū)、較高溫區(qū)、高溫區(qū)和極高溫區(qū).由于不同時(shí)期地表溫度分布不同，3個(gè)年份溫度區(qū)分級(jí)對(duì)應(yīng)的地表溫度范圍也不同，結(jié)果如表1所示.

圖2 不同時(shí)期天津地區(qū)地表溫度分布圖Fig.2 Distribution of land surface temperature in three dates of Tianjin area

表1 標(biāo)準(zhǔn)差法分類(lèi)表Tab.1 Standard deviation classification ℃

根據(jù)標(biāo)準(zhǔn)差法分類(lèi)結(jié)果，對(duì)LST影像進(jìn)行重采樣（Reclassify），得到各溫度區(qū)的分布狀況，如圖3所示，并由此分別計(jì)算不同時(shí)期7類(lèi)溫度區(qū)所占的面積比例，統(tǒng)計(jì)結(jié)果如表2所示.

圖3 研究區(qū)不同時(shí)期的各溫度區(qū)分布圖Fig.3 Distribution of temperature classes in three dates of study area

表2 天津地區(qū)3個(gè)時(shí)期的溫度區(qū)面積比例表Tab.2 Proportion of temperature areas in three dates of study area %

2.2 地表溫度空間分布特征

對(duì)比天津城區(qū)不同時(shí)期的地表溫度分布圖可以發(fā)現(xiàn)，隨著天津城市化進(jìn)程的不斷加快和天津?yàn)I海新區(qū)的快速發(fā)展，地表溫度空間分布發(fā)生了明顯變化，本研究對(duì)研究區(qū)域不同時(shí)期的地表溫度區(qū)分布圖（圖3）和不同溫度區(qū)面積比例（表2）進(jìn)行分析.

從1993年的地表溫度分布圖可以看出，研究區(qū)域中代表高溫的高溫區(qū)和極高溫區(qū)主要分布在天津市中心城區(qū)海河兩側(cè)的沿河一帶，且主要集中在紅橋區(qū)、河北區(qū)、南開(kāi)區(qū)、和平區(qū)和河?xùn)|區(qū)相接壤的地區(qū)，這些區(qū)域主要為天津的老城區(qū)，在20世紀(jì)90年代人口密度較大.塘沽區(qū)沿海河下游處LST略高，但是面積較小，而且比中心城區(qū)溫度要低.同時(shí)，從圖3中可以看出，代表低溫的低溫區(qū)和極低溫區(qū)主要分布在漢沽區(qū)和塘沽區(qū).研究區(qū)域的大部分地區(qū)在1993年屬于較低溫區(qū)、中溫區(qū)和較高溫區(qū)，這3類(lèi)區(qū)域的面積達(dá)到總面積的85.06%，代表高溫的高溫區(qū)和極高溫區(qū)的面積比例只有6.69%，代表低溫的低溫區(qū)和極低溫區(qū)比例為8.24%.此外，從1993年的地表溫度分布中還可以發(fā)現(xiàn)，中心城區(qū)的西北部（主要是北辰區(qū)范圍內(nèi)）和濱海新區(qū)的大港區(qū)范圍內(nèi)地表溫度也較高，這兩地在1993年屬于郊區(qū)，土地類(lèi)型以農(nóng)業(yè)為主，地表溫度本應(yīng)較低，但反演溫度卻與市內(nèi)六區(qū)相接近，這是因?yàn)?月份正值天津冬小麥?zhǔn)斋@季節(jié)，地表植被覆蓋條件較差，大部分地區(qū)在這個(gè)時(shí)期主要為裸地，而且北辰區(qū)、西青區(qū)和大港區(qū)是天津主要的冬小麥種植區(qū)域，因此，1993年6月15日反演所得的這些區(qū)域的地表溫度也較高.

從2001年的地表溫度分布圖可以看出，研究區(qū)的高溫區(qū)主要分布在中心城區(qū)的市內(nèi)六區(qū)和濱海新區(qū)的塘沽一帶.其中，市內(nèi)六區(qū)極高溫區(qū)的范圍較1993年明顯變大，且呈輻射狀向周?chē)l(fā)展.同時(shí)，從圖3中可以發(fā)現(xiàn)，在市內(nèi)六區(qū)中河西區(qū)和南開(kāi)區(qū)范圍內(nèi)存在非常明顯的低溫區(qū)域，這些區(qū)域主要是公園綠地，如水上公園的水域面積較大，植被覆蓋率高，因此成為中心城區(qū)中的低溫區(qū).濱海新區(qū)中塘沽區(qū)、大港城區(qū)和漢沽城區(qū)的地表溫度均呈升高趨勢(shì)溫度區(qū)分布圖顯示塘沽區(qū)中的極高溫區(qū)主要分布在海河入海口處，且塘沽區(qū)高溫區(qū)的面積大于大港區(qū)和漢沽區(qū).與1993年相比，漢沽區(qū)也出現(xiàn)了明顯的高溫區(qū).此外，在市內(nèi)六區(qū)與塘沽區(qū)之間出現(xiàn)了帶狀高溫帶.統(tǒng)計(jì)結(jié)果顯示，2001年研究區(qū)域內(nèi)代表高溫的高溫區(qū)和極高溫區(qū)的面積比例為8.49%，較1993年有所增加.

2009年的地表溫度分布圖顯示，研究區(qū)域內(nèi)高溫區(qū)的分布范圍進(jìn)一步擴(kuò)大.除了水體覆蓋和密度較高的綠地外，市內(nèi)六區(qū)基本都屬于高溫區(qū)和極高溫區(qū)，且向四周擴(kuò)散的范圍和強(qiáng)度較大，中心城區(qū)中4個(gè)近郊區(qū)高溫區(qū)和極高溫區(qū)的面積有所增大，濱海新區(qū)高溫區(qū)的范圍和強(qiáng)度也都明顯變大，特別是塘沽區(qū)內(nèi)，高溫區(qū)的范圍明顯擴(kuò)大.同時(shí)，塘沽區(qū)與中心城區(qū)間的高溫帶也更加顯著，范圍逐漸變寬，這表明隨著中心城區(qū)和濱海新區(qū)之間各種交通設(shè)施的快速發(fā)展，沿途地區(qū)的城市化過(guò)程加快，城鎮(zhèn)建設(shè)用地面積增加，導(dǎo)致地表溫度逐漸增加.在高溫區(qū)面積不斷增加的同時(shí)，研究區(qū)域中低溫區(qū)和極低溫區(qū)的面積進(jìn)一步減少.統(tǒng)計(jì)結(jié)果顯示，2009年研究區(qū)域內(nèi)高溫區(qū)和極高溫區(qū)的面積比例為11.38%，低溫區(qū)和極低溫區(qū)的比例為1.58%.

2.3 地表溫度時(shí)間分布特征

為了進(jìn)一步分析天津城區(qū)和濱海新區(qū)1993年至2009年地表溫度的變化規(guī)律，分別對(duì)中心城區(qū)塘沽區(qū)、大港區(qū)和漢沽區(qū)的地表溫度等級(jí)進(jìn)行了面積統(tǒng)計(jì)，結(jié)果如表3所示.

表3 不同地區(qū)不同時(shí)期各溫度區(qū)面積比例Tab.3 Proportion of temperature areas in three dates of different areas %

由表3可以看出，中心城區(qū)中低溫區(qū)和極低溫區(qū)的面積比例由1993年的1.86%下降到2001年的1.41%，2009年則下降到0.8%，可見(jiàn)中心城區(qū)低溫范圍在這17年中呈減小趨勢(shì)，并且2001年至2009年間的減小趨勢(shì)更為顯著.同時(shí)，高溫區(qū)和極高溫區(qū)的面積比例在1993年為9.04%，2001年增加達(dá)到10.18%，之后仍然逐年增加，在2009年達(dá)到17.39%，增幅十分明顯.其中極高溫區(qū)從無(wú)到有，面積從小到大，變化也十分明顯，這表明天津中心城區(qū)的地表溫度在不斷增加.

濱海新區(qū)中，塘沽區(qū)地表溫度增加顯著.由表3可知，1993年極低溫區(qū)和低溫區(qū)的面積比例為26.87%，極高溫區(qū)的面積比例為0.00%，高溫區(qū)的面積比例也僅為0.06%；2001年該區(qū)極低溫區(qū)和低溫區(qū)的面積比例有所下降，同時(shí)高溫區(qū)和極高溫區(qū)的面積大幅增加，面積比例達(dá)到9.19%；2009年塘沽區(qū)高溫區(qū)和極高溫區(qū)的面積比例繼續(xù)增加，同時(shí)低溫區(qū)和極低溫區(qū)的面積比例迅速下降.漢沽區(qū)從1993年到2009年也呈現(xiàn)出高溫區(qū)面積比例逐漸增加，低溫區(qū)面積比例逐漸減小的趨勢(shì)，同時(shí)還表現(xiàn)出極高溫區(qū)從無(wú)到有的特征.大港區(qū)相對(duì)于其他區(qū)域，情況有所不同，如1993年6月15日的遙感影像中，LST分布異常，這對(duì)該區(qū)1993年的各溫度區(qū)分布比例具有直接影響.

3 結(jié)論

本研究根據(jù)1993年到2009年的遙感影像數(shù)據(jù)，利用單窗算法反演出不同時(shí)期天津地區(qū)的地表溫度分布圖，并借助ArcGIS的空間分析功能對(duì)研究結(jié)果進(jìn)行標(biāo)準(zhǔn)差分類(lèi)，計(jì)算出不同時(shí)期各溫度區(qū)所占面積的百分比，從而對(duì)LST的空間分布及其變化趨勢(shì)進(jìn)行分析.研究結(jié)果表明：1993年到2009年期間，地表溫度相對(duì)較高的區(qū)域由以市內(nèi)六區(qū)為中心逐漸發(fā)展成以市內(nèi)六區(qū)及塘沽區(qū)為中心.這種分布格局的變化同時(shí)也帶動(dòng)了市內(nèi)六區(qū)和塘沽區(qū)之間地區(qū)的地表溫度變化，且影響范圍逐漸擴(kuò)大，強(qiáng)度逐漸增強(qiáng).雖然受天津市地理位置、氣候條件和農(nóng)作物成熟收割等因素的影響，一些地區(qū)存在地表溫度異常區(qū)，但通過(guò)3年的地表溫度等級(jí)圖和溫度等級(jí)統(tǒng)計(jì)表可以看出，天津中心城區(qū)和濱海新區(qū)高溫區(qū)強(qiáng)度從弱到強(qiáng)，范圍從小到大.地表溫度分布格局的變化反映了天津市空間布局的調(diào)整，即由最初的“一條扁擔(dān)挑兩頭”到以中心城區(qū)和濱海新區(qū)核心區(qū)為主副中心構(gòu)建的“一主一副”雙中心布局，再到“雙港雙城”的中心城市布局結(jié)構(gòu)［18］.

［1］ 趙靜，宮阿都，徐捷.國(guó)內(nèi)外城市熱島研究進(jìn)展［J］.自然災(zāi)害學(xué)報(bào)，2008，17（6）：64－67.

［2］ 周志民.城市熱島遙感研究進(jìn)展［J］.安徽農(nóng)業(yè)科學(xué)，2011，39（1）：606－607，614.

［3］ 陳云浩，李京，李曉兵.城市空間熱環(huán)境遙感分析—格局、過(guò)程、模擬與影響［M］.北京：科學(xué)出版社，2003.

［4］ 段金龍，宋軒，張學(xué)雷.基于RS的鄭州市城市熱島效應(yīng)時(shí)空演變［J］.應(yīng)用生態(tài)學(xué)報(bào)，2011，22（1）：165－170.

［5］ 董妍，李星敏，楊艷超，等.西安城市熱島的時(shí)空分布特征［J］干旱區(qū)資源與環(huán)境，2011，25（8）：107－112.

［6］ 覃志豪，ZHANG M H，ARNON K，等.用陸地衛(wèi)星TM6數(shù)據(jù)演算地表溫度的單窗算法［J］.地理學(xué)報(bào)，2001，56（4）：456－466

［7］ JIM?NEZ－MU?OZ J C，SOBRINO J A.A generalized single channel method for retrieving land surface temperature from remote sensing data［J］.Journal of Geophysical Research 2003，108（22）：4688－4695.

［8］ 黃妙芬，邢旭峰，王培娟，等.利用LANDSAT／TM熱紅外通道反演地表溫度的三種方法比較［J］.干旱區(qū)地理，2006，29（1）：132－137.

［9］ 白潔，劉紹民，扈光.針對(duì)TM／ETM＋遙感數(shù)據(jù)的地表溫度反演與驗(yàn)證［J］.農(nóng)業(yè)工程學(xué)報(bào)，2008，24（9）：148－154.

［10］ 杜嘉，張柏，宋開(kāi)山，等.基于Landsat－5 TM的洪河濕地地表溫度估算方法對(duì)比研究［J］.遙感技術(shù)與應(yīng)用，2009，24（3）：312－319.

［11］ 郭軍，李明財(cái)，劉德義.近40年來(lái)城市化對(duì)天津地區(qū)氣溫的影響［J］.生態(tài)環(huán)境學(xué)報(bào)，2009，18（1）：29－34.

［12］ 韓素芹，郭軍，黃歲樑，等.天津城市熱島效應(yīng)演變特征研究［J］.生態(tài)環(huán)境，2007，16（2）：280－284.

［13］ 程晨，蔡喆，閆維，等.基于Landsat TM／ETM＋的天津城區(qū)及濱海新區(qū)熱島效應(yīng)時(shí)空變化研究［J］.自然資源學(xué)報(bào)2010，25（10）：1727－1737.

［14］ 申雙和，趙小艷，楊沈斌，等.利用ASTER數(shù)據(jù)分析南京城市地表溫度分布［J］.應(yīng)用氣象學(xué)報(bào)，2009，20（4）：458－464

［15］ 岳文澤，徐麗華，徐建華.20世紀(jì)90年代上海市熱環(huán)境變化及其社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)力［J］.生態(tài)學(xué)報(bào)，2010，30（1）：155－164

［16］ 梁益同，陳正洪，夏智宏.基于RS和GIS的武漢城市熱島效應(yīng)年代演變及其機(jī)理分析［J］.長(zhǎng)江流域資源與環(huán)境2010，19（8）：914－918.

［17］ 覃志豪，李文娟，徐斌，等.陸地衛(wèi)星TM6波段范圍內(nèi)地表比輻射率的估計(jì)［J］.國(guó)土資源遙感，2004（3）：28－42.

［18］ 尹海林.國(guó)家戰(zhàn)略背景下的天津市空間發(fā)展戰(zhàn)略規(guī)劃［J］.時(shí)代建筑，2010（5）：6－9.

Study on temporal and spatial variation of land surface temperature based on Landsat TM／ETM＋in Tianjin

SUNYan－ling，GUOPeng，GAOXiao－yan

（College of Urban and Environmental Science，Tianjin Normal University，Tianjin 300387，China）

The spatial and temporal variation of land surface temperature（LST）in central city and Binhai new area of Tianjin was studied with Landsat TM／ETM＋.Using mono－window algorithm，the LST was retrieved from Landsat TM／ETM＋thermal band data from three dates in 1993，2001 and 2009.Then the LST was classified into distribution maps by means of standard deviation，and the proportion of temperature class was also calculated and plotted.The results of analysis of maps and plots showed that high temperature area mainly focused in central city，Binhai new area and areas along the roads between the two areas.From 1993 to 2009，high temperature area was becoming larger and harder.The effect of high LST on central city was aggravated，and the effect of high LST on Binhai new area started from scratch，from less to more and experienced obvious changes.The results showed that in some kind of degree，the change of LST distribution reflected the adjustment of urban layout structure made by government.

land surface temperature；Landsat TM／ETM＋；Tianjin；mono－window algorithm

P463.3

A

1671－1114（2012）01－0048－06

2011－10－08

國(guó)家自然科學(xué)基金資助項(xiàng)目（41001022）；中國(guó)科學(xué)院知識(shí)創(chuàng)新工程重要方向資助項(xiàng)目（KZCX2－EW－202）

孫艷玲（1977—），女，講師，主要從事全球變化、資源環(huán)境和3S應(yīng)用方面的研究.

（責(zé)任編校 紀(jì)翠榮）