FPGA implementation of real-time SENSE reconstruction using pre-scan and Emaps sensitivities

Sensitivity Encoding (SENSE) is a widely used technique in Parallel Magnetic Resonance Imaging (MRI) to reduce scan time. Reconfigurable hardware based architecture for SENSE can potentially provide image reconstruction with much less computation time. Application specific hardware platform for SENSE may dramatically increase the power efficiency of the system and can decrease the execution time to obtain MR images. A new implementation of SENSE on Field Programmable Gate Array (FPGA) is presented in this study, which provides real-time SENSE reconstruction right on the receiver coil data acquisition system with no need to transfer the raw data to the MRI server, thereby minimizing the transmission noise and memory usage. The proposed SENSE architecture can reconstruct MR images using receiver coil sensitivity maps obtained using pre-scan and eigenvector (E-maps) methods. The results show that the proposed system consumes remarkably less computation time for SENSE reconstruction, i.e., 0.164 ms @ 200 MHz, while maintaining the quality of the reconstructed images with good mean SNR (29 + dB), less RMSE (< 5 × 10^− 2) and comparable artefact power (< 9 × 10^− 4) to conventional SENSE reconstruction. A comparison of the center line profiles of the reconstructed and reference images also indicates a good quality of the reconstructed images. Furthermore, the results indicate that the proposed architectural design can prove to be a significant tool for SENSE reconstruction in modern MRI scanners and its low power consumption feature can be remarkable for portable MRI scanners.     

3D sensitivity encoded ellipsoidal MR spectroscopic imaging of gliomas at 3T

Purpose

The goal of this study was to implement time efficient data acquisition and reconstruction methods for 3D magnetic resonance spectroscopic imaging (MRSI) of gliomas at a field strength of 3T using parallel imaging techniques.

Methods

The point spread functions, signal to noise ratio (SNR), spatial resolution, metabolite intensity distributions and Cho:NAA ratio of 3D ellipsoidal, 3D sensitivity encoding (SENSE) and 3D combined ellipsoidal and SENSE (e-SENSE) k-space sampling schemes were compared with conventional k-space data acquisition methods.

Results

The 3D SENSE and e-SENSE methods resulted in similar spectral patterns as the conventional MRSI methods. The Cho:NAA ratios were highly correlated (P<.05 for SENSE and P<.001 for e-SENSE) with the ellipsoidal method and all methods exhibited significantly different spectral patterns in tumor regions compared to normal appearing white matter. The geometry factors ranged between 1.2 and 1.3 for both the SENSE and e-SENSE spectra. When corrected for these factors and for differences in data acquisition times, the empirical SNRs were similar to values expected based upon theoretical grounds. The effective spatial resolution of the SENSE spectra was estimated to be same as the corresponding fully sampled k-space data, while the spectra acquired with ellipsoidal and e-SENSE k-space samplings were estimated to have a 2.36–2.47-fold loss in spatial resolution due to the differences in their point spread functions.

Conclusion

The 3D SENSE method retained the same spatial resolution as full k-space sampling but with a 4-fold reduction in scan time and an acquisition time of 9.28 min. The 3D e-SENSE method had a similar spatial resolution as the corresponding ellipsoidal sampling with a scan time of 4:36 min. Both parallel imaging methods provided clinically interpretable spectra with volumetric coverage and adequate SNR for evaluating Cho, Cr and NAA.     

Design and Demonstration of Four-Channel Received Coil Arrays for Vertical-Field MRI

The layout of radio-frequency received coils is related to signal-to-noise ratio (SNR) in magnetic resonance imaging (MRI). In this paper, different structures of four-channel received coil arrays for vertical-field MRI are constructed and optimized by establishing the relationship between coil geometry and SNR to achieve a high SNR and a uniform SNR distribution in the region of interest (ROI). Then, the SNR distributions of three optimized configurations, including rectangular loops, non-definite shape surface coils, and solenoid loops as the main unit, are simulated and compared. The four-channel coil of solenoid loops as the main unit has been found to have the best performance with the highest mean SNR in the ROI when imaging without acceleration. In addition, g-factor and 2D SENSE SNR in yoz-plane are simply analyzed, which show a sharp decrease in SNR for all the coils. Finally, all the coils are manufactured and operated at a 0.5 T permanent magnet MRI system with phantom and joint imaging experiments. Using pixel-by-pixel manner to evaluate SNR map, the experimental results are consistent with the simulation results, while parallel imaging experiment results show that the major consideration in low field MRI is the improvement of SNR value and uniformity rather than that of the imaging speed. As different constructions of four-channel received coils are investigated, we have found the most effective configuration with high and uniform SNR for vertical-field MRI.     

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.     

Localized high-resolution DTI of the human midbrain using single-shot EPI,parallel imaging,and outer-volume suppression at 7 T

Localized high-resolution diffusion tensor images (DTI) from the midbrain were obtained using reduced field-of-view (rFOV) methods combined with SENSE parallel imaging and single-shot echo planar (EPI) acquisitions at 7 T. This combination aimed to diminish sensitivities of DTI to motion, susceptibility variations, and EPI artifacts at ultra-high field. Outer-volume suppression (OVS) was applied in DTI acquisitions at 2- and 1-mm² resolutions, b = 1000 s/mm², and six diffusion directions, resulting in scans of 7- and 14-min durations. Mean apparent diffusion coefficient (ADC) and fractional anisotropy (FA) values were measured in various fiber tract locations at the two resolutions and compared. Geometric distortion and signal-to-noise ratio (SNR) were additionally measured and compared for reduced-FOV and full-FOV DTI scans. Up to an eight-fold data reduction was achieved using DTI-OVS with SENSE at 1 mm², and geometric distortion was halved. The localization of fiber tracts was improved, enabling targeted FA and ADC measurements. Significant differences in diffusion properties were observed between resolutions for a number of regions suggesting that FA values are impacted by partial volume effects even at a 2-mm² resolution. The combined SENSE DTI-OVS approach allows large reductions in DTI data acquisition and provides improved quality for high-resolution diffusion studies of the human brain.     

Potential advantage of higher-order modes of birdcage coil for parallel imaging

The potential advantages of the higher-order resonant modes of a low-pass birdcage coil for parallel imaging are investigated using FDTD simulations in a spherical phantom. Better parallel imaging performance can be achieved in axial planes than in sagittal planes. If more modes are employed, the average g-factor (gmean) and maximum g-factor (gmax) will be improved for a specific acceleration factor (R). G-factor performance at 3 T (gmean=1.79, gmax=3.55) and 7 T (gmean=1.60, gmax=2.45) can be achieved even with an acceleration factor as high as four when all six order modes of 12-rung low-pass birdcage coil are incorporated for imaging on a mid-axial plane. For a specific number of channels, the optimum combination of corresponding modes can be obtained for different acceleration factors. Based on the g-factor and SNR performance, the total degenerate multi-mode birdcage coil with six order resonant modes has better homogeneous coverage and SENSE performance than the 8-element phased array coil, although requiring fewer channels. In addition, the dielectric effects at high field can improve the parallel imaging performance.     

An optimization method for designing SENSE imaging RF coil arrays

An optimization method in RF coil array design for SENSE imaging is described. Using this method the optimized RF coil geometries can be calculated numerically given the required SENSE imaging performance. Although this method can be applied to optimize the RF coil arrays for both 1D and 2D SENSE imaging, to demonstrate the potential applications of this method, we designed RF coil arrays for 2D SENSE imaging and compared their performance by simulation. An optimized 4-channel receive-only RF coil array designed for 2D SENSE imaging was implemented and tested to demonstrate the feasibility of the proposed technique. Imaging results showed reasonable agreement with the simulations, thus the method can be applied to RF coil array designs for SENSE imaging when optimum imaging performance is desired.     

Multichannel transceiver dual-tuned RF coil for proton/sodium MR imaging of knee cartilage at 3 T

Sodium magnetic resonance (MR) imaging is a promising technique for detecting changes of proteoglycan (PG) content in cartilage associated with knee osteoarthritis. Despite its potential clinical benefit, sodium MR imaging in vivo is challenging because of intrinsically low sodium concentration and low MR signal sensitivity. Some of the challenges in sodium MR imaging may be eliminated by the use of a high-sensitivity radiofrequency (RF) coil, specifically, a dual-tuned (DT) proton/sodium RF coil which facilitates the co-registration of sodium and proton MR images and the evaluation of both physiochemical and structural properties of knee cartilage. Nevertheless, implementation of a DT proton/sodium RF coil is technically difficult because of the coupling effect between the coil elements (particularly at high field) and the required compact design with improved coil sensitivity. In this study, we applied a multitransceiver RF coil design to develop a DT proton/sodium coil for knee cartilage imaging at 3 T. With the new design, the size of the coil was minimized, and a high signal-to-noise ratio (SNR) was achieved. DT coil exhibited high levels of reflection S11 (～-21 dB) and transmission coefficient S12 (～-19 dB) for both the proton and sodium coils. High SNR (range 27-38) and contrast-to-noise ratio (CNR) (range 15-21) were achieved in sodium MR imaging of knee cartilage in vivo at 3-mm(3) isotropic resolution. This DT coil performance was comparable to that measured using a sodium-only birdcage coil (SNR of 28 and CNR of 20). Clinical evaluation of the DT coil on four normal subjects demonstrated a consistent acquisition of high-resolution proton images and measurement of relative sodium concentrations of knee cartilages without repositioning of the subjects during the same MR scanning session.     

A 16-channel loop array for in vivo macaque whole-brain imaging at 3 T

Non-human primates (NHPs) are vital models for neuroscience research. These animals have been widely used in behavioral, electrophysiological, molecular, and more recently, multimodal neuroimaging and neuro-engineering studies. Several RF coil arrays have been designed for functional, high-resolution brain magnetic resonance imaging (MRI), but few have been designed to accommodate multimodal devices. In the present study, a 16-channel array coil was constructed for brain imaging of macaques at 3 Tesla (3 T). To construct this coil, a close-fitting helmet-shaped form was designed to host 16 coil loops for whole-brain coverage. This assembly is mountable onto stereotaxic head frame bars, and the coil functions while the monkey is in the sphinx position with a clear line of vision of stimuli presented from outside of the MRI system. In addition, 4 openings were allocated in the coil housing, allowing multimodal devices to directly access visual cortical regions such as V1-V4 and MT. Coil performance was evaluated in an anesthetized macaque by quantifying and comparing signal-to-noise ratios (SNRs), noise correlations, and g-factor maps to a vendor-supplied human pediatric coil frequently used for NHP MRI. The result from in vivo experiments showed that the NHP coil was well-decoupled, had higher SNRs in cortical regions, and improved data acquisition acceleration capability compared with a vendor-supplied human pediatric coil that has been frequently used in macaque MRI studies. Furthermore, whole-brain anatomic imaging, diffusion tensor imaging and functional brain imaging have also been conducted: the details of brain anatomical structure, such as cerebellum and brainstem, can be clearly visualized in T2-SPACE images; b0 SNR calculated from b0 maps was higher than the human pediatric coil in all regions of interest (ROIs); the time-course SNR (tSNR) map calculated for GRE-EPI images demonstrates that the presented coil can be used for high-resolution functional imaging at 3 T.     

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.     

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.     

Comparison of different methods for combining phase-contrast images obtained with multiple coils

The ability to determine coil sensitivities implies that a method optimized in terms of maximized signal-to-noise ratio (SNR) can be applied to the combination of multiple coil images. An optimization of SNR subsequently results in a minimized variance in quantitative velocity measurements using phase-contrast imaging. When coil sensitivities are unknown, the weighted mean method, utilizing the square of the signal magnitude as weights, is suitable for combination of multiple phase images. In this study, the optimized method using estimated coil sensitivities was compared to the weighted mean method both theoretically and experimentally. It is shown that absence of noise correlation between the different coil images implies no difference between the methods regarding the variance of the phase. In the practical situation, noise correlation does exist, implying an opportunity for further reduction of phase variance using the optimized method. In vitro and in vivo studies showed, however, no significant difference between the two methods studied.     

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.     

Parallel magnetic resonance imaging using coils with localized sensitivities

The purpose of this study was to present clinical examples and illustrate the inefficiencies of a conventional reconstruction using a commercially available phased array coil with localized sensitivities. Five patients were imaged at 1.5 T using a cardiac-synchronized gadolinium-enhanced acquisition and a commercially available four-element phased array coil. Four unique sets of images were reconstructed from the acquired k-space data: (a) sum-of-squares image using four elements of the coil; localized sum-of-squares images from the (b) anterior coils and (c) posterior coils and a (c) local reconstruction. Images were analyzed for artifacts and usable field-of-view. Conventional image reconstruction produced images with fold-over artifacts in all cases spanning a portion of the image (mean 90 mm; range 36-126 mm). The local reconstruction removed fold-over artifacts and resulted in an effective increase in the field-of-view (mean 50%; range 20-70%). Commercially available phased array coils do not always have overlapping sensitivities. Fold-over artifacts can be removed using an alternate reconstruction method. When assessing the advantages of parallel imaging techniques, gains achieved using techniques such as SENSE and SMASH should be gauged against the acquisition time of the localized method rather than the conventional sum-of-squares method.     

APIR4EMC: Autocalibrated parallel imaging reconstruction for extended multi-contrast imaging

PurposeTo improve image quality of multi-contrast imaging with the proposed Autocalibrated Parallel Imaging Reconstruction for Extended Multi-Contrast Imaging (APIR4EMC).MethodsAPIR4EMC reconstructs multi-contrast images in an autocalibrated parallel imaging reconstruction framework by adding contrasts as virtual coils. Compensation of signal evolution along the echo train of different contrasts is performed to improve signal prediction for missing samples. As a proof of concept, we performed prospectively accelerated phantom and in-vivo brain acquisitions with T1, T1-fat saturated (Fatsat), T2, PD, and FLAIR contrasts. The k-space sampling patterns of these acquisitions were jointly optimized. Images were jointly reconstructed with the proposed APIR4EMC method as well as individually with GRAPPA. Root mean square error (RMSE) to fully sampled reference images and g-factor maps were computed for both methods in the phantom experiment. Visual evaluation was performed in the in-vivo experiment.ResultsCompared to GRAPPA, APIR4EMC reduced artifacts and improved SNR of the reconstructed images in the phantom acquisitions. Quantitatively, APIR4EMC substantially reduced noise amplification (g-factor) as well as RMSE compared to GRAPPA. Signal evolution compensation reduced artifacts. In the in-vivo experiments, 1 mm³ isotropic 3D images with contrasts of T1, T1-Fatsat, T2, PD, and FLAIR were acquired in as little as 7.5 min with the acceleration factor of 9. Reconstruction quality was consistent with the phantom results.ConclusionCompared to single contrast reconstruction with GRAPPA, APIR4EMC reduces artifacts and noise amplification in accelerated multi-contrast imaging.     

Optimizing Image Reconstruction in SENSE Using GPU

Parallel magnetic resonance imaging (MRI) (pMRI) uses multiple receiver coils to reduce the MRI scan time. To accelerate the data acquisition process in MRI, less amount of data is acquired from the scanner which leads to artifacts in the reconstructed images. SENSitivity Encoding (SENSE) is a reconstruction algorithm in pMRI to remove aliasing artifacts from the undersampled multi coil data and recovers fully sampled images. The main limitation of SENSE is computing inverse of the encoding matrix. This work proposes the inversion of encoding matrix using Jacobi singular value decomposition (SVD) algorithm for image reconstruction on GPUs to accelerate the reconstruction process. The performance of Jacobi SVD is compared with Gauss–Jordan algorithm. The simulations are performed on two datasets (brain and cardiac) with acceleration factors 2, 4, 6 and 8. The results show that the graphics processing unit (GPU) provides a speed up to 21.6 times as compared to CPU reconstruction. Jacobi SVD algorithm performs better in terms of acceleration in reconstructions on GPUs as compared to Gauss–Jordan method. The proposed algorithm is suitable for any number of coils and acceleration factors for SENSE reconstruction on real time processing systems.     

A flexible 9-channel coil array for fast 3D MR thermometry in MR-guided high-intensity focused ultrasound (HIFU) studies on rabbits at 3 T

Signal-to-noise ratio (SNR) is a critical factor in MR-guided high-intensity focused ultrasound (HIFU) for local heating, which can affect the accuracy of temperature measurement. In order to achieve high SNR and higher temporal resolution, dedicated coil arrays for MR-guided HIFU applications need to be developed. In this work, a flexible 9-channel coil array was designed, and constructed at 3 T to achieve fast temperature mapping for MR-guided HIFU applications on rabbit leg muscle. Coil performance was evaluated for SNR, and parallel imaging capability by in-vivo studies. Compared to a commercially available 4-channel flexible coil array, the dedicated 9-channel coil array has a much higher SNR, with at least a 2.6-fold increment in the region of interest (ROI). The inverse g-factors maps demonstrated that the dedicated 9-channel coil array has a better parallel imaging capability than the Flex Small 4. With accelerations normal to the array direction, both coil arrays showed much higher g-factors than those of accelerations along the array direction. Room temperature mapping was implemented to evaluate the temperature measurement accuracy by in-vivo experiments. The precisions of the 9-channel coil, ±0.18 °C for un-acceleration and ± 0.56 °C for acceleration at R = 2 × 2, both improved by an order of magnitude than these of the 4-channel coil, which were ± 1.45 °C for un-acceleration and ± 3.52 °C for acceleration at R = 2 × 2. In the fast temperature imaging on the rabbit leg muscle with heating, a high temporal resolution of 3.3 s with a temperature measurement precision of ±0.56 °C has been achieved using the dedicated 9-channel coil. This study demonstrates that the dedicated 9-channel coil array for rabbit leg imaging provides improved performance in SNR, parallel imaging capability, and the accuracy of temperature measurement compared to a commercial 4-channel coil, and it also achieves fast temperature mapping in practical MR-guided HIFU applications.     

Helmholtz-pair transmit coil with integrated receive array for high-resolution MRI of trabecular bone in the distal tibia at 7T

A Helmholtz-pair local transmit RF coil with an integrated four-element receive array RF coil and foot immobilization platform was designed and constructed for imaging the distal tibia in a whole-body 7T MRI scanner. Simulations and measurements of the B(1) field distribution of the transmit coil are described, along with SAR considerations for operation at 7T. Results of imaging the trabecular bone of three volunteers at 1.5T, 3T and 7T are presented, using identical 1.5T and 3T versions of the 7T four-element receive array. The spatially registered images reveal improved visibility for individual trabeculae and show average gains in SNR of 2.8× and 4.9× for imaging at 7T compared to 3T and 1.5T, respectively. The results thus display an approximately linear dependence of SNR with field strength and enable the practical utility of 7T scanners for micro-MRI of trabecular bone.