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

Compartir:

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 *

Este sitio usa Akismet para reducir el spam. Aprende cómo se procesan los datos de tus comentarios.

Carrito

Últimas publicaciones

  • Simulación 3D de inundaciones en ArcGIS
  • Descargar guía Pix4Dmapper en PDF
  • ¿Cómo seleccionar el mejor método de interpolación?
  • ¿Cuál es la Primera ley de la Geografía de Tobler?
  • Visualización avanzada en 3D usando QGIS
  • Evaluando parámetros de suelos interpolados en QGIS 3
  • Cuadriculas de latitud, longitud y sistema de coordenadas
  • Interpolación espacial en QGIS 3
  • Meridiano de Greenwich (meridiano principal)
  • ¿Qué son las proyecciones de mapas, y por qué engañan al ojo humano?
  • Ventajas y desventajas de los mapas de papel
  • Norte Magnético vs Norte Geográfico (verdadero)

Suscríbete al blog

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

Únete a otros 623 suscriptores