基于遙感影像的2018—2020年肥城市土地利用變化分析

摘 要：隨著經(jīng)濟的快速發(fā)展和城鎮(zhèn)化進程的不斷推進，肥城市的土地利用格局發(fā)生了較大變化。精準把握肥城市土地利用變化特征，對優(yōu)化土地利用空間格局具有重要意義?；?018—2020年肥城市的Landsat 8 OLI_TIRS遙感影像數(shù)據(jù)，采用土地利用動態(tài)度、土地利用類型轉(zhuǎn)移矩陣、重心遷移模型等方法對肥城市土地利用變化狀況進行分析。結(jié)果表明：2018—2020年，肥城市土地利用類型以耕地、城鄉(xiāng)建設(shè)用地、林地為主；2018—2020年，肥城市土地利用類型變化、時空變化相對較快，綜合土地利用動態(tài)度為8.73%，其中土地利用動態(tài)度最大的地類是城鄉(xiāng)建設(shè)用地，最小的是草地；2018—2020年，肥城市城鄉(xiāng)建設(shè)用地、林地、水域、未利用地面積擴大，而耕地、草地面積縮小。

關(guān)鍵詞：GIS；土地利用變化；肥城市

中圖分類號：F321.1；P237 文獻標志碼：A 文章編號：1674-7909（2024）5-122-5

DOI：10.19345/j.cnki.1674-7909.2024.05.028

0 引言

土地利用是指人類開展的所有有目的的開發(fā)利用土地資源的行為［1-2］。分析土地利用的時空變化是各領(lǐng)域的重要研究內(nèi)容，能為土地資源合理利用提供重要的科學依據(jù)［3］。目前，諸多學者已對土地利用變化進行了研究。陳伊多等［4］基于各期Landsat TM/ETM遙感影像，采用土地利用動態(tài)度、轉(zhuǎn)化趨勢等方法分析了西藏自治區(qū)土地利用/覆被特征。馮應(yīng)斌等［5］綜合運用空間自相關(guān)模型、地理探測器等，分析了貴州省耕地顯性邊際化的總體態(tài)勢、空間集聚效應(yīng)及相關(guān)影響因素。廖美玉等［6］采用Mann-Kendall、滑動t檢驗、標準差橢圓和重心遷移模型等方法，探究了近40年黃河流域土地利用/覆被時空變化特征。張晨基等［7］在采用土地利用動態(tài)度及土地利用轉(zhuǎn)移矩陣方法的基礎(chǔ)上，疊合DEM高程數(shù)據(jù)，對新疆土地利用時空變化進行研究。安彬［8］采用土地利用程度綜合指數(shù)、土地利用信息圖譜等方法，研究了1985—2020年漢江生態(tài)經(jīng)濟帶土地利用/覆蓋演變特征。以上研究多從大的時間和空間尺度上對土地利用進行分析，不利于捕捉小的時間和空間尺度上的土地利用時空變化?；诖耍P者以肥城市為研究區(qū)域，基于2018—2020年的土地利用狀況，探究該區(qū)域的時空變化特征，以期為肥城市土地利用空間規(guī)劃提供理論依據(jù)。

1 研究區(qū)概況

肥城市隸屬于山東省泰安市，地處山東中部、泰山西麓，位于北緯35°53′～36°19′，東經(jīng)116°28′～116°59′；地勢北高南低，由東北向西南傾斜；中部隆起地帶為丘陵山地，北部是康匯平原，南部是汶陽平原，自然形成了平原、山地、丘陵等多種地貌形態(tài)［9］；屬溫帶季風氣候區(qū)，四季分明，雨熱同期，光照充足，氣候溫暖，年平均降水量為659 mm。肥城市礦產(chǎn)資源豐富，礦業(yè)勘探開發(fā)歷史悠久，探明資源量大。目前，肥城市已發(fā)現(xiàn)礦產(chǎn)29種，占山東省已發(fā)現(xiàn)礦產(chǎn)總數(shù)的47.6%。

2 研究數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源

2018—2020年肥城市土地利用遙感監(jiān)測數(shù)據(jù)來自地理空間數(shù)據(jù)云。各年份土地利用類型數(shù)據(jù)通過對各期Landsat 8 OLI_TIRS 遙感影像（分辨率為30 m）人工目視解譯生成，利用ENVI 5.6裁剪、輻射定標、大氣校正等預(yù)處理，在ArcGIS 10.2和Excel中進行處理統(tǒng)計。

2.2 研究方法

2.2.1 土地利用動態(tài)度

土地利用動態(tài)度反映一定時期內(nèi)某一區(qū)域不同土地利用類型的數(shù)量變化，一般分為單一土地利用動態(tài)度和綜合土地利用動態(tài)度［10］。

2.2.1.1 單一土地利用動態(tài)度

單一土地利用動態(tài)度表示一定時期內(nèi)某一區(qū)域某一種土地利用類型的數(shù)量變化情況［11］，其表達式見式（1）。

式（1）中：Ki為研究區(qū)域內(nèi)第i類土地利用類型在a～b時段內(nèi)的土地利用動態(tài)度；Uia、Uib分別為a、b時段第i類土地類型的面積；Ta、Tb分別為研究初期時間與研究末期時間。

Ki值越大，表示該土地利用類型轉(zhuǎn)化為其他類型土地的數(shù)量越多，即在研究時段內(nèi)相對變化程度較大；反之，轉(zhuǎn)化數(shù)量越少，在研究時段內(nèi)相對變化程度較小。

2.2.1.2 綜合土地利用動態(tài)度

綜合土地利用動態(tài)度表示一定時期內(nèi)某一區(qū)域不同土地利用類型之間的數(shù)量變化情況［12］，其表達式見式（2）。

式（2）中：LC為一定時期內(nèi)綜合土地利用動態(tài)度；LUi為研究初期第i類土地利用類型的面積；ΔLUi為研究期內(nèi)第i類土地利用類型轉(zhuǎn)化為非i類土地利用類型面積的絕對值；T為研究時段。

LC值越大，表示該區(qū)域土地利用類型的變化程度越大；反之，表示該區(qū)域土地利用類型的變化程度越小。

2.2.2 土地利用類型轉(zhuǎn)移矩陣

土地利用類型轉(zhuǎn)移矩陣以矩陣的形式將2個不同時期土地利用類型之間互相轉(zhuǎn)換的數(shù)量關(guān)系展現(xiàn)出來，是具體反映土地利用變化的結(jié)構(gòu)特征和各類型之間轉(zhuǎn)移方向的重要方法［13］。其表達式見式（3）。

式（3）中：Sij表示從研究初期到研究末期土地類型由i轉(zhuǎn)化為j的土地面積；n表示土地利用類型的數(shù)量；i、j分別為研究初期、研究末期的土地利用類型。

2.2.3 重心遷移模型

重心遷移模型用于反映研究區(qū)域內(nèi)各土地利用類型在時空分布上的動態(tài)變化，可以精確把握土地利用類型的空間遷移趨勢［14］。其表達式見式（4）。

式（4）中：Xt和Yt分別為t時期某種土地利用類型質(zhì)心的經(jīng)度和緯度坐標；Cti為t時期某種土地利用類型第i個斑塊的面積；Xti和Yti分別為t時期土地利用類型第i個斑塊的幾何重心經(jīng)度和緯度坐標。

3 結(jié)果分析

3.1 土地利用分類結(jié)果

研究利用ENVI 5.6與ArcGIS 10.2，根據(jù)土地利用分級系統(tǒng)，對2018年、2019年、2020年肥城市3期土地利用遙感影像進行目視解譯，得到3個時間節(jié)點的土地利用類型分類格局。

肥城市土地利用類型統(tǒng)計情況見表1。由表1可知，在肥城市的土地利用構(gòu)成中，耕地、城鄉(xiāng)建設(shè)用地和林地是三大主導(dǎo)類型，這3種類型的土地總面積占據(jù)了全市土地總面積的90%以上。在所有土地利用類型中，耕地面積所占比重一直最大但呈逐年遞減的趨勢，城鄉(xiāng)建設(shè)用地面積逐年增加，林地面積整體呈上升趨勢，草地、水域及未利用地面積占比較小且波動較小。

具體到空間分布上，肥城市的耕地面積最大，主要分布在肥城市北部的康匯平原和肥城市南部的汶陽平原，而在肥城市中部的丘陵地帶則有相對零星分布；城鄉(xiāng)建設(shè)用地在肥城市北部地區(qū)更為密集，而在肥城市南部較為分散；林地主要集中分布在肥城市北部山區(qū)和肥城市中部丘陵地帶；草地多分布在肥城市北部山麓和肥城市中部的丘陵地帶；水域則主要分布在肥城市南部的大汶河、肥城市中部的康王河及肥城市北部的礦區(qū)塌陷區(qū)；未利用地在肥城市極為少見，只在肥城市北部山區(qū)的某些區(qū)域有少量分布。

3.2 土地利用時間變化

2018—2020年肥城市土地利用動態(tài)度見表2。由表2可知，2018—2020年，土地利用動態(tài)度最大的是城鄉(xiāng)建設(shè)用地，最小的是草地，說明在此期間城鄉(xiāng)建設(shè)用地增加的幅度和草地減少的幅度最為明顯；耕地的土地利用動態(tài)度小于０，表明耕地面積逐年減少，這種現(xiàn)象可能是由城市化進程加快造成的；林地、水域、草地、未利用地的土地利用動態(tài)度均從小于0到大于0，表明其面積均為先減少后增加，2019年后這4類土地利用類型面積增大，可能與人類工農(nóng)業(yè)活動、出行減少有關(guān)；城鄉(xiāng)建設(shè)用地的土地利用動態(tài)度均大于0，表明城鄉(xiāng)建設(shè)用地的土地利用面積持續(xù)增加，說明2018—2020年肥城市經(jīng)濟發(fā)展速度較快、人類活動強度加劇。2018—2020年肥城市的綜合土地利用動態(tài)度為8.73%，綜合土地利用動態(tài)度大于2.50%，說明土地利用變化為快速變動型［15］，整體土地利用類型數(shù)量變化相對較快。

3.3 土地利用空間變化

3.3.1 土地利用轉(zhuǎn)移矩陣

2018—2019年、2019—2020年肥城市土地利用類型轉(zhuǎn)移矩陣分別見表3和表4。由表3和表4可知，2018—2020年，耕地的流出面積最多，流出面積為184.81 km2，主要流向了城鄉(xiāng)建設(shè)用地和林地；林地的流出面積次之，流出面積為47.13 km2，主要流向了耕地。流入面積最多的是城鄉(xiāng)建設(shè)用地，流入面積為110.18 km2，主要來源于耕地和林地；林地的流入面積次之，流入面積為85.13 km2，主要來源于耕地和草地。

總體而言，肥城市土地利用變化表現(xiàn)為城鄉(xiāng)建設(shè)用地、林地和水域的面積增加，而耕地、草地和未利用地的面積減少。這一變化主要是因為城市化進程加快，城鄉(xiāng)邊界不斷擴展，城市化水平逐漸提高。同時，退耕還林政策的實施和礦區(qū)塌陷區(qū)的植被恢復(fù)，尤其是速生楊的種植，也是林地面積增加的重要原因；礦區(qū)塌陷形成的水體導(dǎo)致水域面積增加。這些變化反映了肥城市土地利用的動態(tài)調(diào)整既是城市化發(fā)展的必然結(jié)果，也是生態(tài)恢復(fù)和環(huán)境治理政策實施的具體體現(xiàn)。

3.3.2 重心遷移

2018—2020年肥城市土地利用類型的重心變化情況見表5。由表5可知，肥城市耕地重心向肥城市東南方向遷移，城鄉(xiāng)建設(shè)用地重心向肥城市西北方向遷移，二者遷移方向相反，主要是因為城市擴張占用大量耕地資源、消耗大量水資源，西北部地形平坦且有河流經(jīng)過，適宜城市建設(shè)，因此城市向西北方向遷移。

4 結(jié)束語

筆者根據(jù)地理空間數(shù)據(jù)云的遙感監(jiān)測數(shù)據(jù)，采用土地利用動態(tài)度、土地利用類型轉(zhuǎn)移矩陣、重心遷移模型等方法，并借助ArcGIS10.2、ENVI5.6和Excel，對2018—2020年肥城市土地利用類型時空變化特征進行研究，得出以下結(jié)論。

①2018—2020年，肥城市土地利用類型以耕地、城鄉(xiāng)建設(shè)用地、林地為主；肥城市的綜合土地利用動態(tài)度為8.73%，整體土地利用類型數(shù)量變化相對較快，多數(shù)土地利用類型的單一動態(tài)度數(shù)量變化相對劇烈，與城市快速擴張有關(guān)。因此，政府部門應(yīng)重視土地利用空間規(guī)劃，合理規(guī)劃利用土地，確保城市可持續(xù)發(fā)展。

②2018—2020年，肥城市城鄉(xiāng)建設(shè)用地、林地、水域、未利用地面積擴大，耕地、草地面積縮?。桓刂匦南蚍食鞘袞|南方向遷移，城鄉(xiāng)建設(shè)用地重心向肥城市西北方向遷移。因此，政府部門應(yīng)加強對非農(nóng)業(yè)用地侵占基本農(nóng)田的管理，減少耕地非農(nóng)化、非糧化，確保耕地面積合理穩(wěn)定。

此研究以縣級市為例，從多種土地利用類型的綜合利用到單一土地利用類型的綜合利用的變化，綜合研究了肥城市土地利用類型特征，彌補了從大的時間和空間尺度上研究的局限性與片面性。但由于數(shù)據(jù)的可獲取性有限，僅研究了部分現(xiàn)象，沒有具體分析肥城市土地利用類型變化的驅(qū)動因子，并且直接利用遙感影像數(shù)據(jù)進行監(jiān)督分類，土地利用分類精度可能會存在一定的誤差。因此，在未來的研究中，研究者應(yīng)完善對肥城市土地利用類型驅(qū)動因子的分析，提高分類精度，以得到更精確的研究結(jié)果。

參考文獻：

［1］何英彬，陳佑啟.土地利用/覆蓋變化研究綜述［J］.中國農(nóng)業(yè)資源與區(qū)劃，2004（2）：61-65.

［2］陳佑啟，楊鵬.國際上土地利用/土地覆蓋變化研究的新進展［J］.經(jīng)濟地理，2001（1）：95-100.

［3］陳小蘭，胡娉婷.基于GIS的東莞市土地利用時空變化規(guī)律和特征分析［J］.廣西城鎮(zhèn)建設(shè)，2021（7）：19-21.

［4］陳伊多，楊慶媛.西藏自治區(qū)土地利用/覆被變化時空演變特征及驅(qū)動因素［J］.水土保持學報，2022，36（5）：173-180.

［5］馮應(yīng)斌，楊溪源.基于土地利用變化的貴州省耕地顯性邊際化時空演變及其影響因素［J］.水土保持通報，2023，43（6）：263-271.

［6］廖美玉，方秀琴，蔣心遠，等.黃河流域近40年土地利用/覆被時空變化特征［J/OL］.水土保持學報，1-14［2024-04-19］.https：//doi.org/10.13870/j.cnki.stbcxb.2024.02.030.

［7］張晨基，張永福.基于DEM的新疆土地利用變化時空格局［J］.中南農(nóng)業(yè)科技，2024，45（1）：116-121，126.

［8］安彬.1985—2020年漢江生態(tài)經(jīng)濟帶土地利用/覆蓋演變特征分析［J］.安康學院學報，2023，35（6）：16-23.

［9］周攀清.魯西南地區(qū)村落喪葬習俗探析：以肥城市村落調(diào)查為例［J］.學理論，2014（15）：164-165.

［10］郭健斌，張英，張志偉，等.基于地理探測器的藏東南高山峽谷區(qū)土地利用變化及其驅(qū)動機制：以西藏林芝市為例［J］.中國農(nóng)業(yè)大學學報，2023，28（4）：210-226.

［11］王秀蘭，包玉海.土地利用動態(tài)變化研究方法探討［J］.地理科學進展，1999（1）：83-89.

［12］陳伊多，楊慶媛.西藏自治區(qū)土地利用/覆被變化時空演變特征及驅(qū)動因素［J］.水土保持學報，2022，36（5）：173-180.

［13］任斐鵬，江源，熊興，等.東江流域近20年土地利用變化的時空差異特征分析［J］.資源科學，2011，33（1）：143-152.

［14］張瑞，杜國明，張樹文.1986—2020年東北典型黑土區(qū)耕地資源時空變化及其驅(qū)動因素［J］.資源科學，2023，45（5）：939-950.

［15］薛晉聰，陶志梅.天津市各區(qū)土地利用程度差異及動態(tài)變化分析［J］.合作經(jīng)濟與科技，2019（8）：8-11.

作者簡介：石舒心（2003—），女，本科生，研究方向：GIS。