Ops! My bad! Sorry... Thanks Sören for clarifying this issue...
--- #2735: t.rast.mapcalc input problem --------------------------+---------------------------- Reporter: leohardtke | Owner: grass-dev@… Type: defect | Status: reopened Priority: normal | Milestone: 7.0.1 Component: Temporal | Version: svn-trunk Resolution: | Keywords: t.rast.mapcalc CPU: Unspecified | Platform: Linux --------------------------+---------------------------- Changes (by huhabla): * status: closed => reopened * resolution: fixed => Comment: This bug was not fixed in t.rast.mapcalc. A similar bug was fixed in t.rast.extract. However, fixing the bug in t.rast.mapcalc is more difficult since it uses a simple search and replace approach to substitute the input STRDS with map names and is not aware of what a space time dataset is. If the names of input and output space time datasets are include each other, then the module does not work properly. Choosing a specific order of the input STRDS in the input parameter (i.e: input=ndvi_smooth,QA_mask,ndvi) may reduce the risk of wrong substitution. However, please use t.rast.algebra instead. This module has other issues but is able to correctly identify space time datasets. Ticket URL: <https://trac.osgeo.org/grass/ticket/2735#comment:3> --- Vero 2015-10-06 18:16 GMT-03:00 Leonardo Hardtke <[email protected]>: > Nice! > El 06/10/2015 16:17, "Veronica Andreo" <[email protected]> escribió: > >> Hola Leo, >> >> The problem with basenames was fixed in r66308 ( >> https://trac.osgeo.org/grass/changeset/66308). I've been using basename >> equal to strds name and works as expected. >> >> Enjoy! >> Vero >> >> >> El sep 9, 2015 3:38 PM, "Leonardo Hardtke" <[email protected]> >> escribió: >> >>> Thaks >>> If it is a bug, I should fill a bug report so that some developer can >>> fix it! >>> >>> Cheers >>> >>> 2015-09-09 15:17 GMT-03:00 Veronica Andreo <[email protected]>: >>> >>>> Hi Leo, >>>> >>>> Yes, I had the same problem some time ago and the solution was to >>>> change the name, as you mentioned... >>>> >>>> Maybe a bug? Dunno :) >>>> >>>> Cheers, >>>> Vero >>>> >>>> El sep 9, 2015 3:05 PM, "Leonardo Hardtke" <[email protected]> >>>> escribió: >>>> >>>>> Dear GRASS users (and developers?!), >>>>> >>>>> I'm facing a problem with an expression in the t.rast.mapcalc >>>>> module... >>>>> >>>>> t.rast.mapcalc inputs=QA_mask,ndvi,ndvi_smooth_spacetime >>>>> expression="if(QA_mask==0,ndvi,ndvi_smooth_spacetime)" output=ndvi_fixed >>>>> basename=ndvi_fixed >>>>> It throws the following error: >>>>> Starting temporal sampling... >>>>> Starting mapcalc computation... >>>>> Invalid map >>>>> <MOD13Q1.A2009017.h12v12.005.2009035110230_01_4326_smooth_spacetime_0> >>>>> Parse error >>>>> ERROR: parse error >>>>> ERROR: Error while mapcalc computation >>>>> >>>>> The QA_mask maps are: >>>>> QA_mask_01|PERMANENT|2009-01-01 00:00:00|2009-01-17 00:00:00 >>>>> QA_mask_02|PERMANENT|2009-01-17 00:00:00|2009-02-02 00:00:00 >>>>> ...... >>>>> >>>>> The ndvi_smooth_spacetime maps are: >>>>> ndvi_smooth_spacetime_0|PERMANENT|2009-01-17 00:00:00|2009-02-02 >>>>> 00:00:00 >>>>> ndvi_smooth_spacetime_1|PERMANENT|2009-02-02 00:00:00|2009-02-18 >>>>> 00:00:00 >>>>> ........... >>>>> >>>>> The ndvi maps are >>>>> MOD13Q1.A2009001.h12v12.005.2009019145228_01_4326|PERMANENT|2009-01-01 >>>>> 00:00:00|2009-01-17 00:00:00 >>>>> MOD13Q1.A2009017.h12v12.005.2009035110230_01_4326|PERMANENT|2009-01-17 >>>>> 00:00:00|2009-02-02 00:00:00 >>>>> M >>>>> >>>>> As you see, there is no such a map like.. >>>>> MOD13Q1.A2009017.h12v12.005.2009035110230_01_4326_smooth_spacetime_0 >>>>> >>>>> I guess the module is mixing ndvi (a temporal dataset) with >>>>> _smooth_spacetime... the second part of the name.... >>>>> >>>>> I know it work ok by renaming the ndvi_smooth_spacetime by something >>>>> else... but... is this a bug or shoul I use more '`" to protect names !? >>>>> >>>>> Thanks! >>>>> -- >>>>> Dr. Leonardo A. Hardtke >>>>> Laboratorio de Teledetección y S.I.G. >>>>> Centro Nacional Patagónico (CONICET) >>>>> Bvd. Brown 2825, 9120 >>>>> Puerto Madryn, Chubut, Argentina >>>>> >>>>> _______________________________________________ >>>>> grass-user mailing list >>>>> [email protected] >>>>> http://lists.osgeo.org/mailman/listinfo/grass-user >>>>> >>>> >>> >>> >>> -- >>> Dr. Leonardo A. Hardtke >>> Laboratorio de Teledetección y S.I.G. >>> Centro Nacional Patagónico (CONICET) >>> Bvd. Brown 2825, 9120 >>> Puerto Madryn, Chubut, Argentina >>> >>
_______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
