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 <[email protected]>: > 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 _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
