Intersection of a Polygonal and a Circular Domain

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).

Currently the only supported class is "owin" from package spatstat.geom.

Usage

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

# S3 method for class 'owin'
intersectPolyCircle(object, center, radius, npoly = 32, ...)
<!-- %\method{intersectPolyCircle}{SpatialPolygons}(object, center, radius, npoly = 32, ...) -->

Arguments

object

a polygonal domain of one of the supported classes.

center,radius,npoly

see discpoly.

...

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

letterR <- surveillance:::LETTERR  # an "owin" (internally used for checks)
plot(letterR, axes = TRUE)
plot(intersectPolyCircle(letterR, center = c(-1, 2), radius = 2),
     add = TRUE, col = 4, lwd = 3)