Abstract
An attribute opening is an idempotent, anti-extensive and increasing operator, which removes from an image connected components which do not fulfil a given criterion. When the increasingness property is dropped, we obtain a—more general—attribute thinning. In this paper, we propose efficient grey scale thinnings based on geodesic attributes.
Given that the geodesic diameter is time consuming, we propose a new geodesic attribute, the barycentric diameter to speed up the computation time. Then, we give the theoretical error bound between these two attributes, and we note that in practice, the barycentric diameter gives very similar results in comparison with the geodesic diameter. Finally, we present the algorithm with further optimisations, to obtain a 60× speed up.
We illustrate the use of these thinnings in automated non-destructive material inspection: the detection of cracks. We discuss the advantages of these operators over other methods such as path openings or the supremum of openings with segments.
Similar content being viewed by others
Notes
The normalisation factors are such that in a two-dimensional Euclidian space, the geodesic elongation of a disk is equal to 1.
Database available in www.imageprocessingplace.com/root_files_V3/image_databases.htm.
A demonstration version is available: http://cmm.ensmp.fr/~morard/DemoGeoThinnings.html.
Image from the Institute for Molecular Virology. University of Wisconsin—Madison: http://www.biochem.wisc.edu/faculty/inman/empics/dna-prot.htm.
We made an experiment on more than 100 natural images showing that a pixel is processed less than two times in average.
References
Beucher, S.: Numerical residues. Image Vis. Comput. 25(4), 405–415 (2007)
Beucher, S., Blosseville, J.M., Lenoir, F.: Traffic spatial measurements using video image processing. Intelligent Robots and Computer Vision. Proc. SPIE 848, 648–655 (1987)
Braga-Neto, U.M.: Alternating sequential filters by adaptive-neighborhood structuring functions. In: Mathematical Morphology and Its Applications to Image and Signal Processing, vol. 5, pp. 139–146. Maragos, Schafer and Butt, Atlanta (1996)
Breen, E.J., Jones, R.: Attribute openings, thinnings, and granulometries. Comput. Vis. Image Underst. 64(3), 377–389 (1996)
Grazzini, J., Soille, P.: Improved morphological interpolation of elevation contour data with generalised geodesic propagations. Lect. Notes Comput. Sci. 4673, 742–750 (2007). doi:10.1007/978-3-540-74272-2_92
Heijmans, H., Buckley, M., Talbot, H.: Path openings and closings. J. Math. Imaging Vis. 22(2), 107–119 (2005)
Hendriks, C.L.L.: Constrained and dimensionality-independent path openings. IEEE Trans. Image Process. 19(6), 1587–1595 (2010)
Hernández, J., Marcotegui, B.: Ultimate attribute opening segmentation with shape information. In: Mathematical Morphology and Its Application to Signal and Image Processing, pp. 205–214 (2009)
Lantuejoul, C., Beucher, S.: On the use of the geodesic metric in image analysis. J. Microsc. 121(1), 39–49 (1981)
Lantuéjoul, C., Maisonneuve, F.: Geodesic methods in quantitative image analysis. Pattern Recognit. 17(2), 177–187 (1984)
Maragos, P.: Pattern spectrum and multiscale shape representation. IEEE Trans. Pattern Anal. Mach. Intell. 11(7), 701–716 (1989)
Maragos, P., Ziff, R.D.: Threshold superposition in morphological image analysis systems. IEEE Trans. Pattern Anal. Mach. Intell. 12(5), 498–504 (1990)
Matheron, G.: Random Sets and Integral Geometry, vol. 1. Wiley, New York (1975)
Morard, V., Decencière, E., Dokládal, P.: Geodesic attributes thinnings and thickenings. In: Mathematical Morphology and Its Applications to Image and Signal Processing, vol. 6671, pp. 200–211. Springer, Berlin (2011)
Morard, V., Decencière, E., Dokládal, P.: Region growing structuring elements and new operators based on their shape. In: Signal and Image Processing, pp. 78–85. ACTA Press, Calgary (2011)
Najman, L., Couprie, M.: Building the component tree in quasi-linear time. IEEE Trans. Image Process. 15(11), 3531–3539 (2006)
Ngan, N., Dokládalová, E., Akil, M., Contou-Carrère, F.: Fast and efficient FPGA implementation of connected operators. J. Syst. Archit. 57(8), 778–789 (2011)
Nilsson, S.: Radix sorting and searching. Ph.D. thesis, Department of Computer Science, Lund University, Lund, Sweden (1996)
Ouzounis, G.K., Wilkinson, M.H.F.: Mask-based second-generation connectivity and attribute filters. IEEE Trans. Pattern Anal. Mach. Intell. 29(6), 990–1004 (2007)
Ouzounis, G.K., Wilkinson, M.H.F.: Hyperconnected attribute filters based on k-flat zones. IEEE Trans. Pattern Anal. Mach. Intell. 33(2), 224–239 (2011)
Parra-Denis, E.: Analyse morphologique 3D de particules de forme complexes: application aux intermétalliques dans les alliages d’aluminium. Ph.D. thesis, Université Jean Monnet de Saint-Etienne (2007)
Perlin, K.: Course in advanced image synthesis. In: ACM SIGGRAPH Conference, vol. 18 (1984)
Salembier, P., Oliveras, A., Garrido, L.: Motion connected operators for image sequences. In: European Signal Processing Conference, vol. 96, pp. 1083–1086 (1996)
Salembier, P., Oliveras, A., Garrido, L.: Antiextensive connected operators for image and sequence processing. IEEE Trans. Image Process. 7, 555–570 (1998)
Schmitt, M.: Des algorithmes morphologiques à l’intelligence artificielle. Ph.D. thesis, Ecole Nationale Supérieure des Mines de Paris (1989)
Serra, J.: Image Analysis and Mathematical Morphology, vol. 1. Academic Press, London (1982)
Serra, J.: Image Analysis and Mathematical Morphology, vol. 2. Academic Press, London (1988)
Serra, J., Vincent, L.: An overview of morphological filtering. Circuits Syst. Signal Process. 11(1), 47–108 (1992)
Soille, P.: Generalized geodesic distances applied to interpolation and shape description. In: Mathematical Morphology and Its Applications to Image Processing, pp. 193–200 (1994)
Soille, P.: Generalized geodesy via geodesic time. Pattern Recognit. Lett. 15(12), 1235–1240 (1994)
Soille, P.: Morphological carving. Pattern Recognit. Lett. 25(5), 543–550 (2004). doi:10.1016/j.patrec.2003.12.007
Talbot, H., Appleton, B.: Efficient complete and incomplete path openings and closings. Image Vis. Comput. 25(4), 416–425 (2007)
Urbach, E.R., Wilkinson, M.H.F.: Shape-only granulometries and grey-scale shape filters. In: Proc. Int. Symp. Math. Morphology, vol. 2002, pp. 305–314 (2002)
Urbach, E.R., Roerdink, J.B.T.M., Wilkinson, M.H.F.: Connected shape-size pattern spectra for rotation and scale-invariant classification of gray-scale images. IEEE Trans. Pattern Anal. Mach. Intell. 29(2), 272–285 (2007)
Vincent, L.: Grayscale area openings and closings, their efficient implementation and applications. In: First Workshop on Mathematical Morphology and its Applications to Signal Processing, pp. 22–27 (1993)
Wilkinson, M.H.F., Westenberg, M.A.: Shape preserving filament enhancement filtering. In: Medical Image Computing and Computer-Assisted Intervention. Lecture Notes in Computer Science, vol. 2208, pp. 770–777 (2001)
Acknowledgements
This work was made possible thanks to the support of the “Pôle ASTech” and the “Pôle Nucléaire de Bourgogne”, and has been financed by the French “Département de Seine et Marne”.
The authors are grateful to Ms Raviart, working in the “Centre des matériaux, MINES ParisTech”, for help with the optical microscope.
The authors also wish to thank the anonymous reviewers for their comments and suggestions.
Author information
Authors and Affiliations
Corresponding author
Appendix: Random Connected Component Models
Appendix: Random Connected Component Models
Many strategies exist to generate random binary shapes. The results of the comparison between the geodesic diameter and its approximations are linked to the method used to generate these shapes. Therefore, we use five different methods to have a high variety of objects. For each method, 5 realisations are presented in Fig. 15. The size of the support of these random shapes is a 500 by 500 pixels square.
Convex Shape
A random number of points (between 10 and 100) are randomly and uniformly picked on D. The final connected component is the convex hull of these points.
Pixel Aggregation
This method is used to generate relatively dense objects, which are almost convex. The set is initialised with a single point. At each iteration, a randomly chosen neighbour point is added to the set. The procedure is iterated a random number of times.
Ball Aggregation
This method uses the same process as the pixel aggregation method, except that we iteratively aggregate a ball instead of a point to the set. The ball radius is chosen randomly between 5 and 40 pixels, for each ball. The generated shapes are much more complex than the shapes generated using the pixel aggregation method.
Random Walk
We start from a ball in the centre of the domain. Then, we use a Brownian motion to choose the next location of the ball. At each iteration, the radius of the ball is chosen randomly.
Perlin Noise
Perlin noise [22] is a procedural texture primitive. It has a pseudo random appearance that is highly controllable and multiscale. Figure 14 provides an illustration of a realisation of this noise. By thresholding this image, we get a set of objects and we select the biggest CC of this set (Fig. 14). Some resulting objects can be very smooth, whereas others can have a high tortuosity.
Rights and permissions
About this article
Cite this article
Morard, V., Decencière, E. & Dokládal, P. Efficient Geodesic Attribute Thinnings Based on the Barycentric Diameter. J Math Imaging Vis 46, 128–142 (2013). https://doi.org/10.1007/s10851-012-0374-7
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10851-012-0374-7