基于網(wǎng)格單元的烏魯木齊市土地覆被/利用時(shí)空變化

買(mǎi)買(mǎi)提江·買(mǎi)提尼亞孜，阿里木江·卡斯木

（1. 新疆農(nóng)業(yè)大學(xué)管理學(xué)院，烏魯木齊 830052；2. 新疆師范大學(xué)地理科學(xué)與旅游學(xué)院，烏魯木齊 830054；3. 新疆師范大學(xué)絲綢之路經(jīng)濟(jì)帶城鎮(zhèn)化發(fā)展研究中心，烏魯木齊 830054）

0 引 言

隨著城市化進(jìn)程在全球范圍內(nèi)的快速推進(jìn)，全世界已有50%的人口居住在城市，而且這種趨勢(shì)還將繼續(xù)[1]。伴隨著城市化的發(fā)展，城區(qū)土地覆被/土地利用也發(fā)生了變化，這對(duì)生態(tài)環(huán)境、氣候以及社會(huì)經(jīng)濟(jì)等造成了一定的影響[2-5]。因而，有關(guān)城區(qū)土地覆被/土地利用動(dòng)態(tài)及相關(guān)問(wèn)題的研究已成為一大社會(huì)焦點(diǎn)。近年來(lái)，遙感、空間統(tǒng)計(jì)及地理信息技術(shù)已被用于土地方面的研究，這些技術(shù)的日益成熟推動(dòng)土地動(dòng)態(tài)等相關(guān)研究不斷深入[6-8]。

學(xué)者們對(duì)城市土地覆被/利用變化（land use and land cover change，LUCC）帶來(lái)的一系列問(wèn)題進(jìn)行了大量的研究，涉及城市LUCC的時(shí)空演變[9-10]、驅(qū)動(dòng)力分析[11-12]、情景模擬和預(yù)測(cè)[13-14]以及生態(tài)環(huán)境效應(yīng)等[3,15]。當(dāng)前對(duì)城市LUCC方面的研究大多是以分析研究區(qū)內(nèi)的整體動(dòng)態(tài)變化為主[10]，而對(duì)其內(nèi)部各子區(qū)域土地利用時(shí)空變化的定量分析相對(duì)缺乏[6-7]。在進(jìn)行關(guān)于城市LUCC驅(qū)動(dòng)力研究時(shí)，主要進(jìn)行LUCC總量與各自然人文因子的總量間的相關(guān)分析，以相關(guān)系數(shù)較高的為主要驅(qū)動(dòng)力[11-12]，這雖然可以解釋土地總體變化的原因，但卻無(wú)法定量分析研究區(qū)內(nèi)不同區(qū)域LUCC的空間異質(zhì)性[7]。網(wǎng)格單元（grid cell，GC）方法可以定量分析 LUCC在時(shí)間和空間上的總體及局部的精細(xì)變化，并且能夠通過(guò)將相關(guān)的各類自然及人文因子轉(zhuǎn)入網(wǎng)格單元內(nèi)，定量分析研究區(qū)內(nèi)不同區(qū)域LUCC的驅(qū)動(dòng)力[6-7]。目前，中國(guó)學(xué)者還未將網(wǎng)格單元法用于研究LUCC，國(guó)外學(xué)者Bagan等[6]在研究東京城市增長(zhǎng)及其驅(qū)動(dòng)力時(shí)引入GC分析了LUCC的時(shí)空變化，其后 Maimaitijiang等[7]在一定程度上對(duì)此方法進(jìn)行了改進(jìn)，并在LUCC時(shí)空變化及驅(qū)動(dòng)力研究中進(jìn)行了應(yīng)用?；谵D(zhuǎn)移矩陣與變化檢測(cè)(change detection)的方法可以定性分析研究區(qū)內(nèi)某區(qū)域各種用地類型的增減情況，但不能定量表征其變化情況，而基于GC的方法則可定量分析研究區(qū)內(nèi)任一區(qū)域土地利用的增減情況，這種特性為探討研究區(qū)內(nèi)不同區(qū)域和方向上各類用地的變化及其與相關(guān)自然人文因子空間關(guān)聯(lián)分析奠定了基礎(chǔ)，而轉(zhuǎn)移矩陣與變化檢測(cè)的方法則無(wú)法實(shí)現(xiàn)[7]。

本文引入GC的思路和方法，基于遙感數(shù)據(jù)，對(duì)烏魯木齊市1990年以來(lái)LUCC的時(shí)空變化進(jìn)行定量分析，以期為當(dāng)前烏魯木齊城市發(fā)展提供借鑒，同時(shí)還可更好地為城市未來(lái)的發(fā)展進(jìn)行預(yù)測(cè)。

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

1.1 研究區(qū)概況

烏魯木齊市地處中國(guó)西北新疆中部地區(qū)的天山中段北麓、準(zhǔn)噶爾盆地南緣，位于歐亞大陸腹地，是世界上距離海洋最遠(yuǎn)的內(nèi)陸城市。烏魯木齊屬中溫帶半干旱區(qū)，受高壓控制和地形作用影響，四季都會(huì)出現(xiàn)地形逆溫層和下沉逆溫層，是典型的大陸性氣候特征。烏魯木齊市是新疆的經(jīng)濟(jì)、政治、文化、科教、金融和交通運(yùn)輸中心，是第 2座亞歐大陸橋中國(guó)西部橋頭堡和向西開(kāi)放的重要門(mén)戶，是中國(guó)“一帶一路”戰(zhàn)略絲綢之路經(jīng)濟(jì)帶核心區(qū)的首府城市，近年來(lái)隨著烏魯木齊市經(jīng)濟(jì)迅猛發(fā)展，城市化進(jìn)程不斷加快，土地利用在深度和廣度方面發(fā)生了顯著的變化[16]。本文選擇LUCC變化較大，人文干擾強(qiáng)度高的城市建成區(qū)（市區(qū)）為研究區(qū)（圖1），市區(qū)是一個(gè)形似狹長(zhǎng)葫蘆狀的河谷盆地，建成區(qū)東、西、南3面環(huán)山，地勢(shì)東南高、西北低，研究區(qū)南北長(zhǎng)約48 km，東西寬約36 km，面積為1 716.68 km2。

圖1 研究區(qū)Landsat TM影像圖Fig.1 Landsat TM image of study area

1.2 數(shù)據(jù)來(lái)源

本研究采用5幅Landsat TM5影像：1）1990年1幅(成像時(shí)間為8月6日，軌道號(hào)為142/30)；2）因2000年和2011年1幅影像不能完全涵蓋研究區(qū)，因此每期選擇2幅影像，即2000年2幅（成像時(shí)間均為8月8日，軌道號(hào)分別為143/29、143/30）、2011年2幅（成像時(shí)間為8月11日，軌道號(hào)分別為143/29、143/30）。5幅影像均為同一月份數(shù)據(jù)，可在一定程度上減少或剔除LUCC季節(jié)性變化的影響，同時(shí)數(shù)據(jù)均來(lái)自Landsat TM5傳感器，減少因傳感器不同而造成的誤差[17]。數(shù)據(jù)處理過(guò)程中，首先以烏魯木齊市l(wèi)∶100 000地形圖為基準(zhǔn)來(lái)進(jìn)行地面控制點(diǎn)的選取，再將研究區(qū)影像幾何校正，將誤差控制在l個(gè)像元內(nèi)。本文選取應(yīng)用廣泛且不受區(qū)域特點(diǎn)及目標(biāo)類型影響的成像時(shí)間均為同天的 6S大氣校正模型，其成像時(shí)晴朗無(wú)云，可消除大氣對(duì)影像產(chǎn)生的影響。

2 研究方法

2.1 土地覆被/土地利用分類

城市區(qū)域“異物同譜，同物異譜”現(xiàn)象較普遍（例如城市周邊耕地和城市草地，以及建設(shè)用地和裸地等具有相似的光譜特性），因此基于低、中分辨率的多光譜遙感影像對(duì)城市用地信息進(jìn)行提取和分類時(shí)會(huì)有一定難度[18]，尤其是在干旱區(qū)城市，其裸地分布較為廣泛，建設(shè)用地和裸地在影像分類結(jié)果中混雜的現(xiàn)象較普遍[19]。本文不僅要研究土地利用的總體變化，更側(cè)重土地利用在研究區(qū)內(nèi)部不同區(qū)域空間的精細(xì)變化，因此對(duì)用地分類的精度要求非常高[6-7,20]。本文在充分考慮烏魯木齊市土地利用特征以及各用地類型在影像上的區(qū)分度的基礎(chǔ)上，將11個(gè)土地利用類型合并為城市建設(shè)用地、綠地、水體和裸地（包括沒(méi)有綠色作物的耕地）4個(gè)一級(jí)土地利用單元[21-22]（表 1），以通過(guò)減少類別數(shù)量來(lái)提高分類的精度[17]。

表1 烏魯木齊市土地覆被/利用分類Table 1 Land cover/land use classification of Urumqi city

為了提高信息提取的精度，通過(guò)ENVI5.1進(jìn)行波段運(yùn)算，分別計(jì)算表征植被信息的土壤調(diào)整植被指數(shù)（soil adjusted vegetation index，SAVI）[23]、表征建筑用地信息的歸一化建筑指數(shù)（normalized difference built-up index，NDBI）[24]和表征水體的改進(jìn)的歸一化差異水體指數(shù)（modified normalized differences water index，MNDWI）[25]。本文將這3個(gè)指數(shù)作為影像R、G、B 3個(gè)波段進(jìn)行疊加，采用監(jiān)督分類方法中的最大似然分類法對(duì)影像進(jìn)行分類，經(jīng)精度驗(yàn)證，得知用戶精度和生產(chǎn)者精度均達(dá)到85%以上，符合分類精度 82%的要求[26]，其中歷史結(jié)果的精讀評(píng)價(jià)首先根據(jù)當(dāng)年的烏魯木齊市土地利用圖和目視解譯在ENVI里通過(guò)confusion matrix工具對(duì)感興趣區(qū)隨機(jī)采樣（樣本數(shù)量為300）進(jìn)行精度驗(yàn)證，然后根據(jù)實(shí)地考察結(jié)果對(duì)最終分類結(jié)果進(jìn)行精度分析，3期的Kappa系數(shù)均達(dá)到0.9以上。

2.2 網(wǎng)格單元方法

基于ArcGIS平臺(tái)，在研究區(qū)邊界范圍內(nèi)建立矢量網(wǎng)格。根據(jù)研究的需要確定網(wǎng)格單元大小，將基于遙感數(shù)據(jù)分類獲得的專題圖通過(guò)ArcGIS軟件的連接功能分配轉(zhuǎn)入對(duì)應(yīng)位置的網(wǎng)格單元內(nèi)，并計(jì)算每個(gè)網(wǎng)格單元內(nèi)各類用地占單元總面積的百分比（圖 2）。在向每個(gè)網(wǎng)格單元分配轉(zhuǎn)入各類用地或其他空間要素的過(guò)程中，存在跨越網(wǎng)格單元邊界的像元如何分配的問(wèn)題（圖 2）。傳統(tǒng)的方法是優(yōu)先分配給相應(yīng)的所占面積大的網(wǎng)格單元[6]，這樣會(huì)造成分配到網(wǎng)格單元的用地面積之和大于網(wǎng)格單元自身面積，不符合實(shí)際情況[7]。本文將遙感數(shù)據(jù)分類獲得的專題圖轉(zhuǎn)換為矢量圖，與已建立的矢量網(wǎng)格相交，使跨越網(wǎng)格邊界的斑塊被網(wǎng)格分割成多個(gè)小塊，計(jì)算每個(gè)網(wǎng)格內(nèi)每類用地被分割后的小塊面積，并求和，然后計(jì)算該用地在相應(yīng)網(wǎng)格內(nèi)所占的比例：

式中Pt為用地類型t在每個(gè)網(wǎng)格內(nèi)的面積比例，%；At為每個(gè)網(wǎng)格內(nèi)用地類型t被分割后的每個(gè)小塊的面積；n為每個(gè)網(wǎng)格內(nèi)用地類型t被分割后的數(shù)量；Ag是每個(gè)網(wǎng)格的面積。

綜合考慮研究區(qū)面積、研究目的以及計(jì)算機(jī)運(yùn)行效率等因素，通過(guò)反復(fù)試驗(yàn)，將每個(gè)矢量網(wǎng)格單元的大小確定為250 m×250 m，總計(jì)27 743個(gè)網(wǎng)格單元。

圖2 土地覆被/利用數(shù)據(jù)格網(wǎng)化過(guò)程示意圖Fig.2 Diagram of land use and land cover data gridding procedure

3 結(jié)果與分析

3.1 土地覆被/利用總體變化分析

影像分類結(jié)果及基于分類結(jié)果統(tǒng)計(jì)獲得各類用地面積及比例變化如圖3和表2所示。由圖3可知，1990年以來(lái)裸地減少、建設(shè)用地呈增加趨勢(shì)。與2000年相比，2011年綠地和水體明顯增加。從表中可以看出，1990—2011年烏魯木齊市建設(shè)用地增加較快，從1990年157.37 km2增加到2011年的444.89 km2，比例也由9.17%增加到25.92%，這可能與過(guò)去20多年間城市人口的快速增加和社會(huì)經(jīng)濟(jì)的快速發(fā)展具有緊密的聯(lián)系；其次是裸地減少較多，比例由1990年的75.13%減少到2011年的49.79%；綠地面積變化呈現(xiàn)一定的波動(dòng)，先減少后增加，需要指出的是，由于耕地上有綠色作物時(shí)，在影像上歸為綠地，而沒(méi)有綠色作物時(shí)則歸為裸地，且城市北部（實(shí)際以耕地為主）在2011年綠地增加較多，說(shuō)明除了降水等自然因素影響外，耕地的作物種植或收割情況等人為因素也有重要的影響；最后，研究區(qū)內(nèi)水體的面積呈緩慢增加的趨勢(shì)，比例由0.71%增加到0.84%。

3.2 基于網(wǎng)格單元的各用地類型時(shí)空變化分析

基于網(wǎng)格單元法建立了每類用地變化情況的矢量網(wǎng)格圖（每個(gè)網(wǎng)格大小為250 m×250 m），每個(gè)網(wǎng)格的屬性值表示1990—2011年相應(yīng)用地類型占相應(yīng)網(wǎng)格面積比例的變化情況?；诰W(wǎng)格單元法的建設(shè)用地1990—2011年的時(shí)空變化圖4所示。

圖4表明，過(guò)去20 a間，除極少數(shù)區(qū)域外，建設(shè)用地在所有區(qū)域幾乎都表現(xiàn)出增加的趨勢(shì)，主要是在 1990年建成區(qū)的基礎(chǔ)上向東北、北面、西北以及東南方向蔓延發(fā)展，同時(shí)在過(guò)去的20多年間，建設(shè)用地主要以外延式擴(kuò)張為主，但也有內(nèi)部挖潛的集約式增長(zhǎng)；而少數(shù)區(qū)域如老城區(qū)建設(shè)用地減少可能是由于對(duì)過(guò)去高密度“城中村、棚戶區(qū)”改造為“低密度、高綠地率”的建筑區(qū)等因素引起。

圖3 1990—2011年烏魯木齊市土地覆被/利用變化Fig.3 Land use and land cover change of Urumqi city from year of 1990 to 2011

表2 烏魯木齊市1990—2011年土地覆被/利用面積及比例Table 2 Area and percentage of land use and land cover of Urumqi city from year of 1990 to 2011

從圖 4中還可以看出，雖然建設(shè)用地的擴(kuò)展在不同方向、不同區(qū)域有很大的不同，這可能由很多因素綜合作用形成。1）從自然條件來(lái)看，烏魯木齊市主城區(qū)地勢(shì)高、坡度大且東、西、南3面環(huán)山（圖1），不利于作為建設(shè)用地；2）從政策方面來(lái)看，近幾年國(guó)家和政府出臺(tái)并實(shí)施的區(qū)域一體化戰(zhàn)略措施，如 “烏昌一體化”進(jìn)程（即烏魯木齊市與西北方向的昌吉市一體化發(fā)展）、“烏昌石”城市圈、“兵地融合”發(fā)展等[27]，促進(jìn)了烏魯木齊與西北方向的昌吉市、石河子市，與北部的兵團(tuán)城市五家渠市、東北方向的阜康市等的聯(lián)系和融合發(fā)展；行政區(qū)劃調(diào)整，例如2007年位于烏魯木齊東北方向的米泉市并入首府成為米東區(qū)，同時(shí)推動(dòng)建設(shè)用地向東北方向的擴(kuò)展；同時(shí)，烏魯木齊市在近幾年規(guī)劃中還提出了“南控北擴(kuò)、西延?xùn)|進(jìn)”發(fā)展戰(zhàn)略，城市用地實(shí)行“南控、北擴(kuò)、先兩延、后東進(jìn)”的原則，這在一定程度上決定了烏魯木齊城市建設(shè)用地的發(fā)展趨勢(shì)和方向。3）除了自然地形、政策規(guī)劃等因素外，城市主要路網(wǎng)和對(duì)外交通對(duì)建設(shè)用地在不同方向的增長(zhǎng)和擴(kuò)展具有重要影響，城市道路網(wǎng)絡(luò)是一個(gè)城市在地域空間上的“骨骼”、基礎(chǔ)支撐以及城市用地?cái)U(kuò)張的一項(xiàng)重要驅(qū)動(dòng)因素，與城市用地?cái)U(kuò)張具有密切聯(lián)系[28-29]。烏魯木齊市1990年道路長(zhǎng)度為402 km，2011年達(dá)到了1 695 km，是1990年的4.22倍[30]。將城市道路系統(tǒng)圖和城市建設(shè)用地圖對(duì)比發(fā)現(xiàn)，建設(shè)用地的增加區(qū)域與主干路網(wǎng)在空間具有一定的一致性，主干路對(duì)城市擴(kuò)展的“牽引”特性明顯，尤其是城市對(duì)外交通以及外環(huán)路沿線建設(shè)用地增加明顯，如吐烏大高速路、S111省道、烏昌快速、連霍高速、G314和G216國(guó)道等對(duì)城市建設(shè)用地空間向東北、西北和東南方向延展起到了一定的推動(dòng)和促進(jìn)作用。

圖4 1990—2011年烏魯木齊市建設(shè)用地、綠地和裸地的時(shí)空 變化與道路系統(tǒng)Fig.4 Spatial and temporal change in construction land, green land and bare land of Urumqi city from 1990 to 2011 and road system

過(guò)去20 a間，在建成區(qū)內(nèi)外綠地既有增加的區(qū)域，但也有不少減少的區(qū)域，如圖4c所示。從建成區(qū)范圍內(nèi)來(lái)看，綠地減少的區(qū)域較多，區(qū)域⑧除自然因素外，主要是由建設(shè)用地?cái)U(kuò)展占用草地和林地所致；但建成區(qū)內(nèi)部有很多區(qū)域表現(xiàn)出綠地的增加，而且有幾處為大面積明顯的增加，這與烏魯木齊市近年來(lái)一直實(shí)施的“園林城市建設(shè)、荒山綠化工程”戰(zhàn)略密切相關(guān)，烏魯木齊市在1990年僅有4個(gè)公園，到2011年增加到了26個(gè)，面積從1990年的191 hm2增加到2011年的2 636 hm2，建成區(qū)綠化率從21.8%提高到36.16%[30]，這使得區(qū)域內(nèi)尤其是建成區(qū)范圍內(nèi)很多區(qū)域綠地面積大規(guī)模增加，如圖4c中綠地增加明顯的區(qū)域：①雅瑪里克山、②水磨溝山區(qū)、③蜘蛛山、④紅光山等均為烏魯木齊市重點(diǎn)綠化和改造的區(qū)域。2）從建成區(qū)外部的區(qū)域來(lái)看，研究區(qū)北部（區(qū)域⑤）、西北部（區(qū)域⑤）以及西南部（區(qū)域⑥）主要為耕地，這些區(qū)域均表現(xiàn)出綠地的明顯增減的特性，除了自然因素的影響外，1990年和2011年2期耕地上是否有作物或者被收割的狀態(tài)（沒(méi)有作物或被收割后耕地在影像上表現(xiàn)為裸地，有綠色作物時(shí)表現(xiàn)為綠地）對(duì)這些區(qū)域綠地的增減具有重要影響；研究區(qū)東南部（區(qū)域⑦）綠地表現(xiàn)出明顯增加的特性，這些區(qū)域海拔較高，主要為覆蓋稀疏植被或無(wú)植被的山體，且尚未進(jìn)行人為綠化，即人為干擾的作用非常微弱，因此可以判定該區(qū)域綠地增多主要為不同時(shí)期自然因素如降水和氣溫的不同造成（圖5）。

裸地是干旱區(qū)城市地表覆被的重要特征。圖4d顯示了基于網(wǎng)格單元分析的裸地在1990—2011年的時(shí)空變化情況，可以看出，過(guò)去20 a間，建成區(qū)內(nèi)外裸地既有增加的區(qū)域，但也有不少減少的區(qū)域。在1990—2011年建成區(qū)之間的區(qū)域顯示裸地明顯減少的特性，說(shuō)明過(guò)去20a間城區(qū)向外圍擴(kuò)展的過(guò)程中，占用了大量的裸地（建設(shè)用地占用或綠化后成為綠地）；從 2011年建成區(qū)外部來(lái)看，研究區(qū)北部、西北部以及西南部裸地的增減變化特性突出，因?yàn)樵搮^(qū)域主要為耕地，前后 2期耕地上是否有作物或者被收割的狀態(tài)決定這些區(qū)域裸地的增減變化；2011年建成區(qū)外研究區(qū)西南部表現(xiàn)裸地明顯的減少特征，而該區(qū)域是受人為影響較弱的覆被稀疏植被或無(wú)植被的自然山體，因此該區(qū)域裸地的減少主要是由于前后 2期自然因素如降水和氣溫變化引起植被覆蓋變化，從而引起裸地的變化。

圖5 1990—2011年烏魯木齊年均氣溫和降水情況Fig.5 Annual average temperature and precipitation of Urumqi city from 1990 to 2011

3.3 基于網(wǎng)格單元的各類用地變化相關(guān)分析

基于轉(zhuǎn)移矩陣分析不同用地類型的時(shí)空轉(zhuǎn)移和替代關(guān)系是土地覆被/土地利用研究中的重要方法[31-33]，而基于網(wǎng)格單元的可以通過(guò)計(jì)算不同用地類型變化之間的相關(guān)系數(shù)，分析不同用地類型的時(shí)空轉(zhuǎn)移和替代關(guān)系[6-7]。由于研究區(qū)范圍內(nèi)有大量耕地，其耕種狀態(tài)會(huì)影響綠地和裸地的面積，為了更加準(zhǔn)確地反映20 a間各用地類型的相互轉(zhuǎn)移替代關(guān)系，本文基于2011年建成區(qū)范圍界線內(nèi)的10 236個(gè)網(wǎng)格單元計(jì)算各用地類型之間的變化替代關(guān)系（表3），其中r為皮爾遜相關(guān)系數(shù)，P≤0.05，具有95%置信區(qū)間的統(tǒng)計(jì)顯著性。從表3可以看出，3個(gè)變化階段建設(shè)用地與裸地變化均具有很強(qiáng)的負(fù)相關(guān)性（r約-0.7），同時(shí)其與綠地的變化則具有一定的負(fù)相關(guān)性（r為-0.23左右），而建設(shè)用地與水體相關(guān)性十分微弱，說(shuō)明建設(shè)用地增加和擴(kuò)展主要占用和替代了裸地，其次也占用了一定的綠地；同時(shí)，還發(fā)現(xiàn)在 3個(gè)變化階段裸地與綠地變化均具有較強(qiáng)的負(fù)相關(guān)性（相關(guān)系數(shù)為-0.5左右），說(shuō)明裸地減少部分被綠地所代替，也再次印證烏魯木齊市長(zhǎng)期以來(lái)的“荒山綠化”工程的效果。

表3 1990—2011年各類用地變化量之間的相關(guān)性分析Table 3 Correlation analysis between change of each land use in different stages during 1990-2011

為了進(jìn)一步分析20 a間建設(shè)用地與裸地的轉(zhuǎn)移替代關(guān)系，對(duì)兩者的變化進(jìn)行線性回歸分析（圖 6），說(shuō)明過(guò)去20 a間烏魯木齊市裸地的大規(guī)模減少主要是由建設(shè)用地的增加和擴(kuò)展所驅(qū)動(dòng)。

圖6 1990—2011年烏魯木齊市建設(shè)用地和裸地時(shí)空變化的回歸分析Fig.6 Regression analysis on change in construction land area and that of bare land area in Urumqi city during 1990-2011

4 結(jié) 論

本文引入網(wǎng)格單元的思路和方法，基于遙感數(shù)據(jù)，對(duì)烏魯木齊市 1990—2011年土地覆被/土地利用的時(shí)空變化進(jìn)行定量分析，主要結(jié)論如下：

1）烏魯木齊市在過(guò)去20 a間城市快速擴(kuò)展，建設(shè)用地迅速增加，由 1990年 157.37 km2增加到 2011年的444.89 km2，與城市人口和社會(huì)經(jīng)濟(jì)增長(zhǎng)呈現(xiàn)一致趨勢(shì)。烏魯木齊市地形條件、政策、城市規(guī)劃以及道路交通網(wǎng)絡(luò)等因素對(duì)建設(shè)用地在不同區(qū)域和方向的擴(kuò)展起到了重要作用，在空間上既有在老城區(qū)內(nèi)部挖潛的集約式增長(zhǎng)，也伴隨向外圍的蔓延式擴(kuò)展；水體則呈現(xiàn)緩慢增長(zhǎng)的趨勢(shì)。

2）綠地在20 a間總體呈現(xiàn)增加趨勢(shì)，除降水和氣溫等自然因素影響外，政府長(zhǎng)期實(shí)施的“荒山綠化”工程對(duì)綠地的變化影響較大。

3）裸地始終向著持續(xù)減少的趨勢(shì)發(fā)展，比例由1990年的75.13%減少到2011年的49.79%而城市建設(shè)用地?cái)U(kuò)展占用和“荒山綠化”是導(dǎo)致其變化的關(guān)鍵因素。

本文在利用遙感技術(shù)的基礎(chǔ)上引入了網(wǎng)格的方法，使其結(jié)果更加精細(xì)和準(zhǔn)確的反映出近年來(lái)烏魯木齊市的土地利用動(dòng)態(tài)變化?；诰W(wǎng)格單元的數(shù)據(jù)可更好地為當(dāng)前烏魯木齊城市發(fā)展提供借鑒，同時(shí)還可更好地對(duì)城市未來(lái)發(fā)展進(jìn)行預(yù)測(cè)。

[1] Dept of International Economic and Social Affairs. World Urbanization Prospects: The 2014 Revision [M]. New York:United Nations, 2014.

[2] 李曉文，方創(chuàng)琳，黃金川，等. 西北干旱區(qū)城市土地利用變化及其區(qū)域生態(tài)環(huán)境效應(yīng)：以甘肅河西地區(qū)為例[J]. 第四紀(jì)研究，2003，23(3)：280－290.Li Xiaowen, Fang Chuanglin, Huang Jinchuan, et al. 2003.The urban land use transformation and associated effects of eco-environment in Northwest China arid region: A case study of Hexi region, Gansu province[J]. Quaternary Sciences, 2003, 23(3): 280－290. (in Chinese with English abstract)

[3] 何玲，賈啟建，李超，等. 基于生態(tài)系統(tǒng)服務(wù)價(jià)值和生態(tài)安全格局的土地利用格局模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32(3)：275－284.He Ling, Jia Qijian, Li Chao, et al. Land use pattern simulation based on ecosystem service value and ecological security pattern[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(3): 275－284. (in Chinese with English abstract)

[4] 黃迎春，楊伯鋼，張飛舟，等. 基于同類城市建設(shè)目標(biāo)的北京市土地利用結(jié)構(gòu)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32(4)：217－227.Huang Yingchun, Yang Bogang, Zhang Feizhou, et al.Optimization of land use structure in Beijing based on goal of similar city in world [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2016, 32(4): 217－227. (in Chinese with English abstract)

[5] De Freitas M W D, Santos J R D, Alves D S, et al. Land use and land cover change processes in the Upper Uruguay Basin:Linking environmental and socioeconomic variables[J].Landscape Ecology, 2013, 28(2): 311－327.

[6] Bagan H, Yamagata Y. Landsat analysis of urban growth:How Tokyo became the world’s largest megacity during the last 40 years [J]. Remote Sensing of Environment, 2012, 127:210－222.

[7] Maimaitijiang M, Ghulam A, Sandoval J S O. et al. Drivers of land cover and land use changes in St. Louis metropolitan area over the past 40 years characterized by remote sensing and census population data[J]. International Journal of Applied Earth Observation & Geoinformation, 2015, 35(3):161－174.

[8] 阿里木江·卡斯木，唐兵. 基于RS&GIS的新疆城市擴(kuò)展時(shí)空動(dòng)態(tài)變化分析[J]. 冰川凍土，2013，35(4)：1056–1064.Alimujiang·Kasimu, Tang Bing. Analysis of the spatialtemporal dynamic changes of urban expansion in oasis of Xinjiang based on RS and GIS[J]. Journal of Glaciology and Geocryology, 2013, 35(4): 1056－1064. (in Chinese with English abstract)

[9] 茹克亞·薩吾提，阿里木江·卡斯木. 近 20a來(lái)新疆主要城市擴(kuò)展時(shí)空動(dòng)態(tài)變化特征[J]. 干旱區(qū)研究，2015，32(3)：606–613.Rukeya·Sawuti, Alimujiang·Kasimu. Spatial-temporal dynamic characteristics of major cities expansion in Xinjiang in 20 Years [J]. Arid Zone Research, 2015, 32(3): 606－613. (in Chinese with English abstract)

[10] 王冬辰，杜培軍，蘇紅軍，等. 近 20年大同市土地利用/覆蓋遙感變化分析[J]. 干旱區(qū)資源與環(huán)境，2015，29(7)：68–75.Wang Dongchen, Du Peijun, Su Hongjun, et al. Land use/land cover changes in Datong based on remote sensing data[J]. Journal of Arid Land Resources and Environment,2015, 29(7): 68－75. (in Chinese with English abstract)

[11] 劉盛和. 城市土地利用擴(kuò)展的空間模式與動(dòng)力機(jī)制[J]. 地理科學(xué)進(jìn)展，2002，21(1)：43–50.Liu Shenghe. Spatial patterns and dynamic mechanisms of urban land use growth [J]. Progress in Geography, 2002,21(1): 43－50. (in Chinese with English abstract)

[12] 董雋，臧淑英. 大慶市城市土地利用變化的驅(qū)動(dòng)機(jī)制[J].地理研究，2011，30(1)：1121–1128.Dong Jun, Zang Shuying. Driving force mechanism of urban land use changes in Daqing city[J]. Geographical Research,2011, 30(1): 1121–1128. (in Chinese with English abstract)

[13] 吳桂平，曾永年，馮學(xué)智，等. CLUE-S模型的改進(jìn)與土地利用變化動(dòng)態(tài)模擬: 以張家界市永定區(qū)為例[J]. 地理研究，2010，29(3)：460–470.Wu Guiping, Zeng Yongnian, Feng Xuezhi, et al. Dynamic simulation of land use change based on the improved CLUE-S model: A case study of Yongding county,Zhangjiajie[J]. Geographical Research, 2010, 29(3): 460－470. (in Chinese with English abstract)

[14] 韓會(huì)然，楊成鳳，宋金. 北京市土地利用空間格局演化模擬及預(yù)測(cè)[J]. 地理科學(xué)進(jìn)展，2015，34(8)：976–986.Han Huiran, Yang Chengfeng, Song Jin. Simulation and projection of land-use change in Beijing under different scenarios [J]. Progress in Geography, 2015, 34(8): 976－986.(in Chinese with English abstract)

[15] 張楊，劉艷芳，顧漸萍，等. 武漢市土地利用覆被變化與生態(tài)環(huán)境效應(yīng)研究[J]. 地理科學(xué)，2011，31(10)：1280–1285.Zhang Yang, Liu Yanfang, Gu Jianping, et al. Land use /land cover change and its environmental effects in Wuhan city [J].Scientia Geographica Sinica, 2011, 31(10): 1280－1285. (in Chinese with English abstract)

[16] 玉蘇普江·艾麥提，阿里木江·卡斯木，艾克拜爾·買(mǎi)提尼牙孜. 基于遙感的烏魯木齊市景觀格局演變分析[J]. 西北林學(xué)院學(xué)報(bào)，2015，30(2)：172–179.Yusupujiang·Aimaiti, Alimujiang·Kasimu, Aikebaier·Maitiniyazi.Landscape dynamic change analysis of Urumqi based on remote sensing[J]. Journal of Northwest Forestry University,2015, 30(2): 172－179. (in Chinese with English abstract)

[17] Sexton J O, Urban D L, Donohue M J, et al. Long-term land cover dynamics by multi-temporal classification across the Landsat-5 record[J]. Remote Sensing of Environment, 2013,128(2): 246－258.

[18] Lu Dengsheng, Weng Qihao. Use of impervious surface in urban land-use classification [J]. Remote Sensing of Environment, 2006, 102(1): 146－160.

[19] 吳志杰，趙書(shū)河. 基于TM圖像的“增強(qiáng)的指數(shù)型建筑用地指數(shù)”研究[J]. 國(guó)土資源遙感，2012，93(2)：50–55.Wu Zhijie, Zhao Shuhe. A study of enhanced index- based built-up index based on Landsat TM imagery [J]. Remote Sensing For Land & Resources, 2012, 93(2): 50－55. (in Chinese with English abstract)

[20] Ghulam A, Ghulam O, Maimaitijiang M, et al. Remote Sensing based spatial statistics to document tropical rainforest transition pathways [J]. Remote Sensing, 2015, 7:6257－6279.

[21] 黃粵，陳曦，包安明，等. 近15a烏魯木齊市城市用地?cái)U(kuò)展動(dòng)態(tài)及其空間特征研究[J]. 冰川凍土，2006，28(3)：364－370.Huang Yue, Chen Xi, Bao Anming, et al. Urban sprawl pattern and sprawl features of Urumqi city during the last 15 years [J]. Journal of Glaciology and Geocryology, 2006,28(3): 364－370. (in Chinese with English abstract)

[22] 王珊珊，陳曦，段含明，等. 城市地表溫度對(duì)土地利用/覆被變化響應(yīng)的遙感研究：以烏魯木齊為例[J]. 中國(guó)沙漠，2012，32(3)：878－884.Wang Shanshan, Chen Xi, Duan Hanming, et al. Study on response of land surface temperature to land use and land cover change using remote sensing data: A case on Urumqi,China[J]. Journal of Desert Research, 2012, 32(3): 878－884.(in Chinese with English abstract)

[23] Huete A R. A soil-adjusted vegetation index (SAVI) [J].Remote Sensing of Environment, 1988, 25(3): 295－309.

[24] Zha Y, Gao J, Ni S. Use of normalized difference built-up index in automatically mapping urban areas from TM imagery[J]. Remote Sensing, 2003, 24(3): 583－594.

[25] Xu Hanqiu. Modification of normalized difference water index (NDWI) to enhance open water features in remotely sensed imagery[J]. International Journal of Remote Sensing,2006, 27(14): 3025－3033.

[26] Vol N. Accuracy assessment of satellite derived land-cover data[J]. Photogrammetric Engineering & Remote Sensing ,1994, 60 (6797): 691－691.

[27] 支小軍，丁偉. “絲綢之路經(jīng)濟(jì)帶”戰(zhàn)略背景下天山北坡城市群的功能定位與布局[J]. 新疆農(nóng)墾經(jīng)濟(jì)，2014(3)：24－27.

[28] 王雪微，王士君，宋飏，等. 交通要素驅(qū)動(dòng)下的長(zhǎng)春市土地利用時(shí)空變化[J]. 經(jīng)濟(jì)地理，2015，35(4)：156－161.Wang Xuewei, Wang Shijun, Song Yang, et al. Changchun land use spatio-temporal variation under the transportation elements’ driving[J]. Economic Geography, 2015, 35(4): 156－161. (in Chinese with English abstract)

[29] 趙國(guó)梁，鄭新奇，原智遠(yuǎn)，等. 道網(wǎng)密度與城鎮(zhèn)擴(kuò)張時(shí)空變化特征及關(guān)聯(lián)性分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2015，31(12)：220－229.Zhao Guoliang, Zhen Xinqi, Yuan Zhiyuan, et al. Spatial and temporal characteristics and correlation analysis of road networks and urban sprawl[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2015, 31(12): 220－229. (in Chinese with English abstract)

[30] 烏魯木齊統(tǒng)計(jì)局.烏魯木齊市統(tǒng)計(jì)年鑒2012[M]. 北京：中國(guó)統(tǒng)計(jì)出版社，2012.

[31] 史培軍，陳晉. 深圳市土地利用變化機(jī)制分析[J]. 地理學(xué)報(bào)，2000，55(5)：151－160.Shi Peijun, Chen Jin. Land use change mechanism in Shenzhen city[J]. Acta Geographica Sinica, 2000, 55(5): 151－160. (in Chinese with English abstract)

[32] 謝葉偉，劉兆剛，趙軍，等. 基于RS與GIS的典型黑土區(qū)土地利用變化分析: 以海倫市為例[J]. 地理科學(xué)，2010，30(3)：428－434.Xie Yewei, Liu Zhaogang, Zhao Jun, et al. Analysis of land use dynamic in typical black soil region based on GIS and RS:A case study of Hailun County[J]. Scientia Geographica Sinica, 2010, 30(3): 428－434. (in Chinese with English abstract)

[33] 喬偉峰，盛業(yè)華，方斌，等. 基于轉(zhuǎn)移矩陣的高度城市化區(qū)域土地利用演變信息挖掘: 以江蘇省蘇州市為例[J]. 地理研究，2013，32(8)：1497－1507.Qiao Weifeng, Sheng Yehua, Fang Bin, et al. Land use change information mining in highly urbanized area based on transfer matrix: A case study of Suzhou, Jiangsu Province[J].Geographical Research, 2013, 32(8): 1497－1507. (in Chinese with English abstract)