[R-sig-Geo] geoR help with kriging (newbie!)

Mike Smith mik-smith at supanet.com
Thu Jun 9 14:44:43 CEST 2005


Hi

I have recently been experimenting with the use of kriging, primarily
through Goldensoftware's Surfer. Having had a look around at the
packages available for R I see the geoR (amongst several) supports
kriging and so have started experimenting with this.

What I am interested in doing is batch processing ~300 datasets using
ordinary kriging with an exact interpolator (which I believe has no nugget).

As I am not expert with either R or kriging (!!) I was wondering
what the best command line options might be? I have been working
through the geoR intro and have a few queries:

-the command below appears to cover everything I want to do, however do
I need to fit a model to the variogram first and, if so, where is this
specified (using krige?)
kc <- krige.conv(s100, loc = pred.grid, krige = krige.control(obj.m = ml))

-the input file is a data frame with three columns; x and y coordinates
and a z value to interpolate. I would ideally like to tell geoR (using
loc) to take the max extent of my area and, at a specified point spacing,
interpolate a new grid.


Many thanks

mike


-- 
Dr Mike J Smith
Editor-in-Chief: Journal of Maps
http://www.journalofmaps.com

School of Earth Sciences and Geography
Kingston University
Penrhyn Road
Kingston-upon-Thames
Surrey
KT1 2EE UK
tel +44 (0)208 547 7500
fax +44 (0)208 547 7497
OS X/Y:518098/168540
Lat/Long: 51.403603/-0.303575

Personal Webpage:
http://www.kingston.ac.uk/esg/staff/smith.htm

-----BEGIN PGP PUBLIC KEY BLOCK-----
Version: 2.6

mQCPA0AyMQwBbgEEAMDTRNYB5Zi2SbmOtnBULHWIMXFaK+MVSTqGPi66N8s0ZQaJ
bMggYe5ySdltKLgT1FFmJN48IdSzqvzwPAMD5pQEuyLadnBfFvdfJapBcvHTd1pL
h5ikc/3zUXnK/JegioA4bCPypGnhf/1HFx70GRhpSaEOiul+vK1pQRliQJz1ABEB
AAG0KU1pa2UgU21pdGggPG1pY2hhZWwuc21pdGhAa2luZ3N0b24uYWMudWs+iQCV
AwUQQDIxDK1pQRliQJz1AQHn1QQAqA+203XQM7t0GEhIUh4PFEMysT+s2jMGUkNm
eihQnyK+1ljBcG9Vt8DtzOC5s3Pop1G9nbW9vl7OYzSb4E5uQ5L5wFUz0TJ81mu6
V8u/1CsERi+eXcNre+MON+q1GYoerUSd6LTh/kj+4J+MU05zHpSMXVJQsr2wcaxJ
vDM8DHY=
=oIrO
-----END PGP PUBLIC KEY BLOCK-----




More information about the R-sig-Geo mailing list