[R-sig-Geo] Crop / clip a SpatialPixels* object

Mark Payne markpayneatwork at gmail.com
Sat Jun 22 15:24:58 CEST 2013


This should be simple, but I simply can't find an "elegant" solution.

I have a SpatialPixelsDataFrame that I would like to clip (crop) to a
smaller object. For example,

> data(meuse.grid)
> m = SpatialPixelsDataFrame(points = meuse.grid[c("x", "y")], data = meuse.grid)
> bbox(m)
     min    max
x 178440 181560
y 329600 333760

I would like to constrain m, to all y > 331000 and x >180000. What is
the most efficient way to do this (and have the other functionality
e.g bbox and plotting) reflect this reduced data set.


More information about the R-sig-Geo mailing list