Recortar archivos Raster-Grid con #Python

Frecuentemente se suele extraer áreas de un archivo Raster-Grid, a través de una máscara que determina la superficie de recorte, desechando así zonas que no son de interés para un estudio en particular.

En ArcGIS se dispone de la herramienta Extract by Mask ubicada en: (ArcToolbox > Spatial Analyst Tools > Extraction > Extract by Mask). Sin embargo, usando el lenguaje de programación Python y las funcionalidades de ArcGIS es posible automatizar este proceso cuando se requiere automatizar o cortar varios archivos Raster-Grid.

Para realizar la presente práctica descargar el Script extract by mask.

Script python ArcGIS

Para el funcionamiento de este script se requiere Windows / Python / ArcGIS 10. 0 superior.

Los directorios donde se alojan los archivos no deben ser muy largos, y es imprescindible que no haya espacios en los nombres de las carpetas (se recomienda usar un guion bajo ( _ ) para separar una palabra de otra al momento de nombrar las carpetas).

Para ejecutar el script de manera rápida, se debe abrir la consola de comandos pulsando las teclas Windows + R, y escribiendo CMD.

En la ventana que se abre MS-DOS, arrastrar el script ejecutable de Python, dar un espacio y pulsar Enter.

Script python en MS-DOS

En el momento de la ejecución del script el sistema solicita cada parámetro de entrada:

  • Una carpeta con al menos dos archivos Raster-Grid.
  • Un archivo en formato Raster-Grid o Shape, que será la máscara o límite de la zona a extraer, y
  • Una carpeta donde se almacenarán para los nuevos archivos Raster-Grid “recortados” con el script.

Ingresar párametros python

Una vez introducidos los tres parámetros, el sistema procesa los archivos de la carpeta de entrada y en la ventana se arroja el mensaje “Proceso terminado”.

Resultado python

Finalmente revisar el directorio de salida para ver los archivos procesados. Para mayor claridad se recomienda visualizar el siguiente vídeo.

Una contribución de: Antonio Pantoja

franzpc

2 thoughts on “Recortar archivos Raster-Grid con #Python”

  1. jose lopez dice:

    buenas noches amigo, interesante tu blog, hace dias que he tenia un problema con el arcgis y resulta que al momento de terminar el mapa y realizar la grid, esta no me produce una grid me da error y me saca o simplemente no me sale la grid, eso me preocupa ya he realizado 2 mapa y por mala suerte en ninguno de los 2 e podido realizar la grid

    1. franzpc dice:

      Revisa los valores del Grid, puede que simplemente sea un problema de visualización que se puede chequear en las propiedades, pestaña Symbology del grid (percent clip, minimum maximum)

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