基于GIS與景觀指數的黃驊市景觀格局分析

王永 劉春玲 賀鵬 姚翔龍

摘要通過對黃驊市柵格影像的景觀指數計算，獲得10個斑塊類型的景觀指數和9個景觀水平的指數。對黃驊市土地利用類型進行重分類，通過對景觀指數的分析，得出黃驊市農業景觀、荒漠景觀等景觀分布格局。結果表明，黃驊市景觀整體聚集程度較高、多樣性顯著、均勻度較高、破碎化程度較高，形成以農業景觀為景觀基質，荒漠、陸地景觀散碎分布，工業、水域景觀團塊式分布，居住、交通運輸景觀規律分布的景觀格局。

關鍵詞景觀指數；景觀格局；斑塊；黃驊市

中圖分類號TU984.18；TP79文獻標識碼

A文章編號0517-6611（2017）16-0063-05

Landscape Pattern Analysis Based on GIS and Landscape Index in Huanghua City

WANG Yong1， LIU Chunling2， HE Peng2 et al（1.Faculty of Engineering，China University of Geosciences， Beijing 100083； 2. China Land and Resources Aerial Geophysical Remote Sensing Center，Beijing 100083）

AbstractThrough the calculation by the Huanghua City grid image for the landscape index，the landscape indices of 10 patch types and 9 landscape levels were obtained.Land use types in Huanghua City were reclassified. Through the analysis of landscape index，the landscape distribution pattern of agricultural landscape and desert landscape in Huanghua City was obtained. The results showed that Huanghua City landscape pattern have high degree of aggregation， significant diversity， higher degree of fragmentation and high uniformity. The formation of the agricultural landscape as the landscape stroma，desert and land landscape scattered distribution，industrial and water landscape mass distribution， living and transportation regular distribution landscape pattern.

Key words

Landscape index；Landscape pattern；Plaque；Huanghua City

景觀生態學的研究對象涵蓋了從荒野到城市景觀，對景觀要素或景觀格局現狀分析是景觀生態學的重點研究內容之一[1]。而對景觀內各斑塊類型在景觀中分布規律的研究是景觀格局現狀分析的基礎。目前，以土地利用類型所反映的景觀格局幾何特征，結合數值統計的分析方法在景觀格局分析方面較為成熟。在區域景觀格局分析研究領域，專家學者進行了大量研究。曹瑞娜等[2]通過棲霞市各流域信息提取，計算獲得棲霞市各流域景觀指數，結合DEM模型、自然因素以及地區發展模式對棲霞市進行景觀格局分析。劉飛等[3]通過不同斑塊類型指數的選取計算，分析了貴州習水縣景觀空間特征。寧雅楠等[4]通過景觀指數法結合主成分分析法對青龍滿族自治縣進行景觀分區研究。另外，用景觀指數研究景觀格局對指數的選取依賴性較大，部分學者利用統計回歸、因子分析等方法，建立景觀指數之間的相關性、敏感性的關系。何鵬等[5]通過相關分析方法，解決了指數繁多冗余的問題。彭建等[6]通過研究深圳寶安區土地利用分類，探討景觀格局指數隨土地利用分類系統變化的基本規律。由于人類對自然資源不斷開發和謀求土地資源的可持續利用，使得景觀格局隨空間和時間的轉變而不斷改變[1]。因此，通過對景觀格局的分析有助于把握研究區宏觀的生態環境狀況、景觀格局發展趨勢及科學合理地規劃人們的生活環境，實現土地資源和環境的可持續發展。筆者基于GIS與景觀指數，對黃驊市土地利用類型進行了重分類，研究該市的農業景觀、荒漠景觀的分布格局，以期為今后黃驊市土地資源管理及土地整治提供科學依據。

1研究區概況與研究方法

1.1研究區概況

黃驊市地處117°05′～117°40′ E，38°09′～38°39′ N，位于河北省東南部，地處環渤海經濟圈中部位置和環京津樞紐地帶，有河北省東出西聯的出海口和橋頭堡之稱。全市總面積2 411.6 km2，截至2010年總人口為44.35萬。主要為平原和海岸地貌，其中海岸平坦寬闊，屬于淤泥質海岸，屬暖溫帶半濕潤大陸性季風氣候。轄區土壤類型可歸為潮土、鹽土和沼澤土類3類，其中鹽土遍布全市。全市交通發達，自然資源豐富，其中鹽業資源較為著名，是全國的海鹽生產基地之一。

1.2數據來源及處理

研究區基礎數據來源于2016年GF1遙感數據影像，以2009年第二次全國土地調查分類表為依據，基于ArcGis10.3進行目視人機交互解譯，最終獲取解譯精度為1∶5萬，共29類土地利用類型的gdb數據庫格式的文件。然后結合當地景觀特色及發展模式，重新對29類土地利用類型進行歸并分類，分類形成8種一級景觀類型（圖1）。

最終利用8種一級景觀類型進行景觀指數的計算。計算步驟：①將8種景觀類型放在同一數據庫中，利用ArcGis10的功能轉化為20×20大小的geo.tiff格式的柵格圖像。②將該圖像導入到 Fragstats 4.2軟件中，設置相關權重文件并選擇參數。③計算。以斑塊類型指數（class類）和景觀水平指數（landscape類）為主，最終以這2類景觀指數對黃驊市景觀格局進行分析。研究區景觀類型分布如圖2所示。

1.3景觀指數軟件Fragstats 4.2經過多年的發展和改良，界面操作簡單，并能與ArcGis10軟件進行良好的互補，其包含的景觀指數多達上百種，能從各方面反映區域景觀格局分布特征。同時景觀指數之間存在明顯的相關性、功能的重疊性甚至相互矛盾的問題。該研究通過遙感解譯了解黃驊市景觀特征，結合前人在指數選取方面的方法以及全面了解所選指標所有生態意義、公式的前提下[6-8]，選取了以下斑塊類型指數和景觀水平指數（表1、2）。

2結果與分析

2.1研究區整體結構

斑塊水平指數中CA、PLAND、NP、PD、LPI指數以及景觀水平指數中TA、NP、PD、LPI、PR指數是對斑塊類型和景觀整體水平的基本特征進行描述的指數類型。其余指數均為描述景觀空間異質性、復雜多樣性及連通性方面的特征。通過指數計算，分別獲得表3所示的景觀格局分析數據。

經計算可知，由景觀水平指數TA、NP、PD、LPI可知，研究區的TA為241 160.44 km2，除去沿岸淺海景觀外，總面積約226 263.00 km2。NP為15 773塊，PD為6.54個/100 hm2，LPI占景觀總面積的10.61%。由表3可知，參與計算的8類景觀類型中，農業景觀面積最大，為97 092.12 hm2，占研究區總面積的40.26%，其次為水域景觀，為38 245.08 hm2，占研究區總面積的15.86%。由于水域景觀包括水庫、水工建筑、河流、沼澤等地類，且研究區與海相鄰，水資源豐富，所以面積較其他土地類型較多。此外，荒漠景觀面積達36 776.40 hm2，占總面積的15.25%；工業景觀為23 560.44 hm2，占總面積的9.77%；居住景觀為22 208.44 km2，占總面積的9.21%；沿岸淺海景觀、交通運輸景觀及陸地景觀面積分別為14 897.68、5 511.40、2 868.88 hm2，分別占總面積的6.18%、2.29%、1.19%。農業景觀、水域景觀及荒漠景觀占總面積的71.00%以上，說明這3類景觀是研究區的優勢景觀類型。在景觀基質方面，由于農業景觀占40.28%，大范圍地分布在市中心周邊及內陸地區（圖2），且連通性較好可以認為是研究區的景觀基質。斑塊類型面積分布如圖3所示。

在NP、PD方面，居住景觀和荒漠景觀是研究區最多的，分別是4 366、4 324個和1.81、1.79，說明居住景觀和荒漠景觀的分布特征是部分斑塊散碎，點狀分布于研究區。其次為農業景觀、水域景觀、工業景觀，NP分別為2 658、1 560、1 295，分別為1.10、0.65、0.54；農業景觀、水域景觀是景觀類型中面積最大的2類，但NP、PD并不是最大。交通運輸景觀和陸地景觀的NP、PD較少，分別為833、731個和0.35、0.30（表3），說明黃驊市草地、林地分布面積較少。沿岸淺海景觀是基于海洋資源大力開發的背景下設立的一種景觀類型，該景觀類型沿岸分布且聚集在一起，因此NP、PD均為最少。

2.2景觀面積特征的景觀格局

由上述可知，以農業景觀為景觀基質的黃驊市在景觀類型面積分布上呈多種現狀。農業景觀面積最大占40%以上，但NP、PD不是最多，說明農業景觀在黃驊市具有連通性。而且農業景觀的AREA-MN為36.53 hm2，是除沿岸淺海景觀外最大的一類，說明黃驊市存在著大量農業用地，聚集度程度較高的分布在研究區。由圖2可知，農業用地大量分布在整個研究區，尤其是內陸地區農業景觀連通性更好，所占比例更大。水域景觀主要分布在研究區東北部，其中水庫水面占絕大部分面積。水域景觀中存在著研究區中LPI屬性為水庫水面，為10.61%。水域景觀和工業景觀連通性較好，水域景觀和工業景觀分別為24.52和18.19 hm2，但斑塊數量和斑塊密度較小；說明水域景觀和工業景觀呈團聚式分布，其中海鹽產出為研究區特色資源，以采鹽場為主的工業景觀集中分布在研究區東部且形狀較為規則（圖2）。居住景觀和荒漠景觀的LPI不大，但AREA-MN卻是最多，分別為5.09和8.51 hm2，說明居住景觀和荒漠景觀內部破碎程度高，人為干擾現象嚴重。居住景觀主要分布在研究區中部，在研究區邊緣地區也散布著居住斑塊，但斑塊類型分布均勻有序、形狀規則。荒漠景觀主要是以鹽堿地為主的景觀特征，荒漠景觀主要分布在研究區中部和東南部，形狀不規則，景觀內部支離破碎，是整個研究區整治和改造的重點部分。在其他景觀方面，陸地景觀和交通運輸景觀所占面積較少；陸地景觀中以草地林地為主的綠化帶在整個研究區分布較少，交通運輸景觀中交通用地交錯縱橫，聚集程度較低。沿岸淺海景觀的AREA-MN最大，達2 482.95 hm2，這與集中在海岸帶一側斑塊聚集程度較高有關。

45卷16期王 永等基于GIS與景觀指數的黃驊市景觀格局分析

2.3景觀形狀及邊緣特征的景觀格局

LSI是通過計算某一斑塊類型的形狀與相同面積的圓或正方形之間的偏離程度來測量斑塊形狀復雜程度。因此，可通過LSI指數反映黃驊市某一景觀類型整體的空間展布特征。由表3可知，LSI由大到小依次為交通運輸、荒漠、農業、水域、居住、陸地、工業及沿岸淺海景觀。交通運輸景觀和荒漠景觀的LSI分別為97.95和87.76。從數值上看，交通運輸景觀和荒漠景觀在研究區的空間展布不均衡。交通運輸景觀的屬性是長條形、交錯縱橫的，荒漠景觀是內部散碎地分布在整個研究區，2個景觀類型偏離正方形的程度都很大。農業景觀、水域景觀和居住景觀的LSI分別為70.60、61.47和52.28。一般來說，景觀類型面積越大，則LSI指數的數值越高，研究區農業景觀面積最大，景觀內部破碎程度較高，但內陸地區的連通性較好，因此LSI指數不是最高。水域景觀中由于水庫水面等面積大、聚集程度較高的土地類型的存在也使LSI指數較低。由于居住景觀主要分布在市中心及其周邊地帶，雖然部分斑塊散碎分布在研究區周邊地帶，但所占面積較小，因此LSI指數較小。陸地景觀、工業景觀占景觀面積不大且呈團聚式分布， LSI指數較低，由于沿岸淺海景觀只分布于海岸一側，形狀面積最接近相同類型的正方形，LSI指數最低，為3.62。研究區的LSI指數分布見圖4a。

面積加權的分形維數CFRAC-AM是描述斑塊類型的形狀指標，與LSI側重于面積對比不同的是FRAC從斑塊類型的整體形狀上進行對比。交通運輸景觀和水域景觀的FRAC-AM，由長條分布的交通道路和蔓延曲折形狀的河流水面使FRAC-AM分別為1.43和1.31。其他景觀中，按FRAC-AM越簡單，人為干擾程度就越大的理論，人為干擾程度由強到低分別為工業、居住、沿岸淺海、陸地、荒漠、農業景觀類型，除沿岸淺海景觀外，人工干擾強度次序基本符合實際。沿岸淺海景觀在海岸線一側形狀規則，遙感解譯時人為因素較大，所以FRAC-AM較高。陸地景觀面積最小，但FRAC-AM較高，說明FRAC-AM在描述非優勢種等異質性方面較為突出。

2.4景觀聚集度的景觀格局

AI是基于同一類景觀斑塊間公共邊界長度來計算，當同一類景觀斑塊像元間不存在公共邊界時，該類的聚集程度最低；而當同一景觀類型中所有像元值存在的公共邊界達到最大值時，具有最大的聚集指數。研究區各景觀AI均較高（表3），由高到低依次為沿岸淺海、工業、農業、水域、居住、荒漠、陸地、交通運輸景觀；沿岸淺海景觀由于其天然屬性，AI最高，達99.57%，一方面說明干擾程度較小，另一方面沿岸淺海景觀聚集連通程度遠高于其他景觀類型。由于工業景觀是單一的采鹽場，大部分采鹽場的集中分布是AI較高的原因。也說明黃驊市作為海鹽生產基地之一，大面積海鹽集中生產有利于管理和運作。農業景觀由于大范圍的分布與市中心及內陸地區，是整個研究區的景觀基質，在景觀內部被其他景觀類型填充，因此AI僅次于沿岸淺海景觀和工業景觀，為95.53%。水域景觀、居住景觀、荒漠景觀、陸地景觀的AI分別是93.81%、93.11%、90.94%及86.24%，說明各景觀類型中主體斑塊聚集程度較高，但仍有大量散碎的斑塊類型分布在研究區各處，如居住景觀、荒漠景觀和陸地景觀的散碎性，同時也使這3類景觀的聚集度低于其他聚集程度較高的幾類。交通運輸景觀根據其自身屬性AI為 73.79%，通過參考對比交通運輸景觀的聚集度，同樣可以得到其他幾類景觀類型的聚集度和散碎程度。景觀水平的AI達到88.15%，說明在景觀整體水平具有較高的聚集程度和連貫性。

2.5景觀多樣性及蔓延度的景觀格局

SHDI指數在景觀級別上等于各斑塊類型的面積比乘以其值的自然對數之后的和的負值。當SHDI為0時，表明整個景觀僅由一個斑塊組成，無多樣性；SHEI等于1.00，表明各斑塊類型均勻分布，有最大多樣性；SHDI增大，說明斑塊類型增加或斑塊類型在景觀中呈均衡化趨勢。該研究表明，研究區SHDI為1.70，通過對比其他研究的SHDI，如曹瑞娜等[2]基于多流域計算出的SHDI分布圖及其地類分析以及曹偉等[7]基于宜興市每個城鎮的SHDI計算統計，發現該研究區的SHDI略大于基于多流域計算的最大值（1.60），屬于景觀類型豐富，斑塊均勻程度較為平衡。說明黃驊市土地利用類型豐富，斑塊類型復雜多樣，景觀破碎化程度較高，多樣性顯著以及景觀異質性不明顯。SHEI是指除以給定景觀豐度下的最大可能多樣性。研究區的SHEI數值為0.82，數值較高，表明黃驊市各土地利用類型均勻分布，斑塊破碎化程度較為突出，同時也說明黃驊市不存在優勢度非常明顯的景觀類型。

CONTAG是表征景觀格局重要的指數之一，通過景觀中斑塊的團聚程度和延展趨勢描述景觀格局的蔓延度。由表4可知，研究區的CONTAG為52.02%，表明研究區無某種優勢的景觀類型在整個研究范圍形成較好的連通性，即使農業景觀也只是在研究區內陸形成較好的連通性。研究區中的景觀類型是以多種斑塊類型互相嵌套的密集格局，同時也說明研究區破碎化程度較高。

3結論與討論

（1）黃驊市的景觀格局構成按照比例從大到小依次為農業景觀、水域景觀、荒漠景觀、工業景觀、居住景觀、沿岸淺海景觀、交通運輸景觀及陸地景觀，總面積達2 411.6 km2。農業、水域、荒漠景觀構成黃驊市景觀格局的主體特征，占總面積的71.00%以上。農業景觀在黃驊市內陸地區呈大面積、多斑塊、高聚集度的分布，占研究區總面積的40.28%，成為黃驊市的景觀基質。以水庫水面和坑塘水面為主的水域景觀集中分布在黃驊市東北部，斑塊連通性較好，斑塊破碎程度較低。以鹽堿地為主的荒漠景觀是黃驊市需要大力整治的一類景觀類型，大量存在于黃驊市中南部，面積大，破碎程度高，屬黃驊市重點治理土類。工業景觀人工干擾強烈，呈團塊式分布。居住景觀分布較廣，斑塊數量最多，聚集程度較低，形狀規則。以城市綠化為主的陸地景觀整體面積較少，散碎分布，人為干擾程度強烈。

（2）從景觀格局整體來看，研究區PR為8，具有較高的景觀聚集性和均勻度，景觀種類豐富，蔓延度一般。但荒漠景觀突出，分布較為散亂，影響黃驊市的發展進程。農業景觀內部破碎也較為嚴重，邊緣特征復雜，陸地景觀面積較少，有必要加強城市綠化建設。

（3）黃驊市景觀格局比較清晰，下一步可以通過對比黃驊市10年之內景觀格局的動態演變，提出黃驊市景觀格局的演變方向及其優缺點，為黃驊市更好的土地資源管理以及土地整治提供更好的建議意見。

參考文獻

[1]

陳利頂，孫然好，劉海蓮.城市景觀格局演變的生態環境效應研究進展[J].生態學報，2013，33（4）：1042-1050.

[2] 曹瑞娜，齊偉，李樂，等.基于流域的山區景觀格局分析和分區研究：以山東省棲霞市為例[J].中國生態農業學報，2014，22（7）：859-865.

[3] 劉飛，謝雙喜，謝鵬，等.貴州習水縣土地利用景觀格局特征分析[J].天津農業科學，2015，21（12）：84-88.

[4] 寧雅楠，李貝，楊偉州，等.基于主成分分析法的土地利用景觀分區研究：以青龍滿族自治縣為例[J].中國農業資源與區劃，2016，37（2）：22-28.

[5] 何鵬，張會儒.常用景觀指數的因子分析和篩選方法研究[J].林業科學研究，2009，22（4）：470-474.

[6] 彭建，王仰麟，張源，等.土地利用分類對景觀格局指數的影響[J].地理學報，2006，61（2）：157-168.

[7] 曹偉，周生路，吳紹華，等.基于粗糙集與突變級數法的土地利用景觀分區研究[J].地理科學，2011，31（4）：421-426.

[8] 李秀珍，布仁倉，常禹，等.景觀格局指標對不同景觀格局的反應[J].生態學報，2004，24（1）：123-134.