Monthly Average Solar Maps (R.Sun) in Grass

Question

I'm trying to create monthly average maps as described in R.Sun in Grass. I have my DSM. In reading the page on OSGEO it mentions that it is possible without showing how. Ideas? I assume I will need to write a script.

We need to update the answers. Things are easier now in GRASS7 within it's temporal framework. — Nikos Alexandris, Mar 29 '15 at 13:15

score 4 · Answer 1 · answered Jun 13 '11 at 18:38

I guess that you need to loop over days in a month and then average. Here's attempt at a bash script. The ideas is to calculate the individual solar radiation, and then store the filenames. You can process those with r.series to calculate the monthly average.

#!/bin/bash
mdoys=( 31 59 90 120 151 181 212 243 273 304 334 365 ) 
months=( 'January' 'February' 'March' 'April' 'May' 'June' 'July' 'August' 'September' 'October' 'November' 'December' )
this_month=0
files=""
for ((doy=1;doy<365;doy++)) do

  r.sun -s elev=elevation.dem slop=slope asp=aspect \
        beam_read=beam_${doy} diff_rad=diff_${doy} day=${doy}

  if [[ ${doy} -lt ${mdoys[${this_month}]} ]] ;  then
   files=${files},beam_${doy}
  fi
  if [[ ${doy} -eq ${mdoys[${this_month}]} ]] ; then
   files=${files},beam_${doy}
   echo ${files} > rasters_${months[$this_month]}
   this_month=$((this_month+1))
   # Use r.series to process the files in rasters_January, for example
  fi
done

BTW, I haven't tested it on GRASS, proceed with caution! :) The logic seems to hold, though! — Jose, Jun 13 '11 at 23:20
Just saw this. +1 for indexing in bash. Looks like Python, a bit :-). — Nikos Alexandris, Mar 29 '15 at 13:13

score 4 · Answer 2 · answered Jun 27 '11 at 19:09

4

Here a link to a plain shell script (rather than the bash solution of Jose's comment): http://www.grassbook.org/examples_menu2nd.php

Get there: rsun_crop.sh

answered Jun 27 '11 at 19:09

markusN

12,976
31
48

Monthly Average Solar Maps (R.Sun) in Grass

2 Answers2

Linked