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@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.pebe...@wwu.de
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
