基于rs和gis的北碚区土地利用变化及景观格局分析.docxVIP

基于rs和gis的北碚区土地利用变化及景观格局分析

地球表层系统最重要的景观特征是土地使用和土壤覆盖。土地利用/覆被变化不仅客观记录了人类改变地球表面特征的空间格局,而且还再现了地球表层景观的时空变化过程;近20年来北碚区随着工业化、城镇化的不断发展及人口的快速增长,人地矛盾、人与空间的不协调状态日益严重,在经济发展的同时,也给生态环境带来了一定的负面影响。

使用常规的土地资源调查方法,获取数据的周期长精度差,加上每年土地利用/覆被状况都在变化,往往使耗资巨大的调查结果难以反应当地土地资源的现状,而应用遥感技术则可以快速获得土地利用/覆被的动态信息。并且遥感、地理信息系统和计算机技术已广泛应用于景观格局的定量化研究。地理信息系统使快速、准确地处理大规模空间数据成为可能,它与遥感技术相结合大大拓宽了人类的视觉范围,为景观水平上的监测与评估提供了有效的手段。

1土地使用的变化对景观结构的影响

1.1东北部附作区

北碚区隶属重庆市,地处重庆核心区的西北面嘉陵江下游。地跨东经106°18′～106°56′,北纬29°39′～30°03′,幅员面积755km2,属亚热带湿润;冬暖、春早、秋迟多绵雨;热量丰富,雨量充沛,年平均温度18.3℃,年降水量1105.4mm;耕作制度以一年二熟为主,适宜种植柑桔果木及蚕桑、茶叶等多种经济作物;境内自然条件复杂,植被层次丰富,种类繁多。

1.2试验数据和方法

本文使用的遥感影像为LandsatTM影像和CBERS-2影像,一帧为1988年8月29日,另一帧为2008年9月28日。两帧图像的时间间隔20年,在月份上相差一个月,能够满足研究的需要。另外还有研究区域的行政区划图和1∶5万的地形图及相关年份的土地利用数据。具体工作流程见图1所示。

研究方法:1)土地利用/覆被的研究方法。利用ERDAS软件,将两个时段的影像均以北碚区1∶50000UTM投影的地形图为基准进行纠正,均方差误差(RMS)控制在0.5个象元内,用多项式变换模型进行几何校正。同时进行矢量图到栅格图的转换,利用ERDAS软件的掩膜技术得到研究区的TM和CBERS影像。接着通过监督分类得到研究区域1988年和2008年的土地利用/覆被类型图。根据最新《土地利用现状分类》规定的土地利用分类系统,结合研究区域的实际情况,将研究区域的土地利用/覆被分为水域及水利设施用地、耕地、住宅用地和工矿用地、林地和其他用地5种用地类型。最后把栅格图转换成矢量图,在ARCMAP中结合实地考察和GoogleEarth提供的高分辨图像,进行纠正和分析。

2)景观格局研究方法。景观格局包括空间格局和时间格局,只要空间格局清楚了,时间格局就不难理解了。景观空间格局主要指大小和形状不一的景观斑块在空间上的排列,它是景观异质性的重要表现,同时又是各种生态过程在不同尺度上作用的结果。景观格局变化包括景观形状分析、景观分维数分析、景观破碎度分析等内容。其中景观形状指数(Landscapeshapeindex)是通过计算某一景观类型的斑块周长与同面积圆形的周长之比,来测定该景观类型的复杂程度,计算公式LSI=PiπA???√LSΙ=ΡiπA。式中LSI为景观形状指数。Pi为某一景观类型斑块周长,A为斑块面积;景观分维数(Landscapefractaldimention)主要揭示斑块及斑块组成的景观的形状和面积大小之间的相互关系。计算公式Fd=2ln(P/k)/ln(A),式中Fd为分维数,k为常数,本文k取值4;景观破碎度(Landscapefragmentationindex)指景观被分割的破碎化程度,它与自然资源的保护密切相关。计算公式为LCI=ni/Ai。式中:LCI为景观破碎度;ni为某一景观的斑块数;Ai为某一景观的总面积。

1.3土地利用的数量变化

1988、2008年的土地利用/覆被影像分类的整体精度分别为88.7%和92.6%。显然,这些数据有合理的精度,能够满足研究的需要,分类结果如封二彩图2所示。

利用ARCMAP的空间分析工具ZonalArea得到研究区域的土地利用转移矩阵,从而得到1988—2008年间每一类别的土地利用变化的面积数据(见表1)及土地利用类型占有率图(见图3)。

从上述图表可看出土地利用的数量一直在变化。其中变化幅度依次是耕地林地居住、工矿用地其他用地水域及水利设施用地,耕地和其他用地减少的同时,与之相反的是林地、居住、工矿用地和水体的增加,自1999年国家实行“退耕还林”、“耕地换园地”等政策以来,北碚区政府鼓励实行生态退耕,坡度在25°以上的山坡耕地,拍卖或者承包给当地的农民并大力发展林果,使得农民种林种园的积极性提高,林地大大提高。居民、工矿用地的增加、耕地和其他用地的减少与近20年中国经