erase {raster} | R Documentation |
Erase parts of a Spatial* objects with another Spatial* object
## S4 method for signature 'SpatialPolygons,SpatialPolygons' erase(x, y, ...)
x |
Spatial* object |
y |
Spatial* object |
... |
Additional arguments (none) |
Spatial*
Robert J. Hijmans
if (require(rgdal) & require(rgeos)) { p <- shapefile(system.file("external/lux.shp", package="raster")) b <- as(extent(6, 6.4, 49.75, 50), 'SpatialPolygons') projection(b) <- projection(p) e <- erase(p, b) plot(e) }