ArcGIS Desktopを使用して、1つのシェープファイルから重複するポリゴンをカウントしていますか？

16,400ポリゴンのシェープファイルがあります。各ポリゴンは、全世界の鳥種の拡大を示しています。

次に、重なり合うポリゴンを数える必要があります。ユニオンとディゾルブを試してみました（ユニオンを数える）が、ユニオンがあまりにも多くのポリゴンで機能しない。

その後、大陸を切り取ろうとしましたが、ポリゴンの数が非常に多いため、これも機能しません。さらに、私もこの方法を試しましたが、成功しませんでした。

したがって、16400のポリゴンが1つのシェープファイルにある場合、重複するポリゴンをカウントする方法があるかどうかを尋ねています。

私は10.0を使用しており、10.2を使用できます。ArcPyソリューションも素晴らしいです。

現在、フィッシュネットを作成し、16400ポリゴンでshpの行を反復処理し、ポリゴンがこのセルにあり、次の行（ポリゴン）を取得する場合、フィッシュネットセルの値フィールドに1を書き込みます。これもフィッシュネットのセル数に含まれる場合+1。

しかし、これが良い解決策であるかどうか、そしてそれを実現する方法はわかりません。または、このアプローチを使用するにはRを学ぶ必要があります。

結果：オーバーラップしているポリゴンから新しいポリゴンができたシェイプと、オーバーラップがカウントされるフィールドになります。

したがって、最終的には、同じ場所にいくつの鳥種が見つかったかを確認できるシェープファイルがあるはずです。

カスタムの[ ポリゴンを数える]ツールを使用することをお勧めします。

説明：

このサンプルには、1つのツール、[重複ポリゴンを数える]を備えたツールボックスが含まれています。重複するポリゴンを含むフィーチャクラスまたはレイヤーを指定して、重複を削除した新しいフィーチャクラスと、重複するポリゴンの数を含むJoin_Countフィールドを出力します。

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つのオーバーラップがあるポリゴンのカウントのみを返します。

非常に簡単な方法は次のとおりです。

1. 組合自体にシェープファイル。
2. マルチパート出力を単一パートに変換します。
3. 使用加入空間的（ARE_IDENTICAL_TOマッチオプションを使用）重複をカウントするためのツールを。
4. join_countフィールドを使用してシンボル表示します。

あなたはこの方法を試したと思います：ArcGIS Desktopでのポリゴンの重なりのカウントとラスタライズ？

16,400ポリゴンはそれほど多くありません。ただし、1つの潜在的な解決策は、通常の空間結合を単純に行うことです。でArcMap toolbox, > Analysis Tools -> Overlap -> Spatial Join。

targetとjoin機能の両方を同じデータセットに設定し、出力を指定します。残りの設定はそのままにします。

しばらくすると、「結合カウント」列を含むシェープファイルが返されます。これから1を引きます（明らかに各フィーチャがそれ自体と交差する必要があるため）。これは、各ポリゴンの「オーバーラップ」（実際には交差する）の数になります。

私はそれを上演しました

「重複ポリゴンのカウント」ツールをダウンロードして試してみました。うまくいくかもしれませんが、非常に長い時間がかかります（おそらくファイルサイズのためですが、私の入力FCには5,000レコードしかありませんでした）。

そのツールが実行されるのを待っている間、別のArcMapウィンドウを開いて、ほんの数ステップで必要なものを取得しました。1）空間結合-ターゲットおよび結合フィーチャーと同じフィーチャークラスを使用し、「1対多に結合」オプションを選択します。2）ディゾルブ-最後のステップの出力を使用します。「TARGET_FID」をディゾルブフィールドとして使用し、統計の場合は、「Join_Count」フィールドを合計するか、「JOIN_FID」フィールドをカウントできます。3）ステップ2からの出力ファイルでは、統計情報フィールド（「SUM_Join_Count」、または「COUNT_JOIN_FID」）からフィールド計算機減算1を使用します-各フィーチャがそれ自体と交差するため。

この方法は、[ポリゴンの数を数える]ツールではなく使用することをお勧めします。このJoin-> Dissolveメソッドを開始する約5分前にCOPツールの実行を開始しました。「ポリゴンのカウント」ツールが終了する前にこれを書き込むのに十分な時間のある結果が得られました。

お役に立てれば！