[GRASS-user] Interpreting r.watershed accumulation map

Laurent C. lrntct at gmail.com
Mon Oct 17 10:31:04 PDT 2016


Hello Sören,

The problem occurs when Itzï is calling:
"tgis.SpaceTimeRasterDataset(name).is_in_db()"

tgis.init() is called just before. "name" is a a formatted map or
strds ID (i.e including mapset)

It is difficult to understand what is happening because I can't
reproduce the problem.

Regards,
Laurent


2016-10-16 16:01 GMT-05:00 Sören Gebbert <soerengebbert at googlemail.com>:
> Hi,
> i don't know why this happens. It may be a bug in the temporal framework or
> in itzi. Testing the code in the temporal framework seems to work, but more
> testing is needed.
> However, the traceback is fine, but there must be an error message that will
> tell you what mapset is missing and what mapsets are in the temporal
> framework search path.
>
> Something like:
> """
> You have no permission to "
>        "access mapset <%(mapset)s>, or "
>        "mapset <%(mapset)s> has no temporal database. "
>        "Accessable mapsets are: <%(mapsets)s>
>
> """
>
> Best regards
> Soeren
>
>
>
> 2016-10-16 22:03 GMT+02:00 Rich Shepard <rshepard at appl-ecosys.com>:
>>
>> On Sat, 15 Oct 2016, Sören Gebbert wrote:
>>
>>> Please try to put all mapsets that you use in your analysis into the
>>> mapset search path. The temporal framework will only access strds in
>>> mapsets that are located in the current search path.
>>
>>
>> Soeren,
>>
>>   Current mapset is 'analyses':
>>
>>> g.mapset -p
>>
>> analyses
>>
>>   Some maps might be in other mapsets so ...
>>
>>> g.mapsets -p
>>
>> Accessible mapsets:
>> analyses PERMANENT features soils topography precipitation
>>
>> which is all the mapsets in the location.
>>
>>   The file, itzi-param.txt, contains:
>>
>> [time]
>> duration = 2:00:00
>> record_step = 00:05:00
>>
>> [input]
>> dem = be_cropped at topography
>> friction = n
>> rain = rain
>>
>> [output]
>> prefix = pop
>> values = h, wse, v, vdir, boundaries
>>
>> [statistics]
>> stats_file = pop.csv
>>
>>   And this is the traceback:
>>
>>
>>> itzi run itzi-param.txt
>>
>> Starting simulation of itzi-param.txt... WARNING: Error during execution:
>> Traceback (most recent call last):
>>   File "/usr/lib/python2.7/site-packages/itzi/itzi.py", line 183, in
>> sim_runner_worker
>>     sim_runner.run()
>>   File "/usr/lib/python2.7/site-packages/itzi/itzi.py", line 90, in run
>>     sim_param=self.conf.sim_param)
>>   File "/usr/lib/python2.7/site-packages/itzi/simulation.py", line 77, in
>> __init__
>>     self.gis.read(self.in_map_names)
>>   File "/usr/lib/python2.7/site-packages/itzi/gis.py", line 187, in read
>>     elif self.name_is_stds(self.format_id(map_name)):
>>   File "/usr/lib/python2.7/site-packages/itzi/gis.py", line 138, in
>> name_is_stds
>>     if tgis.SpaceTimeRasterDataset(name).is_in_db():
>>   File
>> "/usr/local/grass-7.3.svn/etc/python/grass/temporal/abstract_dataset.py",
>> line 368, in is_in_db
>>     return self.base.is_in_db(dbif)
>>   File "/usr/local/grass-7.3.svn/etc/python/grass/temporal/base.py", line
>> 315, in is_in_db
>>     dbif.execute(sql, mapset=self.mapset)
>>   File "/usr/local/grass-7.3.svn/etc/python/grass/temporal/core.py", line
>> 981, in execute
>>     self._create_mapset_error_message(mapset)))
>>   File
>> "/usr/local/grass-7.3.svn/etc/python/grass/pygrass/messages/__init__.py",
>> line 271, in fatal
>>     sys.exit(1)
>> SystemExit: 1
>>
>> Simulations complete. Elapsed times: itzi-param.txt: 0:00:01 Total:
>> 0:00:01 Average: 0:00:01
>>
>>   Please suggest what I should try modifying.
>>
>> Thanks,
>>
>>
>> Rich
>> _______________________________________________
>> grass-user mailing list
>> grass-user at lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/grass-user
>
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user


More information about the grass-user mailing list