- •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
10 Applications of Medical Image Processing in the Diagnosis and Treatment... |
229 |
Fig. 10.2 Post-operative X-rays single rod showing anterior (left) and dual rod posterior (right) implants for scoliosis correction
conservative treatment in which the patient is asked to wear an orthotic brace which attempts to exert corrective forces on the deformed spine. Surgical approaches to scoliosis correction (Fig. 10.2) involve attachment of implants to the spine to restore a more normal curvature. Successful surgical approaches typically achieve a 60% reduction of the deformity and prevent further progression; however, there are risks of complication and further deformity progression after surgery which mandate ongoing imaging to monitor the corrected spine.
10.2 Imaging Modalities Used for Spinal Deformity Assessment
Many spinal deformities are visible just by looking at a patient. However, due to differences in body fat levels and bone structure between patients, they cannot be accurately assessed by visual examination of the patient’s appearance. Medical imaging, therefore, plays a key role both in the monitoring of spinal deformity progression before treatment, and in assessing treatment outcomes. There are four medical imaging modalities relevant to the assessment of spinal deformities.
230 |
C. Adam and G. Dougherty |
Planar and biplanar radiography Plane radiographs (X-rays) are the gold standard for imaging spinal deformities, and are used in spine clinics worldwide. A relatively new technology developed in France, the EOS system (EOS imaging, Paris, France), uses bi-planar radiography to simultaneously obtain coronal and sagittal radiographs of standing scoliosis patients with low radiation dose.
Computed Tomography (CT) CT is currently not used routinely for clinical assessment of scoliosis due to its greater cost and higher radiation dose than plane X-rays [1, 2]. However, low dose pre-operative CT is clinically indicated in endoscopic or keyhole scoliosis surgery for planning implant placement [3], and advances in scanner technology now allow CT with much lower radiation doses than previously possible [4]. In the research context, a number of groups have used CT to assess spinal deformities (in particular vertebral rotation in the axial plane) due to the 3D information provided [5–10].
Magnetic Resonance (MR) MR imaging is sometimes used clinically to detect soft tissue abnormalities (particularly the presence of syringomyelia in scoliosis patients), but clinical use of MR for spinal deformities is generally limited. Several research studies have used MR to investigate scoliosis [11–15]. MR is a useful research tool because of the absence of ionizing radiation, although the ability of MR to define bony anatomy is limited and thus image processing of MR datasets is often labour-intensive.
Back surface topography Various optical surface topography systems have been used to directly visualise the cosmetic effects of spinal deformity by assessing back shape in 3D [16–24]. Back surface topography does not involve ionising radiation, and so is useful for functional evaluations involving repeated assessments [25]. However, the relationship between spinal deformity and back shape is complicated by factors such as body positioning, trunk rotation, body build and fat folds [26].
10.2.1 Current Clinical Practice: The Cobb Angle
Currently, the accepted measure for clinical assessment of scoliosis is the Cobb angle [27]. The Cobb angle is measured on plane radiographs by drawing a line through the superior endplate of the superior end vertebra of a scoliotic curve, and another line through the inferior endplate of the inferior-most vertebra of the same scoliotic curve, and then measuring the angle between these lines (Fig. 10.3). Clinically, many Cobb measurements are still performed manually using pencil and ruler on hardcopy X-ray films, but PACS systems (viz. computer networks) are increasingly used which allow manual Cobb measurements to be performed digitally by clinicians on the computer screen. As well as being used to assess scoliosis in the coronal plane, the Cobb angle is used on sagittal plane radiographs to assess thoracic kyphosis and lumbar lordosis.
10 Applications of Medical Image Processing in the Diagnosis and Treatment... |
231 |
Fig. 10.3 Schematic view of a scoliotic deformity showing measurement of the Cobb angle
α
Cobb angle
Although widely used for scoliosis assessment for its simplicity, the Cobb angle has several shortcomings. First, numerous studies have shown that the interand intra-observer measurement variabilities associated with the Cobb angle are high. If the vertebral endplates on a plane X-ray appear blurred due to a forward or backward tilt, considerable inter-observer errors can be introduced in the Cobb method as a result of the difficulty in selecting the endplate orientation. Such errors are even more pronounced in the presence of contour changes resulting from osteoporosis [28, 29]. Figure 10.4 shows a summary of Cobb variability studies performed between 1982 and 2005, indicating that the 95% confidence interval for the difference between two measurements by the same observer is around 5–7◦, and for two measurements by different observers is around 6–8◦. These measurement errors are large enough to make the difference between a diagnosis of progression (requiring treatment) or stability, and so introduce uncertainty into the assessment and treatment process.
Second, because of its simplicity, the Cobb angle omits potentially useful information about the shape of a scoliotic spine. Specifically, the Cobb angle cannot differentiate between a large scoliotic curve which may span 8 or 9 vertebral levels, and a small scoliotic curve which may only span 2 or 3 vertebral levels but that has the same endplate angulation due to its severity (Fig. 10.5).
10.2.2 An Alternative: The Ferguson Angle
Prior to the adoption of the Cobb angle as standard practice by the Scoliosis Research Society in 1966, a number of other scoliosis measurement techniques
232
intra-observer variability |
(degrees) |
95% CI for |
|
C. Adam and G. Dougherty
12
10
8
6
4
2
0
Fig. 10.4 Ninety five percent confidence intervals for intra-observer measurement variability using the Cobb angle from previous studies (figure based on [30])
Fig. 10.5 Simplified representation of three different scoliotic curves with increasing radius and a greater number of vertebral levels in the curve, but the
same Cobb angle |
α |
|
|
|
Cobb angle |
had been proposed, including the Ferguson angle [31]. The Ferguson angle requires identification of three landmark points, the geometric centres of the upper, apical (i.e. the most laterally deviated) and lower vertebrae in a scoliotic curve (Fig. 10.6). Not only does the Ferguson angle take into account the position of the apical vertebra, which the Cobb angle does not, it is less influenced by changes in the shape of the vertebrae [32].