[R-sig-Geo] universal kriging in gstat
Edzer J. Pebesma
e.pebesma at geo.uu.nl
Mon Jul 3 17:40:27 CEST 2006
Alí Santacruz wrote:
> Dear list members,
> I am trying to perform universal kriging in gstat, and according to the help
> of 'krige' function, the names of the coordinates of 'newdata' argument must
> be the same as defined in 'locations'.
I suspect it is time to update this documentation; If you're using sp
objects, the names of coordinates does not play a role anymore, as
we can retrieve them using the coordinates() method.
Does this help?
> I am creating an SpatialPixels object for the prediction locations, using
> the next code:
>> pts<-spsample(myShapePoly, n=3000, "regular")
>> gridded(pts) <- TRUE class(pts)
>  "SpatialPixels"
>  "sp"
> And x1 and x2 are the names of the coordinates of the prediction locations
> So, to perform universal kriging, I have to change the names of the
> coordinates of my observed data to x1 and x2, but I would like to know if I
> can change the name of the coordinates of the prediction grid ('pts').
> I appreciate any comment. Sincerely,
> Alí M. Santacruz
> M.Sc. Geomatics, Student
> National University of Colombia
> Bogotá D.C.
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
More information about the R-sig-Geo