Hi Hanlie,<div><br></div><div>I've this before with 11 gauge stations, also using TRMM data. You'll probably want to write a script and use r.mask to mask the interpolation map to your grid cell, then do r.sum on your masked area. After that just divide by the area to get the average rainfall. Given the small area and limited number of gauges this shouldn't take very long to compute for 2 years of data. You should be able to fairly easily output your desired time series to a text file.</div>
<div><br></div><div>Best,</div><div>Nick<br><br><div class="gmail_quote">On Mon, Aug 15, 2011 at 3:00 PM, Hanlie Pretorius <span dir="ltr"><<a href="mailto:hanlie.pretorius@gmail.com">hanlie.pretorius@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Hi,<br>
<br>
I have rainfall measurements from 6 gauges that I want to interpolate<br>
to an areal value (a 'surface'), so that I can compare the the<br>
interpolated gauge values to a satellite rainfall estimate that covers<br>
a grid cell of 28kmx28km. Two of the gauges are outside, but close to<br>
the border of the grid cell. Therefore, I also need to clip the<br>
interpolated surface to the grid cell and to get the average of the<br>
surface value in this clipped surface.<br>
<br>
A search on Google revealed that v.vol.rst is probably a good choice<br>
for interpolating the rainfall because it takes elevation into<br>
account.<br>
<br>
However, for each rain gauge I have 730 values representing a daily<br>
measurement over two years. As output from GRASS, I need a text file<br>
with the interpolated rainfall values for each day in my time series.<br>
So, I was wondering if there is an 'easy' way to get my output without<br>
creating 730 vectors and rasters?<br>
<br>
Regards<br>
Hanlie<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>
</blockquote></div><br></div>