Improved matrix inversion in image plane parallel MRI

A new 3D parallel magnetic resonance imaging (MRI) method named Generalized Unaliasing Incorporating Support constraint and sensitivity Encoding (GUISE) is presented. GUISE allows direct image recovery from arbitrary Cartesian k-space trajectories. However, periodic k-space sampling patterns are considered for reconstruction efficiency. Image recovery methods such as 2D SENSE (SENSitivity Encoding) and 2D CAIPIRINHA (Controlled Aliasing In Parallel Imaging Results IN Higher Acceleration) are special instances of GUISE where specific restrictions are placed on the k-space sampling patterns used. It is shown that the sampling pattern has large impacts on the image reconstruction error due to noise. An efficient sampling pattern design method that incorporates prior knowledge of object support and coil sensitivity profile is proposed. It requires no experimental trials and could be used in clinical imaging. Comparison of the proposed sampling pattern design method with 2D SENSE and 2D CAIPIRINHA are made based on both simulation and experiment results. It is seen that this new adaptive sampling pattern design method results in a lower noise level in reconstructions due to better exploitation of the coil sensitivity variation and object support constraint. In addition, elimination of the non-object region from reconstruction potentially allows an acceleration factor higher than the number of receiver coils used.     

NMR at very low fields

Although nuclear magnetic resonance in low fields around or below the Earth's magnetic field is almost as old as nuclear magnetic resonance itself, the recent years have experienced a revival of this technique that is opposed to the common trend towards higher and higher fields. The background of this development is the expectation that the low-field domain may open a new window for the study of molecular structure and dynamics. Here, we will give an overview on the specific features in the low-field domain, both from the technical and from the physical point of view. In addition, we present a short passage on the option of magnetic resonance imaging in fields of the micro-Tesla range.     

Noise estimation in parallel MRI: GRAPPA and SENSE

Parallel imaging methods allow to increase the acquisition rate via subsampled acquisitions of the k-space. SENSE and GRAPPA are the most popular reconstruction methods proposed in order to suppress the artifacts created by this subsampling. The reconstruction process carried out by both methods yields to a variance of noise value which is dependent on the position within the final image. Hence, the traditional noise estimation methods – based on a single noise level for the whole image – fail. In this paper we propose a novel methodology to estimate the spatial dependent pattern of the variance of noise in SENSE and GRAPPA reconstructed images. In both cases, some additional information must be known beforehand: the sensitivity maps of each receiver coil in the SENSE case and the reconstruction coefficients for GRAPPA.     

MRI with an atomic magnetometer suitable for practical imaging applications

Conventionally implemented MRI is performed in a strong magnetic field, typically >1 T. The high fields, however, can lead to many limitations. To overcome these limitations, ultra-low field (ULF) [or microtesla] MRI systems have been proposed and implemented. To-date such systems rely on low-Tc Superconducting Quantum Interference Devices (SQUIDs) leading to the requirement of cryogens. In this letter, we report ULF–MRI obtained with a non-cryogenic atomic magnetometer. This demonstration creates opportunities for developing inexpensive and widely applicable MRI scanners.     

基于同时多层激发和并行成像的心脏磁共振电影成像

磁共振成像（MRI）无创无害、对比度多、可以任意剖面成像的特点特别适合用于心脏成像,却因扫描时间长限制了其在临床上的应用.为了解决心脏磁共振电影成像屏气扫描时间过长的问题,该文提出了一种基于同时多层激发的多倍加速心脏磁共振电影成像及其影像重建的方法,该方法将相位调制多层激发（CAIPIRINHA）技术与并行加速（PPA）技术相结合,运用到分段采集心脏电影成像序列中,实现了在相位编码方向和选层方向的四倍加速,并使用改进的SENSE/GRAPPA算法对图像进行重建.分别在水模以及人体上进行了实验,将加速序列图像与不加速序列图像进行对比,结果验证了重建算法的有效性,表明该方法可以在保障图像质量以及准确测量心脏功能的前提下成倍节省扫描时间.     

Evaluation of motion effects on parallel MR imaging with precalibration

Several parallel imaging techniques such as SMASH, SENSE, k-space inherited parallel acquisition (KIPA) and others use reference (calibration) scans to find the parameters required for image reconstruction. Reference data is used to estimate coil sensitivity profiles for image domain techniques such as SENSE or reconstruction coefficients for k-space domain methods such as SMASH and KIPA. Any motion between the reference and accelerated imaging scans can make the reconstruction coefficients determined from the reference scan data suboptimal, resulting in an artifactual reconstruction. This work aims at comparing the effects of motion on the performance of three parallel imaging methods: SENSE, variable-density SENSE and KIPA, which all require one or more reference scans for calibration.     

Truncation effects in SENSE reconstruction

Finite sampling is an important practical issue in Fourier imaging systems. Although data truncation effects are well understood in conventional Fourier imaging where a single uniform receiver channel is used for data acquisition, this issue is not yet fully addressed in parallel imaging where an array of nonuniform receiver channels is used for sensitivity encoding to enable sub-Nyquist sampling of k-space. This article presents a systematic analysis of the problem by comparing the truncation effects in parallel imaging with those in conventional Fourier imaging. Specifically, it derives a convolution kernel function to characterize the truncation effects, which is shown to be approximately equal to that associated with the conventional Fourier imaging scheme. This article also describes a set of conditions under which significant differences between the truncation effects in parallel imaging and conventional Fourier imaging occur. The results should provide useful insight into interpreting and reducing data truncation effects in parallel imaging.     

Correction of concomitant gradient artifacts in experimental microtesla MRI

Magnetic resonance imaging (MRI) suffers from artifacts caused by concomitant gradients when the product of the magnetic field gradient and the dimension of the sample becomes comparable to the static magnetic field. To investigate and correct for these artifacts at very low magnetic fields, we have acquired MR images of a 165-mm phantom in a 66-microT field using gradients up to 350 microT/m. We prepolarize the protons in a field of about 100 mT, apply a spin-echo pulse sequence, and detect the precessing spins using a superconducting gradiometer coupled to a superconducting quantum interference device (SQUID). Distortion and blurring are readily apparent at the edges of the images; by comparing the experimental images to computer simulations, we show that concomitant gradients cause these artifacts. We develop a non-perturbative, post-acquisition phase correction algorithm that eliminates the effects of concomitant gradients in both the simulated and the experimental images. This algorithm assumes that the switching time of the phase-encoding gradient is long compared to the spin precession period. In a second technique, we demonstrate that raising the precession field during phase encoding can also eliminate blurring caused by concomitant phase-encoding gradients; this technique enables one to correct concomitant gradient artifacts even when the detector has a restricted bandwidth that sets an upper limit on the precession frequency. In particular, the combination of phase correction and precession field cycling should allow one to add MRI capabilities to existing 300-channel SQUID systems used to detect neuronal currents in the brain because frequency encoding could be performed within the 1-2 kHz bandwidth of the readout system.     

SNR and functional sensitivity of BOLD and perfusion-based fMRI using arterial spin labeling with spiral SENSE at 3 T

Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) studies using parallel imaging to reduce the readout window have reported a loss in temporal signal-to-noise ratio (SNR) that is less than would be expected given a purely thermal noise model. In this study, the impact of parallel imaging on the noise components and functional sensitivity of both BOLD and perfusion-based fMRI data was investigated. Dual-echo arterial spin labeling data were acquired on five subjects using sensitivity encoding (SENSE), at reduction factors (R) of 1, 2 and 3. Direct recording of cardiac and respiratory activity during data acquisition enabled the retrospective removal of physiological noise. The temporal SNR of the perfusion time series closely followed the thermal noise prediction of a √R loss in SNR as the readout window was shortened, with temporal SNR values (relative to the R=1 data) of 0.72 and 0.56 for the R=2 and R=3 data, respectively, after accounting for physiological noise. However, the BOLD temporal SNR decreased more slowly than predicted even after accounting for physiological noise, with relative temporal SNR values of 0.80 and 0.63 for the R=2 and R=3 data, respectively. Spectral analysis revealed that the BOLD trends were dominated by low-frequency fluctuations, which were not dominant in the perfusion data due to signal processing differences. The functional sensitivity, assessed using mean F values over activated regions of interest (ROIs), followed the temporal SNR trends for the BOLD data. However, results for the perfusion data were more dependent on the threshold used for ROI selection, most likely due to the inherently low SNR of functional perfusion data.     

Magnetic resonance imaging at frequencies below 1 kHz

Within the magnetic resonance imaging (MRI) community the trend is going to higher and higher magnetic fields, ranging from 1.5 T to 7 T, corresponding to Larmor frequencies of 63.8–298 MHz. Since for high-field MRI the magnetization increases with the applied magnetic field, the signal-to-noise-ratio increases as well, thus enabling higher image resolutions. On the other hand, MRI is possible also at ultra-low magnetic fields, as was shown by different groups. The goal of our development was to reach a Larmor frequency range of the low-field MRI system corresponding to the frequency range of human brain activities ranging from near zero-frequency (near-DC) to over 1 kHz. Here, first 2D MRI images of phantoms taken at Larmor frequencies of 100 Hz and 731 Hz will be shown and discussed. These frequencies are examples of brain activity triggered by electrostimulation of the median nerve. The method will allow the magnetic fields of the brain currents to influence the magnetic resonance image, and thus lead to a direct functional imaging modality of neuronal currents.     

MRI gradient fields increase brain mannitol space

Following nephrectomy and intravenous injection of tritiated mannitol, adult male rats were exposed to magnetic resonance imaging (MRI) procedures at 1.5 T, 0.5 T, and 0.3 T. Compared to rats similarly handled but not exposed to MRI procedures, brain mannitol concentration, expressed as a percentage of mean body concentration, was significantly increased at 0.3 T and 0.5 T but not at 1.5 T. At 0.3 T, exposure to gradient-field fluctuations used for imaging increased brain mannitol concentration, but exposures to static main field and pulsed radiofrequency energies did not. Increased brain mannitol associated with gradient-field flux may reflect increased blood-brain barrier permeability or blood volume in brain. MRI effects on brain mannitol space are of uncertain clinical significance, but are consistent with prior evidence of an MRI-induced increase of brain capillary endothelial cell transport observed with horseradish peroxidase. Further studies are needed to confirm these findings and to explore the processes underlying changes in mannitol distribution related to MRI.     

Quantification of microporosity in fruit by MRI at various magnetic fields: comparison with X-ray microtomography

Microstructure determines the mechanical and transport properties of fruit tissues. One important characteristic of the microstructure is the relative volume fraction of gas-filled intercellular spaces, i.e., the tissue microporosity. Quantification of this microporosity is fundamental for investigating the relationship between gas transfer and various disorders in fruit.     

Peculiarities of rf SQUID response in finite magnetic fields

Single-layer washer-type high-T_c YBa₂Cu₃O_7−x rf SQUIDs with grain-boundary Josephson junctions, as well as low-T_c Nb rf SQUIDs with Nb–Al₂O₃–Nb tunnel junctions, have been investigated in finite magnetic fields. It was shown experimentally that the suppression of the critical current of the Josephson junction due to the magnetic field leads to a modulation of the amplitude of the SQUID output signal. The role of the “unwanted” junction in high-T_c rf SQUIDs, which is formed by the grain boundary running through the washer of the SQUIDs on bicrystal substrates, has also been clarified. The drop of the SQUID signal at a finite magnetic field is originated by the penetration of the magnetic field into the unwanted junction. Based on these results, a direct radio-frequency method for the determination of the first critical field H_c1 for long Josephson junctions has been developed.     

Intracranial microvascular imaging at 7 T MRI with transceiver RF coils

Purpose

To investigate intracranial microvascular images with transceiver radio-frequency (RF) coils at ultra-high field 7 T magnetic resonance imaging (MRI).

Materials and methods

We designed several types of RF coils for the study of 7 T magnetic resonance angiography and analyzed quantitatively each coil's performance in terms of the signal-to-noise ratio (SNR) profiles to evaluate the usefulness of RF coils for microvascular imaging applications. We also obtained the microvascular images with different resolutions and parallel imaging technique.

Results

The overlapped 6-channel (ch) transceiver coil exhibited the highest performance for angiographic imaging. Although other multi-channel coils, such as 4- or 8-ch, were also suitable for fast imaging, these coils performed poorly in homogeneity or SNR for angiographic imaging. Furthermore, the 8-ch coil was poor in SNR at the center of the brain, while it had the highest SNR at the periphery.

Conclusion

The present study has demonstrated that the overlapped 6-ch coil with large-size loop coils provided the best performance for microvascular imaging or angiography with the ultra-high-field 7 T MRI, mainly because of its long penetration depth together with high SNR.     

Calculated signal-to-noise ratio of MRI detected with SQUIDs and Faraday detectors in fields from 10 microT to 1.5 T

We examine the calculated signal-to-noise ratio (SNR) achievable with different MRI detection modalities in precession fields ranging from 10 microT to 1.5 T. In particular, we compare traditional Faraday detectors with both tuned and untuned detectors based on superconducting quantum interference devices (SQUIDs). We derive general expressions for the magnetic field noise due to the samples and the detectors, and then calculate the SNR achievable for a specific geometry with each modality with and without prepolarization. We show that each of the three modalities is superior in one of the three field ranges. SQUID-based detection is superior to conventional Faraday detection for MRI in precession fields below 250 mT for a 65 mm diameter surface coil placed a distance of 25 mm from the voxel of interest embedded in a cylinder of tissue 50 mm tall and of radius 50 mm. This crossover field, however, is sensitive to the geometry.     

并行MRI图像重建算法比较及软件实现

首先介绍了不加速的并行MRI图像重建方法,然后对加速的并行MRI的4种图像重建算法进行了比较,得出结论：加速因子相同时,重建质量上,GRAPPA和SENSE的重建质量最好,SMASH的重建质量次之, PILS算法对线圈位置要求极高,重建质量最差;重建速度上,SMASH的重建速度最快,其次是SENSE和PILS,GRAPPA的重建速度最慢. 当加速因子变大时,所有算法重建质量都变差. 最后介绍了算法实现软件,该软件可以读入原始数据,显示数据采集轨迹,计算线圈灵敏度,选择图像重建方法,分析和比较重建图像质量. 该软件为我国在MRI成像领域提供了一个学习和进一步研究图像重建算法的有力工具.     

Reducing SAR in parallel excitation using variable-density spirals: a simulation-based study

Parallel excitation using multiple transmit channels has emerged as an effective method to shorten multidimensional spatially selective radiofrequency (RF) pulses, which have a number of important applications, including B1 field inhomogeneity correction in high-field MRI. The specific absorption rate (SAR) is a primary concern in high-field MRI, where wavelength effects can lead to local peaks in SAR. In parallel excitation, the subjects are exposed to RF pulses from multiple coils, which makes the SAR problem more complex to analyze, yet potentially enables greater freedom in designing RF pulses with lower SAR. Parallel-excitation techniques typically employ either Cartesian or constant-density (CD) spiral trajectories. In this article, variable-density (VD) spiral trajectories are explored as a means for SAR reduction in parallel-excitation pulse design. Numerical simulations were conducted to study the effects of CD and VD spirals on parallel excitation. Specifically, the electromagnetic fields of a four-channel transmit head coil with a three-dimensional head model at 4.7 T were simulated using a finite-difference time domain method. The parallel RF pulses were designed and the resulting excitation patterns were generated using a Bloch simulator. The SAR distributions due to CD and VD spirals were evaluated quantitatively. The simulation results show that, for the same pulse duration, parallel excitation with VD spirals can achieve a lower SAR compared to CD spirals for parallel excitation. VD spirals also resulted in reduced artifact power in the excitation patterns. This gain came with slight, but noticeable, degrading of the spatial resolution of the resulting excitation patterns.     

Parallel imaging with 3D TPI trajectory: SNR and acceleration benefits

Three-dimensional (3D) twisted projection imaging (TPI) trajectory has a unique advantage in sodium (²³Na) imaging on clinical MRI scanners at 1.5 or 3 T, generating a high signal-to-noise ratio (SNR) with a short acquisition time (∼10 min). Parallel imaging with an array of coil elements transits SNR benefits from small coil elements to acquisition efficiency by sampling partial k-space. This study investigates the feasibility of parallel sodium imaging with emphases on SNR and acceleration benefits provided by the 3D TPI trajectory. Computer simulations were used to find available acceleration factors and noise amplification. Human head studies were performed on clinical 1.5/3-T scanners with four-element coil arrays to verify simulation outcomes. In in vivo studies, proton (¹H) data, however, were acquired for concept–proof purpose. The sensitivity encoding (SENSE) method with the conjugate gradient algorithm was used to reconstruct images from accelerated TPI-SENSE data sets. Self-calibration was employed to estimate coil sensitivities. Noise amplification in TPI-SENSE was evaluated using multiple noise trials. It was found that the acceleration factor was as high as 5.53 (corresponding to acceleration number 2×3, ring-by-rotation), with a small image error of 6.9% when TPI projections were reduced in both polar (ring) and azimuthal (rotation) directions. The average noise amplification was as low as 98.7%, or 27% lower than Cartesian SENSE at that acceleration factor. The 3D nature of both TPI trajectory and coil sensitivities might be responsible for the high acceleration and low noise amplification. Consequently, TPI-SENSE may have potential advantages for parallel sodium imaging.