K-Functions on Fibers
Joint work with
Hans Stephensen, Jon Sporring and Stefan sommer (DIKU)
Anne-Marie Svane and Rasmus Waagepetersen (AAU)
Stine Hasselholt (AU)

Content of talk
- Introduce point processes and Ripley's K-function
- Shape valued processes
- Currents metric and fiber K-function
- Result on myelin fibers



Point processes
A point process X on a metric space S is a measurable map from some probability space (Ω,F,P) into the space of locally finite subsets of S.
Usually S=Rn

Ripley's K-function for point processes X on Rn is defined by
K(r)=λ1 E[ #events within distance r from point]
where λ is the density (number per unit area) of events.

Ripley’s K-function is a well-known statistical tool for
- analyzing second order moment
structure of point processes - measuring deviance from complete
spatial randomness in point sets.
Ripley's K-function
Ripley’s K-function is a well-known statistical tool for
- analyzing second order moment
structure of point processes - measuring deviance from complete
spatial randomness in point sets.
Ripley's K-function

r
K(r)

Goal
Define a k-function for point processes where every point is a shape, for example fibers in R3
Issues
- Distance measure
- Well posed K-function
- Uniform shape process
Space of Currents
A piece-wise smooth curve γ is characterized by the path-integral over vector fields w
where τ(x) is the tangent of γ at x.
Vγ(w)=∫γw(x)tτ(x)dλ(x)

Thus, we can think of γ as a functional
Vγ∈W∗={f:W→R∣linear and cont}.
This representation is called a current and formally the space of m-currents is the dual space of the differential m-forms.
By considering W as a RKHS, we get the inner product on W∗
⟨Vγ1,Vγ2⟩=∫γ1∫γ2τγ1(x)tK(x,y)τγ2(x)dλ(x)λ(y)
where K is the kernel of W and the induced distance measure
dc(Vγ1,Vγ2)=∣∣Vγ1−Vγ2∣∣
Space of Currents
The K-function for fiber processes X observed in window W is
where c(.) denotes center point of curve, ν(S0) measure on centered curves and ∣W∣ size of window.

The K-function
- s shape parameter
- t spatial parameter
K(s,t)=∣W∣ν(S0)1 ∑ ∑ 1[∣∣c(γ)−c(γ′)∣∣≤t,dc(γ,γ′)≤s]
γ∈X:c(γ)∈W
γ′=γ∈X
Uniform shape process






Uniformly
rotated lines:
Uniformly
rotated spirals:
Brownian motions:

Myelin fibers
Future Work
- Improve uniform shape process comparison
- Review benefits from using currents metric
- Steel fibers in cement
Thank you for watching!
K-functions for Currents
By pernilleehh
K-functions for Currents
- 408