Skip to content
Prev 9769 / 29559 Next

convert from spatstat to sp

There is; if your question is about converting an object of class im
(spatstat) to an object of class SpatialGridDataFrame (sp), it is in
package maptools:

library(spatstat)
whitenoise <- im(matrix(rnorm(1200),30,40))
library(maptools)
x = as(as(whitenoise, "SpatialGridDataFrame"), "SpatialPixelsDataFrame")

if your polygons are in SpatialPolygons object pol, then

sel = overlay(x, pol)
x.sel = x[!is.na(sel),]
spplot(x.sel)

should get you there.
On 10/29/2010 09:28 PM, Peter Larson wrote: