鄔 丹，鄒 凱，楊 波

（1.湖南師范大學(xué) 資源與環(huán)境科學(xué)學(xué)院，湖南 長(zhǎng)沙 410081；2.益陽(yáng)市桃江縣國(guó)土資源局，湖南 益陽(yáng) 413400）

基于RS與GIS的長(zhǎng)沙市城市熱島影響因子分析

鄔 丹1，鄒 凱2，楊 波1

（1.湖南師范大學(xué) 資源與環(huán)境科學(xué)學(xué)院，湖南 長(zhǎng)沙 410081；2.益陽(yáng)市桃江縣國(guó)土資源局，湖南 益陽(yáng) 413400）

基于RS和GIS技術(shù)，以Landsat-7 ETM+為數(shù)據(jù)源,綜合利用遙感熱紅外影像進(jìn)行地表溫度反演與熱島強(qiáng)度提取。結(jié)果表明，長(zhǎng)沙下墊面熱島效應(yīng)的空間分布與延展和城市建成區(qū)的外圍基本相符，人口密集、建筑物較多之地?zé)釐u效應(yīng)明顯，河流、湖泊及公園等地?zé)釐u相對(duì)較弱。結(jié)合緩沖區(qū)分析、疊置分析等技術(shù)，提取熱島最強(qiáng)與最弱區(qū)域周圍500 m范圍內(nèi)的城市用地類型，分析其相關(guān)關(guān)系，得出的結(jié)論是:水體與綠地具有緩解城市熱島的作用，城市建設(shè)用地與裸地加劇城市熱島效應(yīng)。

長(zhǎng)沙市；RS；GIS；熱島效應(yīng)；影響因子

1 研究區(qū)及數(shù)據(jù)處理

長(zhǎng)沙市位于北緯27°51' ～28°40'，東經(jīng)111°53'～114°5'。全市土地總面積1.18萬(wàn)km2，其中市區(qū)面積1 938 km2。到2010年，長(zhǎng)沙市建成區(qū)面積已達(dá)320 km2，全市城市化率為67.7%[1-9]。

本文采用2006-11-01長(zhǎng)沙市Landsat-7 ETM+數(shù)據(jù)（軌道號(hào)：123/41），圖像質(zhì)量好，清晰無(wú)云，并經(jīng)過(guò)幾何精校正（RMS＜1）。

圖1 長(zhǎng)沙市土地利用分類結(jié)果圖

參考中國(guó)科學(xué)院的中國(guó)土地資源分類系統(tǒng)[10]，結(jié)合長(zhǎng)沙市土地利用的特點(diǎn)，將其分為5大類(如圖1）：林地，城市建設(shè)用地，水體，裸地，耕地。利用最大似然比分類法進(jìn)行土地利用分類，精度達(dá)到本次研究的要求，如表1。

2 研究方法

2.1 地表溫度反演與熱島強(qiáng)度提取

首先，將像元灰度值DN轉(zhuǎn)化為輻射亮度值：

表1 長(zhǎng)沙市2006年土地利用分類誤差矩陣/像元數(shù)/個(gè)

式中，La為輻射亮度，單位是該波段灰度值分別為255和1時(shí)所對(duì)應(yīng)的光譜輻射率。第六波段光譜輻射率低增益和高增益分別是單位將其代入式（1），計(jì)算出輻射亮度La：

式中，La為輻射亮度；T1為亮度溫度。對(duì)于ETM+數(shù)據(jù)由于亮度溫度和真實(shí)地表溫度差異較大（＞5 ℃）[11]，應(yīng)將T1轉(zhuǎn)化為真實(shí)地表溫度：

式中，T2為真實(shí)地表溫度；ε為地表發(fā)射率。根據(jù)不同地物比輻射率值的差異，利用NDVI[12]和MNDWI[13,14]（改進(jìn)的歸一化水體指數(shù)）構(gòu)造決策分類樹(shù)，將土地利用類型分為植被、水體、裸地和建設(shè)用地（圖2）。其中水體比輻率為0.995[15]，植被比輻射率為0.986[15]，農(nóng)閑時(shí)期的耕地沒(méi)有植被覆蓋，其比輻射率與裸地和建筑用地相似，近似看作0.970[16-18]。通過(guò)亮溫-地表溫度轉(zhuǎn)換，溫度反演精度在2 ℃左右[15]。

圖2 決策樹(shù)分類

熱島強(qiáng)度反映的是鄉(xiāng)村與城市之間溫度的相對(duì)概念，以平均溫度為因子能較好地反映熱島強(qiáng)度[19]：

式中，UHI為熱島強(qiáng)度；T2為地表溫度；T為平均溫度。

2.2 緩沖區(qū)分析與疊置分析

緩沖區(qū)分析即給定一個(gè)空間對(duì)象或集合后確定它們的鄰域,鄰域的大小由鄰域半徑R決定，d為歐氏距離或其他距離。對(duì)于給定對(duì)象A[20]：

疊置分析是將代表不同主題的數(shù)據(jù)層面進(jìn)行疊置，產(chǎn)生的新數(shù)據(jù)層綜合了原來(lái)2個(gè)或多個(gè)層面要素所具有的屬性[20]。

3 結(jié)果分析

3.1 長(zhǎng)沙市熱島情況分析

依據(jù)式（1）～（5）得到長(zhǎng)沙市溫度反演圖（圖3），最高溫度22.9 ℃，最低溫16.2 ℃，平均溫度18.5 ℃。據(jù)式（6）得到長(zhǎng)沙市熱島強(qiáng)度，將其均分為10個(gè)等級(jí)（圖 4），1～10代表熱島強(qiáng)度逐漸增大。

1）長(zhǎng)沙下墊面熱島效應(yīng)的空間分布和城市建成區(qū)的外圍基本相符。

2）城市熱島的熱場(chǎng)強(qiáng)度不均勻,其峰值往往是具有高強(qiáng)度熱源的工業(yè)區(qū)或大型公共場(chǎng)所，如長(zhǎng)沙高橋大市場(chǎng)、長(zhǎng)沙火車站、汽車站、五一廣場(chǎng)商業(yè)區(qū)等人口數(shù)量大且密集區(qū)域，熱輻射較大，熱島效應(yīng)較明顯。

3）湘江上的沙洲以及江邊的沙灘由于砂石裸露，在太陽(yáng)輻射下升溫快，在反演圖上亮度值明顯較高。此外像黎托與繞城高速和長(zhǎng)沙大道相鄰，發(fā)展建設(shè)過(guò)程中裸地面積廣，熱島效應(yīng)也較明顯。

4）岳麓山、南郊公園、烈士公園各類綠化較好地的亮度值相對(duì)較低，熱輻射較小，在圖4中基本屬于1、2級(jí)，城市熱島效應(yīng)不明顯。

5）由圖4可看出，水域周圍一定范圍內(nèi)溫度相對(duì)偏低，如湘江、望月湖、烈士公園的年嘉湖等都有此現(xiàn)象，其周圍地段熱島效應(yīng)不明顯，屬于熱輻射較小區(qū)。

圖3 長(zhǎng)沙市LS-7ETM+地表溫度反演結(jié)果圖

圖4 長(zhǎng)沙市熱島強(qiáng)度分級(jí)圖

3.2 影響因子分析

太陽(yáng)對(duì)城市下墊面的輻射強(qiáng)度是均勻的，但因?yàn)橄聣|面的結(jié)構(gòu)和熱力學(xué)性質(zhì)不同而導(dǎo)致了區(qū)域亮溫的差異[21,22]。因此下墊面性質(zhì)差異是影響城市熱島的最主要因素。根據(jù)長(zhǎng)沙市實(shí)際情況，分別對(duì)長(zhǎng)沙城市熱島最強(qiáng)/最弱區(qū)（圖4中的第10級(jí)和第1級(jí)）建立半徑為500 m的緩沖區(qū)，與長(zhǎng)沙市土地利用圖進(jìn)行疊置。分別對(duì)熱島最強(qiáng)區(qū)域與最弱區(qū)域進(jìn)行緩沖分析與疊置分析（表1、表2），得出如下結(jié)論：

城市建設(shè)用地：圖5a中城市建設(shè)用地遠(yuǎn)遠(yuǎn)大于圖5b, 其所占總用地比例高達(dá)43.2%。說(shuō)明城市建設(shè)用地加劇了城市熱島。長(zhǎng)沙市城中高大的建筑物、磚石、水泥路面吸收太陽(yáng)輻射強(qiáng)，同時(shí)由于城市人車數(shù)量多，開(kāi)發(fā)建設(shè)較密集，城市上空有大氣污染覆蓋層吸收地面長(zhǎng)波輻射,地面不易冷卻，熱島效應(yīng)明顯。

裸地：長(zhǎng)沙市的裸地一般是指枯水期的河灘和裸露的地面，圖5a中裸地占了緩沖區(qū)內(nèi)總用地的11.97%。在道路寬且多的地區(qū)熱傳導(dǎo)與對(duì)流效應(yīng)明顯，由于白天吸收熱量多，晚上放熱快，使得近地面空氣被加熱，為城市熱島提供了基礎(chǔ)。

林地：空間布局上形成一定規(guī)模的綠地能更好地緩解城市熱島,即綠地的規(guī)模效應(yīng)[23]。圖5b中林地所占比例高達(dá)33.53%，谷山、岳麓山等聚集的樹(shù)林對(duì)城市熱島效應(yīng)的緩解作用好,其周圍氣溫低于城區(qū)，說(shuō)明林地對(duì)城市熱島效應(yīng)具有緩解作用。

耕地：耕地一般屬于綠地，但根據(jù)實(shí)驗(yàn)結(jié)果圖5a耕地所占比例略高于圖5b。由于實(shí)驗(yàn)采用的是11月份的數(shù)據(jù)，屬于農(nóng)閑期，收割后的耕地基本沒(méi)有植被覆蓋，使得其熱力性質(zhì)接近于裸地，這個(gè)時(shí)期的耕地對(duì)于緩解城市熱島效應(yīng)作用不明顯。

水體：由于長(zhǎng)沙市內(nèi)水域面積小，圖5a與圖5b中水體比例相差不大，在圖5b中可明顯看出水體周圍的城市用地溫度相對(duì)較低，城區(qū)內(nèi)的多個(gè)低溫中心主要分布在城市中幾個(gè)較大的湖區(qū)及其周邊地區(qū)，如咸嘉湖、年嘉湖等。水體可在局部區(qū)域中改善氣候狀況，緩解城市熱島。

此外，城市熱島效應(yīng)是多種因素共同作用的結(jié)果，根據(jù)城市盛行風(fēng)向以及城市地形地勢(shì)，合理設(shè)置城市街道可以形成良好的通風(fēng)走廊，有利于熱量的散失；城市工業(yè)園區(qū)的污染物、煙塵排放量多，合理布局工業(yè)園區(qū)，控制大氣污染，也能緩解城市熱島效應(yīng)。

圖5 長(zhǎng)沙市UHI最強(qiáng)與最弱緩沖區(qū)內(nèi)土地利用類型圖

表2 長(zhǎng)沙市熱島最強(qiáng)緩沖區(qū)土地利用類型像元統(tǒng)計(jì)

表3 長(zhǎng)沙市熱島最弱緩沖區(qū)土地利用類型像元統(tǒng)計(jì)

4 結(jié) 語(yǔ)

本文基于LS-7ETM+數(shù)據(jù)，采用遙感與GIS技術(shù)，對(duì)長(zhǎng)沙市熱島分布情況和影響因子進(jìn)行分析，結(jié)論如下：

1）長(zhǎng)沙市熱島分布情況和城市建成區(qū)較為一致且與其下墊面的不同用地類型有很大關(guān)聯(lián)，建設(shè)用地及裸地的熱島效應(yīng)強(qiáng)，水體與綠地對(duì)城市熱島有緩解作用。

2）利用地表參數(shù)反演能有效地研究城市熱島效應(yīng)。通過(guò)緩沖區(qū)與疊置分析能較準(zhǔn)確直觀地反映城市熱島的影響因子。

日后研究可從以下幾個(gè)方面展開(kāi)：

1）對(duì)城市熱島的研究是根據(jù)地面溫度來(lái)分析的，在有風(fēng)情況下，地表溫度與大氣溫度有一定差別，可以嘗試尋求更好的解決方法。

2）影響長(zhǎng)沙市熱島的因子很多，且影響過(guò)程較為復(fù)雜，日后可多角度深入探討，以求更準(zhǔn)確地進(jìn)行描述。

[1] 陳云浩,李京,李曉兵.城市空間熱環(huán)境[M].北京:科學(xué)出版社,2003

[2] 李志乾,鞏彩蘭,胡勇,等.城市熱島遙感研究進(jìn)展[J].遙感信息,2009(8):100-105

[3] 周紅妹,周成虎,葛偉強(qiáng),等.基于遙感和GIS的城市熱場(chǎng)分布規(guī)律研究[J].地理學(xué)報(bào),2001,56(2):189-196

[4] 楊波,劉宇鵬,馬蘇，等.基于改進(jìn)TAE算法的城市熱島信息提取技術(shù)[J].地理研究,2011,30(5):2-9

[5] 張佳華,侯英雨,李貴才,等.北京城市及周邊熱島日變化及季節(jié)特征的衛(wèi)星遙感研究與影響因子分析[J].中國(guó)科學(xué)D輯: 地球科學(xué),2005, 35 (增刊): 187-194

[6] 蘇偉忠,楊英寶,楊桂山.南京市熱場(chǎng)分布特征及其與土地利用/覆被關(guān)系研究[J].地理科學(xué)，2005,25(6):697-703

[7] 丁金才,周紅妹,葉其欣.從上海市熱島演變看城市綠化的重要意義[J].氣象,2002,28(2): 22-24

[8] 卓嘎,普布次仁,劉偉東.西藏拉薩市熱島效應(yīng)及其影響因子分析[J].中國(guó)環(huán)境監(jiān)測(cè), 2012(26):61-67

[9] 長(zhǎng)沙市政府.長(zhǎng)沙行政區(qū)劃[EB/OL].http://www.xzqh.org/ html/list/217.html,2012-02-28

[10] 李孝芳,中國(guó)土地資源分類原則和系統(tǒng)的探討[J].自然資源,1981(2):1-6

[11] 覃志豪,Li W J, Zhang M H,等.單窗算法的大氣參數(shù)估計(jì)方法[J].國(guó)土資源遙感,2003(56):37-43

[12] 張?jiān)聟?趙志強(qiáng),李雙成,等.基于SPOT NDVI 的華北北部地表植被覆蓋變化趨勢(shì)[J].地理研究,2008,27(4): 746-754

[13] 徐涵秋.利用改進(jìn)的歸一化差異水體指數(shù)(MNDWI)提取水體信息的研究[J].遙感學(xué)報(bào),2005,9(5):589-595

[14] 徐涵秋.基于譜間特征和歸一化指數(shù)分析的城市建筑用地信息提取[J].地理研究,2005,24(2):311-320

[15] 覃志豪,李文娟,徐斌,等.陸地衛(wèi)星TM6波段范圍內(nèi)地表比輻射率的估計(jì)[J].國(guó)土資源遙感,2004(3):28-32

[16] Nerry F, Labed J, Stoll M P. Spectral Properties of Land Surfaces in the Thermal Infrared. 1.Laboratory Measurements of Absolute Spectral Emissivity Signatures[J]. Journal of Geophysical Research,1990(95):7 027-7 044

[17] Salisbury J W, D'aria D M. Emissivity of Terrestrial Materials in the 8-14 μm Atmospheric Window [J]. Remote Sensing of Environment, 1992(42): 83-106.

[18] Zhang J. Spatial and Statistical Analysis of Thermal Satellite Imagery for Extraction of Coal Fire Related Anomalies [D]. Austr：Technical University Vienna, 2004

[19] 劉宇鵬,楊波,陳崇.基于遙感的長(zhǎng)沙市城市熱島效應(yīng)時(shí)空分析[J].遙感信息,2011,15(6): 73-78

[20] 湯國(guó)安,楊昕.地理信息系統(tǒng)[M]. 北京:科學(xué)出版設(shè)，2011

[21] 陳云浩,王潔,李曉兵.夏季城市熱場(chǎng)的衛(wèi)星遙感分析[J].國(guó)土資源遙感,2002(4) :55-59

[22] 劉宇鵬,楊波,陳崇.U-TAE算法提取城市熱島信息的空間尺度分析[J].地理科學(xué)進(jìn)展,2011,30(9) :1 118-1 124

[23] 王勇,李發(fā)斌,李何超.RS與GIS支持下城市熱島效應(yīng)與綠地空間相關(guān)性研究[J].環(huán)境科學(xué)研究，2008,21(4):81-87

P237.9

B

1672-4623（2015）01-0065-04

10.3969/j.issn.1672-4623.2015.01.022

鄔丹，碩士，研究方向?yàn)檫b感與GIS應(yīng)用。

2014-01-21。

項(xiàng)目來(lái)源：國(guó)家自然科學(xué)基金面上資助項(xiàng)目（41171342）；省高校創(chuàng)新平臺(tái)開(kāi)放基金資助項(xiàng)目（10K042）；湖南省重點(diǎn)學(xué)科資助項(xiàng)目。