[R-sig-Geo] error message in function kriege

[Edzer Pebesma]

Maybe you could try:

 rich.kriged = krige(N~1, rich, as(r, "SpatialGridDataFrame"),
    model = rich.fit)

and, if needed, convert rich.kriged back to raster.

On 11/05/2010 07:17 PM, Manuel Spínola wrote:
> Dear list members,
> 
> I am doing a kriging and I got this error message.
> 
>  > rich.kriged = krige(N~1, rich, r, model = rich.fit)
> Error en proj4string(newdata) :
>    proj4string only works for class(es extending) Spatial
> 
> "rich" is SpatialPointsDataFrame
> "r" is a projected raster
> 
> Any help will be appreciatedr?
> 
> Best,
> 
> Manuel
> 
> 
> 
> 
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics      e.pebesma at wwu.de