土地利用与植被覆盖变化研究-以西安市灞桥区为例

《土地利用与植被覆盖变化研究》课程考核

土地利用与植被覆盖变化研究

——以西安市灞桥区为例

学生姓名：JZ-W

学生学号：251020313

土地利用与植被覆盖变化研究

——以西安市灞桥区为例

学生姓名：JZ-W 学生学号：251020313 摘要本文以西安市灞桥区2000年和2009年两期SOPT影像数据基础，利用ENVI软件处理数据，进行研究区域的土地类型划分，进而通过计算土地利用动态度、土地分类指数、土地利用类型转移矩阵、土地利用程度综合指数等量，分析研究区土地利用及覆被变化基本情况，并对变化结果进行解释，对未来土地利用及覆被变化进行合理预测。

关键词灞桥区，SPOT，土地利用，覆被变化

1 绪论

1.1 研究目的意义

此次研究目的在于通过实践练习，处理遥感影像数据，理解LUCC研究一般过程，理解LUCC内涵，熟练ENVI等相关专业软件的使用，提升科研能力。

2 研究内容和技术路线

2.1 研究内容

利用西安市灞桥区2000年和2009年两期分辨率为10米的SPOT影像数据，进行土地利用分类；计算并分析土地利用动态度、土地分类指数、土地利用类型转移矩阵、土地利用程度综合指数的变化特点与变化量；分析该区土地利用及覆被变化规律；对该区土地利用及覆被变化进行解释；利用两期SPOT 影像数据，对未来10年的区域土地利用类型分布进行马尔可夫预测。

2.2 技术路线

对影像进行监督分类，并对分类结果进行检验，结果满意后再计算衡量土地利用及覆被变化的各项指标，分析各个结果数据得出有关信息，对计算结果进行解释，对未来土地利用类型进行预测。

3 研究区域概况与数据来源

3.1 研究区域概况

灞桥区地处陕西关中盆地中部，西安城东部，系西安市辖区之一。地理坐标为东经108°59′～109°16′，北纬34°10′～34°27′。东与临潼区、蓝田县接壤，西与雁塔区、新城区、未央区相连，南与长安区为邻，北以渭河与高陵县相望。距市中心5公里，全区总面积332平方公里，辖9个街道办事处，33个社区，226个行政村，人口50余万，因境内遗存始建于隋代的古灞桥而得名。

灞桥区以渭河冲积平原为主，具有山、坡、川、滩、塬的多样性地貌特征，包括三个类型区，其中北部为渭河冲积平原区，东部为低山丘陵区，东南部为台塬区。白鹿塬区域地质构造属于渭河断陷形成的阶梯式下降复式地堑，阶地台塬沟堑内以砂孵石、砾石及土状堆积物为主，并有丰厚的新生代沉积。台塬上区地形平缓，地势北高南低，以农业用地为主，地形由南向北依次递减至沟沿，沟顶海拔600m～750m，沟底海拔160m～600m之间，长度约3公里，坡度约30°左右，地面切割成的沟谷长度约为1000～2000m/km2。土壤类型以黄垩土、白垩土、褐色土以及少量的红色土，有机质丰富，保墒性能好，但是透气性差。

灞桥区属于暖温带半湿润大陆性季风区，光、热、水、气、土等自然条件优越。年平均降雨量635mm，日照时间为2026～2719小时，年平均气温为12～13.3℃，极端最高气温为41.7℃，绝对最低气温为－20.6℃。年大于10℃的有效积温为3650～4325℃，年无霜期202～208天。年初霜日期为11月1日以后，终霜期为4月1日以后。风向频率（静风频率）为30%左右，平均风速为2m/S。年平均相对湿度为70%左右。全年适于旅游的天数为290天，其中最佳的旅游时节为3月中旬至6月底、8月初至11月中旬，约210天。

3.2 数据来源

本次研究所用的是2000年和2009年两期遥感数据，该数据由指导老师提供，该数据为10米分辨率的SPOT影像数据，其投影为UTM，大地基准为WGS-84。出于本次研究的目的意义不是对目标区域的具体情况得出可靠研究结果，而是对研究一般过程的把握，所以没有对影像进行投影变换，几何校正等预处理工作。

4 研究过程

4.1 土地利用分类

土地利用分类是进行LUCC研究的基础性工作，本次研究选择的分类方法是监督分类，在分类之前要进行训练样本选择。根据研究区面积及影像特征，笔者把灞桥区土地利用分成四类，分别是：建设用地，耕地，水体，林地。在选择训练区样本时，四类样本数量选择在20个左右。选择完样本之后，利用最大似然法进行分类，并利用混淆矩阵来检验分类精度，由于影像质量较差等原因，分类精度不是很好。2000年及2009年分类结果如下图所示。

4.2 土地利用动态度计算

本研究仅涉及单一类型的土地利用动态度。土地利用动态度分析主要采用土地利用动态度研究土地利用变化模式。土地利用类型动态度指是某研究区一定时间范围内某种土地利用类型的数量变化情况，其表达式为：

K=(U b -U α)/U α×1/T×100% （1）

式中：K 为研究时段内某一土地利用类型动态度；U a 、U b 分别为研究期初及研究期末某一种土地利用类型的数量；T 为研究时段长，当T 的时段设定为年时，K 的值就是该研究区某种土地利用类型年变化率。

应用土地利用动态度分析土地利用类型的动态变化，可以真实反映区域土地利用/覆盖变化中土地利用类型变化的剧烈程度。

本研究中，土地类型数量用像元数目表示，U b 代表2009年的土地利用类型大小，U a 代表2000年的土地利用类型大小，时间T 设定为10年，计算结果如下表所示。

表格 1 灞桥区土地利用动态度

4.3 土地分类指数计算

为了考察土地利用类型变化情况，可定义各类土地利用指数，如垦殖指数、植被指数、水域指数、林地指数、城市化指数，通过每类指数的变化可定量表达该区某一类型土地利用的利用程度与变化趋势。通过考察垦殖指数、植被指数、水域指数、林地指数和城市化指数的时空变化，可以定向了解各个地区土地利用与土地覆被中的耕地、草地、林地、水域的变化以及城市的扩展状况。土地各个分类指数定义为：

∑∑==≤⨯=n

i i n

i i A a A a I 1

1

)(100)( （2）

式中：I ——分析区域的土地利用分类指数；

a i ——分析区域内i 类型土地利用所占的土地面积大小； A ——分析区域土地总面积； n ——土地利用分类的数目。

根据各个土地利用分类指数的定义，可以定义土地分类指数变化模型为:

10011

⨯⎪⎪⎪⎪

⎭

⎫ ⎝⎛-=-=∆∑∑==-A a a I I I n

i ia n i ib a b a

b （3） %100)1(⨯⨯∆=--I dI a b a b （4）

式中：Ia 、Ib ——分别为a 时间和b 时间的一定区域的土地利用分类指数； a

b I -∆——在时间段a 与b 之间的土地分类指数变化量；

a

b dI -——与t 时间段对应的土地分类指数变化率；

T ——时间段b －a 。

由于数据量和有关信息量（如实际土地类型种数）限制，不能进行灞桥区土地分类指数计算。 4.4 土地利用类型转移矩阵计算

单纯的土地类型面积的增减并不能较好地反映各类型间的转换情况，即新增部分由哪些土地资源转换而来，减少部分又去向哪里。土地利用状态转移矩阵可全面而又具体地刻画区域土地利用变化的结构