[GRASS-user] problem with t.rast.mapcalc (grass71 daily ubuntu 14.04)

Wed Oct 7 06:54:18 PDT 2015

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 <leohardtke at gmail.com>:

> Nice!
> El 06/10/2015 16:17, "Veronica Andreo" <veroandreo at gmail.com> 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" <leohardtke at gmail.com>
>> 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 <veroandreo at gmail.com>:
>>>
>>>> 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" <leohardtke at gmail.com>
>>>> 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
>>>>> grass-user at lists.osgeo.org
>>>>> 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
>>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20151007/dfeae817/attachment.html>