Cortar shapefiles con Python y QGIS

Este procedimiento permite cortar múltiples shapefiles usando otro de menor extensión espacial (máscara), mediante la ventana de procesamiento de Python en la aplicación QGIS (Consola de Python).

clip phyton qgis
Para desplegar la Consola de Python en QGIS, hay dos métodos: desde el menú Complementos > Consola de Python, o mediante el botón que tiene un ícono similar al de Python:
consola de python en qgis
Cuando aparece la consola de Python, el área de trabajo se divide en tres. En la parte central está el cursor esperando la entrada de código (>>>), bien sea escribiendo directamente o copiando el script desde algún gestor de texto. Y en la parte superior se ve el desarrollo de las órdenes introducidas (descargar los archivos de práctica).

En este caso, por ejemplo, se pegó todo el script que hace posible el procesamiento, se pulsó ENTER, y la consola está esperando ( … ) la introducción de más código o la introducción de otro ENTER para ejecutar.
Insertar comandos Python en QGIS

Se introduce entonces otro ENTER, y el programa empieza el procesamiento. Si todo es correcto aparecen mensajes OUTPUP indicando la dirección de los archivos cortados mediante la máscara. Si hay algún error aparecen ventanas de advertencia o mensajes en color rojo. En este caso, el proceso fue exitoso:
Importar comandos Python
En el espacio de respuesta de los procesos (parte superior de la Consola de Python), se ve tres líneas de color azul, que son las que cada usuario deberá cambiar y adaptar según los directorios donde se encuentren sus datos. Esto se debe hacer previamente a la introducción en la consola, puede ser en un gestor de texto como lo es Notepadd ++. Es preferible que los nombres de las carpetas no contengan espacios. Notar también que esos directorios debes estar entre comillas, tal como se ve en la imagen:
Resultados clip python qgis
Una contribución de: Antonio Pantoja

franzpc

2 thoughts on “Cortar shapefiles con Python y QGIS”

  1. Antonio Pantoja dice:

    Un problema que se puede presentar, es que los directorios no sean reconocidos. En cuyo caso se debe anteponer la letra r, de tal manera que queden así:
    DirectorioShapes = r”D:\_Franz\03_cortar_shapes\datosFranz\ArchivosShape”
    ShapeMascara = r”D:\_Franz\03_cortar_shapes\datosFranz\ArchivosShape\Santiago_de_cali.shp”
    DirectorioSalida = r”D:\_Franz\03_cortar_shapes\datosFranz\ArchivosShapeCortados”

  2. Edgar dice:

    Hola Antonio. Muchas gracias por el post.

    El problema que me surge es que me dice que ‘global features is not defined’

    ¿Qué puede ser?

    Saludos

Deja un comentario

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *

Carrito

Últimas publicaciones

  • Crear automáticamente un mapa de uso del suelo en QGIS
  • Dibujando en QGIS con mayor precisión que AutoCAD
  • Identificando y cuantificando cambios del uso de la tierra en línea
  • Monitoreo ambiental online utilizando un Geoportal
  • Manual PDF Prácticas de QGIS
  • Descargar ArcGIS Pro
  • Libro PDF: Teledetección Espacial
  • Diferencias entre un DSM, DEM & DTM
  • ¿Qué son las distancias geodésicas?
  • El GeoAmor mueve millones
  • Limitaciones o desafíos de los SIG
  • ¿Cómo calcular 18 índices radiométricos en un solo paso?

Suscríbete al blog

Introduce tu correo electrónico para recibir las últimas publicaciones.

Únete a otros 418 suscriptores