Measurements of magnetic field variations in the human brain using a 3D-FT multiple gradient echo technique

EEG-Informed fMRI: A Review of Data Analysis Methods

The simultaneous acquisition of electroencephalography (EEG) with functional magnetic resonance imaging (fMRI) is a very promising non-invasive technique for the study of human brain function. Despite continuous improvements, it remains a challenging technique, and a standard methodology for data analysis is yet to be established. Here we review the methodologies that are currently available to address the challenges at each step of the data analysis pipeline. We start by surveying methods for pre-processing both EEG and fMRI data. On the EEG side, we focus on the correction for several MR-induced artifacts, particularly the gradient and pulse artifacts, as well as other sources of EEG artifacts. On the fMRI side, we consider image artifacts induced by the presence of EEG hardware inside the MR scanner, and the contamination of the fMRI signal by physiological noise of non-neuronal origin, including a review of several approaches to model and remove it. We then provide an overview of the approaches specifically employed for the integration of EEG and fMRI when using EEG to predict the blood oxygenation level dependent (BOLD) fMRI signal, the so-called EEG-informed fMRI integration strategy, the most commonly used strategy in EEG-fMRI research. Finally, we systematically review methods used for the extraction of EEG features reflecting neuronal phenomena of interest.

An illustrated comparison of processing methods for phase MRI and QSM: removal of background field contributions from sources outside the region of interest

The elimination of so-called background fields is an essential step in phase MRI and quantitative susceptibility mapping (QSM). Background fields, which are caused by sources outside the region of interest (ROI), are often one to two orders of magnitude stronger than tissue-related field variations from within the ROI, hampering quantitative interpretation of field maps. This paper reviews the current literature on background elimination algorithms for QSM and provides insights into similarities and differences between the many algorithms proposed. We discuss the basic theoretical foundations and derive fundamental limitations of background field elimination. Copyright © 2016 John Wiley & Sons, Ltd.

Foundations of MRI phase imaging and processing for Quantitative Susceptibility Mapping (QSM)

Quantitative Susceptibility Mapping (QSM) is a novel MRI based technique that relies on estimates of the magnetic field distribution in the tissue under examination. Several sophisticated data processing steps are required to extract the magnetic field distribution from raw MRI phase measurements. The objective of this review article is to provide a general overview and to discuss several underlying assumptions and limitations of the pre-processing steps that need to be applied to MRI phase data before the final field-to-source inversion can be performed. Beginning with the fundamental relation between MRI signal and tissue magnetic susceptibility this review covers the reconstruction of magnetic field maps from multi-channel phase images, background field correction, and provides an overview of state of the art QSM solution strategies.

Accelerated echo planar J-resolved spectroscopic imaging in prostate cancer: a pilot validation of non-linear reconstruction using total variation and maximum entropy

The overlap of metabolites is a major limitation in one-dimensional (1D) spectral-based single-voxel MRS and multivoxel-based MRSI. By combining echo planar spectroscopic imaging (EPSI) with a two-dimensional (2D) J-resolved spectroscopic (JPRESS) sequence, 2D spectra can be recorded in multiple locations in a single slice of prostate using four-dimensional (4D) echo planar J-resolved spectroscopic imaging (EP-JRESI). The goal of the present work was to validate two different non-linear reconstruction methods independently using compressed sensing-based 4D EP-JRESI in prostate cancer (PCa): maximum entropy (MaxEnt) and total variation (TV). Twenty-two patients with PCa with a mean age of 63.8 years (range, 46-79 years) were investigated in this study. A 4D non-uniformly undersampled (NUS) EP-JRESI sequence was implemented on a Siemens 3-T MRI scanner. The NUS data were reconstructed using two non-linear reconstruction methods, namely MaxEnt and TV. Using both TV and MaxEnt reconstruction methods, the following observations were made in cancerous compared with non-cancerous locations: (i) higher mean (choline + creatine)/citrate metabolite ratios; (ii) increased levels of (choline + creatine)/spermine and (choline + creatine)/myo-inositol; and (iii) decreased levels of (choline + creatine)/(glutamine + glutamate). We have shown that it is possible to accelerate the 4D EP-JRESI sequence by four times and that the data can be reliably reconstructed using the TV and MaxEnt methods. The total acquisition duration was less than 13 min and we were able to detect and quantify several metabolites.

Quantitative T2, T2*, and T2' MR imaging in patients with ischemic leukoaraiosis might detect microstructural changes and cortical hypoxia

INTRODUCTION

Quantitative MRI with T2, T2*, and T2' mapping has been shown to non-invasively depict microstructural changes (T2) and oxygenation status (T2* and T2') that are invisible on conventional MRI. Therefore, we aimed to assess whether T2 and T2' quantification detects cerebral (micro-)structural damage and chronic hypoxia in lesions and in normal appearing white matter (WM) and gray matter (GM) of patients with ischemic leukoaraiosis (IL). Measurements were complemented by the assessment of the cerebral blood flow (CBF) and the degree of GM and WM atrophy.

METHODS

Eighteen patients with IL and 18 age-matched healthy controls were included. High-resolution, motion-corrected T2, T2*, and T2' mapping, CBF mapping (pulsed arterial spin labeling, PASL), and segmentation of GM and WM were used to depict specific changes in both groups. All parameters were compared between patients and healthy controls, using t testing. Values of p < 0.05 were accepted as statistically significant.

RESULTS

Patients showed significantly increased T2 in lesions (p < 0.01) and in unaffected WM (p = 0.045) as well as significantly increased T2* in lesions (p = 0.003). A significant decrease of T2' was detected in patients in unaffected WM (p = 0.027), while no T2' changes were observed in GM (p = 0.13). Both unaffected WM and GM were significantly decreased in volume in the patient-group (p < 0.01). No differences of PASL-based CBF could be shown.

CONCLUSION

Non-invasive quantitative MRI with T2, T2*, and T2' mapping might be used to detect subtle structural and metabolic changes in IL. Assessing the grade of microstructural damage and hypoxia might be helpful to monitor disease progression and to perform risk assessment.

Phase-corrected bipolar gradients in multi-echo gradient-echo sequences for quantitative susceptibility mapping

OBJECTIVE

Large echo spacing of unipolar readout gradients in current multi-echo gradient-echo (GRE) sequences for mapping fields in quantitative susceptibility mapping (QSM) can be reduced using bipolar readout gradients thereby improving acquisition efficiency.

MATERIALS AND METHODS

Phase discrepancies between odd and even echoes in the bipolar readout gradients caused by non-ideal gradient behaviors were measured, modeled as polynomials in space and corrected for accordingly in field mapping. The bipolar approach for multi-echo GRE field mapping was compared with the unipolar approach for QSM.

RESULTS

The odd-even-echo phase discrepancies were approximately constant along the phase encoding direction and linear along the readout and slice-selection directions. A simple linear phase correction in all three spatial directions was shown to enable accurate QSM of the human brain using a bipolar multi-echo GRE sequence. Bipolar multi-echo acquisition provides QSM in good quantitative agreement with unipolar acquisition while also reducing noise.

CONCLUSION

With a linear phase correction between odd-even echoes, bipolar readout gradients can be used in multi-echo GRE sequences for QSM.

Slow-MAS NMR: A new technology for in vivo metabolomic studies

Obtaining detailed in vivo metabolic information has been identified as key elements of better understanding the efficacy and toxicity of new therapies. A new nuclear magnetic resonance (NMR) technology called LOCMAT is reported in this paper that yields substantially increased spectral resolution in spatially localized in vivo H NMR metabolite spectra, as illustrated by measurements in the liver of a live mouse. LOCMAT promises to significantly enhance the utility of NMR spectroscopy for biomedical research.:

High-resolution echo-planar spectroscopic imaging of the human calf

Background

This study exploits the speed benefits of echo-planar spectroscopic imaging (EPSI) to acquire lipid spectra of skeletal muscle. The main purpose was to develop a high-resolution EPSI technique for clinical MR scanner, to visualise the bulk magnetic susceptibility (BMS) shifts of extra-myocellular lipid (EMCL) spectral lines, and to investigate the feasibility of this method for the assessment of intra-myocellular (IMCL) lipids.

Methods

The study group consisted of six healthy volunteers. A two dimensional EPSI sequence with point-resolved spectroscopy (PRESS) spatial localization was implemented on a 3T clinical MR scanner. Measurements were performed by means of 64×64 spatial matrix and nominal voxel size 3×3×15 mm³. The total net measurement time was 3 min 12 sec for non-water-suppressed (1 acquisition) and 12 min 48 sec for water-suppressed scans (4 acquisitions).

Results

Spectra of the human calf had a very good signal-to-noise ratio and linewidths sufficient to differentiate IMCL resonances from EMCL. The use of a large spatial matrix reduces inter-voxel signal contamination of the strong EMCL signals. Small voxels enabled visualisation of the methylene EMCL spectral line splitting and their BMS shifts up to 0.5 ppm relative to the correspondent IMCL line. The mean soleus muscle IMCL content of our six volunteers was 0.30±0.10 vol% (range 0.18–0.46) or 3.6±1.2 mmol/kg wet weight (range: 2.1–5.4).

Conclusion

This study demonstrates that high-spatial resolution PRESS EPSI of the muscle lipids is feasible on standard clinical scanners.

Age-related changes of cerebral autoregulation: new insights with quantitative T2'-mapping and pulsed arterial spin-labeling MR imaging

BACKGROUND AND PURPOSE

Cerebral perfusion and O(2) metabolism are affected by physiologic age-related changes. High-resolution motion-corrected quantitative T2'-imaging and PASL were used to evaluate differences in deoxygenated hemoglobin and CBF of the gray matter between young and elderly healthy subjects. Further combined T2'-imaging and PASL were investigated breathing room air and 100% O(2) to evaluate age-related changes in cerebral autoregulation.

MATERIALS AND METHODS

Twenty-two healthy volunteers 60-88 years of age were studied. Two scans of high-resolution motion-corrected T2'-imaging and PASL-MR imaging were obtained while subjects were either breathing room air or breathing 100% O(2). Manual and automated regions of interest were placed in the cerebral GM to extract values from the corresponding maps. Results were compared with those of a group of young healthy subjects previously scanned with the identical protocol as that used in the present study.

RESULTS

There was a significant decrease of cortical CBF (P < .001) and cortical T2' values (P < .001) between young and elderly healthy subjects. In both groups, T2' remained unchanged under hyperoxia compared with normoxia. Only in the younger but not in the elderly group could a significant (P = .02) hyperoxic-induced decrease of the CBF be shown.

CONCLUSIONS

T2'-mapping and PASL in the cerebral cortex of healthy subjects revealed a significant decrease of deoxygenated hemoglobin and of CBF with age. The constant deoxyHb level breathing 100% O(2) compared with normoxia in young and elderly GM suggests an age-appropriate cerebral autoregulation. At the younger age, hyperoxic-induced CBF decrease may protect the brain from hyperoxemia.

Bevacizumab impairs oxidative energy metabolism and shows antitumoral effects in recurrent glioblastomas: a 31P/1H MRSI and quantitative magnetic resonance imaging study

Bevacizumab shows unprecedented rates of response in recurrent glioblastomas (GBM), but the detailed mechanisms are still unclear. We employed in vivo magnetic resonance spectroscopic imaging (MRSI) and quantitative magnetic resonance imaging to investigate whether bevacizumab alters oxygen and energy metabolism and whether this effect has antitumoral activity in recurrent GBM. (31)P and (1)H MRSI, apparent diffusion coefficient (ADC), and high-resolution T2 and T2' mapping (indirect marker of oxygen extraction) were investigated in 16 patients with recurrent GBM at 3 Tesla before and 1.5-2 months after initiation of therapy with bevacizumab. Changes of metabolite concentrations and of the quantitative values in the tumor and normal appearing brain tissue were calculated. The Wilcoxon signed-ranks test was used to evaluate differences for tumor/edema versus control as well as changes before versus after commencement of therapy. Survival analyses were performed for significant parameters. Tumor T2', pH, ADC, and T2 decreased significantly in patients responding to bevacizumab therapy (n = 10). Patients with at least 25% T2' decrease during treatment showed longer progression-free and overall survival durations. Levels of high-energy metabolites were lower at baseline; these persisted under therapy. Glycerophosphoethanolamine as catabolic phospholipid metabolite increased in responders. The MRSI data support the hypothesis that bevacizumab induces relative tumor hypoxia (T2' decrease) and affects energy homeostasis in recurrent GBM, suggesting that bevacizumab impairs vascular function. The antiangiogenic effect of bevacizumab is predictive of better outcome and seems to induce antitumoral activity in the responding GBMs.

Quantitative T*2-mapping based on multi-slice multiple gradient echo flash imaging: retrospective correction for subject motion effects

Numerous clinical and research applications for quantitative mapping of the effective transverse relaxation time T*(2) have been described. Subject motion can severely deteriorate the quality and accuracy of results. A correction method for T*(2) maps acquired with multi-slice multiple gradient echo FLASH imaging is presented, based on acquisition repetition with reduced spatial resolution (and consequently reduced acquisition time) and weighted averaging of both data sets, choosing weighting factors individually for each k-space line to reduce the influence of motion. In detail, the procedure is based on the fact that motion artifacts reduce the correlation between acquired and exponentially fitted data. A target data set is constructed in image space, choosing the data yielding best correlation from the two acquired data sets. The k-space representation of the target is subsequently approximated as linear combination of original raw data, yielding the required weighting factors. As this method only requires a single acquisition repetition with reduced spatial resolution, it can be employed on any clinical system offering a suitable sequence with export of modulus and phase images. Experimental results show that the method works well for sparse motion, but fails for strong motion affecting the same k-space lines in both acquisitions.

A fast B1-mapping method for the correction and normalization of magnetization transfer ratio maps at 3 T

In neuroimaging, there is increasing interest in magnetization transfer (MT) techniques which yield information about bound water protons. One of the main applications is the investigation of the myelin integrity in the central nervous system (CNS). However, several problems may arise, in particular at high magnetic field strengths: B1 inhomogeneities may yield deviations of the MT saturation angle and thus non-uniformities of the measured MT ratio (MTR). This effect can be corrected for but requires in general additional time consuming B1 mapping. Furthermore, increased values of the specific absorption rate (SAR) may require a reduction of the saturation angle for individual subjects, impairing comparability of results. In this work, a B1 mapping method based on magnetization-prepared FLASH with slice selective preparation and excitation pulses and correction for relaxation effects is presented, yielding B1 maps with whole brain coverage, an in-plane resolution of 4 mm, a slice thickness of 3 mm, and a clinically acceptable duration of 46 s. The method is tested both in vitro and in vivo and applied in a subsequent in vivo study to show that MTR values in human brain tissue depend approximately linearly on the preparation angle, with a slope similar to values reported for 1.5 T. Calibration data and B1 maps are applied to B1 inhomogeneity corrections of MTR maps. Subsequently, it is shown that B1-corrected MTR maps acquired at reduced preparation angles due to individual SAR restrictions can be normalized, allowing for a direct comparison with maps acquired at the full angle.

Rapid single-scan T2*-mapping using exponential excitation pulses and image-based correction for linear background gradients

A method for fast quantitative T(2)* mapping based on multiple gradient-echo (multi-GE) imaging with correction for static magnetic field inhomogeneities is described, using an exponential excitation pulse. Field gradient maps are obtained from the phase information and modulus data are subsequently corrected, allowing for simple monoexponential T(2)* fitting. Echoes with long echo times suffering from major signal losses due to field inhomogeneities are excluded from the analysis. The acquisition time for a matrix size of 256 x 256, 1 mm in-plane resolution, and 2 mm slice thickness amounts to 15 s per slice. An additional correction for in-plane field gradients further improves accuracy. Phantom experiments show that the method provides accurate T(2)* values for field gradients up to 200 microT/m; for gradients up to 300 microT/m errors do not exceed 15%. In vivo T(2)* values acquired on healthy volunteers at 3T are in excellent agreement with results from the literature.

Noninvasive monitoring of brain temperature during mild hypothermia

The main purpose of this study was to verify the feasibility of brain temperature mapping with high-spatial- and reduced-spectral-resolution magnetic resonance spectroscopic imaging (MRSI). A secondary goal was to determine the temperature coefficient of water chemical shift in the brain with and without internal spectral reference. The accuracy of the proposed MRSI method was verified using a water and vegetable oil phantom. Selective decrease of the brain temperature of pigs was induced by intranasal cooling. Temperature reductions between 2 degrees C and 4 degrees C were achieved within 20 min. The relative changes in temperature during the cooling process were monitored using MRSI. The reference temperature was measured with MR-compatible fiber-optic probes. Single-voxel (1)H MRS was used for measurement of absolute brain temperature at baseline and at the end of cooling. The temperature coefficient of the water chemical shift of brain tissue measured by MRSI without internal reference was -0.0192+/-0.0019 ppm/degrees C. The temperature coefficients of the water chemical shift relative to N-acetylaspartate, choline-containing compounds and creatine were -0.0096+/-0.0009, -0.0083+/-0.0007 and -0.0091+/-0.0011 ppm/degrees C, respectively. The results of this study indicate that MRSI with high spatial and reduced spectral resolutions is a reliable tool for monitoring long-term temperature changes in the brain.

Reduction of susceptibility-induced signal losses in multi-gradient-echo images: application to improved visualization of the subthalamic nucleus

T2-weighted gradient echo (GE) images yield good contrast of iron-rich structures like the subthalamic nuclei due to microscopic susceptibility induced field gradients, providing landmarks for the exact placement of deep brain stimulation electrodes in Parkinson's disease treatment. An additional advantage is the low radio frequency (RF) exposure of GE sequences. However, T2-weighted images are also sensitive to macroscopic field inhomogeneities, resulting in signal losses, in particular in orbitofrontal and temporal brain areas, limiting anatomical information from these areas. In this work, an image correction method for multi-echo GE data based on evaluation of phase information for field gradient mapping is presented and tested in vivo on a 3 Tesla whole body MR scanner. In a first step, theoretical signal losses are calculated from the gradient maps and a pixelwise image intensity correction is performed. In a second step, intensity corrected images acquired at different echo times TE are combined using optimized weighting factors: in areas not affected by macroscopic field inhomogeneities, data acquired at long TE are weighted more strongly to achieve the contrast required. For large field gradients, data acquired at short TE are favored to avoid signal losses. When compared to the original data sets acquired at different TE and the respective intensity corrected data sets, the resulting combined data sets feature reduced signal losses in areas with major field gradients, while intensity profiles and a contrast-to-noise (CNR) analysis between subthalamic nucleus, red nucleus and the surrounding white matter demonstrate good contrast in deep brain areas.

Effects of simultaneous EEG recording on MRI data quality at 1.5, 3 and 7 tesla

Although the focus of attention on data degradation during simultaneous MRI/EEG recording has to date largely been upon EEG artefacts, the presence of the conducting wires and electrodes of the EEG recording system also causes some degradation of MRI data quality. This may result from magnetic susceptibility effects which lead to signal drop-out and image distortion, as well as the perturbation of the radiofrequency fields, which can cause local signal changes and a global reduction in the signal to noise ratio (SNR) of magnetic resonance images. Here, we quantify the effect of commercially available 32 and 64 electrode caps on the quality of MR images obtained in scanners operating at magnetic fields of 1.5, 3 and 7 T, via the use of MR-based, field-mapping techniques and analysis of the SNR in echo planar image time series. The electrodes are shown to be the dominant source of magnetic field inhomogeneity, although the localised nature of the field perturbation that they produce means that the effect on the signal intensity from the brain is not significant. In the particular EEG caps investigated here, RF inhomogeneity linked to the longer ECG and EOG leads causes some reduction in the signal intensity in images obtained at 3 and 7 T. Measurements of the standard deviation of white matter signal in EPI time series indicates that the introduction of the EEG cap produces a small reduction in the image signal to noise ratio, which increases with the number of electrodes used.

Fast, accurate, and precise mapping of the RF field in vivo using the 180 degrees signal null

RF field B1 nonuniformity is the largest cause of error in the quantitative measurement of many clinically relevant parameters in MR images and spectra. Knowledge of the absolute flip angles at every region will improve the accuracy and precision of such parameters. This method uses the 180 degrees signal null to construct a flip angle map of the entire brain in less than 4 min, independent of T1, T2, and proton density. Three spoiled gradient echo volume acquisitions of the whole brain were made with three different flip angles. The optimum choice of flip angles was determined to be 145 degrees, 180 degrees, and 215 degrees. Linear regression analysis was used to determine the nominal (system calibrated) flip angle required for a signal null at every pixel and thence determine the absolute flip angle at that location. The experiment utilizes an existing MR sequence supplied by the scanner manufacturer. The technique is validated experimentally and a theoretical investigation into the optimum experimental parameters is presented.

Accelerated parallel imaging for functional imaging of the human brain

Accelerated parallel imaging (PI) techniques have recently been applied to functional imaging experiments of the human brain in order to improve the performance of commonly used single-shot techniques like echo-planar imaging (EPI). Potential benefits of PI-fMRI include the reduction of geometrical distortions due to off-resonance signals, the reduction of signal-loss in areas with substantial signal inhomogeneity, increases of the spatial and temporal resolution of the fMRI experiment and reduction of gradient acoustic noise. Although PI generally leads to a substantial decrease in image signal-to-noise ratio (SNR), its effect on the temporal stability of the signal, which ultimately determines fMRI performance, is only partially determined by image SNR. Therefore, the penalty for using PI is generally not as severe as the SNR reduction. The majority of problems related to single-shot techniques become more severe at an increased magnetic field strength, making PI an important tool in achieving the full potential of fMRI at high field.

Quantification of intramyocellular lipids in obese subjects using spectroscopic imaging with high spatial resolution

Quantification of intramyocellular lipids (IMCL) in obese subjects by single-voxel spectroscopy (SVS) or conventional spectroscopic imaging (SI) often fails due to overlap of IMCL spectral lines by extramyocellular lipids (EMCL), and signal contamination from subcutaneous fat and bone marrow. This study demonstrates that these problems can be solved by high-resolution SI with 128 phase-encoding steps and a read gradient during acquisition. The small voxels obtained in this way facilitated differentiation between EMCL and IMCL. This method offers the possibility of studying different muscle groups and the variation of lipids within one muscle.

Lipid content in the musculature of the lower leg: evaluation with high-resolution spectroscopic imaging

A novel spectroscopic imaging method with high spectral and spatial resolution was developed for the specific goal of assessing muscle fat. Sensitivity to the methylene and methyl protons of fatty acids was improved by the use of a binomial 1 1 excitation pulse instead of the standard radiofrequency (RF) pulse. Acceptable measurement time is achieved by using a narrow spectral bandwidth (6 ppm). The spectral resolution is sufficient to resolve extramyocellular (EMCL) and intramyocellular (IMCL) lipids. A post-detection data processing scheme that permits correction of spectral artifacts caused by chemical shifts, spectral line aliasing, and magnetic field inhomogeneities is suggested. The lipid content in different lower leg muscles was evaluated. Muscle fiber orientation was taken into account in assessing quantities of EMCL and IMCL. The proposed technique allows small amounts of inhomogeneously distributed muscle lipids to be quantified.

High-resolution1H NMR spectroscopy in a live mouse subjected to 1.5 Hz magic angle spinning

It is demonstrated that the resolution of the (1)H NMR metabolite spectrum in a live mouse can be significantly enhanced by an ultraslow magic angle spinning of the animal combined with a modified phase-corrected magic angle turning (PHORMAT) pulse sequence. Proton NMR spectra were measured of the torso and the top part of the belly of a female BALBc mouse in a 2 T field while spinning the animal at a speed of 1.5 Hz. It was found that even in this relatively low field, with PHORMAT an isotropic spectrum is obtained with line widths that are a factor of 4.6 smaller than those obtained in a stationary mouse. It is concluded that in vivo PHORMAT has the potential to significantly increase the utility of (1)H NMR spectroscopy for biochemical and biomedical animal research.

Magnetic susceptibility quantification for arbitrarily shaped objects in inhomogeneous fields

Magnetic susceptibility measurement has wide-ranging applications in MR technical development and medical applications. A general susceptibility quantitation method for objects of arbitrary shapes in inhomogeneous magnetic fields is presented in this study. Based on the mean value properties of magnetic fields, the polarizing magnetic field at the location of interest inside an object can be exactly obtained in situ from the field values on a spherical surface enclosing the object. With numerical computation of the self-demagnetizing field and correction of contact shifts, magnetic susceptibilities were quantitatively measured for CuSO(4) phantoms based on their MR gradient echo phase maps.

A modified fuzzy clustering algorithm for operator independent brain tissue classification of dual echo MR images

Methods for brain tissue classification or segmentation of structural magnetic resonance imaging (MRI) data should ideally be independent of human operators for reasons of reliability and tractability. An algorithm is described for fully automated segmentation of dual echo, fast spin-echo MRI data. The method is used to assign fuzzy-membership values for each of four tissue classes (gray matter, white matter, cerebrospinal fluid and dura) to each voxel based on partition of a two dimensional feature space. Fuzzy clustering is modified for this application in two ways. First, a two component normal mixture model is initially fitted to the thresholded feature space to identify exemplary gray and white matter voxels. These exemplary data protect subsequently estimated cluster means against the tendency of unmodified fuzzy clustering to equalize the number of voxels in each class. Second, fuzzy clustering is implemented in a moving window scheme that accommodates reduced image contrast at the axial extremes of the transmitting/receiving coil. MRI data acquired from 5 normal volunteers were used to identify stable values for three arbitrary parameters of the algorithm: feature space threshold, relative weight of exemplary gray and white matter voxels, and moving window size. The modified algorithm incorporating these parameter values was then used to classify data from simulated images of the brain, validating the use of fuzzy-membership values as estimates of partial volume. Gray:white matter ratios were estimated from 20 twenty normal volunteers (mean age 32.8 years). Processing time for each three-dimensional image was approximately 30 min on a 170 MHz workstation. Mean cerebral gray and white matter volumes estimated from these automatically segmented images were very similar to comparable results previously obtained by operator dependent methods, but without their inherent unreliability.

Chemical shift artifact-free microscopy: spectroscopic microimaging of the human skin

A spectroscopic imaging technique with high spatial resolution was used for the study of human skin in vivo. The measurements were performed using a whole-body magnetic resonance system (1.5 T) with standard gradients and a standard 8-cm diameter circular surface coil. A decisive gain in signal-to-noise ratio was achieved by reducing the receiver bandwidth of the imaging system to values less than +/-5 kHz. The chemical shift misregistration was eliminated by post-detection data processing. The method was tested on different kinds of skin, on the foot sole and head. Water, fat, and chemical shift artifact-free images were obtained with resolution 0.107 x 0.143 mm in plane and slice thickness 1 mm. A major advantage of the spectroscopic imaging procedure is that the pulse sequence can be optimized for the maximum signal-to-noise ratio. There is no need for special modification of the sequence to circumvent the chemical shift artifacts (water, fat suppression, etc.).

Magnetic resonance spectroscopic imaging for visualization and correction of distortions in MRI: high precision applications in neurosurgery

We present a method for the quantification and correction of geometrical/intensity distortions of magnetic resonance images predominantly caused by bulk magnetic susceptibility shifts due to susceptibility heterogeneities of measured biologic tissues and shape of the object under investigation. The method includes precise and fast measurements of the static magnetic-field distribution inside the measured object and automated data processing. Magnetic-field deviations in the range (-2.4; 2.6) ppm were found in the human brain at B0 = 1.5 T. For routinely used imaging parameters, with a read gradient strength of about approximately 1 mT/m, the magnetic-field perturbations in the human brain can cause geometrical distortions up to +/-4 mm and intensity changes up to +/-50%. MR images corrected by the described method are suitable for planning high precision applications in neurosurgery.

Design of practical T2-selective RF excitation (TELEX) pulses

Traditional T2-based imaging techniques are geared toward imaging long-T2 species. Traditional techniques are, therefore, not optimal in clinical situations where the information of interest lies in the short-T2 species. T2-selective RF excitation (TELEX) is a technique for obtaining a T2-based contrast that highlights short-T2 values while suppressing long-T2 values-opposite to traditional T2 contrast. Previously, TELEX has been demonstrated qualitatively to highlight only very short-T2 values (T2 approximately 0.001 s). When applied to longer T2 values (T2 > or = 0.01 s), TELEX becomes sensitive to deltaB0 non-uniformities. This restricts its application to problems in which the T2 of interest is very short. In this study, TELEX is characterized quantitatively. Furthermore, a bandwidth broadening scheme is developed that reduces the deltaB0 sensitivity of TELEX. This permits the technique to be applied to longer T2 values. The capabilities and limitations of a practical implementation of TELEX are discussed.

Chemical shift artifact-free imaging: a new option in MRI?

A high-speed proton spectroscopic imaging method with high spatial resolution was used for obtaining water, fat, and chemical shift artifact-free images on a 1.5 T MR scanner. The technique is based on a fast radiofrequency (RF) spoiled gradient-echo sequence. The chemical shift information is encoded by incrementing the echo time in a series of image records. Suppression of water or fat signals is not used. The technique does not require a highly homogeneous magnetic field. Spectroscopic images of a human volunteer were compared with corresponding conventional images obtained using the short inversion time inversion recovery (STIR) and the selective partial inversion recovery (SPIR) methods. The results demonstrate that it is possible to produce images entirely free from chemical shift artifacts using only a few chemical shift encoding steps. The technique also produces pure water and fat images which are significantly better than those produced by using the conventional methods STIR and selective partial inversion recovery. The described method appears to be promising for routine clinical applications because it can be fully automated.

Multiple-echo proton spectroscopic imaging using time domain parametric spectral analysis

A multiple-echo MR spectroscopic imaging (MRSI) method is presented that enables improved metabolite imaging in the presence of local field inhomogeneities and measurement of transverse relaxation parameters. Short echo spacing is used to maximize signal energy from inhomogeneously line-broadened resonances, and time domain parametric spectral analysis of the entire echo train is used to obtain sufficient spectral resolution from the shortened sampling periods. Optimal sequence parameters for 1H MRSI are determined by computer simulation, and performance is compared with conventional single-echo acquisition using phantom studies at a field strength of 4.7 T. A preliminary example for use at 1.5 T is also presented using phantom and human brain MRSI studies. This technique is shown to offer improved performance relative to single-echo MRSI for imaging of metabolites with shortened T2* values due to the presence of local field inhomogeneities. Additional advantages are the intrinsic measurement of metabolite T2 values and determination of metabolite integrals without T2 weighting, thereby facilitating quantitative metabolite imaging.

1H-spectroscopic imaging with read gradient during acquisition in inhomogeneous fields: analysis, measurement strategy, and data processing

The proton magnetic resonance spectroscopic imaging techniques that use read gradient during acquisition produce proton spectra with high spatial and moderately high spectroscopic resolution in a reasonable time for in vivo applications. These techniques suffer mainly from the spatial and spectral distortions caused by the convolution of spectral/spatial information (chemical-shift artifacts) and from the spectral shifts caused by static magnetic field inhomogeneities. The investigators analyze the chemical-shift artifacts in the presence of nonnegligible static magnetic field inhomogeneities and propose a postdetection processing scheme to correct for such effects. Spectral artifacts caused by chemical shifts, spectral line overlapping, streak broadening, and magnetic field inhomogeneities are discussed. The postdetection data processing scheme is demonstrated on measurements of a phantom as well as a human leg.

An in vivo automated shimming method taking into account shim current constraints

Many in vivo imaging techniques require magnetic field homogeneity in the volume of interest. Shim coils of the second and third order spherical harmonics have been used successfully to compensate for complicated field variations caused by the human anatomy itself. The available currents of these coils are invariably limited. In this note we demonstrate that these limits significantly affect the optimal shim condition. We propose an automated in vivo shimming method for arbitrary volumes of interest using 3-dimensional (3D) field maps. This method is a modification of previous works using least-squares criteria. The main difference is that a constrained optimization is performed in vivo under the current limits of the shim coils, which improved the field homogeneity significantly over simple truncations of the least-squares solutions. This shimming method was used with head scans of five normal volunteers on a 4.0 tesla scanner. A fast double-echo sequence was used to obtain field maps, and a new field uniformity measure was derived for this method. The field mapping sequence was tested against a standard single-echo Dixon sequence used by previous investigators, and the stability of the shimming method was tested by repeated studies on the same subject.