Compute the cumulative integral of an image over increasing radial distances from the origin.
radcumint(X, ..., origin, Xname, result = c("fv", "im"))An object of class "fv" or "im",
with the same coordinate units as X.
A pixel image (object of class "im")
with numerical or logical values.
Ignored.
Optional. Origin about which the rotations should be performed.
Either a numeric vector or a character string as described
in the help for shift.owin.
Optional name for X to be used in the
function labels.
Character string specifying the kind of result required: either a function object or a pixel image.
Adrian Baddeley Adrian.Baddeley@curtin.edu.au, Rolf Turner r.turner@auckland.ac.nz and Ege Rubak rubak@math.aau.dk.
This command computes, for each possible distance \(r\), the integral of the pixel values lying inside the disc of radius \(r\) centred at the origin.
If result="fv" (the default) the result is a function
object f of class "fv". For each value of radius \(r\),
the function value f(r) is the integral of X
over the disc of radius \(r\).
If result="im" the result is a pixel image, with the same
dimensions as X. At a given pixel, the result is
equal to f(r) where r is the distance from the given
pixel to the origin. That is, at any given pixel, the resulting value
is the integral of X over the disc
centred at the origin whose boundary passes through the given pixel.
rotmean, spatialcdf
D <- density(redwood)
plot(radcumint(D))
plot(radcumint(D, result="im"))
Run the code above in your browser using DataLab