An HTML attachment was scrubbed... URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110317/d86a3873/attachment.html>
using IDW with anisotropy
4 messages · Edzer Pebesma, smurray at worldbank.org, Paul Hiemstra
Are you interpolating in 3D (x,y,t)?
On 03/17/2011 04:21 PM, smurray at worldbank.org wrote:
hello, I have a timseries dataset of monthly precip data from rain guages. The coverage is relatively sparse (approx 30 stations covering the country) and there is strong anisotropy in the data - at the monthly timestep data is strongly correlated in the E-W direction. Is there any way to use modify the IDW interpolation in gstat, to incorporate a directional bias in the distance weight? many thanks for any suggestions
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebesma at wwu.de
An HTML attachment was scrubbed... URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110317/dee42776/attachment.html> -------------- next part -------------- A non-text attachment was scrubbed... Name: graycol.gif Type: image/gif Size: 105 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110317/dee42776/attachment.gif> -------------- next part -------------- A non-text attachment was scrubbed... Name: pic08723.gif Type: image/gif Size: 1255 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110317/dee42776/attachment-0001.gif> -------------- next part -------------- A non-text attachment was scrubbed... Name: ecblank.gif Type: image/gif Size: 45 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110317/dee42776/attachment-0002.gif> -------------- next part -------------- A non-text attachment was scrubbed... Name: 2A748798.gif Type: image/gif Size: 821 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110317/dee42776/attachment-0003.gif>
Hi, You could also opt for kriging, which supports anisotropy. The intamap package has a function to detect anisotropy, in case you would like to automate the interpolation. cheers, Paul
On 03/18/2011 01:41 AM, smurray at worldbank.org wrote:
Thanks for the reply. I was planning to generate a surface for each time slice independently, so just x ,y. Siobhan Inactive hide details for Re: [R-sig-Geo] using IDW with anisotropyRe: [R-sig-Geo] using IDW with anisotropy *Re: [R-sig-Geo] using IDW with anisotropy* *Edzer Pebesma* to: r-sig-geo 03/17/2011 02:03 PM Sent by: *r-sig-geo-bounces at r-project.org* ------------------------------------------------------------------------ Are you interpolating in 3D (x,y,t)? On 03/17/2011 04:21 PM, smurray at worldbank.org wrote:
hello, I have a timseries dataset of monthly precip data from rain guages.
The coverage
is relatively sparse (approx 30 stations covering the country) and
there is
strong anisotropy in the data - at the monthly timestep data is
strongly
correlated in the E-W direction. Is there any way to use modify the IDW interpolation in gstat, to incorporate a directional bias in the
distance weight?
many thanks for any suggestions
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
-- Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebesma at wwu.de
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo _______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Paul Hiemstra, MSc Global Climate Division Royal Netherlands Meteorological Institute (KNMI) Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39 P.O. Box 201 | 3730 AE | De Bilt tel: +31 30 2206 494 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770 -------------- next part -------------- An HTML attachment was scrubbed... URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110318/c4b98c40/attachment.html> -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: image/gif Size: 105 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110318/c4b98c40/attachment.gif> -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: image/gif Size: 45 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110318/c4b98c40/attachment-0001.gif> -------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: image/gif Size: 821 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20110318/c4b98c40/attachment-0002.gif>