Shapefile aus Lat / Long-Werten mit ArcPy erstellen? [geschlossen]

10

Wie kann ich mit Python in ArcGIS 10 ein Shapefile erstellen?

Ich habe lat & long.

Daraus benötige ich Python-Code, der in ArcGIS Desktop 10 ein Shapefile erstellt.

arcgis-10.0 arcpy shapefile Mehul
quelle

Antworten:

15

Zum Erstellen von Punkten:

ptList =[[20.000,43.000],[25.500, 45.085],[26.574, 46.025], [28.131, 48.124]]
pt = arcpy.Point()
ptGeoms = []
for p in ptList:
    pt.X = p[0]
    pt.Y = p[1]
    ptGeoms.append(arcpy.PointGeometry(pt))

arcpy.CopyFeatures_management(ptGeoms, r"C:\Temp\test.shp")

Es wird eine Nachricht wie folgt zurückgegeben:

<Result 'C:\\Temp\\test.shp'>
Aragon
quelle
6

Eine andere Möglichkeit wäre, nur die vorhandenen Arcpy-Geoverarbeitungswerkzeuge zu verwenden (siehe Code unten).

   # Import arcpy module
import arcpy


# Local variables:
table_dbf = "C:\\temp\\table.dbf"
table_Layer2 = "table_Layer2"
point3_shp = "C:\\temp\\point3.shp"

# Process: Make XY Event Layer
arcpy.MakeXYEventLayer_management(table_dbf, "x_coord", "y_coord", table_Layer2, "", "")

# Process: Copy Features
arcpy.CopyFeatures_management(table_Layer2, point3_shp, "", "0", "0", "0")

mxd = arcpy.mapping.MapDocument(r"C:\temp\Untitled.mxd")
df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]
addLayer = arcpy.mapping.Layer(point3_shp)
arcpy.mapping.AddLayer(df, addLayer, "BOTTOM")
mxd.saveACopy(r"C:\temp\Untitled1.mxd")
Artwork21
quelle
6

Sie können ein Shapefile in Python mit dem Werkzeug Feature-Class erstellen erstellen . Am Ende der Seite befindet sich ein Beispiel.

Um das Shapefile mit Ihren Lat & Long-Daten zu füllen, können Sie einen Insert Cursor verwenden .

Vielleicht können Sie Ihre Lat & Long-Daten als Liste in Python laden und dann das Array durchlaufen, das die Zeilen Ihres neuen Shapefiles mit dem Einfügecursor füllt.

Eine Python-Liste von Koordinaten kann folgendermaßen erstellt werden:

latLonList = [[40.000,-75.000],[39.998,-75.432],[39.981,-75.343]]

Gehen Sie dann folgendermaßen vor, um die Koordinaten in der Liste zu durchlaufen (und sie beispielsweise auszudrucken):

for coord in latLonList:
    print "lat: " + str(coord[0])
    print "lon: " + str(coord[1])

Informationen zum Hinzufügen einer Ebene zu einer MXD-Datei finden Sie unter Hinzufügen eines Shapefiles oder einer Feature-Class als Ebene in ArcGIS Desktop mit Python / ArcPy?

Casey
quelle