[GRASS-user] Aggregating daily maps in relative strds per month

Nikos Alexandris nik at nikosalexandris.net
Mon May 18 10:51:58 PDT 2015


Nikos Alexandris:

> Having absolute timestamps would make life easier with the
> fantastic temporal framework. I'll re-do again my workflow and decide
> upon a year. After all, I did extract the following for a specific year
> (don't remember which one it was now :-?):

I think it was 2009. So, I am trying to go 'absolute':

# list maps
g.list raster pattern=global_rad_zero_[0-9][0-9][0-9]

global_rad_zero_001
..
global_rad_zero_365

# erase timestamps
for Map in `g.list raster pattern=global_rad_zero_[0-9][0-9][0-9]` ;do
r.timestamp map=${Map} date=none ;done

# create an absolute strds
t.create out=global_rad_zero title="Clearsky global irradiation"
description="Clearsky global irradiation on a horizontal surface raster
maps
[Wh.m-2.day-1]" --o

# register maps (which have no timestamp!)
t.register -i input=global_rad_zero maps=`g.list raster pattern=global_rad_zero_[0-9][0-9][0-9] sep=,` start="2009-01-01" increment='1 days' --o

0..ERROR: invalid literal for int() with base 10: '2009-01-01'

Why?  Nikos


More information about the grass-user mailing list