spatstat (version 1.23-1)

density.psp: Kernel Smoothing of Line Segment Pattern

Description

Compute a kernel smoothed intensity function from a line segment pattern.

Usage

## S3 method for class 'psp':
density(x, sigma, \dots, edge=TRUE)

Arguments

x

Line segment pattern (object of class "psp") to be smoothed.

sigma

Standard deviation of isotropic Gaussian smoothing kernel.

...

Extra arguments passed to as.mask which determine the resolution of the resulting image.

edge

Logical flag indicating whether to apply edge correction.

Value

• A pixel image (object of class "im").

Details

This is a method for the generic function density.

A kernel estimate of the intensity of the line segment pattern is computed. The result is the convolution of the isotropic Gaussian kernel, of standard deviation sigma, with the line segments. Computation is performed analytically. If edge=TRUE this result is adjusted for edge effects by dividing it by the convolution of the same Gaussian kernel with the observation window.

See Also

psp.object, im.object, density

Examples

L <- psp(runif(20),runif(20),runif(20),runif(20), window=owin())
  D <- density(L, sigma=0.03)
  plot(D, main="density(L)")
  plot(L, add=TRUE)