Escribir geometrías

Al utilizar los cursores insertar y actualizar, los scripts pueden crear nuevas entidades en una clase de entidad o actualizar las existentes. Un script puede definir una entidad al crear un objeto de Point, completar sus propiedades y colocarlo en un Array. Después, ese conjunto puede utilizarse para establecer la geometría de una entidad utilizando las clases de geometría Polygon, Polyline, PointGeometry o Multipoint.

import arcpy
fc = "c:/data/gdb.gdb/roads"
cursor = arcpy.da.InsertCursor(fc, ["SHAPE@"])
array = arcpy.Array([arcpy.Point(-77.4349451, 37.5408265),
                     arcpy.Point(-78.6384349, 35.7780943)])
spatial_reference = arcpy.SpatialReference(4326)
polyline = arcpy.Polyline(array, spatial_reference)
cursor.insertRow([polyline])

Según se mostró anteriormente, una sola parte de la geometría se define mediante un conjunto de puntos. De igual manera, se puede crear una entidad multiparte a partir de un conjunto de conjuntos de puntos según se muestra más adelante utilizando el mismo cursor.

first_part = arcpy.Array([arcpy.Point(-77.4349451, 37.5408265),                          arcpy.Point(-78.6384349, 35.7780943)]) second_part = arcpy.Array([arcpy.Point(-79.7910143, 36.0786785),                           arcpy.Point(-80.8546435, 35.2315402)])
array = arcpy.Array([first_part, second_part]) spatial_reference = arcpy.SpatialReference(4326) multipart_feature = arcpy.Polyline(array, spatial_reference)
cursor.insertRow([multipart_feature])

Al escribir entidades de punto, se utiliza solamente un solo objeto de punto para establecer la geometría de una entidad de puntos. También se pueden crear puntos más fácilmente (y eficientemente) utilizando el token SHAPE@XY (y los tokens SHAPE@M y SHAPE@Z, según sea necesario).

import arcpy
# fc is a point feature class
fc = "c:/data/gdb.gdb/stops"
cursor = arcpy.da.InsertCursor(fc, ["SHAPE@XY"])
xy = (5997594.4753, 2069901.75682)
cursor.insertRow([xy])

Todas las geometrías se validan antes de que se escriban en una clase de entidad. Problemas como una incorrecta orientación del anillo y polígonos que se intersecan a sí mismos, entre otros, se corrigen cuando se simplifica la geometría antes de su inserción.

El siguiente ejemplo muestra cómo leer un conjunto de coordenadas (definido por coords_list) que contiene una serie de coordenadas lineales y utilizarlo para crear una nueva clase de entidad.

# Create a new line feature class using a text file of coordinates.
#   Each coordinate entry is semicolon delimited in the format of ID;X;Y
import arcpy
import os
# List of coordinates (ID, X, Y)
coords_list = [[1, -61845879.0968, 45047635.4861], 
               [1, -3976119.96791, 46073695.0451],
               [1, 1154177.8272, -25134838.3511],
               [1, -62051091.0086, -26160897.9101],
               [2, 17365918.8598, 44431999.7507],
               [2, 39939229.1582, 45252847.3979],
               [2, 41170500.6291, 27194199.1591],
               [2, 17981554.5952, 27809834.8945],
               [3, 15519011.6535, 11598093.8619],
               [3, 52046731.9547, 13034577.2446],
               [3, 52867579.6019, -16105514.2317],
               [3, 17160706.948, -16515938.0553]]
# The output feature class to be created
outFC = arcpy.GetParameterAsText(0)
# Get the template feature class
template = arcpy.GetParameterAsText(1)
cur = None
try:
    # Create the output feature class
    arcpy.CreateFeatureclass_management(os.path.dirname(outFC),
                                        os.path.basename(outFC), 
                                        "POLYLINE", template)
    # Access spatial reference of template to define spatial
    # reference of geometries
    spatial_reference = arcpy.Describe(template).spatialReference
    # Open an insert cursor for the new feature class
    cur = arcpy.da.InsertCursor(outFC, ["SHAPE@"])
    # Create an array object needed to create features
    array = arcpy.Array()
    # Initialize a variable for keeping track of a feature's ID.
    ID = -1
    for coords in coords_list: 
        if ID == -1:
            ID = coords[0]
        # Add the point to the feature's array of points
        #   If the ID has changed, create a new feature
        if ID != coords[0]:
            cur.insertRow([arcpy.Polyline(array)])
            array.removeAll()
        array.add(arcpy.Point(coords[1], coords[2], ID=coords[0]))
        ID = coords[0]
    # Add the last feature
    polyline = arcpy.Polyline(array, spatial_reference)
    cur.insertRow([polyline])
except Exception as e:
    print(e)
finally:
    # Cleanup the cursor if necessary
    if cur:
        del cur

Es posible crear polígonos multiparte, entidades de polilínea y poligonales con anillos interiores mediante un conjunto de conjuntos y pasar esto a clases de Polygon y Polyline.

Creación de geometría a partir de listas de coordenadas

También es posible crear una geometría a partir de una lista de coordenadas. Este planteamiento puede ofrecer mejoras de rendimiento, ya que evita la sobrecarga de creación de objetos de geometría. Sin embargo, está limitada solamente a entidades que sean de parte simple y, en el caso de los polígonos, sin anillos interiores. Todas las coordenadas deben estar en las unidades de la referencia espacial de la clase de entidad.

En el siguiente ejemplo, se crea una sola entidad poligonal a partir de una lista de pares x,y.

import arcpy
import os
coordinates = [(-117.2000424, 34.0555514), 
               (-117.2000788, 34.0592066), 
               (-117.1957315, 34.0592309), 
               (-117.1956951, 34.0556001)]
# Create a feature class with a spatial reference of GCS WGS 1984
result = arcpy.management.CreateFeatureclass(
    arcpy.env.scratchGDB, 
    "esri_square", "POLYGON", spatial_reference=4326)
feature_class = result[0]
# Write feature to new feature class
with arcpy.da.InsertCursor(feature_class, ['SHAPE@']) as cursor:
    cursor.insertRow([coordinates])