Le portail francophone de la géomatique

Bonjour,

J'essaie de réaliser un calcul avec l'outil r.mapcalc pour obtenir un raster de points de rupture de pente. Pour cela, je veux calculer la différence (en valeur absolue) de pente entre un point et les points de son voisinage (voisinage 3x3). Si cette différence est supérieure à une valeur seuil, le point est considéré comme point de rupture de pente.

Comme décrit dans le manuel de l'outil, je crée la fonction me permettant d'obtenir ces points, et je rentre donc la ligne de commande suivante dans la console SHELL :

r.mapcalc rs = if(abs(s@t-s@t[-1,-1]>10)||abs(s@t-s@t[-1,0]>10)||abs(s@t-s@t[-1,1]>10)||abs(s@t-s@t[0,1]>10)||abs(s@t-s@t[1,1]>10)||abs(s@t-s@t[1,0]>10)||abs(s@t-s@t[1,-1]>10)||abs(s@t-s@t[0,-1]>10)

rs: raster en sortie
s : raster de pente calculé via param.scale ou slope.aspect

Le message d'erreur suivant apparaît : "n'est pas reconnu en tant que commande interne ou externe, un programme exécutable ou un fichier de commandes". J'imagine que la syntaxe n'est pas bonne, mais je n'arrive pas à trouver mon erreur.

Quelqu'un peut-il m'aider ?

Merci d'avance

Arnaud