spatstat.geom (version 3.6-1)

pixellate.psp: Convert Line Segment Pattern to Pixel Image

Description

Converts a line segment pattern to a pixel image by measuring the length or number of lines intersecting each pixel.

Usage

# S3 method for psp
pixellate(x, W=NULL, ..., weights = NULL,
                        what=c("length", "number", "indicator"),
                        DivideByPixelArea=FALSE)

Value

A pixel image (object of class "im") with numeric values.

Details

This function converts a line segment pattern to a pixel image by computing, for each pixel, the total length of intersection between the pixel and the line segments. Alternatively it can count the number of line segments intersecting each pixel.

This is a method for the generic function pixellate for the class of line segment patterns.

The pixel raster is determined by W and the optional arguments .... If W is missing or NULL, it defaults to the window containing x. Then W is converted to a binary pixel mask using as.mask. The arguments ... are passed to as.mask to control the pixel resolution.

If weights are given and what="length", then the length of the intersection between line segment i and pixel j is multiplied by weights[i] before the lengths are summed for each pixel.

See Also

pixellate, as.mask, psp2mask.

Use psp2mask if you only want to know which pixels are intersected by lines.

Examples

  X <- psp(runif(10),runif(10), runif(10), runif(10), window=owin())
  plot(pixellate(X))
  plot(X, add=TRUE)
  sum(lengths_psp(X))
  sum(pixellate(X))
  plot(pixellate(X, what="n"))