[GRASS-user] Temporal framework: calculating annual 5-day extremes
soerengebbert at googlemail.com
Thu Apr 6 02:31:24 PDT 2017
2017-04-06 11:09 GMT+02:00 RichardCooper <richtcooper at hotmail.com>:
> 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 year into
> a single raster.
What you need is a temporal moving window with the size of 5 days to
compute for each day the 5 day aggregate of the future.
You can convert your time series into a voxel dataset (3d raster) and
use r3.mapcalc with the neighbor index operator map[x][y][z] (if i
agg_map3d = map3d + map3d + ... map3d
Or you use t.rast.algebra  and the temporal neighbor operator strds[t]:
agg_strds = prec_strds + prec_strds + ... prec_strds
> To do this in t.rast.accumulate, I can see how to set a temporal cycle of 1
> year (cycle= "12 months"), but not sure how to specify such a rolling sum
> calculation of 5 days as described above. The default method is the 'mean'
> as indicated in r.series.accumulation? I'm not too sure how the accumulation
> is applied in the module.
> Best regards,
> View this message in context: http://osgeo-org.1560.x6.nabble.com/Temporal-framework-calculating-annual-5-day-extremes-tp5316014p5316076.html
> Sent from the Grass - Users mailing list archive at Nabble.com.
> grass-user mailing list
> grass-user at lists.osgeo.org
More information about the grass-user