[R-sig-Geo] universal kriging in gstat

Alí Santacruz amsantac at hotmail.com
Mon Jul 3 18:07:11 CEST 2006

Thanks Edzer,

Yes, I found that the answer was very simple: as you say, we can use the 
coordinates() method.

Best regards,

Alí Santacruz

>From: "Edzer J. Pebesma" <e.pebesma at geo.uu.nl>
>To: Alí Santacruz <amsantac at hotmail.com>
>CC: r-sig-geo at stat.math.ethz.ch
>Subject: Re: [R-sig-Geo] universal kriging in gstat
>Date: Mon, 03 Jul 2006 17:40:27 +0200
>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)
>>[1] "SpatialPixels"
>>[1] "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

MSN Amor: busca tu ½ naranja http://latam.msn.com/amor/

More information about the R-sig-Geo mailing list