<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 13, 2014 at 4:46 PM, Tyler Smith <span dir="ltr"><<a href="mailto:tyler@plantarum.ca" target="_blank">tyler@plantarum.ca</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":1n0" class="a3s" style="overflow:hidden">I'm trying to generate solar irradiance data for a set of points (which<br>


requires r.horizon and r.sun output). My data covers most of eastern<br>
North America at 20m by 20m resolution. This is beyond the capacity of<br>
my computer to process. However, I don't need complete maps for this<br>
area. I only need the beam values for ca. 10,000 point locations. I<br>
expect it would be possible to do this on my laptop, but I can't figure<br>
out how to limit the analysis to these spots. Is it possible? Does<br>
anyone have any suggestions?<br>
</div></blockquote></div><br></div><div class="gmail_extra">I'm not sure if I see all the details but if you can script in Python or Bash you can do something like iteration over all points using PyGRASS or v.out.ascii and textual processing (the later applies also for external files with points) and then set computation region using g.region around each point to sufficient extent. Then you compute the r.horizon and r.sun just in the specified (computational) region and store the results as maps and as text (using e.g. r.univar or v.what.rast). Does this make sense?<br>

</div></div>