The test with meuse.grid works:
#copied from file:SpatialGridDataFrame-class.html
coordinates(meuse.grid) = c("x", "y") # promote toSpatialPointsDataFrame
gridded(meuse.grid) <- TRUE # promote to SpatialPixelsDataFrame
x = as(meuse.grid, "SpatialGridDataFrame") # creates the full grid
#Now me:
>
I think that if you cast to SpatialPixels first, the problem will go away,
because the all-NA cases are dropped with their coordinates. Could you try the
same example with meuse.grid and see how it goes?
Roger
---
Roger Bivand, NHH, Helleveien 30, N-5045 Bergen, [EMAIL PROTECTED]
___
