Skip to contents

[DEFUNC] This is the same as "sector pair correlation function" in Stoyan 1991. Direction defined with a unit vector. Estimated in a cone of angle theta and different radii.

Usage

pcf_sector(
  x,
  r,
  theta,
  unit,
  epsilon,
  h,
  stoyan = 0.15,
  correction = "translation",
  n_dir = 12
)

Arguments

x

pp, list with $x~coordinates $bbox~bounding box

r

Ranges

theta

Direction sector central angles, list of ncol(x$x)-components.

unit

Alternative to theta, a matrix of direction vectors.

epsilon

Direction sector width angles

h

width of epanechnicov kernel

stoyan

if h not given, use h=stoyan/sqrt(lambda)

correction

"none" or translation. Translation only for rectangle box.

n_dir

If unit missing, compute on a grid with resolution n_dir

Details

Otherwise same as default pcf in spatstat, but now with an angle.