No subject

Amy Keeley akeeley at ncdc.noaa.gov
Thu Jan 19 10:48:38 EST 1995


Can anyone give me an idea
as to why the following MAPCALC equation will
not work?
_________________________________________________________
Mapset <akeeley> in Location <uk> 
GRASS 4.1 > r.mapcalc

MY INPUT:
>>mapcalc> sum='if((min("u861002.dat","u861003.dat","u861004.dat",\
>>...> "u861005.dat","u861006.dat","u861007.dat")),(((9.0/24.0)*\
>>...> "@u861002.dat")+"@u861003.dat"+"@u861004.dat"+"@u861005.dat"\
>>...> +"@u861006.dat"+((15.0/24.0)*"@u861007.dat"+1.499)))

PROGRAM RESPONSE:
>>you have confused me
___________________________________________________________

Thanks,
Amy Keeley
akeeley at ncdc.noaa.gov

p.s. "u86100?.dat" are raster files.  What I'm
trying to do is:  [if(x,a)...a if x not zero, 0 otherwise]
|  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |  |
V  V  V  V  V  V  V  V  V  V  V  V  V  V  V  V  V 
So if the smallest value out of all of the raster files is
equal to a zero, the value in output map 'sum' is zero.  If the smallest
value out of all of the raster files if not a zero, the 
values are summed, with a simple multiplication done on the
first and last values, to produce the map 'sum'.



More information about the grass-user mailing list