1987—2018年麗江壩區(qū)景觀格局變化軌跡分析*

田瀟然,王錦,余哲修

(1.西南林業(yè)大學(xué) 園林園藝學(xué)院，云南 昆明 650224；2.昆明理工大學(xué) 建筑與城市規(guī)劃學(xué)院，云南 昆明 650500；3.北京林業(yè)大學(xué) 林學(xué)院，北京 100083)

景觀變化是人類社會(huì)經(jīng)濟(jì)活動(dòng)與地域環(huán)境間相互干擾、相互作用的復(fù)雜動(dòng)態(tài)過(guò)程[1]。產(chǎn)業(yè)結(jié)構(gòu)的轉(zhuǎn)變、新興經(jīng)濟(jì)的進(jìn)入、地方政策等都能促使人類對(duì)土地的利用類型發(fā)生轉(zhuǎn)變，直接影響景觀格局的空間結(jié)構(gòu)[2-3]。為了能夠更為清晰地揭示景觀格局的動(dòng)態(tài)變化特征，將變化軌跡分析模型與土地利用/覆被、景觀格局相結(jié)合進(jìn)行分析，在揭示地域景觀的時(shí)空變化特征的研究中具有較強(qiáng)的優(yōu)勢(shì)。

針對(duì)景觀變化軌跡，學(xué)者們進(jìn)行了大量的深入探討，研究?jī)?nèi)容包括自然生態(tài)系統(tǒng)、城鎮(zhèn)、河流流域等種類眾多、不同尺度的景觀類型。汪東川等[4]通過(guò)對(duì)1992—2008年的天水藉河流域的土地利用/覆被變化軌跡分析，探討土地利用驅(qū)動(dòng)力與水土保持生態(tài)工程措施間的關(guān)聯(lián)性。劉桂林等[5]以塔里木河下游斷流區(qū)為研究對(duì)象，分析得到不同植被的變化軌跡。黃勇等[6]對(duì)環(huán)渤海地區(qū)的景觀格局變化軌跡進(jìn)行分析，結(jié)果顯示耕地是受人類干擾最為劇烈的地類，也是向人工表面轉(zhuǎn)換最大的地類。程琳琳等[7]構(gòu)建了1987—2017年寧夏鹽池縣的景觀格局演變軌跡，為長(zhǎng)期的生態(tài)恢復(fù)檢測(cè)提供支持。孔祥倫等[8]利用土地利用/覆被變化軌跡模型分析1986—2016年間黃河三角洲景觀格局的時(shí)空變化特征，指出人類活動(dòng)是土地利用變化過(guò)程中的主導(dǎo)因素。麗江壩區(qū)內(nèi)的大研古鎮(zhèn)、束河古鎮(zhèn)和白沙古鎮(zhèn)共同列入世界物質(zhì)文化遺產(chǎn)名錄,其豐富的民族旅游資源和多年地方政策的扶持使得當(dāng)?shù)氐穆糜伍_(kāi)發(fā)模式日趨完善，游客量逐年遞增[9-10]。隨著大量游客的涌入，壩區(qū)內(nèi)的服務(wù)業(yè)、房地產(chǎn)業(yè)、城鎮(zhèn)化進(jìn)程也在旅游業(yè)的推動(dòng)下得到了大力發(fā)展[11-12]。麗江壩區(qū)內(nèi)的生產(chǎn)方式從早年間以農(nóng)業(yè)為主逐步向旅游業(yè)過(guò)度的過(guò)程中帶動(dòng)了土地利用類型、景觀格局結(jié)構(gòu)的轉(zhuǎn)變，這使得麗江壩區(qū)成為了景觀格局變化研究的優(yōu)秀案例。

綜上所述，將變化軌跡模型應(yīng)用在土地利用的分析中，能夠更清晰地從時(shí)間和空間2個(gè)維度上得到景觀格局的變化規(guī)律和趨勢(shì)[13]。同時(shí)，關(guān)于麗江地區(qū)的研究多集中于麗江市城市核心區(qū)的旅游業(yè)發(fā)展[14]、古鎮(zhèn)形態(tài)變化等領(lǐng)域[15]，鮮有運(yùn)用定量分析方法對(duì)麗江景觀格局變化進(jìn)行探討的研究[16]。本研究利用變化軌跡模型分析麗江壩區(qū)32 a的土地利用及景觀格局的時(shí)空變化過(guò)程，麗江壩區(qū)在申遺成功及開(kāi)展旅游業(yè)前后過(guò)程中景觀格局的變化趨勢(shì)，全面探析不同景觀類型在研究時(shí)段內(nèi)的空間結(jié)構(gòu)特征，彌補(bǔ)麗江地區(qū)景觀格局研究的不足，為麗江地區(qū)地域景觀的保護(hù)、土地合理開(kāi)發(fā)利用提供有益的參考，并為政府決策提供科學(xué)依據(jù)。

1 研究區(qū)概況

麗江為云南省的地級(jí)市，地處滇西北高原，青藏高原東南部。本研究區(qū)為麗江壩區(qū)，壩區(qū)內(nèi)包括麗江市下轄的古城區(qū)以及玉龍納西縣的白沙鎮(zhèn)和黃山鎮(zhèn)的核心部分。研究區(qū)位于100°10′0″～100°20′0″E和26°45′0″～27°05′0″N之間，壩區(qū)總面積為194.39 km2，壩區(qū)四周環(huán)山，壩區(qū)內(nèi)地勢(shì)較為平坦，年平均溫度14.1 ℃，年降雨量1 205 mm，屬于高原山地季風(fēng)氣候，研究區(qū)地理位置見(jiàn)圖1。

圖1 研究區(qū)的地理位置Fig.1 The map of location of the study area

2 研究方法

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

研究影像來(lái)源于美國(guó)地質(zhì)勘探局(United States Geological Survey， USGS)，選擇1987—2018年共32期的Landsat影像，研究區(qū)成像時(shí)間均為植被生長(zhǎng)季，該區(qū)域未被云層覆蓋。分別對(duì)每期影像進(jìn)行大氣校正，影像裁剪等預(yù)處理。影像信息見(jiàn)表1。

表1 遙感影像信息

2.2 NDVI變化趨勢(shì)及時(shí)間節(jié)點(diǎn)確定

(1)NDVI變化趨勢(shì) 歸一化植被指數(shù)(normalized difference vegetation index，NDVI)除了被廣泛運(yùn)用于植物生長(zhǎng)狀態(tài)和植被覆蓋變化研究中，同時(shí)在NDVI值(-1～1)中，植被大于0.1、裸土在0～0.1之間、水體小于0，在一定程度上能夠反應(yīng)出研究區(qū)的土地利用類型的綜合情況[17-19]。

對(duì)麗江壩區(qū)32 a的NDVI值統(tǒng)計(jì)得到NDVI時(shí)序曲線，運(yùn)用局部加權(quán)回歸法(locally weighted regression，LOWESS)對(duì)曲線平滑處理[20]。由圖2可知，1987至1991年NDVI呈下降趨勢(shì)，1991至1994年為上升狀態(tài)，1994至2010年呈現(xiàn)一個(gè)總體下降的趨勢(shì)，而在2010到2018年則緩慢上升。NDVI年變化趨勢(shì)出現(xiàn)了1987、1991、1994、2010和2018五個(gè)主要的突變點(diǎn)。

圖2 1987—2018年NDVI指數(shù)平均值的變化趨勢(shì)Fig.2 The trend chart of average of NDVIduring 1987—2018

(2)研究區(qū)主要時(shí)間節(jié)點(diǎn)社會(huì)發(fā)展情況 1985年，國(guó)家將麗江列為對(duì)外開(kāi)發(fā)的乙類地區(qū)[21]。在1985年11月《麗江縣城總體規(guī)劃》中，1983—1990年為近期規(guī)劃階段，麗江縣城以發(fā)展農(nóng)林牧產(chǎn)品加工工業(yè)為主，積極為開(kāi)展旅游業(yè)創(chuàng)造條件。1991年為麗江納西族自治縣成立30周年，同年向云南省人民政府申報(bào)《關(guān)于遷建麗江民用飛機(jī)場(chǎng)的報(bào)告》，1994年麗江客運(yùn)公司投入運(yùn)營(yíng)，并批準(zhǔn)了實(shí)施大研古城“五四三二一”工程， 1997年12月，麗江壩區(qū)內(nèi)的大研古鎮(zhèn)、龍泉村(現(xiàn)束河古鎮(zhèn))和白沙鎮(zhèn)共同被列入世界文化遺產(chǎn)名錄[22]。2003年麗江撤地設(shè)市，壩區(qū)內(nèi)的基礎(chǔ)設(shè)施建設(shè)加快，城鎮(zhèn)的建設(shè)重心不斷向南移。2009年底麗江城市建設(shè)實(shí)現(xiàn)“三級(jí)跳”，旅游業(yè)受“昆明99世博會(huì)”的影響開(kāi)始迅猛發(fā)展。2010年麗江城市建成區(qū)綠地率達(dá)到36.78%，麗江市入選2010年中國(guó)10大最佳休閑城市[23]，同年麗江接待國(guó)內(nèi)外游客909.97×104人次[24]。2018年，旅游業(yè)得到了進(jìn)一步發(fā)展，年旅游業(yè)總收入高達(dá)998.45×108元[10]。

(3)研究時(shí)間節(jié)點(diǎn)確定 結(jié)合NDVI時(shí)序曲線與麗江壩區(qū)內(nèi)的社會(huì)與經(jīng)濟(jì)的發(fā)展、城市總體規(guī)劃，以及政策的發(fā)布等影響景觀格局變化的因素，研究時(shí)間節(jié)點(diǎn)確定為1987、1991、1994、2010和2018年。

2.3 土地利用類型分類

根據(jù)《土地利用現(xiàn)狀分類(2010—2017)》[25]以及麗江壩區(qū)的實(shí)際情況將研究區(qū)土地利用類型分為6類：1-耕地、2-林地、3-建設(shè)用地、4-草地、5-水體、6-未利用地。

對(duì)5期影像進(jìn)行監(jiān)督分類，分類結(jié)果結(jié)合同期的Google Earth影像進(jìn)行驗(yàn)證。5期影像解譯總體精度分別為：78.95%、82.88%、85.17%、89.26%、90.93%，Kappa系數(shù)分別為：0.734 7、0.778 2、0.810 6、0.862 9、0.887 4，解譯精度滿足研究需求，景觀分類結(jié)果見(jiàn)圖3。

圖3 各時(shí)間節(jié)點(diǎn)麗江壩區(qū)土地利用覆蓋情況Fig.3 The change of land use/cover in the Lijiangdam area of each time nodes

2.4 土地利用/覆被變化軌跡模型

對(duì)麗江壩區(qū)各時(shí)間節(jié)點(diǎn)的6種地類進(jìn)行變化軌跡分析，最多可以得到65=7776種變化軌跡，實(shí)際分析得到2 630種變化軌跡代碼。運(yùn)用柵格疊加計(jì)算將不同的柵格狀態(tài)標(biāo)記在圖譜中，分析得到土地利用類型、景觀類型等的空間變化特征[26]。

研究在對(duì)5期影像監(jiān)督分類的基礎(chǔ)上對(duì)分類結(jié)果進(jìn)行重分類和代數(shù)運(yùn)算，得到土地利用/覆被變化軌跡結(jié)果。當(dāng)土地利用類型小于10類時(shí)，變化軌跡計(jì)算公式[27]為：

Tij=(G1)ij×10n-1+(G2)ij×10n-2+…+(Gn)ij×10n-n

式中：Tij為軌跡分析計(jì)算結(jié)果中第i行、第j列的柵格軌跡代碼，表示土地利用類型或景觀類型的變化過(guò)程；n為時(shí)間節(jié)點(diǎn)的個(gè)數(shù)；(G1)ij、(G2)ij、…(Gn)ij為某一柵格在各個(gè)時(shí)間節(jié)點(diǎn)所對(duì)應(yīng)的土地利用類型代碼。

2.5 景觀格局指數(shù)計(jì)算

對(duì)5期影像的景觀格局進(jìn)行分析，分析32 a間麗江壩區(qū)景觀格局的時(shí)空變化規(guī)律特征。在參考前人的研究成果[28-30]的基礎(chǔ)上結(jié)合研究目的選擇斑塊類型所占景觀面積的比例(PLAND)、最大斑塊指數(shù)(LPI)、斑塊數(shù)量(NP)、聚集度指數(shù)(AI)、景觀形狀指數(shù)(LSI)、面積加權(quán)平均斑塊分維數(shù)(AWMPFD)，利用Fragstats軟件在類型水平上(class level)計(jì)算景觀格局的變化值。

3 結(jié)果與分析

3.1 麗江壩區(qū)土地利用/覆被變化特征分析

麗江壩區(qū)32 a間6種地類的動(dòng)態(tài)變化見(jiàn)表2。從表2可以看出建設(shè)用地為增幅最大且為持續(xù)增加的用地類型，其中1987—1994年緩慢增長(zhǎng)，1994—2018年則快速增長(zhǎng)，整體增長(zhǎng)達(dá)到356.45%，其面積共增加4 629.24 hm2。水體面積在整個(gè)研究期間的變化幅度非常微弱，研究期間增幅為49.45%，面積增加117.54 hm2。耕地、林地、草地、未利用地則都呈現(xiàn)總體下降的趨勢(shì)，其中耕地出現(xiàn)了先增長(zhǎng)后降低的特征，耕地面積共減少了2 175.30 hm2，下降幅度為29.87%。草地、未利用地、林地的面積總體下降幅度分別為40.31%、22.82%、9.26%。草地面積呈現(xiàn)先減少后增加的特征，1987—2010年持續(xù)下降，2010—2018年出現(xiàn)增長(zhǎng)。未利用地面積總體呈現(xiàn)出緩慢波動(dòng)的變化特征，1987—1991年增長(zhǎng)4.47%，1991—1994年下降10.05%，1994—2010年再增長(zhǎng)6.81%，2010—2018年較快速下降23.10%。

表2 1987—2018麗江壩區(qū)土地利用/覆被變化

3.1.1 土地利用/覆被變化軌跡過(guò)程分析

研究將每種軌跡代碼對(duì)應(yīng)的柵格數(shù)進(jìn)行排序，為了能夠清晰地針對(duì)主要變化軌跡進(jìn)行信息提取，篩選出柵格數(shù)超過(guò)800的軌跡代碼，再排除代碼為11111、22222、33333、44444、55555的5種未變化的代碼，最終獲得28種主要變化軌跡，見(jiàn)表3。

表3 主要軌跡代碼、面積及其所占比例

從表3中可以看出，28種主要變化軌跡中10種變化軌跡(11113、11133、11163、14133、41113、41133、43333、44133、66633、66663)都是由耕地、草地、未利用地轉(zhuǎn)向建設(shè)用地的，共占總變化面積的23.11%，其中由耕地轉(zhuǎn)向建設(shè)用地的面積占總變化面積的15.15%，耕地為壩區(qū)內(nèi)城鎮(zhèn)化建設(shè)的主要征用地類。除了耕地外，占總變化面積4.53%的草地和3.43%的未利用地轉(zhuǎn)向建設(shè)用地。耕地除了被征用為建設(shè)用地外，代碼為11114、11116、11166、16166的4種軌跡顯示由耕地轉(zhuǎn)向草地和未利用地，占總變化面積的7.40%。部分耕地和林地短暫地轉(zhuǎn)向未利用地、草地的現(xiàn)象，后又恢復(fù)為耕地或林地。值得注意的是有460.71 hm2的未利用地轉(zhuǎn)向了林地，占總變化面積的3.57%；750.69 hm2的未利用地轉(zhuǎn)向耕地，占總變化面積的5.82%。由此看出，在麗江壩區(qū)快速城鎮(zhèn)化過(guò)程中未利用地不僅用于城鎮(zhèn)的建設(shè)用地，也有少部分用于補(bǔ)充耕地的不足。

3.1.2 土地利用/覆被時(shí)空變化分析

土地利用/覆被主要變化軌跡時(shí)空分布見(jiàn)圖4。

圖4 1987—2018年麗江壩區(qū)土地利用/覆被 主要變化軌跡及其空間分布Fig.4 The spatial distribution of land use/cover main changetrajectory in Lijiang dam area during 1987—2018

由圖4可知，1994—2010年間核心城區(qū)北、西、南側(cè)的耕地首先被征用為城市的建設(shè)用地，2010—2018年間西側(cè)、北側(cè)的耕地繼續(xù)轉(zhuǎn)為建設(shè)用地，城區(qū)東側(cè)新增建設(shè)用地與耕地形成相互交錯(cuò)的形態(tài)。

在城市向四周擴(kuò)展的進(jìn)程中建設(shè)用地間隙中零星耕地轉(zhuǎn)為未利用地。玉龍新城位于城區(qū)南側(cè)，是由耕地和未利用地轉(zhuǎn)入的建設(shè)用地建設(shè)而成，與主城區(qū)之間的耕地相間隔。2010—2018年間，壩區(qū)北部村落周邊、中部城市周邊以及南部耕地中都有分散的耕地轉(zhuǎn)為未利用地。

通過(guò)變化軌跡分析可知，土地利用/覆被時(shí)空變化的主要特征表現(xiàn)為：建設(shè)用地首先呈弧形緊密包圍城市核心區(qū)，后向東側(cè)發(fā)展與耕地形成交錯(cuò)的松散布局形態(tài)，核心區(qū)周邊的耕地變化顯著，原分散的村落、城鎮(zhèn)逐漸被新城市包圍或半包圍。由于麗江壩區(qū)為南北向長(zhǎng)、東西向窄的山間盆地，地形的限制使城市核心區(qū)由放射狀逐漸向近長(zhǎng)方形的形態(tài)發(fā)展；壩區(qū)中部土地利用/覆被變化劇烈，北部的未利用地以及南部的耕地變化并不顯著。

3.2 麗江壩區(qū)景觀格局變化軌跡

麗江壩區(qū)景觀格局變化軌跡見(jiàn)表4。從斑塊類型所占景觀面積的比例(PLAND)計(jì)算結(jié)果可以看出，耕地轉(zhuǎn)向建設(shè)用地、未利用地轉(zhuǎn)向耕地、未利用地轉(zhuǎn)向建設(shè)用地(11113、11133、66661、66663)3類變化軌跡在整個(gè)景觀中占比較高，PLAND值分別為：2.47、1.58、1.57、0.72，這說(shuō)明研究區(qū)內(nèi)發(fā)生深刻變化的主要斑塊類型都與人類活動(dòng)密切相關(guān)。在人類的干擾下壩區(qū)內(nèi)景觀格局的破碎度和景觀異質(zhì)性明顯加劇，斑塊數(shù)量(NP)大于500的變化軌跡有11113、11114、41113、66661。最大斑塊指數(shù)(LPI)值大于0.1的有11113、11133、66661，說(shuō)明主要變化軌跡中的斑塊面積都較小，其中11113的LPI值最大，為0.22，由耕地轉(zhuǎn)向建設(shè)用地以及未利用地轉(zhuǎn)向耕地的斑塊為整個(gè)壩區(qū)內(nèi)的優(yōu)勢(shì)斑塊類型。LPI為0.01的有14133、41113、41133、44133、66664、66665，主要是零星耕地轉(zhuǎn)為草地再轉(zhuǎn)為建設(shè)用地、草地轉(zhuǎn)向耕地再轉(zhuǎn)為建設(shè)用地、未利用地轉(zhuǎn)向草地或水體的變化類型。

表4 麗江壩區(qū)的景觀格局變化軌跡

從景觀形狀指數(shù)(LSI)計(jì)算結(jié)果中可知，由未利用地轉(zhuǎn)向水體的變化類型(66665)LSI值最小(10.88)，其形狀較為整齊。而受人類活動(dòng)干擾大的變化軌跡11113的最大LSI值為32.97，其形狀最為復(fù)雜，其次為11114和66661，LSI值為30.36、30.04。面積加權(quán)平均斑塊分維數(shù)(AWMPFD)中11133最大，為1.16，空間形狀最為復(fù)雜，其次為11113、22622，AWMPFD值分別為1.15、1.13。從LSI和AWMPFD的計(jì)算結(jié)果中可以看出人類活動(dòng)與斑塊形狀的復(fù)雜程度密切相關(guān)。

從聚集度(AI)計(jì)算結(jié)果可知，11133、11113的AI值為74.88、69.72，聚集度較高，這些變化軌跡集中于城市核心區(qū)周圍成片的建設(shè)用地，66622的AI值為66.70，集中在玉龍新城東南側(cè)。較高的聚集度以及較大的面積，說(shuō)明變化軌跡是人類大規(guī)模開(kāi)發(fā)活動(dòng)造成的。草地轉(zhuǎn)向建設(shè)用地的斑塊類型(41133、44133、14133)為主要變化軌跡中最分散的變化軌跡。

綜上所述，促使景觀格局快速變化的最主要驅(qū)動(dòng)軌跡集中在受到旅游城鎮(zhèn)化深刻影響的城市核心區(qū)周圍，斑塊形狀的復(fù)雜性以及景觀的破碎性都是受到人類開(kāi)發(fā)利用影響而形成的。麗江壩區(qū)景觀32 a的變化中人類活動(dòng)成為了驅(qū)動(dòng)景觀變化的最主要因素。

4 討論與結(jié)論

(1)在1987—2018年的變化過(guò)程中，麗江壩區(qū)內(nèi)的建設(shè)用地面積增幅最大(356.45%)，面積為4 629.24 hm2，其主要來(lái)源為征用城市核心區(qū)及村落周邊的耕地、草地和未利用地。耕地呈現(xiàn)出先升后降的變化特征，1987—1994年間壩區(qū)的生產(chǎn)方式以農(nóng)耕為主， 1994—2018年期間壩區(qū)旅游業(yè)的興起以及麗江市政府辦公區(qū)南遷造成耕地大面積流失，共減少2 175.30 hm2，降幅為29.87%，為壩區(qū)內(nèi)面積減少最多的景觀類型。

(2)壩區(qū)內(nèi)土地利用/覆被變化過(guò)程中受到人類活動(dòng)干擾的影響深刻，以城市核心區(qū)周邊的變化最為劇烈，壩區(qū)北部和南部變化并不顯著。通過(guò)對(duì)麗江壩區(qū)土地利用/覆被軌跡變化分析，有占總面積66.31%的土地發(fā)生了變化，共有12 889.53 hm2，城區(qū)人口增長(zhǎng)、政策扶持旅游業(yè)、旅游經(jīng)濟(jì)發(fā)展、游客涌入、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、旅游城鎮(zhèn)化等人類活動(dòng)都成為土地利用/覆被快速變化的主要驅(qū)動(dòng)因素。

(3)通過(guò)對(duì)壩區(qū)的景觀格局變化軌跡分析，土地轉(zhuǎn)變?yōu)榻ㄔO(shè)用地的變化軌跡中景觀斑塊所呈現(xiàn)出的較高的破碎度、不規(guī)則形狀，這都與人類活動(dòng)干擾高度相關(guān)。城市核心區(qū)周邊由耕地轉(zhuǎn)向建設(shè)用地的景觀斑塊最為集中、連接度高，成為景觀格局變化中的優(yōu)勢(shì)斑塊類型；由草地轉(zhuǎn)向建設(shè)用地的變化軌跡斑塊遠(yuǎn)離城市核心區(qū)、聚集度最低，分散于人類活動(dòng)干擾較低的耕地間。在人類的社會(huì)、經(jīng)濟(jì)活動(dòng)影響下城市核心區(qū)以及其周邊的破碎度、聚集度、景觀格局空間結(jié)構(gòu)的復(fù)雜程度均高于遠(yuǎn)離核心區(qū)的北部和南部地區(qū)。

由旅游服務(wù)業(yè)推動(dòng)區(qū)域城鎮(zhèn)化以及地方產(chǎn)業(yè)結(jié)構(gòu)調(diào)整從而導(dǎo)致景觀格局發(fā)生深刻變化[13]，麗江壩區(qū)乃至麗江市區(qū)成為了一個(gè)由社會(huì)經(jīng)濟(jì)發(fā)展引起景觀格局變化的景觀研究關(guān)注地[16，21-23]。本研究通過(guò)運(yùn)用軌跡分析法與土地利用/覆被、景觀格局相結(jié)合的方法從時(shí)間、空間維度探索麗江壩區(qū)32 a間土地利用的變化過(guò)程，較為全面地揭示出景觀變化的特征。經(jīng)過(guò)對(duì)主要變化軌跡的分析，能夠清晰地梳理出各土地類型的變化方向以及景觀格局的變化過(guò)程，在篩選引起景觀變化的主要驅(qū)動(dòng)因子方面有較為明顯的優(yōu)越性。

景觀格局變化一直以來(lái)都是景觀研究的熱點(diǎn)議題，而揭示其發(fā)展變化的方法眾多[31-32]。本研究對(duì)麗江壩區(qū)的景觀時(shí)空變化過(guò)程和特征進(jìn)行研究，初步探討了影響景觀變化的主要驅(qū)動(dòng)因素，但在景觀破碎化、空間異質(zhì)性趨勢(shì)以及快速旅游城鎮(zhèn)化對(duì)景觀生態(tài)系統(tǒng)的穩(wěn)定性風(fēng)險(xiǎn)評(píng)價(jià)等方面的研究還需進(jìn)一步深入探索。