<div dir="ltr"><br><br>On Tue, Sep 25, 2018 at 2:06 PM Thomas Adams <<a href="mailto:tea3rd@gmail.com">tea3rd@gmail.com</a>> wrote:<br>><br>> Hello all,<br>><br><div>> I think there is sufficient evidence to recommend using elevation as a secondary variable for spatial interpolation of precipitation <br></div><div><br></div><div>as an alternative to kriging, TPS interpolation also allows to use covariables, e.g. with v.surf.tps or r.resamp.tps (both are addons)</div><div><br></div><div>Markus M<br></div><div><br></div><div>> using, for instance, Universal Kriging or spatial drift, or a method such as PRISM in the USA (<a href="http://www.prism.oregonstate.edu/">http://www.prism.oregonstate.edu/</a>) -- PRISM estimates are used widely in the US by the National Weather Service and others. See:</div>><br>> <a href="ftp://140.115.123.47/download/%A6h%C5%DC%B6q%A6a%BD%E8%B2%CE%ADp%BE%C7/01_Geostatistical%20approaches%20for%20incorporating%20elevation%20into%20the%20spatial%20interpolation%20of%20rainfall.pdf">ftp://140.115.123.47/download/%A6h%C5%DC%B6q%A6a%BD%E8%B2%CE%ADp%BE%C7/01_Geostatistical%20approaches%20for%20incorporating%20elevation%20into%20the%20spatial%20interpolation%20of%20rainfall.pdf</a><br>><br>> <a href="https://www.sciencedirect.com/science/article/pii/S0022169418303548">https://www.sciencedirect.com/science/article/pii/S0022169418303548</a><br>><br>> <a href="https://www.hydrol-earth-syst-sci.net/10/197/2006/hess-10-197-2006.pdf">https://www.hydrol-earth-syst-sci.net/10/197/2006/hess-10-197-2006.pdf</a><br>><br>> But, there is this, e.g., which supports Micha's statement:<br>><br>> <a href="ftp://140.115.123.47/download/%A6h%C5%DC%B6q%A6a%BD%E8%B2%CE%ADp%BE%C7/01_Geostatistical%20approaches%20for%20incorporating%20elevation%20into%20the%20spatial%20interpolation%20of%20rainfall.pdf">ftp://140.115.123.47/download/%A6h%C5%DC%B6q%A6a%BD%E8%B2%CE%ADp%BE%C7/01_Geostatistical%20approaches%20for%20incorporating%20elevation%20into%20the%20spatial%20interpolation%20of%20rainfall.pdf</a><br>><br>> Best,<br>> Tom<br>><br>> On Mon, Sep 24, 2018 at 2:19 PM Micha Silver <<a href="mailto:tsvibar@gmail.com">tsvibar@gmail.com</a>> wrote:<br>>><br>>> Hi Rich<br>>><br>>><br>>> On 09/24/2018 06:32 PM, Rich Shepard wrote:<br>>><br>>>   I want to determine whether GRASS or R is best suited to<br>>> interpolating/extrapolating annual mean precipitation data from 58 reporting<br>>> stations (unevenly distributed within the county) across the county. Some<br>>> flavor of kriging would be applied to these data to illustrate a general<br>>> long-term pattern of rainfall.<br>>><br>>>   Some years ago there was both a static display of a chemical constituent<br>>> in a river reach, a "heat map", (and an automation of temporal changes, if I<br>>> correctly recall) and I'm not finding this in the web site galleries.<br>>><br>>>   While elevation could be included as an explanatory variable using<br>>> regression kriging my purpose is to illustrate county-wide mean annual rainfall<br>>> distribution over a 13 year period, not to interpolate values for specific,<br>>> unsampled locations.<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 I<br>>> have the rgrass7 package working well; it imported the GRASS county boundary<br>>> map which I converted to a SpatialPolygonDataFrame.)<br>>><br>>><br>>> The guidelines that I follow include:<br>>><br>>> Rainfall interpolation (of any kind) should be done only for long time periods = at least a full season. since you are looking at 13 years of data then this requirement is fulfilled.<br>>> The rules of thumb for "how many points" for kriging interpolation usually says > 30. So again you are fine.<br>>> The GRASS modules offer only ordinary kriging AFAIK, which might be appropriate in this case. But if you want to use Kriging with External Drift with the elevation as the secondary "trend" variable, I think you'll need to go with with R. Having said that, I think that "the jury is still out" on whether elevation improves the interpolation or not. <br>>><br>>> Regards, Micha<br>>>  <br>>><br>>> Rich<br>>><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><br>>><br>>><br>>> -- <br>>> Micha Silver<br>>> Ben Gurion Univ.<br>>> Sde Boker, Remote Sensing Lab<br>>> cell: +972-523-665918<br>>><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><br>><br>><br>><br>> --<br>> Thomas E Adams, III<br>> 1724 Sage Lane<br>> Blacksburg, VA 24060<br>> <a href="mailto:tea3rd@gmail.com">tea3rd@gmail.com</a> (personal)<br>> <a href="mailto:tea@terrapredictions.org">tea@terrapredictions.org</a> (work)<br>><br>> 1 (513) 739-9512 (cell)<br>><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>