許冬凡 張加龍 鮑 瑞 廖 易 馮亞飛

（1. 西南林業(yè)大學(xué)林學(xué)院，云南 昆明 650233；2. 昆明市信息中心，云南 昆明 650500）

土地利用/覆蓋變化（LUCC）越來(lái)越被認(rèn)為是全球環(huán)境變化與可持續(xù)發(fā)展研究的重要內(nèi)容[1]。土地利用/覆蓋變化對(duì)區(qū)域生態(tài)環(huán)境產(chǎn)生重要影響，并以累積的方式影響全球環(huán)境變化[2]。Landsat衛(wèi)星系列圖像則是研究土地覆蓋類(lèi)型和變化的重要數(shù)據(jù)源[3]。傳統(tǒng)的使用遙感分析軟件進(jìn)行土地覆蓋分類(lèi)的工程量過(guò)于巨大，使用Google Earth Engine（GEE）云平臺(tái)是一個(gè)更好的選擇。近年來(lái)，GEE 平臺(tái)被廣泛使用于基于時(shí)間序列數(shù)據(jù)的土地利用/覆蓋變化研究中，而長(zhǎng)時(shí)間序列的區(qū)域研究對(duì)于研究該區(qū)域的發(fā)展特征、發(fā)展模式具有非常重要的意義[4]。縱觀已有的基于GEE 平臺(tái)的進(jìn)行土地覆蓋分類(lèi)的研究，研究對(duì)象主要集中在大型城市群、經(jīng)濟(jì)發(fā)達(dá)地區(qū)[5-7]，對(duì)于西南地區(qū)的研究較少。而基于長(zhǎng)時(shí)間序列遙感影像數(shù)據(jù)對(duì)于滇池流域進(jìn)行土地覆蓋分析則未見(jiàn)報(bào)道。有關(guān)滇池流域的研究多集中于水環(huán)境、景觀生態(tài)、流域治理等方面[8-10]，對(duì)于滇池流域制作長(zhǎng)時(shí)間序列的土地覆蓋數(shù)據(jù)是必要的。

建設(shè)用地的擴(kuò)展能直觀體現(xiàn)一個(gè)區(qū)域的城市化水平[11]，研究城鎮(zhèn)用地?cái)U(kuò)展的驅(qū)動(dòng)機(jī)制不僅有助于深入理解城市擴(kuò)展的規(guī)律及方向，而且對(duì)于推動(dòng)區(qū)域協(xié)調(diào)發(fā)展具有重要意義[12]。驅(qū)動(dòng)機(jī)制著重分析區(qū)域土地利用變化過(guò)程的主要自然和社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)因素及其驅(qū)動(dòng)機(jī)理[13]。在自然系統(tǒng)中，氣候、土壤、水文等被認(rèn)為是主要的驅(qū)動(dòng)力；在社會(huì)系統(tǒng)中，驅(qū)動(dòng)力通常有人口變化、貧富狀況、技術(shù)進(jìn)步、經(jīng)濟(jì)增長(zhǎng)、政治經(jīng)濟(jì)結(jié)構(gòu)等[14]。目前針對(duì)建設(shè)用地變化驅(qū)動(dòng)力的研究多以定性分析為主，定量分析較少。在方法上多采用相關(guān)性分析[15]、主成分分析[16]、Logistic 回歸分析[17]等方法。在驅(qū)動(dòng)力研究尺度上多以整個(gè)研究區(qū)為對(duì)象進(jìn)行驅(qū)動(dòng)力分析，但研究區(qū)內(nèi)部各行政區(qū)的自然環(huán)境特點(diǎn)、區(qū)域發(fā)展方針不同，以行政區(qū)劃為研究單元來(lái)分析土地覆蓋變化特征更為科學(xué)。因此，本研究基于GEE 云平臺(tái)，以滇池流域?yàn)檠芯繀^(qū)，基于長(zhǎng)時(shí)間序列的Landsat 數(shù)據(jù)進(jìn)行土地覆蓋變化分析，以行政區(qū)劃為單元來(lái)研究土地覆蓋變化的特點(diǎn)，并結(jié)合多項(xiàng)驅(qū)動(dòng)因子來(lái)研究流域內(nèi)不同時(shí)段的建設(shè)用地?cái)U(kuò)張的驅(qū)動(dòng)力。以期揭示滇池流域的發(fā)展特征，為后續(xù)學(xué)者研究滇池流域的發(fā)展提供科學(xué)數(shù)據(jù)支撐，為實(shí)現(xiàn)流域的可持續(xù)發(fā)展和自然資源的永續(xù)利用提供科學(xué)參考。

1 材料與方法

1.1 研究區(qū)概況

滇池流域?qū)倏偯娣e為2 920 km2，包含了昆明市轄區(qū)范圍，內(nèi)有五華區(qū)、官渡區(qū)、盤(pán)龍區(qū)、西山區(qū)、呈貢區(qū)、晉寧縣和嵩明縣。滇池流域占云南省總國(guó)土面積的0.7%，聚集著云南省7.6%的人口，創(chuàng)造高達(dá)云南省24%的GDP[18]。隨著昆明市的城市化進(jìn)程加快，與解放初期相比，滇池水面面積已減少至少13 km2，林地面積減少至少200 km2[9]。

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

采用的遙感影像為L(zhǎng)andsat 系列衛(wèi)星影像，涵蓋了從1988—2018 年的全部滇池流域影像，共采用63 景影像，所有Landsat 數(shù)據(jù)處理均基于GEE云平臺(tái)，通過(guò)代碼編寫(xiě)進(jìn)行處理分析，詳細(xì)介紹請(qǐng)參考官方網(wǎng)站（https://developers.google.com/earthengine/datasets/）。研究所采用的氣象數(shù)據(jù)均來(lái)源于中國(guó)氣象數(shù)據(jù)網(wǎng)（http://data.cma.cn），主要包含年平均氣溫、年平均降水、年平均降水、年日照時(shí)長(zhǎng)、年極端氣溫、年極端降水、年均氣壓等數(shù)據(jù)。

研究所采用的社會(huì)經(jīng)濟(jì)數(shù)據(jù)均來(lái)源于云南省統(tǒng)計(jì)局所公布的《云南省統(tǒng)計(jì)年鑒》[19]（1989—2019 年），主要包括地區(qū)生產(chǎn)總值、戶籍人口、城鎮(zhèn)人口比重、工農(nóng)業(yè)產(chǎn)值、各產(chǎn)業(yè)從業(yè)人數(shù)等。

1.3 基于GEE 平臺(tái)的監(jiān)督分類(lèi)

監(jiān)督分類(lèi)就是根據(jù)影像所展示的不同地物，繪制感興趣區(qū)（ROI）作為訓(xùn)練樣本，再通過(guò)選取分類(lèi)器進(jìn)行分類(lèi)。本研究基于GEE 平臺(tái)，采用Python 與JavaScript 語(yǔ)言，使用支持向量機(jī)的監(jiān)督分類(lèi)方法[7]。參考中科院土地覆蓋分類(lèi)標(biāo)準(zhǔn)，將滇池流域內(nèi)所有土地分成林地、草地、農(nóng)地、水體、建設(shè)用地、其他未利用這6 類(lèi)。在得到多期土地覆蓋的分類(lèi)結(jié)果后，使用檢驗(yàn)樣本，采取基于混淆矩陣的精度評(píng)價(jià)方法對(duì)于分類(lèi)結(jié)果進(jìn)行精度評(píng)價(jià)。

分類(lèi)的樣本需要保證均勻分布于研究區(qū)各處，各樣本類(lèi)之間的樣本分離度要控制在1.5 以上。本次實(shí)驗(yàn)的時(shí)間尺度為1988—2018 年共31年，每一年的訓(xùn)練樣本每類(lèi)選取30 個(gè)，一共180個(gè)訓(xùn)練樣本；檢驗(yàn)樣本每類(lèi)也選取30 個(gè)，一共180 個(gè)檢驗(yàn)樣本。

1.4 土地覆蓋變化測(cè)度分析

土地覆蓋變化動(dòng)態(tài)度[20]可以反映不同地類(lèi)的綜合動(dòng)態(tài)變化情況。

1.5 驅(qū)動(dòng)因子選取

在多種要素共同構(gòu)成形成驅(qū)動(dòng)力體系的同時(shí)，各個(gè)驅(qū)動(dòng)力因子之間存在著多種關(guān)聯(lián)模式[22]，因此要研究滇池流域建設(shè)用地的變化驅(qū)動(dòng)力，應(yīng)將多種要素關(guān)聯(lián)納入驅(qū)動(dòng)體系。選取的指標(biāo)如表1 所示，所有數(shù)據(jù)均來(lái)源于《云南省統(tǒng)計(jì)年鑒》[19]。

表1 驅(qū)動(dòng)力指標(biāo)體系Table 1 The driving force index system

1.6 驅(qū)動(dòng)力分析方法

本次研究選取社會(huì)經(jīng)濟(jì)數(shù)據(jù)來(lái)分析滇池流域土地覆蓋變化的驅(qū)動(dòng)力，使用熵值法以消除確定權(quán)重過(guò)程中的人為主觀因素，得出各個(gè)因子的驅(qū)動(dòng)力貢獻(xiàn)權(quán)重，再通過(guò)灰色關(guān)聯(lián)法[23]計(jì)算經(jīng)熵值法篩選后的驅(qū)動(dòng)因子與滇池流域建設(shè)用地變化之間的關(guān)聯(lián)度。

2 結(jié)果與分析

2.1 滇池流域土地覆蓋變化

2.1.1 分類(lèi)結(jié)果精度評(píng)價(jià)

使用基于混淆矩陣的方法進(jìn)行精度評(píng)價(jià)，得到各年份分類(lèi)結(jié)果的總體精度與Kappa 系數(shù)，各年份的分類(lèi)精度評(píng)價(jià)結(jié)果見(jiàn)表2。1992 年、2003年與2011 年由于影像質(zhì)量較差，故分類(lèi)精度稍低，但多年來(lái)的分類(lèi)結(jié)果均保持了良好的精度。

表2 精度評(píng)價(jià)表Table 2 Accuracy evaluation

2.1.2 土地覆蓋變化特征分析

滇池流域四面環(huán)山，林地主要分布于山地區(qū)域，環(huán)繞流域四周，是滇池流域內(nèi)面積占比最大的土地覆蓋類(lèi)型；草地主要集中于流域的東部、北部，雜亂分布于林地之間，為滇池流域占比第二多的土地覆蓋類(lèi)型；滇池坐落于流域中部，除滇池外仍有較多小微水體與大型水庫(kù)分布于流域內(nèi)；建設(shè)用地主要分布在昆明主城區(qū)、呈貢區(qū)、南部的晉寧縣和北部的嵩明縣，其面積增長(zhǎng)較為迅速；農(nóng)地多數(shù)環(huán)滇池而分布，多年來(lái)面積不斷減少，現(xiàn)存農(nóng)地主要集中在呈貢區(qū)、晉寧縣；流域內(nèi)的其他未利用地主要表現(xiàn)為裸地，多分布于建設(shè)用地周邊。

1988—2018 年的滇池流域土地覆蓋類(lèi)型總體構(gòu)成如圖1 所示，長(zhǎng)期以來(lái)面積持續(xù)減少的為草地、農(nóng)地、水體這3 類(lèi)；而持續(xù)增長(zhǎng)的為林地、建設(shè)用地、其他未利用地這3 類(lèi)。

圖1 土地覆蓋類(lèi)型結(jié)構(gòu)圖Fig. 1 Plot of land cover type structure

統(tǒng)計(jì)滇池流域各個(gè)時(shí)期的土地覆蓋數(shù)據(jù)，基于土地覆蓋變化動(dòng)態(tài)度、土地覆蓋年變化率得出1988—2018 年內(nèi)滇池流域各用地類(lèi)型的變化情況，見(jiàn)表3～4。

表3 各用地類(lèi)型面積變化Table 3 Changes of each land typekm2

表4 各用地類(lèi)型動(dòng)態(tài)變化Table 4 Dynamic changes of various land types%

2.1.3 土地覆蓋變化流向分析

多年來(lái)，滇池流域的6 種土地覆蓋類(lèi)型的類(lèi)型變化過(guò)程較為劇烈?；?988 年、1998 年、2008 年、2018 年的土地覆蓋分類(lèi)結(jié)果制作轉(zhuǎn)移矩陣弦圖，見(jiàn)圖2。1988—1998 年，林地變化程度較大，由林地轉(zhuǎn)化為其他未利用地轉(zhuǎn)的面積較大，同時(shí)草地、農(nóng)地、其他未利用地也有較多面積轉(zhuǎn)為林地；1998—2008 年，由林地、草地、水體、其他未利用地轉(zhuǎn)為建設(shè)用地的面積較多；2008—2018 年，由林地、草地、水體轉(zhuǎn)化為建設(shè)用地的面積較多；縱觀1988—2018 年的土地覆蓋流轉(zhuǎn)，建設(shè)用地為最大的轉(zhuǎn)入類(lèi)型，林地、草地的轉(zhuǎn)出面積較多。

圖2 土地覆蓋轉(zhuǎn)移矩陣弦圖Fig. 2 Chord plot of land cover transition matrix

2.1.4 土地覆蓋變化區(qū)劃分析

根據(jù)1988 年、1998 年、2008 年、2018 年的土地覆蓋結(jié)果，按照流域內(nèi)的行政區(qū)劃進(jìn)行分區(qū)統(tǒng)計(jì)，如圖3 所示。由于水體與其他未利用地的面積變化幅度較小，所以未對(duì)這2 種地類(lèi)進(jìn)行區(qū)劃分析。

圖3 各地類(lèi)行政區(qū)劃統(tǒng)計(jì)圖Fig. 3 Statistical maps of land cover types by various regions

由圖3 可知，各區(qū)縣的草地面積均在1988 年為最多，逐年減少，尤其是晉寧與嵩明在1988—2008 年間，草地面積減少了一半以上，而五華區(qū)、西山區(qū)、盤(pán)龍區(qū)的草地面積減少幅度較小。

盤(pán)龍區(qū)、五華區(qū)的林地面積多年來(lái)的增長(zhǎng)幅度較小，其余區(qū)縣的林地面積均持續(xù)增加，1988—2008 年是林地增長(zhǎng)面積最快的階段，晉寧縣域嵩明縣的林地面積在這一時(shí)段內(nèi)得以翻番。農(nóng)地整體趨勢(shì)呈現(xiàn)先增后減，五華區(qū)與盤(pán)龍區(qū)的農(nóng)地面積多年來(lái)始終維持在較低水平，其余區(qū)縣在1998—2008 年間的農(nóng)地面積達(dá)到最大，而后劇烈減少。建設(shè)用地為持續(xù)增長(zhǎng)較為顯著的地類(lèi)，五華區(qū)、官渡區(qū)、西山區(qū)、盤(pán)龍區(qū)作為昆明市的老城區(qū)，多年來(lái)的建設(shè)用地平穩(wěn)增長(zhǎng)，而呈貢區(qū)、晉寧縣、嵩明縣近年來(lái)城區(qū)擴(kuò)張速度持續(xù)增快。

根據(jù)土地利用程度模型[21]，對(duì)不同的土地覆蓋類(lèi)型進(jìn)行對(duì)應(yīng)的賦值：其他未利用地為1，林地、草地、水體為2，耕地為3，建設(shè)用地為4，進(jìn)而對(duì)滇池流域內(nèi)的7 個(gè)行政區(qū)分別計(jì)算土地利用程度綜合指數(shù)，結(jié)果見(jiàn)表5。

表5 各行政區(qū)土地利用程度綜合指數(shù)表Table 5 Comprehensive index of land utilization degree in each region

1988—1998 年屬于主城的五華區(qū)、西山區(qū)、官渡區(qū)、盤(pán)龍區(qū)這四區(qū)由于建設(shè)用地面積較多，具有較高的土地利用程度綜合指數(shù)，而晉寧縣由于農(nóng)地較多，此項(xiàng)指數(shù)也偏大；1998—2008 年，呈貢區(qū)與晉寧縣由于建設(shè)用地的增加，導(dǎo)致土地利用程度綜合指數(shù)增長(zhǎng)較快；2008—2018 年，僅五華區(qū)、官渡區(qū)與呈貢區(qū)的土地利用程度綜合指數(shù)增加明顯，其余各區(qū)縣未見(jiàn)顯著變化。1988—2018 年內(nèi)，官渡區(qū)、五華區(qū)與呈貢區(qū)的土地覆蓋變化較為劇烈，而嵩明、晉寧的變化稍微緩和。

1988—1998 年各區(qū)縣土地覆蓋變化進(jìn)程較為平緩，而1998—2008 年，各區(qū)縣均發(fā)生較大的變化，是滇池流域發(fā)展最為迅猛的一個(gè)階段，2008—2018 年的土地覆蓋變化進(jìn)程稍有減緩。

2.2 滇池流域建設(shè)用地變化驅(qū)動(dòng)力分析

滇池流域的城市規(guī)模較為龐大，近年來(lái)人口不斷增長(zhǎng)、經(jīng)濟(jì)水平逐年提升，定量研究其城市化進(jìn)程的驅(qū)動(dòng)因素具有重要意義。所以本研究選取滇池流域的建設(shè)用地為研究對(duì)象，進(jìn)行驅(qū)動(dòng)分析。

2.2.1 基于熵值法的驅(qū)動(dòng)因子篩選

本研究基于長(zhǎng)時(shí)間序列的土地覆蓋數(shù)據(jù)和對(duì)應(yīng)年份的驅(qū)動(dòng)因子，采用熵值法篩選出1988—2018 年對(duì)土地覆蓋變化具有較強(qiáng)關(guān)聯(lián)性的因子，篩選結(jié)果見(jiàn)表6，前5 項(xiàng)指標(biāo)的權(quán)重系數(shù)顯著高于其他指標(biāo)，后4 項(xiàng)指標(biāo)的權(quán)重系數(shù)則過(guò)低。信息熵值越大代表此項(xiàng)數(shù)據(jù)效用不足，去掉權(quán)重系數(shù)小于5%的指標(biāo)，將剩余11 項(xiàng)指標(biāo)代入灰色關(guān)聯(lián)法模型進(jìn)行關(guān)聯(lián)度計(jì)算。

表6 熵值法計(jì)算權(quán)重結(jié)果Table 6 The weight results calculated by entropy method

2.2.2 基于灰色關(guān)聯(lián)法的驅(qū)動(dòng)因子分類(lèi)

篩選后的11 項(xiàng)指標(biāo)即為1988—2018 年與滇池流域的土地覆蓋變化具有較強(qiáng)關(guān)聯(lián)性的驅(qū)動(dòng)因子。將滇池流域劃分為3 個(gè)時(shí)間段：1988—1998年、1998—2008 年、2008—2018 年，將對(duì)應(yīng)年份的土地覆蓋數(shù)據(jù)及篩選后的驅(qū)動(dòng)因子代入模型計(jì)算灰色關(guān)聯(lián)度，結(jié)果見(jiàn)表7。

表7 灰色關(guān)聯(lián)度結(jié)果Table 7 Result of grey association

將關(guān)聯(lián)度指標(biāo)按數(shù)值大小區(qū)分為3 個(gè)關(guān)聯(lián)等級(jí)：0～0.5 為低關(guān)聯(lián)，0.5～0.8 為中關(guān)聯(lián)，0.8～1.0 為高關(guān)聯(lián)。按關(guān)聯(lián)等級(jí)對(duì)灰色關(guān)聯(lián)度的結(jié)果進(jìn)行區(qū)分，1988—1998 年對(duì)滇池流域建設(shè)用地的增長(zhǎng)起高關(guān)聯(lián)的因子為：人均道路鋪裝面積和社會(huì)銷(xiāo)售品銷(xiāo)售總額；1998—2008 年對(duì)滇池流域建設(shè)用地的增長(zhǎng)起高關(guān)聯(lián)的因子為：第三產(chǎn)業(yè)從業(yè)人數(shù)比、人均道路鋪裝面積和社會(huì)銷(xiāo)售品銷(xiāo)售總額；2008—2018 年對(duì)滇池流域建設(shè)用地的增長(zhǎng)起高關(guān)聯(lián)的因子為：人均道路鋪裝面積、高等學(xué)校數(shù)量和客運(yùn)總量。1988—2018 年，人均道路鋪裝面積因子持續(xù)表現(xiàn)為高關(guān)聯(lián)，社會(huì)銷(xiāo)售品銷(xiāo)售總額因子與客運(yùn)總量因子有2 個(gè)階段表現(xiàn)為高關(guān)聯(lián)，而衛(wèi)生醫(yī)療機(jī)構(gòu)數(shù)、總?cè)丝跀?shù)、自然科學(xué)研究機(jī)構(gòu)數(shù)、國(guó)內(nèi)生產(chǎn)總值這4 項(xiàng)因子的多年來(lái)的關(guān)聯(lián)度較低。

按表1 的驅(qū)動(dòng)因子分類(lèi)來(lái)計(jì)算各年度驅(qū)動(dòng)因子的關(guān)聯(lián)度總值，1988—1998 年的驅(qū)動(dòng)力排序應(yīng)為：經(jīng)濟(jì)驅(qū)動(dòng)（2.082）＞人口驅(qū)動(dòng)（1.748）＞基礎(chǔ)建設(shè)驅(qū)動(dòng)（1.668）＞教育、衛(wèi)生條件驅(qū)動(dòng)（1.119）；1998—2008 年內(nèi)的驅(qū)動(dòng)力排序應(yīng)為：人口驅(qū)動(dòng)（2.347）＞經(jīng)濟(jì)驅(qū)動(dòng)（1.944）＞基礎(chǔ)建設(shè)驅(qū)動(dòng)（1.729）＞教育、衛(wèi)生條件驅(qū)動(dòng)（1.641）；2008—2018 年內(nèi)的驅(qū)動(dòng)力排序應(yīng)為：教育、衛(wèi)生條件驅(qū)動(dòng)（2.034）＞基礎(chǔ)建設(shè)驅(qū)動(dòng)（1.789）＞經(jīng)濟(jì)驅(qū)動(dòng)（1.514）＞人口驅(qū)動(dòng)（1.452）。1988—2008 年，經(jīng)濟(jì)驅(qū)動(dòng)因素與人口驅(qū)動(dòng)因素的貢獻(xiàn)比重較高，而2008—2018 年，教育、衛(wèi)生條件驅(qū)動(dòng)因素與基礎(chǔ)建設(shè)驅(qū)動(dòng)因素的貢獻(xiàn)比重較高。

3 結(jié)論與討論

本次研究基于GEE 云平臺(tái)，采用1988—2018年的Landsat 遙感影像對(duì)滇池流域進(jìn)行了土地覆蓋分類(lèi)，結(jié)合動(dòng)態(tài)度、年變化率、土地利用程度綜合指數(shù)等指標(biāo)，按行政區(qū)劃的不同來(lái)分析滇池流域多年來(lái)的土地覆蓋變化特征，并且結(jié)合多項(xiàng)驅(qū)動(dòng)因子，采用熵值法與灰色關(guān)聯(lián)法來(lái)分析不同時(shí)期內(nèi)的建設(shè)用地變化驅(qū)動(dòng)力。結(jié)果表明：

1）1988—2018 年，面積持續(xù)減少的為草地、農(nóng)地、水體這3 類(lèi)；而持續(xù)增長(zhǎng)的為林地、建設(shè)用地、其他未利用地這3 類(lèi)；滇池流域的主要土地覆蓋類(lèi)型為林地、草地與水體，三者面積總占比長(zhǎng)期超過(guò)60%；1988—2018 年，林地面積增長(zhǎng)最多，建設(shè)用地次之，面積減少最多的地類(lèi)為草地。

2）1988—2018 年，滇池流域的土地覆蓋變化過(guò)程較為劇烈，建設(shè)用地為最大的轉(zhuǎn)入類(lèi)型，林地、草地的轉(zhuǎn)出面積較多；流域內(nèi)土地覆蓋的主要流轉(zhuǎn)方式為草地轉(zhuǎn)為林地、建設(shè)用地，農(nóng)地轉(zhuǎn)化為建設(shè)用地、其他未利用地，其他未利用地轉(zhuǎn)化為草地、林地。

3）滇池流域在1988—1998 年的各區(qū)縣土地覆蓋變化進(jìn)程較為平緩，而1998—2008 年，各區(qū)縣均發(fā)生較大的變化，是滇池流域發(fā)展最為迅猛的階段，2008—2018 年內(nèi)的土地覆蓋變化進(jìn)程稍有減緩；1988—2018 年，官渡區(qū)、五華區(qū)與呈貢區(qū)的土地覆蓋變化進(jìn)程較為劇烈，而嵩明、晉寧的變化進(jìn)程相對(duì)緩和；晉寧、嵩明2 個(gè)縣的林地、草地總占比較大，呈貢與晉寧的農(nóng)地總占比較大。

4）各因子在不同時(shí)期內(nèi)對(duì)于滇池流域建設(shè)用地變化的驅(qū)動(dòng)作用也不同；1988—2018 年，人均道路鋪裝面積因子持續(xù)表現(xiàn)為高關(guān)聯(lián)，社會(huì)銷(xiāo)售品銷(xiāo)售總額因子與客運(yùn)總量因子有兩個(gè)階段表現(xiàn)為高關(guān)聯(lián)；而衛(wèi)生醫(yī)療機(jī)構(gòu)數(shù)、總?cè)丝跀?shù)、國(guó)內(nèi)生產(chǎn)總值、自然科學(xué)研究機(jī)構(gòu)數(shù)這4 項(xiàng)因子有2 個(gè)階段表現(xiàn)為低關(guān)聯(lián)；1988—1998 年、1998—2008 年，經(jīng)濟(jì)驅(qū)動(dòng)因素與人口驅(qū)動(dòng)因素為主要的驅(qū)動(dòng)因素，基礎(chǔ)建設(shè)驅(qū)動(dòng)因素與教育、衛(wèi)生條件驅(qū)動(dòng)因素貢獻(xiàn)雖然較低，但是穩(wěn)步增長(zhǎng)；2008—2018 年，教育、衛(wèi)生條件驅(qū)動(dòng)因素與基礎(chǔ)建設(shè)驅(qū)動(dòng)因素的貢獻(xiàn)比較高，而人口驅(qū)動(dòng)因素與經(jīng)濟(jì)驅(qū)動(dòng)因素的貢獻(xiàn)比重下降較多。

滇池作為全國(guó)第六大淡水湖、最大高原淡水湖泊，其流域內(nèi)的土地覆蓋變化進(jìn)程較為獨(dú)特。在以往的滇池流域土地利用/覆蓋變化研究中，部分學(xué)者從景觀格局[8]、地形梯度[24]等角度進(jìn)行深入探討，但未能考慮到滇池流域內(nèi)的各行政區(qū)具有較大差異。本次研究基于GEE 平臺(tái)，獲取了較高精度的長(zhǎng)時(shí)間序列的土地覆蓋結(jié)果，可為滇池流域的景觀格局優(yōu)化、生態(tài)安全建設(shè)、生態(tài)系統(tǒng)保護(hù)等相關(guān)研究提供基礎(chǔ)數(shù)據(jù)，并揭示流域內(nèi)各行政區(qū)不同的發(fā)展規(guī)律，為地方發(fā)展提供有效的科學(xué)依據(jù)。

在建設(shè)用地變化驅(qū)動(dòng)力的研究上，各行政區(qū)均有部分地區(qū)處于滇池流域之外，這導(dǎo)致數(shù)據(jù)具有不一致性，所以本研究的驅(qū)動(dòng)因子體系未能按行政區(qū)的不同來(lái)構(gòu)建，這是本研究的不足之處。但本研究仍然按照不同時(shí)間段，對(duì)不同時(shí)期內(nèi)滇池流域的建設(shè)用地進(jìn)行了驅(qū)動(dòng)力分析。在今后的研究中，將結(jié)合多模型，全面、系統(tǒng)、深入地研究如何構(gòu)建更科學(xué)的區(qū)域土地覆蓋變化驅(qū)動(dòng)體系，以期提高研究的科學(xué)性。