Dear sig-geos,
is there any way of converting from vector
(imported with rgdal from a shp) to raster
within R?
Thanks!
Agus
Yes, using the classes and methods in sp -- you'd first have to create a
grid e.g. with spsample, sampling type regular or makegrid. Next, if you
want to have the polygon values at grid cell centres, overlay the
polygons with the grid, using the overlay method.
> library(sp)
> data(meuse.riv)
> meuse.sr = SpatialPolygons(list(Polygons(list(Polygon(meuse.riv)),
>
> plot(meuse.sr)
> points(spsample(meuse.sr, n = 1000, "regular"), pch = 3)
> class(meuse.sr)
[1] "SpatialPolygons"
attr(,"package")
[1] "sp"
> pts = spsample(meuse.sr, n = 1000, "regular")
> class(pts)
[1] "SpatialPoints"
attr(,"package")
[1] "sp"
> gridded(pts)=T
> class(pts)
[1] "SpatialPixels"
attr(,"package")
[1] "sp"
--
Edzer