Lifespan Trajectories of White Matter Changes in Rhesus Monkeys.

Progress in neurodevelopmental brain research has been achieved through the use of animal models. Such models not only help understanding biological changes that govern brain development, maturation and aging, but are also essential for identifying possible mechanisms of neurodevelopmental and age-related chronic disorders, and to evaluate possible interventions with potential relevance to human disease. Genetic relationship of rhesus monkeys to humans makes those animals a great candidate for such models. With the typical lifespan of 25 years, they undergo cognitive maturation and aging that is similar to this observed in humans. Quantitative structural neuroimaging has been proposed as one of the candidate in vivo biomarkers for tracking white matter brain maturation and aging. While lifespan trajectories of white matter changes have been mapped in humans, such knowledge is not available for nonhuman primates. Here, we analyze and model lifespan trajectories of white matter microstructure using in vivo diffusion imaging in a sample of 44 rhesus monkeys. We report quantitative parameters (including slopes and peaks) of lifespan trajectories for 8 individual white matter tracts. We show different trajectories for cellular and extracellular microstructural imaging components that are associated with white matter maturation and aging, and discuss similarities and differences between those in humans and rhesus monkeys, the importance of our findings, and future directions for the field. Significance Statement: Quantitative structural neuroimaging has been proposed as one of the candidate in vivo biomarkers for tracking brain maturation and aging. While lifespan trajectories of structural white matter changes have been mapped in humans, such knowledge is not available for rhesus monkeys. We present here results of the analysis and modeling of the lifespan trajectories of white matter microstructure using in vivo diffusion imaging in a sample of 44 rhesus monkeys (age 4-27). We report and anatomically map lifespan changes related to cellular and extracellular microstructural components that are associated with white matter maturation and aging.

Inter-site and inter-scanner diffusion MRI data harmonization.

We propose a novel method to harmonize diffusion MRI data acquired from multiple sites and scanners, which is imperative for joint analysis of the data to significantly increase sample size and statistical power of neuroimaging studies. Our method incorporates the following main novelties: i) we take into account the scanner-dependent spatial variability of the diffusion signal in different parts of the brain; ii) our method is independent of compartmental modeling of diffusion (e.g., tensor, and intra/extra cellular compartments) and the acquired signal itself is corrected for scanner related differences; and iii) inter-subject variability as measured by the coefficient of variation is maintained at each site. We represent the signal in a basis of spherical harmonics and compute several rotation invariant spherical harmonic features to estimate a region and tissue specific linear mapping between the signal from different sites (and scanners). We validate our method on diffusion data acquired from seven different sites (including two GE, three Philips, and two Siemens scanners) on a group of age-matched healthy subjects. Since the extracted rotation invariant spherical harmonic features depend on the accuracy of the brain parcellation provided by Freesurfer, we propose a feature based refinement of the original parcellation such that it better characterizes the anatomy and provides robust linear mappings to harmonize the dMRI data. We demonstrate the efficacy of our method by statistically comparing diffusion measures such as fractional anisotropy, mean diffusivity and generalized fractional anisotropy across multiple sites before and after data harmonization. We also show results using tract-based spatial statistics before and after harmonization for independent validation of the proposed methodology. Our experimental results demonstrate that, for nearly identical acquisition protocol across sites, scanner-specific differences can be accurately removed using the proposed method.

Multi-shell diffusion signal recovery from sparse measurements.

For accurate estimation of the ensemble average diffusion propagator (EAP), traditional multi-shell diffusion imaging (MSDI) approaches require acquisition of diffusion signals for a range of b-values. However, this makes the acquisition time too long for several types of patients, making it difficult to use in a clinical setting. In this work, we propose a new method for the reconstruction of diffusion signals in the entire q-space from highly undersampled sets of MSDI data, thus reducing the scan time significantly. In particular, to sparsely represent the diffusion signal over multiple q-shells, we propose a novel extension to the framework of spherical ridgelets by accurately modeling the monotonically decreasing radial component of the diffusion signal. Further, we enforce the reconstructed signal to have smooth spatial regularity in the brain, by minimizing the total variation (TV) norm. We combine these requirements into a novel cost function and derive an optimal solution using the Alternating Directions Method of Multipliers (ADMM) algorithm. We use a physical phantom data set with known fiber crossing angle of 45° to determine the optimal number of measurements (gradient directions and b-values) needed for accurate signal recovery. We compare our technique with a state-of-the-art sparse reconstruction method (i.e., the SHORE method of Cheng et al. (2010)) in terms of angular error in estimating the crossing angle, incorrect number of peaks detected, normalized mean squared error in signal recovery as well as error in estimating the return-to-origin probability (RTOP). Finally, we also demonstrate the behavior of the proposed technique on human in vivo data sets. Based on these experiments, we conclude that using the proposed algorithm, at least 60 measurements (spread over three b-value shells) are needed for proper recovery of MSDI data in the entire q-space.

Gray matter alterations in early aging: a diffusion magnetic resonance imaging study.

Many studies have observed altered neurofunctional and structural organization in the aging brain. These observations from functional neuroimaging studies show a shift in brain activity from the posterior to the anterior regions with aging (PASA model), as well as a decrease in cortical thickness, which is more pronounced in the frontal lobe followed by the parietal, occipital, and temporal lobes (retrogenesis model). However, very little work has been done using diffusion MRI (dMRI) with respect to examining the structural tissue alterations underlying these neurofunctional changes in the gray matter. Thus, for the first time, we propose to examine gray matter changes using diffusion MRI in the context of aging. In this work, we propose a novel dMRI based measure of gray matter "heterogeneity" that elucidates these functional and structural models (PASA and retrogenesis) of aging from the viewpoint of diffusion MRI. In a cohort of 85 subjects (all males, ages 15-55 years), we show very high correlation between age and "heterogeneity" (a measure of structural layout of tissue in a region-of-interest) in specific brain regions. We examine gray matter alterations by grouping brain regions into anatomical lobes as well as functional zones. Our findings from dMRI data connects the functional and structural domains and confirms the "retrogenesis" hypothesis of gray matter alterations while lending support to the neurofunctional PASA model of aging in addition to showing the preservation of paralimbic areas during healthy aging.

Sexually dimorphic white matter geometry abnormalities in adolescent onset schizophrenia.

The normal human brain is characterized by a pattern of gross anatomical asymmetry. This pattern, known as the "torque", is associated with a sexual dimorphism: The male brain tends to be more asymmetric than that of the female. This fact, along with well-known sex differences in brain development (faster in females) and onset of psychosis (earlier with worse outcome in males), has led to the theory that schizophrenia is a disorder in which sex-dependent abnormalities in the development of brain torque, the correlate of the capacity for language, cause alterations in interhemispheric connectivity, which are causally related to psychosis (Crow TJ, Paez P, Chance SE. 2007. Callosal misconnectivity and the sex difference in psychosis. Int Rev Psychiatry. 19(4):449-457.). To provide evidence toward this theory, we analyze the geometry of interhemispheric white matter connections in adolescent-onset schizophrenia, with a particular focus on sex, using a recently introduced framework for white matter geometry computation in diffusion tensor imaging data (Savadjiev P, Kindlmann GL, Bouix S, Shenton ME, Westin CF. 2010. Local white geometry from diffusion tensor gradients. Neuroimage. 49(4):3175-3186.). Our results reveal a pattern of sex-dependent white matter geometry abnormalities that conform to the predictions of Crow's torque theory and correlate with the severity of patients' symptoms. To the best of our knowledge, this is the first study to associate geometrical differences in white matter connectivity with torque in schizophrenia.

Human middle longitudinal fascicle: variations in patterns of anatomical connections.

Based on high-resolution diffusion tensor magnetic resonance imaging (DTI) tractographic analyses in 39 healthy adult subjects, we derived patterns of connections and measures of volume and biophysical parameters, such as fractional anisotropy (FA) for the human middle longitudinal fascicle (MdLF). Compared to previous studies, we found that the cortical connections of the MdLF in humans appear to go beyond the superior temporal (STG) and angular (AG) gyri, extending to the temporal pole (TP), superior parietal lobule (SPL), supramarginal gyrus, precuneus and the occipital lobe (including the cuneus and lateral occipital areas). Importantly, the MdLF showed a striking lateralized pattern with predominant connections between the TP, STG and AG on the left and TP, STG and SPL on the right hemisphere. In light of the results of the present study, and of the known functional role of the cortical areas interconnected by the MdLF, we suggested that this fiber pathway might be related to language, high order auditory association, visuospatial and attention functions.

DIFFEOMORPHIC POINT SET REGISTRATION USING NON-STATIONARY MIXTURE MODELS.

This paper investigates a diffeomorphic point-set registration based on non-stationary mixture models. The goal is to improve the non-linear registration of anatomical structures by representing each point as a general non-stationary kernel that provides information about the shape of that point. Our framework generalizes work done by others that use stationary models. We achieve this by integrating the shape at each point when calculating the point-set similarity and transforming it according to the calculated deformation. We also restrict the non-rigid transform to the space of symmetric diffeomorphisms. Our algorithm is validated in synthetic and human datasets in two different applications: fiber bundle and lung airways registration. Our results shows that non-stationary mixture models are superior to Gaussian mixture models and methods that do not take into account the shape of each point.

A review of magnetic resonance imaging and diffusion tensor imaging findings in mild traumatic brain injury.

Mild traumatic brain injury (mTBI), also referred to as concussion, remains a controversial diagnosis because the brain often appears quite normal on conventional computed tomography (CT) and magnetic resonance imaging (MRI) scans. Such conventional tools, however, do not adequately depict brain injury in mTBI because they are not sensitive to detecting diffuse axonal injuries (DAI), also described as traumatic axonal injuries (TAI), the major brain injuries in mTBI. Furthermore, for the 15 to 30 % of those diagnosed with mTBI on the basis of cognitive and clinical symptoms, i.e., the "miserable minority," the cognitive and physical symptoms do not resolve following the first 3 months post-injury. Instead, they persist, and in some cases lead to long-term disability. The explanation given for these chronic symptoms, i.e., postconcussive syndrome, particularly in cases where there is no discernible radiological evidence for brain injury, has led some to posit a psychogenic origin. Such attributions are made all the easier since both posttraumatic stress disorder (PTSD) and depression are frequently co-morbid with mTBI. The challenge is thus to use neuroimaging tools that are sensitive to DAI/TAI, such as diffusion tensor imaging (DTI), in order to detect brain injuries in mTBI. Of note here, recent advances in neuroimaging techniques, such as DTI, make it possible to characterize better extant brain abnormalities in mTBI. These advances may lead to the development of biomarkers of injury, as well as to staging of reorganization and reversal of white matter changes following injury, and to the ability to track and to characterize changes in brain injury over time. Such tools will likely be used in future research to evaluate treatment efficacy, given their enhanced sensitivity to alterations in the brain. In this article we review the incidence of mTBI and the importance of characterizing this patient population using objective radiological measures. Evidence is presented for detecting brain abnormalities in mTBI based on studies that use advanced neuroimaging techniques. Taken together, these findings suggest that more sensitive neuroimaging tools improve the detection of brain abnormalities (i.e., diagnosis) in mTBI. These tools will likely also provide important information relevant to outcome (prognosis), as well as play an important role in longitudinal studies that are needed to understand the dynamic nature of brain injury in mTBI. Additionally, summary tables of MRI and DTI findings are included. We believe that the enhanced sensitivity of newer and more advanced neuroimaging techniques for identifying areas of brain damage in mTBI will be important for documenting the biological basis of postconcussive symptoms, which are likely associated with subtle brain alterations, alterations that have heretofore gone undetected due to the lack of sensitivity of earlier neuroimaging techniques. Nonetheless, it is noteworthy to point out that detecting brain abnormalities in mTBI does not mean that other disorders of a more psychogenic origin are not co-morbid with mTBI and equally important to treat. They arguably are. The controversy of psychogenic versus physiogenic, however, is not productive because the psychogenic view does not carefully consider the limitations of conventional neuroimaging techniques in detecting subtle brain injuries in mTBI, and the physiogenic view does not carefully consider the fact that PTSD and depression, and other co-morbid conditions, may be present in those suffering from mTBI. Finally, we end with a discussion of future directions in research that will lead to the improved care of patients diagnosed with mTBI.

Sparse multi-shell diffusion imaging.

Diffusion magnetic resonance imaging (dMRI) is an important tool that allows non-invasive investigation of neural architecture of the brain. The data obtained from these in-vivo scans provides important information about the integrity and connectivity of neural fiber bundles in the brain. A multi-shell imaging (MSI) scan can be of great value in the study of several psychiatric and neurological disorders, yet its usability has been limited due to the long acquisition times required. A typical MSI scan involves acquiring a large number of gradient directions for the 2 (or more) spherical shells (several b-values), making the acquisition time significantly long for clinical application. In this work, we propose to use results from the theory of compressive sampling and determine the minimum number of gradient directions required to attain signal reconstruction similar to a traditional MSI scan. In particular, we propose a generalization of the single shell spherical ridgelets basis for sparse representation of multi shell signals. We demonstrate its efficacy on several synthetic and in-vivo data sets and perform quantitative comparisons with solid spherical harmonics based representation. Our preliminary results show that around 20-24 directions per shell are enough for robustly recovering the diffusion propagator.

A diffusion tensor imaging study of the anterior limb of the internal capsule in schizophrenia.

Frontal-subcortical cognitive and limbic feedback loops modulate higher cognitive functioning. The final step in these feedback loops is the thalamo-cortical projection through the anterior limb of the internal capsule (AL-IC). Using diffusion tensor imaging (DTI), we evaluated abnormalities in the AL-IC fiber tract in schizophrenia. Participants comprised 16 chronic schizophrenia patients and 19 male, normal controls, who were group matched for handedness, age, and parental socioeconomic status, and underwent DTI on a 1.5 Tesla GE system. We measured the diffusion indices, fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD), and manually segmented, based on FA maps, AL-IC volume, normalized for intracranial contents (ICC). The results showed a significant reduction in the ICC-corrected volume of the AL-IC in schizophrenia, but did not show diffusion measure group differences in the AL-IC in FA, MD, RD or AD. In addition, in the schizophrenia patients, AL-IC FA correlated positively with performance on measures of spatial and verbal declarative/episodic memory, and right AL-IC ICC-corrected volume correlated positively with more perseverative responses on the Wisconsin Card Sort Test (WCST). We found a reduction in AL-IC ICC-corrected volume in schizophrenia, without FA, MD, RD or AD group differences, implicating the presence of a structural abnormality in schizophrenia in this subcortical white matter region which contains important cognitive, and limbic feedback pathways that modulate prefrontal cortical function. Despite not demonstrating a group difference in FA, we found that AL-IC FA was a good predictor of spatial and verbal declarative/episodic memory performance in schizophrenia.

Sequential anisotropic multichannel Wiener filtering with Rician bias correction applied to 3D regularization of DWI data.

It has been shown that the tensor calculation is very sensitive to the presence of noise in the acquired images, yielding to very low quality Diffusion Tensor Images (DTI) data. Recent investigations have shown that the noise present in the Diffusion Weighted Images (DWI) causes bias effects on the DTI data which cannot be corrected if the noise characteristic is not taken into account. One possible solution is to increase the minimum number of acquired measurements (which is 7) to several tens (or even several hundreds). This has the disadvantage of increasing the acquisition time by one (or two) orders of magnitude, making the process inconvenient for a clinical setting. We here proposed a turn-around procedure for which the number of acquisitions is maintained but, the DWI data are filtered prior to determining the DTI. We show a significant reduction on the DTI bias by means of a simple and fast procedure which is based on linear filtering; well-known drawbacks of such filters are circumvented by means of anisotropic neighborhoods and sequential application of the filter itself. Information of the first order probability density function of the raw data, namely, the Rice distribution, is also included. Results are shown both for synthetic and real datasets. Some error measurements are determined in the synthetic experiments, showing how the proposed scheme is able to reduce them. It is worth noting a 50% increase in the linear component for real DTI data, meaning that the bias in the DTI is considerably reduced. A novel fiber smoothness measure is defined to evaluate the resulting tractography for real DWI data. Our findings show that after filtering, fibers are considerably smoother on the average. Execution times are very low as compared to other reported approaches which allows for a real-time implementation.

Diffusion tractography of the fornix in schizophrenia.

BACKGROUND: White matter fiber tracts, especially those interconnecting the frontal and temporal lobes, are likely implicated in pathophysiology of schizophrenia. Very few studies, however, have focused on the fornix, a compact bundle of white matter fibers, projecting from the hippocampus to the septum, anterior nucleus of the thalamus and the mamillary bodies. Diffusion Tensor Imaging (DTI), and a new post-processing method, fiber tractography, provides a unique opportunity to visualize and to quantify entire trajectories of fiber bundles, such as the fornix, in vivo. We applied these techniques to quantify fornix diffusion anisotropy in schizophrenia. METHODS: DTI images were used to evaluate the left and the right fornix in 36 male patients diagnosed with chronic schizophrenia and 35 male healthy individuals, group matched on age, parental socioeconomic status, and handedness. Regions of interest were drawn manually, blind to group membership, to guide tractography, and fractional anisotropy (FA), a measure of fiber integrity, was calculated and averaged over the entire tract for each subject. The Doors and People test (DPT) was used to evaluate visual and verbal memory, combined recall and combined recognition. RESULTS: Analysis of variance was performed and findings demonstrated a difference between patients with schizophrenia and controls for fornix FA (p=0.006). Protected post-hoc independent sample t-tests demonstrated a bilateral FA decrease in schizophrenia, compared with control subjects (left side: p=0.048; right side p=0.006). Higher fornix FA was statistically significantly correlated with DPT and measures of combined visual memory (r=0.554, p=0.026), combined verbal memory (r=0.647, p=0.007), combined recall (r=0.516, p=0.041), and combined recognition (r=0.710, p=0.002) for the control group. No such statistically significant correlations were found in the patient group. CONCLUSIONS: Our findings show the utility of applying DTI and tractography to study white matter fiber tracts in vivo in schizophrenia. Specifically, we observed a bilateral disruption in fornix integrity in schizophrenia, thus broadening our understanding of the pathophysiology of this disease.

Signal LMMSE estimation from multiple samples in MRI and DT-MRI.

A method to estimate the magnitude MR data from several noisy samples is presented. It is based on the Linear Minimum Mean Squared Error (LMMSE) estimator for the Rician noise model when several scanning repetitions are available. This method gives a closed-form analytical solution that takes into account the probability distribution of the data as well as the existing level of noise, showing a better performance than methods such as the average or the median.

A method for clustering white matter fiber tracts.

BACKGROUND/PURPOSE: Despite its potential for visualizing white matter fiber tracts in vivo, diffusion tensor tractography has found only limited applications in clinical research in which specific anatomic connections between distant regions need to be evaluated. We introduce a robust method for fiber clustering that guides the separation of anatomically distinct fiber tracts and enables further estimation of anatomic connectivity between distant brain regions. METHODS: Line scanning diffusion tensor images (LSDTI) were acquired on a 1.5T magnet. Regions of interest for several anatomically distinct fiber tracts were manually drawn; then, white matter tractography was performed by using the Runge-Kutta method to interpolate paths (fiber traces) following the major directions of diffusion, in which traces were seeded only within the defined regions of interest. Next, a fully automatic procedure was applied to fiber traces, grouping them according to a pairwise similarity function that takes into account the shapes of the fibers and their spatial locations. RESULTS: We demonstrated the ability of the clustering algorithm to separate several fiber tracts which are otherwise difficult to define (left and right fornix, uncinate fasciculus and inferior occipitofrontal fasciculus, and corpus callosum fibers). CONCLUSION: This method successfully delineates fiber tracts that can be further analyzed for clinical research purposes. Hypotheses regarding specific fiber connections and their abnormalities in various neuropsychiatric disorders can now be tested.

DTI and MTR abnormalities in schizophrenia: analysis of white matter integrity.

Diffusion tensor imaging (DTI) studies in schizophrenia demonstrate lower anisotropic diffusion within white matter due either to loss of coherence of white matter fiber tracts, to changes in the number and/or density of interconnecting fiber tracts, or to changes in myelination, although methodology as well as localization of such changes differ between studies. The aim of this study is to localize and to specify further DTI abnormalities in schizophrenia by combining DTI with magnetization transfer imaging (MTI), a technique sensitive to myelin and axonal alterations in order to increase specificity of DTI findings. 21 chronic schizophrenics and 26 controls were scanned using Line-Scan-Diffusion-Imaging and T1-weighted techniques with and without a saturation pulse (MT). Diffusion information was used to normalize co-registered maps of fractional anisotropy (FA) and magnetization transfer ratio (MTR) to a study-specific template, using the multi-channel daemon algorithm, designed specifically to deal with multidirectional tensor information. Diffusion anisotropy was decreased in schizophrenia in the following brain regions: the fornix, the corpus callosum, bilaterally in the cingulum bundle, bilaterally in the superior occipito-frontal fasciculus, bilaterally in the internal capsule, in the right inferior occipito-frontal fasciculus and the left arcuate fasciculus. MTR maps demonstrated changes in the corpus callosum, fornix, right internal capsule, and the superior occipito-frontal fasciculus bilaterally; however, no changes were noted in the anterior cingulum bundle, the left internal capsule, the arcuate fasciculus, or inferior occipito-frontal fasciculus. In addition, the right posterior cingulum bundle showed MTR but not FA changes in schizophrenia. These findings suggest that, while some of the diffusion abnormalities in schizophrenia are likely due to abnormal coherence, or organization of the fiber tracts, some of these abnormalities may, in fact, be attributed to or coincide with myelin/axonal disruption.

Processing and visualization for diffusion tensor MRI.

This paper presents processing and visualization techniques for Diffusion Tensor Magnetic Resonance Imaging (DT-MRI). In DT-MRI, each voxel is assigned a tensor that describes local water diffusion. The geometric nature of diffusion tensors enables us to quantitatively characterize the local structure in tissues such as bone, muscle, and white matter of the brain. This makes DT-MRI an interesting modality for image analysis. In this paper we present a novel analytical solution to the Stejskal-Tanner diffusion equation system whereby a dual tensor basis, derived from the diffusion sensitizing gradient configuration, eliminates the need to solve this equation for each voxel. We further describe decomposition of the diffusion tensor based on its symmetrical properties, which in turn describe the geometry of the diffusion ellipsoid. A simple anisotropy measure follows naturally from this analysis. We describe how the geometry or shape of the tensor can be visualized using a coloring scheme based on the derived shape measures. In addition, we demonstrate that human brain tensor data when filtered can effectively describe macrostructural diffusion, which is important in the assessment of fiber-tract organization. We also describe how white matter pathways can be monitored with the methods introduced in this paper. DT-MRI tractography is useful for demonstrating neural connectivity (in vivo) in healthy and diseased brain tissue.

Nonrigid registration of 3D tensor medical data.

New medical imaging modalities offering multi-valued data, such as phase contrast MRA and diffusion tensor MRI, require general representations for the development of automated algorithms. In this paper we propose a unified framework for the registration of medical volumetric multi-valued data using local matching. The paper extends the usual concept of similarity between two pieces of data to be matched, commonly used with scalar (intensity) data, to the general tensor case. Our approach to registration is based on a multiresolution scheme, where the deformation field estimated in a coarser level is propagated to provide an initial deformation in the next finer one. In each level, local matching of areas with a high degree of local structure and subsequent interpolation are performed. Consequently, we provide an algorithm to assess the amount of structure in generic multi-valued data by means of gradient and correlation computations. The interpolation step is carried out by means of the Kriging estimator, which provides a novel framework for the interpolation of sparse vector fields in medical applications. The feasibility of the approach is illustrated by results on synthetic and clinical data.

CURVES: curve evolution for vessel segmentation.

The vasculature is of utmost importance in neurosurgery. Direct visualization of images acquired with current imaging modalities, however, cannot provide a spatial representation of small vessels. These vessels, and their branches which show considerable variations, are most important in planning and performing neurosurgical procedures. In planning they provide information on where the lesion draws its blood supply and where it drains. During surgery the vessels serve as landmarks and guidelines to the lesion. The more minute the information is, the more precise the navigation and localization of computer guided procedures. Beyond neurosurgery and neurological study, vascular information is also crucial in cardiovascular surgery, diagnosis, and research. This paper addresses the problem of automatic segmentation of complicated curvilinear structures in three-dimensional imagery, with the primary application of segmenting vasculature in magnetic resonance angiography (MRA) images. The method presented is based on recent curve and surface evolution work in the computer vision community which models the object boundary as a manifold that evolves iteratively to minimize an energy criterion. This energy criterion is based both on intensity values in the image and on local smoothness properties of the object boundary, which is the vessel wall in this application. In particular, the method handles curves evolving in 3D, in contrast with previous work that has dealt with curves in 2D and surfaces in 3D. Results are presented on cerebral and aortic MRA data as well as lung computed tomography (CT) data.

Three-dimensional adaptive filtering in magnetic resonance angiography.

In order to enhance 3D image data from magnetic resonance angiography (MRA), a novel method based on the theory of multidimensional adaptive filtering has been developed. The purpose of the technique is to suppress image noise while enhancing important structures. The method is based on local structure estimation using six 3D orientation selective filters, followed by an adaptive filtering step controlled by the local structure information. The complete filtering procedure requires approximately 3 minutes of computational time on a standard workstation for a 256 x 256 x 64 data set. The method has been evaluated using a mathematical vessel model and in vivo MRA data (both phase contrast and time of flight (TOF)). 3D adaptive filtering results in a better delineation of small blood vessels and efficiently reduces the high-frequency noise. Depending on the data acquisition and the original data type, contrast-to-noise ratio (CNR) improvements of up to 179% (8.9 dB) were observed. 3D adaptive filtering may provide an alternative to prolonging the scan time or using contrast agents in MRA when the CNR is low.