Hi all, I'm trying to figure out why coercing a SpatialGrid directly to a SpatialPoints doesn't result in the same thing as coercing it via an intermediate SpatialPixels. Here's a minimal example:
require(sp) grd <- SpatialGrid(GridTopology(c(0,0), c(10,10), c(4,4))) as(as(grd, "SpatialPixels"), "SpatialPoints") # result == what I expect as(grd, "SpatialPoints") # result != expected Is this behavior as expected? Maybe I am doing it wrong? Many thanks, David -- David Holstius PhD Student in Environmental Health Sciences UC Berkeley School of Public Health _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
