Abstract
SUMMARY: Diffusion tensor imaging is a valuable measure in clinical settings to assess diagnosis and prognosis of neonatal brain development. However, obtaining reliable images is not straightforward because of the tissue characteristics of the neonatal brain and the high likelihood of motion artifacts. In this review, we present guidelines on how to acquire DTI data of the neonatal brain and recommend high-quality data acquisition and processing as an essential means to obtain accurate and robust parametric maps. Sudden head movements are problematic for DTI in neonates, and these may lead to incorrect values. We describe strategies to minimize the corrupting effects both in terms of acquisition (eg, more gradient directions) and postprocessing (eg, tensor estimation methods). In addition, tools are described that can help assess whether a dataset is of sufficient quality for further assessment.
ABBREVIATIONS:
- MD
- mean diffusivity
- FA
- fractional anisotropy
- RESTORE
- robust estimation of tensors by outlier rejection
- TBSS
- tract-based spatial statistics
Both premature birth and complications around term birth are risk factors for brain injury and subsequent neurodevelopmental impairment. However, this injury often remains “silent” long after the threshold of irreversible injury has been crossed. The most important challenge is to define early proxy measurements of long-term neurodevelopmental outcome that will enable intervention in the early stages of the still adaptive developing human brain. MR imaging is widely used to monitor development and injury of the neonate brain and to predict neurodevelopmental outcome1⇓⇓⇓⇓–6 (Fig 1). Advanced quantitative MR imaging techniques such as DTI7⇓⇓⇓–11 can detect subtle differences between normal and abnormal tissue. DTI has become invaluable for the assessment of brain development in preterm infants because it enables detection of white matter maturation during premyelination.1,3,12⇓–14 Several DTI studies have revealed that abnormal white matter maturation is correlated to neuromotor and neurocognitive performances to childhood.15⇓⇓–18 DTI measurements can therefore provide early biomarkers to be targeted in neuroprotective intervention trials.
Although diffusion MR imaging of the neonatal brain is gaining popularity, time constraints often impede acquisition of high quality data. The total time in the scanner should be kept to a minimum, not only for ethical reasons but also to prevent hemodynamic instability. Moreover, a longer acquisition time increases the chance of motion. The feasibility of an MR imaging study and the optimal settings for the MR imaging scan are determined by the available time in combination with MR hardware and the requirements for obtaining reliable and meaningful data.
From the information above, it may be clear that acquisition, processing, and interpretation of the images with this advanced MR imaging technique are not as straightforward as with the conventional MR imaging sequences. Recent reviews have covered the general concerns within the DTI pipeline that can severely influence the measured results and therefore the study outcome.11,19,20 Many of the specifics of diffusion MR imaging acquisition in adult brains are also applicable to imaging the neonatal brain. However, the latter presents challenges related to size, tissue composition, and motion and therefore requires specific acquisition settings. Targeted acquisition schemes and postprocessing methods are necessary to account for motion and to obtain reliable DTI parameter maps.
In this review, we present guidelines on how to acquire diffusion MR imaging data of the neonatal brain. Furthermore, we recommend high-quality data acquisition and processing as essential means to obtain accurate and robust parametric maps. We will first describe some general topics related to DTI acquisition. Next, we broadly describe the specific concerns and acquisition settings related to DTI in the neonate. Finally, we discuss the main steps related to DTI acquisition and processing and the ways they influence the neonatal DTI data quality.
DTI of the Neonate
Imaging the neonatal brain is more difficult than imaging the adult brain. The neonate head is much smaller, and the brain is developing at a fast rate, has a high water content, and is unmyelinated.21 MR imaging acquisition settings therefore must be adjusted to neonatal MR imaging relaxation times and diffusion coefficients. Additionally, there is a greater likelihood of motion artifacts in the images because neonates have higher heart and breathing rates and will not lie motionless because sedation is usually not given in this situation. Fig 2 shows examples of motion artifacts. Thus, the optimal acquisition strategies established for the adult brain22 are not necessarily applicable to the neonatal brain because diffusivities, anisotropy, T2, and SNR will differ. Additionally, optimal acquisition depends on multiple, interrelated parameters, including but not limited to, TR, TE, acquisition matrix, field of view, section thickness, acquisition time, number of averages, maximum b-value, number of b-values, number of non-DWIs, and number of gradient directions.
Setup for Imaging Neonates
Neonatal MR imaging is challenging because many neonates are receiving respiratory support and are vulnerable to hemodynamic instability.23⇓–25 A recent study showed that despite of the use of a dedicated guideline, adverse events during neonatal MR imaging scanning were common.26 Because of the increased application of neonatal MR imaging, the provision of a safe environment during the MR imaging procedures is of great concern. An MR-compatible incubator with a specialized neonatal head coil allows the infant to be kept stable and safe during transport and MR imaging scan. A small-size neonatal head coil typically offers better SNR than the larger adult head coils.27 With improved SNR, the scan time can be reduced and/or the spatial resolution for the conventional scans can be improved.
Two common strategies are used to limit motion. One is sedation, which in neonates can compromise breathing and therefore must be performed with great caution. The second strategy is promoting sleep by placing the infant in the incubator in a comfortable and secure way.28 This will generally help, but it takes some practice to comfort the child and prevent it from waking up when imaging starts.29 The infant's perception of acquisition noise can be reduced with the use of moldable earplugs, neonatal earmuffs, and/or an acoustic hood.30
Resolution, TR, TE
A striking difference between the neonatal and adult brain is the inverse contrast between white and gray matter on both the T1- and T2-weighted images.31 This is caused by the high water content and the unmyelinated white matter of the neonate. T1 and T2 relaxation times quickly decrease over the first year of life, and the exact pattern is different for the different brain structures. Typical T2 relaxation times are approximately 300 ms at <30 weeks' gestational age, approximately 200 ms at term age,32 and approximately 100 ms at adulthood.33 T1 relaxation times for white matter are approximately 1600 ms at term age and approximately 500 ms at 2 years of age.31 For DTI measurements, a longer TR can be needed to reduce the effect of the high T1 relaxation time (ie, saturation and therefore reduced SNR and/or bias to other proton pools). Fortunately, the long T2 relaxation time does favor the DTI measurement because there is more signal left at a similar TE. However, the DTI estimates for different ages or brain structures must be interpreted carefully because they can be influenced by the varying T1 and T2 relaxation times, possibly caused by altered relative contributions of the different water compartments (eg, intracellular versus extracellular).
All structures in the neonatal brain are smaller than in the adult brain, thereby requiring a higher resolution. Typical in-plane resolutions range from 0.6 × 0.6 mm2 to 2 × 2 mm2, with a section thickness of 1.9–3.0 mm (Table 1), which is still relatively large; therefore images should be interpreted cautiously.34 Additionally, the in-plane resolution is often obtained with a low acquisition matrix that is zero-filled to a higher matrix; the images therefore are smoother and boundaries are less sharp. The imaging resolution and the partial volume effects determine what structures can be resolved. For fiber-tracking purposes, the resolution is ideally isotropic and contiguous sections are needed without section gap. The use of isotropic voxels prevents preferential averaging of fiber orientation along a certain axis.11
In addition, the resolution and SNR are related because higher resolution means lower SNR. Therefore, SNR constraints can limit the resolution. Higher resolutions are obtained by decreasing the field of view or increasing the matrix size. However, both options have their limitations because the field of view must fit the whole image volume, and higher matrix size introduces artifacts and increases scan time.
General DTI Acquisition Considerations
Generally, a single-shot EPI sequence allows fast acquisition but, unfortunately, is prone to artifacts that affect the data acquisition and the processing.11,19,20,35 The most critical artifacts are image distortions caused by magnetic susceptibility effects and eddy currents. Among the solutions proposed to reduce the susceptibility artifacts are parallel imaging and B0 field correction. Eddy currents can be reduced by applying bipolar gradients and dual-echo or twice-refocused sequences.36 For a more comprehensive and detailed description of these artifacts and other problems, the authors recommend the recent review by Jones and Cercignani19 or other reviews regarding DTI.11,20,35
Diffusion MR imaging is intrinsically a low SNR technique, and the SNR can significantly influence the reliability of each tensor estimate. For example, in adult brains, low SNR levels reduce the accuracy of the diffusion estimates (eg, increased FA) and decrease precision (eg, larger standard deviations).37 The SNR dependence is influenced by the underlying diffusivities and FA values and by the acquisition settings such as b-value.37,38
The magnetic field strength affects the quality of the images both positively and negatively. At 3T, the SNR is generally higher than at 1.5T. However, the images acquired at 3T are more sensitive to susceptibility artifacts, which are prone around air-tissue boundaries. Therefore, the field strength must be considered carefully when setting up a protocol.
The performance of the gradient systems is an important determinant of DTI image quality. Stronger gradients allow for higher diffusion weighting within a shorter time, thereby reducing TE and improving SNR and reducing artifacts. Gradient systems with advanced eddy-current compensation are preferred because these limit image distortions.
Hardware and software of the MR imaging systems differ between manufacturers and are constantly being upgraded. Therefore, imaging protocols cannot always be exchanged between sites, and what is achievable at one scanner might not be possible at another. Additionally, even equally equipped MR imaging systems can introduce a bias.39
The coil should adequately cover the volume of interest and exhibit sufficient SNR. Preferably, a multichannel coil should be used with the option to use parallel imaging. Parallel imaging techniques can shorten the EPI readout, thereby reducing imaging artifacts and TE.40
Diffusion Weighting
The b-value indicates the amount of diffusion weighting that is applied. The optimal b-value depends on the tensor information of interest but is generally approximately 1.09/ADC.22 For clinical adult brain DTI, a b-value of 1000 s/mm2 is common practice, which is a concession between a longer optimal b-value and the need for shorter echo times to ensure enough SNR. As mentioned, the neonatal brain contains more water and is less myelineated than the adult brain. Therefore, the neonatal MD is higher than the adult MD (2.0 × 10−3 mm2/s versus 0.7 × 10−3 mm2/s) and the optimal b-value should be lower to accurately estimate the diffusion tensor. Currently, there is no consensus on the b-value in neonates, and b-values range between 600–1100 s/mm2 (Table 1). Simulation studies are needed to determine the optimal b-value, just as has been performed for adult brains.
At least 6 diffusion gradient directions and 1 non-DWI are needed to estimate the diffusion tensor. However, more directions increase the accuracy and precision of the diffusion estimates. Tensor estimates with 30 directions are statistically rotationally invariant.41 Increasing the number of directions also entails increasing the number of b = 0 acquisitions because 1 b = 0 image is needed for every 8–10 DWIs.22 The optimal distribution of the gradients is a uniform distribution along the surface of a sphere,22 which minimizes the orientation dependency of the tensor estimates. Additionally, the order of the gradients ideally should be uniformly distributed in case the acquisition is corrupted by motion or interrupted acquisition.42,43 Dubois et al43 have proposed optimal diffusion gradient orientation schemes that allow calculation of the diffusion tensor estimates with the use of partial datasets.
The higher likelihood of motion in the neonates requires a higher number of gradient directions. The relative contribution of a corrupted image is less if more gradient directions are acquired. Therefore, it may be preferable to increase the number of gradient directions at the cost of reducing TR (and thereby SNR) to keep the total scan time approximately equal. In addition, we recommend the use of extra b = 0 images because these images also can be affected by intersection and intrasection motion, and good-quality b = 0 images are needed to correctly calculate the tensor estimates.
Motion
Subject-related motion can affect the resulting parameter maps. The motion can be evident between images or within a single image. During the long acquisition, adult subjects tend to slowly move their heads even if they are instructed to lay still and the head is secured by padding.44 A slow displacement and/or rotation causes a misalignment between images. This can be adjusted for by image registration. The motion within a single image is caused by cardiac pulsation or sudden large amplitude subject-related motion. The latter occurs rarely in adults (1.4%)45 but is a severe problem in neonates. Fig 2 shows the effects of severe motion, cardiac pulsation, and the large differences in signal intensities throughout the image volume.
Cardiac pulsation causes both a local deformation of the tissue and additional signal loss. The tensor estimates are then less reliable. Because neonates have much faster heart rates than adults, the likelihood that the diffusion sensitization is obtained during a pulsation is larger, leading to a higher likelihood of artifacts. The effects of cardiac pulsation can be reduced by triggering the acquisition on the cardiac cycle; however, scan times become longer because the effective TR should be at least 5 times the T1 of the tissue. For neonates, cardiac gating seems possible only by monitoring the ECG because the delay time for the pulse wave on the pulse oximeter is too long (249 ms for adults46) compared with the heart beat (approximately every 350 ms).
Subject motion is a major concern for neonatal imaging because neonates tend to move about even when sedated. Intrasection motion from a sudden shake of the head may corrupt the image and result in miscalculation of the DTI estimates. A typical motion-corrupted image is characterized by severe signal loss; this is caused by tissue displacement during the diffusion encoding, which also causes signal loss next to the dephasing caused by diffusion. This severe signal loss cannot be recovered, and an apparent high displacement would be calculated. Improper dealing with the motion artifacts may result in biased MD and high FA values.
Only a few studies have reported data rejection on the grounds of severe artifacts. Rejection was necessary in 15–60% of the subjects,15,47⇓–49 and the incidence of less noticeable artifacts will be larger. In a pilot study of 27 preterm neonates, we investigated the statistical outliers of the tensor regression (see “DTI Quality and Pathology” section) and showed that data of 60% of the subjects were corrupted by motion (defined as >10 sections with >30% outliers) (Table 2).50 This high occurrence of motion can have devastating results. Unfortunately, the typical tensor regression method used by MR imaging vendors and commercial processing tools is insufficient with corrupted data, and it appears that many DTI users are unaware of the motion-related problems.
The solutions to prevent corrupted images are limited. Obviously, the movement should be minimized a priori by placing the child comfortably and supporting the head with pads. The acquisition strategies to reduce the signal loss are limited because the sudden gross motion occurs during the time the diffusion encoding is performed. Therefore, navigator echoes or altered readout sequences appear not to be useful. However, shortening the diffusion time by stronger gradients or lower b-value might reduce the effects. Fortunately, there are strategies to decrease the detrimental effects of corrupted images. These include the oversampling of gradient directions, the removal of corrupted images, and/or the use of more advanced tensor regression methods.
Temperature
Another issue one must be aware of is that diffusion is temperature-dependent. A 1° difference in temperature leads to a 2–3% difference in ADC.51,52 It will therefore be necessary to maintain the child at a constant temperature during scanning. Also, consistent core temperatures are necessary in studying group differences or performing a longitudinal study because the differences in MD or FA between ages, groups, and therapies are small.
Processing
The basic steps in data postprocessing are motion/distortion correction, estimation of the diffusion tensor, and computation of the parametric maps. The data should be inspected before and after each step to ensure good image quality and absence of artifacts. Postprocessing results in data reduction, and it can be difficult to spot erroneous results on the parametric maps.
Motion/Distortion Correction
Registration and correction of the images is necessary to adjust for slight variation in brain position and to correct for the eddy current. A quick check for misregistration between images is to examine the FA maps for rims of high anisotropy at the edges of the brain or at the interface between CSF and white matter.
Because most available registration tools are optimized for the contrast in adult brain, one should check whether the registration software is adequate for the neonatal brain, in view of its different contrast. In addition, most tools perform their registration on the whole image volume and it is therefore logical that section drop-outs caused by motion severely alter the registration, especially in the case of affine registration. A final step is to also correct the diffusion encoding gradients with the same rotational correction as the images.44 Neglecting to do so will induce a mismatch between the actual diffusion weighting and the expected diffusion weighting and therefore result in errors in the estimates of the diffusion tensor and the fiber orientation.
It is of paramount importance to ensure that motion correction, which is often necessary in neonatal scanning, does not negatively affect the DTI dataset. Intersection movement is problematic for volume registration, and the b = 0 volume should be carefully checked on misalignment between sections. Large signal drop-outs are problematic with through-plane motion correction because thereafter, the incorrect signal intensities affect 2 sections, resulting in possible incorrect tensor estimations. Zhou et al47 investigated the use of local texture features to identify and reject outlier images automatically before estimating the diffusion tensor. Their method is fast and removes sections that are corrupted and cannot be used for tensor estimation, resulting in more accurate data. Implementation of fast detection techniques that are based on image characteristics will improve data quality and might offer possibilities to repeat corrupted gradient directions while scanning.53
Estimation of the Tensor
Several tensor estimation methods have been developed to estimate the diffusion tensor. Each of the approaches is based on different principles, and generally speed and accuracy are inversely proportional: methods range from fast but less accurate, to slow but more accurate.9,49,54 The main problem in estimating the tensor is the presence of data outliers caused by motion, pulsation, artifacts, noise, and so on. Depending on the method used, these outliers can have a significant impact on the resulting eigenvalues and eigenvectors. The linear least-squares method is commonly used by vendors and is the default setting for commonly used DTI postprocessing software such as FSL (http://www.fmrib.ox.ac.uk/fsl)55; however, this method proved to be the least appropriate to estimate the tensor because the method assumptions are restrictive and physically implausible.56 Fig 3 shows the effect of the tensor estimation method for a dataset with few outliers and a dataset corrupted as the result of motion. The difference between the methods is clearly visible, with the ordinary least-squares method providing obviously erroneous results, whereas the weighted least-squares method, which only takes a fraction more computational time, provides more accurate results.
Therefore, for neonatal DTI data, with its high likelihood of corrupted data, we must use proper tensor estimation procedures to obtain reliable diagnostic images or research data. RESTORE54 is a method that detects and removes outlier data in the tensor estimation. This method requires additional gradient encoding directions to eliminate erroneous data. For neonatal DTI measurements, Morris et al49 investigated the use of removing outliers based on the RESTORE algorithm before motion correction because the latter can result in averaging outliers with uncorrupted voxels. With this new method, they found an increased sensitivity to outliers, which is important as outliers correlated strongly with subject movement.
Quality Assurance
It is of paramount importance to ensure that all the steps from data acquisition to data analysis are correct and that the parameter maps are of sufficient quality before analyzing the data and drawing conclusions. Good-quality data are of paramount importance for the neonatal population because the chances of artifacts and corrupted data are larger than for the average adult population, and artifacts will bias the quantification and increase the variation. Data of sufficient quality will provide quantitative measures with low variation.
Quality assurance begins by carefully checking the b = 0 images for accuracy and absence of artifacts and intersection movement. Intersection movement causes severe problems with the image registration, resulting in erroneous results; therefore, those image volumes must be excluded from further analysis and hence our previous suggestion to increase the number of b = 0 volumes. Thereafter, the DWIs are inspected for large signal dropouts, which are easily observed by through-plane projections (Fig 2C). Image misalignments caused by motion or eddy currents can be visualized by a high FA around the rim of the brain or by calculating the standard deviations across the DW images, in which large SD at the rims indicates misregistration. This is also a good method to check whether images were registered accurately or if registration issues persist or even are introduced. Artifacts in the data can also be spotted by locating the areas where FA is >1. By definition, this should not be plausible but occurs when the DW signal intensity is larger than the b = 0 signal intensity.
The tensor estimation model also enables possibilities to spot artifacts by means of the residuals of the tensor fit and the calculation of data outliers. Residuals represent the difference between the fit and the data signal, with high residuals indicating either a poor estimation model or underlying data artifacts. The average residual maps as depicted in Fig 4 are helpful to detect and locate artifacts, assuming that high residuals are caused by artifacts and that good image data result in a homogeneous residuals map.20 High residuals are generally found in areas with partial volume effects (rims) and in regions with low signal intensities,20 and therefore one should extra carefully interpret results from those regions.
In addition, the residuals are used to detect signal outliers. The location and extent of outliers is a functional tool to investigate the DTI quality. The example in Fig 5 shows that 2 DTI volumes have a high percentage of outliers, which corresponds to corrupted images. We have found that determining the volumes and sections with high percentage outliers is an excellent method to assess data corruption caused by gross motion.50 However, labeling a point as an outlier also is dependent on the presence or absence of other (larger) outliers. Indeed, one expects a more or less even distribution of outliers with high-quality data, whereas corrupted sections cause a significant increase in outliers (Fig 6). Further research is needed to determine the relation between the distribution of outliers and the resulting parameter maps, the needed processing steps, and suitability of a dataset for further analysis.
Data Analysis
After the computation of the tensor estimates, the parametric maps can be used for further data analysis to assess the white matter integrity or structure. Commonly used techniques are region of interest selection, voxel-based analysis, tract-based spatial statistics, and fiber tractography. Region of interest selection is used to compare white matter properties within a specific region of the brain. The lower number of comparisons between subjects increases the statistical power, and small differences between subjects may be found. The technique is also independent of individual variations in shape and size of the brain. Major limitations include reliance on the accuracy of the region of interest placements, the likelihood of missing differences in regions that are not included, the time-consuming nature when the regions of interest are manually drawn, and the observer dependency.
Whole-brain analysis can be performed with voxel-based analyses, in which all voxels are analyzed in terms of a tensor estimate. Because of individual variation in size and shape of the brain, all individual brains must first be aligned to a template (normalization) before they can be compared.57,58 Normalization is the main challenge in voxel-based analysis because inaccurate normalization results in incorrect comparison of individuals.59 This is particularly problematic in patients with cerebral injuries where normalization results in distortion of the brains. TBSS overcomes this limitation by aligning the FA images from multiple subjects to a “mean FA skeleton.” With this technique, only voxels that are at the center of tracts common to all the individuals in the population are included. Because voxels with poor alignment are excluded, one of the limitations of this strategy is that it may miss variations in the periphery of white matter tracts. Other limitations and concerns related to TBSS have been addressed by others.60,61
Fiber tractography enables delineation and comparison of the orientation and direction of a specific white matter tract between individuals. This method is independent of variation in brain organization. In the deterministic approach,62 the algorithm moves in the direction of the principal eigenvector (ε1); this is assumed to be parallel to the dominant fiber orientation in each voxel. The main limitation of deterministic tractography is uncertainty about the reliability of a reconstructed trajectory. Probabilistic tractography takes the uncertainty in ε1 into account by propagating a distribution of possible orientations from the seed point.63 Although probabilistic tractography seems more reliable, it remains dependent on the placement of the seed points, which is manually done.3,12,35,64⇓–66
Although these analytic measurements provide in vivo information about the orientation, organization, and microstructural properties of the cerebral white matter, DTI analysis remains a proxy technique and not a direct visualization of the cerebral white matter. Moreover, because tensor estimations are subject to error, we emphasize the importance of appropriate acquisition settings with regard to gestational age, a quality check before postprocessing, reliable postprocessing techniques, and good knowledge of neuro anatomy.35
DTI Quality and Pathology
Investigation of DTI results with pathology brings out the importance of accurate data acquisition, processing, and data interpretation. Fig 7 shows 2 examples of DTI measurements: a patient with a middle cerebral artery infarct and a patient diagnosed with massive cerebral infarction, both scanned within days after the injury. Large portions of the brain are affected, and these consequently have lower MD values as the result of edema. These datasets appear to be of lesser quality in terms of outliers and residuals than the nonpathologic datasets without motion corruption. However, a closer look reveals that the outliers are randomly distributed in time and that the residuals are high in the areas with low MD values. This is expected because the diffusivities in the pathologic areas are lower and therefore the b-value might not be optimal to accurately determine MD and FA values in these areas. Although data quality assurance is not directly a vital issue for clinical practice in which the main goal is to visualize the areas with altered MD and FA values, it is imperative for the analysis of the DTI values to exclude any bias in DTI indices.
Conclusions
High-quality, quantitative data are essential to ensure reliable and meaningful imaging findings. Both data quality and analysis will improve with the use of protocols and image-processing tools specifically designed for neonatal imaging, seeing that tissue composition and occurrence of artifacts in neonates are significantly different from those in the adult population. Optimized acquisition and targeted processing will increase the diagnostic and prognostic value of the scans but will also considerably improve intervention or outcome studies because a reduced data spread will lead to stronger correlations, will require fewer subjects, and will allow earlier determination of valuable intervention methods.
REFERENCES
- 1.↵
- 2.↵
- 3.↵
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.↵
- 64.↵
- 65.↵
- 66.↵
- 67.
- 68.
- 69.
- 70.
- 71.
- 72.
- 73.
- 74.
- 75.
- 76.
- 77.
- 78.
- 79.
- 80.
- Received October 3, 2012.
- Accepted after revision November 15, 2012.
- © 2013 by American Journal of Neuroradiology