Erstellen eines quadratischen Puffers um ein Punkt-Feature mit ArcGIS for Desktop?

Ich möchte einen quadratischen Puffer aus einem Punkt-Feature erstellen, verstehe aber den Code, der darin enthalten ist, nicht.

Ähnliche Fragen wurden auf der Website forums.esri gestellt, aber das war vor über 10 Jahren und es funktionierte nicht, als ich den Code ausprobierte.

Wie erstelle ich einen quadratischen Puffer aus einem Punkt-Feature?

Probieren Sie diese Schritte mit ArcMap 10 aus:

1. Puffer Ihre Punktmerkmal (ArcToolbox> Analyse Tools> Nähe> Buffer). Stellen Sie sicher, dass im Feld Lineareinheit der richtige Abstand ausgewählt ist.
2. Geben Sie Ihre neu erstellten Puffer in das Werkzeug Feature-Hüllkurve in Polygon ein (Datenverwaltungstools> Features> Feature-Hüllkurve in Polygon). Stellen Sie sicher, dass das Kontrollkästchen "Mehrteilige Features erstellen" aktiviert ist, wenn Sie mehrere Punkte haben.

Für eine Python-Lösung:

Verwenden von SearchCursor und InsertCursor zum Erstellen quadratischer Puffer

Eine mögliche Lösung besteht darin, Ihre "normalen" runden Puffer mit dem Standard-ESRI-Pufferwerkzeug mit einem beliebigen Radius zu erstellen und dann ein Feature-Envelope-To-Polygon für die resultierende Feature-Class von Puffern auszuführen. Dadurch wird ein quadratisches Hüllkurven-Feature um die Ausdehnung jedes Features erstellt. Feature Envelope to Polygon befindet sich unter Data Management> Features. Das Model Builder-Modell würde etwa so aussehen:

Da das Skript, das am Ende von Aarons Code verlinkt ist, nur für quadratische Puffer verwendet werden kann und nicht das neuere arcpy.da-Modul verwendet, habe ich ein Skript geschrieben, mit dem Rechteckpuffer erstellt werden können. Bei einem Datensatz mit 10.000 zufälligen Punkten dauerte es 10 Sekunden:

import os, arcpy

point_FC = arcpy.GetParameterAsText(0)
w = float(arcpy.GetParameterAsText(1))
h = float(arcpy.GetParameterAsText(2))
output_FC = arcpy.GetParameterAsText(3)

def rect(coord, w, h):
        #Given XY coordinates and rectangle dimensions,
        #return a polygon object of a rectangle centered about the point
        x,y = coord
        w *= 0.5
        h *= 0.5
        xmin,xmax = x-w, x+w
        ymin,ymax = y-h, y+h
        poly = ((xmin, ymax), (xmax, ymax), (xmax, ymin), (xmin, ymin))
        return arcpy.Polygon(arcpy.Array(arcpy.Point(*p) for p in poly))

#Create output feature class.
spatref = arcpy.Describe(point_FC).spatialReference
folder, base = os.path.split(output_FC)
arcpy.CreateFeatureclass_management(folder, base, "POLYGON", spatial_reference=spatref)

#Get field object for every field in input except OID and Shape.
fields = [f for f in arcpy.ListFields(point_FC) if f.type not in ("OID", "Geometry")]

for field in fields:
         arcpy.AddField_management(output_FC, field.name, field.type, field.precision,
                                   field.scale, field.length, field.aliasName,
                                   field.isNullable, field.required, field.domain)

#Get field names to be inputted to cursors.
#Need SHAPE@XY token to read point coords and SHAPE@ token to write polygon coords.
fnames = [f.name for f in fields]
fields_in = fnames[::]
fields_out = fnames[::]
fields_in.append("SHAPE@XY")
fields_out.append("SHAPE@")

#Create buffers and write attributes to output FC, if any.
count = int(arcpy.GetCount_management(point_FC)[0])
arcpy.SetProgressor("step", "Buffering...", 0, count, 1)
with arcpy.da.SearchCursor(point_FC, fields_in) as Scursor,  arcpy.da.InsertCursor(output_FC, fields_out) as Icursor:
        for i,row_in in enumerate(Scursor):                
                #"Convert" point to rectangle
                arcpy.SetProgressorPosition(i)
                feature = list(row_in)
                feature[-1] = rect(feature[-1], w, h)                
                Icursor.insertRow(feature)

Angenommen, Sie verwenden ArcObjects (verwenden Sie die Tags, um die Sprache und die API anzugeben, die Sie verwenden), können Sie IEnvelope.Expandeinen quadratischen Puffer aus der Hüllkurve eines Punkts erstellen, wie im folgenden Beispiel dargestellt: Abrufen aller Features aus der Punktsuche im GeoFeatureLayer-Snippet

ESRI.ArcGIS.Geometry.IEnvelope envelope = point.Envelope;
envelope.Expand(searchTolerance, searchTolerance, false);

Als Alternative zu Aarons Antwort können Sie für Benutzer ohne Advanced-Lizenz das Tool Minimum Bounding Geometry verwenden. Schritte unten (modifiziert von Aaron):

1. Puffer Ihre Punktmerkmal (ArcToolbox> Analyse Tools> Nähe> Buffer). Stellen Sie sicher, dass im Feld Lineareinheit der richtige Abstand ausgewählt ist.
2. Geben Sie Ihre neu erstellten Puffer in das Werkzeug Minimale Begrenzungsgeometrie ein (Datenverwaltungswerkzeuge> Funktionen> Minimale Begrenzungsgeometrie). Verwenden Sie 'RECTANGLE_BY_AREA' oder 'RECTANGLE_BY_WIDTH', die anderen Optionen sind nur mit einer erweiterten Lizenz verfügbar.

BEARBEITEN : Mit dieser Option können Sie die Ausrichtung der resultierenden quadratischen Puffer nicht steuern, ohne die Option "ENVELOPE" (für die eine erweiterte Lizenz erforderlich ist) zu verwenden. Durch Aktivieren der Option "Geometrieeigenschaften als Attribute zur Ausgabe hinzufügen (optional)" wird der resultierende Versatz als "MBG_Orientation" in der Ausgabe-Feature-Class aufgezeichnet. Auf diese Weise können Sie die Features bei Bedarf zurück in die Mitte drehen. Weitere Informationen finden Sie unter Drehen von Polygonen nach Wert aus der Attributtabelle mit ArcPy ?. für eine mögliche Lösung dafür.

Diese Site beschreibt, wie man csv mit Hilfe von JavaScript und js2shapefile von geographiclib in quadratische, rechteckige oder kreisförmige Puffer konvertiert.

Sie können nachsehen, ob dies Ihr Problem löst.