Rodrigo de Luis-Garcia1, Miguel Angel Tola-Arribas2, Claudio Delrieux3, and Carlos Alberola-Lopez1
1Universidad de Valladolid, Valladolid, Spain, 2Hospital Universitario Rio Hortega, Valladolid, Spain, 3Universidad Nacional del Sur, Bahia Blanca, Argentina
Synopsis
Simple global measures describing the complexity of the
white matter architecture can provide useful information when analyzing
diffusion MRI data, and can be even capable of finding statistical differences
between groups. We propose the use of the fractal dimension of the FA maps for
that purpose, and illustrate its potential on a dataset composed of elderly
subjects and patients from three different stages of Alzheimer’s disease.PURPOSE AND MOTIVATION
This
abstract is focused on the use of the fractal dimension as a global descriptor
of the white matter. The fractal nature of diffusion MRI images has been
studied before
1,2. However, to the best of our knowledge, fractal descriptors
have not been employed as markers of disease.
Simple
fractal measures, such as the fractal dimension as estimated with the box counting
technique, can provide simple yet powerful procedure to obtain information
about the complexity of the white matter architecture.
METHODS
Four
groups of subjects from an Alzheimer study were analyzed, containing a healthy
control group (group A: N = 17, age= 74.5 ± 3.5y), patients with mild cog-
nitive impairment (group B: N = 13, age= 76.3 ± 1.1y), patients with mild
Alzheimer’s disease (group C: N = 19, age= 76.1 ± 2.7y) and patients with
moderate Alzheimer’s disease (group D: N = 7, age= 76.6±1.4y). Differences in
age were not significant between the cohorts. Patients were diagnosed according
to NINCDS-ADRDA Alzheimer's Criteria.
Diffusion
weighted images were acquired in a GE Signa 1.5 T MRI unit at QDiagnóstica,
Valladolid, Spain. The parameters of the acquisition protocol were the
following: 25 gradient directions, one baseline volume, b = 1000 s/mm2, 1.015 ×
1.015 × 3 mm3 of voxel size,TR = 13,000 ms, TE = 85.5 ms, 256 × 256 matrix, NEX
= 2 and 39 slices covering the entire brain.
After
preprocessing, including the removal of non-brain structures such as the skull,
diffusion tensors were estimated using a least squares method [27]. From the
tensor volume, Fractional Anisotropy (FA), Mean Diffusivity (MD) and Radial
Diffusivity (RD) maps were computed.
The
FA maps were afterwards slightly eroded in order to remove possible outliers. From
them, binary maps can be obtained simply by thresholding the FA maps.
The
Hausdoff fractal dimension can describe how much a certain pattern changes when
the scale at which it is measured also changes. While the most simple objects
can have an integer Hausdoff dimension (1 for a line, 2 for a square, 3 for a
cube), more complex objects can have non-integer Hausdorff dimensions. Among
the many techniques for the calculation of fractal properties, the box counting
method, which approximates the Hausdorff fractal dimension, is the most
commonly employed. Using this method, the (3D) space is partitioned in equal
boxes of size r. Then, N(r) is the number of boxes of size r that contain at
least a non-zero voxel. The
estimation of the fractal dimension, FD, is performed by computing the slope of
N(r), when plotted in a double logarithmic scale.
The
notion of fractal dimension can also be extended to gray-level images. In this
case, N(r) is the mean value inside each box, instead of the number of non-zero
voxels.
FD
values for the binarized FA maps (using a threshold of FA=0.3) and for the
gray-level FA maps were computed for all subjects, together with mean values
over the white matter of the FA, MD and RD maps. A one-way Anova test was
performed to investigate whether the four groups belong to the same
distribution. When they did not, bilateral t-tests were applied to check for pairwise
differences.
RESULTS
Figure
2 collects the p-values corresponding to the Anova and t-tests carried out.
Significant differences were found for the FD over the binarized (FD 0.3) and
gray-level FA (FA gray) maps, while no significant differences were found using
the mean FA values. There is extensive literature indicating that MD and RD are
more powerful descriptors of the changes within the white matter in Alzheimer’s
disesase and, accordingly, significant differences were also found for the mean
values of these maps.
With regard to the pairwise comparisons, the FD 0.3
and, to a lesser extent, the FD gray showed a considerable capacity to
differenciate subjects at different stages of Alzheimer’s disease. Notably, both
FD measures found significant differences between groups C and D, while mean MD
and mean RD discovered sifnificant differences between groups B and C (RD was
close to statistical significance). Although further investigation is needed,
this is a possible indication of different mechanisms of neurodegeneration
taking place at different stages of Alzheimer’s disease and thus affecting
different diffusion properties as measured by diffusion MRI.
CONCLUSION
Fractal
dimension of FA maps is a simple yet powerful method for providing global
descriptors of the white matter architecture. In the case of a group study on
Alzheimer’s disease, FD was able to reveal significant differences between
subjects at different stages of the disease.
Acknowledgements
The authors acknowledge the Ministerio de Ciencia e Innovación of Spain for research
grant TEC2013-44194-P.References
1.
P. Katsaloulis, P. Verganelakis, A. Provata, Fractal dimension and lacunarity
of tractography images of the human brain, Fractals 17(02): 181-189, 2009.
2.
P. Katsaloulis, A. Ghosh, et al, Fractality in the neuron axonal topography of
the human brain based on 3-D diffusion MRI, The European Physical Journal B 85:
150, 2012.