本文介绍在ArcMap 软件中,基于矢量面要素集 ,计算在其中每一个面区域 内,遥感影像的像元个数、平均值、总和等统计值 ,并将统计信息附加 到矢量图层的属性表中的方法。
首先,明确一下本文的需求。现在有一个矢量面要素集 ,其中包含了多个分散的面要素,如下图所示。
其中,其属性表中包含了每一个要素的名称,以及其他相关信息,如下图所示。
此外,我们还有同一地区的一景遥感影像,在本文中就以当地的NDVI数据为例;如下图所示。
我们希望实现的就是,为每一个面要素 ,计算其区域 内NDVI的平均值、累计值(求和),并统计区域内像元的数量。
知道了需求,即可开始操作。其中,我们依次选择"Spatial Analyst Tools.tbx "→"Zonal "→"Zonal Statistics as Table"工具,如下图所示。
随后,在弹出的窗口中,配置相应的参数。其中,第一个参数就是需要进行统计的矢量要素集 ;第二个参数则是统计这些统计值时,对矢量面加以区分的字段 的名称------我们可以选择其属性表中,任何一个不含有相同值 的字段(或者直接选择FID
这个字段就行);第三个参数是需要被统计的遥感影像;第四个参数则是要导出的结果表的路径与名称。此外,还可以在最后一个参数中选择我们所需要的统计指标。
执行上述工具,将得到一个结果表格。在表格上右键,并选择"Open"选项,如下图所示。
可以看到,其按照我们所选的字段,分别对每一个矢量面中,遥感影像的不同统计值加以了计算,如下图所示。
但是,因为我这里是基于FID
这个字段来区分不同的面的,但是光看一个FID
我也不知道这个面是哪一个面,所以希望可以将带有区域名称的矢量数据的属性表(也就是本文开头的第二张图片),与我们上述工具计算得到的结果表匹配起来。
这个操作也很简单。在矢量要素集上右键,选择"Joins and Relates "→"Join"选项,如下图所示。
在其中,选择匹配所依据的关键字,我这里就选择FID
这个字段了;同时,选择好我们得到的结果表。
接下来,即可开始匹配。如果匹配过程中,出现如下图所示的提示,直接选择"Yes"即可。
完成匹配后,可以看到我们此时的表格中,就既有每一个区域的名称,也有其对应的统计数值了;如下图所示。
但是,此时这个匹配只是暂时的,我们还需要将匹配后的结果图层导出,才可以保持匹配结果。
依然在矢量要素集上右键,选择"Data "→"Export Data"选项,如下图所示。
注意选择导出的范围为"All features",并配置好导出图层的路径与文件名称,如下图所示。
接下来,即可得到属性表带有匹配后结果的图层了,如下图所示。
至此,大功告成。
欢迎关注:疯狂学习GIS