Wie teile ich Punkte mit ArcGIS Desktop nach Polygon?

Ich suche nach einer Möglichkeit, Punkte durch eine County-Polygon-Feature-Class zu teilen. Das angehängte Bild zeigt zwei Feature-Classes: eine aus Punkten und die andere eine County-Polygon-Feature-Class. Das Teilungswerkzeug wäre ideal, wenn es Punkte zulässt. Ich würde eine saubere Ausgabe ohne zusätzliche Felder bevorzugen, die den Punktattributen hinzugefügt werden. Außerdem möchte ich den County-Namen aus der Polygon-Feature-Class haben, die die neuen Punkt-Feature-Classes definiert (z. B. Kiowa, Clark, Comanche). Ich freue mich über Lösungen und Ratschläge.

In diesem Beispiel sollte das Endprodukt aus drei Punktmerkmalsklassen mit den Namen "Kiowa", "Clark" und "Comanche" bestehen, die aus einer größeren Punktmerkmalsklasse hergestellt wurden. Eine automatisierte Lösung wäre ideal, da ich viele, viele FCs mit zusammengeführten Punkten in Dutzenden von Landkreisen habe, mit denen ich arbeiten kann (ungefähr über den Bundesstaat KS).

Versuchen Sie, nachdem Sie Spatial Join oder Intersect verwendet haben, um das County-Namensattribut für jeden Punkt abzurufen , das in der Geoverarbeitungsmodell- und Skript-Tool-Galerie verfügbare Werkzeug " Ebene nach Attributen teilen" von Dan Patterson zu verwenden .

Alternativ können Sie ModelBuilder verwenden, um dies mithilfe eines anderen Ansatzes zu automatisieren, bei dem Ebene nach Standort auswählen ausgewählt wird (klicken Sie auf das Miniaturbild, um das vollständige Bild anzuzeigen).

Aus diesem Thread in den ESRI-Foren entfernt: Stapelauswahl nach Speicherort und Exportieren von Shapefiles

Viel einfacher mit dem Befehl intersect.
Schalten Sie unerwünschte Felder aus.
1. Schneiden Sie mit einer neuen Ausgabe.
2. Verbinden Sie die Ausgabe mit der ursprünglichen OID.
3. Berechnen Sie ein neues Feld mit dem Namen cnty.

In ettools finden Sie einige erweiterte Versionen der räumlichen Verknüpfung und Aufteilung nach Standort.

Nach dem, was ich gelesen habe, würde ich einfach eine Überlagerung (ArcToolbox-Analysetools - Überlagerung - und entweder Identität, Schnittmenge oder räumliche Verknüpfung) Ihrer Punkte mit einem County-Polygon-Dataset ausführen. In ArcMap können Sie die in der County-Ebene angezeigten Felder auf das Feld mit dem County-Namen beschränken, und nur dieses Feld wird dem Ergebnis hinzugefügt. Wenn selbst dieses hinzugefügte Feld zu groß ist, berechnen Sie seine Werte in das übereinstimmende Feld in der Punktetabelle und löschen Sie das hinzugefügte Feld.

Mit Python können Sie eine "Ebene nach Ort auswählen" schleifen, um die Punkte zu erfassen, die die einzelnen Counties überlappen, und sie dann als Shapefile speichern. Es wäre ungefähr so:

import arcgisscripting
# starts geoprocessing
gp = arcgisscripting.create()
gp.OverWriteOutput = 1

# Variable iniciation
points = u"Path for your point shape"
counties = u"Path for counties shape"
outDir = u"path for output directory"

#Load points as a layer
gp.MakeFeatureLayer(points, "points")

# Go county by county
rows = gp.searchcursor(counties)
row = rows.next()

# loop County by county overlap
while row.countyName: #adapt to your countie table of attributes
    # Make a layer from the feature class
    gp.MakeFeatureLayer(counties, "counties lyr", "[countieName]='" + row.countieName + "'")
    # Select all points that intersect the current countie polygon
    gp.SelectLayerByLocation(points, "intersect", "counties_lyr", 0, "new_selection")
    outSHP = outDir + row.counties + u".shp"
    gp.CopyFeatures_management(points, outSHP)
    print(outSHP)
    row = rows.next()
 # End of loop

 del rows, row, gp

 #END

Ich habe es nicht getestet. Sie müssen es wahrscheinlich an Ihren Fall anpassen.

Bearbeiten: Select_analysis von CopyFeatures_management gemäß den Empfehlungen von @ blah238 geändert.

Es gibt jetzt ein eLearning-Video-Tutorial (nicht kostenlos, aber relativ günstig) mit dem Titel Verwenden von ArcPy-Cursorn zum Teilen nach Attributen und Speicherorten, in dem Folgendes gezeigt wird:

Schreiben Sie ein Python-Skript, um eine Feature-Class mithilfe einer Polygon-Feature-Class in eine Feature-Class pro Feature aufzuteilen, das sich an den einzelnen Polygonpositionen befindet

Ich denke, es geht um die Anforderungen dieser Frage.

Haftungsausschluss : Ich bin der Autor und Moderator dieses eLearning-Video-Tutorials