
- •Biological and Medical Physics, Biomedical Engineering
- •Medical Image Processing
- •Preface
- •Contents
- •Contributors
- •1.1 Medical Image Processing
- •1.2 Techniques
- •1.3 Applications
- •1.4 The Contribution of This Book
- •References
- •2.1 Introduction
- •2.2 MATLAB and DIPimage
- •2.2.1 The Basics
- •2.2.2 Interactive Examination of an Image
- •2.2.3 Filtering and Measuring
- •2.2.4 Scripting
- •2.3 Cervical Cancer and the Pap Smear
- •2.4 An Interactive, Partial History of Automated Cervical Cytology
- •2.5 The Future of Automated Cytology
- •2.6 Conclusions
- •References
- •3.1 The Need for Seed-Driven Segmentation
- •3.1.1 Image Analysis and Computer Vision
- •3.1.2 Objects Are Semantically Consistent
- •3.1.3 A Separation of Powers
- •3.1.4 Desirable Properties of Seeded Segmentation Methods
- •3.2 A Review of Segmentation Techniques
- •3.2.1 Pixel Selection
- •3.2.2 Contour Tracking
- •3.2.3 Statistical Methods
- •3.2.4 Continuous Optimization Methods
- •3.2.4.1 Active Contours
- •3.2.4.2 Level Sets
- •3.2.4.3 Geodesic Active Contours
- •3.2.5 Graph-Based Methods
- •3.2.5.1 Graph Cuts
- •3.2.5.2 Random Walkers
- •3.2.5.3 Watershed
- •3.2.6 Generic Models for Segmentation
- •3.2.6.1 Continuous Models
- •3.2.6.2 Hierarchical Models
- •3.2.6.3 Combinations
- •3.3 A Unifying Framework for Discrete Seeded Segmentation
- •3.3.1 Discrete Optimization
- •3.3.2 A Unifying Framework
- •3.3.3 Power Watershed
- •3.4 Globally Optimum Continuous Segmentation Methods
- •3.4.1 Dealing with Noise and Artifacts
- •3.4.2 Globally Optimal Geodesic Active Contour
- •3.4.3 Maximal Continuous Flows and Total Variation
- •3.5 Comparison and Discussion
- •3.6 Conclusion and Future Work
- •References
- •4.1 Introduction
- •4.2 Deformable Models
- •4.2.1 Point-Based Snake
- •4.2.1.1 User Constraint Energy
- •4.2.1.2 Snake Optimization Method
- •4.2.2 Parametric Deformable Models
- •4.2.3 Geometric Deformable Models (Active Contours)
- •4.2.3.1 Curve Evolution
- •4.2.3.2 Level Set Concept
- •4.2.3.3 Geodesic Active Contour
- •4.2.3.4 Chan–Vese Deformable Model
- •4.3 Comparison of Deformable Models
- •4.4 Applications
- •4.4.1 Bone Surface Extraction from Ultrasound
- •4.4.2 Spinal Cord Segmentation
- •4.4.2.1 Spinal Cord Measurements
- •4.4.2.2 Segmentation Using Geodesic Active Contour
- •4.5 Conclusion
- •References
- •5.1 Introduction
- •5.2 Imaging Body Fat
- •5.3 Image Artifacts and Their Impact on Segmentation
- •5.3.1 Partial Volume Effect
- •5.3.2 Intensity Inhomogeneities
- •5.4 Overview of Segmentation Techniques Used to Isolate Fat
- •5.4.1 Thresholding
- •5.4.2 Selecting the Optimum Threshold
- •5.4.3 Gaussian Mixture Model
- •5.4.4 Region Growing
- •5.4.5 Adaptive Thresholding
- •5.4.6 Segmentation Using Overlapping Mosaics
- •5.6 Conclusions
- •References
- •6.1 Introduction
- •6.2 Clinical Context
- •6.3 Vessel Segmentation
- •6.3.1 Survey of Vessel Segmentation Methods
- •6.3.1.1 General Overview
- •6.3.1.2 Region-Growing Methods
- •6.3.1.3 Differential Analysis
- •6.3.1.4 Model-Based Filtering
- •6.3.1.5 Deformable Models
- •6.3.1.6 Statistical Approaches
- •6.3.1.7 Path Finding
- •6.3.1.8 Tracking Methods
- •6.3.1.9 Mathematical Morphology Methods
- •6.3.1.10 Hybrid Methods
- •6.4 Vessel Modeling
- •6.4.1 Motivation
- •6.4.1.1 Context
- •6.4.1.2 Usefulness
- •6.4.2 Deterministic Atlases
- •6.4.2.1 Pioneering Works
- •6.4.2.2 Graph-Based and Geometric Atlases
- •6.4.3 Statistical Atlases
- •6.4.3.1 Anatomical Variability Handling
- •6.4.3.2 Recent Works
- •References
- •7.1 Introduction
- •7.2 Linear Structure Detection Methods
- •7.3.1 CCM for Imaging Diabetic Peripheral Neuropathy
- •7.3.2 CCM Image Characteristics and Noise Artifacts
- •7.4.1 Foreground and Background Adaptive Models
- •7.4.2 Local Orientation and Parameter Estimation
- •7.4.3 Separation of Nerve Fiber and Background Responses
- •7.4.4 Postprocessing the Enhanced-Contrast Image
- •7.5 Quantitative Analysis and Evaluation of Linear Structure Detection Methods
- •7.5.1 Methodology of Evaluation
- •7.5.2 Database and Experiment Setup
- •7.5.3 Nerve Fiber Detection Comparison Results
- •7.5.4 Evaluation of Clinical Utility
- •7.6 Conclusion
- •References
- •8.1 Introduction
- •8.2 Methods
- •8.2.1 Linear Feature Detection by MDNMS
- •8.2.2 Check Intensities Within 1D Window
- •8.2.3 Finding Features Next to Each Other
- •8.2.4 Gap Linking for Linear Features
- •8.2.5 Quantifying Branching Structures
- •8.3 Linear Feature Detection on GPUs
- •8.3.1 Overview of GPUs and Execution Models
- •8.3.2 Linear Feature Detection Performance Analysis
- •8.3.3 Parallel MDNMS on GPUs
- •8.3.5 Results for GPU Linear Feature Detection
- •8.4.1 Architecture and Implementation
- •8.4.2 HCA-Vision Features
- •8.4.3 Linear Feature Detection and Analysis Results
- •8.5 Selected Applications
- •8.5.1 Neurite Tracing for Drug Discovery and Functional Genomics
- •8.5.2 Using Linear Features to Quantify Astrocyte Morphology
- •8.5.3 Separating Adjacent Bacteria Under Phase Contrast Microscopy
- •8.6 Perspectives and Conclusions
- •References
- •9.1 Introduction
- •9.2 Bone Imaging Modalities
- •9.2.1 X-Ray Projection Imaging
- •9.2.2 Computed Tomography
- •9.2.3 Magnetic Resonance Imaging
- •9.2.4 Ultrasound Imaging
- •9.3 Quantifying the Microarchitecture of Trabecular Bone
- •9.3.1 Bone Morphometric Quantities
- •9.3.2 Texture Analysis
- •9.3.3 Frequency-Domain Methods
- •9.3.4 Use of Fractal Dimension Estimators for Texture Analysis
- •9.3.4.1 Frequency-Domain Estimation of the Fractal Dimension
- •9.3.4.2 Lacunarity
- •9.3.4.3 Lacunarity Parameters
- •9.3.5 Computer Modeling of Biomechanical Properties
- •9.4 Trends in Imaging of Bone
- •References
- •10.1 Introduction
- •10.1.1 Adolescent Idiopathic Scoliosis
- •10.2 Imaging Modalities Used for Spinal Deformity Assessment
- •10.2.1 Current Clinical Practice: The Cobb Angle
- •10.2.2 An Alternative: The Ferguson Angle
- •10.3 Image Processing Methods
- •10.3.1 Previous Studies
- •10.3.2 Discrete and Continuum Functions for Spinal Curvature
- •10.3.3 Tortuosity
- •10.4 Assessment of Image Processing Methods
- •10.4.1 Patient Dataset and Image Processing
- •10.4.2 Results and Discussion
- •10.5 Summary
- •References
- •11.1 Introduction
- •11.2 Retinal Imaging
- •11.2.1 Features of a Retinal Image
- •11.2.2 The Reason for Automated Retinal Analysis
- •11.2.3 Acquisition of Retinal Images
- •11.3 Preprocessing of Retinal Images
- •11.4 Lesion Based Detection
- •11.4.1 Matched Filtering for Blood Vessel Segmentation
- •11.4.2 Morphological Operators in Retinal Imaging
- •11.5 Global Analysis of Retinal Vessel Patterns
- •11.6 Conclusion
- •References
- •12.1 Introduction
- •12.1.1 The Progression of Diabetic Retinopathy
- •12.2 Automated Detection of Diabetic Retinopathy
- •12.2.1 Automated Detection of Microaneurysms
- •12.3 Image Databases
- •12.4 Tortuosity
- •12.4.1 Tortuosity Metrics
- •12.5 Tracing Retinal Vessels
- •12.5.1 NeuronJ
- •12.5.2 Other Software Packages
- •12.6 Experimental Results and Discussion
- •12.7 Summary and Future Work
- •References
- •13.1 Introduction
- •13.2 Volumetric Image Visualization Methods
- •13.2.1 Multiplanar Reformation (2D slicing)
- •13.2.2 Surface-Based Rendering
- •13.2.3 Volumetric Rendering
- •13.3 Volume Rendering Principles
- •13.3.1 Optical Models
- •13.3.2 Color and Opacity Mapping
- •13.3.2.2 Transfer Function
- •13.3.3 Composition
- •13.3.4 Volume Illumination and Illustration
- •13.4 Software-Based Raycasting
- •13.4.1 Applications and Improvements
- •13.5 Splatting Algorithms
- •13.5.1 Performance Analysis
- •13.5.2 Applications and Improvements
- •13.6 Shell Rendering
- •13.6.1 Application and Improvements
- •13.7 Texture Mapping
- •13.7.1 Performance Analysis
- •13.7.2 Applications
- •13.7.3 Improvements
- •13.7.3.1 Shading Inclusion
- •13.7.3.2 Empty Space Skipping
- •13.8 Discussion and Outlook
- •References
- •14.1 Introduction
- •14.1.1 Magnetic Resonance Imaging
- •14.1.2 Compressed Sensing
- •14.1.3 The Role of Prior Knowledge
- •14.2 Sparsity in MRI Images
- •14.2.1 Characteristics of MR Images (Prior Knowledge)
- •14.2.2 Choice of Transform
- •14.2.3 Use of Data Ordering
- •14.3 Theory of Compressed Sensing
- •14.3.1 Data Acquisition
- •14.3.2 Signal Recovery
- •14.4 Progress in Sparse Sampling for MRI
- •14.4.1 Review of Results from the Literature
- •14.4.2 Results from Our Work
- •14.4.2.1 PECS
- •14.4.2.2 SENSECS
- •14.4.2.3 PECS Applied to CE-MRA
- •14.5 Prospects for Future Developments
- •References
- •15.1 Introduction
- •15.2 Acquisition of DT Images
- •15.2.1 Fundamentals of DTI
- •15.2.2 The Pulsed Field Gradient Spin Echo (PFGSE) Method
- •15.2.3 Diffusion Imaging Sequences
- •15.2.4 Example: Anisotropic Diffusion of Water in the Eye Lens
- •15.2.5 Data Acquisition
- •15.3 Digital Processing of DT Images
- •15.3.2 Diagonalization of the DT
- •15.3.3 Gradient Calibration Factors
- •15.3.4 Sorting Bias
- •15.3.5 Fractional Anisotropy
- •15.3.6 Other Anisotropy Metrics
- •15.4 Applications of DTI to Articular Cartilage
- •15.4.1 Bovine AC
- •15.4.2 Human AC
- •References
- •Index

7 Detecting and Analyzing Linear Structures in Biomedical Images: A Case Study... |
163 |
Fig. 7.13 Boxplots showing the NFL scores for each of the NDS groups calculated manually (left) and automatically (right)
(3.68–33.91) for the manual analysis to (1.22–20.03) for the automated analysis. ANOVA analysis results in a p-value for discrimination among these groups which is slightly higher for the automated than the manual analysis, though both are highly significant (p ≈ 0) (Table 7.2).
7.6 Conclusion
The analysis of CCM images requires the identification of fiber-like structures with low contrast in noisy images. This is a requirement shared by a number of imaging applications in biology, medicine, and other fields. A number of methods have been applied in these applications, and we have compared some of these, and more generic methods, with a dual-model detection algorithm devised for this study. The comparison used a large set of images with manual ground-truth. In terms of both error-rates (pixel misclassification) and signal-to-noise ratio, the dual model achieved the highest performance. It seems reasonable to propose that this filter is likely to prove equally useful in applications of a similar nature.
The question of the clinical utility of the method was also addressed. The evaluation has shown that the automatic analysis is consistent with the manual ground-truth with a correlation of r = 0.93. Similarity in grouping control and patient subjects between manual and automated analysis was also achieved with (p ≈ 0). Therefore, we conclude that automated analysis of corneal nerve fibers is a much quicker and potentially more reliable practical alternative to manual analysis due to its consistency and immunity to the inter/intra-observer variability. These prosperities will help deliver CCM from a research tool to a practical and potentially large scale clinical technique for the assessment of neuropathic severity.
164 |
M.A. Dabbah et al. |
Acknowledgments The authors thank Dr. Mitra Tavakoli and Mr. Ioannis Petropoulos for capturing the CCM images and providing the manually delineated ground-truth. The authors also thank the Juvenile Diabetes Research Foundation (JDRF) for the financial support under the scholar grant 17–2008–1031.
References
1.Malik, R.A., Kallinikos, P., Abbott, C.A., van Schie, C.H.M., Morgan, P., Efron, N., Boulton, A.J.M.: Corneal confocal microscopy: A non-invasive surrogate of nerve fibre damage and repair in diabetic patients. Diabetologia 46, 683–688 (2003)
2.Kallinikos, P., Berbanu, M., O’Donnell, C., Boulton, A., Efron, N., Malik, R.: Corneal nerve tortuosity in diabetic patients with neuropathy. Invest. Ophthalmol. Vis. Sci. 45, 418–422 (2004)
3.Hossain, P., Sachdev, A., Malik, R.A.: Early detection of diabetic peripheral neuropathy with corneal confocal microscopy. Lancet 366, 1340–1343 (2005)
4.Ruggeri, A., Scarpa, F., Grisan, E.: Analysis of corneal images for the recognition of nerve structures. In: IEEE Conference of the Engineering in Medicine and Biology Society (EMBS),
pp.4739–4742, September 2006
5.Dabbah, M.A., Graham, J., Tavakoli, M., Petropoulos, Y., Malik, R.A.: Nerve fibre extraction in confocal corneal microscopy images for human diabetic neuropathy detection using gabor filters. In: Medical Image Understanding and Analysis (MIUA), pp. 254–258, July 2009
6.Dixon, R.N., Taylor, C.J.: Automated asbestos fibre counting. Mach. Aided Image Anal. 178–185 (1979)
7.Bryson, N., Dixon, R.N., Hunter, J.J., Taylor, C.J.: Contextual classification of cracks. Image Vis. Comput. 12, 149–154 (1994)
8.Zwiggelaar, R., Astley, S., Boggis, C., Taylor, C.: Linear structures in mammographic images: Detection and classification. IEEE Trans. Med. Imag. 23, 1077–1086 (2004)
9.Jain, A.K., Farrokhnia, F.: Unsupervised texture segmentation using Gabor filters. Pattern Recogn. 24, 1167–1186 (1991)
10.Lorenz, C., Carlsen, I., Buzug, T.: Fassnacht, C., Weese, J.: Multi-scale line segmentation with automatic estimation of width, contrast and tangential direction in 2D and 3D medical images. In: CVRMed-MRCAS’97, pp. 233–242 (1997)
11.Sato, Y., Nakajima, S., Atsumi, H., Koller, T., Gerig, G., Yoshida, S., Kikinis, R.: 3D multiscale line filter for segmentation and visualization of curvilinear structures in medical images. In: CVRMed-MRCAS’97, pp. 213–222 (1997)
12.Frangi, A.F., Niessen, W.J., Vincken, K.L., Viergever, M.A.: Multiscale vessel enhancement filtering. In: Medical Image Computing and Computer-Assisted Interventation (MICCAI),
pp.130–137, July 1998
13.Krissian, K., Malandain, G., Ayache, N., Vaillant, R., Trousset, Y.: Model-based detection of tubular structures in 3D images. Comput. Vis. Image Underst. 80, 130–171 (2000)
14.Kingsbury, N.: Complex wavelets for shift invariant analysis and filtering of signals. Appl. Comput. Harmonic Anal. 10, 234–253 (2001)
15.Kingsbury, N.: Rotation-invariant local feature matching with complex wavelets. In: European Conference on Signal Processing (EUSIPCO), Florence, pp. 4–8, 2006
16.Sadeghzadeh, R., Berks, M., Astley, S., Taylor, C.: Detection of retinal blood vessels using complex wavelet transforms and random forest classification. In: Proceedings of Medical Image Understanding and Analysis (MIUA), pp. 127–131 (2010)
17.Chen, Z., Berks, M.: Astley, S., Taylor, C., Classification of linear structures in mammograms using random forests. In: Digital Mammography, pp. 153–160 (2010)
18.Felsberg, M., Sommer, G.: The Monogenic Signal. IEEE Trans. Signal Process. 49, 3136–3144 (2001)
7 Detecting and Analyzing Linear Structures in Biomedical Images: A Case Study... |
165 |
19.Pan, X.B., Brady, M., Highnam, R., Declerck, J.: The use of multi-scale monogenic signal on structure orientation identification and segmentation. In: Digital Mammography, pp. 601–608 (2006)
20.Ali, R., Gooding, M., Christlieb, M., Brady, M.: Advanced phase-based segmentation of multiple cells from brightfield microscopy images. In: IEEE International Symposium on Biomedical Imaging: From Nano to Macro (ISBI), pp. 181–184 (2008)
21.Minsky, M.: Microscopy apparatus US Patent 3013467. U. S. P. Office, Ed. US, 1961
22.Dyck, P., Norell, J., Tritschler, H., Schuette, K., Samigullin, R., Ziegler, D., Bastyr, E., Litchy, W., O’Brien, P.: Challenges in design of multicenter trials: end points assessed longitudinally for change and monotonicity. Diabetes Care 30, 2619–2625 (2007)
23.Umapathi, T., Tan, W., Loke, S., Soon, P., Tavintharan, S., Chan, Y.: Intraepidermal nerve fiber density as a marker of early diabetic neuropathy. Muscle Nerve 35, 591–598 (2007)
24.Loseth, S., Stalberg, E., Jorde, R., Mellgren, S.: Early diabetic neuropathy: thermal thresholds
and intraepidermal nerve fibre density in patients with normal nerve conduction studies.
J. Neurol. 255, 1197–1202 (2008)
25.Malik, R., Tesfaye, S., Newrick, P., Walker, D., Rajbhandari, S., Siddique, I., Sharma, A., Boulton, A., King, R., Thomas, P., Ward, J.: Sural nerve pathology in diabetic patients with minimal but progressive neuropathy. Diabetologia 48, 578–585 (2005)
26.Malik, R., Veves, A., Walker, D., Siddique, I., Lye, R., Schady, W., Boulton, A.: Sural nerve fibre pathology in diabetic patients with mild neuropathy: relationship to pain, quantitative sensory testing and peripheral nerve electrophysiology. Acta Neuropathol. (Berl.) 101, 367–374 (2001)
27.Novella, S., Inzucchi, S., Goldstein, J.: The frequency of undiagnosed diabetes and impaired glucose tolerance in patients with idiopathic sensory neuropathy. Muscle Nerve 24, 1229–1231 (2001)
28.Singleton, J., Smith, A., Bromberg, M.: Increased prevalence of impaired glucose tolerance in patients with painful sensory neuropathy. Diabetes Care 24, 1448–1453 (2001)
29.Sumner, C., Sheth, S., Griffin, J., Cornblath, D., Polydefkis, M.: The spectrum of neuropathy in diabetes and impaired glucose tolerance. Neurology 60, 108–111 (2003)
30.Quattrini, C., Tavakoli, M., Jeziorska, M., Kallinikos, P., Tesfaye, S., Finnigan, J., Marshall, A., Boulton, A.J.M., Efron, N., Malik, R.A.: Surrogate markers of small fiber damage in human diabetic neuropathy. Diabetes 56, 2148–2154 (2007)
31.Tavakoli, M., Kallinikos, P.A., Efron, N., Boulton, A.J.M., Malik, R.A.: Corneal sensitivity is reduced and relates to the severity of neuropathy in patients with diabetes. Diabetes Care 30, 1895–1897 (2007)
32.Mehra, S., Tavakoli, M., Kallinikos, P.A., Efron, N., Boulton, A.J.M., Augustine, T., Malik, R.A.: Corneal confocal microscopy detects early nerve regeneration after pancreas transplantation in patients with type 1 diabetes. Diabetes Care 30, 2608–2612 (2007)
33.Tavakoli, M., Marshall, A., Thompson, L., Kenny, M., Waldek, S., Efron, N., Malik, R.A.: Corneal confocal microscopy: A novel noninvasive means to diagnose neuropathy in patients with Fabry disease. Muscle and Nerve 40, 976–984 (2009)
34.Tavakoli, M., Marshall, A., Pitceathly, R., Fadavi, H., Gow, D., Roberts, M.E., Efron, N., Boulton, A.J.M., Malik, R.A.: Corneal confocal microscopy: A novel means to detect nerve fibre damage in idiopathic small fibre neuropathy. Exp. Neurol. 223, 245–250 (2010)
35.Tavakoli, M., Quattrini, C., Abbott, C., Kallinikos, P., Marshall, A., Finnigan, J., Morgan, P., Efron, N., Boulton, A.J.M., Malik, R.A.: Corneal confocal microscopy: A novel non-invasive test to diagnose and stratify the severity of human diabetic neuropathy. Diabetes Care 33, 1792–1797 (2010)
36.Dougherty, G., Johnson, M.J., Wiers, M.D.: Measurement of retinal vascular tortuosity and its application to retinal pathologies. Med. Biol. Eng. Comput. 48, 87–95 (2010)
37.Dabbah, M.A., Graham, J., Tavakoli, M., Petropoulos, Y., Malik, R.A.: Dual-model automatic detection of nerve-fibres in corneal confocal microscopy images. In: The International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI), vol. 6361, pp. 300–307 (2010)
166 |
M.A. Dabbah et al. |
38.Daugman, J.G.: Two-dimensional spectral analysis of cortical receptive field profiles. Vis. Res. 20, 847–856 (1980)
39.Rao, A.R.: A taxonomy for texture description and identification. Springer, New York (1990)
40.Kass, M., Witkin, A.: Analyzing oriented patterns. Comp. Vis. Graph. Image Process. 37, 362–385 (1987)
41.Hong, L., Wan, Y., Jain, A.: Fingerprint image enhancement: algorithm and performance evaluation. IEEE Trans. Pattern Anal. Mach. Intell. 20, 777–789 (1998)
42.Zhang, T.Y., Suen, C.Y.: A fast parallel algorithm for thinning digital patterns. Commun. ACM 27, 236–239 (1984)
43.Rutovitz, D.: Pattern recognition. J. R. Stat. Soc. A (General) 129, 504–530 (1966)
44.Kingsbury, N.: Dual-tree complex wavelet transform pack. http://www-sigproc.eng.cam.ac. uk/\ ngk/. Accessed June 2002
45.Kovesi, P.: An implementation of Felsberg’s monogenic filters. http://www.csse.uwa.edu. au1pk/research/matlabfns/. Accessed August 2005
46.Kroon, D.J., Schrijver, M.: Hessian based Frangi Vesselness filter. http://www.mathworks.co. uk/. Accessed October 2009
47.Abbott, C.A., Carrington, A.L., Ashe, H., Bath, S., Every, L.C., Griffiths, J., Hann, A.W., Hussein, A., Jackson, N., Johnson, K.E., Ryder, C.H., Torkington, R., Ross, E.R.E.V., Whalley, A.M., Widdows, P., Williamson, S., Boulton, A.J.M.: The North-West diabetes foot care study: Incidence of, and risk factors for, new diabetic foot ulceration in a community-based patient cohort. Diabetic Med. 19, 377–384 (2002)