Dear r-sig-geo, In the course of my PhD I have written a R package that performs automatic interpolation (kriging) and might be interesting for other people. I work on real-time automatic interpolation of radioactivity levels. It depends heavily on the gstat pacakge written by Edzer (Pebesma, he is my supervisor). If people are interested in this package, it can be downloaded from http://intamap.geo.uu.nl/~paul/Downloads.html. If you have any comments, complaints or suggestions please do not hesitate to e-mail me.
Some sample code: library(automap) data(meuse) coordinates(meuse) = ~x+y data(meuse.grid) gridded(meuse.grid) = ~x+y kr.ok = autoKrige(log(zinc)~1, meuse, meuse.grid) plot(kr.ok) kr.uk = autoKrige(log(zinc)~sqrt(dist), meuse, meuse.grid) plot(kr.uk) cheers, Paul -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +31302535773 Fax: +31302531145 http://intamap.geo.uu.nl/~paul _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
