Verschieben / Versetzen von Punktpositionen mit ArcPy oder ModelBuilder?

Ich habe eine Reihe nicht georeferenzierter CAD-Ebenen (siehe diese Frage ) mit Textanmerkungsfunktionen. Ich habe ein Modell zum Konvertieren des Texts in Punkte erstellt, aber nach dem Konvertieren der Anmerkung in eine Punkt-Feature-Class sehe ich, dass die CAD-Textankerpunkte nicht mit der Mitte des CAD-Textes übereinstimmen (wo die Punkte hingehören).

Daher möchte ich ein Feature programmgesteuert (mit ArcPy oder ModelBuilder) relativ zu seiner aktuellen Position (Delta x, y ) [verschieben] , indem ich einen gemessenen X, Y-Wert verwende, den ich bereitstellen werde.

Auf diese Weise könnte ich die GIS-Punkte anstelle des versetzten CAD-Ankerpunkts wieder dorthin verschieben, wo sie hingehören.

Wie kann ich diese Aufgabe erfüllen?

@PolyGeo hat mit SHAPE @ XY IN 10.1 eine hervorragende Antwort gegeben , aber derzeit verwende ich 10.0. Irgendwelche 10.0 Ideen?

Dieser Code sollte dies mit dem SHAPE @ XY-Token tun, das mit arcpy.da.UpdateCursor in ArcGIS 10.1 geliefert wurde.

import arcpy
# Set some variables
fc = r"C:\temp\test.gdb\testFC"
fc2 = r"C:\temp\test.gdb\testFCcopy"
xOffset = 0.001
yOffset = 0.001
# Code to make a copy which will have its coordinates moved (and can be compared with original)
if arcpy.Exists(fc2):
    arcpy.Delete_management(fc2)
arcpy.Copy_management(fc,fc2)
# Perform the move
with arcpy.da.UpdateCursor(fc2, ["SHAPE@XY"]) as cursor:
    for row in cursor:
        cursor.updateRow([[row[0][0] + xOffset,row[0][1] + yOffset]])

Das hier verwendete Codierungsmuster stammt von ArcPy Café .

Ich verdanke @ artist21, dass er mich zu meiner endgültigen Lösung geführt hat. Im Online-Hilfeartikel von ArcGIS 10.0 mit dem Titel " Feldbeispiele berechnen ", der unter den Unterkategorien " Codebeispiele - Geometrie " und " Verschieben Sie für eine Punkt-Feature-Class die x-Koordinate jedes Punkts um 100 " aufgeführt ist, habe ich tatsächlich ein fast vollständiges Skript gefunden.

Das letzte Skript, das ich im ModelBuilder-Tool "Feld berechnen" verwendet habe, war:

Ausdruck:

shiftXYCoordinates(!SHAPE!,%ShiftX%,%ShiftY%)

Dabei sind ShiftX und ShiftY Variablen (als Parameter), die auf der ModelBuilder-Zeichenfläche definiert sind.

Ausdruckstyp:

PYTHON_9.3

Codeblock:

def shiftXYCoordinates(shape,x_shift,y_shift):
   point = shape.getPart(0)
   point.X += float(x_shift)
   point.Y += float(y_shift)
   return point

Da alle Modelle an einem ausgewählten Satz arbeiten, sollten Sie dies auch als generisches Tool erstellen können, das in Verbindung mit anderen Modellen / Tools in anderen Modelbuilder-Sitzungen funktioniert. Das sehr einfache Modell, das ich erstellt habe (als "Plugin" für andere Modelle zum Verschieben von Koordinatenwerten), sieht folgendermaßen aus. Auf diese Weise kann ich die Verschiebung pro Auswahlsatz steuern (wie in anderen Modellen definiert):

Es hat wie ein Zauber funktioniert, ich danke Ihnen allen für Ihre Beiträge!

Sie können dieses Feldrechnerskript auch verwenden , um Feature-Positionen zu verschieben:

def XYsetVALUE( shape, X_value, Y_value): 
  myMoveX = 0.001
  myMoveY = 0.001
  point = shape.getPart(0) 
  point.X = X_value + myMoveX
  point.Y = Y_value + myMoveY
  return point 

XYsetVALUE (! SHAPE!,! X_COORD!,! Y_COORD!)

Mit der obigen Funktion können Sie eine zusätzliche Methode zum Berechnen von Feldern in Ihr Modell aufnehmen.

Ich habe die Lösung so angepasst, dass Punkte in eine bestimmte Richtung (Winkel) und eine bestimmte Entfernung verschoben / verschoben werden.

Sieht aus wie:

def shiftXYCoordinates(shape,angle,distance):
point = shape.getPart(0)
point.Y += distance * math.cos(math.radians(angle))
point.X += distance * math.sin(math.radians(angle))
return point

und werden wie shiftXYCoordinates (! SHAPE!,! Angle !, 5000) aufgerufen, wenn Sie ein Feld "Winkel" für Ihre Punktmerkmale haben (oder natürlich mit einer Konstanten). Der Winkel sollte in Dezimalgraden angegeben werden. 0 verschiebt sich nach oben, 90 nach rechts usw. Ich habe sie erhalten, nachdem ich Strip-Map-Index-Features erstellt und diese in Punkte konvertiert habe.

Stellen Sie außerdem sicher, dass Sie vor dem Ausführen den Feldnamen „Form“ auswählen :)

(Lösung in ArcMap 10.0 SP5 getestet)

Wie Sie sehen können, ist es in 10.1 viel einfacher, wenn Sie Zugriff auf Cursortoken erhalten.

import arcpy
# Code to move features in copy of same dataset
fc = r"C:\temp\test.gdb\testFC"
fc2 = r"C:\temp\test.gdb\testFCcopy"
xOffset = 0.001
yOffset = 0.001
if arcpy.Exists(fc2):
    arcpy.Delete_management(fc2)
arcpy.Copy_management(fc, fc2)

shape = arcpy.Describe(fc2).ShapeFieldName

cursor = arcpy.UpdateCursor(fc2)
for row in cursor:    
    point = row.getValue(shape).getPart()
    row.setValue(shape, arcpy.Point(point.X + xOffset, point.Y + yOffset))
    cursor.updateRow(row) 

del point, row, cursor

Dies funktioniert für 10.0:

# Featureclass here
FC = r'featureclass'

fcount = 0
shapefield = arcpy.Describe(FC).shapeFieldName
featureUpdate = arcpy.UpdateCursor(FC)
for f in featureUpdate:
    # Hard coded shifts but easy enough to set up a lookup function if needed
    sLon = 0.001
    sLat = 0.001
    # Optional but I like to count to see where it is at in the process
    if fcount % 1000 == 0:
        print('Updating feature %s...' %(fcount))
    # Get the original value
    cF = f.getValue(shapefield)
    cPNT = cF.getPart()
    # Create a new point with the shifted value
    sPNT = arcpy.Point(cPNT.X - sLon, cPNT.Y - sLAT)
    # Set the shapefield to the new point and update feature
    f.setValue(shapefield, sPNT)
    featureUpdate.updateRow(f)
    fcount += 1
del featureUpdate