ArcGIS计算矢量要素集中每一个面的遥感影像平均值、最大值等统计指标

  本文介绍在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