plot.owin {spatstat}R Documentation

Plot a Spatial Window

Description

Plot a two-dimensional window of observation for a spatial point pattern

Usage

 plot.owin(x, main, add=FALSE, ..., box=TRUE, edge=0.04)

Arguments

x The window to be plotted. An object of class owin, or data which can be converted into this format by as.owin().
main text to be displayed as a title above the plot.
add logical flag: if TRUE, draw the window in the current plot; if FALSE, generate a new plot.
... extra arguments passed to the generic plot function.
box logical flag; if TRUE, plot the enclosing rectangular box
edge nonnegative number; the plotting region will have coordinate limits that are 1 + edge times as large as the limits of the rectangular box that encloses the pattern.

Details

This is the plot method for the class owin. The action is to plot the boundary of the window on the current plot device, using equal scales on the x and y axes.

If the window x is of type "rectangle" or "polygonal", the boundary of the window is plotted as a polygon or series of polygons. If x is of type "mask" the discrete raster approximation of the window is displayed as a binary image (white inside the window, black outside).

Graphical parameters controlling the display (e.g. setting the colours) may be passed directly via the ... arguments, or indirectly reset using spatstat.options. The parameter par.binary of spatstat.options controls the display when x is of type "mask".

Value

none.

Author(s)

Adrian Baddeley adrian@maths.uwa.edu.au http://www.maths.uwa.edu.au/~adrian/ and Rolf Turner rolf@math.unb.ca http://www.math.unb.ca/~rolf

See Also

owin.object, plot.ppp, spatstat.options

Examples

## Not run: 
  # rectangular window
   data(nztrees)
   plot(nztrees$window)
   abline(v=148, lty=2)

  # polygonal window
  plot(c(0,1),c(0,1),type="n")
  bdry <- locator()
  # click the vertices of a polygon (anticlockwise)
  w <- owin(c(0,1), c(0,1), poly=bdry)
  plot(w)

  # binary mask
  we <- erode.owin(w, 0.05, FALSE)
  plot(we)
  spatstat.options(par.binary=list(col=grey(c(0.5,1))))
  plot(we)
## End(Not run)

[Package spatstat version 1.8-8 Index]