基于RS和GIS的百色市西林縣耕地景觀格局演變及其驅動因子分析

童新華張福文韋燕飛

（1.廣西師范學院地理科學與規劃學院，廣西 南寧 530001；2.廣西師范學院國土資源與測繪學院，廣西 南寧 530001）

基于RS和GIS的百色市西林縣耕地景觀格局演變及其驅動因子分析

童新華1張福文1韋燕飛2

（1.廣西師范學院地理科學與規劃學院，廣西 南寧 530001；2.廣西師范學院國土資源與測繪學院，廣西 南寧 530001）

利用Erdas軟件對獲取的2005年和2012年西林縣TM影像數據進行遙感解譯，并用ARCGIS軟件對圖像進行數字化，柵格格式轉換處理，最終獲取百色市西林縣的耕地格局；選擇相關景觀指數使用 fragstas3.4軟件結合相關的數學分析方法，對西林縣的耕地空間分布特征進行定量分析，研究其變化的驅動因素。結果顯示：西林縣2012年的耕地景觀比2005年更顯離散化分布，耕地景觀破碎度增加，耕地的數量有所下降，而耕地景觀的整體形狀顯得更有規則化。

耕地；景觀指數；格局演變；驅動因子

國家統計局公布的數據顯示，2013年我國糧食總產量達到 60193.5萬噸，人口總數 13.6億人，人均糧食產量為442.60kg，還遠遠低于發達國家的標準。依據2013年我國國土資源公報，截止2012年底我國耕地面積為20.27億畝。從數據上看全國耕地總面積雖然大于國家提出的18億畝耕地，但由于我國是一個多山區的國家，耕地資源分布極不均勻，特別是西部地區耕地數量少且耕地的質量水平總體較低；另外耕地資源利用方式粗放，糧食生產成本較高、產量提升空間小，耕地保護形勢依然嚴峻。

廣西壯族自治區是一個多山地少平原的地區，2012年底廣西壯族自治區總人口為5240萬人，2012年廣西土地利用年度變更調查成果顯示2012年底廣西的耕地面積為6624.59萬畝，人均耕地面積為 0.0843公頃，低于全國的平均水平0.093公頃[1]；耕地資源總量較少，經濟、科技相對比較落后，對耕地資源的利用更是處于比較粗淺的層次，耕地資源少及利用不合理的雙重制約，更加加劇了快速增長的糧食需求與糧食產量不足的矛盾。因此切實保護耕地，充分合理利用耕地資源成為了人們生活的重中之重。

1 研究區概況

西林縣位于廣西壯族自治區最西端，地處桂、滇、黔三省（區）結合部，介于東經 104°29′～105°36′，北緯24°11′～24°44′之間。屬亞熱帶季風氣候，年內最冷月與最熱月平均氣溫差值較小，夏無酷暑、冬無嚴寒。全境以土山為主，占全縣總面積的95%。土壤以山地紅壤為主，其他類型還包括水稻土、黃壤土、石灰巖土、沖積土等四類。土壤質地多屬壤質，土壤結構旱地以團粒結構為主，水田以小塊狀結構為主，土壤耕層疏松。2012年末西林縣土地總面積為299725.99公頃，占廣西壯族自治區國土面積的1.26%，2012年末西林耕地面積為22018.55公頃，占廣西壯族自治區總耕地的0.5%。

圖1 西林縣位置示意圖

2 研究方法

2.1 數據來源

研究數據來源包括：由統計年鑒記錄的2012年西林縣的自然、經濟及人口數據；西林縣2005年和2012年的TM遙感影像數據；西林縣2012年度土地利用年度變更調查數據及部分相關的土地利用規劃數據。

2.2 數據處理及信息獲取

參考全國土地利用分類體系結合研究區域的景觀類型特點[2]，西林縣耕地景觀包括水田、旱地、水澆地、望天田、菜地和其他農用地。在ERDARS9.1軟件的支持下，對2005年和2012年西林縣兩期TM遙感影像進行幾何校正與配準，輻射校正，投影變換，空間增強和裁剪等預處理。在通過野外實地調查獲得研究區域實際經驗的基礎上，將土地利用類型分為耕地和非耕地兩類，并結合2005年和2012年西林縣土地利用現狀圖及野外調查數據對TM影像進行遙感解譯，從而得到西林縣耕地空間分布圖。選用隨機采樣方法在各地取 100個像元，經實地野外調查形成誤差評價表（表1），從表1可以看出2005年和2012年得解譯結果總體精度分別為0.865和0.875，kappa系數分別為0.73和0.75，均滿足最低容許判別精度0.7的要求。

表1 解譯誤差評價表

將解譯結果導入ARCGIS軟件進行數字化，在西林縣土地利用現狀矢量數據基礎上，對數據進行預處理，包括耕地數據的提取，屬性的剔除與屬性的增加，數據的裁剪，拓撲錯誤的改正，同時建立空間數據庫；再利用ARCGIS軟件進行拓撲處理，柵格轉化，空間統計，最后將耕地景觀數據轉化成GRID格式，使用Fragstats3.4軟件進行景觀指數的計算。

2.3 景觀指數選取和計算

景觀格局包括景觀組成單元的類型、數目以及大小不一的景觀斑塊在空間上的配置，是景觀異質性的具體體現[3]。景觀格局特征可以在 3個層次上分析：單個斑塊、由若干單個斑塊類型的整個組成的斑塊類型和包括若干斑塊類型的整個景觀鑲嵌體。相應的景觀格局指數也存在三個水平上不同的指數:斑塊水平指數、類型水平指數和景觀水平指數。斑塊水平指數包括斑塊的面積、形狀、邊界特征等[4]；類型水平有斑塊的平均面積、平均形狀指數、斑塊密度等；景觀水平指數除了各種類型指數外還有各種多樣性指數、聚集度指數、距離指數和生境破碎化指數等。本文選取以下幾個指數進行分析。

（1）平均斑塊最近距離，即斑塊之間距離的大小，可以直接反映景觀要素分布離散或團聚程度，一般來說ENN值大，反映出同類型斑塊間相隔距離遠，分布較離散；反之，說明同類型斑塊間相距近，呈團聚分布[5]。其計算公式為：

式中，hij為景觀中每一個斑塊與其最近鄰體的距離，N’是具有鄰體斑塊的總數。

（2）)斑塊密度（PD）

該指數是斑塊個數與景觀總面積的比值，即每平方千米（100hm2）的斑塊數。值愈大，景觀破碎化程度愈高，空間異質性愈大。它反映了景觀空間結構的復雜性，根據這一指數可以比較耕地景觀破碎化程度。斑塊密度表達式為：

N表示研究區斑塊總數或景觀中某類型景觀的斑塊個數，A表示各研究區域的總面積或景觀的總面積[6]。

（3）最大斑塊指數（LPI）等于某一斑塊類型中的最大斑塊占據整個景觀面積的比例。該指數有助于確定景觀中的優勢類型，其值的大小決定著景觀中的優勢種、內部種的豐度等生態特征；其值的變化可以改變干擾的強度和頻率，反映人類活動的方向和強弱[7]。計算公式為：

其中aij是斑塊的面積，A為包括景觀內部背景在內的總面積。

（4）平均形狀斑塊指數（SHAPE＿MN）

該指數是度量景觀斑塊邊界的規則程度、斑塊形狀復雜性的重要指標。該指數的值越小說明景觀斑塊的形狀越接近最為簡單的方形，值越大說明斑塊形狀越復雜越沒有規則[8]。計算公式為：

式中：m為斑塊類型，p為斑塊的周長，a為斑塊的面積，N為景觀中的斑塊總數。

（5）平均斑塊分維數（FRAC＿MN）是景觀組中單個斑塊的分維數以面積為基準的加權平均值，用來衡量斑塊邊界的復雜程度，值越大說明斑塊邊界的形狀越復雜，它的取值范圍是1≤FRAC＿MN≤2[9]。計算公式為：

其中aij，pij分別是斑塊的面積和周長。m為斑塊類型，n為斑塊的數量。

（6）聚合度指數（AI）是基于同類型斑塊公共邊界長度來計算，當類型斑塊不存在公共邊界時，該類型的聚合度最低，而當類型斑塊的公共邊界達到最大值時，該類型的聚合度最高。可以反映景觀組分的空間配置特征，度量斑塊類型的非隨機性或聚合程度。計算公式為

其中gij是相應景觀類型的相似鄰接斑塊數量[10]。

3 西林縣耕地景觀格局變化結果與分析

利用RS和GIS方法處理西林縣2005年和2012年兩期TM遙感圖像，并結合土地利用現狀數據建立耕地空間數據庫，并Fragstats3.4軟件進行景觀指數的計算，得出西林縣耕地景觀格局各種分布特征。

3.1 耕地景觀單元特征變化

西林縣耕地景觀斑塊數從2005年的17370減少到17294，斑塊總面積減少了 24.3961hm22，年均減少 3.4852 hm22。另外斑塊最小面積減小，斑塊的平均面積增大。斑塊總面積的減少及斑塊最小面積的減小反映了耕地景觀呈退化的趨勢，這一方面與國家實施的退耕還林政策有關[11]；另一方面也受快速的城鎮化及工業發展等人類的經濟建設擠占耕地資源的影響。

表2 耕地景觀單元特征指數

3.2 耕地景觀斑塊數量特征

西林縣 2005年、2012年耕地斑塊面積在 1hm2以下和1-5hm2的斑塊占的比重較大，在 5hm2以下的分別達到了95.64%和96%，而10hm2以上的斑塊數均不到總數的2%，體現了該縣耕地較為破碎零散，不易于管理和有效利用耕地資源。1hm2以下的斑塊比重由71.42%增加到72.32%，1hm2以上的比重均呈下降趨勢，說明斑塊破碎程度進一步增大，反映耕地資源退化與流失。

圖2 年度斑塊面積等級比重圖

3.3 耕地景觀指數變化

耕地景觀指數變化表明：西林縣2012年耕地景觀平均斑塊最近距離（ENN＿MN）大于 2005年，同時斑塊密度（PD）減小，聚合度（AI）降低，說明耕地景觀的破碎度增加，耕地的空間分布日趨離散，集中度越來越來小。2012年最大斑塊指數（LPI）由2005年得1.9578上升到2012年的1.9600，說明面積較大的耕地斑塊數量有所增加，連片的耕地分布變得更加緊密。2012年的耕地平均斑塊形狀指數（SHAPE＿MN）和平均斑塊分維數（FRAC＿MN）均比2005有所降低，表明耕地景觀邊界日益規則化，形狀的復雜程度下降。整體表明，隨著土地利用總體規劃的實施及基本農田劃定工作的實行，西林的耕地景觀整體的復雜度降低變得日趨有規則，較之以前更方便耕作和管理；但隨著快速的經濟社會發展和城鎮化，人地矛盾的進一步加劇，耕地受到的破壞和侵占也進一步加重，往后保護耕地的工作有待進一步的加強。

表3 耕地景觀格局指數

3.4 驅動因子分析

通過2005年和2012年兩個時段的耕地景觀格局對比分析，可以得出西林縣的耕地景觀數量有所下降，空間分布呈離散化趨勢，耕地景觀的內部結構不穩定，存在退化的現象。結合近年來的西林縣統計年鑒和社會經濟實際數據，把影響耕地景觀格局演變的因素分為內驅動因子和外驅動因子。內驅動因子主要包括地形、地貌、區位、水熱、光照等；外驅動因子主要包括工業化、城鎮化、耕地政策等人為的干擾活動。

3.4.1 內驅動因子

西林縣地處云貴高原向廣西丘陵過渡地帶，全境以土山為主，占全縣總面積的95%[12]。地貌類型復雜，地勢起伏不平，谷深坡陡，耕地主要沿著河流兩岸零散分布；另外該地區巖溶地貌發育，水土流失發生和發展的潛在危險性很大，為治理水土流失不得不實行退耕還林政策，其中 1996-2005年，退耕還林面積就達到了 9586.73公頃。可見復雜地貌的類型及該地區生態安全的需求是決定該地區耕地離散化分布及耕地資源量少的主要內驅動因子。

3.4.2 外驅動因子

該地區以山地丘陵為主，適合開發建設的用地較少，同時為治理水土流失需保護森林植被，為了發展和生存不得不占用部分耕地。但隨著社會經濟的快速發展及人口的增加，城鎮化及工業化的同步推進，人們的生產生活用地的需求量也急劇增加，使得耕地保護布局與城鎮化、工業化布局時有重疊，耕地保護與經濟社會發展沖突日益增多，人地矛盾日趨加重。人類生存，生產需求對耕地產生的強烈干擾活動是驅動該區域耕地景觀格局發生變化的主要外在因子。

4 結語

運用RS技術和GIS技術結合景觀格局分析方法，對西林縣耕地景觀格局的變化研究及驅動因子分析表明，該縣2005-2012年間的耕地景觀格局破碎程度有所增加，主要表現為斑塊的最大面積減小，面積較小的斑塊數量增加；景觀結構復雜程度有所降低，主要表現為斑塊邊界日益規則化；景觀數量有所下降，主要表現為2012年的耕地數量比2005年的少。耕地景觀格局的變化一方面受由該地區的地貌以丘陵土山為主，耕地資源總量少，同時為治理水土流失保證生態安全實行退耕還林政策驅動變化；另一方面，經濟社會的發展，人口的增加，“一要吃飯，二要建設”兩者矛盾的激化也是驅動耕地景觀格局變化的重要因子。

[1] 譚莉梅.河北省太行山區域耕地資源空間分布特征研究[J].中國生態農業學報,2012,18(4):872-875.

[2] 于柳媚.右江河谷土地利用與景觀格局變化的空間分異特征研究[J].自然科學學報,2014:81-86.

[3] 陳張麗,吳志峰.基于遙感和 GIS的廣州市天河區水域景觀演變及其驅動因子分析[J].國土與自然資源研究,2012,(6):55-57.

[4] 黃琯,吳娜.基于景觀單元的四川德陽市旌陽區城市綠地系統景觀生態規劃研究[D].雅安:四川農業大學,2013.

[5] 謝軍飛,李延明,郭佳.北京城市公共綠地景觀格局分析[J].西北林學院學報,2011,26(3):194-198.

[6] 王玉學.張雪松.基于景觀生態理論的土地整治分區研究[D].武漢:華中師范大學,2012.

[7] 米仕洪.基于景觀格局的成都市溫江區綠地網格構建[D].雅安:四川農業大學,2012.

[8] 張永民,鐘林生.松嫩平原典型鹽堿化地區土地利用格局變化研究[J].干旱區資源與環境,2008,9(22):124-128.

[9] 肖捷穎.土地利用與景觀格局變化的空間分異特征研究[J].中國生態農業學報,2010,18(2):416-421.

[10] Lichen Liu,Xiaoying Li.Quantitative analysis of urban expansion using RS and GIS,a case study in Lanzhou [J].Journal of Urban Planning and Development,2012.

[11] 馬慧.農村居民點用地空間特征分析[J].國土資源信息技術,2014,(2):27-31.

[12] 廣西壯族自治區統計局.廣西統計年鑒-2012[M].北京:中國統計出版社,2012.

Based on RS and GIS in the Baise Xilin farmland landscape pattern evolution and its driving factors analysis

Using Erdas software to obtain the xilin county in 2005 and 2005 TM remote sensing image data, and interpret the image by ARCGIS software digital, grid format conversion processing, eventually taking baise distribution of farmland landscape xilin；Select the relevant landscape index using fragstas3.4 software combined with related mathematical analysis method, carries on the quantitative analysis on the arable land, the spatial distribution characteristics of xilin, to study the changes of drivers. Results show that the xilin of farmland landscape in 2012 more than 2005 discrete distribution of farmland landscape fragmentation increases, the amount of cultivated land decreased, while the overall shape of the farmland landscape more rules.

arable land;landscape index;pattern of evolution;driving factors

S29

A

1008-1151(2015)02-0052-03

2015-01-10

童新華（1964－），男，廣西百色人，廣西師范學院地理科學與規劃學院副教授，從事國土信息技術研究。