translational_correction: Translational edge-correction

In myllym/marksummary: Calculate summary functions for marked point patterns

Description

Calculates a vector of edge correction weights for given point pairs using the translational edge-correction. Notice that the inverse has already been taken in the returned vector. Each matrix element e_1, 2^-1 has value 1.0 / ((win_x - abs(x1 - x2)) * (win_y - abs(y1 - y2))) for a rectangular window. For a circular window, it is 1.0 / ( 2.0 * R^2 * acos(d/(2*R)) - d/2 * sqrt(4*R^2 - d^2) ) where d is the distance between (x1, y1) and (x2, y2) and R is the radius of the disc.

Usage

translational_correction(window, x1, y1, x2, y2)

Arguments

window	A window object from owin where the points reside. The window has to have type "rectangle".
x1	A vector containing the x-coordinates of the first points of the pairs. The order of the elements has to correspond with y1.
y1	A vector containing the y-coordinates of the first points of the pairs. The order of the elements has to correspond with x1.
x2	A vector containing the x-coordinates of the second points of the pairs. The order of the elements has to correspond with y2.
y2	A vector containing the y-coordinates of the second points of the pairs. The order of the elements has to correspond with x2.

Value

A vector containing the edge correction coefficients. Notice that the inverse has already been taken. The order of the pairs is the same as in ((x1, y1), (x2, y2)).

References

[1] J. Illian, A. Penttinen, H. Stoyan, and D. Stoyan, Statistical Analysis and Modelling of Spatial Point Patterns, 1st ed. p. 188. John Wiley & Sons, Ltd, 2008.

myllym/marksummary documentation built on May 23, 2019, 11:59 a.m.