Mapa Solar promedio mensual (R.Sun) en pasto

Question

Estoy tratando de crear mapas medios mensuales como se describe en R.Sun en la hierba. Tengo mi DSM. En la lectura de la página de OSGEO menciona que es posible sin mostrar cómo. ¿Ideas? Supongo que tendrá que escribir un guión.

Answer 1

Aquí un enlace a un simple shell script (en lugar de la solución de bash del comentario de José): http://www.grassbook.org/examples_menu2nd.php

Llegar: rsun_crop.sh

Answer 2

Supongo que necesitará en bucle durante días en un mes y medio entonces. Aquí es un script bash. Las ideas es calcular la radiación solar individual y luego almacenar los nombres de archivo. Puede procesarlos r.series para calcular el promedio mensual.

#!/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