GIS-空间分析（10）

提取分析

在 GIS 项目中，我们经常需要裁剪、拼接从 internet 上下载的栅格数据。使其适合于研究区域。裁剪可以通过指定分析掩膜或由 x 和 y 的最小、最大值定义矩形的研究区域范围，然后使用较大的栅格作为输入栅格。镶嵌（mosaic）可以将多个输入栅格拼接成一个栅格。如果输入栅格出现重叠，GIS 软件一般都会提供在重叠区域填充像元值的选项。例如 ArcGIS允许用户选择第一个输入栅格的数据或由所有输入栅格的混合数据作为重叠区的值；如果输入栅格之间存在小缝隙，有种方法是用邻域平均运算来填充缺失的数值。
提取分析是指从一个现有栅格提取数据生成一个新的栅格，提取栅格数据的条件可以是一个数据集、图形对象或查询表达式。如果数据集是点要素图层，提取工具提取该点位置的像元值，并把该像元值赋予到要素属性表的新字段中。如果数据集是一个栅格或一个多边形要素图层，提取工具提取由栅格或者多边形图层定义范围内的像元值，区域以外的像元被赋予 no data。
用于栅格数据提取的图形对象可以是一个矩形、一组点、一个圆或一个多边形。该图形对象以 x 和 y 坐标值输入。例如，一个圆可以用一对 x 和 y 坐标值为其圆心，一个长度值为其半径输入。使用这种方法，可提取距震中 50 米半径范围内的高程数据，或者提取一组气象站的高程数据。
使用“按属性提取”可以生成一个像元值符合查询表达式的新栅格。例如，我们可以创建一个落在特定高程分区的新栅格。在输出栅格中，这一高程分区以外的像元被赋予no data。

任务 1：按属性、形状或位置提取 
提取分析工具可用于根据像元的属性或其空间位置从栅格中提取像元的子集。也可以获取特定位置的像元值作为点要素类中的属性或表。
基于像元的属性或空间位置提取像元值到一个新栅格的工具包括以下几种：
按照属性值提取像元（按属性提取）可通过一个 where 子句来完成。例如，在您的分析中可能需要从高程栅格中提取高程大于 100 米的像元，如图所示。

                                                         

按照像元空间位置的几何提取像元时，要求像元组必须位于指定几何形状的内部或外部。按照指定位置提取像元时，需要根据像元的 x，y 点位置来识别像元的位置（用点提取），或者通过使用掩膜或栅格数据来识别像元的位置（按掩膜提取）。
1）实验数据：“\chp2\2.6 其他的栅格数据操作\data\Stowe.gdb”。
2）实验步骤：
①在 ArcToolbox 中双击【Spatial Analyst 工具】-【提取分析】-【按属性提取】，打开【按属性提取】对话框，如图 所示，其参数说明见表所示。

②在【按属性提取】对话框中，输入【输入栅格】数据，指定【输出栅格】的保存路径和名称。
③单击【Where 子句】右边的 按钮，在弹出的【查询构造器】对话框中输入用于选择栅格像元子集的逻辑表达式，如图所示。双击“VALUE”，点击【>=】按钮，然后单击【获取唯一值】，选择“3146”完成逻辑表达式，单击【确定】按钮返回

                                        

                                                    

④单击【确定】按钮，完成操作，高程提取结果如图 所示，黑色部分为提取结果

         

                                 

实验说明：
如果 Where 子句的求值结果是 True，则将为该像元位置返回初始输入值。如果其求值结果是 False，则将为像元位置指定 NoData。 
任务2：将像元值提取到点要素 
用于指定将像元值提取到属性表或常规表中的像元位置的工具包括以下各项：
通过点要素类识别的像元值可以记录为新输出要素类的属性（值提取至点）。通过点要素类识别的像元值可以追加到要素类的属性表中（多值提取至点）。所识别位置（栅格和要素）的像元值可记录在表中（采样）。
1）实验数据：“\chp2\2.6 其他的栅格数据操作\data\Stowe.gdb”。
2）实验步骤：
①在 ArcToolbox 中双击【Spatial Analyst 工具】-【提取分析】-【值提取至点】，打开【值提取至点】对话框，如图 所示，其参数说明见表所示。

②在【值提取至点】对话框中，输入【输入点要素】和【输入栅格】数据，指定【输出点要素】的保存路径和名称。
③点击【确定】按钮，完成操作

                     

                    

实验说明：
1）输入点要素类中的所有字段均将包括在输出点要素类中。
2）输出要素类将添加一个名为 RASTERVALU 的新字段，该字段不能已存在于输入要素的属性表中。否则，将发生错误。对多波段栅格使用值提取至点时，RASTERVALU 字段将包含输入栅格的最后一个波段的值。要从多个栅格或某个多波段栅格数据集提取值，请使用多值提取至点工具。对于属性表的 RASTERVALU 字段，若赋值栅格中的像元值为 NoData，
则该像元的字段值为 -9999。
3）如果输入栅格为浮点型，则生成的输出点数据集将仅包含输入要素数据的属性和像元值（由插值选项确定）。
4）添加来自输入栅格的属性时，如果输出点要素数据集是 shapefile，则输入栅格中不能存在名称长度超过 10 个字符的字段。如果存在此类字段，则必须在运行该工具之前对其进行重命名