Conteggio dei poligoni sovrapposti da un Shapefile usando ArcGIS Desktop?

Ho uno shapefile con 16.400 poligoni. Ogni poligono mostra l'espansione di una specie di uccelli per tutto il mondo.

Ora devo contare i poligoni sovrapposti. L'ho provato con l'unione e lo scioglimento (conta l'unione), ma l'unione non funziona per così tanti poligoni.

Poi ho provato a tagliare i continenti, ma anche questo non funziona a causa dell'enorme numero di poligoni. Inoltre ho provato questo metodo , anche senza successo.

Quindi ti sto chiedendo che esiste un modo per contare poligoni sovrapposti se 16400 poligoni sono in un file di forma?

Sto lavorando con 10.0 e posso lavorare con 10.2. Anche una soluzione ArcPy è meravigliosa.

Al momento sto pensando di creare una rete da pesca e iterare sopra le file dello shp con i poligoni 16400 e scrivere 1 su un campo valore di una cella a rete se il poligono si trova in questa cella e quindi prendere la riga successiva (poligono) e se questo è anche nel conteggio delle celle a rete +1.

Ma non so se questa è una buona soluzione e come realizzarla. Oppure devo imparare R per usare questo approccio.

Il risultato: dovrebbe essere una forma in cui si hanno nuovi poligoni tra quelli sovrapposti e un campo in cui vengono conteggiate le sovrapposizioni.

Quindi alla fine dovrebbe esserci uno shapefile in cui puoi vedere quante specie di uccelli si trovano nello stesso posto.

Consiglierei di utilizzare lo strumento Conta sovrapposizione funzionalità (analisi) .

Genera funzioni sovrapposte planarizzate dalle funzioni di input. Il conteggio delle funzioni sovrapposte viene scritto nelle funzioni di output.

Utilizzando token geometrici arcpici, potresti provare qualcosa del genere:

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))

Per i dati di esempio sopra riportati, il codice restituirà i seguenti conteggi di sovrapposizione:

Il codice così com'è restituirà i conteggi solo per i poligoni che hanno almeno una sovrapposizione.

Un metodo molto semplice è:

1. Unione shapefile con se stesso;
2. Converti output multipart in single part ;
3. Utilizzare lo strumento di unione spaziale per contare le sovrapposizioni (utilizzare l'opzione di corrispondenza ARE_IDENTICAL_TO);
4. Simbolizza usando il join_countcampo.

Immagino tu abbia provato questo metodo: il conteggio e la rasterizzazione dei poligoni si sovrappongono in ArcGIS Desktop?

16.400 poligoni non sono così tanti. Tuttavia, una potenziale soluzione è semplicemente fare un Join spaziale regolare. Nel ArcMap toolbox, > Analysis Tools -> Overlap -> Spatial Join.

Impostare entrambe le funzioni targete joinsullo stesso set di dati e specificare un output. Lascia il resto delle impostazioni.

Dopo qualche istante dovresti recuperare uno shapefile che contiene una colonna "conteggio dei join". Sottrai 1 da questo (come ovviamente ogni caratteristica dovrebbe intersecarsi) e quello dovrebbe essere il numero di "sovrapposizioni" (effettivamente interseca) per ogni poligono.

L'ho appena suonato

Ho scaricato e provato lo strumento "Conta poligoni sovrapposti". Potrebbe funzionare, ma richiede molto tempo (probabilmente a causa della dimensione del file, ma il mio input FC aveva solo <5.000 record).

Mentre aspettavo che lo strumento funzionasse, ho aperto un'altra finestra di ArcMap e ci sono voluti solo un paio di passaggi rapidi per ottenere quello che volevo. 1) Join spaziale: utilizza la stessa classe di funzionalità Target e Join e seleziona l'opzione "Unisci da uno a molti". 2) Dissolve: utilizza l'output dell'ultimo passaggio. Utilizza "TARGET_FID" come campo di dissolvenza e per le statistiche puoi SOMMARE il campo "Join_Count" o COUNT il campo "JOIN_FID". 3) Nel file di output dal passaggio 2, utilizzare la calcolatrice di campo sottrai 1 dal campo stats ("SUM_Join_Count" o "COUNT_JOIN_FID") - poiché ogni funzione si interseca da sola.

Suggerisco di usare questo metodo con lo strumento "Conta poligoni sovrapposti". Ho iniziato a eseguire lo strumento COP ~ 5 minuti prima di iniziare questo metodo Join-> Dissolve e mi ha dato il risultato con abbastanza tempo per scriverlo prima ancora che lo strumento "Count Overlapping Polygon" fosse finito.

Spero che sia di aiuto!