erase: Erase parts of a SpatialPolygons* or SpatialLines* object....

Description Usage Arguments Value Author(s) See Also Examples

Description

Erase parts of a SpatialPolygons* or SpatialLines* object with a SpatialPolygons* object

Usage

1
2
3
4
## S4 method for signature 'SpatialPolygons,SpatialPolygons'
erase(x, y, ...)
## S4 method for signature 'SpatialLines,SpatialPolygons'
erase(x, y, ...)

Arguments

x

SpatialPolygons or SpatialLines object

y

SpatialPolygons object

...

Additional arguments (none)

Value

Spatial*

Author(s)

Robert J. Hijmans

See Also

The equivalent for raster data is mask

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
if (require(rgdal) & require(rgeos)) {
# erase parts of polygons with other polygons
	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)
	
	
# erase parts of lines with polygons	
	r <- raster(extent(p) +c(-.1,.1,-.1,.1), crs=crs(p))
	start <- xyFromCell(r, cellFromCol(r, 1))
	end <- xyFromCell(r, cellFromCol(r, ncol(r)))
	lines <- do.call(spLines, lapply(1:10, function(i)rbind(start[i,], end[i,])))
	crs(lines) <- crs(p)
	
	e2 <- erase(lines, p)
	plot(p)
	lines(lines, col='blue', lwd=4, lty=3)
	lines(e2, col='red', lwd=2)
}

raster documentation built on Nov. 17, 2017, 5:51 a.m.