Top Banner
Integral Geometry Descriptors for Characterizing Emphysema and Lung Fibrosis in HRCT Images Michal Charemza 1 , Elke Th¨ onnes 1,2 , Abhir Bhalerao 3 , David Parr 4 1 Centre for Scientific Computing, University of Warwick, UK, [email protected]. 2 Department of Statistics, University of Warwick, UK, [email protected]. 3 Department of Computer Science, University of Warwick, UK, [email protected]. 4 Department of Radiology, University Hospital Coventry and Warwickshire NHS Trust, UK, [email protected]. Integral geometry descriptors (Minkowski functionals) are used to characterize local textural properties of lung parenchyma from HRCT images. They describe the morphology & topol- ogy of 2D and 3D binary structures, and have been shown to be effective in describing properties of complex and disordered media. We propose their use for detecting and grading emphysema and fibrosis in HRCT images. We present some illustrative results and make proposals for their use in a larger validation study. Abstract Defined on sets in 3D space. They are proportional to [4]: Volume V Surface Area S Mean Breadth B : a measure of average width Euler-Poincar´ e Characteristic χ: Connectivity number = # connected components -# tunnels +# cavities Minkowski Functionals: Definition Fundamental role over certain functions of binary images All Real-valued, additive, isometry invariant, continuous func- tions on a compact 3D binary image A are of form 3 i=0 α i W i (A) where the α i are constants and the W i are the Minkowski functionals. Are themselves additive: for K 1 , K 2 compact parts of a 3D image W i (K 1 K 2 )= W i (K 1 )+ W i (K 2 ) - W i (K 1 K 2 ). (1) Are themselves isometry invariant: do not change their value under rotation or translation. Can be extended to the interior of compact set K by W i (K )=(-1) 3+i+dim K (W i ) (2) Minkowski Functionals: Properties Finding functionals reduces to counting cells, faces, edges and vertices of image seen as collection of cubic voxels. The Minkowski functionals for a cube, rectangle, line and vertex all known analytically. = can use (1) and (2) repeatedly to find: V = n 3 n 3 : # voxels S = -6n 3 +2n 2 n 2 : # faces B =3n 3 /2 - n 2 + n 1 /2 n 1 : # edges χ = -n 3 + n 2 - n 1 + n 0 n 0 : # vertices To find n i We use the Algorithm of Equa- tions [1], where the faces, edges and vertices of the voxel (light red) counted assuming the 13 adjacent voxels (black) already counted. To describe heterogeneity functionals are found on win- dows of fixed size over the 3D image Finding Functionals Models, left-to-right: healthy, fibrotic and honeycomb. Based on Boolean model and Voronoi tessellation. [4] Supervised kNN classification of MFs Projection of MFs onto first PCA modes Results: Synthetic Model Functionals found on thresholded in-vivo lung HRCT scan. Projection of 4 MFs onto 3 PCA principal modes Results: Lung Data Promising results: Minkowski functionals do appear to separate tissue types. Larger validation study required. More advanced probabilistic model development needed. Work concurrent to, but independent from, Boehm et al. [2, 3]. Conclusion [1] I. Blasquez and J.-F. Poiraudeau. Efficient processing of Minkowski functionals on a 3D binary image using binary decision diagrams. Journal of MSCG, 11(1), 2003. [2] H. Boehm, C. Fink, U. Attenberger, C. Becker, J. Behr, and M. Reiser. Automated Classification of Normal and Pathologic Pulmonary Tissue by Topological Texture Features Extracted rom Multi-detector CT in 3D. European Radiology, 18, 2008. DOI: 10.1007/s00330-008-1082-y. [3] H. Boem, C. Fink, C. Becker, and M. Reiser. Automated Characterization of Normal and Pathologic Lung Tissue by Topological Texture Analysis of Multi-Detector CT. In Giger and Karssemeijer, editors, Medical Imaging 2007: Computer- Aided Diagnosis, Proceedings of SPIE Vol 6514, page DOI: 10.117/12.702697, 2007. [4] D. Stoyan, W. Kendall, and J. Mecke. Stochastic Geometry and its Application. Probability and Statistics. John Wiley & Sons, New York, 1995.
1

Charemza 08 Poster MICCAI Integral Geometry Descriptors for Characterizing Emphysema and Lung Fibrosis in HRCT Images

Feb 05, 2016

Download

Documents

Poster
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Charemza 08 Poster MICCAI Integral Geometry Descriptors for Characterizing Emphysema and Lung Fibrosis in HRCT Images

Integral Geometry Descriptors for Characterizing Emphysema and Lung Fibrosis in HRCT ImagesMichal Charemza1, Elke Thonnes1,2, Abhir Bhalerao3, David Parr41Centre for Scientific Computing, University of Warwick, UK, [email protected]. 2Department of Statistics, University of Warwick, UK, [email protected]. 3 Department of Computer Science, University of Warwick, UK, [email protected]. 4Department of Radiology, University Hospital Coventry and Warwickshire NHS Trust, UK, [email protected].

Integral geometry descriptors (Minkowski functionals) areused to characterize local textural properties of lung parenchymafrom HRCT images. They describe the morphology & topol-ogy of 2D and 3D binary structures, and have been shown to beeffective in describing properties of complex and disordered media.

We propose their use for detecting and grading emphysema andfibrosis in HRCT images. We present some illustrative results andmake proposals for their use in a larger validation study.

Abstract

Defined on sets in 3D space. They are proportional to [4]:

• Volume V

• Surface Area S

• Mean Breadth B: a measure of average width

• Euler-Poincare Characteristic χ: Connectivity number =# connected components −# tunnels +# cavities

Minkowski Functionals: Definition

Fundamental role over certain functions of binary images

• All Real-valued, additive, isometry invariant, continuous func-tions on a compact 3D binary image A are of form

3∑i=0

αiWi(A)

where the αi are constants and the Wi are the Minkowskifunctionals.

• Are themselves additive: for K1, K2 compact parts of a 3Dimage

Wi(K1 ∪K2) = Wi(K1) + Wi(K2)−Wi(K1 ∩K2). (1)

• Are themselves isometry invariant: do not change theirvalue under rotation or translation.

• Can be extended to the interior of compact set K by

Wi(K◦) = (−1)3+i+dim K(Wi) (2)

Minkowski Functionals: Properties

Finding functionals reduces to counting cells, faces, edges andvertices of image seen as collection of cubic voxels.

• The Minkowski functionals for a cube, rectangle, line and vertexall known analytically.

• =⇒ can use (1) and (2) repeatedly to find:

V = n3 n3: # voxels

S = −6n3 + 2n2 n2: # faces

B = 3n3/2− n2 + n1/2 n1: # edges

χ = −n3 + n2 − n1 + n0 n0: # vertices

• To find ni We use the Algorithm of Equa-tions [1], where the faces, edges and verticesof the voxel (light red) counted assuming the13 adjacent voxels (black) already counted.

• To describe heterogeneity functionals are found on win-dows of fixed size over the 3D image

Finding Functionals

Models, left-to-right: healthy, fibrotic and honeycomb.

• Based on Boolean model and Voronoi tessellation. [4]

Supervised kNN classification of MFs

Projection of MFs onto first PCA modes

Results: Synthetic Model

Functionals found on thresholded in-vivo lung HRCT scan.

Projection of 4 MFs onto3 PCA principal modes

Results: Lung Data

Promising results: Minkowski functionals do appear to separatetissue types.

• Larger validation study required.

• More advanced probabilistic model development needed.

• Work concurrent to, but independent from, Boehm et al. [2, 3].

Conclusion

[1] I. Blasquez and J.-F. Poiraudeau. Efficient processing of Minkowski functionals on a 3D binary image using binary decisiondiagrams. Journal of MSCG, 11(1), 2003.

[2] H. Boehm, C. Fink, U. Attenberger, C. Becker, J. Behr, and M. Reiser. Automated Classification of Normal and PathologicPulmonary Tissue by Topological Texture Features Extracted rom Multi-detector CT in 3D. European Radiology, 18, 2008.DOI: 10.1007/s00330-008-1082-y.

[3] H. Boem, C. Fink, C. Becker, and M. Reiser. Automated Characterization of Normal and Pathologic Lung Tissue byTopological Texture Analysis of Multi-Detector CT. In Giger and Karssemeijer, editors, Medical Imaging 2007: Computer-Aided Diagnosis, Proceedings of SPIE Vol 6514, page DOI: 10.117/12.702697, 2007.

[4] D. Stoyan, W. Kendall, and J. Mecke. Stochastic Geometry and its Application. Probability and Statistics. John Wiley &Sons, New York, 1995.