<div dir="ltr"><br><br>On Thu, Oct 11, 2018 at 4:45 PM Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br>><br>> 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><div>> having difficulty understanding the output.</div><div><br></div><div>With only 58 points, you will not get a lot of spatial detail, no matter what interpolation method you use. For example, using 58 points to create a surface with e.g. 1 million points is a waste of time: there is simply not enough spatial detail in the input. You can check, as you have done, by changing the resolution of the current region for the output and comparing the results. The only solution for more spatial detail is to get more input points.</div><div><br></div><div>Markus M<br></div><div><br></div><div>> I suspect that because my data</div>> 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">grass-user@lists.osgeo.org</a><br>> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a></div>