Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-07 Thread RichardCooper
With a manual work around I was able to use t.rast.algebra, but needed to both manually create a relative STRDS and then register the output rasters from t.rast.algebra into the latter. The error encountered in using t.rast.algebra (I'm not sure if me or a bug) appears to be that an absolute

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-07 Thread Sören Gebbert
Patching t.rast.to.rast3 to support file input requires the modification of r.to.rast3 to implement file support. This may not happen in the near future. I would suggest to use t.rast.algebra or to compute the 5day rainfall extrema in 5 year chunks using t.rast.to.rast3 + r3.mapcalc. Use

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-07 Thread Moritz Lennert
On 07/04/17 12:45, RichardCooper wrote: Many thanks for the suggestions. I'd really like to use the voxel approach but am getting the following error on running t.rast.to.rast3 to create a voxel. I've increased the hard and soft limits for open files to 4 on my system, but I still am

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-07 Thread RichardCooper
I should add that t.rast.to.rast3 works for a STRDS containing 1000 layers, but if more than 5000 layers I get the abovementioned error after increasing the open file limits (to 40) on my notebook. The number of open files on the system: lsof | wc -l 156036 -- View this message in

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-07 Thread RichardCooper
Many thanks for the suggestions. I'd really like to use the voxel approach but am getting the following error on running t.rast.to.rast3 to create a voxel. I've increased the hard and soft limits for open files to 4 on my system, but I still am unable to convert a space time raster dataset

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-06 Thread Sören Gebbert
Hi, 2017-04-06 11:09 GMT+02:00 RichardCooper : > I have a time series of rainfall data, and for each year I want to calculate > the five-day period with maximum rainfall. So I would need to calculate the > sum of day1 to day5, then day2 to day6, then day3 to day7 etc for

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-06 Thread Mira Kattwinkel
Hi Richard t.rast.aggregate might be the function you are looking for. It has the method 'max'. For the cycling you might need to somehow loop over the data with different starting days. Best, Mira On 06/04/17 11:09, RichardCooper wrote: I have a time series of rainfall data, and for

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-06 Thread RichardCooper
I have a time series of rainfall data, and for each year I want to calculate the five-day period with maximum rainfall. So I would need to calculate the sum of day1 to day5, then day2 to day6, then day3 to day7 etc for the whole year, and then output the maximum grid cell 5-day values for each

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-06 Thread Sören Gebbert
Hi, i am not sure if i understand the problem correctly. However, you can use t.rast.accumulate [1] to create the rolling sum for an arbitrary interval (5 days, one year?) and then use t.rast.aggregate to compute the yearly maximum value time series based on the time series output of

Re: [GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-06 Thread Sitansu Pattnaik
r.series for the computation and g.mlist for selecting the rasters is how I have used in the past. Sitansu On Thu, Apr 6, 2017 at 10:01 AM, RichardCooper wrote: > Hi, > > I have a data set containing multiple years of daily raster layers and > would > like to aggregate

[GRASS-user] Temporal framework: calculating annual 5-day extremes

2017-04-05 Thread RichardCooper
Hi, I have a data set containing multiple years of daily raster layers and would like to aggregate and output annual raster layers of 5-day extremes (maxima). Essentially, for each grid cell, I need to calculate rolling 5-day sums of each year and then find the annual max of the latter sums, and