1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
\name{union}
\docType{methods}
\alias{union}
\alias{union,Extent,Extent-method}
\alias{union,SpatialPolygons,SpatialPolygons-method}
\alias{union,SpatialPolygons,missing-method}
\alias{union,SpatialPoints,SpatialPoints-method}
\alias{union,SpatialLines,SpatialLines-method}
\title{
Union Extent or SpatialPolygons* objects
}
\description{
Extent objects: Objects are combined into their union. See \code{\link{crop}} and \code{\link{extend}} to union a Raster object with an Extent object.
Two SpatialPolygons* objects. Overlapping polygons (between layers, not within layers) are intersected, other spatial objects are appended. Tabular attributes are joined. See \code{\link{bind}} if you want to combine polygons without intersection.
Single SpatialPolygons* object. Overlapping polygons are intersected. Original attributes are lost. New attributes allow for determining how many, and which, polygons overlapped.
Union for SpatialLines and SpatialPoints simply combines the two data sets; without any geometric intersections. This is equivalent to \code{\link{bind}}.
}
\usage{
\S4method{union}{Extent,Extent}(x, y)
\S4method{union}{SpatialPolygons,SpatialPolygons}(x, y)
\S4method{union}{SpatialPolygons,missing}(x, y)
\S4method{union}{SpatialLines,SpatialLines}(x, y)
\S4method{union}{SpatialPoints,SpatialPoints}(x, y)
}
\arguments{
\item{x}{Extent or SpatialPolygons* object}
\item{y}{Same as \code{x} or missing}
}
\value{
Extent or SpatialPolygons object
}
\seealso{
\code{\link[raster]{intersect}, \link[raster]{extent}, \link[raster]{setExtent}}
\code{\link[sp]{merge}} for merging a data.frame with attributes of Spatial objects
and \code{\link{+,SpatialPolygons,SpatialPolygons-method}} for an algebraic notation
}
\examples{
e1 <- extent(-10, 10, -20, 20)
e2 <- extent(0, 20, -40, 5)
union(e1, e2)
#SpatialPolygons
p <- shapefile(system.file("external/lux.shp", package="raster"))
p0 <- aggregate(p)
b <- as(extent(6, 6.4, 49.75, 50), 'SpatialPolygons')
crs(b) <- crs(p)
u <- union(p0, b)
plot(u, col=2:4)
}
\keyword{methods}
\keyword{spatial}
|
