Cortar archivos Raster-TIF con un script de Python en QGIS

En esta entrada se muestra el proceso para realizar una extracción por mascara (“cortar”) de varios archivos raster-TIF con un archivo shape, ejecutando un script en la consola de Python de QGIS.

El script puede ejecutarse también con el CMD (ventana del Símbolo del sistema) de Windows, pero en este post únicamente se mostrará cómo realizarlo directamente el QGIS.

Los programas para el funcionamiento de este script son: Windows / QGIS.

En la imagen que se ve a continuación, en el recuadro rojo están las tres líneas donde el usuario debe modificar los directorios correspondientes a sus propios datos. Nótese que dicho directorio se puede copiar directamente desde el navegador de Windows, y se pega en el script, anteponiendo la letra r en las comillas.

Por ejemplo, el directorio D:_Franz\04_cortarTIF_QGIS\datosFranz\TIF_originales
Quedaría así:
CarpetaTIForiginales = r”D:_Franz\04_cortarTIF_QGIS\datosFranz\TIF_originales”

Las otras partes del script se deben dejar sin modificaciones, incluso las primeras 11 líneas se pueden copiar hasta a la consola de Python de QGIS, y no habrá problema pues son líneas comentadas ya que tienen el símbolo # de comentarios.

Una vez modificados los directorios del script, se puede copiar a la consola de Python de QGIS, la cual se abre desde el menú Complementos > Consola de Python, o mediante el botón que tiene un ícono similar al de Python:

Luego se copia todas las líneas de código del script hasta la ventana de entrada, que se señala en verde en la imagen a continuación:

Luego se introduce uno o dos ENTER para que se ejecute. Si todo está en orden, los colores de las líneas que quedan en la parte superior tendrán color negro, azul o verde, si resultan líneas en color rojo, indica que se presentó algún error, por lo que se debe revisar el procedimiento y las modificaciones realizadas al script original, o bien, copiar las líneas de código nuevamente.

Es necesario también revisar la carpeta de los archivos de salida (raster-TIF cortados) y eliminar los archivos que se hayan generado, antes de ejecutar nuevamente el script.

Una contribución de: Antonio Pantoja

franzpc

2 thoughts on “Cortar archivos Raster-TIF con un script de Python en QGIS”

  1. Antonio Pantoja dice:

    Interesante post pero no termino de ver la necesidad de meter las cosas en un script de Python cuando se puede correr el comando de GDAL directamente desde la consola (que es en el fondo a lo que llama Python) y te ahorras incluso abrir QGIS.
    😉

    1. Antonio Pantoja C dice:

      Efectivamente, el método mostrado en este post es tan solo un más de los múltiples posibles.
      Gracias por comentar
      Saludos
      Antonio Pantoja C

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 422 suscriptores