第11章 栅格数据的空间分析

ArcGIS10地理信息系统教程教学PPT牟乃夏Email：mounaixia微博：新浪博客：QQ：153807203京东：,第11章栅格数据的空间分析,本章主要内容11.1栅格数据的基础知识11.2数据分析的环境设置11.3密度分析11.4距离分析11.5提取分析11.6栅格插值11.7重分类11.8条件分析与栅格计算器11.9太阳辐射分析11.10表面分析11.11统计分析11.12实例,3,第11章栅格数据的空间分析,基于栅格数据的空间分析在ArcGIS中占有重要地位，空间建模的基本过程也是通过栅格数据的空间分析进行的。本章讲述栅格数据的基础知识、空间分析的环境设置、密度分析、距离分析、提取分析、栅格插值、重分类、条件分析、太阳辐射分析、表面分析、统计分析等基本概念和操作方面。,4,11.1栅格数据的基础知识,栅格数据是由按行和列组织的单元矩阵组成，每个单元都包含一个信息值。栅格数据一般分为两类：专题数据和图像数据。,5,11.1.1栅格数据的组成,1.单个栅格数据由栅格单元组成，单元是特定区域的方块，所有单元大小相同。2.值每个单元被分配一个特定的值以标识或描述单元归属的类或组，或所描述现象的大小或数量。3.分区和区域具有相同值得任意两个或多个像元属于同一个分区。,6,11.1.2栅格数据的应用,1.栅格数据用作底图栅格数据通常被用做其他要素图层的背景来显示场景。2.栅格数据用做表面地图栅格数据非常适合表示沿表面连续变化的数据。3.栅格数据用做主题地图表示主题的栅格数据可通过分析其他数据获得。4.栅格数据用做要素的属性栅格结构的数字照片、扫描文档还可以用做地理对象的属性。,7,11.2.1为分析结果指定磁盘位置,分析结果的缺省位置是系统的临时目录，可以为分析结果指定新的存放位置，操作步骤如下：1）打开【环境设置】对话框，单击【工作空间】标签，如下图所示。2）输入【当前工作空间】、【临时工作空间】的存放路径。3）单击【确定】按钮，完成设置。,分区和区域,8,11.2数据分析的环境设置,ArcGIS10的环境设置有四个级别，分别为应用程序级别设置、工具级别设置、模型级别设置和模型流程级别设置。以第一种方法为例，说明设置步骤。在ArcMap窗口中单击模型构建器窗口按钮，打开【模型】对话框。单击【模型】-【模型属性】，打开【模型属性】对话框。切换到【环境】选项卡，选中要设置的环境旁边的复选框，单击【值】按钮，打【环境设置】对话框进行设置。,11.2.2使用分析掩膜,分析掩膜标识了分析过程中需要考虑到的分析单元即分析范围。设置分析掩膜分两步：首先创建分析掩膜，可通过提取分析工具创建；然后在【环境设置】对话框中的【栅格分析】标签中指定，使之能应用于后续的分析。掩膜提取结果如左图所示。,9,掩膜提取结果,10,11.2.3选择坐标系统,在ArcGIS的空间分析中，可以指定分析结果的坐标系统，操作步骤如下：1）在【环境设置】对话框中单击【输入坐标系】标签。2）在【输入坐标系】下拉框中选择坐标系统。通常情况下，分析结果将使用第一个输入栅格数据集的坐标系统。3）单击【确定】按钮，完成设置。,11.2.4设置分析结果的范围,在栅格数据的空间分析中，分析范围由所用的实用工具决定。操作步骤如下：1）在【环境设置】对话框中单击【处理范围】标签，如下图所示。2）在【范围】下拉框中选择空间分析的处理范围。3）单击【确定】按钮，完成设置。,11,分析结果的范围设置,12,11.3密度分析,密度分析是指根据输入的要素数据集计算整个区域的数据集状况，从而产生一个连续的密度表面。通过计算密度，将每个采样点的值散布到整个研究区域，并获得输出栅格中每个像元的密度值。下图为人口分布的密度表面图。,人口分布的密度表面图,11.3.1核密度分析,核密度分析用于计算要素在其周围邻域中的密度，既可计算点要素的密度也可计算线要素的密度。在核密度分析中，落入搜素区域内的点具有不同的权重，靠近格网搜索中心的点会被赋予较大的权重。核密度分析对话框如左图所示。,13,【核密度分析】对话框,14,11.3.2线密度分析,线密度操作步骤如下：1)在ArcToolbox中双击【SpatialAnalyst工具】-【密度分析】-【线密度分析】，打开【线密度分析】对话框。2）在【线密度分析】对话框中，输入【输入折线要素】【Population字段】数据，指定【输出栅格】的保存路径和名称。3）在【输入像元大小】和【搜索半径】文件框中输入输出栅格数据集的单元大小和密度计算的搜索半径。4）单击【确定】按钮，完成密度图制作。,15,11.3.3点密度分析,点密度分析用于计算每个输出栅格像元周围点要素密度。从理论上讲，每个栅格像元中心的周围都定义了一个邻域，将落入邻域内的样本点的Population字段值相加，然后除以邻域面积，即得到点要素的密度。,16,11.4距离分析,距离分析是指根据每一个栅格相距其最邻近要素的距离分析结果，得到每一个栅格与其临近源的相互关系。通过距离分析，便于人们对资源进行合理的配置和利用。距离分析的两种主要方法是欧氏距离工具盒成本加权距离工具。距离分析的两个重要概念是源和成本。,17,11.4.1源和成本的概念,1.源源是距离分析中的目标或目的地。源是一些离散的点、线、面要素，要素可以相邻，但属性必须不同。源可以是栅格数据，也可以是矢量数据。2.成本成本是指到达目标、目的地的话费，如金钱、时间等。,18,11.4.2欧氏距离,欧氏距离工具根据直线距离描述每个像元与一个源或一组源的关系。欧氏距离工具有三种：欧氏距离、欧氏方向、欧氏分配。下图为某区域的欧氏距离图。,欧氏距离数据,19,11.4.3成本距离,成本距离操作步骤1）在ArcToolbox中双击【SpatialAnalyst工具】-【距离分析】-【成本距离】，打开【成本距离】对话框。2）在【成本距离】对话框中，输入【输入栅格数据或要素源数据】、【输入成本栅格数据】，指定【输出距离栅格数据】的保存路径和名称。3）在【最大距离】文本框中输入积累成本值不能超过的阈值。4）单击【确定】按钮，生成成本距离数据。,20,11.4.4成本路径,成本路径工具用于确定目标点与源点之间的最低成本路径。成本路径的计算过程中，出发地可以是要素点，也可以是区域要素。所以存在三种成本路径的计算方法，即1)EACH-CELL.为每个区域中的没个栅格单元寻找一条成本最低的路径。2）EACH-ZONE.为每个区域寻找一条成本最低路径。3）BEST-SINGLE.在所有的区域中寻找一条成本最低路径。,21,11.4.5最小成本廊道,廊道是不同于两侧基质的狭长地带，具有通道和阻隔的双重作用。如果从廊道栅格中选择值小于最大积累距离的所有像元，则输出栅格为不超过制定成本的一片带状像元。生成的输出阈值即可视为像元的最小成本廊道，下图为廊道分析结果。,廊道分析结果,22,11.4.6成本分配,具体操作步骤如下：1）在ArcToolbox中双击【SpatialAnalyst】-【距离分析】-【成本分配】，打开【成本分配】对话框。2）在【成本分配】对话框中，输入【输入栅格数据或要素源数据】、【源字段】以及【输入成本栅格数据】数据，指定【输出分配栅格数据】的保存路径和名称。3）在【最大距离】文本框中输入积累成本值不能超过的阈值4）单击确定按钮，完成操作。,23,11.5提取分析,提取分析工具用于获取感兴趣的栅格单元，如提坡度大于10%的所有像元。栅格数据的提取分析大致可以分为两种：一类是按属性、形状或者位置提取像元子集；另一类是将像元值提取到点要素，并将这些值记录到点要素类的属性表中。,24,11.5.1按属性、形状或位置提取,1、按属性提取利用【按属性提取】工具可将满足指定属性查询条件的像元提取到心输出栅格中，如提取高程大于1000m的所有像元，或提取土地利用类型中属性为商业用地的所有像元。2、按形状提取3、按位置提取基于空间位置提取像元有两种方法：一种是通过定义一组感兴趣的坐标点，从某栅格中提取规定的像元；另一种是使用眼膜来确定要提取的像元。,25,11.5.2将像元值提取到点要素,1、将像元素提取到新点要素数据集中利用【值提取至点】工具，可以基于一组点要素提取栅格的像元素，并将这些值记录到输出要素类的属性表中，同时为输出要素类添加一个存储新值的名为RASTERVALU的新字段。2、将像元值追加到现有点要素数据集中利用【多值提取至点】工具，可以在点要素类的指定位置提取一个或者多个栅格中的像元值，并将值记录到点要素类的属性表中。,26,11.6.1插值的概念,在区域研究过程中，要获得区域内每个点的数据时非常困难的。一般情况下只采集研究区域内的部分数据，这些数据以离散点的形式存在，只有在采样点上才有准确的数值，未采样点上都没有数值。然而，在实际应用中却经常需要用到某些未采样点的值，此时就需要将已知样本点的值按照一定方法扩散开来，给其他的点赋予一个合理的预测值，这就是插值。,27,11.6.2插值方法,1、反距离权重法2、样条函数法3、克里金法4、自然邻域法,28,11.7重分类,简单的说，重分类就是对原有栅格像元值重新分类从而得到一组新值并输出。所有重分类方法适合区域中的每个像元。也就是说，当对现有值应用某替代值时，所有重分类方法都可以将该替代值应用到原始区域的各个像元。重分类方法不会仅对输入区域的一部分应用替代值。重分类工具包括重分类、查找表、分割、使用表和ASCII文件中分类等。,29,11.7.1重分类,在实际应用中，进行重分类的原因一般有新值替代、将值组合到一起、按相同等级对一组栅格的值进行重分类、将特定值设定为NoData或者为NoData像元设置某个值四种。1、新值替代2、将值组合到一起3、按相同等级对一组栅格的值进行重分类4、将特定值设置为NoData或者为NoData像元设置为某个值,11.7.2查找表,操作步骤：1）在ArcToolbox中双击【SpatialAnalyst工具】-【重分类】-【查找表】，打开【查找表】对话框。2）在【查找表】对话框，输入【输入栅格】和【查找字段】数据，指定【输出栅格】的保存路径和名称。3）单击【确定】按钮，完成操作。下图为查找结果属性表,30,查找结果属性表,31,11.7.3分割,分割工具是按一系列相等的间隔来划分值得整个范围，或者通过将各个像元数量划分到一定数量的组中并保证每组分到的像元数量相等来进行划分。例如，如果输入栅格中值得范围为1-200，而要分割的间隔数为10，则输出栅格的值将介于1-10.输入栅格种植介于1-20的像元将指定为1，值介于21-40的像元将指定为2，以此类推。,32,11.7.4使用表重分类,使用表重分类通过使用重映射表和重分类表将单个值、一定范围内的值、字符串或NoData映射为其他值或NoData。重映射表可以是ASCII文件或INFO表，由两部分组成：第一部分是要重分类的特定像元值，第二部分是像元重分类后的输出值，如INFO表。ASCII重映射表由注释、关键字和赋值语句组成。,33,11.7.5使用ASCII文件重分类,操作步骤如下：1）在ArcToolbox中双击【SpatialAnalyst工具】-【重分类】-【使用ASCII文件重分类】，打开【使用ASCII文件重分类】对话框。2）在【使用ASCII文件重分类】对话框，输入【输入栅格】和【查找字段】数据，指定【输出栅格】的保存路径和名称。3）【将缺失值更改为NoData】为可选项，若选中则栅格像元中未在重映射表中出现或重分类的值被重分类为NoData.4）单击【确定】按钮，完成操作。,34,11.8.1条件分析,条件分析工具包括天剑函数工具、选取函数工具和设为空函数工具。条件函数根据像元值在指定的条件语句中的真假来控制每个像元的输出值。选取函数根据位置栅格数据上每个像元的值来确定输出栅格上的相应位置将使用哪一个输出栅格的值。,35,11.8.2栅格计算器,栅格计算式栅格数据空间分析中进行数据处理和分析最为常用的方法。利用栅格计算器，除了可以方便地完成基于数学运算符、基于数学函数的栅格运算，还支持调用ArcGIS自带的栅格数据空间分析函数。栅格计算器的几种用途：1）简单算术运算2）数学函数运算3）空间分析函数运算,36,11.9太阳辐射分析,太阳辐射是地球上各种物理过程和生物过程的主要能量来源。入射太阳辐射穿过大气层时会发生改变，由于地形和表面要素影响又进一步改变，最后在地球表面被拦截成直射部分、散射部分和反射部分。,37,11.9.1太阳辐射的基本概念,太阳辐射分析中常用的几个基本概念有：视域、太阳图和星空图等。1、视域视域是从某特定位置观看天空时，整个天空可见或遮挡的栅格数据表达。2、太阳图太阳图用来表示太阳轨迹，即太阳随时间变化而产生的位置变化。3、星空图散射辐射是云、粒子等大气成分分散光线的结果。,38,11.9.2太阳辐射的实现,下面对单个点的太阳辐射情况进行分析。操作步骤：（1）在ArcToolbox中双击【SpatialAnalyst工具】【太阳辐射】【太阳辐射点】，打开【太阳辐射点】对话框。（2）在【太阳辐射点】对话框中，填入【输入栅格】、【输入点要素或表】数据，指定【输出总辐射要素】的保存路径和名称。,（3）【高度偏移（可选）】：可选项，输入要执行计算的DEM表面上的高度（以米为单位），高度偏差将应用到所有输入位置。（4）在【纬度（可选）】和【天空大小/分辨率（可选）】文本框中填入位置区域的纬度和天空大小。（5）在【时间配置（可选）】下拉框中选择用于计算太阳辐射的时间配置（时段）。,39,40,11.10表面分析,表面分析是为了获得原始数据中暗含的空间特征信息，如等值线、坡度、坡向、山体阴影等。ArcGIS表面分析的主要功能有：从表面获取坡度和坡向信息、创建等值线、分析表面的可视性、从表面计算山体的阴影、确定坡面线的高度等。,11.10.1坡向,操作步骤：（1）在ArcToolbox中双击【SpatialAnalyst工具】【表面分析】【坡向】选项，打开【坡向】对话框。（2）在【坡向】对话框中，填入【输入栅格】数据，指定【输出栅格】的保存路径和名称。（3）单击【确定】，完成操作。右图为某区域的坡向图。,41,坡向结果图,42,11.10.2等值线,等值线就是将表面上相邻的等值点（诸如高程、温度、降水、大气压力等）连接起来的线。等值线的集合常被称为等值线图，但也可拥有特定的术语称谓。例如表示压力的称为等压线图，表示温度的称为等温线图，表示高程的等高线图是最常使用的等值线图。等值线的分布显示出整个表面上值的变化情况，等值线越密，表面值得变化越大，反之越小。,43,11.10.3填挖方,操作步骤：（1）在ArcToolbox中双击【SpatialAnalyst工具】【表面分析】【填挖方】选项，打开【填挖方】对话框。（2）在【填挖方】对话框中，填入【输入填/挖之前的栅格表面】、【输入填/挖之后的栅格表面】数据，指定【输出栅格】的保存路径和名称。（3）在【Z因子（可选）】文本框中填入Z因子。（4）单击【确定】，完成操作。,44,11.10.4山体阴影,山体阴影通过考虑照明源的角度和阴影，根据表面栅格创建晕渲地貌。它根据假想的照明光源对高程栅格图运用山影函数，计算每个单元以及相关邻域单元的照明值，很好地表达地形的立体形态，而且可以提取地形遮蔽信息。在创建山体阴影图时，主要考虑太阳方位角和太阳高度。,11.10.5坡度,（1）在ArcToolbox中双击【SpatialAnalyst工具】【表面分析】【坡度】，打开【坡度】对话框。（2）在【坡度】对话框中，填入【输入栅格】数据，指定【输出栅格】的保存路径和名称。（3）【输出测量单位（可选）】：可选项，选择坡度的表示方法，有两种情况。（4）在【Z因子（可选）】文本框中填入Z因子。（5）单击【确定】，完成操作。左图为某区域的坡度图。,45,坡向图,46,11.10.6曲率,地面曲率是对地形表面上一点扭曲变化程度的定量化度量因子。曲率计算的输出结果为每个像元的表面曲率，该值通过将该像元与八个相邻像元拟合而得。曲率输出栅格的单位以及可选输出剖面曲线栅格和输出平面曲线栅格的单位是z单位的百分之一（1/100）。某山区（平缓地貌）的全部三个输出栅格的合理期望值介于-0.5至0.5之间；如果山势较为陡峭崎岖（极端地貌），那么期望值介于-4至4之间。,47,11.10.7可见性分析,操作步骤：（1）在ArcToolbox中双击【SpatialAnalyst工具】【表面分析】【视域】，打开【视域】对话框。（2）