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

2 comentarios en “Recortar archivos Raster-Grid con #Python”

  1. 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

    • 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)

Los comentarios están cerrados.