spatstat.geom (version 3.3-6)

crossdist.pp3: Pairwise distances between two different three-dimensional point patterns

Description

Computes the distances between pairs of points taken from two different three-dimensional point patterns.

Usage

# S3 method for pp3
crossdist(X, Y, ..., periodic=FALSE, squared=FALSE)

Value

A matrix whose [i,j] entry is the distance from the i-th point in X

to the j-th point in Y.

Arguments

X,Y

Point patterns in three dimensions (objects of class "pp3").

...

Ignored.

periodic

Logical. Specifies whether to apply a periodic edge correction.

squared

Logical. If squared=TRUE, the squared distances are returned instead (this computation is faster).

Author

Adrian Baddeley based on code for two dimensions by Pavel Grabarnik .

Details

Given two point patterns in three-dimensional space, this function computes the Euclidean distance from each point in the first pattern to each point in the second pattern, and returns a matrix containing these distances.

This is a method for the generic function crossdist for three-dimensional point patterns (objects of class "pp3").

This function expects two point patterns X and Y, and returns the matrix whose [i,j] entry is the distance from X[i] to Y[j].

Alternatively if periodic=TRUE, then provided the windows containing X and Y are identical and are rectangular, then the distances will be computed in the `periodic' sense (also known as `torus' distance): opposite edges of the rectangle are regarded as equivalent. This is meaningless if the window is not a rectangle.

See Also

crossdist, pairdist, nndist, G3est

Examples

if(require(spatstat.random)) {
   X <- runifpoint3(20)
   Y <- runifpoint3(30)
} else {
   X <- osteo$pts[[1]]
   Y <- osteo$pts[[2]]
   Y <- Y[domain(X)]
}
   d <- crossdist(X, Y)
   d <- crossdist(X, Y, periodic=TRUE)