基于TM影像的青海省瑪多縣土地利用及景觀分析

張朋濤+殷青軍

摘要：通過遙感圖像預處理、人工解譯1990年、2000年、2007年和2010年四期TM影像得到土地利用數據，求得土地利用轉移矩陣和各種生態環境景觀指數。通過分析四期土地利用數據和總體景觀指數之間的差異性，指出由自然和人為原因引起的生態環境的變化：瑪多縣的生態環境已經得到較為明顯的改善。

關鍵詞：遙感數據；瑪多縣；土地利用；景觀分析

1 引言

瑪多縣是國內人類生存環境最惡劣的地區之一，生態環境也相對比較脆弱。但其位于“三江源”國家級自然保護區[1]核心腹地，生態地位十分重要，是青藏高原重要的生態屏障，也是黃河中下游地區經濟社會發展的重要生態功能平衡區。但是隨著人類活動影響的增加，當地生態環境破壞也較為嚴重，政府部門也開始生態恢復重建工作，而利用遙感和GIS技術作為監測技術，分析土地利用與覆蓋變化是重要的研究方向[2-3]。

2 研究區域概況

瑪多縣位于東經96°55′～99°20′，北緯33°50′～35°40′，東部與海南藏族自治州興海縣和果洛藏族自治州瑪沁縣相毗連，北依積石山支脈布青山與海西蒙古族藏族自治州都蘭縣相鄰，西靠玉樹藏族自治州曲麻萊縣，西南以巴顏喀拉山為界與玉樹藏族自治州稱多縣接壤，南與四川省石渠縣和果洛州達日縣相連。瑪多縣是青海省平均海拔最高的縣，全縣總面積24477平方公里，平均海拔4500米以上，縣城駐地海拔4276米，自然條件惡劣，高寒缺氧，環境嚴酷，年均氣溫-4℃，年日照時數4442.2小時，年均降水量303.9毫米。

3 研究內容：

3.1數據處理

所用的數據為1990年8月、2000年10月、2007年8月和2010年7月四期的TM遙感影像，經過幾何糾正、投影變換、建立統一坐標系統等預處理之后，采用人工解譯的方法，對四期影像數據進行解譯分類處理。以上處理在遙感軟件ERDAS和ARCGIS下進行。

3.2土地利用分類

土地利用分類參考國家土地分類系統，結合外業調查。地形地貌和土壤類型等資料，人工解譯對四期遙感影像進行土地利用分類。對最終解譯成果做出土地利用專題地圖，見圖1。

圖1 各年份土地利用圖

3.3景觀指數

景觀指數是指能高度濃縮景觀格局信息，反映其結構組成和空間配置某些方面特征的簡單定量指標[4]。景觀指數有兩個層次[5]，即總體指數和要素指數。本次研究根據解譯出的土地利用數據計算了5種總體景觀指數，如多樣性指數、優勢度指數、分維數、均勻度指數、破碎度指數描述總體景觀的要素特征。

4 結果分析

4.1 土地利用變化分析

根據四期土地利用解譯數據求得土地利用轉移矩陣，見表1、表2、表3。從1990年到2000年變化最大的是湖泊，共有35.18km?面積發生變化，其中32.25 km?變為濕地、1.04 km?變為中覆蓋草地、0.27 km?變為低覆蓋草地、1.62 km?變為灘地；但需要注意的是有2.71 km?高覆蓋草地和1.26 km?濕地變為湖泊。從2000年到2007年變化最大的是濕地，有41.66 km?面積發生變化，其中主要變化為34.06 km?變為湖泊，另外有0.43 km?變為高覆蓋草地、7.17 km?變為水庫坑塘；湖泊面積變化也相對較大，有12.65 km?變為濕地、0.03 km?變為高覆蓋草地、0.57 km?變為中覆蓋草地、3.74 km?變為低覆蓋草地、0.33 km?變為書庫坑塘、0.55 km?變為灘地；需要注意的是有5.32 km?的高覆蓋草地變為中覆蓋草地。從2007年到2010年變化最大的是湖泊面積，增加了55.59 km?，變化來源為濕地44.42 km?、高覆蓋草地0.74 km?、中覆蓋草地5.63km?、低覆蓋草地2.44 km?、河流1.00 km?、灘地1.70 km?；另外面積變化較大的是中覆蓋草地和低覆蓋草地，有40.69 km?中覆蓋草地變為高覆蓋草地、3.79 km?低覆蓋草地變為高覆蓋草地。從1990年到2010年總體上變化比較大的動態變化是有46.15 km?中覆蓋草地變為高覆蓋草地、41.26 km?濕地變為湖泊。總的看來，湖泊和濕地之間的轉化較多，高覆蓋草地和中低覆蓋草地之間的轉化也較多，有一部分濕地變為湖泊，也有一部分中低覆蓋草地變為高覆蓋草地；另外一個顯著變化是從2000年到2007年水庫坑塘這一土地利用類型從無到有，其他土地利用類型面積變化不大。研究區域內主要土地利用類型為草地、湖泊和濕地，在四期土地利用數據中，高覆蓋草地面積分別占總面積的33.19%、33.19%、33.17%、33.34%，中覆蓋草地面積分別占總面積的26.27%、26.26%、26.28%、26.09%，低覆蓋草地面積分別占總面積的19.10%、19.13%、19.09%、19.06%，湖泊面積分別占總面積的10.57%、10.70%、10.58%、10.39。

表 1 1990—2000年土地利用轉移矩陣 km?

有林地 灌木林 高覆蓋草地 中覆蓋草地 低覆蓋草地 河流 湖泊 灘地 城鎮居民點 農村居民點 工程用地 沙地 濕地 裸巖石礫地 總計

（km?）

有林地 0.03 0.03

灌木林地 15.78 15.78

高覆蓋草地 8120.51 2.71 8123.22

中覆蓋草地 4.61 6426.54 6431.15

低覆蓋草地 0.33 4675.33 4675.66

河流 264.33 264.33

湖泊 1.04 0.27 1615.12 1.62 32.25 1650.30

灘地 45.74 45.74endprint

城鎮居民點 1.03 1.03

農村居民點 1.14 1.14

工程用地 6.58 6.58

沙地 185.53 185.53

濕地 1.26 2585.88 2587.14

裸巖石礫地 489.67 489.67

總計 0.03 15.78 8125.12 6427.91 4675.60 264.33 1619.09 47.36 1.03 1.14 6.58 185.53 2618.13 489.67 24477.30

表 2 2000—2007年土地利用轉移矩陣 km?

有林地 灌木林地 高覆蓋草地 中覆蓋草地 低覆蓋草地 河流 湖泊 水庫、坑塘 灘地 城鎮居民點 農村居民點 工程用地 沙地 濕地 裸巖石礫地 總計

（km?）

有林地 0.03 0.03

灌木林地 15.78 15.78

高覆蓋草地 8119.44 5.32 0.36 8125.12

中覆蓋草地 6426.30 0.27 1.34 6427.91

低覆蓋草地 4669.11 0.12 6.38 4675.60

河流 262.48 0.36 1.49 264.32

湖泊 0.03 0.57 3.74 1601.22 0.33 0.55 12.65 1619.09

灘地 0.11 0.15 1.65 45.45 47.36

城鎮居民點 1.03 1.03

農村居民點 1.14 1.14

工程用地 6.58 6.58

沙地 185.53 185.53

濕地 0.43 34.06 7.17 2576.46 2618.13

裸巖石礫地 489.67 489.67

總計 0.03 15.78 8119.90 6432.31 4672.85 262.62 1636.39 18.36 45.99 1.03 1.14 6.58 185.53 2589.12 489.67 24477.30

表 3 2007—2010年土地利用轉移矩陣 km?

有林地 灌木林地 高覆蓋草地 中覆蓋草地 低覆蓋草地 河流 湖泊 水庫、坑塘 灘地 城鎮居民點 農村居民點 工程用地 沙地 濕地 裸巖石礫地 總計

（km?）

有林地 0.03 0.03

灌木林地 15.78 15.78

高覆蓋草地 8118.34 0.74 0.83 8119.90

中覆蓋草地 40.69 6385.99 5.63 6432.31

低覆蓋草地 3.79 4666.62 2.44 4672.85

河流 261.62 1.00 262.62

湖泊 1636.03 0.35 1636.38

水庫、坑塘 18.36 18.36

灘地 1.70 44.29 46.00

城鎮居民點 1.03 1.03

農村居民點 1.14 1.14

工程用地 6.58 6.58

沙地 185.53 185.53

濕地 44.42 2544.70 2589.12

裸巖石礫地 489.67 489.67

總計 0.03 15.78 8162.82 6385.99 4666.62 261.62 1691.97 18.36 44.29 1.03 1.14 6.58 185.53 2545.88 489.67 24477.30

4.2景觀指數分析

根據土地利用解譯數據分析了體現總體特征的總體景觀指數。從表4中數據可以看出多樣性指數略微增大，整體來看呈增長趨勢，但是2000年相對于1990年是從1.8463減小到1.8459；分維數總體上稍有減小；優勢度指數除了2000年小于1990年，整體上一直在增大；破碎度指數整體上是減小的趨勢，但2010年由207年的0.3726增加到0.3733；均勻度指數整體上也是逐漸減小的趨勢。

表4 總體景觀指數

多樣性

指數 分維

數 優勢度

指數 破碎度

指數 均勻度

指數

1990 1.8463 1.5472 0.3473 0.3731 0.9105

2000 1.8459 1.5463 0.3469 0.3728 0.9104

2007 1.8471 1.5467 0.3485 0.3726 0.9111

2010 1.8468 1.5456 0.3493 0.3733 0.8969

5 結論

根據土地利用變化分析，可知在所研究的時間段內瑪多縣高覆蓋草地面積是一個先減小再增加的變化過程，而中覆蓋草地和低覆蓋草地面積則是先增加再減小的變化過程；湖泊面積同樣是先減少后增加的變化過程，而濕地面積則是先增加再減少的變化過程。在研究時段內瑪多縣生態環境整體上是在得到逐步的改善，也驗證了2000年成立三江源自然保護區之后所做的大量生態恢復工作是有效的。分維數的減小、優勢度的增加和均勻度的減小說明了瑪多縣的整體景觀類型向著更好的方向發展，破碎度指數的變化也驗證了這一點，但是2010破碎度指數的增加則表明人類活動對整個研究區域的影響在變大。結合總體景觀指數來看，當地居民也在不破壞生態環境的基礎上保證經濟的發展。從整個區域的土地利用動態變化來看，瑪多縣的生態環境已經得到較為明顯的改善。

參考文獻：

[1] 王紅星.三江源自然保護區[J].黑龍江水專學報.2000，27（4）：61.

[2] 趙庚星，李玉環，虛春達.遙感和GIS支持的土地利用動態監測研究：以黃河三角洲墾利縣為例[J].應用生態學報，2000，11（4）：573-576.

[3] 王長耀，駱成鳳，齊述華，等.NDVI-Ts空間全國土地利用覆蓋分類方法研究[J].遙感學報，2005，9（1）：93-99.

[4] 李遠方.廢黃河三角洲的演變[J].地理研究，1991，10（4）：29-39.

[5] 鄔建國.景觀生態學——格局、過程、尺度與等級[M].北京：高等教育出版社，2000.

作者簡介：張朋濤，1988年生，青海師范大學生命與地理科學學院在讀碩士研究生endprint