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

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: https://gis.stackexchange.com/questions/46954/clip-spatial-object-to-bounding-box-in-r/.

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

See also

Other geo: bbox_scale(), geo_bb_matrix(), geo_bb(), mapshape_available(), mapshape(), quadrant(), reproject()

Examples

if (FALSE) {
cb <- sf::st_buffer(cents_sf[8, ], dist = 0.012)
plot(cents_sf$geometry)
plot(cb, add = TRUE)
clipped <- gclip(cents_sf, cb)
plot(clipped, add = TRUE)
clipped$avslope # gclip also returns the data attribute
points(clipped)
points(cents[cb, ], col = "red") # note difference
gclip(cents_sf, cb)
}