Dynamic magnetic resonance imaging of swallowing and laryngeal motion using parallel imaging at 3 T

Object

To evaluate the feasibility of an optimized MRI protocol based on high field imaging at 3 T in combination with accelerated data acquisition by parallel imaging for the analysis of oropharyngeal and laryngeal function.

Materials and Methods

Fast 2D gradient echo (GRE) MRI with different spatial resolutions (1.7×2.7 and 1.1×1.5 mm²) and image update rates (4 and 10 frames per second) was employed to assess pharyngeal movements and visualize swallowing via tracking of an oral contrast bolus (blueberry juice). In a study with 10 normal volunteers, image quality was semi-quantitatively graded by three independent observers with respect to the delineation of anatomical detail and depiction of oropharynx and larynx function. Additionally, the feasibility of the technique for the visualization of pathological pre- and post-surgical oropharynx and larynx function was evaluated in a patient with inspiratory stridor.

Results

Image grading demonstrated the feasibility of dynamic MRI for the assessment of normal oropharynx and larynx anatomy and function. Superior image quality (P<.05) was found for data acquisition with four frames per second and higher spatial resolution. In the patient, dynamic MRI detected pathological hypermobility of the epiglottis resulting in airway obstruction. Additional post-surgical MRI for one clinical case revealed morphological changes of the epiglottis and improved function, i.e., absence of airway obstruction and normal swallowing.

Conclusion

Results of the volunteer study demonstrated the feasibility of dynamic MRI at 3 T for the visualization of the oropharynx and larynx function during breathing, movements of the tongue and swallowing. Future studies are necessary to evaluate its clinical value compared to existing modalities based on endoscopy or radiographic techniques.     

Data consistency criterion for selecting parameters for k-space-based reconstruction in parallel imaging

k-space-based reconstruction in parallel imaging depends on the reconstruction kernel setting, including its support. An optimal choice of the kernel depends on the calibration data, coil geometry and signal-to-noise ratio, as well as the criterion used. In this work, data consistency, imposed by the shift invariance requirement of the kernel, is introduced as a goodness measure of k-space-based reconstruction in parallel imaging and demonstrated. Data consistency error (DCE) is calculated as the sum of squared difference between the acquired signals and their estimates obtained based on the interpolation of the estimated missing data. A resemblance between DCE and the mean square error in the reconstructed image was found, demonstrating DCE's potential as a metric for comparing or choosing reconstructions. When used for selecting the kernel support for generalized autocalibrating partially parallel acquisition (GRAPPA) reconstruction and the set of frames for calibration as well as the kernel support in temporal GRAPPA reconstruction, DCE led to improved images over existing methods. Data consistency error is efficient to evaluate, robust for selecting reconstruction parameters and suitable for characterizing and optimizing k-space-based reconstruction in parallel imaging.     

SNR-optimized myocardial perfusion imaging using parallel acquisition for effective density-weighted saturation recovery imaging

The concept of density-weighted imaging and parallel acquisition for effective density-weighted (PLANED) imaging was transferred to saturation recovery (SR) sequences, in order to increase the SNR in first-pass myocardial perfusion imaging. Filtering in combination with density-weighted imaging allows SNR-optimized data weighting and the free choice of the corresponding spatial response function (SRF) simultaneously. This method was evaluated in simulations and applied successfully to phantom and in vivo first-pass myocardial perfusion studies. Unfiltered, Cartesian sampled images were compared to images acquired with SR-PLANED, which has been adjusted to result in an identical SRF as the Cartesian imaging. SNR-optimized SR-PLANED imaging improved the SNR up to 15% without changing acquisition time, the SRF or the field of view (FOV). The presented method provides high image quality and optimized SNR for first-pass myocardial perfusion imaging.     

Ictal imaging using functional magnetic resonance

Magnetic resonance imaging (MRI) can now provide maps of human brain function with high spatial and temporal resolution. This noninvasive technique can also map the coritical activation that occurs during focal seizures, as demonstrated here by the results obtained using a conventional 1.5 T clinical MRI system for the investigation of a 4-year-old boy suffering from frequent partial motor seizures of his right side. FLASH images (TE = 60 ms) were acquired every 10 s over a period of 25 min, and activation images derived by subtracting baseline images from images obtained during clinical seizures. Functional MRI revealed sequential activation associated with specific gyri within the left hemisphere with each of five consecutive clinical seizures, and also during a period that was not associated with a detectable clinical seizure. The activated regions included gyri that were structurally abnormal. These results demonstrate (a) that functional MRI can potentially provide new insights into the dynamic events that occur in the epileptic brain and their relationship to brain structure; and (b) that there is the possibility of obtaining similar information in the absence of clinical seizures, suggesting the potential for studies in patients with interictal electrical disturbances.     

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

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

High resolution dynamic cardiac MRI using partial separability of spatiotemporal signals with a novel sampling scheme

The partial separability (PS) of spatiotemporal signals has been exploited to accelerate dynamic cardiac MRI by sampling two datasets (training and imaging datasets) without breath-holding or ECG triggering. According to the theory of partially separable functions, the wider the range of spatial frequency components covered by the training dataset, the more accurate the temporal constraint imposed by the PS model. Therefore, it is necessary to develop a new sampling scheme for the PS model in order to cover a wider range of spatial frequency components. In this paper, we propose the use of radial sampling trajectories for collecting the training dataset and Cartesian sampling trajectories for collecting the imaging dataset. In vivo high resolution cardiac MRI experiments demonstrate that the proposed data sampling scheme can significantly improve the image quality. The image quality using the PS model with the proposed sampling scheme is comparable to that of a commercial method using retrospective cardiac gating and breath-holding. The success of this study demonstrates great potential for high-quality, high resolution dynamic cardiac MRI without ECG gating or breath-holding through use of the PS model and the novel data sampling scheme.     

Influence of multichannel combination, parallel imaging and other reconstruction techniques on MRI noise characteristics

The statistical properties of background noise such as its standard deviation and mean value are frequently used to estimate the original noise level of the acquired data. This requires the knowledge of the statistical intensity distribution of the background signal, that is, the probability density of the occurrence of a certain signal intensity. The influence of many new MRI techniques and, in particular, of various parallel-imaging methods on the noise statistics has neither been rigorously investigated nor experimentally demonstrated yet. In this study, the statistical distribution of background noise was analyzed for MR acquisitions with a single-channel and a 32-channel coil, with sum-of-squares (SoS) and spatial-matched-filter (SMF) data combination, with and without parallel imaging using k-space and image-domain algorithms, with real-part and conventional magnitude reconstruction and with several reconstruction filters. Depending on the imaging technique, the background noise could be described by a Rayleigh distribution, a noncentral chi-distribution or the positive half of a Gaussian distribution. In particular, the noise characteristics of SoS-reconstructed multichannel acquisitions (with k-space-based parallel imaging or without parallel imaging) differ substantially from those with image-domain parallel imaging or SMF combination. These effects must be taken into account if mean values or standard deviations of background noise are employed for data analysis such as determination of local noise levels. Assuming a Rayleigh distribution as in conventional MR images or a noncentral chi-distribution for all multichannel acquisitions is invalid in general and may lead to erroneous estimates of the signal-to-noise ratio or the contrast-to-noise ratio.     

High-field magnetic resonance imaging using solenoid radiofrequency coils

High-resolution magnetic resonance imaging using dedicated high-field radiofrequency micro-coils at 16.4 T (700 MHz) was investigated. Specific solenoid coils primarily using silver and copper as conductors with enamel and polyurethane coatings were built to establish which coil configuration produces the best image. Image quality was quantified using signal-to-noise ratio and signal variation over regions of interest. Benchmarking was conducted using 5-mm diameter coils, as this size is comparable to an established coil of the same size. Our 1.4-mm-diameter coils were compared directly to each other, from which we deduce performance as a function of conductor material and coating. A variety of materials and conductor coatings allowed us to choose an optimal design, which we used to image a kidney section at 10-micron resolution. We applied zero-fill extrapolation to achieve 5-micron resolution.     

Surface-based functional magnetic resonance imaging analysis of partial brain echo planar imaging data at 1.5 T

Surface-based functional magnetic resonance imaging (fMRI) analysis is more sensitive and accurate than volume-based analysis for detecting neural activation. However, these advantages are less important in practical fMRI experiments with commonly used 1.5-T magnetic resonance devices because of the resolution gap between the echo planar imaging data and the cortical surface models. We expected high-resolution segmented partial brain echo planar imaging (EPI) data to overcome this problem, and the activation patterns of the high-resolution data could be different from the low-resolution data. For the practical applications of surface-based fMRI analysis using segmented EPI techniques, the effects of some important factors (e.g., activation patterns, registration and local distortions) should be intensively evaluated because the results of surface-based fMRI analyses could be influenced by them. In this study, we demonstrated the difference between activations detected from low-resolution EPI data, which were covering whole brain, and high-resolution segmented EPI data covering partial brain by volume- and surface-based analysis methods. First, we compared the activation maps of low- and high-resolution EPI datasets detected by volume- and surface-based analyses, with the spatial patterns of activation clusters, and analyzed the distributions of activations in occipital lobes. We also analyzed the high-resolution EPI data covering motor areas and fusiform gyri of human brain, and presented the differences of activations detected by volume- and surface-based methods.     

Composite MR image reconstruction and unaliasing for general trajectories using neural networks

In rapid parallel magnetic resonance imaging, the problem of image reconstruction is challenging. Here, a novel image reconstruction technique for data acquired along any general trajectory in neural network framework, called “Composite Reconstruction And Unaliasing using Neural Networks” (CRAUNN), is proposed. CRAUNN is based on the observation that the nature of aliasing remains unchanged whether the undersampled acquisition contains only low frequencies or includes high frequencies too. Here, the transformation needed to reconstruct the alias-free image from the aliased coil images is learnt, using acquisitions consisting of densely sampled low frequencies. Neural networks are made use of as machine learning tools to learn the transformation, in order to obtain the desired alias-free image for actual acquisitions containing sparsely sampled low as well as high frequencies. CRAUNN operates in the image domain and does not require explicit coil sensitivity estimation. It is also independent of the sampling trajectory used, and could be applied to arbitrary trajectories as well. As a pilot trial, the technique is first applied to Cartesian trajectory-sampled data. Experiments performed using radial and spiral trajectories on real and synthetic data, illustrate the performance of the method. The reconstruction errors depend on the acceleration factor as well as the sampling trajectory. It is found that higher acceleration factors can be obtained when radial trajectories are used. Comparisons against existing techniques are presented. CRAUNN has been found to perform on par with the state-of-the-art techniques. Acceleration factors of up to 4, 6 and 4 are achieved in Cartesian, radial and spiral cases, respectively.     

Sectoral sampling in centric-scan SPRITE magnetic resonance imaging

A new approach to the construction of k-space trajectories for centric-scan SPRITE in both 2D and 3D is presented. All benefits of previous SPRITE methods are retained, most importantly the ability to image objects with short T*(2). This new approach gives more flexibility in the choice of number of interleaves with points more evenly distributed across k-space. All these improvements positively contribute to image quality and resolution, which can be also traded off against experimental speed. Sectoral sampling will have significant benefits for magnetisation preparation contrast imaging.     

Improved simultaneous multislice magnetic resonance imaging using total variation regularization

Controlled aliasing in parallel imaging results in higher acceleration(CAIPIRINHA) for simultaneous multislice imaging has been proposed recently, which combines multiband excitation and phase cycling techniques to reduce scan time and improve subsequent imaging reconstruction. In this work, the total variation(TV) regularization method is used to further improve CAIPIRINHA. The TV regularization uses an edge-preserving prior, which establishes a relationship between neighboring pixels for image reconstruction. It reduces artifacts and suppresses noise amplification simultaneously.The results are presented with a standard eight-channel head coil with an acceleration factor of 4, where the TV-regularized CAIPIRINHA generates an improved reconstruction as compared with a typical nonregularized CAIPIRINHA.     

Labeling of macrophages using bacterial magnetosomes and their characterization by magnetic resonance imaging

This work investigated macrophages labeled with magnetosomes for the possible detection of inflammations by MR molecular imaging. Pure magnetosomes and macrophages containing magnetosomes were analyzed using a clinical 1.5 T MR-scanner. Relaxivities of magnetosomes and relaxation rates of cells containing magnetosomes were determined. Peritonitis was induced in two mice. T₁, T₂ and T₂* weighted images were acquired following injection of the probes. Pure magnetosomes and labeled cells showed slight effects on T₁, but strong effects on T₂ and T₂* images. Labeled macrophages were located with magnetic resonance imaging (MRI) in the colon area, thus demonstrating the feasibility of the proposed approach.     

Head and body imaging by hydrogen nuclear magnetic resonance

A hydrogen (¹H) nuclear magnetic resonance (NMR) imaging study of the normal head, thorax, and limbs is reported. The images are 10 to 15 mm thick transverse slices obtained in 2 to 4 min using a two-dimensional Fourier transform technique. Spatial resolution in the imaging plane is about 2 mm, enabling the optic nerve and many small blood vessels to be observed. Thorax scans show details of the cardiac chambers, aorta wall, and lungs without artefacts arising from physiological motion.     

Sub-optical resolution of single spins using magnetic resonance imaging at room temperature in diamond

There has been much recent interest in extending the technique of magnetic resonance imaging (MRI) down to the level of single spins with sub-optical wavelength resolution. However, the signal to noise ratio for images of individual spins is usually low and this necessitates long acquisition times and low temperatures to achieve high resolution. An exception to this is the nitrogen-vacancy (NV) color center in diamond whose spin state can be detected optically at room temperature. Here we apply MRI to magnetically equivalent NV spins and demonstrate fully resolved spectra with resolution well below the optical wavelength of the readout light. In addition, using a microwave version of MRI we achieved a resolution that is 1/270 in size of the coplanar striplines, which define the effective wavelength of the microwaves that were used to excite the transition. This technique can eventually be extended to imaging of large numbers of NVs in a confocal spot and possibly to image nearby dark spins via their mutual magnetic interaction with the NV spin.     

基于深度神经网络的时空编码磁共振成像超分辨率重建方法

单扫描时空编码磁共振成像是一种新型超快速磁共振成像技术,它对磁场不均匀和化学位移伪影有较强的抵抗性,但是其固有的空间分辨率较低,因此通常需要进行超分辨率重建,以在不增加采样点数的情况下提高时空编码磁共振图像的空间分辨率.然而,现有的重建方法存在迭代求解时间长、重建结果有混叠伪影残留等问题.为此,本文提出了一种基于深度神经网络的单扫描时空编码磁共振成像超分辨率重建方法.该方法采用模拟样本训练深度神经网络,再利用训练好的网络模型对实际采样信号进行重建.数值模拟、水模和活体鼠脑的实验结果表明,该方法能快速重建出无残留混叠伪影、纹理信息清楚的超分辨率时空编码磁共振图像.适当增加训练样本数量以及在训练样本中加入适当的随机噪声水平,有助于改善重建效果.     

Cardiac magnetic resonance imaging in arrhythmogenic right ventricular cardiomyopathy: correlation to the QRS dispersion

The aim of the study was to evaluate the relationship between the presence of right ventricular abnormalities detected by cardiac magnetic resonance (CMR) and QRS dispersion, the strongest independent predictor of sudden death in ARVC. A consecutive series of 40 patients from a single institution were recruited with a clinical diagnosis of ARVC based on the diagnostic criteria. All patients underwent systematic clinical evaluation, including history and examination, electrocardiography, 24-h Holter monitor, chest radiography, echocardiography and CMR examination and were divided into two groups according to the QRS dispersion: group I, QRS dispersion ≥ 40 ms; group II, QRS dispersion < 40 ms. The relationship between the characteristic parameters of CMR image and QRS dispersion were analyzed in two groups. There were significant differences in QRS dispersion (57±14 ms vs. 26±11 ms), right ventricular end-diastolic diameter (57±10 mm vs. 48±11 mm, P=.012), right ventricular end-systolic diameter (52±10 mm vs. 44±11 mm, P=.010), right ventricular end-diastolic volume (260±105 ml vs. 180±66 ml, P=.006), right ventricular end-systolic volume (222±98 ml vs. 148±61 ml, P=.006) and myocardial fibrosis detection rate (74% vs. 38%, P=.024) between two groups. For all patients with ARVC, QRS dispersion and right ventricular end-diastolic volume (r= 0.66, P<.001), right ventricular end-systolic volume (r= 0.67, P<.001), right ventricular outflow tract area (r= 0.68, P<.001) showed a moderate positive correlation. Right ventricular outflow tract area, right ventricular end-diastolic volume and end-systolic volume detected by CMR in patients with ARVC were positively correlated to the extent of QRS dispersion (≥ 40 ms), the strongest independent predictor of sudden cardiac death.