Skip to content
Prev 12753 / 29559 Next

grid overlay

Hi Josh,

Am 09.09.2011 18:16, schrieb Josh Tewksbury:
Have a look at ?expand.grid to see how you can get that easier.
Yes it is:

library(sp)
coordinates(gridA) <- 1:2; gridded(gridA) <- TRUE; fullgrid(gridA) <- TRUE
coordinates(gridB) <- 1:2; gridded(gridB) <- TRUE; fullgrid(gridB) <- TRUE
gridA <- SpatialGridDataFrame(gridA, over(gridA, gridB))

But there might be more appropriate means in package raster.

Regards,
Tom