重慶市榮昌區(qū)土地利用變化與驅(qū)動(dòng)力研究

付偉 周寶同 田雨

摘要：為了解重慶市榮昌撤縣設(shè)區(qū)后的土地利用現(xiàn)狀和存在問(wèn)題，協(xié)調(diào)城市發(fā)展與生態(tài)保護(hù)之間的關(guān)系，基于2000年、2003年、2006年、2009年、2012年及2015年的6期遙感影像，借助遙感（RS）技術(shù)和地理信息系統(tǒng)（GIS）技術(shù)，分析榮昌區(qū)2000—2015年土地利用變化特征，并采用灰色關(guān)聯(lián)法和主成分分析法定量分析土地利用變化驅(qū)動(dòng)因素。結(jié)果表明，2000—2015年，榮昌區(qū)土地利用特征變化主要表現(xiàn)為耕地、草地、未利用地大致呈現(xiàn)減少趨勢(shì)，林地、建設(shè)用地、水域大致呈現(xiàn)增長(zhǎng)趨勢(shì)，未利用地大部分轉(zhuǎn)變?yōu)榻ㄔO(shè)用地;2000—2006年土地利用變化量為負(fù)值，2006—2015年土地利用變化量為正值，土地利用整體由調(diào)整期步入發(fā)展期;土地利用破碎程度降低，優(yōu)勢(shì)土地類型減少，土地利用景觀類型逐漸分布均勻;高關(guān)聯(lián)度指標(biāo)中城鎮(zhèn)化率、公路里程、人均地區(qū)生產(chǎn)總值等與第一主成分關(guān)聯(lián)性高，常住人口與第二主成分關(guān)聯(lián)性高，榮昌區(qū)土地利用驅(qū)動(dòng)力在研究時(shí)段以社會(huì)經(jīng)濟(jì)、交通運(yùn)輸和生態(tài)環(huán)境等因素為主。

關(guān)鍵詞：土地利用;土地利用驅(qū)動(dòng)力;榮昌區(qū);重慶;景觀指數(shù);灰色關(guān)聯(lián)度

中圖分類號(hào)：F301.2?? 文獻(xiàn)標(biāo)志碼： A? 文章編號(hào)：1002-1302（2020）06-0214-09

土地利用/覆蓋變化（LUCC）是自然環(huán)境、時(shí)空關(guān)系和經(jīng)濟(jì)技術(shù)等結(jié)構(gòu)和行為對(duì)土地利用方式與目的改變的體現(xiàn)[1]，對(duì)全球和環(huán)境變化有著重要的影響[2]，一直是地理學(xué)、生態(tài)學(xué)等相關(guān)科學(xué)領(lǐng)域研究的熱點(diǎn)[3-6]。隨著社會(huì)的發(fā)展，人類正以空前的幅度和規(guī)模影響著全球環(huán)境，造成全球人口、資源與環(huán)境問(wèn)題的日益突出[7]。因此，只有科學(xué)合理地分區(qū)利用土地，才能夠保證土地資源的節(jié)約與集約利用，從而做到各類土地的高效率利用[8]。

土地利用變化的驅(qū)動(dòng)力及驅(qū)動(dòng)機(jī)制是土地利用變化研究的重要問(wèn)題[9]。驅(qū)動(dòng)力研究是土地利用研究的核心問(wèn)題之一。在這方面，以往主要是單因子分析，目前則更多地綜合考慮多因子及其相互作用[10]。近些年，國(guó)外學(xué)者的做法是通過(guò)獲取大量的土地利用變化歷史與現(xiàn)狀數(shù)據(jù)，采用多種數(shù)學(xué)統(tǒng)計(jì)分析方法，探討引起土地利用變化的相關(guān)主導(dǎo)因子，建立土地利用變化定性的概念模型與數(shù)學(xué)模型[11]，如系統(tǒng)動(dòng)力學(xué)模型（system dynamics models）、細(xì)胞自動(dòng)機(jī)模型（cellular automata）以及混合模型（hybrid/integrated models）[12]，這些模型在一定程度上模擬了LUCC可能的原因。然而，這些模型也有一定的缺點(diǎn)，細(xì)胞自動(dòng)機(jī)模型可以模擬LUCC中生物物理方面的問(wèn)題，對(duì)人類決策的模擬是一個(gè)弱點(diǎn);系統(tǒng)動(dòng)力學(xué)模型能夠全面考慮LUCC過(guò)程中的驅(qū)動(dòng)因素，但缺乏處理空間因素的能力。國(guó)內(nèi)一些學(xué)者嘗試用Logistic回歸模型分析土地利用變化驅(qū)動(dòng)力，如謝花林運(yùn)用Logistic回歸模型分析了生態(tài)功能區(qū)內(nèi)土地利用變化的影響因素及其作用機(jī)制[13]。本研究主要通過(guò)分析土地利用3個(gè)單一因子變化，嘗試采用灰色關(guān)聯(lián)度和主成分分析方法，探究重慶市榮昌區(qū)土地利用的主要影響因素和驅(qū)動(dòng)力，以期為榮昌區(qū)未來(lái)的土地資源管理政策和制定合理的土地利用規(guī)劃提供決策依據(jù)。

1 研究地區(qū)與研究方法

1.1 研究區(qū)概況

榮昌區(qū)位于重慶市與成都市的交匯處，被稱為“渝西明珠”，是重慶市的西大門，位于105°17′～105°44′E，29°15′～29°41′N，處于成渝兩大城市的經(jīng)濟(jì)輻射范圍之內(nèi)，周圍有瀘州、內(nèi)江、自貢、宜賓、永川、江津等多個(gè)中等城市。榮昌區(qū)境內(nèi)以平壩、丘陵為主，地勢(shì)起伏不大，平均海拔為300～400 m，屬于亞熱帶季風(fēng)性氣候，全區(qū)年平均氣溫達(dá)到 17.7 ℃，年降水量為1 118 mm，年日照時(shí)數(shù)為 1 282 h，年太陽(yáng)輻射總量為3 497.84 MJ/m2。全區(qū)總面積為 1 075.37 km2，占重慶市總面積的 1.31%，轄15個(gè)鎮(zhèn)和6個(gè)街道。2015年，榮昌區(qū)的常住人口為 70.10 萬(wàn)人，城鎮(zhèn)化率為49.26%，國(guó)內(nèi)生產(chǎn)總值（GDP）為329.87億元。

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

本研究所采用的土地利用數(shù)據(jù)來(lái)源于國(guó)家地理空間數(shù)據(jù)云的2000年、2006年、2009年 Landsat-5 TM和2012年、2015年Landsat-8 OLI遙感影像，同時(shí)搜集了榮昌區(qū)2000—2015年各種社會(huì)經(jīng)濟(jì)指標(biāo)、行政界限、研究期土地利用現(xiàn)狀數(shù)據(jù)等資料。

在ENVI 5.1軟件的支持下，對(duì)榮昌區(qū)6期遙感影像按行政邊界裁剪后，采用幾何精校正、大氣校正以及圖像鑲嵌、裁剪和增強(qiáng)的方式得到了6期遙感解譯圖。根據(jù)GB/T 21010—2007《土地利用現(xiàn)狀分類》和中國(guó)科學(xué)院土地利用遙感檢測(cè)分類系統(tǒng)[14]，采取監(jiān)督與非監(jiān)督分類法和目視解譯相結(jié)合的方法進(jìn)行分類，將榮昌區(qū)的土地利用類型劃分為6類：林地、耕地、草地、水域、建設(shè)用地、未利用地，生成榮昌區(qū)2000—2015年6期土地利用圖（圖1至圖6）。為保證分類結(jié)果的準(zhǔn)確性，通過(guò)Kappa系數(shù)法對(duì)榮昌區(qū)6期遙感解譯分類結(jié)果進(jìn)行精度評(píng)價(jià)，得到Kappa系數(shù)依次為0.826、0.814、0.806、0.833、0.815、0.822，評(píng)價(jià)結(jié)果精度均在80%以上，說(shuō)明可以用于后續(xù)研究。

1.3 研究方法

1.3.1 土地利用變化動(dòng)態(tài)度指數(shù) 單一土地利用類型動(dòng)態(tài)度表達(dá)的是某一研究區(qū)一定時(shí)間范圍內(nèi)某種土地利用類型的數(shù)量變化情況，在本研究中表達(dá)為年變化率，其表達(dá)式[15-16]為式中：K為研究時(shí)段內(nèi)某類土地利用類型的年變化率;Ua、Ub分別為研究期初和研究期末某類土地利用類型的面積;T為研究年段長(zhǎng)度。

1.3.2 土地利用變化程度指數(shù) 土地利用程度反映土地利用的廣度和深度， 反映人類活動(dòng)和自然環(huán)境對(duì)土地利用的綜合作用。土地利用程度變化量可表達(dá)為式中：La和Lb分別為研究初期和研究末期的土地利用程度指數(shù);Ai為第i級(jí)土地利用程度分級(jí)指數(shù);Cib和Cia分別為b時(shí)間和a時(shí)間第i級(jí)土地利用程度所占區(qū)域面積百分比。若ΔLb-a>0，則表明該地區(qū)土地利用處于發(fā)展期;若ΔLb-a<0，則表明該地區(qū)土地利用處于調(diào)整期或衰退期。

1.3.3 土地利用變化景觀指數(shù) 在景觀生態(tài)學(xué)中，

景觀指數(shù)的變化是分析景觀要素動(dòng)態(tài)變化的重要指標(biāo)[17]。土地利用是由各種不同類型土地利用面積也就是斑塊組成，空間特征和時(shí)間特征顯著。本研究從景觀水平上選取最大斑塊指數(shù)（LPI）、斑塊個(gè)數(shù)（NP）、景觀形狀指數(shù)（LSI）、香農(nóng)多樣性指數(shù)（SHDI）、景觀均勻度指數(shù)（SHEI）等5個(gè)景觀指數(shù)來(lái)研究榮昌區(qū)土地利用變化情況，具體指標(biāo)見(jiàn)表1。

1.3.4 灰色關(guān)聯(lián)分析法和主成分分析法 如果系統(tǒng)結(jié)構(gòu)不明確，作用原理難以明確闡述，但對(duì)系統(tǒng)的最后結(jié)果總能得到一些資料和信息，這種信息部分明確、部分不明確的抽象系統(tǒng)稱之為灰色系統(tǒng)[18]。本研究共選取常住人口、城鎮(zhèn)化率、公路里程、地區(qū)生產(chǎn)總值、工業(yè)總產(chǎn)值、人均地區(qū)生產(chǎn)總值、畜牧業(yè)總產(chǎn)值、森林覆蓋率等27個(gè)指標(biāo)來(lái)建立榮昌區(qū)土地利用驅(qū)動(dòng)力體系，由于指標(biāo)較多樣本較小，在進(jìn)行灰色關(guān)聯(lián)度分析和主成分分析之前，對(duì)榮昌區(qū)2000—2015年土地利用變化程度和這27個(gè)指標(biāo)進(jìn)行常規(guī)擬合，然后對(duì)這27個(gè)指標(biāo)和土地利用綜合指數(shù)進(jìn)行灰色關(guān)聯(lián)定量分析。本研究將榮昌區(qū)土地利用變化綜合指數(shù)設(shè)為參考序列，即母序列，將27個(gè)指標(biāo)設(shè)為比較序列，即子序列。根據(jù)步驟計(jì)算得出27個(gè)驅(qū)動(dòng)因子變量的灰色關(guān)聯(lián)度，從而得到灰色關(guān)聯(lián)序。

通過(guò)灰色關(guān)聯(lián)法分析得到27個(gè)指標(biāo)與土地利用變化綜合指數(shù)的關(guān)聯(lián)度（R），通過(guò)觀察灰色關(guān)聯(lián)度結(jié)果并綜合多方面考慮，挑選出10個(gè)關(guān)聯(lián)度高且典型的驅(qū)動(dòng)指標(biāo)來(lái)進(jìn)行下一步的主成分分析。

2 結(jié)果與分析

2.1 土地利用動(dòng)態(tài)度變化

由表2可知，在研究期的15年間，耕地、草地、未利用地的面積呈現(xiàn)減少趨勢(shì)，林地、水域、建設(shè)用地面積總體呈增加趨勢(shì)。由于退耕還林還湖政策的開展，林地面積大幅度增加，變化率達(dá)11.08%，面積增加 11 958.53 hm2，耕地面積減少 14 924.42 hm2。受城鎮(zhèn)化、社會(huì)經(jīng)濟(jì)水平提高等因素的影響，建設(shè)用地面積變化較明顯，增加 4 888.98 hm2。草地、水域、未利用地面積較小，其中草地和未利用地面積分別減少701.24、1 231.90 hm2，水域面積增加144.35 hm2。榮昌區(qū)各類型土地面積比例都有所變化，但一直以耕地、林地和建設(shè)用地為主，其中耕地在這3種類型土地面積中占主導(dǎo)地位，這主要和榮昌區(qū)工業(yè)化水平較低、以第一產(chǎn)業(yè)生產(chǎn)為主等有關(guān)。

從表3可以看出，2000—2015年林地的年變化率最高，為 11.48%;建設(shè)用地的年變化率為 2.53%;未利用地的減少速度最快，年變化率為 -6.32%;草地的變化速度較未利用地慢，但年變化

率達(dá)到 -5.28%;耕地和水域的變化與其他地類相比不是特別明顯，耕地的年變化率為-1.21%，水域的年變化率為1.30%。林地增加的主要原因除了退耕還林外，還和榮昌區(qū)積極實(shí)施荒山造林、水系綠化和道路綠化有關(guān)。建設(shè)用地一直處于增長(zhǎng)狀態(tài)，2006—2009年變化最快，年變化率為7.37%，說(shuō)明這一階段是榮昌區(qū)社會(huì)經(jīng)濟(jì)發(fā)展較好較快的階段。未利用地的年變化率在5個(gè)階段均為負(fù)值，說(shuō)明未利用地一直在減少，2006—2009年快速減少，年變化率為-30.55%，建設(shè)用地在這一階段增長(zhǎng)最快，而未利用地在這一階段快速減少，說(shuō)明未利用地大部分轉(zhuǎn)為了建設(shè)用地，服務(wù)了社會(huì)經(jīng)濟(jì)的發(fā)展。

2.2 土地利用程度變化

由表4可知，榮昌區(qū)2000—2015年5個(gè)階段土地利用變化量有正有負(fù)，2000—2003年ΔLb-a為 -5.35;2003—2006年ΔLb-a為-0.81;2006—2009年ΔLb-a為1.95;2009—2012年ΔLb-a為0.64;2012—2015年ΔLb-a為0.88。同樣，榮昌區(qū)土地利用程度變化率也有正有負(fù)，2000—2003年為 -1.82%;2003—2006年為-0.28%;2006—2009年為0.66%;2009—2012年為0.22%;2012—2015年為0.30%。榮昌區(qū)2000—2003、2003—2006年的土地利用程度變化量和土地利用程度變化率的值均為負(fù)，說(shuō)明這2個(gè)階段榮昌區(qū)處于調(diào)整期，人類活動(dòng)對(duì)土地利用的程度較低，對(duì)土地來(lái)說(shuō)，土地的負(fù)載較小，自然狀態(tài)較好;后3個(gè)階段土地利用程度變化量和土地利用程度變化率均為正值，說(shuō)明這3個(gè)階段榮昌區(qū)處于發(fā)展期，人類活動(dòng)對(duì)土地的利用程度在提高，土地的產(chǎn)值不斷被挖掘，土地的負(fù)荷較大。2000—2015年榮昌區(qū)的土地利用程度變化量為 -2.69，雖然后3個(gè)階段土地利用程度變化量為正值，但數(shù)值較小，導(dǎo)致榮昌區(qū)這15年以-0.91%的速度在調(diào)整發(fā)展，從后3個(gè)階段的土地利用程度變化量和土地利用程度變化率來(lái)看，榮昌區(qū)會(huì)慢慢從調(diào)整期步入到發(fā)展期，這也和榮昌區(qū)近些年城鎮(zhèn)化水平提高、工業(yè)化不斷發(fā)展、土地得到了充分利用相符合。

2.3 土地利用景觀變化

從表5可以看出，各景觀指數(shù)在2000—2015年間都發(fā)生了變化。斑塊個(gè)數(shù)（NP）逐漸下降，在2012—2015年間下降速度最快，說(shuō)明榮昌區(qū)土地利用景觀破碎度逐漸降低，土地利用類型分布逐漸集中;最大斑塊指數(shù)（LPI）逐漸降低并趨向平穩(wěn)，說(shuō)明榮昌區(qū)內(nèi)的優(yōu)勢(shì)土地利用類型面積在逐漸減少，根據(jù)土地利用數(shù)量分析可知，耕地是榮昌區(qū)內(nèi)的優(yōu)勢(shì)土地利用類型，說(shuō)明耕地面積在逐漸減少，土地利用類型在逐漸變得豐富;景觀形狀指數(shù)（LSI）在2012—2015年間下降速度最快，但最高值和最低值之間差距不大，說(shuō)明榮昌區(qū)的土地利用類型面積變化不大，規(guī)則程度較穩(wěn)定，受人為管控制約較強(qiáng);香農(nóng)多樣性指數(shù)（SHDI）和景觀均勻度指數(shù)（SHEI）在這15年均呈現(xiàn)緩慢上升趨勢(shì)，說(shuō)明榮昌區(qū)土地利用景觀類型在逐漸變得豐富，土地利用景觀類型分布逐漸均勻。

2.4 灰色關(guān)聯(lián)和主成分分析

本研究將選取的27個(gè)指標(biāo)分別與土地利用程度進(jìn)行常規(guī)擬合，從擬合程度的高低來(lái)判斷驅(qū)動(dòng)指標(biāo)對(duì)土地利用程度的影響，將擬合結(jié)果r2>0.70的定為高度擬合程度指標(biāo)，將r2<0.30的定為低度擬合程度指標(biāo)，將0.30≤r2≤0.70的指標(biāo)定為中度擬合程度指標(biāo)（r2的取值為0～1，越接近1說(shuō)明擬合程度越好）。如表6所示，高度擬合指標(biāo)包括常住人口、城鎮(zhèn)化率、公路里程、人均地區(qū)生產(chǎn)總值、油料產(chǎn)量、地區(qū)生產(chǎn)總值、工業(yè)總產(chǎn)值等，這些指標(biāo)均與土地利用程度指數(shù)有密切的相關(guān)性，說(shuō)明在榮昌區(qū)撤縣設(shè)區(qū)后，人口及城鎮(zhèn)化率因素、社會(huì)經(jīng)濟(jì)因素、農(nóng)村發(fā)展因素、科技文化因素和生態(tài)因素對(duì)土地利用程度影響較大;中度擬合指標(biāo)有豬肉產(chǎn)量、農(nóng)用化肥施用量、普通中學(xué)在校學(xué)生數(shù)、森林覆蓋率、工業(yè)固體廢物排放量等，說(shuō)明這幾個(gè)指標(biāo)對(duì)榮昌區(qū)土地利用程度變化的影響一般。低度擬合指標(biāo)有失業(yè)人員登記數(shù)、糧食產(chǎn)量、建筑企業(yè)數(shù)、農(nóng)作物播種面積等，這些指標(biāo)對(duì)土地利用程度的影響較小。

由土地利用動(dòng)態(tài)度、土地利用程度指數(shù)和土地利用景觀變化指數(shù)構(gòu)成土地利用變化綜合指數(shù)，該綜合指數(shù)對(duì)3個(gè)指標(biāo)用變異系數(shù)法設(shè)置權(quán)重并通過(guò)加權(quán)求和法算出，將綜合指數(shù)設(shè)為參考序列即母序列，將27個(gè)指標(biāo)設(shè)為比較序列即子序列。根據(jù)步驟計(jì)算得出27個(gè)驅(qū)動(dòng)因子變量的灰色關(guān)聯(lián)度，從而得到灰色關(guān)聯(lián)序。27個(gè)驅(qū)動(dòng)指標(biāo)和土地利用變化綜合指數(shù)的關(guān)聯(lián)度具體情況見(jiàn)表7。

通過(guò)觀察得到的27個(gè)指標(biāo)與土地利用變化綜合指數(shù)的關(guān)聯(lián)度，發(fā)現(xiàn)10個(gè)關(guān)聯(lián)度高且典型的驅(qū)動(dòng)指標(biāo)為X1常住人口（0.999 3）、X2城鎮(zhèn)化率（0.983 9）、X3公路里程（0.927 0）、X9農(nóng)林牧漁總產(chǎn)值（0.963 6）、X14畜牧業(yè)總產(chǎn)值（0.964 9）、X15人均地區(qū)生產(chǎn)總值（0.903 9）、X19金融機(jī)構(gòu)人民幣存款余額（0.917 3）、X22社會(huì)消費(fèi)品零售總額（0.946 7）、X24工業(yè)SO2排放量（0.996 5）、X26人均園林綠地面積（0.943 4）。由以上結(jié)果可以看出，榮昌區(qū)撤縣設(shè)區(qū)后，人口增長(zhǎng)和城鎮(zhèn)化發(fā)展都促使榮昌區(qū)土地利用發(fā)生轉(zhuǎn)變，這與榮昌區(qū)常住人口從2000年的65.6萬(wàn)人增加到2015年的70.1萬(wàn)人，戶籍人口數(shù)從2000年的80.64萬(wàn)人增加到2015年的84.48萬(wàn)人，榮昌區(qū)的城鎮(zhèn)化率從2000年的21.1%增加到2015年的49.26%密切相關(guān);地區(qū)生產(chǎn)總值、人均地區(qū)生產(chǎn)總值的增長(zhǎng)直接說(shuō)明一個(gè)地區(qū)的經(jīng)濟(jì)發(fā)展水平，隨著榮昌區(qū)社會(huì)經(jīng)濟(jì)水平的發(fā)展，榮昌區(qū)人民的生活水平明顯提高，人們對(duì)環(huán)境質(zhì)量要求越來(lái)越高，社會(huì)生活環(huán)境得到顯著改善;在2000—2015年這15年里，榮昌區(qū)公路里程不斷延長(zhǎng)，說(shuō)明榮昌區(qū)的交通運(yùn)輸條件在不斷改善，榮昌區(qū)的通達(dá)性和可進(jìn)入性明顯得到了提高; 榮昌區(qū)的特色農(nóng)牧業(yè)發(fā)展也促使了土地利用方式的轉(zhuǎn)變。

以這10個(gè)指標(biāo)作為變量，借助SPSS 21.0進(jìn)行分析研究，從表8可以得知，X2城鎮(zhèn)化率、X3公路里程、X9農(nóng)林牧漁總產(chǎn)值、X14畜牧業(yè)、X15人均地區(qū)生產(chǎn)總值、X19金融機(jī)構(gòu)人民幣存款余額、X22社會(huì)消費(fèi)品零售總額、X26人均園林綠地面積和第1主成分的相關(guān)系數(shù)較高;X1常住人口和第2主成分的關(guān)聯(lián)度較高，由此可將第1主成分確定為社會(huì)經(jīng)濟(jì)、基礎(chǔ)設(shè)施因素，人均園林綠地面積可反映榮昌區(qū)的生態(tài)環(huán)境，第2主成分和常住人口相關(guān)性密切，由此看來(lái)人口情況對(duì)榮昌區(qū)土地利用影響較大。

3 結(jié)論

本研究先借助遙感（RS）和地理信息系統(tǒng)（GIS）技術(shù)獲得遙感解譯數(shù)據(jù)，通過(guò)對(duì)土地利用變化動(dòng)態(tài)度、土地利用程度，景觀變化度等單一指標(biāo)進(jìn)行分析，結(jié)果表明， 從數(shù)量變化上看，在研究階段期間，6類土地利用類型均發(fā)生了不同程度的變化。耕地、草地、未利用地的面積呈減少趨勢(shì)，其中耕地的面積減少的最多，2000—2015年間共減少 14 924.4 hm2;林地、建設(shè)用地、水域的面積呈增加趨勢(shì)，林地的面積增加最多，共增加11 958.53 hm2，建設(shè)用地面積共增加4 888.98 hm2;研究期間，榮昌區(qū)土地利用方式以耕地、林地、建設(shè)用地為主，未利用地主要轉(zhuǎn)變?yōu)榻ㄔO(shè)用地。從土地利用程度變化上看，在2000—2015年間，土地利用程度變化量和變化率都小于0，土地利用程度總體屬于中上水平;研究前期榮昌區(qū)工業(yè)化基礎(chǔ)薄弱，土地利用處于調(diào)整期，人類活動(dòng)對(duì)土地利用程度較低，隨著榮昌區(qū)社會(huì)經(jīng)濟(jì)水平的提高，基礎(chǔ)設(shè)施的完善，后期土地利用逐漸步入發(fā)展期。從土地利用景觀變化上看，榮昌區(qū)斑塊個(gè)數(shù)、最大斑塊指數(shù)、景觀形狀指數(shù)逐漸下降，景觀多樣性和景觀均勻度緩慢上升，榮昌區(qū)土地利用格局逐漸均衡發(fā)展，草地、水域、未利用地所占比例較小，主要是耕地、林地、建設(shè)用地均衡發(fā)展。

榮昌區(qū)撤縣設(shè)區(qū)后，城市化水平不斷推進(jìn)，土地利用方式發(fā)生了很大轉(zhuǎn)變。在對(duì)27個(gè)驅(qū)動(dòng)指標(biāo)中10個(gè)關(guān)聯(lián)度指標(biāo)的分析中可以看出，城鎮(zhèn)化率水平的提高，迫使區(qū)域土地利用發(fā)生動(dòng)態(tài)調(diào)整，城鎮(zhèn)用地的擴(kuò)張，耕地和未利用地等自然、半自然景觀不斷轉(zhuǎn)變?yōu)榻ㄔO(shè)用地等，給生態(tài)環(huán)境造成了較大的壓力;社會(huì)經(jīng)濟(jì)水平的提高，人均地區(qū)生產(chǎn)總值的增長(zhǎng)使人們對(duì)環(huán)境提出了更高的要求，耕地一部分轉(zhuǎn)化為城市綠地，另外退耕還林政策使得耕地面積進(jìn)一步減少，將影響到區(qū)域糧食產(chǎn)量。因此，應(yīng)正確處理好發(fā)展用地、耕地保護(hù)和生態(tài)保護(hù)之間的關(guān)系，合理保護(hù)耕地，加大對(duì)未利用地的開發(fā)利用，同時(shí)注意對(duì)水域等生態(tài)景觀的保護(hù);明確榮昌區(qū)產(chǎn)業(yè)發(fā)展方向和用地標(biāo)準(zhǔn)，注重產(chǎn)業(yè)布局的規(guī)?；图夯?，節(jié)約集約利用土地，發(fā)揮第一特色產(chǎn)業(yè)，注重與第二、第三產(chǎn)業(yè)協(xié)調(diào)發(fā)展，優(yōu)化利用土地資源。

本研究通過(guò)研究土地利用變化的3個(gè)指數(shù)，分析榮昌區(qū)2000—2015年的土地利用變化情況，利用灰色關(guān)聯(lián)和主成分分析法來(lái)探究土地利用驅(qū)動(dòng)力因素，對(duì)榮昌區(qū)的土地利用及轉(zhuǎn)變作了系統(tǒng)的分析，對(duì)于榮昌區(qū)未來(lái)的土地規(guī)劃和決策具有一定的借鑒作用，但文中土地利用變更數(shù)據(jù)的獲取以及分析指標(biāo)選取的合理性還需要進(jìn)一步研究。由于數(shù)據(jù)獲取的有限性和遙感影像本身分辨率的問(wèn)題，最終解譯結(jié)果可能會(huì)存在差異，因此，盡可能地結(jié)合多種解譯方法從而提高解譯精度是必要的。

參考文獻(xiàn)：

[1]Overmars K P，Verburg P H.Analysis of land use drivers at the watershed and household level：linking two paradigms at the Philippine forest fringe[J]. International Journal of Geographical Information Science，2005，19（2）：125-152.

[2]侯 鵬，王 橋，王昌佐，等. 流域土地利用/土地覆被變化的生態(tài)效應(yīng)[J]. 地理研究，2011，30（11）：2092-2098.

[3]秦 鵬，董玉祥. 廣州市土地利用空間分布影響因素的粒度效應(yīng)分析[J]. 資源科學(xué)，2013，35（11）：2239-2247.

[4]張 舟，吳次芳，譚 榮. 生態(tài)系統(tǒng)服務(wù)價(jià)值在土地利用變化研究中的應(yīng)用：瓶頸和展望[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2013，24（2）：556-562.

[5]劉金勇，孔繁花，尹海偉，等. 濟(jì)南市土地利用變化及其對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2013，24（5）：1231-1236.

[6]劉紀(jì)遠(yuǎn)，匡文慧，張?jiān)鱿椋? 20世紀(jì)80年代末以來(lái)中國(guó)土地利用變化的基本特征與空間格局[J]. 地理學(xué)報(bào)，2014，69（1）：3-14.

[7]馬 靜. 基于GIS的土地利用/覆被變化及其社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)力研究——以張掖市為例[D]. 蘭州：西北師范大學(xué)，2008.

[8]王 炯，許月明，郭 慶. 基于聚類分析法的保定市土地利用分區(qū)及建議[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃，2011，32（2）：63-67.

[9]李秀彬. 全球環(huán)境變化研究的核心領(lǐng)域——土地利用/土地覆被變化的國(guó)際研究動(dòng)向[J]. 地理學(xué)報(bào)，1996，51（6）：553-557.

[10]黃秉維. 論地球系統(tǒng)科學(xué)與可持續(xù)發(fā)展戰(zhàn)略科學(xué)基礎(chǔ)（Ⅰ）[J]. 地理學(xué)報(bào)，1996，63（4）：350-354.

[11]Bi cˇík I，Jele cˇek L，Sˇtěpánek V. Land-use changes and their social driving forces in Czechia in the 19th and 20th centuries[J]. Land Use Police，2001，18（1）：65-73.

[12]韓超峰，陳仲新. LUCC驅(qū)動(dòng)力模型研究綜述[J]. 中國(guó)農(nóng)學(xué)通報(bào)，2008，24（4）：365-368.

[13]謝花林. 基于Logistic回歸模型的區(qū)域生態(tài)用地演變影響因素分析——以京津冀地區(qū)為例[J]. 資源科學(xué)，2011，33（11）：2063-2070.

[14]史 舟，梁宗正，楊媛媛，等. 農(nóng)業(yè)遙感研究現(xiàn)狀與展望[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2015，46（2）：247-260.

[15]李忠峰，王一謀，馮毓蓀，等. 基于RS與GIS的榆林地區(qū)土地利用變化分析[J]. 水土保持學(xué)報(bào)，2003，17（2）：97-99.

[16]王美玲，邴龍飛，郗鳳明，等. 老工業(yè)搬遷區(qū)土地利用變化時(shí)空特征及其驅(qū)動(dòng)力——以沈陽(yáng)市鐵西老工業(yè)區(qū)為例[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2013，24（7）：1969-1979.

[17]陸華麗. RS和GIS支持下辛集市土地利用/土地覆蓋變化監(jiān)測(cè)研究[D]. 保定：河北農(nóng)業(yè)大學(xué)，2008.

[18]張軍良. 利用灰色控制系統(tǒng)模型進(jìn)行預(yù)測(cè)分析研究[J]. 計(jì)算機(jī)與網(wǎng)絡(luò)，2012，38（11）：69-71.王錦杰，陳 昊，張 瑩，等. 基于植被健康指數(shù)的2001—2018年間江蘇省農(nóng)業(yè)干旱時(shí)空分析[J]. 江蘇農(nóng)業(yè)科學(xué)，2020，48（6）：223-231.