raster (version 2.5-8)

erase: Erase parts of a Spatial* object

Description

Erase parts of a Spatial* objects with another Spatial* object

Usage

# S4 method for SpatialPolygons,SpatialPolygons
erase(x, y, ...)

Arguments

x

Spatial* object

y

Spatial* object

...

Additional arguments (none)

Value

Spatial*

Examples

# NOT RUN {
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)
}
# }

Run the code above in your browser using DataCamp Workspace