<div dir="ltr"><div>Hello Rich,</div><div><br></div><div>Few (if any) of the spatial interpolation algorithms require "evenly spaced" data. Ultimately, you want spatial estimates that produce the least bias and where the errors are gaussian distributed, which can only be determined using some process of bootstrapping where stations are removed from the spatial interpolation process and are then used to evaluate error at those locations; and then, repeated, with a new set of stations. This should be done with several spatial interpolation algorithms that are generally known to produce the best estimates.</div><div><br></div><div>If you are trying to produce estimates over an entire county, you really need stations well outside the county boundary to avoid having problems near the county periphry (interior of the county boundary). GRASS and R play really well together through spgrass7, so that open the use of many R based geostatistical packages.</div><div><br></div><div>Tom<br></div><div><br></div><br><div class="gmail_quote"><div dir="ltr">On Thu, Oct 11, 2018 at 10:45 AM Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On Mon, 24 Sep 2018, Rich Shepard wrote:<br>
<br>
> Please provide some thoughts on the work flow to do this within GRASS. I'm<br>
> digging into the gstat docs to get a sense of how to do this within R (and<br>
> I have the rgrass7 package working well; it imported the GRASS county<br>
> boundary map which I converted to a SpatialPolygonDataFrame.)<br>
<br>
   I have run v.surf.tps on the point map of 58 reporting stations and am<br>
having difficulty understanding the output. I suspect that because my data<br>
does not consist of evenly spaced sparsed points, I've accepted the default<br>
thinning of 1.5, and I don't (yet) have a DEM covering the entire county the<br>
output is not as usable as it could be. (I've changed the region's<br>
resolution to 100m.)<br>
<br>
   Regardless, I'm asking for input from those who have experience with the<br>
geostatistical modules on how to determine the most appropriate one for a<br>
given data set. I see that v.krige[1], v.kriging, and v.surf.tps are available<br>
and my lack of experiences with all of these hinders my making an informed<br>
decision.<br>
<br>
   Pointers to relevant documentation, as well as other advice, will be much<br>
appreciated.<br>
<br>
TIA,<br>
<br>
Rich<br>
<br>
[1] The rpy2 dependency for v.krige is not yet available for R-3.5.1.<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a></blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><br><div><br></div></div></div></div></div></div></div>