Archivo *.bat para determinar sucesivamente NDVI, SAVI, LAI y εNB en cónsola de GRASS-QGIS con Windows XP

En el artículo anterior propuse un script de bash para determinar sucesivamente NDVI, SAVI, LAI y εNB en cónsola de GRASS-QGIS usando Linux Debian. Debido a que ocasionalmente uso Windows (XP y Seven), tenía curiosidad en saber si el script anterior, modificado como archivo *.bat, podía ser ejecutado en este sistema operativo sin problemas. Como tengo montadas las particiones NTFS de XP y Seven en Debian con privilegios de lectura/escritura, lo que hice fue copiar el dataset empleado en Debian (calabozo2) en Windows XP (C:\Documents and Settings\home_usuario\grassdata\calabozo2) y posteriormente ejecutar QGIS en Windows XP; cargando el dataset calabozo2 en ambiente de GRASS. La ejecución de la cónsola de GRASS tiene la apariencia que se observa en la siguiente imagen:

Como tenía mucho tiempo sin crear y ejecutar archivos *.bat, después de una ligera revisión de literatura, me di cuenta que éste podía funcionar (sólo con un cambio en la primera línea):

cls
echo "Calculando NDVI ..."
r.mapcalc "NDVI=(banda4-banda3)/(banda4+banda3+0.)"
echo "Hecho!"
echo "Calculando SAVI con L=0.1 ..."
r.mapcalc "SAVI_0.1=((banda4-banda3)/(banda4+banda3+0.1))*(1+0.1)"
echo "Hecho!"
echo "Calculando LAI..."
r.mapcalc "LAI=if(SAVI_0.1<0.1,0.00001,(if(0.1<SAVI_0.1 && SAVI_0.1<0.687,-log((0.69-SAVI_0.1)/0.59)/0.91,if(SAVI_0.1>0.687,6,0))))"
echo "Hecho!"
echo "Calculando eNBf ..."
r.mapcalc "eNB=if(LAI<3 && NDVI>0.,0.97+0.0033*LAI,if(LAI>=3 && NDVI>0.,0.98))"
r.mapcalc "eNBf=if(eNB==0,0.99,eNB)"
echo "Hecho!

Lo anterior lo edité como texto plano en C:\Archivos de programa\Quantum GIS Wroclaw y lo grabé como script.bat. Para ejecutarlo sólo bastó script [Enter] en cónsola de GRASS. El resultado fue exitoso; aunque la salida no tan elegante como en Linux (tampoco tengo interés en mejorarla):

Los resultados fueron idénticos a los de Linux; tal como se observa en la muestra de la imagen siguiente:

Esta entrada fue publicada en GRASS, QGIS, SIG, Software Libre. Guarda el enlace permanente.

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s