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

[Lionel Hertzog]

In this example it is pretty straightforward using the subset function:

m1<-subset(m,y>331000 & x>180000)
bbox(m1)

See this link for another way to do this using rgeos:
http://stackoverflow.com/questions/13982773/crop-for-spatialpolygonsdataframe

Cheers,
Lionel


On 22/06/2013 15:24, Mark Payne wrote:
> 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
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo