GIS-空间分析（6）

局域运算

局域运算是一个像元接一个像元运算，建立栅格数据分析的核心。局域运算由单个或多个输入栅格生成一个新的栅格，新栅格的像元值可以由输入与输出栅格的关系函数计算得到，或通过分类表对其赋值。

1 单一栅格的局域运算

假定以单一栅格为源数据，基于输入栅格的像元值，局域运算通过空间数学函数计算输出栅格的每个像元值，如表所示

例如，浮点型栅格转换成为整型栅格的过程就是一种简单的局域运算。它用取整（integer）函数逐个像元地进行取整运算。把用百分数表示的坡度栅格转换成为用度数表示的坡度栅格也是一种局域运算，表达式为[Slope_d]=57.296*arctan([slope_p]/100)，如图 所示

                                      

任务1：把高程值从米转成千米

将两个栅格的值逐个像元地相除，或将一个栅格的值除以一个常数，如图 所示

                                     

1）实验数据：“chp2\2.2 局域运算\data\emidalat”；
2）实验步骤：
①在 ArcToolbox 中双击【空间分析】-【数学分析】-【除】，打开【除】对话框，如图 所示，其参数说明见表 所示

                                

②输入【输入栅格数据 1】、【输入栅格数据 2】指定输入和输出栅格的保存路径和名称。
③单击【确定】按钮，完成操作。
实验说明： 
1）使用此工具时输入的顺序很重要。
2）当一个数除以零时，输出结果为 NoData。
3）输入的数据类型将决定输出的数据类型：如果两个输入均为整数，则将执行整数除法，并且输出结果为整数。例如，如果 3 除以 2，则输出为 1。如果两个输入中的任何一个为浮点型，则将执行浮点型除法，并且输出结果也为浮点型值。例如，如果 3 除以 2.0，则输出为 1.5。在“地图代数”中，此工具的等价运算符为“/”。

2 重新分类

重新分类是通过分类生成一个新的栅格数据的局域运算方法，也称再编码，或通过查找表的转换。重新分类方法有两种：第一种方法是一对一改变，即输入栅格中的一个像元值在输出栅格中被赋予一个新值。例如，在输出栅格中，将土地利用栅格数据中的灌溉农地像元赋值为 1。第二种方法是在输入栅格中对一系列像元值赋予新值。例如，在人口密度栅格中，将人口密度为 0~25 人/每平方公里的像元对应的输出栅格赋值为 1，等等。整型的栅格数据可以利用以上两种方法中的任意一种进行重新分类，但是，浮点型栅格数据则只能利用第二种方法进行重新分类。进行重新分类的目的。第一，创建简化的栅格数据。例如，栅格数据可以用 1 代表坡度段 0.0~10.0%，用 2 代表坡度段 10.0~20.0%，依次类推，以此来取代一系列连续的坡度值。第二，生成包含唯一类别或数值的新栅格。例如，10.0%~20.0%的坡度范围。第三，生成表示输入栅格像元值排序结果的新栅格。例如，重新分类结果可表示 1~5 的适宜性排序，1 为最不适宜，5 为最适宜。

任务1：重分类 
1）实验数据：“chp2\2.2 局域运算\data\重分类.tif”。
2）实验步骤：
①在 ArcToolbox 中双击【Spatial Analyst 工具】-【重分类】-【重分类】，打开【重分类】对话框，如图 2-12 所示，其参数说明见表 2-4 所示。
②在【重分类】对话框中，输入【输入栅格】和【重分类】字段数据，指定输出栅格的保存路径和名称。
③单击【分类】按钮，弹出【分类】对话框，如图 2-13 所示。在【方法】下拉框中选择分类方法，在【类别】下拉框中选择分组数，单击【确定】按钮完成。
④若要对每个旧值赋予一个唯一的新值，则单击【唯一】按钮。若要添加新的条目，单击【添加条目】按钮，反之则单击【删除条目】按钮。此外，还可以对新值取反，以及设置新值的精度等

⑤【将缺失值更改为 NoData】为可选项，若选中则栅格像元中未在重映射表中出现或
重分类的值被分类为 NoData。
⑥单击【确定】按钮，完成操作。

                                   

                                   

实验说明： 
1）输入栅格必须具有有效的统计数据。
2）如果要对值的范围重新分类，除两个输入范围的边界外，范围不应重叠。在发生重叠的位置，较低输入范围的最大值将包含在取值范围中，而较高输入范围的最小值将不包含在取值范围中。举例来说，如果指定了两个范围，比如将值 1 到 5 重新分类为 100，将 5 到10 重新分类为 200，则小于或等于 5 的输入值将指定给输出值 100，而大于 5 的输入值（如5.1）则会指定给 200。
3）如果使用工具对话框，则可通过保存【按钮】来存储重映射表以供以后使用。使用【加载】按钮可打开先前通过【保存】按钮创建的重映射表。建议只加载先前用重分类工具保存的表。表的格式是特定的，且必须含有字段 FROM、TO、OUT 和 MAPPING。
4）重分类的重映射表经过修改后，如果选择新的输入栅格，则这些值将不会更新。如果重分类不适合新的栅格，则可按如下方式重新初始化默认重分类：
 编辑重分类字段的值或为其选择新值
 重新启动该工具
5）该工具包含一个用来管理小数位数处理方式的精度控件

3 多个栅格的局域运算

多个栅格的局域运算也涉及图层合成，地图叠置或叠加等操作，由于可以用多个栅格图层进行运算，所以局域运算相当于基于矢量的地图叠置操作。
许多局域运算都同时用多个输入栅格，而非仅用单一输入栅格。除了可用于独立栅格的数学公式外，其他的基于输入栅格的像元值或其频率的度量也都可存储于输出栅格。最大值、最小值、值域、总和、平均值、中值和标准差等统计值都是应用于数值型栅格的度量。例如，显示局域运算计算三个输入栅格数据的平均值的例子。如果输入栅格的一个像元为 no data，则输出栅格中该像元也为 no data，如图所示

                                                                    

其他适于数值型或者类别型数据栅格的度量，诸如众数、少数和唯一值数目等统计值。对于每个像元，众数表示输出频率最高的像元值，少数则为输出频率最低的像元值，类别型栅格则输出不同像元值的数目。例如，图 2-15 显示由三个输入栅格获取的众数统计值的输出栅格。

                                                                    

还有一些局域运算不涉及统计或计算。称为合并的一种局域运算，是将一个独特输出值赋予输入值的每种独特组合。假设一个坡度栅格有三种像元值（0%-20%、20%-40%和大于40%），一个坡向栅格有四种像元值（北、东、南和西）。合并运算生成的输出栅格，每种
坡度和坡向的独特组合有一个值，如 1 代表坡度大于 40%和南坡，2 代表坡度为 20%-40%和南坡，依次类推，如图 所示

                                                                     

任务 1：局部分析

局部分析工具包括合并、查找输入列表中满足指定条件的值的出现次数、查找输入列表中满足指定条件的值和查找输入列表中满足指定条件的位置
合并：接受多个栅格并为输入值的各种唯一组合指定一个新值。每个输入的原始像元值将被记录在输出栅格的属性表中。组合中出现NoData 数据时，输出数据为 NoData。附加项将被添加到输出栅格的属性表中，每个输入栅格对应一个附加项，如图所示

                                          

查找输入列表中满足指定条件的值的出现次数：局部分析中有等于频数、大于频数和小于频数三个工具，可以用来确定每个像元从栅格列表得到的输入值等于、大于或小于指定值输入的次数。例如，已知特定的像元位置，考虑赋值像元为 2，计算三个输入栅格的输出结
果，如图 所示。

           

查找输入列表中满足指定条件的值：可以根据输入栅格中哪个值或哪些值满足指定条件来指定像元的输出值。
频数取值工具：可以为每个位置确定输入值中第 n 个出现频率最高的值。例如，如果指定的频数取值为 2，那么对于每个像元，应该输出输入栅格列表中出现频率第二高的值，如图 所示。
注意：①在频数取值的过程中，只要某输入栅格上的某个位置的像元是 NoData，该位置就将接收 NoData 输出；②如果某个特定位置出现的值都相同，则输出该值；③如果某个特定位置没有出现频率第二高的值或出现频率第二高的值不唯一，则输出 NoData 值；④如果某个特定位置出现的值各不相同则输出 NoData 值。等级工具可以为输入栅格集中的值逐个像元地进行等级排定，返回哪个等级将由等级输
入栅格的值来确定，如图 所示。

        

查找输入列表中满足指定条件的位置：可以用来指定含有特定像元值的栅格的位置，包括最高位置和最低位置。
最高位置：逐个像元来确定一组栅格中具有最大值的栅格的位置，如图 所示，输入三个栅格数据，依次返回每个栅格像元最大值所在的栅格数据编号

                                      

最低位置：基于单个像元确定一组栅格数据中具有最小值的栅格的位置，如图所示

                         

像元统计：根据多个栅格数据计算每个像元的统计数据。可用的统计数据有：众数、最大值、均值、中位数、最小值、少数、范围、标准差、总和及变异度，如图所示。

                     

1）实验数据：“chp2\2.2 局域运算\data\emidalat.tif”。
2）实验步骤：
①在 ArcToolbox 中单击【Spatial Analyst 工具】-【局部分析】-【像元统计数据】，打开像元统计对话框，如图 2-23 所示，其参数说明见表 所示；
②在【输入栅格数据或常量值】中选择一个图层，点击加号按钮将其加入数据列表框，或者单击浏览按钮从磁盘中选择要使用的栅格数据；
③在【输出栅格】中为输出结果指定目录及名称；
④在【叠加统计】中选择统计类型：MEAN；
⑤单击【确定】，完成操作