[GRASS-dev] Re: [GRASS GIS] #498: r.sun2 commissioning trials
GRASS GIS
trac at osgeo.org
Wed Jan 25 22:27:37 EST 2012
#498: r.sun2 commissioning trials
--------------------+-------------------------------------------------------
Reporter: hamish | Owner: hamish
Type: defect | Status: assigned
Priority: major | Milestone: 6.4.2
Component: Raster | Version: svn-develbranch6
Keywords: r.sun | Platform: All
Cpu: All |
--------------------+-------------------------------------------------------
Comment(by hamish):
Replying to [comment:19 dylan]:
> The output looks correct when slope and aspect are specified.
One thing I notice about sa_diff_rad.015.jpg is that the diffuse output
there seems to incorporate the effects of direct shading, while
diff_rad.015.jpg does not.
I would ~ guess that the amount of diffuse light reaching a grid cell
would be directly proportional to the visible sky. (stable throughout the
day; should be independent of aspect (sa_diff.015.jpg is correlated with
aspect); perhaps we should try a new r.surf.volcano surface and see how
much diffuse light gets into a deep caldera or well, with the only visible
sky a small amount directly above)
Ultimately I think a formal verification of the model using (multiple)
real light meters is going to be necessary, but to get the deployment
design right we're going to have to become experts in the algorithm. (e.g.
if a light meter is shaded by a reflective wall, its beam will be 0 but
will it's diffuse be increased vs a black wall? how many bounces does
reflected light follow? just one?)
...
here's another pseudo-data test-case, I'll post results once it finishes.
{{{
#spearfish (further north than NC so more defined shadows)
g.region -d
r.mapcalc "undulates = (2 + sin( row() * 2 ) + cos( col() * 2 )) * 500"
r.colors undulates color=bcyr
nviz undulates
DAY=355
time r.sun -s elevin=undulates day=$DAY step=0.05 \
beam_rad=rad_test.und.355.beam \
diff_rad=rad_test.und.355.diff \
refl_rad=rad_test.und.355.refl
}}}
Hamish
--
Ticket URL: <https://trac.osgeo.org/grass/ticket/498#comment:22>
GRASS GIS <http://grass.osgeo.org>
More information about the grass-dev
mailing list