Draw the Stienen diagram of a point pattern, or compute the region covered by the Stienen diagram.
stienen(X, …, bg = "grey", border = list(bg = NULL))
stienenSet(X, edge=TRUE)Point pattern (object of class "ppp").
Arguments passed to plot.ppp
to control the plot.
Fill colour for circles.
Either a list of arguments passed to plot.ppp
to control the display of circles at the border of the diagram,
or the value FALSE
indicating that the border circles should not be plotted.
Logical value indicating whether to include the circles at the border of the diagram.
The plotting function stienen returns NULL.
The return value of stienenSet is a window (object of class
"owin").
The Stienen diagram of a point pattern (Stienen, 1982) is formed by drawing a circle around each point of the pattern, with diameter equal to the nearest-neighbour distance for that point. These circles do not overlap. If two points are nearest neighbours of each other, then the corresponding circles touch.
stienenSet(X) computes the union of these circles and
returns it as a window (object of class "owin").
stienen(X) generates a plot of the Stienen diagram of
the point pattern X. By default, circles are shaded in grey
if they lie inside the window of X, and are not shaded
otherwise.
Stienen, H. (1982) Die Vergroeberung von Karbiden in reinen Eisen-Kohlenstoff Staehlen. Dissertation, RWTH Aachen.
# NOT RUN {
Y <- stienenSet(cells)
stienen(redwood)
stienen(redwood, border=list(bg=NULL, lwd=2, cols="red"))
# }
Run the code above in your browser using DataLab