Thanks Markus
Would it be more productive for me to get and compile GRASS 7.0? Given that
most of the addons which I am interested in are built into this version. On
Win7.
ta
--
View this message in context:
Ciao grass users.
I have three maps A, B, C.
I would use r.mapcalc to combine them according to this logic:
Put the values of A and if A is null you put the values of B and if B is
null you put the values of C.
I have tried in various ways without success. .. now I'm trying with a
double IF...
I think you almost got there.
the IF condition would be this:
r.mapcalc output = if((isnull(a)),b,if(isnull(b),c),a)
so if A is null, use B, otherwise test if B is null. if it is, use C,
otherwise use A.
(note: this is from the top of my mind, I didn't tried it)
Hope this helps.
Carlos
*r.patch* - Creates a composite raster map layer by using known category
values from one (or more) map layer(s) to fill in areas of no data in
another map layer.
The GRASS program *r.patch* allows the user to build a new raster map the
size and resolution of the current region by assigning
Hi Carlos and thanks.
A and B do not overlap, and I should take the values of A and B also
(always). Then.in areas that are empty, I would then put the values of
C.
I tried withr.mapcalc output = if((isnull(a)),b,if(isnull(b),a),c)
but I lose the values of C
--
View
Hi Sylvain and thanks.
I had not thought of r.patch and I think that works (but I have to
carefully examine the result)
I did so
r.patch input = a, b, c = output output_patch
Now I want to do it with r.mapcalc:)
--
View this message in context:
W dniu 16.01.2011 18:29, Gabriele N. pisze:
A and B do not overlap, and I should take the values of A and B also
(always). Then.in areas that are empty, I would then put the values of
C.
I tried withr.mapcalc output = if((isnull(a)),b,if(isnull(b),a),c)
but I lose the values of C
Hello all.
I've been trying to convert some of my bash scripts to python, but
some issues appeared:
1 - it seems that http://grass.osgeo.org/programming6/pythonlib.html is offline
2 - dealing with maps names. if I use find_file, I end up with the
full name of the map (that is, map name plus
just to add to the last email:
grray() isnot working.
this is what I get, with any raster:
a = garray.array()
a
array([[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.],
[ 0., 0., 0., ..., 0., 0., 0.],
...,
[ 0., 0., 0., ..., 0., 0.,