The Annals of Applied Statistics

Fibre-generated point processes and fields of orientations

Bryony J. Hill, Wilfrid S. Kendall, and Elke Thönnes

Full-text: Open access


This paper introduces a new approach to analyzing spatial point data clustered along or around a system of curves or “fibres.” Such data arise in catalogues of galaxy locations, recorded locations of earthquakes, aerial images of minefields and pore patterns on fingerprints. Finding the underlying curvilinear structure of these point-pattern data sets may not only facilitate a better understanding of how they arise but also aid reconstruction of missing data. We base the space of fibres on the set of integral lines of an orientation field. Using an empirical Bayes approach, we estimate the field of orientations from anisotropic features of the data. We then sample from the posterior distribution of fibres, exploring models with different numbers of clusters, fitting fibres to the clusters as we proceed. The Bayesian approach permits inference on various properties of the clusters and associated fibres, and the results perform well on a number of very different curvilinear structures.

Article information

Ann. Appl. Stat., Volume 6, Number 3 (2012), 994-1020.

First available in Project Euclid: 31 August 2012

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Markov chain Monte Carlo spatial birth–death process earthquakes empirical Bayes fibre processes field of orientations fingerprints spatial point processes tensor fields


Hill, Bryony J.; Kendall, Wilfrid S.; Thönnes, Elke. Fibre-generated point processes and fields of orientations. Ann. Appl. Stat. 6 (2012), no. 3, 994--1020. doi:10.1214/12-AOAS553.

Export citation


  • Adler, R. J. and Taylor, J. E. (2007). Random Fields and Geometry. Springer, New York.
  • Arsigny, V., Fillard, P., Pennec, X. and Ayache, N. (2006). Log-Euclidean metrics for fast and simple calculus on diffusion tensors. Magn. Reson. Med. 56 411–421.
  • August, J. and Zucker, S. W. (2003). Sketches with curvature: The curve indicator random field and Markov processes. IEEE Trans. Pattern Anal. Mach. Intell. 25 387–400.
  • Barrow, J. D., Bhavsar, S. P. and Sonoda, D. H. (1985). Minimal spanning trees, filaments and galaxy clustering. Royal Astronomical Society, Monthly Notices 216 17–35.
  • Brooks, S. P. and Roberts, G. O. (1998). Convergence assessment techniques for Markov chain Monte Carlo. Statist. Comput. 8 319–335.
  • Dryden, I. L., Koloydenko, A. and Zhou, D. (2009). Non-Euclidean statistics for covariance matrices, with applications to diffusion tensor imaging. Ann. Appl. Stat. 3 1102–1123.
  • Genovese, C. R., Perone-Pacifico, M., Verdinelli, I. and Wasserman, L. (2009). On the path density of a gradient field. Ann. Statist. 37 3236–3271.
  • Hill, B. J. (2011). An orientation field approach to modelling fibre-generated spatial point processes. Ph.D. thesis, Univ. Warwick. Available at
  • Illian, J., Penttinen, A., Stoyan, H. and Stoyan, D. (2008). Statistical Analysis and Modelling of Spatial Point Patterns. Wiley, New York.
  • Kaspi, H. and Mandelbaum, A. (1994). On Harris recurrence in continuous time. Math. Oper. Res. 19 211–222.
  • Le Bihan, D. L., Mangin, J. F., Poupon, C., Clark, C. A., Pappata, S., Molko, N. and Chabriat, H. (2001). Diffusion tensor imaging: Concepts and applications. J. Magn. Reson. Imaging 13 534–546.
  • Li, C., Sun, X., Zou, K., Yang, H., Huang, X., Wang, Y., Lui, S., Li, D., Zou, L. and Chen, H. (2007). Voxel based analysis of DTI in depression patients. International Journal of Magnetic Resonance Imaging 1 43–48.
  • Martínez, V. J. and Saar, E. (2002). Statistics of the Galaxy Distribution. Chapman & Hall/CRC, Boca Raton, FL.
  • Møller, J. and Waagepetersen, R. P. (2004). Statistical Inference and Simulation for Spatial Point Processes. Monographs on Statistics and Applied Probability 100. Chapman & Hall, London.
  • Pearl, J. (1988). Probabilistic Reasoning in Intelligent Systems: Networks of Plausible Inference. Morgan Kaufmann, San Mateo.
  • Stanford, D. C. and Raftery, A. E. (2000). Finding curvilinear features in spatial point patterns: Principal curve clustering with noise. IEEE Transactions on Pattern Analysis and Machine Intelligence 22 601–609.
  • Stoica, R. S., Martínez, V. J. and Saar, E. (2007). A three-dimensional object point process for detection of cosmic filaments. Appl. Statist. 56 459–477.
  • Stoica, R. S., Martínez, V. J. and Saar, E. (2010). Filaments in observed and mock galaxy catalogues. Astronomy and Astrophysics 510 1–12.
  • Stoica, R. S., Martínez, V. J., Mateu, J. and Saar, E. (2005). Detection of cosmic filaments using the Candy model. Astronomy and Astrophysics 434 423.
  • Stoyan, D., Kendall, W. S. and Mecke, J. (1995). Stochastic Geometry and Its Applications, 2nd ed. Wiley, New York.
  • Su, J. (2009). A tensor approach to fingerprint analysis. Ph.D. thesis, Univ. Warwick.
  • Su, J., Hill, B. J., Kendall, W. S. and Thönnes, E. (2008). Inference for point processes with unobserved one-dimensional reference structure. CRiSM Working Paper 8-10, Univ. Warwick.
  • Watson, C. (2001). NIST special database 30: Dual resolution images from paired fingerprint cards. National Institute of Standards and Technology, Gaithersburg.