Criando shapefile a partir de valores lat / long usando o ArcPy? [fechadas]

10

Como posso criar um shapefile usando Python no ArcGIS 10?

Eu tenho lat e long.

A partir disso, preciso do código Python, que criará um shapefile no ArcGIS Desktop 10.

arcgis-10.0 arcpy shapefile Mehul
fonte

Respostas:

15

Para criar pontos:

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

Ele retornará uma mensagem como esta:

<Result 'C:\\Temp\\test.shp'>
Aragão
fonte
6

Outra opção seria usar apenas as ferramentas de geoprocessamento existentes, consulte o código abaixo.

   # 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
fonte
6

Você pode criar um shapefile no Python usando a ferramenta Create Feature Class . Há um exemplo na parte inferior da página.

Para preencher o shapefile com seus dados de latência e longitude, você pode usar um Insert Cursor .

Talvez você possa carregar seus dados de latência e longitude como uma lista no Python, e percorrer a matriz que preenche as linhas do seu novo shapefile com o cursor de inserção.

Uma lista python de coordenadas pode ser construída assim:

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

Em seguida, para percorrer as coordenadas na lista (e imprimi-las, por exemplo), faça o seguinte:

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

Para adicionar uma camada a um arquivo mxd, consulte Adicionando shapefile ou classe de recurso como camada no ArcGIS Desktop usando Python / ArcPy?

Casey
fonte