土地利用变化报告.docx
- 文档编号:10982023
- 上传时间:2023-02-24
- 格式:DOCX
- 页数:13
- 大小:754.44KB
土地利用变化报告.docx
《土地利用变化报告.docx》由会员分享,可在线阅读,更多相关《土地利用变化报告.docx(13页珍藏版)》请在冰豆网上搜索。
土地利用变化报告
土地利用变化报告..
土
地
利
用
变
化
报
告
1前言
经过一学期对《土地利用与植被覆盖变化研究》课程的学习,使我对这门课程有了更进一步的认识和了解。
为了巩固所学内容,灵活的将理论与实际相结合,本次研究特选取西安市临潼区2000、2009年的tm影像,利用所学理论知识,结合软件操作,完成土地利用分类、动态度计算、土地利用转移矩阵、土地利用程度综合指数计算和马尔科夫预测等工作,并对研究结果进行分析总结。
本次研究旨在掌握土地利用变化研究的基本流程,从而加深对遥感与土地利用基本理论的理解,着重培养我们分析问题和解决问题的能力。
2研究过程
2.1土地利用分类
土地利用分类是区分土地利用空间地域组成单元的过程。
这种空间地域单元是土地利用的地域组合单位,表现人类对土地利用、改造的方式和成果,反映土地的利用形式和用途(功能)。
分类提供土地覆盖以及土地利用中涉及的人类活动类型等信息。
它也可能有助于环境影响评价,以及潜在的土地利用多样性。
本研究采用监督分类的方法对临潼区2000年和2009年土地利用情况进行分类,以分析研究区的土地利用分布格局。
根据已有的影像资料分析判别,可将研究区地物类别划分为四类:
建设用地、耕地、水体和林地。
在ENVI中利用ROITool来定义训练样本,也就是将感兴趣区当做训练样本。
根据人工经验知识和影像波段组合,可确定影像中白色有尖锐颗粒的呈片状分布地区为建设用地,影响颜色呈浅绿、深绿,并且纹理清晰,分布在城市着农村周围的平缓区域的为耕地,水体的分布呈现出带状或颜色均一的片状,一般为蓝色或浅蓝色,而在影像东南地区的山上则分布着大量的林地,颜色为墨绿,色泽均一或有浅绿色颗粒。
根据上述的解译特征,在ENVI软件中的RegionofInterest中选择土地类别的样本,并按照3、7的规则分为验证样本和训练样本。
如图所示:
图1.样本选取过程
对于选好的训练样本在Options下的ComputeROISeparability中计算可分离度,一般可分离度大于1.8,则表示训练样本的可分离性较好。
本研究所选的分类样本分离度均大于1.6,表明可分离性较好。
根据分类的复杂程度和精度要求,本次研究选择最大似然法分类器。
最大似然法分类器假设每一个波段的每一类统计都呈正态分布,计算给定像元属于某一训练样本的似然度,像元最终被归并到似然度最大的一类当中。
在ENVI软件菜单中选择Classification——Supervised——LikelihoodClassification功能,对2000年和2009年的遥感影像进行土地利用分类,制作成专题地图如下:
图2.2000年、2009年土地利用分布图
执行监督分类后,需要对分类结果进行评价。
用之前选好的验证样本对分类后的影像进行精度评价。
使用ConfusionMatrix工具可以把分类结果的精度显示在一个混淆矩阵,用于比较分类结果和地表真实信息。
ENVI可以使用一幅地表真实图像或地表真实感兴趣区来计算混淆矩阵。
在Classification——PostClassification——ConfusionMatrix——UsingGroundTruthImage下完成精度评定的工作,得到下面的报表:
表12000年和2009年土地利用精度评价表
2000年土地利用精度检验
类型
建设用地
耕地
水体
林地
总计
建设用地
2275
72
0
159
2506
耕地
2
5110
71
253
5436
水体
307
15
2260
6
2588
林地
363
232
0
2612
3207
总计
2947
5429
2331
3030
13737
总精度=89.2262%;Kappa=0.8506
2009年土地利用精度检验
类型
建设用地
耕地
水体
林地
总计
建设用地
1662
164
269
0
2095
耕地
240
3183
0
597
4020
水体
150
0
1072
0
1222
林地
4
602
239
1254
2099
总计
2056
3949
1580
1851
9436
总精度=75.9962%;Kappa=0.6610
2.2土地利用变化分析
土地覆盖主要表达的是土地表层以植被为主要覆盖物的自然类型和状态;土地利用则是人类对土地的利用方式和使用状况。
因此,土地利用变化主要体现在土地利用类型变化、土地利用类型数量变化、土地资源生态背景质量变化、土地利用程度变化及土地利用变化的区域差异等方面。
在全球环境变化研究中,土地利用和土地覆被动态越来越被认为是一个关键而迫切的研究课题。
土地利用转移矩阵可全面而又具体地分析区域土地利用变化的数量结构特征与各用地类型变化的方向,因而在土地利用变化和模拟分析中具有重要意义,并得到了广泛应用。
但是目前对土地利用转移矩阵的应用主要还是集中在直接分析土地利用类型面积的变化信息上。
土地利用转移矩阵来源于系统分析中对系统状态与状态转移的定量描述。
如表2中的土地利用转移矩阵中行表示2000年土地利用类型,列表示2009年土地利用类型。
非对角线上百分比表示2000年-2009年期间一种土地类型转换为另一种土地类型的面积占2009年该种土地利用类型总面积的百分比;对角线上的百分比表示2000年-2009年期间土地利用类型保持不变的面积百分比。
表22000-2009年土地利用转移矩阵
2000-2009年土地利用转移矩阵
面积(Km²)
2000年土地利用类型
建设用地
耕地
水体
林地
总计
2009年土地利用类型
建设用地
63.3855
116.9687
21.0657
20.3223
221.8496
百分比
28.57
52.72
9.50
9.16
100
耕地
70.7224
423.9598
13.7578
85.4306
594.1856
百分比
11.90
71.35
2.32
14.38
100
水体
1.8663
5.1144
4.2806
0.6524
11.9181
百分比
15.66
42.91
35.92
5.47
100
林地
14.2299
27.7788
2.3925
43.4614
87.9535
百分比
16.18
31.58
2.72
49.41
100
总计
150.2041
573.8217
41.4966
149.8667
1447.3822
从表中可以看出,2009年的建设用地有28.57%对比2000年保持不变,有52.72%是由耕地转换来的,有9.5%和9.16%是由水体和林地转换而来的。
而耕地除了保持不变的部分,分别有11.9%、2.32%和14.38%是由建设用地、水体和林地转化而来的。
对水体而言,2009年相比于2000年有42.91%是由耕地转换而来的。
林地有31.58%是由耕地转换而来的,这可能与国家退耕还林的政策有关。
总之,十年来临潼区的土地利用转换较大。
但总体上建设用地和耕地是呈增长趋势,而水体和林地等依靠土地自身再生的资源却呈减少趋势,值得有关部门和广大人民群众关注和保护。
土地利用动态度对比较土地利用变化区域差异和分析、预测土地利用变化趋势,具有十分重要的作用。
利用土地利用动态度模型分析土地利用类型的动态变化,可以反映区域土地利用类型变化的剧烈程
度。
土地资源数量变化可以用土地利用动态度来表示。
单一土地利用类型动态度表达的是某研究区一定时间范围内某种土地利用类型的数量变化情况,其表达式为:
式中,K为研究时段内某一土地利用类型动态度;Ua、Ub分别为研究期初及研究期末某一种土地利用类型的数量;T为研究时段长,当T的时段设定为年时,K的值就是该研究区某种土地利用类型的年变化率。
采用上述公式,我们对西安市临潼区的土地利用类型数量在ENVI软件中进行了统计,并对其动态度进行了计算,得到以下结果:
表3单一土地利用动态度计算
单一土地利用类型动态度
土地类型
2000
2009
年份间隔
动态度
建设用地
1503008
2218496
10
4.76
耕地
5741017
5941856
10
0.35
水体
415464
119181
10
-7.13
林地
1499403
879535
10
-4.13
由表3可知,土地利用的动态度有正数也有负数,但总体而言变化率并不是很大。
水体和林地的动态度为负数,变化较大,表明该土地利用类型减少较快,说明水体和林地开发速度较快;而建设用地和耕地则相反,其动态度均为正数,表明该土地利用类型有增加的趋势,这主要是随着城镇化进程的推进,必须占用土地以适应城镇化的发展。
而耕地的增长主要表现在对自然林地的改造利用上,可能是由于城镇化的发展和人口基数的快速增加,致使耕地面积不断扩张。
土地利用程度主要反映土地利用的广度和深度,它不仅反映了土地利用中土地本身的自然属性,同时也反映了人类因素与自然环境因素的综合效应。
根据刘纪远先生等提出的土地利用程度的综合分析方法。
将土地利用程度按照土地自然综合体在社会因素影响下的自然平衡状态分为若干级,并赋予分级指数,从而给出了某研究区土地利用程度综合指数的定量化表达:
式中,Dj为某研究区域土地利用程度综合指数;Ai为研究区域内第i级土地利用分级指数;Ci为研究区域内第i级土地利用程度分级面积百分比;n为土地利用程度分级数。
按照土地利用程度指数的分级标准,城乡建设,交通水利建设等建设用地的土地利用程度最高为第一级,第二级为耕地、园地以及其他农用地,第三极的是林地、牧草地和水域,土地利用程度最低的是自然保留地为第四级(见表4)。
西安市临潼区内土地面积所占比例最大的是农业用地中的耕地(属于土地人为再生利用级),占土地总面积的60%左右;其次为建设用地(属于土地非再生利用级);然后是未经开发,自然生长的林地(属于土地未利用级);最后是水域(属于土地自身再生利用级)。
表明临潼区的土地资源开发较好。
表4临潼区土地利用分级
类型
土地未利用级
土地自身再生利用级
土地认为再生利用级
土地非再生利用级
分类指数
1
2
3
4
土地利用类型
林地
水体
耕地
建设用地
根据土地利用程度综合指数模型计算的2000年和2009年的土地利用程度综合指数为下表:
表5土地利用程度综合指数
土地类型
2000
2009
像元数
分类等级
面积比
综合指数
像元数
分类等级
面积比
综合指数
建设用地
1503008
4
16.41
279
2218496
4
24.22
304
耕地
5741017
3
62.68
5941856
3
64.87
水体
415464
2
4.53
119181
2
1.30
林地
1499403
1
16.37
879535
1
9.60
总计
9158892
9159068
从上表中可以看出:
临潼区土地利用程度整体处于一个中等偏上水平。
并且2000年到2009年土地利用程度总体变化差异不大。
土地利用程度是人类有意识的一种选择行为活动。
土地利用程度越高,其自然属性的可利用程度就越低,人类对其干扰程度就越大。
通过采用土地利用程度的综合分析方法对临潼区进行研究,表明:
临潼区土地利用处于发展时期,土地利用程度整体处于一个中等偏上水平,仍有很多土地资源有待进一步开发利用,土地利用程度依地形地貌特征差异分布。
土地利用程度主要反映了土地系统中人类活动因素的影响程度,但这种影响程度常常受地理区位、经济政治和交通运输因素的限制。
因此需要大量的数据做进一步的研究分析。
2.3马尔科夫预测
马尔科夫预测过程是根据俄国数学家马尔科夫的随机过程理论提出来的,是一种具有“无后效应”的特殊随机过程。
所谓的“无后效应”是指某随机过程时刻(t+1)的状态S(t+1)的概率只与时刻t的状态S(t)有关,而与时刻t之前的状态概率无关。
这点用于土地利用结构动态变化预测是合适的,因为在一定条件下,土地利用的动态演变具有马尔科夫过程的性质。
具体而言,在通过转移矩阵生成区域土地利用变化转移概率矩阵的基础上,马尔科夫模型被用来预测特定情景下土地利用变化的未来趋势。
(1)转移概率矩阵的确立是成功应用马尔科夫模型的关键,转移概率矩阵可以定量说明土地利用类型之间的相互转化状况,揭示转移速率,其土地利用/覆被类型转移概率矩阵的数学表达式为:
式中:
Pij—土地利用/覆被类型i转化为j的转移概率。
该矩阵的每一项元素需满足以下两个条件:
(2)假设己知土地利用/覆被系统中n=0时的,初始状态
,经过n次转移后的状态
为:
式中:
—初始状态土地利用/覆被概率矩阵经过n步转移后的概率矩阵,经过n步转移后的系统状态概率为
,所以
。
由此可见,系统在任何时刻的状态概率是由初始状态概率和转移概率确定的。
将土地利用/覆被类型划分为一系列相互转化的状态,各状态在系统中所占的份额即每种土地利用类型而积占全部土地利用类型面积的百分比作为各状态的初始概率,从而构成初始状态矩阵A0。
用ENVI软件中的ClassStatistics工具统计临潼区2000年土地利用类型面积,计算得到各覆被类型面积百分比,形成初始状态矩阵,见表6:
表6初始状态矩阵
土地利用类型
A0
建设用地
16.41
耕地
62.68
水体
4.54
林地
16.37
转移概率矩阵的计算:
利用己求出的某时间段内的某地土地利用/覆被类型而积的转移矩阵,求出该时间段内该地土地利用/覆被类型的转移概率,即某地类转化后的各土地利用类型而积占转化前该地类的百分比。
在ENVI中的ChangeDetection中可以计算得到转移概率矩阵如表7:
表72000-2009年的转移概率矩阵
2000
2009
建设用地
耕地
水体
林地
建设用地
0.42172
0.47054
0.01242
0.09468
耕地
0.20374
0.73848
0.000891
0.04839
水体
0.50704
0.33114
0.10303
0.05759
林地
0.13554
0.56976
0.000453
0.28986
利用Matlab软件,以9年为一个步长,根据公式建立马尔科夫预测模型。
由此可计算出今后每个阶段的土地利用/覆被类型的转移概率矩阵和各土地利用/覆被类型所占的面积比例,即可模拟出各类地物之间面积比例的转化情况。
得到2018年的土地利用分布情况,如表8:
表8马尔科夫预测2018年土地利用分布情况
土地利用类型
A0
建设用地
24.23
耕地
64.86
水体
1.32
林地
9.59
3心得
“纸上得来终觉浅,绝知此事要躬行”,《土地利用与植被覆盖变化》的研究实习结束了,经过这次系统全面的学习,让我更深刻的了解了遥感图像处理和土地利用变化监测的的操作流程,也掌握了面对常出现问题的解决方法,同时也让我对这门课程有了更确切的理解,把理论应用于实际,同时加深了对理论部分的重点理解,还提升了动手能力。
当然还有很多不懂的问题需要日后的不断探索和研究,但这次实践给我打开了一扇自主学习的门,让我懂得了凡是不要怕难不要怕不会。
所谓“困难像弹簧,看你强不强,你强它就弱,你弱它就强”,当我们遇到不会或者不清楚的地方,要有效利用周围的资源,网上咨询,询问老师,只要有心,就一定会有解决之法。
总之,这次研究实习让我受益匪浅。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 土地利用 变化 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)