Accelerated magnetic resonance imaging using the sparsity of multi-channel coil images

Joint estimation of coil sensitivities and output image (JSENSE) is a promising approach that improves the reconstruction of parallel magnetic resonance imaging (pMRI). However, when acceleration factor increases, the signal to noise ratio (SNR) of JSENSE reconstruction decreases as quickly as that of the conventional pMRI. Although sparse constraints have been used to improve the JSENSE reconstruction in recent years, these constraints only use the sparsity of the output image, which cannot fully exploit the prior information of pMRI. In this paper, we use the sparsity of coil images, instead of the output image, to exploit more prior information for JSENSE. Numerical simulation, phantom and in vivo experiments demonstrate that the proposed method has better performance than the SparseSENSE method and the constrained JSENSE method using the sparsity of the output image only.     

Normalization of brain magnetic resonance images using histogram even-order derivative analysis

The even-ordered (2nd, 4th and 6th) derivatives of a brain MRI histogram were used to calculate a characteristic value for white matter, which was used to normalize the image intensity scale. Simulated image histograms were used to estimate the methodological error as a function of noise level, and the optimum derivative order was determined for each image type studied (T1-, T2- and density-weighted). The algorithm yielded highly reproducible results when used in conjunction with a threshold-sensitive brain segmentation algorithm. It also proved insensitive to the presence of extra-cranial tissues. This method of histogram analysis could find utility in a variety of applications that demand robust intensity normalization including image registration, brain segmentation, tissue classification and spatial inhomogeneity correction.     

Spatially encoded NMR and the acquisition of 2D magnetic resonance images within a single scan

An approach that enables the acquisition of multidimensional NMR spectra within a single scan has been recently proposed and demonstrated. The present paper explores the applicability of such ultrafast acquisition schemes toward the collection of two-dimensional magnetic resonance imaging (2D MRI) data. It is shown that ideas enabling the application of these spatially encoded schemes within a spectroscopic setting, can be extended in a straightforward manner to pure imaging. Furthermore, the reliance of the original scheme on a spatial encoding and subsequent decoding of the evolution frequencies endows imaging applications with a greater simplicity and flexibility than their spectroscopic counterparts. The new methodology also offers the possibility of implementing the single-scan acquisition of 2D MRI images using sinusoidal gradients, without having to resort to subsequent interpolation procedures or non-linear sampling of the data. Theoretical derivations on the operational principles and imaging characteristics of a number of sequences based on these ideas are derived, and experimentally validated with a series of 2D MRI results collected on a variety of model phantom samples.     

Bilateral filtering of magnetic resonance phase images

High-pass filtering is required for the removal of background field inhomogeneities in magnetic resonance phase images. This high-pass filtering smooths across boundaries between areas with large differences in phase. The most prominent boundary is the surface of the brain where areas with large phase values inside the brain are located close to areas outside the brain where the phase is, on average, zero. Cortical areas, which are of great interest in brain MRI, are therefore often degraded by high-pass filtering. Here, we propose the use of the bilateral filter for the high-pass filtering step. The bilateral filter is essentially a Gaussian filter that stops smoothing at boundaries. We show that the bilateral filter improves image quality at the brain's surface, without sacrificing contrast within the brain.     

Unbiased segmentation of diffusion-weighted magnetic resonance images of the brain using iterative clustering

Segmentation of diffusion-weighted echo-planar imaging (DW-EPI) is challenging because of concerns regarding spatial resolution and distortion. Methods commonly used require manual input and often need thresholding measures to segment white matter (WM), gray matter (GM) and cerebrospinal fluid (CSF). This may introduce operator bias and misclassification error. When comparing patients with a diffuse disease process-such as multiple sclerosis (MS)--with healthy controls, although information from all images may be biased due to disease effect, this is more so if the data set employed to perform segmentation is also used as a measured outcome for the study, for example, fractional anisotropy maps. Presented in this work is an unbiased method for segmenting DW-EPI data sets using the b=0 and single-shot inversion recovery EPI into WM, GM and CSF. The method employs an iterative clustering technique to account for partial volume effects and signal variation caused by radiofrequency inhomogeneity. The technique is evaluated with both real and synthetic brain data and results compared with statistical parametric mapping (SPM02). With synthetic brain data, where a gold standard of segmentation exists, the presented method showed less misclassification compared to SPM02. The unbiased method proposed may provide a more accurate methodology of segmentation in the analysis of DWI-EPI images in conditions such as MS.     

Adaptive segmentation of magnetic resonance images with intensity inhomogeneity using level set method

It is a big challenge to segment magnetic resonance (MR) images with intensity inhomogeneity. The widely used segmentation algorithms are region based, which mostly rely on the intensity homogeneity, and could bring inaccurate results. In this paper, we propose a novel region-based active contour model in a variational level set formulation. Based on the fact that intensities in a relatively small local region are separable, a local intensity clustering criterion function is defined. Then, the local function is integrated around the neighborhood center to formulate a global intensity criterion function, which defines the energy term to drive the evolution of the active contour locally. Simultaneously, an intensity fitting term that drives the motion of the active contour globally is added to the energy. In order to segment the image fast and accurately, we utilize a coefficient to make the segmentation adaptive. Finally, the energy is incorporated into a level set formulation with a level set regularization term, and the energy minimization is conducted by a level set evolution process. Experiments on synthetic and real MR images show the effectiveness of our method.     

Characterization of trabecular bone structure from high-resolution magnetic resonance images using fuzzy logic

The purpose of this work was to apply fuzzy logic image processing techniques to characterize the trabecular bone structure with high-resolution magnetic resonance images. Fifteen ex vivo high-resolution magnetic resonance images of specimens of human radii at 1.5 T and 12 in vivo high-resolution magnetic resonance images of the calcanei of peri- and postmenopausal women at 3 T were obtained. Soft segmentation using fuzzy clustering was applied to MR data to obtain fuzzy bone volume fraction maps, which were then analyzed with three-dimensional (3D) fuzzy geometrical parameters and measures of fuzziness. Geometrical parameters included fuzzy perimeter and fuzzy compactness, while measures of fuzziness included linear index of fuzziness, quadratic index of fuzziness, logarithmic fuzzy entropy, and exponential fuzzy entropy. Fuzzy parameters were validated at 1.5 T with 3D structural parameters computed from microcomputed tomography images, which allow the observation of true trabecular bone structure and with apparent MR structural indexes at 1.5 T and 3 T. The validation was statistically performed with the Pearson correlation coefficient as well as with the Bland-Altman method. Bone volume fraction correlation values (r) were up to .99 (P<.001) with good agreements based on Bland-Altman analysis showing that fuzzy clustering is a valid technique to quantify this parameter. Measures of fuzziness also showed consistent correlations to trabecular number parameters (r>.85; P<.001) and good agreements based on Bland-Altman analysis, suggesting that the level of fuzziness in high-resolution magnetic resonance images could be related to the trabecular bone structure.     

Motion estimation of magnetic resonance cardiac images using the Wigner-Ville and hough transforms

Myocardial motion analysis and quantification is of utmost importance for analyzing contractile heart abnormalities and it can be a symptom of a coronary artery disease. A fundamental problem in processing sequences of images is the computation of the optical flow, which is an approximation of the real image motion. This paper presents a new algorithm for optical flow estimation based on a spatiotemporal-frequency (STF) approach. More specifically it relies on the computation of the Wigner-Ville distribution (WVD) and the Hough Transform (HT) of the motion sequences. The latter is a well-known line and shape detection method that is highly robust against incomplete data and noise. The rationale of using the HT in this context is that it provides a value of the displacement field from the STF representation. In addition, a probabilistic approach based on Gaussian mixtures has been implemented in order to improve the accuracy of the motion detection. Experimental results in the case of synthetic sequences are compared with an implementation of the variational technique for local and global motion estimation, where it is shown that the results are accurate and robust to noise degradations. Results obtained with real cardiac magnetic resonance images are presented. The text was submitted by the authors in English.     

3D reconstruction of the brain from magnetic resonance images using a connectivity algorithm

We present high resolution three dimensional (3D) connectivity, surface construction and display algorithms that detect, extract, and display the surface of a brain from contiguous magnetic resonance (MR) images. The algorithms identify the external brain surface and create a 3D image, showing the fissures and surface convolutions of the cerebral hemispheres, cerebellum, and brain stem. Images produced by these algorithms also show the morphology of other soft tissue boundaries such as the cerebral ventricular system and the skin of the patient. For the purposes of 3D reconstruction, our experiments show that T1 weighted images give better contrast between the surface of the brain and the cerebral spinal fluid than T2 weighted images. 3D reconstruction of MR data provides a non-invasive procedure for examination of the brain surface and other anatomical features.     

A soft-segmentation visualization scheme for magnetic resonance images

Prevalent visualization tools exploit gray value distribution in images through modified histogram equalization and matching technique, referred to as the window width/window level-based method, to improve visibility and enhance diagnostic value. The window width/window level tool is extensively used in magnetic resonance (MR) images to highlight tissue boundaries during image interpretation. However, the identification of different regions and distinct boundaries between them based on gray-level distribution and displayed intensity levels is extremely difficult because of the large dynamic range of tissue intensities inherent in MR images. We propose a soft-segmentation visualization scheme to generate pixel partitions from the histogram of MR image data using a connectionist approach and then generate selective visual depictions of pixel partitions using pseudo color based on an appropriate fuzzy membership function. By applying the display scheme in clinical examples in this study, we could demonstrate additional overlapping regions between distinct tissue types in healthy and diseased areas (in the brain) that could help improve the tissue characterization ability of MR images.