[R-sig-Geo] Crop / clip a SpatialPixels* object
Mark Payne
markpayneatwork at gmail.com
Sat Jun 22 15:24:58 CEST 2013
Hi,
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.
Mark
More information about the R-sig-Geo
mailing list