- •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
12 Tortuosity as an Indicator of the Severity of Diabetic Retinopathy |
279 |
line-like structures or segments. We mitigated this by applying a three-point moving average filter to the exported coordinates prior to calculating tortuosities.
NeuronJ is a semiautomated tracing method, as it requires the user to guide the vessel tracing. However, it has major advantages in that there is no need to preprocess the images (e.g., by histogram equalization) nor is segmentation required (i.e., the tracing can be performed on a grayscale image (usually the green plane of an RGB color image) directly).
12.5.2 Other Software Packages
Hessian-based detectors are computationally expensive. HCA-vision is a software platform for the automated detection and analysis of linear features based on multidirectional nonmaximum suppression (MDNMS). HCA-vision has been successfully applied to a number of applications including neurite tracing for drug discovery and functional genomics [60], quantifying astrocyte morphology [61], and separating adjacent bacteria under phase contrast microscopy [62]. It is discussed in detail in Chapter 8. The software can be downloaded after completing a request at http://www.hca-vision.com/product download hca vision.html. The large noise content in our images precluded us from utilizing the program successfully.
Retinal vessel centerline extraction can also be achieved using multiscale matched filters, with a vessel confidence measure defined as a projection of a vector formed from a normalized pixel neighborhood on to a normalized ideal vessel profile [54]. Vessel boundary measures and associated confidences are computed at potential vessel boundaries. A training technique is used to develop a mapping of this vector to a likelihood ratio that measures the “vesselness” at each pixel. Results comparing this vesselness measure to matched filters alone and to measures based on the Hessian of intensities have shown substantial improvements both qualitatively and quantitatively. Binary executables of the code are available at http://www.sofka.com/LRV.html.
This represents a possible route for fully automated tracing of retinal vessels. A typical tracing is shown in Fig. 12.6. This is a promising result, although it was difficult to avoid breaks in the tracings without adding spurious vessels. It may be possible to achieve a more acceptable result with an optimal choice of tracer sensitivity, and judicious preprocessing of the image to remove background variation and reduce noise.
12.6 Experimental Results and Discussion
The grades associated with the Messidor project database (and all other databases) are rather coarse. To explore changes in the vessels with the severity of retinopathy, we need a finer grading scheme. The number of detected microaneurysms has been
280 |
M. Iorga and G. Dougherty |
Fig. 12.6 The result of a fully automated tracing of retinal vessels using Sofka’s program on the same image as in Fig. 12.5. The traced vessels are shown in green, with the bifurcations subsequently colored red
used as a surrogate for the severity of retinopathy in its early stages [12, 13]. This information was available in the marked database (of 82 unique images) supplied by MINES ParisTech.
Figure 12.7 shows that the tortuosity (per cm) increases steadily with the number of manually identified microaneurysms from images in the marked database. The correlation coefficient (r = 0.4159, n = 46) corresponds to a probability, p, for the null hypothesis of 0.004.4 This suggests that tortuosity is related to microaneurysm count, and that the tortuosity increases with the severity of retinopathy.
The number of microaneurysms detected by the Waikato Automated Microaneurysm Detector did not match the numbers manually counted by experts for images from the marked database. Despite the normalization of the images prior to segmentation, the number of microaneurysms detected is still somewhat dependent on the threshold probability used in the Bayes classifier. Useful thresholds normally lie between 10−5 and 10−7 (Cree, private communication). The key is to select a threshold which is high enough to detect all the microaneurysms but small enough not to count noise as aneurysm, but this optimum threshold depends on the particular image. We investigated two methods for finding the optimum value of the threshold.
4http://faculty.vassar.edu/lowry/VassarStats.html.
12 Tortuosity as an Indicator of the Severity of Diabetic Retinopathy |
281 |
Fig. 12.7 The tortuosity/length for 46 vessels from the marked database. The best-fitted line is superimposed
For both methods, we plotted the number of microaneurysms detected in an image for a series of thresholds between 10−5 and 10−7. The higher the threshold, the smaller the number of aneurysms detected. We considered the optimum threshold to be a balance between the microaneurysm count changing too fast with threshold and it hardly changing at all with threshold. In our first method (method 1), we considered this be the value where the plot of microaneurysm counts was furthest from a straight line connecting the two extreme thresholds (Fig. 12.8). In our second method (method 2), we calculated the local curvature of datapoints (by calculating the local tortuosity using five points centered on the datapoint of interest) in the plot of microaneurysm counts vs. threshold, and considered the largest value to indicate the position of the optimum threshold.
We tested both methods on the marked database, for which we know the microaneurysm counts (viz., the “gold standard” counts from manual counting by three experts). Table 12.1 shows the number of microaneurysms in the images, assessed by the various methods. Probability thresholds of 10−5, 10−6, and 10−7 are included, although microaneurysm counts for 100 different thresholds were computed. The correlation coefficient for these fixed values with the “gold standard” values varied between 0.4696 and 0.6805, depending on the value of the threshold. Although it is possible to achieve a correlation of 0.6805 with a fixed threshold, it would be highly unlikely that this particular value of the threshold would be chosen a priori. The correlation coefficient for the number of microaneurysms found by method 1 with the “gold standard” counts was 0.6458, while the correlation coefficient for the number of microaneurysms found by method 2 with the “gold standard” counts was higher at 0.7303. (This corresponds to p < 0.001 (n = 82)
282 |
M. Iorga and G. Dougherty |
Fig. 12.8 (a) Image (MA Originaux 09, from the marked database) and (b) a plot showing how the optimum threshold is obtained from this image using method 1
for the null hypothesis). In the light of these findings, we consider that method 2 is particularly well suited for finding the optimum threshold for use with the Waikato Automated Microaneurysm Detector.
The microaneurysms in the Messidor project database are not so clearly delineated. We detected and counted them with the Waikato Automated Microaneurysm Detector using method 2 (Fig. 12.9). Only images with a microaneurysm count of five or greater were considered, since we had greatest confidence that these corresponded to actual microaneurysms. The correlation
12 Tortuosity as an Indicator of the Severity of Diabetic Retinopathy |
283 |
Table 12.1 The number of microaneurysms detected for the images of the marked database using various thresholds, and our methods 1 and 2, compared with the actual “gold standard” counts
|
“Gold |
|
|
|
|
|
|
standard” |
Threshold |
Threshold |
Threshold |
|
|
Image ID |
counts |
= 10−5 |
= 10−6 |
= 10−7 |
Method 1 |
Method 2 |
1 |
2 |
0 |
1 |
11 |
1 |
2 |
2 |
3 |
2 |
3 |
11 |
3 |
4 |
3 |
2 |
1 |
1 |
4 |
1 |
2 |
4 |
1 |
3 |
4 |
16 |
3 |
4 |
5 |
4 |
2 |
4 |
19 |
4 |
4 |
6 |
16 |
5 |
17 |
34 |
11 |
16 |
7 |
3 |
1 |
6 |
11 |
2 |
7 |
8 |
8 |
3 |
7 |
12 |
4 |
7 |
9 |
9 |
8 |
15 |
36 |
11 |
11 |
10 |
19 |
7 |
20 |
39 |
16 |
16 |
11 |
14 |
4 |
15 |
41 |
11 |
9 |
12 |
3 |
2 |
3 |
11 |
3 |
4 |
13 |
3 |
2 |
4 |
11 |
4 |
5 |
14 |
4 |
0 |
4 |
7 |
1 |
4 |
15 |
1 |
3 |
4 |
16 |
3 |
4 |
16 |
5 |
2 |
4 |
19 |
4 |
4 |
17 |
11 |
3 |
10 |
22 |
10 |
11 |
18 |
16 |
7 |
17 |
29 |
13 |
15 |
19 |
7 |
2 |
8 |
21 |
4 |
4 |
20 |
0 |
1 |
1 |
7 |
1 |
3 |
21 |
2 |
0 |
1 |
4 |
1 |
1 |
22 |
3 |
1 |
6 |
11 |
2 |
7 |
23 |
6 |
2 |
7 |
11 |
4 |
5 |
24 |
7 |
3 |
13 |
25 |
11 |
10 |
25 |
0 |
1 |
1 |
3 |
1 |
3 |
26 |
0 |
0 |
1 |
6 |
0 |
1 |
27 |
6 |
3 |
10 |
39 |
8 |
6 |
28 |
10 |
8 |
17 |
48 |
14 |
14 |
29 |
11 |
0 |
3 |
12 |
2 |
3 |
30 |
5 |
5 |
13 |
20 |
5 |
15 |
31 |
2 |
1 |
7 |
21 |
5 |
6 |
32 |
8 |
4 |
13 |
36 |
9 |
7 |
33 |
7 |
0 |
6 |
13 |
3 |
6 |
34 |
5 |
4 |
7 |
33 |
6 |
8 |
35 |
3 |
1 |
3 |
15 |
1 |
3 |
36 |
0 |
1 |
2 |
43 |
2 |
2 |
37 |
1 |
0 |
6 |
34 |
2 |
2 |
38 |
18 |
5 |
12 |
24 |
10 |
10 |
39 |
18 |
1 |
13 |
34 |
8 |
14 |
40 |
22 |
5 |
17 |
42 |
11 |
17 |
41 |
7 |
3 |
9 |
25 |
9 |
9 |
42 |
17 |
3 |
8 |
31 |
5 |
8 |
43 |
10 |
5 |
18 |
57 |
10 |
10 |
|
|
|
|
|
|
(continued) |
284 |
|
|
|
|
M. Iorga and G. Dougherty |
|
Table 12.1 |
(continued) |
|
|
|
|
|
|
“Gold |
|
|
|
|
|
|
standard” |
Threshold |
Threshold |
Threshold |
|
|
Image ID |
counts |
= 10−5 |
= 10−6 |
= 10−7 |
Method 1 |
Method 2 |
45a |
4 |
2 |
6 |
11 |
7 |
4 |
46 |
4 |
1 |
4 |
22 |
3 |
3 |
47 |
1 |
0 |
1 |
10 |
5 |
2 |
48 |
18 |
4 |
9 |
25 |
1 |
9 |
49 |
2 |
1 |
2 |
14 |
9 |
2 |
50 |
13 |
5 |
15 |
42 |
2 |
12 |
51 |
10 |
4 |
9 |
30 |
12 |
8 |
52 |
7 |
1 |
8 |
23 |
9 |
4 |
53 |
1 |
0 |
1 |
4 |
2 |
1 |
54 |
0 |
1 |
4 |
11 |
1 |
4 |
55 |
3 |
1 |
3 |
6 |
1 |
3 |
56 |
12 |
3 |
5 |
36 |
2 |
5 |
57 |
14 |
3 |
17 |
52 |
5 |
11 |
58 |
8 |
2 |
5 |
8 |
11 |
6 |
59 |
0 |
0 |
2 |
5 |
3 |
2 |
60 |
2 |
2 |
3 |
12 |
2 |
3 |
61 |
4 |
1 |
2 |
9 |
3 |
3 |
62 |
11 |
3 |
11 |
22 |
2 |
9 |
63 |
21 |
4 |
11 |
25 |
9 |
12 |
64 |
5 |
5 |
12 |
27 |
6 |
14 |
65 |
14 |
8 |
21 |
36 |
12 |
13 |
66 |
5 |
1 |
2 |
17 |
16 |
4 |
67 |
2 |
0 |
2 |
17 |
2 |
2 |
68 |
0 |
1 |
2 |
22 |
0 |
2 |
69 |
0 |
0 |
0 |
3 |
2 |
3 |
70 |
2 |
3 |
7 |
24 |
0 |
7 |
71 |
1 |
1 |
2 |
5 |
6 |
2 |
72 |
1 |
1 |
1 |
14 |
1 |
2 |
73 |
13 |
9 |
26 |
61 |
1 |
14 |
74 |
12 |
15 |
22 |
93 |
14 |
20 |
75 |
6 |
3 |
15 |
50 |
22 |
13 |
76 |
0 |
0 |
1 |
4 |
11 |
3 |
77 |
0 |
0 |
1 |
3 |
0 |
1 |
78 |
4 |
10 |
37 |
111 |
1 |
20 |
79 |
2 |
1 |
3 |
7 |
25 |
3 |
80 |
1 |
0 |
1 |
4 |
3 |
1 |
81 |
1 |
0 |
8 |
39 |
1 |
8 |
82 |
15 |
5 |
15 |
32 |
3 |
16 |
83 |
6 |
5 |
15 |
8 |
9 |
7 |
aNo ID 44 exists because the image is identical to ID 43
of microaneurysm count with tortuosity is not as strong (r = 0.2236 (n = 34), corresponding to a probability value for the null hypothesis of 0.1018) as with the marked database. Although we optimized the counting within the Waikato
12 Tortuosity as an Indicator of the Severity of Diabetic Retinopathy |
285 |
Fig. 12.9 The tortuosity/length for 34 vessels from the Messidor project database. The best-fitted line is superimposed
Fig. 12.10 The tortuosity of the first quintile length compared to the tortuosity of the entire vessel, for 46 vessels from the marked database
Automated Microaneurysm Detector, we do not expect the number of counts to be completely accurate. Despite this, the level of correlation supports our earlier finding with the marked database that tortuosity is related to microaneurysm count, and that tortuosity increases with severity of retinopathy.
A complicating factor in the use of vessel tortuosity as an indicator of the severity of retinopathy would be a change in tortuosity along an individual blood vessel. We measured the tortuosity of vessels along the quintiles of their length. Figure 12.10