#!/bin/sh
# Author: Marcello Benigno, May/2013; copyright: GPL >= 2
# Purpose: Create a reclass slope map
# Usage: ./declividade.sh
#-----------------------------------------------------------
#verifica se o GRASS está em execução:
if  [ -z "$GISBASE" ] ; then
    echo "ERRO: O GRASS não está em execução!"
 	exit 1
fi
echo -n "Digite o nome do layer raster (MDT): "
read _MTD
#verifica se o arquivo existe
g.findfile element=cell file="$_MTD" > /dev/null
if [ $? -eq 1 ] ; then
	g.message -e "O layer $_MTD não foi encontrado neste mapset!"
exit 1 
else
#definição da região da análise (dem SRTM):
g.region rast=$_MTD
#criação do mapa de declividade em porcentagem
r.slope.aspect --o elevation=$_MTD slope=slope format=percent
#reclassificação em função das classes de declividade
r.mapcalc <<EOF
class_1 = if(slope <= 2.0, 1, 0)
class_2 = if(slope > 2.0 && slope <= 5.0, 2, 0)
class_3 = if(slope > 5.0 && slope <= 10.0, 3, 0)
class_4 = if(slope > 10.0 && slope <= 15.0, 4, 0)
class_6 = if(slope > 15.0 && slope <= 45.0, 6, 0)
class_7 = if(slope > 45.0 && slope <= 70.0, 7, 0)
class_8 = if(slope > 70.0, 8,0)
slope.reclass = class_1 + class_2 + class_3 + class_4 + class_6 + class_7 + class_8
EOF
#eliminação das classes
g.remove --q rast=class_1,class_2,class_3,class_4,class_6,class_7,class_8
echo "------------------------------------------------" 
echo "O Layer 'slope.reclass' foi criado com sucesso."
echo "------------------------------------------------"
fi
exit 0