楊松勇， 劉士彪， 劉昭華

（江西理工大學(xué)建筑與測(cè)繪工程學(xué)院，江西 贛州341000）

0 引 言

城市化進(jìn)程的加快，不透水面為主的人工景觀逐漸取代城市周邊的植被、水體等自然景觀格局，在人工下墊面比例提高的同時(shí)，由于地表覆被的變化導(dǎo)致了地表熱輻射、傳遞與存儲(chǔ)的一系列改變，最終以城市熱島的形式表現(xiàn)出來(lái)[1].城市化進(jìn)程的逐漸深入，城市熱島效應(yīng)愈演愈烈，導(dǎo)致人類聚居地的氣候環(huán)境、生態(tài)環(huán)境等遭到嚴(yán)重破壞，能源、健康乃至安全問(wèn)題隨之而來(lái)，對(duì)人類社會(huì)的可持續(xù)發(fā)展構(gòu)成嚴(yán)重威脅.

隨著空間信息技術(shù)的發(fā)展，高分辨率的衛(wèi)星熱紅外遙感技術(shù)逐步發(fā)展并趨于完善，利用熱紅外數(shù)據(jù)獲取地表信息、監(jiān)測(cè)城市熱場(chǎng)已成為現(xiàn)實(shí)，且該技術(shù)手段在城市熱島研究中發(fā)揮著越來(lái)越重要的作用[2].目前應(yīng)用遙感技術(shù)在城市熱島研究方面主要包括地表溫度解譯方法的研究[3-6]、城市熱島空間格局與演化研究[7-9]以及城市熱島形成與驅(qū)動(dòng)力研究[10-11]，除此之外，還有城市熱環(huán)境與地表特征之間的關(guān)系研究[12]，主要包括城市熱島與土地利用變化[13-14]以及與不透水面之間[15-16]的關(guān)系研究.例如樊亞鵬等[17]利用四期Landsat影像對(duì)廣州市的城市擴(kuò)展及熱環(huán)境效應(yīng)進(jìn)行研究，并討論了城市熱環(huán)境與建筑指數(shù)和植被指數(shù)之間的定量關(guān)系.Mathew A等[18]利用MODIS數(shù)據(jù)對(duì)印度的艾哈邁達(dá)巴德和齋浦爾的城市熱島效應(yīng)進(jìn)行了研究，得出雖白天觀測(cè)到的溫度較高，但城市熱島現(xiàn)象在夜間更明顯這一結(jié)論，并在齋浦爾進(jìn)行了24小時(shí)的溫度監(jiān)測(cè)，進(jìn)行了地表溫度測(cè)量研究，結(jié)果表明土壤具有快速升溫的特性，植被覆蓋區(qū)在整個(gè)觀測(cè)期間的溫度最低.

然而絕大多數(shù)研究主要集中于大型城市的熱島效應(yīng)，對(duì)于中小城市的熱島變遷過(guò)程及格局研究較少.隨著城市化進(jìn)程的深入，中小城市擴(kuò)張明顯，地類演化劇烈，研究由此引起的熱場(chǎng)分布變化無(wú)疑對(duì)城市的規(guī)劃發(fā)展具有一定的指導(dǎo)作用，研究城市熱島效應(yīng)與下墊面各參數(shù)之間的關(guān)系對(duì)提出減緩熱島效應(yīng)非常具有價(jià)值.文中以2009年、2014年、2017年Landsat影像為數(shù)據(jù)源，ENVI5.3、ArcGIS為軟件平臺(tái)，以石城縣縣城總體規(guī)劃方案（2013—2030）中的縣城規(guī)劃區(qū)域?yàn)檠芯繀^(qū)，采用輻射傳輸方程法對(duì)該地區(qū)進(jìn)行地表溫度反演，此基礎(chǔ)上對(duì)該地區(qū)城市熱島效應(yīng)與各地表參數(shù)進(jìn)行定量化相關(guān)性研究.這對(duì)熱島效應(yīng)的城市的合理規(guī)劃、城市熱島效應(yīng)的緩解具有重要的意義，為石城縣創(chuàng)建森林城市、園林縣城提供參考.

1 研究區(qū)概況及數(shù)據(jù)源

1.1 研究區(qū)概況

石城縣位于江西省贛州市的東北部，石城縣縣城總體規(guī)劃方案（2013—2030）中將石城縣定為贛江源頭生態(tài)宜居宜業(yè)城市、贛東南交通樞紐和商貿(mào)物流中心、贛州市新型工業(yè)城，有著“三區(qū)兩片、一帶兩廊、三軸三心”的規(guī)劃結(jié)構(gòu).石城縣位于中亞熱帶季風(fēng)濕潤(rùn)氣候區(qū)，日照充足，降水充沛，四季分明，全年平均日照時(shí)長(zhǎng)為1825.6 h，年平均溫度和年降水量分別介于 15.5～19.1℃和 1400～1900 mm之間，兩者均呈東高西低、南高北低的空間格局.石城縣屬東南丘陵低山地區(qū)，中部地勢(shì)平坦，東北部群山包圍，西南部丘陵連綿，是一片原生態(tài)的生命綠洲，生態(tài)環(huán)境基礎(chǔ)良好.近些年來(lái)城市化擴(kuò)張對(duì)資源具有嚴(yán)重的破壞，城市熱島效應(yīng)越發(fā)顯著.

1.2 數(shù)據(jù)源及數(shù)據(jù)處理

選用地理空間數(shù)據(jù)云網(wǎng)站(http：//www.gscloud.cn)下載的2009年10月26日、2014年10月8日和2017年11月1日的三期Landsat遙感數(shù)據(jù)作為主要數(shù)據(jù)源，云量較少，影響質(zhì)量較好；影像空間分辨率為30 m，Landsat TM5第6波段（熱紅外波段）分辨率為120 m，Landsat 8第10波段 （熱紅外波段）分辨率為100 m；非遙感數(shù)據(jù)主要包括石城縣航空影像數(shù)據(jù)和1∶500的地形圖.

在ENVI5.3軟件中，對(duì)遙感數(shù)據(jù)進(jìn)行讀取、輻射定標(biāo)、大氣校正等處理；利用石城縣1∶500的規(guī)劃圖對(duì)各期影像進(jìn)行幾何校正；利用石城縣城市矢量圖裁剪出研究區(qū)域遙感影像.研究區(qū)概況（航拍圖）如圖1所示.

圖1 研究區(qū)概況

2 研究方法

2.1 地表溫度反演

地表溫度作為遙感影像反演結(jié)果，是影像在剔除大氣以及傳感器自身的各種影響后，考慮到不同地表比輻射率的情況下得到的地表輻射量所對(duì)應(yīng)的溫度[19].由于遙感影像的最小單元是像元，通常都是混合像元；大氣作為一個(gè)復(fù)雜的系統(tǒng)，很難通過(guò)模型來(lái)進(jìn)行參數(shù)估計(jì)，故而增加了遙感反演地表溫度的難度.地表溫度又稱地面輻射溫度，指的是與地表有相同輻射亮度的黑體溫度，與傳感器在衛(wèi)星高度所觀測(cè)到的熱輻射強(qiáng)度相對(duì)應(yīng)的溫度.地表亮度溫度與地表溫度在數(shù)值上雖不相等，但二者相關(guān)性強(qiáng)，相比地表亮度溫度，地表溫度包含更多大氣信息，在一定程度上能夠模擬地表溫度的空間分布，更有利于城市熱島效應(yīng)的反演與評(píng)估[20].

常用的基于熱紅外數(shù)據(jù)的地表溫度反演方法有：輻射傳輸方程算法、單窗算法和分裂窗算法等，一些學(xué)者對(duì)不同算法的反演精度進(jìn)行了定量對(duì)比分析,結(jié)果表明輻射傳輸方程法精度最高[21-23].輻射傳輸方程法（也稱大氣校正法）是目前地表溫度反演主要算法之一，主要使用大氣模型模擬大氣對(duì)地表熱輻射的影響，包括估計(jì)大氣對(duì)熱輻射傳導(dǎo)的吸收作用以及大氣自身的向上和該方法主要是依靠MODTRAN輻射傳輸代碼和整體算法獲取特定點(diǎn)的透過(guò)率和大氣上下輻射亮度，估計(jì)誤差較小，僅有3%左右[24].測(cè)算過(guò)程主要包括：

1）計(jì)算歸一化植被指數(shù)NDVI

式（1）中：ρNIR表示近紅外波段，ρR表示紅波段.

2）計(jì)算植被覆蓋度

式（2）中：b1表示公式（1）計(jì)算的結(jié)果.

3）計(jì)算地表比輻射率

由于Landsat TM的第6波段與Landsat8 TIRS的第10波段波譜范圍較為相近，處理過(guò)程中均可采用Sobrino[25]提出的NDVI閾值法計(jì)算地表比輻射率，如公式（3）所示：

式（3）中：Pv是植被覆蓋度，用公式（2）計(jì)算的植被覆蓋度結(jié)果.

4）計(jì)算同溫度下黑體輻射亮度值

式（4）中：b1為地表比輻射率結(jié)果，b2為熱紅外波段經(jīng)過(guò)輻射定標(biāo)的圖像，Lu為大氣向上輻射亮度、Ld為大氣向下輻射亮度、t為大氣在熱紅外波段的透過(guò)率.

輻射傳輸方程法大氣參數(shù)值包括大氣上行輻射、大氣下行輻射和大氣透過(guò)率可以在NASA公布的網(wǎng)站查詢（http：//atmcorr.gsfc.nasa.gov），以 2014 年Landsat遙感數(shù)據(jù)為例，輸入成影時(shí)間：2014-10-08 02：45 和中心經(jīng)緯度（Lat：26，Lon：115.8），以及其他相應(yīng)的參數(shù)，得到大氣剖面信息如圖2所示.

圖2 輻射傳輸方程法大氣參數(shù)值

5）計(jì)算地表溫度

根據(jù)普朗克公式的反函數(shù)，地表溫度計(jì)算公式如式（5）所示：

式（5）中：b1為同溫度黑體輻射亮度值，LST為真實(shí)的地表溫度值，單位是℃，參數(shù)K1、K2為普朗克校正常量（表1），可從遙感影像元數(shù)據(jù)中獲取.

表1 傳感器參數(shù)設(shè)置

2.2 熱島效應(yīng)等級(jí)劃分

對(duì)地表溫度進(jìn)行歸一化處理，使得處理后像元數(shù)據(jù)區(qū)間為[0，1]；按公差為0.2的等差級(jí)數(shù)將熱島效應(yīng)密度分割為五個(gè)等級(jí)，其中，0.0～0.2表示強(qiáng)綠島區(qū)，0.2～0.4表示綠島區(qū)，0.4～0.6表示中等區(qū)，0.6～0.8表示熱島區(qū)，0.8～1.0表示強(qiáng)熱島區(qū).歸一化處理方法如公式（6）所示：

式（6）中：N表示熱紅外影像中第i個(gè)像元?dú)w一化地表溫度值；Tmax、Tmin分別表示地表溫度的最大值和最小值；Ti則表示第i個(gè)像元的地表溫度值.

2.3 地表反演溫度驗(yàn)證

對(duì)于地表溫度遙感反演結(jié)果的驗(yàn)證，最普遍的方法就是將Landsat熱紅外影像反演的地表溫度與衛(wèi)星飛過(guò)該地區(qū)上空時(shí)刻的地面測(cè)量(對(duì)應(yīng)像元尺度)溫度值進(jìn)行比較.但由于實(shí)測(cè)值不容易獲取，這里選用MODIS地表溫度合成產(chǎn)品（MODLT1F，精度為1 km）作為標(biāo)準(zhǔn)值來(lái)驗(yàn)證本文的亮溫反演結(jié)果.以2014年數(shù)據(jù)為例，LST反演結(jié)果與同天MOD11A1進(jìn)行對(duì)比，經(jīng)過(guò)數(shù)據(jù)統(tǒng)計(jì)得出，Landsat反演結(jié)果顯示研究區(qū)地表溫度最高值為38.42℃，最低溫度為19.51℃，MODIS地表溫度的最高值為36.54℃，最低溫度為18.82℃，空間分布整體一致性較高.由遙感物理知識(shí)可知，傳感器在接收地表輻射值時(shí)，除地表輻射值外，還包括了大氣輻射的地表反射以及大氣散射等多來(lái)源輻射值，使得地表溫度值普遍高于地表溫度值.

為進(jìn)一步驗(yàn)證反演結(jié)果，這里分別從建筑用地、水體、植被和裸地四類不同地表覆蓋類型中各選取15個(gè)樣本點(diǎn)，并統(tǒng)計(jì)各個(gè)樣本點(diǎn)地表溫度反演值與對(duì)應(yīng)的MODIS地表溫度值，用兩者間差值的絕對(duì)值來(lái)表示兩溫度值的差異.通過(guò)統(tǒng)計(jì)分析可得，兩者差異的平均值為0.94℃，反演總體精度較高.其中不同地表覆蓋類型反演結(jié)果差異較大，其中溫度差異最大值為2.02℃，出現(xiàn)在建筑用地覆蓋類型樣本中；最小值為0.45℃，出現(xiàn)在植被覆蓋類型樣本中，且整體上植被類型溫度反演較好，樣本差異絕對(duì)值均小于1℃.

2.4 地表參數(shù)的計(jì)算與提取

地表參數(shù)是導(dǎo)致城市熱環(huán)境發(fā)生變化的主要因素，比如歸一化建筑用地指數(shù)、土壤調(diào)節(jié)植被指數(shù)、改進(jìn)歸一化水體指數(shù)及裸土指數(shù)等幾個(gè)主要的相關(guān)指數(shù).城市熱島效應(yīng)的形成與建筑用地、水體、植被及裸地之間相互轉(zhuǎn)換緊密相關(guān)[26]，根據(jù)石城縣城區(qū)地表覆蓋情況，主要覆蓋類型包括建筑用地、水體、植被以及部分裸地.為進(jìn)一步研究該區(qū)域地表參數(shù)的變化和地表參數(shù)對(duì)城市熱島的影響，選用石城縣2014年影像數(shù)據(jù)及熱島反演結(jié)果 （云量更少（0.18），能夠更準(zhǔn)確的反應(yīng)熱島效應(yīng)與各地表參數(shù)之間的相關(guān)性），選取能夠反應(yīng)建筑用地、水體、植被及裸地的相關(guān)指數(shù)進(jìn)行研究.

1）歸一化建筑用地指數(shù)NDBI

建筑用地指數(shù) （Normalized Difference Buiding Index，NDBI）是基于歸一植被指數(shù)提出的，根據(jù)不同地物光譜特征差異準(zhǔn)確反應(yīng)建筑用地信息，在地表植被覆蓋較高地區(qū)可提高提取精度，其數(shù)值越大表明建筑用地比例越高.

式（7）中：NDBI表示建筑用地指數(shù)；ρSWIR2表示短波紅外波段反射率；ρNIR表示近紅外波段反射率.

2）土壤調(diào)節(jié)植被指數(shù)SAVI

土壤調(diào)節(jié)植被指數(shù) （Soil Adjusted Vegetation Index，SAVI）是在NDVI公式基礎(chǔ)上引入土壤調(diào)節(jié)因子，以此來(lái)消除土壤背景對(duì)植被光譜信號(hào)的影響.同時(shí)SAVI是最適合于研究城市建成區(qū)，其植被覆蓋率降低至15%.

式（8）中：SAVI表示土壤調(diào)節(jié)植被指數(shù)；ρNIR表示近紅外波段反射率；ρR表示紅光波段反射率；L表示土壤調(diào)節(jié)因子，相關(guān)研究表明，當(dāng)L=0.5時(shí)，可以將土壤的亮度差異降到最低，故這里取L=0.5.

3）改進(jìn)歸一化水體指數(shù)MNDWI

改進(jìn)歸一化差異水體指數(shù) （Modified Normalized Difference Water Index，MNDWI） 采用了歸一化的比值運(yùn)算，能夠去除地形等外界環(huán)境的影響，在抑制植被噪音信息的同時(shí)能夠較好的抑制土壤和建筑用地信息，從而有利于水體信息的提取，特別是提取城鎮(zhèn)范圍的水體；另外，MNDWI容易區(qū)分陰影和水體，解決了水體提取中難消除陰影的難題.

式（9）中：MNDWI表示改進(jìn)的歸一化水體指數(shù)；ρG表示綠光波段反射率；ρSWIR1表示短波熱紅外波段反射率.

4）裸土指數(shù)BSI

在裸土指數(shù)（Bare Soil Index，BSI）圖像上，裸地亮度最高，能夠用于某地區(qū)裸地信息的提取.表達(dá)式如公式（10）所示.

式（10）中：BSI表示裸土指數(shù)；ρSWIR1表示短波熱紅外波段反射率；ρR示紅光波段反射率；ρNIR表示近紅外波段反射率.

通過(guò)公式（7）到公式（10），在 ENVI5.3 軟件中依次進(jìn)行Landsat 8波段運(yùn)算命令，分別提取到4種指數(shù)分布圖（圖 3），建設(shè)用地和裸土的 NDBI、BSI指數(shù)值較高，SAVI、MNDWI指數(shù)值較低；植被的NDBI、MNDWI 指數(shù)值較低，SAVI、BSI 指 數(shù) 值 較高； 水體的 MNDWI指數(shù)值高，NDBI、SAVI、BSI指數(shù)值均較低.

圖3 2014年石城縣4種指數(shù)分布

3 結(jié)果與分析

3.1 熱島效應(yīng)等級(jí)時(shí)空分布

由于3幅影像獲取的時(shí)相不完全一致，不能直接將絕對(duì)地表溫度值進(jìn)行對(duì)比，月份及季節(jié)的差異只影響絕對(duì)溫度值，不影響溫度的分布規(guī)律，故對(duì)各年份影像的地表溫度進(jìn)行正規(guī)化；利用密度分割技術(shù)將正規(guī)化以后的數(shù)據(jù)密度分割成強(qiáng)綠島區(qū)、綠島區(qū)、中等區(qū)、熱島區(qū)和強(qiáng)熱島區(qū)5個(gè)等級(jí)（圖4），提取各等級(jí)分區(qū)面積（表2），進(jìn)一步分析石城縣城市熱島時(shí)空分布變化.

圖4 石城縣城市熱島效應(yīng)等級(jí)影像分布

表2 熱島等級(jí)區(qū)面積及所占比例

由圖4可知，各年研究區(qū)綠島區(qū)及中等區(qū)分布較廣，由縣城中心到郊區(qū)熱島強(qiáng)度逐漸減弱，其中城市熱島主要位于琴江西岸琴江鎮(zhèn)的城西片區(qū)以及縣城南部.強(qiáng)綠島區(qū)分布較少，主要集中在植被覆蓋度較高的區(qū)域，如東城郊野生態(tài)主題公園以及西南部花園森林公園局部；熱島區(qū)主要分布于琴江兩岸地區(qū)，近年來(lái)還集中于縣城南部及東南部地區(qū).由于縣城中心坐落于此，具有密集的城市建筑群，且工業(yè)企業(yè)相對(duì)較多，土地利用以建設(shè)用地為主，包括行政、文化、教育、商業(yè)、體育、工業(yè)、居住等用地；隨著城市的擴(kuò)張及發(fā)展，琴江兩岸、縣城南部建筑群落不斷增加，縣城南部橋村周邊分布了工業(yè)基地、物流倉(cāng)儲(chǔ)基地、商務(wù)設(shè)施及各級(jí)居民地，熱島區(qū)面積有所覆蓋，呈現(xiàn)放射狀由熱島中心向周邊環(huán)境遞減；強(qiáng)熱島區(qū)范圍較小但分布集中，主要位于琴江西岸的城西片區(qū)，以老城區(qū)和工業(yè)區(qū)為主，尤以石城縣工業(yè)園區(qū)之一的古樟工業(yè)園為主，眾多工業(yè)企業(yè)集中于此，加之流動(dòng)人口較多，故而該區(qū)熱島強(qiáng)度較為嚴(yán)重.

由表2可知，各階段性縣城熱島等級(jí)分布差異明顯.2009—2014年，中等區(qū)面積變化最大，覆蓋面積減少了9.19 km2；強(qiáng)熱島區(qū)面積增量最大，東部植被面積增加的同時(shí)，縣城西南部的花園森林公園范圍強(qiáng)綠島區(qū)面積也有所增加；2014—2017年，強(qiáng)綠島區(qū)面積變化最大，減少了4.36 km2，中等區(qū)變化量次之，面積增加了3.97 km2，強(qiáng)熱島區(qū)面積變化量最小.石城縣2009—2017年近10年中，中等區(qū)面積變化量最大，范圍減少了5.22 km2；綠島區(qū)面積增量最大，增加了4.30 km2，近年來(lái)隨著石城縣生態(tài)縣建設(shè)強(qiáng)度的不斷增加，東、西城郊野生態(tài)公園、通天寨等各公園景區(qū)生態(tài)范圍不斷增加；強(qiáng)熱島區(qū)范圍由2009年的0.23 km2增加到0.37 km2，面積增加了0.14 km2.

3.2 熱島效應(yīng)與各參數(shù)指數(shù)關(guān)系的建立

為進(jìn)一步定量分析地表溫度與各指數(shù)之間的關(guān)系，確保統(tǒng)計(jì)分析的客觀性,結(jié)合研究區(qū)域面積的大小，對(duì)2014年的各指數(shù)圖像及熱島溫度反演圖進(jìn)行隨機(jī)采集40個(gè)樣本點(diǎn)，并對(duì)樣本點(diǎn)做統(tǒng)計(jì)和回歸分析（圖5）.

圖5 地表溫度與各指數(shù)相關(guān)性分析

地表溫度與四種指數(shù)回歸分析結(jié)果與散點(diǎn)圖結(jié)果保持一致，地表溫度與SAVI、MNDWI、NDBI和BSI指數(shù)之間存在顯著相關(guān)性.石城縣地表溫度LST與NDBI呈正相關(guān)，相關(guān)系數(shù)r為0.86，NDBI數(shù)值每增加0.1，地表溫度增加2.63℃（圖5（a）），說(shuō)明城市建筑用地面積的增加對(duì)地表溫度具有增強(qiáng)作用.地表溫度LST與SAVI呈負(fù)相關(guān)，相關(guān)系數(shù)為0.92，SAVI數(shù)值每增加0.1，地表溫度減少2.48℃（圖5（b）），可以提高城市綠地覆蓋率，對(duì)熱島效應(yīng)具有緩解作用.地表溫度LST與MNDWI呈負(fù)相關(guān)，相關(guān)系數(shù)為0.77，MNDWI數(shù)值每增加0.1，地表溫度減少 4.26℃（圖 5（c）），說(shuō)明地表水體含量以及植被含水量也是影響溫度變化的重要因素.地表溫度LST與BSI呈正相關(guān)，相關(guān)系數(shù)為0.87，BSI數(shù)值每增加0.1，地表溫度增加3.24℃（圖5（d）），說(shuō)明城市裸土面積等不透水下墊面的增長(zhǎng)會(huì)使城市溫度增高,熱島效應(yīng)加劇.

4 結(jié) 論

地表溫度作為地表物理過(guò)程中的一個(gè)重要參數(shù)，是全球物質(zhì)能量循環(huán)、能量平衡、氣候變化的重要因素，是自然科學(xué)研究的一個(gè)重要領(lǐng)域.隨著遙感技術(shù)的快速發(fā)展，目前利用單一的熱紅外通道數(shù)據(jù)反演地表溫度研究較多，通過(guò)獲取過(guò)境時(shí)刻的熱紅外輻射量，利用一定的反演算法估算地表溫度，進(jìn)而反應(yīng)地表溫度.

文中選用多時(shí)相Landsat遙感數(shù)據(jù)，通過(guò)對(duì)石城縣城市熱島效應(yīng)的測(cè)算以及相關(guān)地表參數(shù)的計(jì)算，分析了石城縣2009年、2014年及2017年城市熱島效應(yīng)各等級(jí)時(shí)空分布變化，及其與相關(guān)參數(shù)指數(shù)之間的相關(guān)性，得出如下結(jié)論：

1）石城縣城市熱島效應(yīng)等級(jí)空間分布分區(qū)顯著，熱島效應(yīng)強(qiáng)度由縣城中心向周邊郊區(qū)逐漸降低；綠島區(qū)及中等區(qū)分布較廣，強(qiáng)綠島區(qū)主要分布于縣城的東城郊野生態(tài)公園及西南部的花園森林公園等植被覆蓋較多的地區(qū)；熱島區(qū)主要集中于中心縣城的琴江兩岸地區(qū)以及縣城南部地區(qū)；強(qiáng)熱島區(qū)主要分布于琴江鎮(zhèn)琴江西岸的老城區(qū)及工業(yè)區(qū).

2）2009—2017年近10年中，研究區(qū)熱島強(qiáng)度等級(jí)空間分布變化顯著.隨著石城縣生態(tài)縣建設(shè)的實(shí)施，綠島區(qū)面積持續(xù)增加，強(qiáng)綠島區(qū)范圍也有所擴(kuò)大；熱島區(qū)面積持續(xù)減少，但空間分布區(qū)域集中化；同時(shí)強(qiáng)熱島區(qū)面積也有所增加，以工業(yè)區(qū)和住宅地為主.

3）城市地表溫度與SAVI和MNDWI指數(shù)呈明顯的線性負(fù)相關(guān)關(guān)系，與NDBI和BSI指數(shù)呈明顯的正相關(guān)關(guān)系.為定量化分析地表參數(shù)之間的關(guān)系，通過(guò)回歸分析擬合地表溫度與四種地表參數(shù)指數(shù)之間的相關(guān)性函數(shù)，求解相關(guān)系數(shù)，表明適當(dāng)提高綠地覆蓋率、增加水域面積，合理規(guī)劃城市道路及建筑物有助于減緩熱島效應(yīng).