Estrarre il numero di vertici in ciascun poligono?

Ho ArcGIS Desktop 10.2 e la mia sfida è come estrarre il numero di vertici in ciascun poligono per tutte le funzionalità come questa:

Ho molti pacchi nella mia classe di caratteristiche e voglio estrarre il numero di vertici per tutte le caratteristiche separatamente, quindi voglio visualizzare le coordinate XY per tutti i vertici.

per ulteriori informazioni, voglio solo convertire vertici per ogni poligono e visualizzare il numero di ogni vertice a partire da 1 numero, quindi se ho un poligono e ha 4 vertici, voglio convertire un poligono in vertici e visualizzo il numero di vertici come questo (1,2,3,4,5), quindi visualizzando xy per ogni vertice, penso che il vero ID sfida in, come convertire tutti i poligoni in vertici e creare ogni numero di vertici a partire da 1 numero.

Il codice seguente combina le altre risposte e aggiunge un po 'per numerare i vertici.

import arcpy
arcpy.env.workspace = "in_memory"
#paths
fc = r"...\polygons"
fc_out = r"...\vertices"
arcpy.MakeFeatureLayer_management(fc, "lyr")
# add fields if needed
for FIELD in ["DRAW_ORDER", "COUNT"]:
    if FIELD not in [field.name for field in arcpy.ListFields(fc)]:
        try:
            arcpy.AddField_management("lyr", FIELD, "SHORT")
        except Exception as e:
            print e
# get the number of points minus overlapping (@dmahr - GSE)
arcpy.CalculateField_management("lyr", "COUNT", "!Shape!.pointCount-!Shape!.partCount", "PYTHON")
# dict to iterate and check count
OIDS = {}
for row in arcpy.da.SearchCursor("lyr", ["OBJECTID", "COUNT"]):
    OIDS[row[0]] = row[1]
del row
# get vertices as points and add XY (@Aaron - GSE)
arcpy.FeatureVerticesToPoints_management("lyr", fc_out)
arcpy.AddXY_management(fc_out)
# start adding a number to the points
for OID in OIDS:
    order_count = 1
    rows = arcpy.da.UpdateCursor(fc_out, ["DRAW_ORDER", "COUNT"], "ORIG_FID = %d"%OID)
    for row in rows:
        # will leave the overlapping as NULL
        if order_count <= OIDS[OID]:
            row[0] = order_count
            rows.updateRow(row)
            order_count += 1
##        # this can set the overlapping to 0 or some unique value (999)
##        else:
##            row[0] = 0
##            rows.updateRow(row)

I punti sono etichettati in ordine di disegno. L'ultimo punto (sotto il primo) non avrà etichetta e può essere cancellato selezionando tutti i punti che hanno valori Null o unici "DRAW_ORDER" se non necessari per la ricostruzione. È possibile utilizzare una query di definizione per rimuovere i punti sovrapposti dalla visualizzazione.

Sono presenti dati XY, ma li lascerò ai tuoi desideri di etichettatura / visualizzazione. Vedi la risposta di Aaron sull'aggiunta di un campo XY per l'etichettatura.

Stavo anche giocando con FeatureClass all'array numpy, ma ho finito prima.

Il modo più semplice per farlo è aggiungere un nuovo campo intero alla tabella degli attributi del livello pacchi. Quindi, eseguire il calcolatore di campo con la seguente espressione:

!Shape!.pointCount-!Shape!.partCount

Il !Shape!.pointCountrestituisce il numero totale di vertici nella funzione. Tuttavia, il primo vertice di ciascuna parte viene ripetuto alla fine, al fine di chiudere la funzione. Per gestire ciò, sottrarre un vertice per ogni parte usando -!Shape!.partCount.

Nota che dovrai usare il parser Python per far funzionare questa espressione.

dmahr ha fornito una buona soluzione per contare i vertici. Per un modo non programmatico di etichettare ciascun punto con i cordoni XY, provare il flusso di lavoro seguente:

1. Caratteristica Vertici a punti
2. Aggiungi due nuovi campi (tipo: doppio) nel nuovo punto FC "X", "Y"
3. Calcola la geometria. Campo clic destro> Calcola geometria ...> Coordinata X del punto (ripetere per il campo Y)
4. Aggiungi un altro campo "XY" (tipo: testo)

5. Calcola il campo "XY" nel calcolatore di campo, dove XY =

   str (! x!) + "," + str (! y!)

6. Funzionalità dell'etichetta. Livello clic destro> Etichette> Campo etichetta: XY

Ciò produce i seguenti risultati:

È inoltre possibile eseguire queste azioni a livello di codice utilizzando explode_to_pointsun cursore di ricerca (come inizio).

Decostruisci una caratteristica nei suoi singoli punti o vertici. Se explode_to_points è impostato su True, una funzione multipunto con cinque punti, ad esempio, è rappresentata da cinque righe.

(Il valore predefinito è falso)

arcpy.da.SearchCursor (in_table, field_names, {where_clause}, {spatial_reference}, {explode_to_points}, {sql_clause})

Se non si desidera calcolare un nuovo campo e si desidera recuperare rapidamente un certo numero di vertici per livello (per scopi di generalizzazione, ad esempio quando si espongono i set di dati sul Web), è possibile creare uno strumento di script personalizzato all'interno di un toolbox o esporre il codice come componente aggiuntivo Python.

Codice dello strumento di script personalizzato:

import arcpy
in_fc = arcpy.GetParameterAsText(0)

features = [feature[0] for feature in arcpy.da.SearchCursor(in_fc,"SHAPE@")]
count_vertices = sum([f.pointCount-f.partCount for f in features])
arcpy.AddMessage("***************************************")
arcpy.AddMessage("Number of vertices in the layer: {0}".format(count_vertices))
arcpy.AddMessage("***************************************")

Codice componente aggiuntivo Python (selezionare un livello nel sommario per contare i vertici):

import arcpy
import pythonaddins

arcpy.env.overwriteOutput = True
mxd = arcpy.mapping.MapDocument("current")
in_fc = pythonaddins.GetSelectedTOCLayerOrDataFrame()

features = [feature[0] for feature in arcpy.da.SearchCursor(in_fc,"SHAPE@")]
count_vertices = sum([f.pointCount-f.partCount for f in features])

pythonaddins.MessageBox(count_vertices, 'Number of vertices in {0}'.format(in_fc.name), 0)