Clinical Evaluation of Scout Accelerated Motion Estimation and Reduction Technique for 3D MR Imaging in the Inpatient and Emergency Department Settings.

BACKGROUND AND PURPOSE: A scout accelerated motion estimation and reduction (SAMER) framework has been developed for efficient retrospective motion correction. The goal of this study was to perform an initial evaluation of SAMER in a series of clinical brain MR imaging examinations. MATERIALS AND METHODS: Ninety-seven patients who underwent MR imaging in the inpatient and emergency department settings were included in the study. SAMER motion correction was retrospectively applied to an accelerated T1-weighted MPRAGE sequence that was included in brain MR imaging examinations performed with and without contrast. Two blinded neuroradiologists graded images with and without SAMER motion correction on a 5-tier motion severity scale (none = 1, minimal = 2, mild = 3, moderate = 4, severe = 5). RESULTS: The median SAMER reconstruction time was 1 minute 47 seconds. SAMER motion correction significantly improved overall motion grades across all examinations (P < .005). Motion artifacts were reduced in 28% of cases, unchanged in 64% of cases, and increased in 8% of cases. SAMER improved motion grades in 100% of moderate motion cases and 75% of severe motion cases. Sixty-nine percent of nondiagnostic motion cases (grades 4 and 5) were considered diagnostic after SAMER motion correction. For cases with minimal or no motion, SAMER had negligible impact on the overall motion grade. For cases with mild, moderate, and severe motion, SAMER improved the motion grade by an average of 0.3 (SD, 0.5), 1.1 (SD, 0.3), and 1.1 (SD, 0.8) grades, respectively. CONCLUSIONS: SAMER improved the diagnostic image quality of clinical brain MR imaging examinations with motion artifacts. The improvement was most pronounced for cases with moderate or severe motion.

Imaging G-Ratio in Multiple Sclerosis Using High-Gradient Diffusion MRI and Macromolecular Tissue Volume.

BACKGROUND AND PURPOSE: Remyelination represents an area of great therapeutic interest in multiple sclerosis but currently lacks a robust imaging marker. The purpose of this study was to use high-gradient diffusion MRI and macromolecular tissue volume imaging to obtain estimates of axonal volume fraction, myelin volume fraction, and the imaging g-ratio in patients with MS and healthy controls and to explore their relationship to neurologic disability in MS. MATERIALS AND METHODS: Thirty individuals with MS (23 relapsing-remitting MS, 7 progressive MS) and 19 age-matched healthy controls were scanned on a 3T MRI scanner equipped with 300 mT/m maximum gradient strength using a comprehensive multishell diffusion MRI protocol. Macromolecular tissue volume imaging was performed to quantify the myelin volume fraction. Diffusion data were fitted to a 3-compartment model of white matter using a spheric mean approach to yield estimates of axonal volume fraction. The imaging g-ratio was calculated from the ratio of myelin volume fraction and axonal volume fraction. Imaging metrics were compared between groups using 2-sided t tests with a Bonferroni correction. RESULTS: The mean g-ratio was significantly elevated in lesions compared with normal-appearing WM (0.74 vs 0.67, P < .001). Axonal volume fraction (0.17 vs 0.23, P < .001) and myelin volume fraction (0.17 vs 0.25, P < .001) were significantly lower in lesions than normal-appearing WM. Myelin volume fraction was lower in normal-appearing WM compared with that in healthy controls (0.25 vs 0.27, P = .009). Disability, as measured by the Expanded Disability Status Scale, was significantly associated with myelin volume fraction (ß = -40.5, P = .001) and axonal volume fraction (ß = -41.0, P = .016) in normal-appearing WM. CONCLUSIONS: The imaging g-ratio may serve as a biomarker for the relative degree of axonal and myelin loss in MS.

g-Ratio weighted imaging of the human spinal cord in vivo.

The fiber g-ratio is defined as the ratio of the inner to the outer diameter of the myelin sheath. This ratio provides a measure of the myelin thickness that complements axon morphology (diameter and density) for assessment of demyelination in diseases such as multiple sclerosis. Previous work has shown that an aggregate g-ratio map can be computed using a formula that combines axon and myelin density measured with quantitative MRI. In this work, we computed g-ratio weighted maps in the cervical spinal cord of nine healthy subjects. We utilized the 300mT/m gradients from the CONNECTOM scanner to estimate the fraction of restricted water (fr) with high accuracy, using the CHARMED model. Myelin density was estimated using the lipid and macromolecular tissue volume (MTV) method, derived from normalized proton density (PD) mapping. The variability across spinal level, laterality and subject were assessed using a three-way ANOVA. The average g-ratio value obtained in the white matter was 0.76+/-0.03, consistent with previous histology work. Coefficients of variation of fr and MTV were respectively 4.3% and 13.7%. fr and myelin density were significantly different across spinal tracts (p=3×10-7 and 0.004 respectively) and were positively correlated in the white matter (r=0.42), suggesting shared microstructural information. The aggregate g-ratio did not show significant differences across tracts (p=0.6). This study suggests that fr and myelin density can be measured in vivo with high precision and that they can be combined to produce a g-ratio-weighted map robust to free water pool contamination from cerebrospinal fluid or veins. Potential applications include the study of early demyelination in multiple sclerosis, and the quantitative assessment of remyelination drugs.

Pushing the limits of in vivo diffusion MRI for the Human Connectome Project.

Perhaps more than any other "-omics" endeavor, the accuracy and level of detail obtained from mapping the major connection pathways in the living human brain with diffusion MRI depend on the capabilities of the imaging technology used. The current tools are remarkable; allowing the formation of an "image" of the water diffusion probability distribution in regions of complex crossing fibers at each of half a million voxels in the brain. Nonetheless our ability to map the connection pathways is limited by the image sensitivity and resolution, and also the contrast and resolution in encoding of the diffusion probability distribution. The goal of our Human Connectome Project (HCP) is to address these limiting factors by re-engineering the scanner from the ground up to optimize the high b-value, high angular resolution diffusion imaging needed for sensitive and accurate mapping of the brain's structural connections. Our efforts were directed based on the relative contributions of each scanner component. The gradient subsection was a major focus since gradient amplitude is central to determining the diffusion contrast, the amount of T2 signal loss, and the blurring of the water PDF over the course of the diffusion time. By implementing a novel 4-port drive geometry and optimizing size and linearity for the brain, we demonstrate a whole-body sized scanner with G(max) = 300 mT/m on each axis capable of the sustained duty cycle needed for diffusion imaging. The system is capable of slewing the gradient at a rate of 200 T/m/s as needed for the EPI image encoding. In order to enhance the efficiency of the diffusion sequence we implemented a FOV shifting approach to Simultaneous MultiSlice (SMS) EPI capable of unaliasing 3 slices excited simultaneously with a modest g-factor penalty allowing us to diffusion encode whole brain volumes with low TR and TE. Finally we combine the multi-slice approach with a compressive sampling reconstruction to sufficiently undersample q-space to achieve a DSI scan in less than 5 min. To augment this accelerated imaging approach we developed a 64-channel, tight-fitting brain array coil and show its performance benefit compared to a commercial 32-channel coil at all locations in the brain for these accelerated acquisitions. The technical challenges of developing the over-all system are discussed as well as results from SNR comparisons, ODF metrics and fiber tracking comparisons. The ultra-high gradients yielded substantial and immediate gains in the sensitivity through reduction of TE and improved signal detection and increased efficiency of the DSI or HARDI acquisition, accuracy and resolution of diffusion tractography, as defined by identification of known structure and fiber crossing.

Effects of image reconstruction on fiber orientation mapping from multichannel diffusion MRI: reducing the noise floor using SENSE.

PURPOSE: To examine the effects of the reconstruction algorithm of magnitude images from multichannel diffusion MRI on fiber orientation estimation. THEORY AND METHODS: It is well established that the method used to combine signals from different coil elements in multichannel MRI can have an impact on the properties of the reconstructed magnitude image. Using a root-sum-of-squares approach results in a magnitude signal that follows an effective noncentral-χ distribution. As a result, the noise floor, the minimum measurable in the absence of any true signal, is elevated. This is particularly relevant for diffusion-weighted MRI, where the signal attenuation is of interest. RESULTS: In this study, we illustrate problems that such image reconstruction characteristics may cause in the estimation of fiber orientations, both for model-based and model-free approaches, when modern 32-channel coils are used. We further propose an alternative image reconstruction method that is based on sensitivity encoding (SENSE) and preserves the Rician nature of the single-channel, magnitude MR signal. We show that for the same k-space data, root-sum-of-squares can cause excessive overfitting and reduced precision in orientation estimation compared with the SENSE-based approach. CONCLUSION: These results highlight the importance of choosing the appropriate image reconstruction method for tractography studies that use multichannel receiver coils for diffusion MRI acquisition.

Improving diffusion MRI using simultaneous multi-slice echo planar imaging.

In diffusion MRI, simultaneous multi-slice single-shot EPI acquisitions have the potential to increase the number of diffusion directions obtained per unit time, allowing more diffusion encoding in high angular resolution diffusion imaging (HARDI) acquisitions. Nonetheless, unaliasing simultaneously acquired, closely spaced slices with parallel imaging methods can be difficult, leading to high g-factor penalties (i.e., lower SNR). The CAIPIRINHA technique was developed to reduce the g-factor in simultaneous multi-slice acquisitions by introducing inter-slice image shifts and thus increase the distance between aliased voxels. Because the CAIPIRINHA technique achieved this by controlling the phase of the RF excitations for each line of k-space, it is not directly applicable to single-shot EPI employed in conventional diffusion imaging. We adopt a recent gradient encoding method, which we termed "blipped-CAIPI", to create the image shifts needed to apply CAIPIRINHA to EPI. Here, we use pseudo-multiple replica SNR and bootstrapping metrics to assess the performance of the blipped-CAIPI method in 3× simultaneous multi-slice diffusion studies. Further, we introduce a novel image reconstruction method to reduce detrimental ghosting artifacts in these acquisitions. We show that data acquisition times for Q-ball and diffusion spectrum imaging (DSI) can be reduced 3-fold with a minor loss in SNR and with similar diffusion results compared to conventional acquisitions.

An implanted 8-channel array coil for high-resolution macaque MRI at 3T.

An 8-channel receive coil array was constructed and implanted adjacent to the skull in a male rhesus monkey in order to improve the sensitivity of (functional) brain imaging. The permanent implant was part of an acrylic headpost assembly and only the coil element loop wires were implanted. The tuning, matching, and preamplifier circuitry was connected via a removable external assembly. Signal-to-noise ratio (SNR) and noise amplification for parallel imaging were compared to single-, 4-, and 8-channel external receive-only coils routinely used for macaque fMRI. In vivo measurements showed significantly improved SNR within the brain for the implanted versus the external coils. Within a region-of-interest covering the cerebral cortex, we observed a 5.4-, 3.6-fold, and 3.4-fold increase in SNR compared to the external single-, 4-, and 8-channel coils, respectively. In the center of the brain, the implanted array maintained a 2.4×, 2.5×, and 2.1× higher SNR, respectively compared to the external coils. The array performance was evaluated for anatomical, diffusion tensor and functional brain imaging. This study suggests that a stable implanted phased-array coil can be used in macaque MRI to substantially increase the spatial resolution for anatomical, diffusion tensor, and functional imaging.

T2* mapping and B0 orientation-dependence at 7 T reveal cyto- and myeloarchitecture organization of the human cortex.

Ultra-high field MRI (≥ 7 T) has recently shown great sensitivity to depict patterns of tissue microarchitecture. Moreover, recent studies have demonstrated a dependency between T2* and orientation of white matter fibers with respect to the main magnetic field B0. In this study we probed the potential of T2* mapping at 7 T to provide new markers of cortical architecture. We acquired multi-echo measurements at 7 T and mapped T2* over the entire cortex of eight healthy individuals using surface-based analysis. B0 dependence was tested by computing the angle θ(z) between the normal of the surface and the direction of B0, then fitting T2*(θ(z)) using model from the literature. Average T2* in the cortex was 32.20 +/- 1.35 ms. Patterns of lower T2* were detected in the sensorimotor, visual and auditory cortices, likely reflecting higher myelin content. Significantly lower T2* was detected in the left hemisphere of the auditory region (p<0.005), suggesting higher myelin content, in accordance with previous investigations. B0 orientation dependence was detected in some areas of the cortex, the strongest being in the primary motor cortex (∆R2*=4.10 Hz). This study demonstrates that quantitative T2* measures at 7 T MRI can reveal patterns of cytoarchitectural organization of the human cortex in vivo and that B0 orientation dependence can probe the coherency and orientation of gray matter fibers in the cortex, shedding light into the potential use of this type of contrast to characterize cyto-/myeloarchitecture and to understand the pathophysiology of diseases associated with changes in iron and/or myelin concentration.

32-channel RF coil optimized for brain and cervical spinal cord at 3 T.

Diffusion and functional magnetic resonance imaging of the spinal cord remain challenging due to the small cross-sectional size of the cord and susceptibility-related distortions. Although partially addressable through parallel imaging, few highly parallel array coils have been implemented for the cervical cord. Here, we developed a 32-channel coil that fully covers the brain and c-spine and characterized its performance in comparison with a commercially available head/neck/spine array. Image and temporal signal-to-noise ratio were, respectively, increased by 2× and 1.8× in the cervical cord. Averaged g-factors at 4× acceleration were lowered by 22% in the brain and by 39% in the spinal cord, enabling 1-mm isotropic R = 4 multi-echo magnetization prepared gradient echo of the full brain and c-spine in 3:20 min. Diffusion imaging of the cord at 0.6 × 0.6 × 5 mm(3) resolution and tractography of the full brain and c-spine at 1.7-mm isotropic resolution were feasible without noticeable distortion. Improvements of this nature potentially enhance numerous basic and clinical research studies focused on spinal and supraspinal regions.

Simulating Magnetic Nanoparticle Behavior in Low-field MRI under Transverse Rotating Fields and Imposed Fluid Flow.

In the presence of alternating-sinusoidal or rotating magnetic fields, magnetic nanoparticles will act to realign their magnetic moment with the applied magnetic field. The realignment is characterized by the nanoparticle's time constant, τ. As the magnetic field frequency is increased, the nanoparticle's magnetic moment lags the applied magnetic field at a constant angle for a given frequency, Ω, in rad/s. Associated with this misalignment is a power dissipation that increases the bulk magnetic fluid's temperature which has been utilized as a method of magnetic nanoparticle hyperthermia, particularly suited for cancer in low-perfusion tissue (e.g., breast) where temperature increases of between 4°C and 7°C above the ambient in vivo temperature cause tumor hyperthermia. This work examines the rise in the magnetic fluid's temperature in the MRI environment which is characterized by a large DC field, B(0). Theoretical analysis and simulation is used to predict the effect of both alternating-sinusoidal and rotating magnetic fields transverse to B(0). Results are presented for the expected temperature increase in small tumors (~1 cm radius) over an appropriate range of magnetic fluid concentrations (0.002 to 0.01 solid volume fraction) and nanoparticle radii (1 to 10 nm). The results indicate that significant heating can take place, even in low-field MRI systems where magnetic fluid saturation is not significant, with careful The goal of this work is to examine, by means of analysis and simulation, the concept of interactive fluid magnetization using the dynamic behavior of superparamagnetic iron oxide nanoparticle suspensions in the MRI environment. In addition to the usual magnetic fields associated with MRI, a rotating magnetic field is applied transverse to the main B(0) field of the MRI. Additional or modified magnetic fields have been previously proposed for hyperthermia and targeted drug delivery within MRI. Analytical predictions and numerical simulations of the transverse rotating magnetic field in the presence of B(0) are investigated to demonstrate the effect of Ω, the rotating field frequency, and the magnetic field amplitude on the fluid suspension magnetization. The transverse magnetization due to the rotating transverse field shows strong dependence on the characteristic time constant of the fluid suspension, τ. The analysis shows that as the rotating field frequency increases so that Ωτ approaches unity, the transverse fluid magnetization vector is significantly non-aligned with the applied rotating field and the magnetization's magnitude is a strong function of the field frequency. In this frequency range, the fluid's transverse magnetization is controlled by the applied field which is determined by the operator. The phenomenon, which is due to the physical rotation of the magnetic nanoparticles in the suspension, is demonstrated analytically when the nanoparticles are present in high concentrations (1 to 3% solid volume fractions) more typical of hyperthermia rather than in clinical imaging applications, and in low MRI field strengths (such as open MRI systems), where the magnetic nanoparticles are not magnetically saturated. The effect of imposed Poiseuille flow in a planar channel geometry and changing nanoparticle concentration is examined. The work represents the first known attempt to analyze the dynamic behavior of magnetic nanoparticles in the MRI environment including the effects of the magnetic nanoparticle spin-velocity. It is shown that the magnitude of the transverse magnetization is a strong function of the rotating transverse field frequency. Interactive fluid magnetization effects are predicted due to non-uniform fluid magnetization in planar Poiseuille flow with high nanoparticle concentrations.

Three dimensional echo-planar imaging at 7 Tesla.

Functional MRI (fMRI) most commonly employs 2D echo-planar imaging (EPI). The advantages for fMRI brought about by the increasingly popular ultra-high field strengths are best exploited in high-resolution acquisitions, but here 2D EPI becomes impractical for several reasons, including the very long volume acquisitions times. In this study at 7 T, a 3D EPI sequence with full parallel and partial Fourier imaging capability along both phase encoding axes was implemented and used to evaluate the sensitivity of 3D and corresponding 2D EPI acquisitions at four different spatial resolutions ranging from small to typical voxel sizes (1.5-3.0 mm isotropic). Whole-brain resting state measurements (N=4) revealed a better, or at least comparable sensitivity of the 3D method for gray and white matter. The larger vulnerability of 3D to physiological effects was outweighed by the much shorter volume TR, which moreover allows whole-brain coverage at high resolution within fully acceptable limits for event-related fMRI: TR was only 3.07 s for 1.5 mm, 1.88 s for 2.0 mm, 1.38 s for 2.5 mm and 1.07 s for 3.0 mm isotropic resolution. In order to investigate the ability to detect and spatially resolve BOLD activation in the visual cortex, functional 3D EPI experiments (N=8) were performed at 1 mm isotropic resolution with parallel imaging acceleration of 3x3, resulting in a TR of only 3.2 s for whole-brain coverage. From our results, and several other practical advantages of 3D over 2D EPI found in the present study, we conclude that 3D EPI provides a useful alternative for whole-brain fMRI at 7 T, not only when high-resolution data are required.

Heating in the MRI environment due to superparamagnetic fluid suspensions in a rotating magnetic field.

In the presence of alternating-sinusoidal or rotating magnetic fields, magnetic nanoparticles will act to realign their magnetic moment with the applied magnetic field. The realignment is characterized by the nanoparticle's time constant, τ. As the magnetic field frequency is increased, the nanoparticle's magnetic moment lags the applied magnetic field at a constant angle for a given frequency, Ω, in rad/s. Associated with this misalignment is a power dissipation that increases the bulk magnetic fluid's temperature which has been utilized as a method of magnetic nanoparticle hyperthermia, particularly suited for cancer in low-perfusion tissue (e.g., breast) where temperature increases of between 4°C and 7°C above the ambient in vivo temperature cause tumor hyperthermia. This work examines the rise in the magnetic fluid's temperature in the MRI environment which is characterized by a large DC field, B(0). Theoretical analysis and simulation is used to predict the effect of both alternating-sinusoidal and rotating magnetic fields transverse to B(0). Results are presented for the expected temperature increase in small tumors (~1 cm radius) over an appropriate range of magnetic fluid concentrations (0.002 to 0.01 solid volume fraction) and nanoparticle radii (1 to 10 nm). The results indicate that significant heating can take place, even in low-field MRI systems where magnetic fluid saturation is not significant, with careful selection of the rotating or sinusoidal field parameters (field frequency and amplitude). The work indicates that it may be feasible to combine low-field MRI with a magnetic hyperthermia system using superparamagnetic iron oxide nanoparticles.

Broadband slab selection with B1+ mitigation at 7T via parallel spectral-spatial excitation.

Chemical shift imaging benefits from signal-to-noise ratio (SNR) and chemical shift dispersion increases at stronger main field such as 7 Tesla, but the associated shorter radiofrequency (RF) wavelengths encountered require B1+ mitigation over both the spatial field of view (FOV) and a specified spectral bandwidth. The bandwidth constraint presents a challenge for previously proposed spatially tailored B1+ mitigation methods, which are based on a type of echovolumnar trajectory referred to as "spokes" or "fast-kz". Although such pulses, in conjunction with parallel excitation methodology, can efficiently mitigate large B1+ inhomogeneities and achieve relatively short pulse durations with slice-selective excitations, they exhibit a narrow-band off-resonance response and may not be suitable for applications that require B1+ mitigation over a large spectral bandwidth. This work outlines a design method for a general parallel spectral-spatial excitation that achieves a target-error minimization simultaneously over a bandwidth of frequencies and a specified spatial-domain. The technique is demonstrated for slab-selective excitation with in-plane B1+ mitigation over a 600-Hz bandwidth. The pulse design method is validated in a water phantom at 7T using an eight-channel transmit array system. The results show significant increases in the pulse's spectral bandwidth, with no additional pulse duration penalty and only a minor tradeoff in spatial B1+ mitigation compared to the standard spoke-based parallel RF design.

Magnitude least squares optimization for parallel radio frequency excitation design demonstrated at 7 Tesla with eight channels.

Spatially tailored radio frequency (RF) excitations accelerated with parallel transmit systems provide the opportunity to create shaped volume excitations or mitigate inhomogeneous B(1) excitation profiles with clinically relevant pulse lengths. While such excitations are often designed as a least-squares optimized approximation to a target magnitude and phase profile, adherence to the target phase profile is usually not important as long as the excitation phase is slowly varying compared with the voxel dimension. In this work, we demonstrate a method for a magnitude least squares optimization of the target magnetization profile for multichannel parallel excitation to improve the magnitude profile and reduce the RF power at the cost of a less uniform phase profile. The method enables the designer to trade off the allowed spatial phase variation for the improvement in magnitude profile and reduction in RF power. We validate the method with simulation studies and demonstrate its performance in fourfold accelerated two-dimensional spiral excitations, as well as for uniform in-plane slice selective parallel excitations using an eight-channel transmit array on a 7T human MRI scanner. The experimental results are in good agreement with the simulations, which show significant improvement in the magnitude profile and reductions in the required RF power while still maintaining negligible intravoxel phase variation.

32-channel 3 Tesla receive-only phased-array head coil with soccer-ball element geometry.

A 32-channel 3T receive-only phased-array head coil was developed for human brain imaging. The helmet-shaped array was designed to closely fit the head with individual overlapping circular elements arranged in patterns of hexagonal and pentagonal symmetry similar to that of a soccer ball. The signal-to-noise ratio (SNR) and noise amplification (g-factor) in accelerated imaging applications were quantitatively evaluated in phantom and human images and compared with commercially available head coils. The 32-channel coil showed SNR gains of up to 3.5-fold in the cortex and 1.4-fold in the corpus callosum compared to a (larger) commercial eight-channel head coil. The experimentally measured g-factor performance of the helmet array showed significant improvement compared to the eight-channel array (peak g-factor 59% and 26% of the eight-channel values for four- and fivefold acceleration). The performance of the arrays is demonstrated in high-resolution and highly accelerated brain images.

3T phased array MRI improves the presurgical evaluation in focal epilepsies: a prospective study.

BACKGROUND: Although detection of concordant lesions on MRI significantly improves postsurgical outcomes in focal epilepsy (FE), many conventional MR studies remain negative. The authors evaluated the role of phased array surface coil studies performed at 3 Tesla (3T PA MRI). METHODS: Forty patients with medically intractable focal epilepsies were prospectively imaged with 3T PA-MRI including high matrix TSE T2, fluid attenuated inversion recovery, and magnetization prepared rapid gradient echo. All patients were considered candidates for epilepsy surgery. 3T PA-MRIs were reviewed by a neuroradiologist experienced in epilepsy imaging with access to clinical information. Findings were compared to reports of prior standard 1.5T MRI epilepsy studies performed at tertiary care centers. RESULTS: Experienced, unblinded review of 3T PA-MRI studies yielded additional diagnostic information in 48% (19/40) compared to routine clinical reads at 1.5T. In 37.5% (15/40), this additional information motivated a change in clinical management. In the subgroup of patients with prior 1.5T MRIs interpreted as normal, 3T PA-MRI resulted in the detection of a new lesion in 65% (15/23). In the subgroup of 15 patients with known lesions, 3T PA-MRI better defined the lesion in 33% (5/15). CONCLUSION: Phased array surface coil studies performed at 3 Tesla read by an experienced unblinded neuroradiologist can improve the presurgical evaluation of patients with focal epilepsy when compared to routine clinical 1.5T studies read at tertiary care centers.

Nonstationary noise estimation in functional MRI.

An important issue in functional MRI analysis is accurate characterisation of the noise processes present in the data. Whilst conventional fMRI noise representations often assume stationarity (or time-invariance) in the noise generating sources, such approaches may serve to suppress important dynamic information about brain function. As an alternative to these fixed temporal assumptions, we present in this paper two time-varying procedures for examining nonstationary noise structure in fMRI data. In the first procedure, we approximate nonstationary behaviour by means of a collection of simple but numerous time-varying parametric models. This is accomplished through the derivation of a locally parametric AutoRegressive (AR) plus drift model which tracks temporal covariance by allowing the model parameters to evolve over time. Before exploring time variation in these parameters, window-widths (bandwidths) that are well suited to the latent time-varying noise structure must be determined. To do this, we employ a bandwidth selection mechanism based on Stein's Unbiased Risk Estimator (SURE) criterion. In the second procedure, we describe the fMRI noise using a nonparametric method based on Functional Data Analysis (FDA). This process generates well-conditioned nonstationary covariance estimates that reflect temporal continuity in the underlying data structure whilst penalizing effective model dimension. We demonstrate both methods on simulated data and investigate the presence of nonstationary noise in resting fMRI data using the whitening capabilities of the locally parametric procedure. We evaluate the comparative behaviour of the stationary and nonstationary AR-based methods on data acquired at 1.5, 3 and 7 T magnetic field strengths and show that incorporation of time variation in the AR parameters leads to an overall decrease in the level of residual structure in the data. The FDA noise modelling technique is formulated within an activation mapping procedure and compared to the SPM (Statistical Parametric Mapping) toolbox on a cognitive face recognition task. Both the SPM and FDA methods show good sensitivity on this task, but we find that inclusion of the nonstationary FDA noise model seems to improve detection power in important task-related medial temporal regions.

Eight-channel phased array coil and detunable TEM volume coil for 7 T brain imaging.

An eight-channel receive-only brain coil and table-top detunable volume transmit coil were developed and tested at 7 T for human imaging. Optimization of this device required attention to sources of interaction between the array elements, between the transmit and receive coils and minimization of common mode currents on the coaxial cables. Circular receive coils (85 mm dia.) were designed on a flexible former to fit tightly around the head and within a 270-mm diameter TEM transmit volume coil. In the near cortex, the array provided a fivefold increase in SNR compared to a TEM transmit-receive coil, a gain larger than that seen in comparable coils at 3 T. The higher SNR gain is likely due to strong dielectric effects, which cause the volume coil to perform poorly in the cortex compared to centrally. The sensitivity and coverage of the array is demonstrated with high-resolution images of the brain cortex.

Comparison of physiological noise at 1.5 T, 3 T and 7 T and optimization of fMRI acquisition parameters.

Previous studies have shown that under some conditions, noise fluctuations in an fMRI time-course are dominated by physiological modulations of the image intensity with secondary contributions from thermal image noise and that these two sources scale differently with signal intensity, susceptibility weighting (TE) and field strength. The SNR of the fMRI time-course was found to be near its asymptotic limit for moderate spatial resolution measurements at 3 T with only marginal gains expected from acquisition at higher field strengths. In this study, we investigate the amplitude of image intensity fluctuations in the fMRI time-course at magnetic field strengths of 1.5 T, 3 T, and 7 T as a function of image resolution, flip angle and TE. The time-course SNR was a similar function of the image SNR regardless of whether the image SNR was modulated by flip angle, image resolution, or field strength. For spatial resolutions typical of those currently used in fMRI (e.g., 3 x 3 x 3 mm(3)), increases in image SNR obtained from 7 T acquisition produced only modest increases in time-course SNR. At this spatial resolution, the ratio of physiological noise to thermal image noise was 0.61, 0.89, and 2.23 for 1.5 T, 3 T, and 7 T. At a resolution of 1 x 1 x 3 mm(3), however, the physiological to thermal noise ratio was 0.34, 0.57, and 0.91 for 1.5 T, 3 T and 7 T for TE near T2*. Thus, by reducing the signal strength using higher image resolution, the ratio of physiologic to image noise could be reduced to a regime where increased sensitivity afforded by higher field strength still translated to improved SNR in the fMRI time-series.

In vivo GABA+ measurement at 1.5T using a PRESS-localized double quantum filter.

A point-resolved spectroscopy (PRESS)-localized double quantum filter was implemented on a 1.5T clinical scanner for the estimation of gamma-amino butyric acid (GABA) concentrations in vivo. Several calibrations were found to be necessary for consistent results to be obtained. The apparent filter yield was approximately 38%; filter strength was sufficient to reduce the singlet metabolite peaks in vivo to below the level of the noise. Metabolite-nulled experiments were performed, which confirmed that significant overlap occurred between macromolecule signals and the GABA resonance at 3.1 ppm. Although the multiplet arm at 2.9 ppm was confirmed to be relatively free of contamination with macromolecules, some contribution from these and from peptides is likely to remain; therefore, the term GABA+ is used. GABA+ concentrations were estimated relative to creatine (Cr) at the same echo time (TE) in a group of controls, studied on two occasions. The GABA+ concentration in 35-ml regions of interest (ROIs) in the occipital lobe was found to be 1.4 +/- 0.2 mM, with scan-rescan repeatability of 38%.