适宜性评价及空间分析.docx
- 文档编号:30280361
- 上传时间:2023-08-13
- 格式:DOCX
- 页数:30
- 大小:2.25MB
适宜性评价及空间分析.docx
《适宜性评价及空间分析.docx》由会员分享,可在线阅读,更多相关《适宜性评价及空间分析.docx(30页珍藏版)》请在冰豆网上搜索。
适宜性评价及空间分析
3、适宜性评价与空间分析
3.1研究目的和意义
土地利用研究对于当代的中国具有非常重大的意义。
改革开放以来,随着经济和社会的快速开展,对土地的开发利用导致了自然资源需求量的急剧增加,引起了土地利用/土地覆被发生了显著的变化(周万村等,2000),由此带来了诸如土地退化、水土流失、植被破坏、水资源短缺等一系列生态环境问题,这引起了国家和社会的高度关注。
从70年代后期开场,国家陆续开展了三北防护林、退耕还林还草、天然林保护等林业生态建立工程。
黄土高原地区作为我国典型的生态环境脆弱区,首当其冲成为生态环境建立的重点地区,也因此成为我国土地利用变化研究的热点地区(李丽娟等,2005;李志等,2007;郝慧梅等,2007)。
XX省X起县地处黄土高原腹地,是典型的水土流失和生态脆弱地区。
1998年以来,全县积极响应党中央“再造一个山川秀美的西北地区〞的伟大号召,在全县范围展开了封禁造林工作,于当年一次性退耕155.5万亩,实现了“封得住、退得下、还得下〞的目标,是全国封得最早、退得最快、面积最大、成效最为显著、群众得到实惠最多的县份之一,成为全国退耕还林的一面旗帜,被国家林业局命名为“全国退耕还林先进县〞。
退耕还林工程的深入实施,使X起的生态环境发生了巨大变化。
土壤年侵蚀模数由1997年的1.53x104t/km2,下降到目前的0.54x104t/km2,年降雨量已由1997年的478.3mm增至582mm,干早、暴雨、冰雹、霜冻等自然灾害明显减少,五级以上的大风已由1997年之前的年均19次降为5次,多年罕见的飞禽走兽重新显现,秀美山川初露端倪,良好的生态链正在形成。
在最新的EOS卫星遥感图片上,一片浓绿的颜色清晰地勾勒出了X起的地貌轮廓,与毗邻的XX、XX、XX地界形成了鲜明的比照。
本文选取X起县为研究区,在GIS平台上,运用退耕还林前后的遥感数据,综合分析退耕还林工程实施以来当地土地利用的时空变化,意义如下:
(l)科学评价和总结X起县退耕还林工程实施的成果,通过对其退耕还林前后土地利用时空变化及其生态效应的研究,反映其在生态建立驱动力作用下的土地利用变化特点,展现其生态建立的巨大成果。
(2)通过研究其土地利用变化时空特点及限制因子,为当地决策部门制定生态建立和土地利用政策、开展土地利用规划提供依据。
(3)由于X起县退耕还林工作非常典型,对其土地利用变化的过程和特点进展深入研究,为黄土高原地区的退耕还林、植树造林等生态建立工作的实施提供了很好的经历和样板,有利于同区的其他县份借鉴。
3.2评价标准及方法
1、评价标准
基于GIS的退耕还林分析:
陡坡地退耕还林还草的核心问题是研究地面坡度与土地利用,特别是耕地的使用状况之间的关系。
故本次研究着重于调查样区内主要土地利用类型(重点是耕地,包括滩旱地、沟旱地、梯旱地和坡地)中的地表坡度组成状况。
经查阅相关文献资料,初步确定本次工作主要任务是为黄土高原地区坡耕地退耕还林还草提供准确可靠的根本地理信息,应用GIS技术可迅速清查土地利用详细状况。
根据多年试验结果,黄土丘陵沟壑区地面坡度在3°,8°,15°,25°,35°是较为明显的土壤侵蚀临界坡度值;3°以下无明显侵蚀,8°以上的坡地细沟、浅沟普遍出现;15°是产沙的突然增大点,小于15°的坡地,土壤侵蚀相对较弱,当坡度超过15°时产沙量突然增大,侵蚀渐趋加剧,坡面发育大量切沟;25°是土壤侵蚀方式的一个转折点,25°以上的坡地开场出现重力侵蚀,并有冲沟发育;35°是黄土状堆积物的临界休止角,大于35的坡地侵蚀方式由以流水作用为主转变为以重力作用为主,崩塌、滑坡、错落等重力侵蚀现象普遍出现,侵蚀量到达最大值。
在黄土丘陵沟壑区的土地利用中,6°~25°是可修筑梯田的地域范围,也是主要的种植区域;小于6°的耕地只要在耕作时采用水保耕作法便可以防止面蚀发生;25°以上的坡地应作为林牧用地,是需进展退耕还林还草的区域。
综上所述,可以得到黄土丘陵沟壑区的坡度共有5级:
<6°,6°~15°,15°~25°,25°~35°,>35°。
退耕还林还草属水土保持中的生物措施,其根本要求是6°~15°耕地作为主要种植区域;15°~25°耕地需加强工程措施以防止水土流失发生;25°以上坡地作为林牧用地区域,需退耕还林还草;35°以上坡地是重点治理区域,必需尽快进展治理。
按照上述要求,根本确定了GIS在本次研究中的任务如下:
〔1〕叠置分析:
叠置分析是提取空间隐含信息的一种手段,其操作为将相关数据进展叠加产生一个新的数据层面,要求被叠加的数据是基于一样坐标系统的一样区域,叠加的结果层将综合被叠加层的要素具有的属性。
叠置分析包括矢量数据和栅格数据的叠置。
在ArcGIS中使用intersect命令进展叠置分析,将耕地图层与坡度图层、坡向图层、高程带图层分别叠加分析,得到耕地的坡度、坡向、高程属性。
进而对各种坡度下的耕地适宜性按照预先设定的指标体系进展评价。
〔2〕进展面积统计。
将ENVI中运用决策树分类方法分好类的耕地图导入到ArcGIS中,首先进展矢量栅格数据的转换操作,然后在属性表中加一个面积属性〔AddField操作〕,进而计算各种类型耕地的面积以便后续的评价。
〔3〕叠加河流、行政区划、道路、居民点等根底地理信息,在ArcGIS中运用地图整饰的方法生成专题图。
2、评价方法
〔1〕叠置分析
Identity操作
(1)在ArcToolbox中选择AnalystTools|Overlay|Identity选项,翻开其对话框。
(2)填入输入图层〔Inputfeatures〕,识别参照图层〔Erasefeature〕,输出图层〔OutputFeatureClass〕和选择需要连接过去的属性字段,在右下角的环境设置〔Environments〕中,可以对输入输出数据的参数进展设置,参数设置见下列图。
(3)单击OK按钮,完成操作,结果见下列图。
Union操作
(1)在ArcToolbox中选择AnalystTools,翻开后选择Overlay中的Union选项,翻开其对话框。
(2)逐个输入要进展合并的图层〔InputFeature〕,按右边的“加号〞来将图层添加进来,在中间“Features〞组合框内的就是要进展合并操作的图层列表,输入要输出的文件的路径和名称〔OutputFeatures〕,同时在下方的属性字段中选择要进展连接的属性字段〔JoinAttributes〕或全部,输出文件的类型,也可以对环境参数进展相关的设置。
(3)单击OK进展合并操作,输出结果。
图1Union操作对话框
Intersect操作
(1)从ArcToolbox中选择AnalystTools,翻开后选择Overlay中的Intersect选项,翻开其对话框。
(2)逐个输入要进展相交的图层〔Inputfeatures〕,按右边的“加号〞来将图层添加进来,在中间“Features〞组合框内的就是要进展相交操作的图层列表,输入要输出的文件的路径和名称〔OutputFeatureClass〕,同时在下方的属性字段中选择要进展连接的属性字段〔JoinAttributes〕或全部,设置输出文件的类型,也可以对环境参数进展相关的设置。
(3)单击OK进展交集操作,输出结果。
图2Intersect操作对话框
Symmetricaldifference操作
(1)在Arctoolbox中选择AnalystTools,翻开后选择Overlay中的SymmetricalDifference选项,翻开其对话框〔见下列图〕。
(2)输入要进展操作的输入图层〔InputFeatures〕,同时在下面输入参照的差值图层〔UpdateFeatures〕,输入要输出的文件的路径和名称〔OutputFeatureClass〕,同时在下方的属性字段中选择要进展连接的属性字段〔JoinAttributes〕或全部,还有容错量,也可以对环境参数进展相关的设置。
(3)单击OK进展交集操作,输出结果如下列图所示。
〔2〕栅格计算器
A、翻开栅格计算器功能
B、选定一个图层后,编辑计算公式
C、单击evaluate即可计算得到结果
〔3〕重分类
新值替代
在SpatialAnalyst的下拉菜单中选择Reclassify;
在Inputraster的下拉菜单中选择需要变更值的图层;
在ReclassField栏的下拉菜单中选择需变更的字段;
在NewValue栏中定位需要改变数值的位置,然后键入新值。
可点击Load按钮导入已经制作好的重映射表,也可以点击Save按钮来保存当前重映射表;
ChangemissingvaluestoNoData:
可选项,将missingvalues值改成无数据〔NoData〕;
Outputraster:
为输出结果指定目录及名称;
点击OK按钮。
旧值合并
在SpatialAnalyst的下拉菜单中选择Reclassify;
在Inputraster的下拉菜单中选择需要重新组合分类的图层;
在ReclassField栏的下拉菜单中选择需用的字段;
在OldValue栏中连续选择要进展归并的数据〔按住Ctrl键可以逐个进展选择,按Shift键可以选择两者之间的全部数据〕;
在选择并弹亮的数据项上点右键;在弹出的菜单中选择GroupEntries;
给新组合成的数据赋予新值。
可点击Load按钮导入已经制作好的重映射表,也可以点击Save按钮来保存当前重映射表;
ChangemissingvaluestoNoData:
可选项,将missingvalues值改成无数据〔NoData〕;
Outputraster:
为输出结果指定目录及名称;
点击OK按钮。
重新分类
在数据的使用过程中经常会因某种需要将数据用一种等级体系来进展分类,或需要将多个栅格数据用统一的等级体系重新归类。
例如,在对洪水灾害进展预测时,需要分析降雨量数据,地形数据和植被数据等。
首先要按照这几种数据对洪灾带来的影响分别将其分为大概10类〔1-10〕。
每种数据中,级别越高其对洪灾的影响度越大。
经过分类后就可以通过这些分类信息进展洪灾模拟。
下面为分类根本操作:
在SpatialAnalyst的下拉菜单中选择Reclassify;
在Inputraster的下拉菜单中选择需要重新分类的图层;
在ReclassField栏的下拉菜单中选择需用的字段;
点击Classify按钮;
在Method栏的下拉菜单中选择一种分类方法,包括Manual〔手工分类〕,EqualInterval〔等间距分类〕,DefinedInterval〔自定义间距分类〕,Quantile〔分位数分类〕,NaturalBreaks〔自然间距分类〕,StandardDeviation〔标准差分类〕;
在Classes栏的下拉菜单中选择进展重分类的类别数目;
点击Classification对话框中的OK按钮;
根据需要,在Reclassify对话框中修改将输出的栅格数据的值,或点击Load按钮导入已经制作好的重映射表,也可以点击Save按钮来保存当前重映射表;
Outputraster:
为输出结果指定目录及名称;
点击OK按钮。
重分类后的图像如下列图所示
3.3操作步骤
1、从DEM数据中提取坡度数据。
2、将提取的坡度数据进展重分类。
黄土丘陵沟壑区的坡度共有5级:
<6°,6°~15°,15°~25°,25°~35°,>35°,按照上述标准将提取的坡度数据进展重分类。
地外表任一点的坡度是指过该点的切平面与水平地面的夹角,表示了地外表在该点的倾斜程度〔葛翠萍等,2021〕。
坡度作为地形的一个根本因子,可以通过调节太阳辐射!
水分等自然要素的分配,间接影响土地利用类型的空间分布。
另一方面,坡度的不同使得人类活动对土地利用方式的干扰程度不同,从而也影响到土地利用方式的形成和变化〔任志远等,2005〕。
在黄土高原地区,坡度是重要的地形指标,与退耕还林等引起土地利用变化的各项措施的实施有着密不可分的联系。
坡度分级确实定:
坡度分级的原那么是能够有效反映研究区坡度组合的特点,本研究以国土资源部?
土地利用更新调查技术规定〔试行〕?
〔2004〕为根底,结合X起县地形地貌的特点,适当参考相关研究的成果〔汤国安等,2006;胡世雄等,1999;赵牡丹等,2002〕,将本区坡度范围定为5级,如下表:
表1X起县坡度分级表
坡度分级
平缓坡
缓坡
斜坡
陡坡
极陡坡
坡度范围〔度〕
0~6
6~15
15~25
25~35
≥35
具体操作步骤如下:
重分类的结果图
3、坡度数据的栅格转换。
将坡度数据转换为矢量格式,以便于与两期土地利用数据进展叠置分析。
具体操作如下:
栅格转矢量结果图
4、叠置分析。
将两期遥感数据进展叠置〔叠置分析中的Identity操作,以便更好的识别出退耕还林的情况〕。
以水平农田为例,具体操作步骤如下:
5、将第4步中的两期遥感影像叠置的结果和坡度图层再进展叠置分析。
Intersect操作的结果图
林地草地两期数据叠加效果图
耕地两期数据叠加效果图
耕地叠加坡度效果图
林草地和坡度Intersect结果图
居民点叠加效果图
6、导出属性表,进展每一级坡度的各类用地变化情况分析。
在Excel中分析不同坡度等级下各类用地的变化情况。
首先将X起县1:
1万的DEM数据在Arcgis10平台上利用立方卷积法进展像元重采样,设定分辨率为30x30m,以使分辨率与土地利用数据保持一致。
然后应用ArcgIs中的外表分析〔surfaceAnalyst〕模块,提取坡度〔slope〕数据,将坡度值按上述分级标准分为5级。
最后,把栅格的坡度分级数据转为矢量数据,与土地利用变化图进展叠加分析,即可得到不同坡度梯度下的土地利用变化数据,如下表所示:
表2不同坡度梯度下的X起县土地利用状况〔hm2〕
坡度范围〔度〕
0~6
6~15
15~25
25~35
≥35
林草地
1991
16202.7
35591.2
84479.5
60445.1
28858.2
2021
20597.5
53127.7
130043.8
92496.2
30302.4
耕地
1991
13731.5
27604.6
57870.7
38997.4
18833.7
2021
7143.2
8114.5
9296.3
4437.8
1332.2
居民点
1991
102.8
77.5
102.2
58.3
19.3
2021
1115.9
765.2
730.4
359.2
194.
水域
1991
294.5
194.1
359.3
289.8
150.3
2021
456
256.4
587.2
369.5
176.2
3.4各类土地利用变化情况分析:
〔1〕林草地
林地在1991年和2021年不同坡度下的分布状况是一致的,其特点为低坡度级的林地分布较少,中高坡度级的林地分布较广,74%以上的林地分布在坡度为15°以上的地区,优势坡度范围均为15°~25°从12年来林地在坡度梯度的变化来看,各坡度级的林地面积均有增加,其中中高坡度级的林地增量最多,到达了总增幅的75.09%,优势坡度级〔15°~25°〕的优势程度进一步增大。
林地变化的这些特征说明了当地开展林业生态建立主要集中在坡度较高的山区。
草地在两个年度不同坡度下的分布特征与林地一致,表现为中高坡度范围的草地分布较广泛,坡度15°以上的草地面积均到达当年草地总面积的77%以上,优势坡度范围也为150~25°。
从1991年到2021年,草地面积整体表达减少的趋势,但是在坡度为35。
以上地区草地面积略有增加,这与当地开展的封山禁牧等生态保护政策因素有关。
其在不同坡度级上的变化趋势根本上与林地呈反向互补关系,所以草地在各坡度梯度的面积减少并没有造成当地生态环境的恶化。
〔2〕耕地
在1997年和2021年,X起县的耕地在不同坡度级均呈倒U型分布,15°~25°坡度级的耕地面积最多,为这两年的优势坡度。
区别在于1991年陡坡以上的耕地面积占当年总耕地面积的比例较大,为36.06%,远高于2021年的19.37%。
这22年间,耕地总体呈减少趋势,各坡度梯度的耕地面积也表现为减少,中高坡度级〔15°以上〕减幅较大,占到耕地减少面积的80%以上,是造本钱区耕地减少的决定因素,由此说明X起县退耕还林工程主要针对陡坡耕地进展退耕。
〔3〕居民点
1991年和2021年,X起县居民点占地面积随坡度的增加大致呈递减的趋势,中低坡度范围〔25°以下〕分布较多,占总建立用地面积的79%以上,优势坡度范围为0°~8°,22年间居民点面积大幅增加,增加的居民点也主要分布在中低坡度。
3.5X起县土地利用变化驱动力分析
3.5.1退耕还林
X起县自1998年开场实施“退耕还林〞工程,是全国150多个退耕还林县(市、区)封得最早、退得最快、面积最大的县份。
其特点为“一次性退耕〞,1999年一次性退耕坡耕地155.5万亩。
由于退耕还林工程的实施,X起县耕地面积尤其是坡耕地的面积大幅减少,林地面积大幅增加,由此造成了相关土地利用类型的剧烈变化,是本区土地利用变化的主要驱动力。
;从空间上来说,退耕还林工程实施的空间分布非常广泛,反映空间分布的多度值排在所有驱动因素的第一位;其对当地土地利用变化的重要程度仅次于荒山造林、人工种草与自然封育,位列各驱动因素的第三位。
3.5.2荒山造林
除退耕还林以外,X起县还在中低覆盖度草地、未利用地上开展了以生态水保林和生态经济林为主的人工造林活动。
生态水保林包括XX杨、刺槐、侧柏、油松等地方性乔木树种;生态经济林包括沙棘、山杏、山桃等具有较高经济价值的灌木。
这些林木的栽植,不仅取得了良好的经济效益,而且取得了巨大的生态效益,使得当地土地利用发生了巨大变化。
3.5.3人工种草与自然封育
1998年5月,X起县作出了实施封山禁牧、大力开展舍饲养羊的决定,并于当年淘汰出栏当地土种山羊23.8万只,顺利实现了整体封禁目标,在全国首开了封山禁牧的先河。
其目的在于通过自然的作用,使植被得以恢复。
除此以外,当地也开展了以牧草为主要草种类型的人工种草植草活动,到2005年底,累计种植人工牧草127万亩,保存面积%万亩,是全国人工种植紫花首藉面积最大的县。
这些措施的实施,引起了当地土地利用的较大变化:
共有722%.7hm2土地因人工种草和自然封育而发生转化,占因生态建立引起的土地利用变化总面积的30.26%,由此使该地的草地覆盖率提高了19.1个百分点;其多度值为26.7,说明此项驱动因素在空间范围的分布也较为广泛;其重要值为56.8,说明人工种草与自然封育因素在当地土地利用变化的过程中也起着很重要的作用。
3.6土地利用的潜力
土地利用是人类为了到达一定的目的或者满足一定的需求开发利用土地资源的具体表现形式,其变化过程表达了人为因素在其中所起的重大作用。
同时,土地既然作为一种自然资源,它的变化与变迁必然受到自然规律的制约。
所以,土地利用变化是一个自然和人文两种因素共同作用的过程。
在X起县,当地为了实现改善生态环境的目标,人为开展了大规模的退耕还林工作,使当地土地利用发生了巨大变化,但是这种变化能否一直持续下去?
或者说这种变化有没有极限?
这是本章需要探讨的问题。
目前,有关土地利用变化预测的研究也很多,但大局部仅仅是从统计学的角度利用各种模型(如马尔科夫模型)进展预测,得出的结果仅仅是在统计学上有意义,并不能够反映人类活动和自然环境对土地利用的限制作用,所以对当地的土地利用政策制定和规划并没有特别实际的指导意义。
本文选取林地和耕地两种土地利用类型,通过分析自然和人为因素中对其起最大限制作用的因子;来综合反映这两种地类变化的潜力和极限,从而为当地今后的退耕还林政策的制定提供依据和参考。
3.7X起县生态建立的启示
X起县土地利用变化的结果说明当地生态建立取得了巨大成绩,其经历应当成为黄土高原其他地区学习和借鉴的范例。
依据本次研究结果,将其生态建立的经历和特点总结为以下两点:
(l)生态建立要遵循自然规律,因地制宜、因时制宜。
由土地利用空间变化分析可知,当地各土地利用类型的分布与变化受自然规律的制约,尤其是坡度、坡向、高程等因素影响着水分、光照的再分配,在生态建立时需要重点考虑。
(2)生态建立政府是主导,但要注意调动普通民众的积极性,以确保生态建立的成果。
无论是退耕还林还是封山禁牧,都需要政府的正确决策和相关配套措施的出台,并且还要积极思考退耕禁牧以后农民和牧民的生计问题、去向问题,只有这样才能保证生态建立不反复。
X起县的实践经历和我们土地利用变化的研究结果均说明,当地将生态建立与新农村建立严密结合,一方面强化退耕还林等生态建立工程的实施,另一方面着力解决农民的收入问题,大力开展新农村建立,使得生态建立与农民生活取得双丰收。
3.8讨论
(l)通过对X起县土地利用时空变化的研究,我们发现局部林草地存在退化现象,陡坡耕地在局部地区也还存在,这是今后当地土地利用规划和改造的重点。
(2)本文主要采用的遥感影像分别为TM和Landsat数据,两种数据源的分辨率有较大差异,这样就有可能因数据源的不同产生误差。
但是,随着遥感技术的开展,数据的分辨率会越来越高,应用高分辨率数据进展土地利用研究也将会是开展趋势,那么如何解决高、低分辨率数据的匹配问题,以及如何定量的评估不同数据源带来的计算误差是今后相关研究的新课题。
(3)由于受到数据源分辨率的影响,为了确保一定的准确率,本文在研究与分析时仅采用了土地利用的一级分类,对二级分类的土地利用类型没有过多涉及,在今后的研究中,应重点探讨一级分类下的不同二级分类之间的转化,这对于分析土地利用的自然和人文驱动力无疑具有重要意义。
(4)在定量计算土地利用变化的生态效应时,运用的是全国的平均数据,众所周知,黄土高原地区是全国较为特殊的区域,全国的平均数据肯定会存在较大误差,所以在今后的研究中,应当探寻适合当地的核算因子,以提高准确度。
(5)通过研究,作者认为X起县退耕还林的做法和经历局部可供黄土高原其他地区借鉴,比方退耕还林的思路、技术及相关配套政策等,但同时其退耕还林是在大量的石油收入保证下才得以顺利开展,对于地区经济较为薄弱的黄土高原其他县份,如何保证退耕还林后的成果,或者协调退耕还林与增加农民收入、促进经济转型等关系,还需要继续探讨。
致谢
这个学期的资源与环境遥感,让我们感觉收获很大,通过对这门课的学习,我觉得我对生命科学这个专业有了更深的了解。
自从黄健教师给我们上完第一节课后,我才明白,如今生命科学已经成为推动社会开展,经济进步十分重要的一门学科。
在这次
参考文献
[1]葛翠萍,赵军,尹升葛,等.典型黑土区地形因子对土地利用变化的影响[J].农业系统科学与综合研究,2021,25(l):
69-73.
[2]任志远,黄青,李唐.XX省生态平安及空间差异定量分析[J].地理学报,2005,60(4):
597-606.
[3]汤国安,宋佳.基于DEM坡
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 适宜 评价 空间 分析
![提示](https://static.bdocx.com/images/bang_tan.gif)