[GRASS-user] t.rast.mapcalc error

Sören Gebbert soerengebbert at googlemail.com
Sun Jul 27 15:13:51 PDT 2014


Hi Veronica,
i have tested the strds/map comparison and it works for me with the
latest grass7.1 build. Can you please check if the max_my_strds map
exists and has valid values and that all maps in the strds exists and
have valid values?

Best regards
Soeren

2014-07-27 23:24 GMT+02:00 Veronica Andreo <veroandreo at gmail.com>:
> Hi Soeren, Markus and list,
>
> I'm having trouble with t.rast.mapcalc... these are the commands I'm using:
>
> # get the max value for each pixel in the strds
> t.rast.series input=my_strds method=maximum output=max_my_strds
>
> # then, i need to compare all values in the strds with the map of maximums
> and whenever it's the same value, save the start_doy, otherwise null
> t.rast.mapcalc -n inputs=my_strds output=date_max_my_strds \
> expression="if(my_strds == max_my_strds,start_doy(),null())"
> basename=date_max_my_strds
>
> # to get a map with dates of maximum values afterwards
> t.rast.series input=date_max_my_strds method=maximum
> output=series_date_max_my_strds
>
> This is the error I'm getting in newly compiled GRASS 7
>
> Start temporal sampling
> Start mapcalc computation
> Invalid map <max_cla_hants25_corrected_1>
> Parse error
> ERROR: parse error
> ERROR: Error while mapcalc computation
>
> However, when i follow the exact same steps, but using an older map of
> maximums (created a month ago or so), it works just fine... that's the only
> difference i've found... this used to work a couple of weeks ago...
>
> Am I doing something wrong? Or is it no longer possible to compare all maps
> in strds with a single map??? Please, help!
>
> Thanks much in advance!
>
> Best,
> Vero


More information about the grass-user mailing list