<p dir="ltr">Hi Nikos!</p>
<p dir="ltr">If the data in fact represent daily data, why don't you register it as absolute time?? In that way, you would just loop over strds...</p>
<p dir="ltr">for STRDS in global_rad_zero beam_rad_zero diff_rad_zero insol_time<br>
do     <br>
t.rast.aggregate input=${STRDS} output=${STRDS}_monthly_average basename=${STRDS}_average granularity="1 month" method=average<br>
done</p>
<p dir="ltr">Cheers,<br>
Vero</p>
<div class="gmail_quot<blockquote class=" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Is there a better way to aggregate a *relative* spatio-temporal raster<br>
data set composed by daily maps (output from r.sun.daily) in monthly<br>
average strds-es?<br>
<br>
At the moment I use the following csv and script:<br>
<br>
--%<---<br>
# month_start_end_day.csv<br>
January|1|31<br>
February|32|59<br>
March|60|90<br>
April|91|120<br>
May|121|151<br>
June|152|181<br>
July|182|212<br>
August|213|243<br>
September|244|273<br>
October|274|304<br>
November|305|334<br>
December|335|365<br>
--->%--<br>
<br>
and<br>
<br>
--%<---<br>
# loop over spatio-temporal raster data sets of interest<br>
for STRDS in global_rad_zero beam_rad_zero diff_rad_zero insol_time<br>
do<br>
    # loop over months (start day, end day, month name) in csv<br>
    for LINE in `cat month_start_end_day.csv`<br>
<br>
    do<br>
        # set positional parameters<br>
        set -- $(echo ${LINE} |tr "|" " ")<br>
<br>
        # aggregate<br>
        t.rast.aggregate input=${STRDS} output=${STRDS}_$1_average basename=${STRDS}_$1_average where="start_time >= $2 AND start_time <= $3" granularity=1 --o<br>
    done<br>
<br>
done<br>
--->%--<br>
<br>
Any tip highly appreciated, Nikos<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</div>