張玉杰，劉新俠，張安兵，王賀封

(1.河北工程大學 資源學院，河北 邯鄲 056038；2.河北省煤炭資源綜合開發(fā)與利用協(xié)同創(chuàng)新中心，河北 邯鄲 056038)

?

基于TM影像的漳河流域土地利用時空演變分析

張玉杰1，2，劉新俠1，張安兵1，2，王賀封1，2

(1.河北工程大學 資源學院，河北 邯鄲 056038；2.河北省煤炭資源綜合開發(fā)與利用協(xié)同創(chuàng)新中心，河北 邯鄲 056038)

摘要：為獲取和進一步探討漳河流域土地利用時空演變規(guī)律，基于1993、2006及2010年3期Landsat TM 遙感影像數(shù)據(jù)和2011年SPOT5高分影像數(shù)據(jù)，利用ArcGIS空間分析功能與ENVI動態(tài)監(jiān)測功能對流域土地利用時空演變特征進行研究。分析結果顯示：1993-2010年，漳河上游區(qū)域下墊面土地利用發(fā)生了巨大改變；耕地和裸地呈減少趨勢，各減少411.23、341.11 km2；林地、城鎮(zhèn)與建設用地和水域用地分別增加了656.8、72.01、23.53 km2；耕地與裸地主要轉入林地，部分耕地轉入城鎮(zhèn)與建設用地。

關鍵詞：TM影像；土地利用/覆被變化；漳河流域

土地利用/覆被變化可引起許多自然現(xiàn)象和生態(tài)環(huán)境的變化。它反應了不同時期人類出于各種目的對土地利用方式的改變。這種變化既包括土地資源數(shù)量、質(zhì)量隨時間的變化，也包括土地利用空間格局的變化及土地利用類型組合方式的變化[1-2]。近年來對土地利用變化的研究成果較多，如朱運海等[3]探討了土地利用/覆被變化(LUCC)自動檢測方法，得出綜合應用上述方法能夠高效、準確的自動監(jiān)測出LUCC信息；HAO Hui mei等[4]對中國北部農(nóng)牧交錯帶地區(qū)的土地利用/覆被變化及引起的生態(tài)環(huán)境效應進行研究；黃寶榮等[5]分析了北京市城鄉(xiāng)結合部11年的土地利用/覆被變化情況，總結了造成該地區(qū)景觀格局演變的主要驅(qū)動因子；陳松林[6]、鄒亞榮[7]等研究了土地利用類型與土壤侵蝕的分布相關關系；陳龍泉[8]等分析了基于Markov-CA的土地利用/土地覆蓋變化動態(tài)模型并進行了驗證；李忠鋒[9]等利用RS與GIS分析榆林地區(qū)土地利用變化情況，得出該地區(qū)土地利用類型發(fā)生了顯著變化；Xiaohe Gu 等[10]基于遙感圖像對北京地區(qū)近30年土地利用/覆被變化進行動態(tài)監(jiān)測并結合驅(qū)動因子進行了針對性分析。漳河上游流域地區(qū)隨著社會經(jīng)濟的快速發(fā)展，人類活動不斷加劇，以及自然降雨的減少，土地利用/土地覆被變化較大，土地類型間轉換頻繁。為探索分析該區(qū)域多年間地表覆被變化狀況，本文利用RS、GIS技術對該流域土地利用類型變化強度進行了分析，得出該區(qū)域土地類型基本變化趨勢。

1 研究區(qū)概況

漳河位于海河流域西南部，屬漳衛(wèi)河水系一級支流，由上游的清漳河和濁漳河匯合而成?？鐣x、冀、豫三省四個地級市的22個縣(市、區(qū))，流域總面積18 284 km2。漳河上游區(qū)域以山區(qū)地貌為主，東部為太行山脈，西部為太岳山脈，兩山之間是長治盆地，面積1 169 km2，是流域內(nèi)主要農(nóng)耕區(qū)。漳河流域四季分明，屬溫帶大陸性氣候。

2 數(shù)據(jù)與處理

2.1數(shù)據(jù)獲取

通過USGS(united states Geological Survey)、中國地理空間數(shù)據(jù)云(Geospatial Data Cloud)獲得美國陸地衛(wèi)星(Landsat-5)拍攝的1993、2006、2010年的Landsat TM5衛(wèi)星圖像(見表1)，該研究采用空間分辨率為30 m的1-5、7等六個波段。選取圖像時，為確保數(shù)據(jù)分析效果，采用了云量低于5%，且直接判讀標志(色調(diào)、色彩、紋理、形狀等)較清晰，植被生長良好、植被生長狀況近似月份的圖像。需要指出的是1993年植被長勢較好的7、8、9月份因圖像云量過大，不能滿足信息提取的要求，本實驗選取了相近的6、10兩個月份的圖像予以代替。輔助驗證數(shù)據(jù)采用由法國發(fā)射的spot-5衛(wèi)星拍攝的2011年5 m分辨率遙感圖像。

表1 獲取LANDSAT-5數(shù)據(jù)

2.2 數(shù)據(jù)處理方法

土地利用動態(tài)度可定量描述區(qū)域土地利用變化的速度，它對描述土地利用變化的區(qū)域差異和預測未來土地利用變化趨勢具有較好的應用。單一土地利用類型動態(tài)度和綜合土地動態(tài)度都是土地利用資源數(shù)量變化模型，本文采用單一土地利用類型動態(tài)度模型來進行研究區(qū)土地利用變化研究。單一土地利用類型動態(tài)度反映的是某一研究區(qū)域一定時期內(nèi)某種土地利用類型的數(shù)量變化情況，其計算公式如下：

式中，K為研究時段內(nèi)某一土地利用類型動態(tài)度；Ua、Ub分別為研究初期及研究末期某一土地利用類型的數(shù)量；T為研究時段長，當T的時段設定為年時， K的值就是該研究區(qū)內(nèi)某種土地利用類型的年變化率，當K為正值時，表示該類土地的年增加率，K為負值時表示該類土地的年減少率；K的絕對值越大說明該類土地年變化率越大。轉移矩陣可全面而又具體地分析區(qū)域土地利用變化的結構特征與各用地類型變化的方向。該方法來源于系統(tǒng)分析中對系統(tǒng)狀態(tài)與狀態(tài)轉移的定量描述，為國際、國內(nèi)所常用[11]。ENVI軟件具備該項分析功能模塊，本項目采用該模塊進行應用分析，得出了該區(qū)域土地利用類型間的轉換關系。

3 漳河上游地區(qū)土地覆蓋變化特征

3.1 土地覆蓋類型圖

利用ENVI軟件進行數(shù)據(jù)預處理，以現(xiàn)場采集和SPOT5目視解譯結果為訓練和驗證樣本，針對該研究區(qū)域選取了600個樣本，并按照2:1的比例進行分配，其中訓練樣本、驗證樣本分別為400、200個。采用監(jiān)督分類方法中的最大似然法(Maximum Likelihood method)進行圖像分類，通過該分類方法分別提取出1993、2006、2010年各類土地覆蓋類型的空間分布信息，如圖1。

表2 漳河上游地區(qū)土地覆蓋類型面積狀況

3.2 土地覆蓋空間格局與動態(tài)變化

(1) 土地覆蓋類型統(tǒng)計分析

由表2可得該地區(qū)土地覆蓋類型的具體變化情況，林地、耕地、裸地是漳河上游流域主要土地利用類型，城鎮(zhèn)與建設用地、水域占地面積較小，兩者占地面積不足2%。1993-2010年，林地、城鎮(zhèn)與建設用地及水域占地面積均呈增加趨勢，但耕地、裸地情況相反；1993-2006年，耕地減少面積最大，減少了356.89 km2，占總面積的比例下降2%；林地、城鎮(zhèn)與建設用地、水域3種類別占地面積均增加，其中林地面積增加量最大，增加了515.25 km2，城鎮(zhèn)與建設用地和水域增加面積分別為39.38、8.32 km2；2006-2010年之間，裸地面積減少量最大，減少135.05 km2；耕地面積減少幅度相對于前14年有所減緩，耕地面積減少54.34 km2。林地、城鎮(zhèn)與建設用地和水域仍保持增加趨勢，分別增加141.58、32.6、15.21 km2。

(2)土地覆蓋類型變化分析

由表3和圖1可知，1993-2006年間漳河上游流域土地覆蓋類型變化情況：林地、城鎮(zhèn)與建設用地、水域占地面積呈增加趨勢，其中林地面積增加515.25 km2，增加面積占林地現(xiàn)狀面積的4.93%，林地面積變化速度36.80 km2/y；城鎮(zhèn)與建設用地面積增加39.38 km2，增加面積占城鎮(zhèn)與建設用地現(xiàn)狀面積的14.52%；年變化速度為2.81 km2/y；水域面積增加8.32 km2，增加面積占水域現(xiàn)狀面積的14.52%，年變化速度達0.594 km2/y；耕地、裸地占地面積呈減少趨勢，其中耕地面積減少356.89 km2，減少面積占耕地現(xiàn)狀面積的8.1%，年變化速度為-25.49 km2/y；裸地面積減少206.06 km2，減少面積占裸地現(xiàn)狀面積的9.03%，年變化速度為-14.72 km2/y。

2006～2010年間，漳河上游流域土地利用變化情況：該時間段內(nèi)各土地類型變化趨勢未變，但變化幅度發(fā)生一定改變：林地面積增加141.58 km2，增加面積占林地現(xiàn)狀面積的1.34%，林地面積變化速度達28.316 km2/y；城鎮(zhèn)與建設用地面積增加32.6 km2，增加面積占城鎮(zhèn)與建設用地現(xiàn)狀面積的10.73%，年變化速度為6.52 km2/y；水域面積增加15.21 km2，增加面積占水域現(xiàn)狀面積的8.86%，年變化速度達3.042 km2/y；耕地面積減少54.34 km2，減少面積占耕地現(xiàn)狀面積的1.25%，年變化速度為-10.868 km2/y；裸地面積減少135.05 km2，減少面積占裸地現(xiàn)狀面積的1.34%，年變化速度為-27.01 km2/y。

表3 1993、2006、2010年漳河上游流域土地覆蓋類型轉移速度

由本文2.2中公式求得漳河上游流域單一土地利用類型年變化率，見圖2。由該圖可知：1993-2006年，城鎮(zhèn)與建設用地的變化率最大，城鎮(zhèn)與建設用地的年增加率為1.213 %，其次為裸地，該地類年減少率為-0.592%；2006-2010年間，耕地、林地年變化率較低，其中耕地年變化率最小為-0.246%；林地年增加率為0.271%；水域年變化率呈上升趨勢，由0.401% 變?yōu)?.944 %；城鎮(zhèn)與建設用地的年變化率仍為最大，并且有上增趨勢，為2.404%；裸地的年減少率由1993-2006年的-0.592%變?yōu)?006-2010年間的-1.184%；表明漳河上游流域土地類型變化幅度較大，城鎮(zhèn)與建設用地尤為突出。

3.3 土地覆蓋類型的時空變化

單純從面積的變化來看，很難揭示土地利用的內(nèi)部結構變化，為更好地揭示各種土地利用類型的內(nèi)部變化特征，采用土地利用轉移矩陣定量說明土地利用類型之間的相互轉化狀況，分析區(qū)域土地利用變化特征[12]。

利用ENVI的動態(tài)監(jiān)測功能中的分類比較法工具—Change Detection Statistics和流程化動態(tài)監(jiān)測工具—Thematic Change Workflow ，對1993、2010年兩期土地覆蓋分類結果進行動態(tài)監(jiān)測：獲得該地區(qū)土地利用的變化類型、面積、百分比等變化矩陣信息(表4)，該矩陣可以反映出漳河上游地區(qū)不同時期土地覆蓋類型的轉變過程；并得到1993-2010年間，漳河上游地區(qū)土地覆蓋類型空間轉換情況(圖3)。由表4和圖3得出，該流域土地覆蓋的時空變化特點如下：

(1)總面積情況：1993年，漳河上游地區(qū)五種地型的占地面積分別為水域148.2 km2，城鎮(zhèn)與建設用地231.83 km2，森林占地面積為9 928.07 km2，耕地4 768.8 km2，裸地 2 488.13 km2；2010年，各地類占地面積分別為，水域171.73 km2，城鎮(zhèn)與建設用地303.81 km2，森林占地面積為10 584.9 km2，耕地4 357.57 km2，裸地 2 147.02 km2。

(2)不同時期同種地型的差別：從1993～2010年，水域面積增加23.53 km2，城鎮(zhèn)與建設用地增加71.98 km2，森林占地面積變化量最大，相比2010年，增加656.83 km2；耕地占地面積和裸地面積均呈減少趨勢，其中裸地減少量最大，減少面積為411.23 km2，耕地面積減少341.11 km2。

(3)土地類別轉移情況：1993-2010年間，耕地面積轉移量最大，面積為2 114.19 km2，其中耕地面積主要轉變?yōu)榱值?，轉移面積量為1 537.34 km2，其次為裸地400.13 km2，少量轉變?yōu)槌擎?zhèn)與建設用地和水域，分別為128.57 、48.15 km2；水域面積轉出75.19 km2，主要轉變?yōu)榱值睾透?，面積分別為35.87、28.56 km2，少量轉變?yōu)槌擎?zhèn)與建設用地和裸地；林地面積轉出2 106.81 km2，其面積轉移量僅次于耕地，林地主要轉變?yōu)槁愕兀娣e為1 183.92 km2，868.09 km2的林地轉入耕地，少量林地面積轉入水域和城鎮(zhèn)與建設用地；139.81 km2的城鎮(zhèn)與建設用地轉入其他土地類別，主要轉變?yōu)楦?，面積為83.88 km2，另外31.56 km2轉入林地，少量轉入水域和裸地；裸地轉出1 938.72 km2，主要轉入林地和耕地，面積分別為1 166.18、715.12 km2，少量轉移到水域11.25 km2、城鎮(zhèn)與建設用地46.17 km2。

表 4 1993～2010年漳河上游流域土地利用類型轉移矩陣(km2)

4 結論

1) 耕地、林地為漳河上游流域土地利用的主要類型。1993-2010年，林地、水域、城鎮(zhèn)與建設用地均呈遞增趨勢，但耕地和裸地呈遞減趨勢。

2) 近18年，漳河上游地區(qū)不同地類轉移頻繁且復雜，土地利用變化強度較大，其中耕地面積大幅度降低，且主要轉化為林地，少部分轉變?yōu)樗w；同時1 695.65 km2其他土地類別轉入耕地，主要來自林地與裸地。表明該時期漳河上游流域耕地轉化較活躍，生態(tài)環(huán)境得到一定改善，土地利用強度增大，實施墾荒種田。

3)1993-2010年，漳河上游區(qū)域林地面積呈增加趨勢，耕地與裸地為主要補充來源，說明該地區(qū)實施退耕還林政策，大力倡導植樹造林、發(fā)展林業(yè)；同時城鎮(zhèn)與建設用地、水域用地呈不同程度增加，說明人類活動加劇，城市化進程加快，人口壓力增大，對該地區(qū)的生態(tài)環(huán)境產(chǎn)生一定壓力；大量水渠、水庫的修建改變該地區(qū)的土

地利用類型，對生物及周圍環(huán)境產(chǎn)生一定影響。在未來建設中可繼續(xù)鼓勵正確的植樹造林、植草，加強該地區(qū)蓄水能力，減緩地下水開采；減緩城市化進程，保護道路及水渠、水庫附近的植被。

參考文獻：

[1]王秀蘭，包玉海.土地利用動態(tài)變化研究方法探討[J].地理科學進展，1999，18(1)：81-87.

[2]傅伯杰，陳利頂，馬克明.黃土丘陵區(qū)小流域土地利用變化對生態(tài)環(huán)境的影響[J].地理學報，1999，54(3)：241-246.

[3]朱運海，張百平，曹銀璇，等.土地利用/覆被變化遙感檢測方法與應用分析[J].地球信息科學，2007(3)：116-122.

[4]HAO HUI MEI.Land use/land cover change(LUCC) and eco-environment response to LUCC in Farming-Pastoral Zone,China[J].Agricultural Sciences in China，2009(1)：91-97.

[5]黃寶榮，張慧智，王學志.城市擴張對北京市城鄉(xiāng)結合部自然和農(nóng)業(yè)景觀的影響——以昌平區(qū)三鎮(zhèn)為例[J].生態(tài)學報，2014(22)：6756-6766.

[6]陳松林.基于GIS的土壤侵蝕與土地利用關系研究[J].福建師范大學學報：自然科版，2000(1)：106-109.

[7]鄒亞榮，張增祥，周全斌，等.基于GIS的土壤侵蝕與土地利用關系分析[J].水土保持研究，2002(1)：67-69.

[8]陳龍泉，鄭海金.基于Markov-CA的土地利用/土地覆蓋變化動態(tài)模型研究[J].測繪信息與工程，2004(1)：36-38.

[9]李忠鋒，王一謀，馮毓蓀，等.基于RS與GIS的榆林地區(qū)土地利用變化分析[J].水土保持學報，2003(2)：97-99.

[10]GU XIAOHE.Dynamic monitoring and driving power analysis of LUCC based on remote sensing in Beijing in recent thirty years[J].SPIE，2013，89(21)：1-9

[11]史培軍.土地利用順蓋變化與生態(tài)安全響應機制度[M].北京：科學出版社，2004.

[12]趙東波，梁 偉，楊勤科，等.陜北黃土丘陵區(qū)近30年來土地利用動態(tài)變化分析[J].水土保持通報，2008，28(2)：22-28.

(責任編輯王利君)

Analysis on the temporal and spatial process of dynamic change of land use and land cover in upstream of zhanghe river basin using Landsat TM

ZHANG Yu-jie1,2, LIU Xin-xia1, ZHANG An-bing1,2, WANG He-feng1,2

(1. College of Resources, Hebei University of Engineering,Hebei Handan 056038,China;2. Collaborative Innovation Center of the Comprehensive Development and Utilization of Coal Resource, Hebei Handan 056038,China)

Abstract：Land use/land cover change in upstream of Zhanghe river basin was affected by natural and human factors. During the past decades, reduced natural rainfall and frequent human activities has an influence on this region, so this research is important. Through analyzing the spatio-temporal characteristics of dynamic change of land use /land cover change in basin, Land use/land cover change for many years will be extracted in zhanghe river basin ,using the method of ENVI change detection and ArcGIS Spatial analyst, based on 1993,2006,2010 images derived from Thematic Mapper (TM) imagery .Moreover,Spot-5 image derived from 2011,will be used for precision evaluation. These results provide an important statistic reference for land use planning and sustainable development in this region. Finally, results indicate that forest, buildings and water was respectively increased by 656.8km2、72.01km2、23.53km2,rather cultivated and barren was decreased by 411.23 km2、341.11 km2,from 1993 to 2010.The input forest area was mainly originated from farmland and barren. Part of the farmland was transformed into urban land.

Key words：landsat TM; land use and land cover; upstream of Zhanghe

中圖分類號：X171

文獻標識碼：A

文章編號：1673-9469(2016)01-0076-05

doi:10.3969/j.issn.1673-9469.2016.01.017

作者簡介：張玉杰(1988-)，女，山東濰坊人，碩士，研究方向為遙感應用。

基金項目：河北省自然科學基金資助項目(D2015402134)；河北省教育廳自然科學研究項目(QN2014184，YQ2013012)

收稿日期：2015-09-26