[R-sig-Geo] error in krig interpolation
Edzer J. Pebesma
e.pebesma at geo.uu.nl
Thu Apr 5 14:30:33 CEST 2007
el1 is already a SpatialPointsDataFrame, so omit the locations argument.
You may try something as simple as
zz <- krige(log(z) ~ 1, el1, mask_SG, efitted)
--
Edzer
temiz wrote:
> when does this error occur ?
>
> > zz <- krige(formula = log(el1$z) ~ 1, locations = ~x+y , data = el1
> ,newdata = mask_SG , model = efitted) Error in function (classes, fdef,
> mtable) :
> unable to find an inherited method for function "coordinates<-",
> for signature "SpatialPointsDataFrame"
>
>
> str(el1) gives:
>
> Formal class 'SpatialPointsDataFrame' [package "sp"] with 5 slots
> ..@ data :'data.frame': 24051 obs. of 4 variables:
> .. ..$ cat: num [1:24051] 1 2 3 4 5 6 7 8 9 10 ...
> .. ..$ x : num [1:24051] 449114 449204 449262 449263 449276 ...
> .. ..$ y : num [1:24051] 4590406 4590375 4590117 4590304 4590214 ...
> .. ..$ z : num [1:24051] 110 110 110 110 110 110 110 110 110 110 ...
> ..@ coords.nrs : num(0)
> ..@ coords : num [1:24051, 1:2] 449114 449204 449262 449263 449276 ...
> .. ..- attr(*, "dimnames")=List of 2
> .. .. ..$ : NULL
> .. .. ..$ : chr [1:2] "coords.x1" "coords.x2"
> ..@ bbox : num [1:2, 1:2] 449083 4584531 456552 4591890
> .. ..- attr(*, "dimnames")=List of 2
> .. .. ..$ : chr [1:2] "coords.x1" "coords.x2"
> .. .. ..$ : chr [1:2] "min" "max"
> ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots
> .. .. ..@ projargs: chr " +proj=utm +zone=36 +ellps=intl +units=m
> +no_defs"
>
>
>
> regards
>
>
More information about the R-sig-Geo
mailing list