stplanr (version 0.0.2)

gclip: Crops spatial object x to the bounding box of spatial object (or matrix) b

Description

This function is a cross between the spatial subsetting funtions such as sp::over(), rgeos::gIntersects() etc, and the cropping functions of raster::crop() and rgeos::gIntersection(). The output is the subset of spatial object a with an outline described by a square bounding box. The utility of such a function is illustrated in the following question: http://gis.stackexchange.com/questions/46954/clip-spatial-object-to-bounding-box-in-r/.

Usage

gclip(shp, bb)

Arguments

shp

The spatial object a to be cropped

bb

the bounding box or spatial object that will be used to crop shp

Examples

library(sp)
data(cents)
bb <- bbox(cents)
cb <- rgeos::gBuffer(cents[8, ], width = 0.012, byid = TRUE)
plot(cents)
plot(cb, add = TRUE)
clipped <- gclip(cents, cb)
row.names(clipped)
clipped$avslope # gclip also returns the data attribute
points(clipped)
points(cents[cb,], col = "red") # note difference