This is a unifying wrapper around functionality of various packages dealing with spatial data. It computes the intersection of a circular domain and a polygonal domain (whose class defines the specific method).

intersectPolyCircle(object, center, radius, ...)

# S3 method for owin
intersectPolyCircle(object, center, radius, npoly = 32, ...)
# S3 method for SpatialPolygons
intersectPolyCircle(object, center, radius, npoly = 32, ...)
# S3 method for gpc.poly
intersectPolyCircle(object, center, radius, npoly = 32,
                   useGEOS = FALSE, ...)

Arguments

object

a polygonal domain of one of the supported classes.

center,radius,npoly

see discpoly.

useGEOS

logical indicating if package rgeos (gIntersection) should be used instead of package gpclib. The latter (default) requires explicit acceptance of gpclib's restricted license via surveillance.options(gpclib=TRUE).

...

potential further arguments (from the generic).

Value

a polygonal domain of the same class as the input object.

Author

Sebastian Meyer

See also

discpoly to generate a polygonal approximation to a disc

Examples

library("spatstat.geom")
plot(letterR)
plot(intersectPolyCircle(letterR, c(3,2), 1), add=TRUE, col=2, lwd=3)