dir_x and dir_y are maps, and as far as I know they can not use in map[x,y]
since it expects relative coordinates for the moving window...even though
your raster values are in [0,1,-1]
Christian.
--------------------------------------------------
From: "Achim Kisseler" <[email protected]>
Sent: Thursday, January 28, 2010 2:13 PM
To: <[email protected]>
Subject: [GRASS-user] r.mapcalc question
Hello grass-users,
I have a problem running a r.mapcalc procedure:
given two maps
dir_x,
dir_y
with 0,1,-1 entries
I want to run:
r.mapcalc "the_map=if(isnull(the_map),the_map[dir_x,dir_y],null())"
in order to make a r.water.outlet on many many areas at the same time.
The error I get is:
syntax error, unexpected NAME, expecting INTEGER or '-'
why does dir_x/dir_y does not return an integer here? Using int(dir_x)
also does not work.
Until now I use a r.mapcalc command with many "if", checking the direction
of every pixel every time. But this is quite slow.
Achim
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
