Anzahl der Eckpunkte in jedem Polygon extrahieren?

Ich habe ArcGIS Desktop 10.2 und meine Herausforderung besteht darin, die Anzahl der Eckpunkte in jedem Polygon für alle Funktionen wie diese zu extrahieren:

Ich habe viele Flurstücke in meiner Feature-Class und möchte die Anzahl der Scheitelpunkte für alle Features separat extrahieren. Anschließend möchte ich die XY-Koordinaten für alle Scheitelpunkte anzeigen.

Für weitere Informationen möchte ich nur Scheitelpunkte für jedes Polygon konvertieren und die Nummer jedes Scheitelpunkts anzeigen. Wenn ich also ein Polygon mit 4 Scheitelpunkten habe, möchte ich ein Polygon in Scheitelpunkte konvertieren und die Anzahl dieser Scheitelpunkte anzeigen (1,2,3,4,5), dann Anzeige von xy für jeden Scheitelpunkt, ich denke, die eigentliche Herausforderung ID, wie alle Polygone in Scheitelpunkte konvertiert und jede Anzahl von Scheitelpunkten ab 1 Nummer.

Der folgende Code kombiniert die anderen Antworten und fügt ein Bit hinzu, um die Eckpunkte zu nummerieren.

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)

Die Punkte sind in Zeichnungsreihenfolge beschriftet. Der letzte Punkt (unter dem ersten) hat keine Beschriftung und kann gelöscht werden, indem alle Punkte mit Nullwerten oder eindeutigen "DRAW_ORDER" -Werten ausgewählt werden, wenn sie nicht für die Rekonstruktion benötigt werden. Mit einer Definitionsabfrage können die überlappenden Punkte aus der Anzeige entfernt werden.

XY-Daten sind vorhanden, aber ich überlasse dies Ihren Beschriftungs- / Anzeigewünschen. Siehe Aarons Antwort zum Hinzufügen eines XY-Felds zur Beschriftung.

Ich spielte auch mit FeatureClass, um Array zu numpieren, aber ich beendete dies zuerst.

Der einfachste Weg, dies zu tun, besteht darin, der Attributtabelle des Parzellen-Layers ein neues Ganzzahlfeld hinzuzufügen. Führen Sie dann den Feldrechner mit dem folgenden Ausdruck aus:

!Shape!.pointCount-!Shape!.partCount

Das !Shape!.pointCountgibt die Gesamtanzahl der Scheitelpunkte in dem Feature zurück. Der erste Scheitelpunkt jedes Teils wird jedoch am Ende wiederholt, um das Feature zu schließen. Subtrahieren Sie dazu für jedes Teil einen Scheitelpunkt mit -!Shape!.partCount.

Beachten Sie, dass Sie den Python-Parser verwenden müssen, damit dieser Ausdruck funktioniert.

dmahr bot eine gute Lösung zum Zählen von Eckpunkten. Verwenden Sie den folgenden Arbeitsablauf, um jeden Punkt nicht programmgesteuert mit den XY-Koordinaten zu kennzeichnen:

1. Feature-Stützpunkte zu Punkten
2. Füge zwei neue Felder (Typ: double) in den neuen Punkt FC "X", "Y" ein
3. Geometrie berechnen. Rechtsklickfeld> Geometrie berechnen ...> X-Koordinate des Punktes (Wiederholung für Y-Feld)
4. Füge ein weiteres Feld "XY" hinzu (Typ: Text)

5. Berechnen Sie das Feld "XY" im Feldrechner, wobei XY =

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

6. Beschriften von Features. Klicken Sie mit der rechten Maustaste auf Ebene> Beschriftungen> Beschriftungsfeld: XY

Dies führt zu folgenden Ergebnissen:

Sie können diese Aktionen auch programmgesteuert explode_to_pointsmit einem Suchcursor ausführen (als Start).

Dekonstruieren Sie ein Feature in seine einzelnen Punkte oder Scheitelpunkte. Wenn explode_to_points auf True gesetzt ist, wird ein Mehrpunkt-Feature mit beispielsweise fünf Punkten durch fünf Zeilen dargestellt.

(Der Standardwert ist False.)

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

Wenn man kein neues Feld berechnen und nur eine Anzahl von Eckpunkten pro Ebene erhalten möchte sehr schnell zurückerhalten (zum Zwecke der Verallgemeinerung, z. B. wenn Sie die Datensätze im Web verfügbar machen), können Sie ein benutzerdefiniertes Skript-Tool in a erstellen Toolbox oder machen Sie den Code als Python-Add-In verfügbar.

Benutzerdefinierter Skript-Tool-Code:

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("***************************************")

Python-Add-In-Code (wählen Sie eine Ebene im Inhaltsverzeichnis aus, um die Scheitelpunkte zu zählen):

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)