使用ArcGIS Desktop从一个Shapefile中计数重叠的多边形？

我有一个具有16,400个多边形的shapefile。每个多边形表示整个世界鸟类的扩展。

现在，我必须计算重叠的多边形。我用并集进行了尝试并溶解（计算了并集），但是并集不适用于那么多的多边形。

然后，我尝试裁剪大洲，但是由于多边形数量巨大，这也无法正常工作。而且我尝试了这种方法，也没有成功。

因此，我问你，如果一个shapefile中有16400个多边形，有没有一种方法可以计算重叠的多边形？

我正在使用10.0，可以使用10.2。ArcPy解决方案也很棒。

目前，我正在考虑创建一个渔网，并用16400个多边形遍历shp的行，如果多边形位于该单元格中，则将1写入渔网单元的value字段，然后取下一行（多边形），然后如果也是在渔网单元中，则为+1。

但是我不知道这是否是一个好的解决方案以及如何实现它。或者我必须学习R才能使用这种方法。

结果：应该是一个形状，在其中重叠的多边形中有新的多边形，在其中计算重叠的字段。

因此，最后应该有一个shapefile，您可以在其中查看在同一位置发现了多少种鸟类。

我建议使用“ 计数重叠特征（分析）”工具。

从输入要素生成平面化的重叠要素。重叠要素的数量被写入输出要素。

使用arcpy几何标记，您可以尝试如下操作：

import os
import arcpy

arcpy.env.workspace = r"" #path to workspace
arcpy.env.overwriteOutput = 1

polygon_fc = r"" #path to polygon fc

base = [row for row in arcpy.da.SearchCursor(polygon_fc,["OID@","SHAPE@"])]
compare = base

overlaps_stats = {}

for b in base:
    for c in compare:
        if b[1].overlaps(c[1]):
            #print "{0} overlaps {1}".format(b[0],c[0])
            if overlap_stats.has_key(b[0]):
                overlap_stats[b[0]].append(c[0])
            else:
                overlap_stats[b[0]] = [c[0]]

for key,value in overlap_stats.iteritems():
    print "Polygon {0}:  Overlaps: {1}.".format(key,len(value))

对于上面的示例数据，代码将返回以下重叠计数：

照原样的代码将仅返回至少有一个重叠的多边形的计数。

一个非常简单的方法是：

1. 联盟与自身shape文件;
2. 将多部分输出转换为单部分 ;
3. 使用空间连接工具计算重叠（使用ARE_IDENTICAL_TO匹配选项）；
4. 使用该join_count字段进行符号化。

我猜您尝试过这种方法：在ArcGIS Desktop中对多边形重叠进行计数和栅格化吗？

16,400个多边形并不多。但是，一种潜在的解决方案是简单地进行常规的空间连接。在ArcMap toolbox, > Analysis Tools -> Overlap -> Spatial Join。

将target和join功能都设置为相同的数据集并指定输出。保留其余设置。

片刻之后，您应该返回一个包含“连接计数”列的shapefile。从中减去1（很明显，每个要素都应相交），这应该是每个多边形的“重叠”数（实际上是相交）。

我只是在

我下载并尝试了“计数重叠多边形”工具。它可能可以工作，但是要花很长时间（可能是因为文件大小，但是我的输入FC仅具有<5,000条记录）。

当我等待该工具运行时，我打开了另一个ArcMap窗口，只花了几个简单的步骤就可以得到我想要的。1）空间联接-使用与“目标”和“联接特征”相同的要素类，然后选择“一对多联接”选项。2）溶解-使用最后一步的输出。使用“ TARGET_FID”作为溶解字段，对于统计信息，您可以将“ Join_Count”字段相加，也可以将“ JOIN_FID”字段相加。3）在第2步的输出文件中，使用字段计算器从stats字段中减去1（“ SUM_Join_Count”或“ COUNT_JOIN_FID”），因为每个要素都相交。

我建议在“计数重叠多边形”工具上使用此方法。在开始此Join-> Dissolve方法之前，我开始运行COP工具约5分钟，它给了我足够的时间来编写结果，甚至还没有完成“ Count Overlapping Polygon”工具的编写。

希望这可以帮助！