曹倩倩,劉瑞芬

基于不同土地利用類型的植被指數(shù)與地表溫度的關(guān)系——以張家口市為例

曹倩倩,劉瑞芬*

河湖生態(tài)修復(fù)與藻類利用湖北省重點(diǎn)實(shí)驗(yàn)室, 湖北工業(yè)大學(xué)土木建筑與環(huán)境學(xué)院, 湖北 武漢 430068

本文以河北省張家口市區(qū)六個(gè)區(qū)作為研究對(duì)象，基于Landsat 8遙感影像、土地利用數(shù)據(jù)和ArcGIS空間分析技術(shù)，獲取歸一化植被指數(shù)()以及陸地表面溫度（）的空間分布，并采用多重比較量化分析六種土地利用類型下、值的差異，最后通過(guò)線性回歸方法分析不同土地利用類型與之間的關(guān)系。結(jié)果表明，張家口市區(qū)的和的空間上值的分布呈相反趨勢(shì)。對(duì)于而言，城鄉(xiāng)工礦居民用地與耕地、草地、林地等土地利用類型存在顯著性差異；對(duì)于而言，耕地與草地、林地、城鄉(xiāng)工礦居民用地等土地利用類型存在顯著性差異。不同土地利用類型下的和呈現(xiàn)出明顯的負(fù)相關(guān)關(guān)系，判定系數(shù)2=0.519。結(jié)果可以為張家口市城市生態(tài)規(guī)劃以及緩解熱島效應(yīng)提供科學(xué)依據(jù)。

土地利用; 地表溫度; 植被指數(shù)

在快速城鎮(zhèn)化的過(guò)程中，水泥和瀝青的硬化人工地表逐漸取代了自然地表，城市土地覆蓋的局部景觀格局的改變使得整個(gè)區(qū)域的生態(tài)環(huán)境狀況發(fā)生了變化[1]，以及人口的快速增長(zhǎng)和工業(yè)的發(fā)展，大氣環(huán)境、地表結(jié)構(gòu)和人為熱排放等多種因素相互作用，從而使得城市內(nèi)部溫度要高于郊區(qū)的溫度，形成了熱島效應(yīng)[2]。植被覆蓋可以通過(guò)蒸騰作用影響著地表的能量平衡[3]，植被的光合作用還可以將光能轉(zhuǎn)化為化學(xué)能[4]，通過(guò)蒸騰作用和光合作用影響著地氣之間的熱交換形式和減小太陽(yáng)輻射效應(yīng)，從而降低地表溫度，進(jìn)而緩解城市熱島效應(yīng)。植被指數(shù)()用來(lái)衡量植被覆蓋度。城市的不同土地利用類型如建設(shè)用地、耕地、林地有不同的植被覆蓋度，對(duì)地表溫度()也有著不同的影響，因此量化城市尺度下不同土地利用類型的植被覆蓋度及植被覆蓋度對(duì)地表溫度的影響可以指導(dǎo)土地利用布局，對(duì)緩解城市熱島效應(yīng)和構(gòu)建生態(tài)友好城市有著重要的意義。

近年來(lái)隨著遙感技術(shù)的發(fā)展，星載遙感傳感器可以以較短的周期和成本快速準(zhǔn)確的獲取地表的下墊面的溫度特征[5]，因此在土地利用及其覆蓋變化與的關(guān)系的研究中得以廣泛應(yīng)用。牟雪潔[6]、彭文甫[7]、劉宇[8]等人從熱島效應(yīng)分布特征的角度探究了研究土地利用與之間的關(guān)系；彭璐[3]、高尚[4]等人使用遙感影像研究的變化值對(duì)變化的影響，量化變化值和變化值之間的線性關(guān)系，研究表明不同的范圍對(duì)的影響不同；Wilson JS[9]、岳文澤[10]、Lo CP[11]、周媛[12]等人運(yùn)用遙感影像和GIS空間分析技術(shù)探討了不同土地利用類型下的、值的差異，以及不同土地利用與之間的關(guān)系，研究均表明、在不同土地利用下存在差異，并且與存在著明顯的負(fù)相關(guān)關(guān)系，該研究為城市規(guī)劃以及綠地系統(tǒng)規(guī)劃提供技術(shù)支撐和依據(jù)。但是這些研究的研究區(qū)域多集中于南方植被茂盛的城市，例如福建省廈門市、平澤島以及上海市,對(duì)于北方植被較稀疏的城市研究比較少。本文以北方城市張家口市區(qū)為例，利用張家口市區(qū)一級(jí)土地利用類型shp文件和Landsat8遙感影像進(jìn)行和地表溫度的反演，量化與的關(guān)系，進(jìn)而探討不同土地利用類型下的與之間的關(guān)系，對(duì)張家口緩解城市熱島效應(yīng)和建成更加生態(tài)友好城市的生態(tài)規(guī)劃提出依據(jù)。

1 材料與方法

1.1 研究區(qū)概況及數(shù)據(jù)來(lái)源

張家口市屬于河北省西北部區(qū)域（如圖1），介于113°50′～116°30′E，39°30′～42°10′N之間。張家口市整體地勢(shì)西北高、東南低，有壩上高原區(qū)、壩下山間盆地區(qū)兩大典型地貌區(qū)，兩種地貌以陰山山脈大馬群山分水嶺為界，屬溫帶大陸性季風(fēng)氣候，有干旱、風(fēng)沙、霜凍和雹災(zāi)等主要災(zāi)害。

本文研究區(qū)域?yàn)閺埣铱谑袇^(qū)，包括橋東區(qū)、橋西區(qū)、宣化區(qū)、崇禮區(qū)、下花園區(qū)、萬(wàn)全區(qū)等六個(gè)區(qū)（如圖1），區(qū)域面積達(dá)6429 km2。土地利用類型主要分為耕地、草地、林地、城鄉(xiāng)工礦居民用地、水域、未利用土地等六大類型，占地比分別為36.5%、32%2、24.49%、6%、1%、0.01%。

圖 1 張家口市區(qū)的地理位置以及各個(gè)分區(qū)

本研究以2015年7月12日河北省張家口市Landsat 8遙感影像作為數(shù)據(jù)源，其空間分辨率為30 m，包括陸地成像儀OLI含9個(gè)波段和熱紅外傳感器TIRS含2個(gè)波段等11個(gè)波段。應(yīng)用Landsat8 OLI的第4、5波段進(jìn)行研究區(qū)的波段運(yùn)算，同時(shí)利用Landsat8 OLI獲得的第4、5波段運(yùn)算得到的和Landsat8 TIRS10反演研究區(qū)的。同時(shí)結(jié)合張家口市2015年的土地利用現(xiàn)狀shp文件，將研究區(qū)的土地利用類型分為耕地、草地、林地、城鄉(xiāng)工礦居民用地、水域、未利用土地等六種類型（如圖2）。從研究區(qū)的土地利用布局來(lái)看，林地主要集中于崇禮區(qū)，草地主要集中在萬(wàn)全區(qū)和橋東區(qū)兩個(gè)區(qū)域，建設(shè)用地大多數(shù)分布于橋東區(qū)、橋西區(qū)、宣化區(qū)，耕地和草地為主要用地類型，面積百分比為68.5%。

圖 2 張家口市區(qū)土地利用類型圖

1.2 植被指數(shù)及地表溫度計(jì)算

歸一化植被指數(shù)提取，有研究表明指標(biāo)可以很好的表征植被覆蓋狀況和地表輻射溫度[11]。通過(guò)綠色植物的紅外波段和近紅外波段反射光譜的差異可以進(jìn)行的計(jì)算。波段運(yùn)算時(shí)，先要對(duì)Landsat8 OLI的遙感影像進(jìn)行輻射定標(biāo)和大氣校正，然后對(duì)經(jīng)過(guò)處理的遙感影像進(jìn)行波段運(yùn)算。的計(jì)算公式如（1）所示。計(jì)算出來(lái)的值介于-1～1，一般認(rèn)為在植物的生長(zhǎng)季節(jié)>0，表示有植被覆蓋，值越大說(shuō)明植被越茂盛，超過(guò)0.5表明植被生長(zhǎng)茂盛[12]。

式中：5，4分別表示紅光波段Band4和近紅外波段Band5的灰度值。

地表溫度反演，本研究利用Landsat8 TIRS10波段所接收到的地面各處的熱輻射值進(jìn)行定量反演，采用胡德勇[13]的改進(jìn)的Landsat8陸地表面溫度計(jì)算的方法：利用ENVI軟件進(jìn)行遙感影像處理，首先分別對(duì)Landsat8 OLI和Landsat8 TIRS10數(shù)據(jù)進(jìn)行輻射定標(biāo)，接著對(duì)Landsat8 OLI數(shù)據(jù)進(jìn)行大氣校正、利用其第4、5波段進(jìn)行和植被覆蓋度計(jì)算（見(jiàn)式2）、再進(jìn)行地表比輻射率的計(jì)算（見(jiàn)式3、4、5），結(jié)合基于Landsat8 TIRS10數(shù)據(jù)進(jìn)行輻射亮度的計(jì)算以及從航天局獲取的大氣剖面參數(shù)，進(jìn)行同溫度下黑體輻射亮度的計(jì)算（見(jiàn)等式6），最后運(yùn)用波段運(yùn)算得出（見(jiàn)式7）。

基于像元二分模型估算植被覆蓋度計(jì)算，NDVI裸土或者建筑表面的值，NDVI為全植被覆蓋區(qū)的值。本研究區(qū)域?qū)嶋H的植被覆蓋狀況將NDVI設(shè)置為0.05，將NDVI設(shè)置為0.7，當(dāng)像元的>NDVI時(shí)，代表全植被覆蓋，P=1；當(dāng)像元的>NDVI時(shí)，代表裸土，P=0。

城鎮(zhèn)地表、自然地表兩類混合像元地表比輻射率計(jì)算。植被、建筑表面、裸土的純凈像元的地表發(fā)射率分別為ε=0.986、ε=0.970、ε=0.972，城鎮(zhèn)地表可大致看作是建筑表面和植被組成的，自然地表可看作是植被和裸土組成的。

地表比輻射率計(jì)算公式中，溫度比率R定義為R=(T/)4，植被R、建筑表面R、裸土R的溫度比率可以使用上式的公式求出[14]：

式中T，T分別代表亮溫和地表溫度（K），代表TIRS波段；L(T)為亮溫為T時(shí)傳感器接收到的輻射能量(W?m-2?sr-1?μm-1)；L(T)為為T時(shí)的輻射能量(W?m-2?sr-1?μm-1)；τ為大氣透過(guò)率（無(wú)量綱）；ε代表地表發(fā)射率（無(wú)量綱）；L↓為大氣下行輻射(W?m-2?sr-1?μm-1)；L↑為大氣上行輻射(W?m-2?sr-1?μm-1)。通過(guò)航天局獲取的數(shù)據(jù)張家口市區(qū)的τ=0.79，L↓=2.88，L↑=1.69。

地表溫度表達(dá)式如上式所示[15]，Landsat8中，通常取1=774.8853，2=1321.0789，再將計(jì)算結(jié)果的單位轉(zhuǎn)化成攝氏度。

1.3 研究方法

研究主要用到的數(shù)學(xué)分析方法有post-hoc多重比較分析、線性回歸分析。多重比較方法是通過(guò)對(duì)總體均值之間的兩兩比較，來(lái)檢驗(yàn)各個(gè)因素總體均值之間相等或者差異的關(guān)系。線性回歸分析是利用數(shù)學(xué)關(guān)系來(lái)確定兩種或者兩種以上的變量之間的相互依賴的定量關(guān)系的統(tǒng)計(jì)分析方法。

本研究中使用post-hoc多重比較分析的方法，分析耕地、草地、林地、城鄉(xiāng)工礦居民用地、水域兩兩土地利用之間總體均值差，分析兩種土地利用類型之間、的差異程度。使用線性回歸分析，分析基于土地利用類型的和之間的線性關(guān)系，判斷對(duì)不同土地利用類型的的影響程度。

構(gòu)建特征剖面，以張家口市政府（114.892，40.774）為中心原點(diǎn)坐標(biāo)，分別向東西南北四個(gè)方向延展至區(qū)域邊界，最終形成東西、南北兩個(gè)方向的軸線（如圖3），在ArcGIS軟件中進(jìn)行數(shù)據(jù)處理，獲得東西、南北兩條軸線上和的值剖面。

圖 3 東西、南北方向特征剖面軸線

2 結(jié)果與分析

2.1 NDVI與LST計(jì)算結(jié)果與剖面變化情況

張家口市區(qū)與空間分布如圖4所示，二者值的空間分布具有相反的趨勢(shì)，值大（?。┑膮^(qū)域一般較?。ù螅?。具體來(lái)說(shuō)，在研究區(qū)域東北部的值相對(duì)較大，相應(yīng)的的值比較??；的較大值出現(xiàn)在區(qū)域中部和南部。同樣的，和空間分布呈相反的趨勢(shì)在同類型研究中均有發(fā)現(xiàn)[16,17]，充分說(shuō)明區(qū)域的植被覆蓋會(huì)直接影響到，主要通過(guò)影響熱輻射、熱動(dòng)力和土壤水分等地表特征來(lái)使得發(fā)生變化[18]。

(a)NDVI空間分布圖；(b)LST空間分布圖

和在東西和南北兩個(gè)方向剖面的值的變化趨勢(shì)情況如圖5。東西方向最大值為0.79，其土地利用類型為林地，最低值為0.09，其土地利用類型為城鄉(xiāng)工礦居民用地。還可觀察到東西方向剖面居中區(qū)域，當(dāng)林地與城鄉(xiāng)工礦居民用地相鄰時(shí)，林地的要明顯高于城鄉(xiāng)工礦居民用地的，如圖5(a)所示。的最高值及次高值分別為48.88 ℃和48.59 ℃，所對(duì)應(yīng)的土地利用類型為耕地和城鄉(xiāng)工礦居民用地，最低值為22.63 ℃，所對(duì)應(yīng)的土地利用類型為林地。張家口市的耕地類型有水田和旱地兩種，種植的作物類型也較多樣，不同的作物類型對(duì)溫度造成的影響也不同，導(dǎo)致東西方向上耕地溫度值（范圍為24.96～48.88 ℃）波動(dòng)比較大，出現(xiàn)耕地溫度偶然值是最高值的現(xiàn)象。南北方向的最大值及最小值分別為0.71和0.07，所對(duì)應(yīng)的土地利用類型為林地和城鄉(xiāng)工礦居民用地，與東西向的極值相同。南北方向的最高值為47.1 ℃，對(duì)應(yīng)的土地利用類型為草地，最低值及次低值分別為30.11 ℃和30.29 ℃，對(duì)應(yīng)的土地利用類型為耕地和林地。草地的出現(xiàn)最高值一方面可能是受到周邊土地利用類型的影響，如圖5(d)右側(cè)區(qū)域，城鄉(xiāng)工礦居民用地周邊的草地較高，另一方面可能草地生長(zhǎng)狀況不太好，出現(xiàn)了相對(duì)裸露的地表，較低就出現(xiàn)草地的溫度最高的現(xiàn)象。

從圖中可以看出，東西和南北方向的與總體呈現(xiàn)出相反的變化趨勢(shì)，增長(zhǎng)時(shí)，減少，在高值區(qū)時(shí)對(duì)應(yīng)的處于低值區(qū)，和的折線圖關(guān)于水平方向?qū)ΨQ。但對(duì)于同一種土地利用類型（如林地），當(dāng)其最大時(shí)對(duì)應(yīng)的不一定最小，為深入分析不同土地利用類型下的、的差異性，需對(duì)研究區(qū)的土地利用類型的和進(jìn)行統(tǒng)計(jì)分析。

(a)東西方向NDVI變化趨勢(shì) (c)南北方向NDVI變化趨勢(shì)(b)東西方向LST變化趨勢(shì) (d)南北方向LST變化趨勢(shì)

2.2 不同土地利用類型的NDVI與LST統(tǒng)計(jì)分析

本文利用ArcGIS軟件對(duì)土地利用類型數(shù)據(jù)與、的柵格圖像分別進(jìn)行疊加，分區(qū)統(tǒng)計(jì)不同土地利用類型下的、的平均值（對(duì)應(yīng)區(qū)域的像元均值）和標(biāo)準(zhǔn)差，統(tǒng)計(jì)結(jié)果如圖6所示。結(jié)果表明，不同土地利用類型下的、的平均值有所差異。林地的平均值最大為0.65，耕地和草地的平均值均為0.46，城鄉(xiāng)工礦居民用地平均值最低為0.29。林地的平均值最小為34.21 ℃，耕地和草地的平均值也較接近，城鄉(xiāng)工礦居民用地平均值最高為40.84 ℃。其中林地的最大對(duì)應(yīng)的最低，城鄉(xiāng)工礦居民用地的最低對(duì)應(yīng)的最大。有研究發(fā)現(xiàn)[19]，林地在所有土地利用類型中，對(duì)溫度的影響最大，與本文結(jié)果一致。

圖 6 不同土地利用類型下的NDVI、LST平均值及標(biāo)準(zhǔn)差

從上文結(jié)果來(lái)看，不同的土地利用類型上的和的平均值存在差異，但是平均值并不能表示這種差異的顯著程度。為深入了解不同土地利用類型之間兩兩的差異顯著程度，采用了TamhanceT2 post-hoc多重比較對(duì)東西、南北兩個(gè)軸向上土地利用類型的、的均值差進(jìn)行分析，檢驗(yàn)結(jié)果如表1，表2所示。對(duì)于來(lái)說(shuō)，耕地與林地、城鄉(xiāng)工礦居民用地，林地與草地、城鄉(xiāng)工礦居民用地，草地與城鄉(xiāng)工礦居民用地等五組土地利用類型下存在顯著性差異。對(duì)于來(lái)說(shuō)，耕地與林地、草地、城鄉(xiāng)工礦居民用地，林地與草地、城鄉(xiāng)工礦居民用地等五組土地利用類型下存在顯著性差異。

多重比較分析的結(jié)果進(jìn)一步證實(shí)，林地與耕地、草地的存在著顯著的差異，出現(xiàn)這種差異可能是由于林地與耕地、草地的葉面積指數(shù)()有所差異。葉面積指數(shù)是指植物所有葉片面積的總和與植物所覆蓋的地面面積的比值，能夠反映截獲的太陽(yáng)輻射、進(jìn)行氣體交換的潛在的葉片[20]。根據(jù)謝軍飛[21]的研究，林地最大值介于3.18～3.81之間，作物的最大值為1.44，草地的最大值為2.36，因此一般而言林地的值是要高于耕地和草地的。又由于與存在相關(guān)關(guān)系[22-23]，Baret[24]認(rèn)為兩者之間的函數(shù)關(guān)系可近似表示為：=1-exp（-K*），消光系數(shù)K取值范圍為0.7～1.4，因此林地和耕地、草地也出現(xiàn)顯著性差異。城鄉(xiāng)工礦居民用地與耕地、林地、草地的均存在著顯著性的差異，這是由于城市化過(guò)程中人口聚集，對(duì)大自然進(jìn)行森林采伐以及地面硬質(zhì)化等，導(dǎo)致城鄉(xiāng)工礦居民用地缺乏綠色植被[25]，所以城鄉(xiāng)工礦居民用地的最低。

不同土地利用類型的出現(xiàn)顯著差異也與有關(guān)。在沒(méi)有水分限制時(shí)，主要表示地表蒸散，包括蒸發(fā)（裸土以及植被冠層截獲水分的蒸發(fā)）和植物的蒸騰，很大程度上有蒸散的能量（潛熱）來(lái)決定[26]。顯熱和潛熱通過(guò)近地層大氣和下墊面能量和水分交換而產(chǎn)生[27]。在植被生長(zhǎng)旺盛的地方，和值高，植被生命活動(dòng)旺盛，蒸騰量大，蒸騰阻力小，潛熱占比增大，地表溫度降低[26]。由于林地的是要高于耕地、草地，且劉鳳山[28]的模擬結(jié)果顯示，森林(=2.23)>農(nóng)田(=1.04)>草地(=0.62)，三者的蒸騰量不同，潛熱占比不同，因此三者的也存在著顯著性的差異。城鄉(xiāng)工礦居民用地與耕地、林地的存在著較大的差異，但是與草地的差異性并不顯著。顯著是因?yàn)槌青l(xiāng)工礦居民用地為人工地表植被較少，其較耕地、林地等自然地表的低，植被的蒸騰量小，地表溫度高。與草地的差異不顯著，有可能是相鄰的城鄉(xiāng)工礦用地與草地進(jìn)行氣體交換、流通，形成局部氣體環(huán)流[29]，所以差異不明顯。

表 1 土地利用類型的NDVI多重比較

表2 土地利用類型的LST多重比較

*平均值差值的顯著性水平為0.0001

2.3 不同土地利用類型的LST與NDVI的定量關(guān)系

張家口市的與的線性回歸分析如圖7所示，總體來(lái)說(shuō)與呈明顯的負(fù)相關(guān)關(guān)系，判定系數(shù)2=0.519。岳文澤對(duì)于上海城市市區(qū)的和的研究，和存在著明顯的負(fù)相關(guān)關(guān)系，判定系數(shù)2=0.497，與本文的研究結(jié)論一致。對(duì)不同的土地利用類型的與的線性回歸分析如表3所示。水域的與為負(fù)相關(guān)關(guān)系，但是其線性關(guān)系不強(qiáng)烈，有研究表明，水體具有特殊性，其和具有正相關(guān)關(guān)系[30,31]，與本文研究結(jié)果相反，這可能是由于本文研究中水域的樣本數(shù)量較少（見(jiàn)表3），導(dǎo)致線性關(guān)系不明顯，水域的和的相關(guān)關(guān)系仍需要進(jìn)一步探究。除水域外，其它四種類型的土地利用的和呈現(xiàn)明顯負(fù)相關(guān)關(guān)系，且線性關(guān)系較顯著，相關(guān)的線性回歸函數(shù)如表3所示。一次項(xiàng)系數(shù)最大的是耕地，最小的是城鄉(xiāng)工礦居民用地，代表著隨著的增加，耕地的降低的最快、城鄉(xiāng)工礦居民用地的降低的最慢，耕地比城鄉(xiāng)工礦居民用地隨著的變化的敏感性比較高。

圖 7 NDVI與LST的線性回歸圖

表 3 不同土地利用類型的回歸函數(shù)分析

3 結(jié)論

本文基于Landsat8遙感影像數(shù)據(jù)和GIS空間分析技術(shù)，分析了張家口市區(qū)和的空間格局，范圍為-0.55～0.94，的范圍為1.02～59.9 ℃，且兩者的空間上值的分布呈相反趨勢(shì)。東西、南北兩條特征剖面線同樣顯示，各方向上的兩者均呈現(xiàn)相反的變化趨勢(shì)，增長(zhǎng)時(shí)減少，高時(shí)對(duì)應(yīng)的低，東西和南北方向的的最高值對(duì)應(yīng)的均是林地，最低值對(duì)應(yīng)的均是城鄉(xiāng)工礦居住用地，且各方向的和的折線圖關(guān)于水平方向?qū)ΨQ。

不同土地利用類型下的、的平均值有所差異，在所有類型中林地的平均值最大為0.65，耕地和草地的平均值均為0.46，城鄉(xiāng)工礦居民用地平均值最低為0.29；林地的最小為34.21 ℃，耕地和草地的也較接近。

不同的土地利用類型下、存在顯著性差異。對(duì)于而言，城鄉(xiāng)工礦居民用地與耕地、草地、林地等土地利用類型存在顯著性差異；對(duì)于而言，耕地與草地、林地、城鄉(xiāng)工礦居民用地等土地利用類型存在顯著性差異。研究區(qū)域除水域之外的與的回歸函數(shù)為顯著的負(fù)相關(guān)關(guān)系。不同的土地利用類型的負(fù)相關(guān)關(guān)系的相關(guān)程度不同，不同的土地利用類型的對(duì)于的影響效果不同。其中，耕地的變化對(duì)影響較大，城鄉(xiāng)工礦居民用地的對(duì)影響比較小。本文研究成果能夠?yàn)榫徑獬鞘袩釐u效應(yīng)提供理論依據(jù)，可以為張家口市的城市生態(tài)規(guī)劃提供科學(xué)依據(jù)。

[1] Wu J, Jelinski DE, Luck M,. Multiscale analysis of landscape heterogeneity: scale variance and pattern metrics[J]. Annals of GIS, 2000,6(1):6-19

[2] 湯君友楊山趙銳.無(wú)錫市城鎮(zhèn)用地變化及其環(huán)境效應(yīng)研究[J].國(guó)土資源遙感,2002(3):16-18

[3] 彭璐,李小梅,傅佳麗,等.福建省平潭島植被覆蓋變化及其對(duì)地表溫度影響的定量研究[J].福建師范大學(xué)學(xué)報(bào)(自 然科學(xué)版),2016,32(1):78-85

[4] 高尚,沙晉明,帥晨.廈門市地表溫度與植被覆蓋關(guān)系定量研究[J].福建師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2019,35(2):14-21

[5] 姚遠(yuǎn),陳曦,錢靜.城市地表熱環(huán)境研究進(jìn)展[J].生態(tài)學(xué)報(bào),2018,38(3):1134-1147

[6] 牟雪潔,趙昕奕.珠三角地區(qū)地表溫度與土地利用類型關(guān)系[J].地理研究,2012,31(9):1589-1597

[7] 彭文甫,周介銘,羅懷良,等.城市土地利用與地面熱效應(yīng)時(shí)空變化特征的關(guān)系——以成都市為例[J].自然資源學(xué)報(bào),2011,26(10):1738-1749

[8] 劉宇,匡耀求,吳志峰,等.不同土地利用類型對(duì)城市地表溫度的影響——以廣東東莞為例[J].地理科學(xué),2006(5):5597-5602

[9] Wilson JS, Clay M, Martin E,. Evaluating environmental influences of zoning in urban ecosystems with remote sensing[J]. Remote Sensing of Environment, 2003,86(3):303-321

[10] 岳文澤,徐建華,徐麗華.基于遙感影像的城市土地利用生態(tài)環(huán)境效應(yīng)研究——以城市熱環(huán)境和植被指數(shù)為例[J]. 生態(tài)學(xué)報(bào),2006,26(5):1450-1460

[11] Lo CP, Quattrochi DA, Luvall JC. Application of high-resolution thermal infrared remote sensing and GIS to assess the urban heat island effect[J]. International Journal of Remote Sensing, 2010,18(2):287-304

[12] 周媛,石鐵矛,胡遠(yuǎn)滿,等.基于城市土地利用類型的地表溫度與植被指數(shù)的關(guān)系[J].生態(tài)學(xué)雜志,2011,30(7):1504-1512

[13] 胡德勇,喬琨,王興玲,等.單窗算法結(jié)合Landsat8熱紅外數(shù)據(jù)反演地表溫度[J].遙感學(xué)報(bào),2015,19(6):964-976

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

[15] 白潔,劉紹民,扈光.針對(duì)TM/ETM+遙感數(shù)據(jù)的地表溫度反演與驗(yàn)證[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(9):148-155

[16] 魏寶成,銀山,宋潔,等.呼和浩特市不同植被指數(shù)與地表溫度的定量遙感關(guān)系[J].水土保持研究,2015,22(5):79-85

[17] 張曉娟,周啟剛,黃麗盺,等.重慶市主城區(qū)地表溫度與植被覆蓋指數(shù)關(guān)系研究[J].土壤通報(bào),2018,49(2):293-302

[18] Weng QH, Lu DS, Schubring J. Estimation of land surface temperature–vegetation abundance relationship[J]. Remote Sensing of Environment, 2004,89(4):467-483

[19] 蔣晶,喬治.北京市土地利用變化對(duì)地表溫度的影響分析[J].遙感信息,2012,27(3):105-111

[20] 陳拉,黃敬峰,王秀珍.不同傳感器的模擬植被指數(shù)對(duì)水稻葉面積指數(shù)的估測(cè)精度和敏感性分析[J].遙感學(xué)報(bào),2008,12(1):143-151

[21] 謝軍飛,許蕊.基于MODIS數(shù)據(jù)的北京植被葉面積指數(shù)的時(shí)空變化[J].生態(tài)學(xué)雜志,2014,33(5):1374-1380

[22] Carlson TN, Ripley DA. On the relation between, fractional vegetation cover, and leaf area index[J]. Remote Sensing of Environment, 1997,62(3):241-252

[23] Nemani RR, Running SW. Testing a theoretical climate-soil-leaf area hydrologic equilibrium of forests using satellite data[J]. Agricultural and Forest Meteorology, 1989(44):245-260

[24] Baret F, Guyot G. Potentials and limits of vegetation indices for LAI and APAR assessment[J]. Remote Sensing of Environment, 1991,35(2-3):161-173

[25] Owen TW, Carlson TN, Gillies RR. An assessment of satellite remotely-sensed land cover parameters in quantitatively describing the climaticeffect of urbanization[J]. Remote Sensing, 1998,m19(9):1663-1681

[26] 韓麗娟,王鵬新,王錦地,等.植被指數(shù)_地表溫度構(gòu)成的特征空間研究[J].中國(guó)科學(xué)(D輯:地球科學(xué)),2005(4):371-377

[27] 賀有為,王秋兵,溫學(xué)發(fā),等.季節(jié)性干旱對(duì)中亞熱帶人工林顯熱和潛熱通量日變化的影響[J].生態(tài)學(xué)報(bào),2011,31(11):3069-3081

[28] 劉鳳山,陶福祿,肖登攀,等.土地利用類型轉(zhuǎn)換對(duì)地表能量平衡和氣候的影響——基于SiB2模型的模擬結(jié)果[J].地 理科學(xué)進(jìn)展,2014,33(6):815-824

[29] 王修信,胡玉梅,劉馨,等.城市草地的小氣候調(diào)節(jié)作用初步研究[J].廣西師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2007(3):23-27

[30] 唐曦,束炯,樂(lè)群.基于遙感的上海城市熱島效應(yīng)與植被的關(guān)系研究[J].華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2008(1):119-128

[31] Yue W, Xu J, Tan W,. The relationship between land surface temperature andwith remote sensing: application to Shanghai Landsat 7 ETM+ data[J]. International Journal of Remote Sensing, 2007,28(15):3205-3226

Relationships between Normalized Vegetation Index and Land Surface Temperature Based on Different Land Use Types——A Case of Zhangjiakou City

CAO Qian-qian, LIU Rui-fen*

430068,

Based on Landsat 8 remote sensing data and ArcGIS spatial analysis technology, the spatial normalized vegetation index () and surface temperature () were studied for Zhangjiakou City in Hebei Province. Different values ofandfor six land use types were characterized and compared by the multiple comparisons method. Finally, the relationship betweenandof different land use types was analyzed by the linear regression method. The results show that the spatial distributions ofandin Zhangjiakou City have opposite trends. For, there are significant differences among urban-rural construction land, cultivated land, grassland, woodland and other land use types. For, there are significant differences among cultivated land, grassland, forest land, and urban-rural construction land. The negative relationships exist betweenandunder different land use types, and the determination coefficient2=0.519. The study provides a scientific basis for urban ecological planning and heat island alleviation in Zhangjiakou City.

Land use; land surface temperature; normalized vegetation index

S181

A

1000-2324(2021)06-0955-09

2021-10-26

2021-11-12

國(guó)家自然科學(xué)基金青年科學(xué)基金資助項(xiàng)目:基于優(yōu)先流控制的綠色屋頂植物-人工基質(zhì)結(jié)構(gòu)優(yōu)化設(shè)計(jì)(51909081);國(guó)家科技重大專項(xiàng)水專項(xiàng)項(xiàng)目(2017ZX07101003-08)

曹倩倩(1997-),女,碩士研究生,主要研究方向:低影響開(kāi)發(fā)技術(shù)及城市規(guī)劃. E-mail:1141617706@qq.com

通訊作者:Author for correspondence. E-mail:ruifen1986@aliyun.com