¿Calcular coordenadas en los vértices de un SHP de polígonos? - Parte 1 🤩🤯🤓 #python

Únete a este canal para acceder a sus beneficios:

Facebook "El Geógrafo"

Blog "El Rincón Geográfico"
________________________________________________________________
Código de Python

import arcpy

# Cambia la ruta según la ubicación de tu shapefile

# Agrega las columnas necesarias para almacenar las coordenadas de los vértices
for i in range(1, 5):
arcpy.AddField_management(shapefile_path, f"Vertice{i}_X", "DOUBLE")
arcpy.AddField_management(shapefile_path, f"Vertice{i}_Y", "DOUBLE")

# Función para extraer las coordenadas X e Y de los vértices de un polígono
def extract_vertices(polygon, num_vertices):
vertices = []
for part in polygon:
for point in part:
if len(vertices) == num_vertices:
break
if len(vertices) == num_vertices:
break
return vertices

# Procesa cada fila del shapefile y guarda las coordenadas en las nuevas columnas
for row in cursor:
polygon = row[0]
vertices = extract_vertices(polygon, 4)

# Guarda las coordenadas en las columnas correspondientes
for i, (x, y) in enumerate(vertices, 1):
row[i] = x
row[i + 4] = y


print("Columnas de coordenadas de vértices agregadas exitosamente.")
# Columnas de coordenadas de vértices agregadas exitosamente.