arcgis实之栅格运算
- 格式:docx
- 大小:843.31 KB
- 文档页数:15
栅格分析
从DEM上提取山顶点,凹陷点
数据
DEM数据,地理信息:KRASOVSKY-1940-Transverse-Mercator
Cellsize 5×5 单位:meter
学校选址
数据
DEM LANDUSE 栅格数据无坐标信息
School rec_sites 矢量数据点类型无坐标信息
设置栅格分析环境
Cellsize设置,输出单元格大小设置
Mask设置,分析范围设置
从Dem上提取山顶点,凹陷点
原理
DEM (数字高程模型)属于场模型(连续的栅格),表达和存储时,使用单元格中间值代替整个栅格的值,可视化时,采用插值(如线性)生成连续的栅格
采用领域的方法,可以用领域中的最高值或最低值代替该单元格的值,当该单元格是山顶点时(相对领域而言),领域计算后的值(取最高值)与原值相同,相反,凹陷点也是如
此
操作
1.加载DEM add date 或直接从catalog 中拖入
使用IDENTIFY工具查看每个单元格时,每个单元格只有一个值,对单元格而言是离散的,但对DEM而言,是连续的
2.环境设置
在10.1中,设置在ARCTOOLS的环境参数设置中完成
常规设置
输出工作空间设置
栅格分析设置
栅格存储设置
关于构建金字塔等,此处默认
领域运算
在10.1中,这个工具是focal Ststistics
在这个工具中,由两部分
NEIGHBERHOOD部分,设置领域类型,大小。类型有:矩形,园,环,楔形,和用户自定义的,大小用户自己定义
Stastistics type 就是领域的功能选择,此处选择max,既将领域最大值赋予单元格
结果预计:计算后的高值和DEM高值相等,低值应该大于等于DEM的低值
栅格运算
使用计算器计算DEM 和LINYU的差值,结果为0的值对应的栅格单元就是山顶点
重分类
利用重分类工具提取栅格为0的点,其他点赋值为NODATE
重分类之后,只保留0值
矢量化显示
将栅格转化为矢量,以便于编辑
结果
可以看出:
这些山顶点均分布在
领域内的相对高值点
如果改变领域的大
小,类型,得到的山
顶点也会发生变化
凹陷点的提取类似,使用modelbuilder提取结果
改变领域大小
10×10 100×100
可以看出,凹陷点分布在低值区域
3×3时,凹陷点和山顶点分布
绿色点为山顶点
红色点为凹陷点
这是由于:该处地势平
坦
注意
当选择同一个领域计算山顶点和凹陷点时,在栅格运算后,最大值或最小绝对值是相同的,因为这是领域的内的最大值减去最小值,或相反
学校选址
原理
学校的选址受到各方面因素的综合影响,并且各因素的影响程度不同,并且在因素内部,各个因子对选址也有不同的影响
为了在arcgis 中实现选址,规定每个因素内部分为10个等级,按对选址的影响程度分级
加载数据
由于没有投影信息,显示,为了计算方便,显示单位设置为meter,这并不影响结果,相对位置并没有改变
环境设置
常规设置
输出工作空间设置
栅格分析设置
栅格存储设置
关于构建金字塔等,此处默认
DEM分析
由于在选址中,DEM并不能直接使用,选址和坡度息息相关,越平坦,越适宜建学校Slope
单位选择度,Z-factor 高程因子,平面单位和高程单位的换算
Slope内部因子分级
使用reclassfy工具,分为10类,采用等间距分类方法,并且坡度越小值越高
Landuse 分析
Landuse 重分类
由于water ,wetland 和grass 不适宜建校,所以重分类是删除这三项
School 分析
分析现有的学校,赋予权重,距离现有学校越远,值越高 为了和其他栅格图像分析,选择距离制图生成栅格图像
距离制图
计算其它栅格到源栅格的距离
Reclassify
按照距离远近分类,越远值越高
Rec_site 分析
同学校分析
距离制图
Reclassify
距离越近,值越高
栅格运算
计算最适宜区域
结果至于在
提取大于8的部分
这就是最适合建校的区域,当然可以根据个人意愿建校
问题
在环境设置中,输出范围设置DEM 和LANDUSE为什么landuse成功,而DEM失败呢?
关于Modelbuilder的构建,如果构建MODEL,在model中,尽是reclassify的设置,即字段的赋值,而这一切尽是在没有字段值的情况下操作,这就要求我们必须事先单独操作过整个步骤,且必须对每一步非常了解,其实,这个moderl 构建中涉及太多属性设置,我认为就不应该建立model