1、空间分析是地理信息系统的主要特征。空间分析能力(特别是对空间隐含信息的提取和传输能力)是地理信息系统区别与一般信息系统的主要方面,也是评价一个地理信息系统的主要指标。空间分析赖以进行的基础是地理空间数据库。空间分析运用的手段包括各种几何的逻辑运算、数理统计分析,代数运算等数学手段。空间分析可以基于矢量数据或栅格数据进行,具体是情况要根据实际需要确定。空间分析步骤根据要进行的空间分析类型的不同,空间分析的步骤会有所不同。通常,所有的空间分析都涉及以下的基本步骤,具体在某个分析中,可以作相应的变化。空间分析的基本步骤: a) 确定问题并建立分析的目标和要满足的条件b) 针对空间问题选择合适的分析工
2、具c) 准备空间操作中要用到的数据。d) 定制一个分析计划然后执行分析操作。e) 显示并评价分析结果空间分析实际上是一个地理建模过程,它涉及:问题的确定、使用哪些空间分析操作、评价数据、以合适的次序执行一系列的空间分析操作、显示及评价分析结果。 实验数据:实验数据包括:Slope1(栅格数据),Landuse (栅格数据), landuse92,r5yield,emidalat街道图层AIOStreets和城市地籍图层:AIOZonecov气温.shp,YNBoundary.shp (云南省的边界)三、实验内容及步骤空间分析模块本章的大部分练习都会用到空间分析扩展模块,要使用“空间分析模块”首
3、先在ArcMap中执行菜单命令,在extensions模块管理窗口中, 将“Spatial Analys”前的检查框打勾。(非常重要,否则无法正常使用空间分析功能) 然后,在ArcMap 工具栏的空白区域点右键,在出现的右键菜单中找到“空间分析”项,点击该项,在ArcMap中显示“空间分析”工具栏。执行“空间分析”工具栏中的菜单命令选项设定与空间分析操作有关的一些参数。这里请在常规选项中设定一个工作目录。因为在空间分析的过程种会产生一些中间结果,默认的情况下这些数据会存储在Windows 系统的临时路径下(C:temp),当设置了工作目录后,这些中间结果就会保存在指定的路径下。空间分析工具栏1
4、. 了解栅格数据在ArcMap中,新建一个地图文档,加载栅格数据:Slope1,在TOC 中右键点击图层Slope1, 查看属性在图层属性对话框中,点击“数据源”选项,可以查看此栅格图层的相关属性及统计信息。打开“空间分析”工具栏,点击图标,查看栅格数据的统计直方图:新建ArcMap地图文档:加载离散栅格数据:Landuse,在TOC中右键点击Landuse ,“打开属性表”查看字段“Count”可以看到每种地类所占栅格单元的数目2. 用任意多边形剪切栅格数据(矢量数据转换为栅格数据)在ArcCatalog下新建一个要素类(要素类型为:多边形),命名为:ClipPoly.shp在ArcMap中
5、,加载栅格数据:Landuse、和ClipPoly.shp打开 编辑器 工具栏,开始编辑 ClipPoly ,根据要剪切的区域,绘制一个任意形状的多边形。打开属性表,修改多边形的字段“ID”的值为1,保存修改,停止编辑。打开 空间分析 工具栏 执行命令:-指定栅格大小:查询要剪切的栅格图层Landuse的栅格大小,这里指定为25指定输出栅格的名称为路径执行命令: -构造表达式:Landuse*polyClip4 ,执行 栅格图层:Landuse 和 用以剪切的栅格 polyClip4 之间的 相乘运算(polyClip4是上一步得到的图层)得到的结果即是以任意多边形剪切的Landuse数据3.
6、 栅格重分类(Raster Reclassify) 通过栅格重分类操作可以将连续栅格数据转换为离散栅格数据在ArcMap中,新建地图文档,加载栅格数据Slope1,打开“空间分析”工具栏,执行菜单命令“重分类”手动添加5个条目(Add Entry),值如上图所示,将坡度栅格重新分为5类:0 8 、8 15 、15 25 、25 35、 35 度以上。4. 栅格计算查询符合条件的栅格(Raster Calculator)找出坡度在25度以下的区域在上一步的基础上进行,执行“空间分析”工具栏上的命令:构造表达式 Slope1=25满足条件的栅格赋值为1,其余的栅格赋值为 05. 面积制表(Tabu
7、late Area)在上一步的基础上进行。加载Landuse92栅格图层,打开ArcToolbox在ArcToolbox中,执行Zonal下的“面积制表”工具按上图所示,指定分区数据和输入栅格数据打开得到的交叉面积数据表,观查其中的记录,理解本操作的意义是什么?6. 分区统计(Zonal Statistic)在ArcMap中新建地图文档,加载栅格图层 r5yield (粮食产区分类图)、栅格Organic(土壤有机质含量分布图)在r5yield 中,根据产量不同分为5个粮食产区打开ArcToolbox,执行定制 在出现的对话框中的 “命令”选项页。在左边栏中,目录列表框中,选择“工具”在右边栏
8、中,命令列表框中,选择“缓冲区向导”拖放“缓冲区向导”图标到菜单通过属性选择STR_NAME=CYPRESS ,从图层AIOStrees中,选择街道名称为CYPRESS的街道执行“工具”菜单中的“缓冲区向导”命令,或点击“缓冲区向导”图标,打开缓冲区向导对话框:通过缓冲区向导,建立所选择街道的0.05 km缓冲区(一个多边形图层),名字命名为:缓冲_AIOStreets得到沿街道“CYPRESS”的50米缓冲区8. 空间关系查询l Select By Location:根据位置选择在上一步的基础上进行,找出与 街道“CYPRESS”的50米缓冲区相交的地块。9. 采样数据的空间内插(Inter
9、polate)空间插值常用于将离散点的测量数据转换为连续的数据曲面,以便与其它空间现象的分布模式进行比较,它包括了空间内插和外推两种算法。空间内插算法是一种通过已知点的数据推求同一区域其它未知点数据的计算方法;空间外推算法则是通过已知区域的数据,推求其它区域数据的方法。数据:气温.shp 中有两个字段 Y01 Y02 记录的是16个气象观测站,2001年和2002年的年平均气温,下面要通过空间内插的方法将点上的数据扩展到连续的空间上,得到气温空间分布图。 YNBoundary.shp 是云南省的边界新建地图文档,加载图层:气温.shp 、YNBoundary, 打开“空间分析”工具栏,执行菜单
10、命令样条在样条函数内插对话框中,按下图所示指定参数确定后,得到如下的气温空间分布图(通过修改图例得到相同的效果)2001年平均气温样条函数空间内插参考以上操作,生成2002年的平均气温空间分布图:2002年平均气温样条函数空间内插执行菜单命令,在general选项页的working栏选任一个工作区,在analysis mask栏选YNBoundary.shp。切换到extent选项页,在analysis extent栏下拉选择same as layer YNBoundary,切换到cell size在analysis cell栏下拉选择maximum of inputs或minimum of
11、inputs,最后“确定”。然后重新进行空间插值(空间分析内插成栅格距离权重倒数,在input point选气温.shp,在z value field选Y02,在search radius type下拉选variable,其余默认,确定),得到如下的结果。10. 栅格单元统计(Cell Statistic)在上一步的基础上进行现在我们要根据2001年和2002年的年平均气温得到多年平均气温空间分布图,打开“空间分析”工具栏,执行菜单命令2001、2002年间平均气温空间内插11. 邻域统计(Neighborhood)邻域分析也称为窗口分析,主要应用于栅格数据模型。地理要素在空间上存在着一定的关联性。对于栅格数据所描述的某项地学要素,其中的(I,J)栅格往往会影响其周围栅格的属性特征。准确而有效地反映这种事物空间上联系的特点,是计算机地学分析的重要任务。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1