Golden-ratio rotated stack-of-stars acquisition for improved volumetric MRI

Efficient 3D cone trajectory design for improved combined angiographic and perfusion imaging using arterial spin labeling

PURPOSE

To improve the spatial resolution and repeatability of a non-contrast MRI technique for simultaneous time resolved 3D angiography and perfusion imaging by developing an efficient 3D cone trajectory design.

METHODS

A novel parameterized 3D cone trajectory design incorporating the 3D golden angle was integrated into 4D combined angiography and perfusion using radial imaging and arterial spin labeling (CAPRIA) to achieve higher spatial resolution and sampling efficiency for both dynamic angiography and perfusion imaging with flexible spatiotemporal resolution. Numerical simulations and physical phantom scanning were used to optimize the cone design. Eight healthy volunteers were scanned to compare the original radial trajectory in 4D CAPRIA with our newly designed cone trajectory. A locally low rank reconstruction method was used to leverage the complementary k-space sampling across time.

RESULTS

The improved sampling in the periphery of k-space obtained with the optimized 3D cone trajectory resulted in improved spatial resolution compared with the radial trajectory in phantom scans. Improved vessel sharpness and perfusion visualization were also achieved in vivo. Less dephasing was observed in the angiograms because of the short TE of our cone trajectory and the improved k-space sampling efficiency also resulted in higher repeatability compared to the original radial approach.

CONCLUSION

The proposed 3D cone trajectory combined with 3D golden angle ordering resulted in improved spatial resolution and image quality for both angiography and perfusion imaging and could potentially benefit other applications that require an efficient sampling scheme with flexible spatial and temporal resolution.

T1-weighted Motion Mitigation in Abdominal MRI: Technical Principles, Clinical Applications, Current Limitations, and Future Prospects

T1-weighted (T1W) pulse sequences are an indispensable component of clinical protocols in abdominal MRI but usually require multiple breath holds (BHs) during the examination, which not all patients can sustain. Patient motion can affect the quality of T1W imaging so that key diagnostic information, such as intrinsic signal intensity and contrast enhancement image patterns, cannot be determined. Patient motion also has a negative impact on examination efficiency, as multiple acquisition attempts prolong the duration of the examination and often remain noncontributory. Techniques for mitigation of motion-related artifacts at T1W imaging include multiple arterial acquisitions within one BH; free breathing with respiratory gating or respiratory triggering; and radial imaging acquisition techniques, such as golden-angle radial k-space acquisition (stack-of-stars). While each of these techniques has inherent strengths and limitations, the selection of a specific motion-mitigation technique is based on several factors, including the clinical task under investigation, downstream technical ramifications, patient condition, and user preference. The authors review the technical principles of free-breathing motion mitigation techniques in abdominal MRI with T1W sequences, offer an overview of the established clinical applications, and outline the existing limitations of these techniques. In addition, practical guidance for abdominal MRI protocol strategies commonly encountered in clinical scenarios involving patients with limited BH abilities is rendered. Future prospects of free-breathing T1W imaging in abdominal MRI are also discussed. ©RSNA, 2024 See the invited commentary by Fraum and An in this issue.

Artificial intelligence for neuro MRI acquisition: a review

OBJECT

To review recent advances of artificial intelligence (AI) in enhancing the efficiency and throughput of the MRI acquisition workflow in neuroimaging, including planning, sequence design, and correction of acquisition artifacts.

MATERIALS AND METHODS

A comprehensive analysis was conducted on recent AI-based methods in neuro MRI acquisition. The study focused on key technological advances, their impact on clinical practice, and potential risks associated with these methods.

RESULTS

The findings indicate that AI-based algorithms have a substantial positive impact on the MRI acquisition process, improving both efficiency and throughput. Specific algorithms were identified as particularly effective in optimizing acquisition steps, with reported improvements in workflow efficiency.

DISCUSSION

The review highlights the transformative potential of AI in neuro MRI acquisition, emphasizing the technological advances and clinical benefits. However, it also discusses potential risks and challenges, suggesting areas for future research to mitigate these concerns and further enhance AI integration in MRI acquisition.

T1 mapping of myocardium in rats using self-gated golden-angle acquisition

PURPOSE

The aim of this study is to design a method of myocardial T1 quantification in small laboratory animals and to investigate the effects of spatiotemporal regularization and the needed acquisition duration.

METHODS

We propose a compressed-sensing approach to T1 quantification based on self-gated inversion-recovery radial two/three-dimensional (2D/3D) golden-angle stack-of-stars acquisition with image reconstruction performed using total-variation spatiotemporal regularization. The method was tested on a phantom and on a healthy rat, as well as on rats in a small myocardium-remodeling study.

RESULTS

The results showed a good match of the T1 estimates with the results obtained using the ground-truth method on a phantom and with the literature values for rats myocardium. The proposed 2D and 3D methods showed significant differences between normal and remodeling myocardium groups for acquisition lengths down to approximately 5 and 15 min, respectively.

CONCLUSIONS

A new 2D and 3D method for quantification of myocardial T1 in rats was proposed. We have shown the capability of both techniques to distinguish between normal and remodeling myocardial tissue. We have shown the effects of image-reconstruction regularization weights and acquisition length on the T1 estimates.

Stochastic optimization of three-dimensional non-Cartesian sampling trajectory

PURPOSE

Optimizing three-dimensional (3D) k-space sampling trajectories is important for efficient MRI yet presents a challenging computational problem. This work proposes a generalized framework for optimizing 3D non-Cartesian sampling patterns via data-driven optimization.

METHODS

We built a differentiable simulation model to enable gradient-based methods for sampling trajectory optimization. The algorithm can simultaneously optimize multiple properties of sampling patterns, including image quality, hardware constraints (maximum slew rate and gradient strength), reduced peripheral nerve stimulation (PNS), and parameter-weighted contrast. The proposed method can either optimize the gradient waveform (spline-based freeform optimization) or optimize properties of given sampling trajectories (such as the rotation angle of radial trajectories). Notably, the method can optimize sampling trajectories synergistically with either model-based or learning-based reconstruction methods. We proposed several strategies to alleviate the severe nonconvexity and huge computation demand posed by the large scale. The corresponding code is available as an open-source toolbox.

RESULTS

We applied the optimized trajectory to multiple applications including structural and functional imaging. In the simulation studies, the image quality of a 3D kooshball trajectory was improved from 0.29 to 0.22 (NRMSE) with Stochastic optimization framework for 3D NOn-Cartesian samPling trajectorY (SNOPY) optimization. In the prospective studies, by optimizing the rotation angles of a stack-of-stars (SOS) trajectory, SNOPY reduced the NRMSE of reconstructed images from 1.19 to 0.97 compared to the best empirical method (RSOS-GR). Optimizing the gradient waveform of a rotational EPI trajectory improved participants' rating of the PNS from "strong" to "mild."

CONCLUSION

SNOPY provides an efficient data-driven and optimization-based method to tailor non-Cartesian sampling trajectories.

Kz-accelerated variable-density stack-of-stars MRI

This work aimed to develop a modified stack-of-stars golden-angle radial sampling scheme with variable-density acceleration along the slice (kz) dimension (referred to as VD-stack-of-stars) and to test this new sampling trajectory with multi-coil compressed sensing reconstruction for rapid motion-robust 3D liver MRI. VD-stack-of-stars sampling implements additional variable-density undersampling along the kz dimension, so that slice resolution (or volumetric coverage) can be increased without prolonging scan time. The new sampling trajectory (with increased slice resolution) was compared with standard stack-of-stars sampling with fully sampled kz (with standard slice resolution) in both non-contrast-enhanced free-breathing liver MRI and dynamic contrast-enhanced MRI (DCE-MRI) of the liver in volunteers. For both sampling trajectories, respiratory motion was extracted from the acquired radial data, and images were reconstructed using motion-compensated (respiratory-resolved or respiratory-weighted) dynamic radial compressed sensing reconstruction techniques. Qualitative image quality assessment (visual assessment by experienced radiologists) and quantitative analysis (as a metric of image sharpness) were performed to compare images acquired using the new and standard stack-of-stars sampling trajectories. Compared to standard stack-of-stars sampling, both non-contrast-enhanced and DCE liver MR images acquired with VD-stack-of-stars sampling presented improved overall image quality, sharper liver edges and increased hepatic vessel clarity in all image planes. The results have suggested that the proposed VD-stack-of-stars sampling scheme can achieve improved performance (increased slice resolution or volumetric coverage with better image quality) over standard stack-of-stars sampling in free-breathing DCE-MRI without increasing scan time. The reformatted coronal and sagittal images with better slice resolution may provide added clinical value.

Real-time MRI motion estimation through an unsupervised k-space-driven deformable registration network (KS-RegNet)

Purpose. Real-time three-dimensional (3D) magnetic resonance (MR) imaging is challenging because of slow MR signal acquisition, leading to highly under-sampled k-space data. Here, we proposed a deep learning-based, k-space-driven deformable registration network (KS-RegNet) for real-time 3D MR imaging. By incorporating prior information, KS-RegNet performs a deformable image registration between a fully-sampled prior image and on-board images acquired from highly-under-sampled k-space data, to generate high-quality on-board images for real-time motion tracking. Methods. KS-RegNet is an end-to-end, unsupervised network consisting of an input data generation block, a subsequent U-Net core block, and following operations to compute data fidelity and regularization losses. The input data involved a fully-sampled, complex-valued prior image, and the k-space data of an on-board, real-time MR image (MRI). From the k-space data, under-sampled real-time MRI was reconstructed by the data generation block to input into the U-Net core. In addition, to train the U-Net core to learn the under-sampling artifacts, the k-space data of the prior image was intentionally under-sampled using the same readout trajectory as the real-time MRI, and reconstructed to serve an additional input. The U-Net core predicted a deformation vector field that deforms the prior MRI to on-board real-time MRI. To avoid adverse effects of quantifying image similarity on the artifacts-ridden images, the data fidelity loss of deformation was evaluated directly in k-space. Results. Compared with Elastix and other deep learning network architectures, KS-RegNet demonstrated better and more stable performance. The average (±s.d.) DICE coefficients of KS-RegNet on a cardiac dataset for the 5- , 9- , and 13-spoke k-space acquisitions were 0.884 ± 0.025, 0.889 ± 0.024, and 0.894 ± 0.022, respectively; and the corresponding average (±s.d.) center-of-mass errors (COMEs) were 1.21 ± 1.09, 1.29 ± 1.22, and 1.01 ± 0.86 mm, respectively. KS-RegNet also provided the best performance on an abdominal dataset. Conclusion. KS-RegNet allows real-time MRI generation with sub-second latency. It enables potential real-time MR-guided soft tissue tracking, tumor localization, and radiotherapy plan adaptation.

Accelerating multi-echo chemical shift encoded water-fat MRI using model-guided deep learning

PURPOSE

To accelerate chemical shift encoded (CSE) water-fat imaging by applying a model-guided deep learning water-fat separation (MGDL-WF) framework to the undersampled k-space data.

METHODS

A model-guided deep learning water-fat separation framework is proposed for the acceleration using Cartesian/radial undersampling data. The proposed MGDL-WF combines the power of CSE water-fat imaging model and data-driven deep learning by jointly using a multi-peak fat model and a modified residual U-net network. The model is used to guide the image reconstruction, and the network is used to capture the artifacts induced by the undersampling. A data consistency layer is used in MGDL-WF to ensure the output images to be consistent with the k-space measurements. A Gauss-Newton iteration algorithm is adapted for the gradient updating of the networks.

RESULTS

Compared with the compressed sensing water-fat separation (CS-WF) algorithm/2-step procedure algorithm, the MGDL-WF increased peak signal-to-noise ratio (PSNR) by 5.31/5.23, 6.11/4.54, and 4.75 dB/1.88 dB with Cartesian sampling, and by 4.13/6.53, 2.90/4.68, and 1.68 dB/3.48 dB with radial sampling, at acceleration rates (R) of 4, 6, and 8, respectively. By using MGDL-WF, radial sampling increased the PSNR by 2.07 dB at R = 8, compared with Cartesian sampling.

CONCLUSIONS

The proposed MGDL-WF enables exploiting features of the water images and fat images from the undersampled multi-echo data, leading to improved performance in the accelerated CSE water-fat imaging. By using MGDL-WF, radial sampling can further improve the image quality with comparable scan time in comparison with Cartesian sampling.

Golden-Angle Radial MRI: Basics, Advances, and Applications

In recent years, golden‐angle radial sampling has received substantial attention and interest in the magnetic resonance imaging (MRI) community, and it has become a popular sampling trajectory for both research and clinical use. However, although the number of relevant techniques and publications has grown rapidly, there is still a lack of a review paper that provides a comprehensive overview and summary of the basics of golden‐angle rotation, the advantages and challenges/limitations of golden‐angle radial sampling, and recommendations in using different types of golden‐angle radial trajectories for MRI applications. Such a review paper is expected to be helpful both for clinicians who are interested in learning the potential benefits of golden‐angle radial sampling and for MRI physicists who are interested in exploring this research direction. The main purpose of this review paper is thus to present an overview and summary about golden‐angle radial MRI sampling. The review consists of three sections. The first section aims to answer basic questions such as: what is a golden angle; how is the golden angle calculated; why is golden‐angle radial sampling useful, and what are its limitations. The second section aims to review more advanced trajectories of golden‐angle radial sampling, including tiny golden‐angle rotation, stack‐of‐stars golden‐angle radial sampling, and three‐dimensional (3D) kooshball golden‐angle radial sampling. Their respective advantages and limitations and potential solutions to address these limitations are also discussed. Finally, the third section reviews MRI applications that can benefit from golden‐angle radial sampling and provides recommendations to readers who are interested in implementing golden‐angle radial trajectories in their MRI studies.

Radial sequences and compressed sensing in pediatric body magnetic resonance imaging

Magnetic resonance imaging (MRI) is often an ideal imaging modality for children of any age for any anatomy and for many pathologies. MRI sequences can be prescribed to produce high-resolution images of anatomical structures, characterize tissue composition, and detect physiological states and organ function. Shortening imaging sequences in any manner possible has been a topic of research and development in MRI since its emergence. Selection of imaging sequence parameters influences more than just the appearance and signal qualities of the imaged tissues; these details along with spatial encoding and data readout steps determine the time it takes to acquire an image. As each piece of image data is acquired and encoded with spatial and temporal information it is stored in k-space. As k-space is filled, either completely or partially, a diagnostic image or physiological data can be reconstructed. Shortening the length of time required for the readout step by efficiently filling k-space using compressed sensing and radial techniques is the subject of this manuscript.

Three dimensional radial echo planar imaging for functional MRI

PURPOSE

To demonstrate a novel 3D radial echo planar imaging (3D REPI) sequence for flexible, rapid, and motion-robust sampling in fMRI.

METHODS

The 3D REPI method expands on the recently described golden angle rotated EPI trajectory using radial batched internal navigator echoes (TURBINE) approach by exploiting the unused perpendicular direction in the EPI readout to form fast analogues of rotated stack of stars or spirals trajectories that cover all 3 dimensions of k-space. An iterative conjugate gradient algorithm with SENSE reconstruction and time-segmented non-uniform fast Fourier transform (FFT) was used for parallel imaging acceleration and to account for the effects of B₀ inhomogeneity. The golden angle rotation allowed for sliding window reconstruction schemes to be applied in brain BOLD fMRI experiments.

RESULTS

Combined whole brain visual and motor fMRI experiments were successfully carried out on a clinical 3T scanner at 2 mm isotropic and 1 × 1 × 2 mm³ resolutions using the 3D REPI design. Improved sampling characteristics and image quality were observed for twisted trajectories at the expense of prolonged readout times and off-resonance effects. The ability to correct for rigid motion correction was also demonstrated.

CONCLUSIONS

3D REPI presents a flexible approach for segmented volumetric fMRI with motion correction and high in-plane spatial resolutions. Improved BOLD fMRI brain activation maps were obtained using a sliding window reconstruction.

Dynamic MRI of swallowing: real-time volumetric imaging at 12 frames per second at 3 T

Objective

Dysphagia or difficulty in swallowing is a potentially hazardous clinical problem that needs regular monitoring. Real-time 2D MRI of swallowing is a promising radiation-free alternative to the current clinical standard: videofluoroscopy. However, aspiration may be missed if it occurs outside this single imaged slice. We therefore aimed to image swallowing in 3D real time at 12 frames per second (fps).

Materials and methods

At 3 T, three 3D real-time MRI acquisition approaches were compared to the 2D acquisition: an aligned stack-of-stars (SOS), and a rotated SOS with a golden-angle increment and with a tiny golden-angle increment. The optimal 3D acquisition was determined by computer simulations and phantom scans. Subsequently, five healthy volunteers were scanned and swallowing parameters were measured.

Results

Although the rotated SOS approaches resulted in better image quality in simulations, in practice, the aligned SOS performed best due to the limited number of slices. The four swallowing phases could be distinguished in 3D real-time MRI, even though the spatial blurring was stronger than in 2D. The swallowing parameters were similar between 2 and 3D.

Conclusion

At a spatial resolution of 2-by-2-by-6 mm with seven slices, swallowing can be imaged in 3D real time at a frame rate of 12 fps.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10334-021-00973-6.

A T1-based correction method for proton resonance frequency shift thermometry in breast tissue

PURPOSE

Develop and evaluate the effectiveness of a T1-based correction method for errors in proton resonant frequency shift thermometry due to non-local field effects caused by heating in fatty breast tissues.

METHODS

Computational models of human breast tissue were created by segmenting MRI data from a healthy human volunteer. MR-guided focused ultrasound (MRgFUS) heating and MR thermometry measurements were simulated in several locations in the heterogeneous segmented breast models. A T1-based correction method for PRF thermometry errors was applied and the maximum positive and negative errors and the root mean squared error (RMSE) in a region around each heating location was evaluated with and without correction. The method uses T1 measurements to estimate the temperature change in fatty tissues and correct for their influence. Experimental data from a heating study in cadaver breast tissue were analyzed, and the expected PRFS error computed.

RESULTS

The simulated MR thermometry had maximum single voxel errors ranging between 10% and 18% when no correction was applied. Applying the correction led to a considerable improvement, lowering the maximum error range to 2%-5%. The 5th to 95th percentile interval of the temperature error distribution was also lowered with correction, from approximately 3.5 to 1°C. This correction worked even when T1 times were uniformly raised or lowered by 5%-10%. The experimental data showed predicted errors of 15%.

CONCLUSIONS

This simulation study demonstrates that the T1-based correction method reduces MR thermometry errors due to non-local effects from heating in fatty tissues, potentially improving the accuracy of thermometry measurements during MRgFUS treatments. The presented correction method is reliant on having a patient-specific 3D model of the breast, and may be limited by the accuracy of the fat temperatures which in turn may be limited by noise or bias present in the T1 measurements.

Tiny golden angle stack-of-stars (tygaSoS) free-breathing functional lung imaging

PURPOSE

MRI of the lung parenchyma is still challenging due to cardiac and respiratory motion, and the low proton density and short T₂*. Clinical feasible MRI methods for functional lung assessment are of great interest. It was the objective of this study to evaluate the potential of combining the ultra-short echo-time stack-of-stars approach with tiny golden angle (tyGASoS) profile ordering for self-gated free-breathing lung imaging.

METHODS

Free-breathing tyGASoS data were acquired in 10 healthy volunteers (3 smoker (S), 7 non-smoker (NS)). Images in different respiratory phases were reconstructed applying an image-based self-gating technique. Resulting image quality and sharpness, and parenchyma visibility were qualitatively scored by three blinded independent reader, and the signal-to-noise ratio (SNR), proton fraction (f_P) and fractional ventilation (FV) quantified.

RESULT

The imaging protocol was well tolerated by all volunteers. Image quality was sufficient for subsequent quantitative analysis in all cases with good to excellent inter-reader reliability. Between expiration (EX) and inspiration (IN) significant differences (p < 0.001) were observed in SNR (EX: 3.73 ± 0.89, IN: 3.14 ± 0.74) and f_P (EX: 0.27 ± 0.09, IN: 0.25 ± 0.08). A significant (p < 0.05) higher f_P (EX/IN: 0.22 ± 0.07/0.21 ± 0.07 (NS), 0.33 ± 0.07/0.30 ± 0.06 (S)) was observed in the smoker group. No significant FV differences resulted between S and NS.

CONCLUSION

The study proves the feasibility of free-breathing tyGASoS for multiphase lung imaging. Changes in f_P may indicate an initial response in the smoker group and as such proves the sensitivity of the proposed technique. A major limitation in FV quantification rises from the large inter-subject variability of breathing patterns and amplitudes, requiring further consideration.

Model-based reconstruction for simultaneous multi-slice T1 mapping using single-shot inversion-recovery radial FLASH

PURPOSE

To develop a single-shot multi-slice T 1 mapping method by combing simultaneous multi-slice (SMS) excitations, single-shot inversion-recovery (IR) radial fast low-angle shot (FLASH), and a nonlinear model-based reconstruction method.

METHODS

SMS excitations are combined with a single-shot IR radial FLASH sequence for data acquisition. A previously developed single-slice calibrationless model-based reconstruction is extended to SMS, formulating the estimation of parameter maps and coil sensitivities from all slices as a single nonlinear inverse problem. Joint-sparsity constraints are further applied to the parameter maps to improve T 1 precision. Validations of the proposed method are performed for a phantom and for the human brain and liver in 6 healthy adult subjects.

RESULTS

Phantom results confirm good T 1 accuracy and precision of the simultaneously acquired multi-slice T 1 maps in comparison to single-slice references. In vivo human brain studies demonstrate the better performance of SMS acquisitions compared to the conventional spoke-interleaved multi-slice acquisition using model-based reconstruction. Aside from good accuracy and precision, the results of 6 healthy subjects in both brain and abdominal studies confirm good repeatability between scan and re-scans. The proposed method can simultaneously acquire T 1 maps for 5 slices of a human brain ( 0.75 × 0.75 × 5 mm 3 ) or 3 slices of the abdomen ( 1.25 × 1.25 × 6 mm 3 ) within 4 seconds.

CONCLUSIONS

The IR SMS radial FLASH acquisition together with a nonlinear model-based reconstruction enable rapid high-resolution multi-slice T 1 mapping with good accuracy, precision, and repeatability.

Improved 3D real-time MRI of speech production

PURPOSE

To provide 3D real-time MRI of speech production with improved spatio-temporal sharpness using randomized, variable-density, stack-of-spiral sampling combined with a 3D spatio-temporally constrained reconstruction.

METHODS

We evaluated five candidate (k, t) sampling strategies using a previously proposed gradient-echo stack-of-spiral sequence and a 3D constrained reconstruction with spatial and temporal penalties. Regularization parameters were chosen by expert readers based on qualitative assessment. We experimentally determined the effect of spiral angle increment and k_z temporal order. The strategy yielding highest image quality was chosen as the proposed method. We evaluated the proposed and original 3D real-time MRI methods in 2 healthy subjects performing speech production tasks that invoke rapid movements of articulators seen in multiple planes, using interleaved 2D real-time MRI as the reference. We quantitatively evaluated tongue boundary sharpness in three locations at two speech rates.

RESULTS

The proposed data-sampling scheme uses a golden-angle spiral increment in the k_x -k_y plane and variable-density, randomized encoding along k_z . It provided a statistically significant improvement in tongue boundary sharpness score (P < .001) in the blade, body, and root of the tongue during normal and 1.5-times speeded speech. Qualitative improvements were substantial during natural speech tasks of alternating high, low tongue postures during vowels. The proposed method was also able to capture complex tongue shapes during fast alveolar consonant segments. Furthermore, the proposed scheme allows flexible retrospective selection of temporal resolution.

CONCLUSION

We have demonstrated improved 3D real-time MRI of speech production using randomized, variable-density, stack-of-spiral sampling with a 3D spatio-temporally constrained reconstruction.

Compressed sensing with signal averaging for improved sensitivity and motion artifact reduction in fluorine-19 MRI

Fluorine-19 (19 F) MRI of injected perfluorocarbon emulsions (PFCs) allows for the non-invasive quantification of inflammation and cell tracking, but suffers from a low signal-to-noise ratio and extended scan time. To address this limitation, we tested the hypotheses that a 19 F MRI pulse sequence that combines a specific undersampling regime with signal averaging has both increased sensitivity and robustness against motion artifacts compared with a non-averaged fully sampled pulse sequence, when both datasets are reconstructed with compressed sensing. As a proof of principle, numerical simulations and phantom experiments were performed on selected variable ranges to characterize the point spread function of undersampling patterns, as well as the vulnerability to noise of undersampling and reconstruction parameters with paired numbers of x signal averages and acceleration factor x (NAx-AFx). The numerical simulations demonstrated that a probability density function that uses 25% of the samples to fully sample the k-space central area allowed for an optimal balance between limited blurring and artifact incoherence. At all investigated noise levels, the Dice similarity coefficient (DSC) strongly depended on the regularization parameters and acceleration factor. In phantoms, the motion robustness of an NA8-AF8 undersampling pattern versus NA1-AF1 was evaluated with simulated and real motion patterns. Differences were assessed with the DSC, which was consistently higher for the NA8-AF8 compared with the NA1-AF1 strategy, for both simulated and real cyclic motion patterns (P < 0.001). Both strategies were validated in vivo in mice (n = 2) injected with perfluoropolyether. Here, the images displayed a sharper delineation of the liver with the NA8-AF8 strategy than with the NA1-AF1 strategy. In conclusion, we validated the hypotheses that in 19 F MRI the combination of undersampling and averaging improves both the sensitivity and the robustness against motion artifacts.

Semiautomatic cerebrovascular territory mapping based on dynamic ASL MR angiography without vessel-encoded labeling

PURPOSE

Characterizing vessel territories can provide crucial information for evaluation of cerebrovascular disorders. In this study, we present a novel postprocessing pipeline for vascular territorial imaging of cerebral arteries based on a noncontrast enhanced time-resolved 4D magnetic resonance angiography (MRA).

METHODS

Eight healthy participants, 1 Moyamoya patient, and 1 arteriovenous malformations patient were recruited. Territorial segmentation and relative blood flow rate calculations of cerebral arteries including left and right middle cerebral arteries and left and right posterior cerebral arteries were carried out based on the 4D MRA-derived arterial arrival time maps of intracranial vessels.

RESULTS

Among healthy young subjects, the average relative blood flow rate values corresponding to left and right middle cerebral arteries and left and right posterior cerebral arteries were 35.9 ± 5.9%, 32.9 ± 7.5%, 15.4 ± 3.8%, and 15.9 ± 2.5%, respectively. Excellent agreement was observed between relative blood flow rate values obtained from the proposed 4D MRA-based method and reference 2D phase contrast MRI. Abnormal cerebral circulations were visualized and quantified on both patients using the developed technique.

CONCLUSION

The vascular territorial imaging technique developed in this study allowed for the generation of territorial maps with user-defined level of details within a clinically feasible scan time, and as such may provide useful information to assess cerebral circulation balance in different pathologies.

Cardiac and Respiratory Self-Gating in Radial MRI Using an Adapted Singular Spectrum Analysis (SSA-FARY)

Cardiac Magnetic Resonance Imaging (MRI) is time-consuming and error-prone. To ease the patient's burden and to increase the efficiency and robustness of cardiac exams, interest in methods based on continuous steady-state acquisition and self-gating has been growing in recent years. Self-gating methods extract the cardiac and respiratory signals from the measurement data and then retrospectively sort the data into cardiac and respiratory phases. Repeated breathholds and synchronization with the heart beat using some external device as required in conventional MRI are then not necessary. In this work, we introduce a novel self-gating method for radially acquired data based on a dimensionality reduction technique for time-series analysis (SSA-FARY). Building on Singular Spectrum Analysis, a zero-padded, time-delayed embedding of the auto-calibration data is analyzed using Principle Component Analysis. We demonstrate the basic functionality of SSA-FARY using numerical simulations and apply it to in-vivo cardiac radial single-slice bSSFP and Simultaneous Multi-Slice radiofrequency-spoiled gradient-echo measurements, as well as to Stack-of-Stars bSSFP measurements. SSA-FARY reliably detects the cardiac and respiratory motion and separates it from noise. We utilize the generated signals for high-dimensional image reconstruction using parallel imaging and compressed sensing with in-plane wavelet and (spatio-)temporal total-variation regularization.

Towards accelerated quantitative sodium MRI at 7 T in the skeletal muscle: Comparison of anisotropic acquisition- and compressed sensing techniques

PURPOSE

To compare three anisotropic acquisition schemes and three compressed sensing (CS) approaches for accelerated tissue sodium concentration (TSC) quantification using ²³Na MRI at 7 T.

MATERIALS AND METHODS

Three anisotropic 3D-radial acquisition sequences were evaluated using simulations, phantom- and in vivo TSC measurements: An anisotropic density-adapted 3D-radial sequence (3DPR-C), a 3D acquisition-weighted density-adapted stack-of-stars sampling scheme (SOS) and a SOS approach with golden-ratio rotation (SOS-GR). Eight healthy volunteers were examined at a 7 Tesla MRI system. TSC measurements of the calf were conducted with a nominal spatial resolution of Δx = (3.0 × 3.0 × 15.0) mm³ and a field of view of (156.0 × 156.0 × 240.0) mm³ for multiple undersampling factors (USF). Three CS reconstructions were evaluated: Total variation CS (TV-CS), 3D dictionary-learning compressed sensing (3D-DLCS) and TV-CS with a block matching prior (TV-BL-CS). Results of the simulations and measurements were compared to a simulated ground truth (GT) or a fully sampled reference measurement (FS), respectively. The deviation of the mean TSC evaluated in multiple ROI (mE_GT/FS) and the normalized root-mean-squared error (NRMSE) for simulations were evaluated for CS and NUFFT reconstructions.

RESULTS

In simulations, the SOS-GR yielded the lowest NRMSE and mE_GT (< 4%) with NUFFT for an acquisition time (TA) of less than 2 min. CS further improved the results. In simulations and measurements, the best TSC quantification results were obtained with 3D-DLCS and SOS-GR (lowest NRMSE, mE_GT < 2.6% in simulations, mE_GT < 10.7% for phantom measurements and mE_FS < 6% in vivo) with an USF = 4.1 (TA < 2 min). TV-CS showed no or only slight improvements to NUFFT. The results of TV-BL-CS were similar to 3D-DLCS.

DISCUSSION

The TA for TSC measurements could be reduced to less than 2 min by using adapted sequences such as SOS-GR and CS reconstruction approaches such as 3D-DLCS or TV-BL-CS, while the quantitative accuracy stays comparable to a fully sampled NUFFT reconstruction (approx. 8 min TA). In future, the lower TA could improve clinical applicability of TSC measurements.

Rotated spiral RARE for high spatial and temporal resolution volumetric arterial spin labeling acquisition

Background:

Arterial Spin Labeling (ASL) MRI can provide quantitative images that are sensitive to both time averaged blood flow and its temporal fluctuations. 3D image acquisitions for ASL are desirable because they are more readily compatible with background suppression to reduce noise, can reduce signal loss and distortion, and provide uniform flow sensitivity across the brain. However, single-shot 3D acquisition for maximal temporal resolution typically involves degradation of image quality through blurring or noise amplification by parallel imaging. Here, we report a new approach to accelerate a common stack of spirals 3D image acquisition by pseudo golden-angle rotation and compressed sensing reconstruction without any degradation of time averaged blood flow images.

Methods:

28 healthy volunteers were imaged at 3T with background-suppressed unbalanced pseudo-continuous ASL combined with a pseudo golden-angle Stack-of-Spirals 3D RARE readout. A fully-sampled perfusion-weighted volume was reconstructed by 3D non-uniform Fast Fourier Transform (nuFFT) followed by sum-of-squares combination of the 32 individual channels. Coil sensitivities were estimated followed by reconstruction of the 39 single-shot volumes using an L₁-wavelet Compressed-Sensing reconstruction. Finally, brain connectivity analyses were performed in regions where BOLD signal suffers from low signal-to-noise ratio and susceptibility artifacts.

Results:

Image quality, assessed with a non-reference 3D blurring metric, of full time averaged blood flow was comparable to a conventional interleaved acquisition. The temporal resolution provided by the acceleration enabled identification and quantification of resting-state networks even in inferior regions such as the amygdala and inferior frontal lobes, where susceptibility artifacts can degrade conventional resting-state fMRI acquisitions.

Conclusion:

This approach can provide measures of blood flow modulations and resting-state networks for free within any research or clinical protocol employing ASL for resting blood flow.

Whole-heart, ungated, free-breathing, cardiac-phase-resolved myocardial perfusion MRI by using Continuous Radial Interleaved simultaneous Multi-slice acquisitions at sPoiled steady-state (CRIMP)

PURPOSE

To develop a whole-heart, free-breathing, non-electrocardiograph (ECG)-gated, cardiac-phase-resolved myocardial perfusion MRI framework (CRIMP; Continuous Radial Interleaved simultaneous Multi-slice acquisitions at sPoiled steady-state) and test its quantification feasibility.

METHODS

CRIMP used interleaved radial simultaneous multi-slice (SMS) slice groups to cover the whole heart in 9 or 12 short-axis slices. The sequence continuously acquired data without magnetization preparation, ECG gating or breath-holding, and captured multiple cardiac phases. Images were reconstructed by a motion-compensated patch-based locally low-rank reconstruction. Bloch simulations were performed to study the signal-to-noise ratio/contrast-to-noise ratio (SNR/CNR) for CRIMP and to study the steady-state signal under motion. Seven patients were scanned with CRIMP at stress and rest to develop the sequence. One human and two dogs were scanned at rest with a dual-bolus method to test the quantification feasibility of CRIMP. The dual-bolus scans were performed using both CRIMP and an ungated radial SMS saturation recovery (SMS-SR) sequence with injection dose = 0.075 mmol/kg to compare the sequences in terms of SNR, cardiac phase resolution and quantitative myocardial blood flow (MBF).

RESULTS

Perfusion images with multiple cardiac phases in all image slices with a temporal resolution of 72 ms/frame were obtained. Simulations and in-vivo acquisitions showed CRIMP kept the inner slices in steady-state regardless of motion. CRIMP outperformed SMS-SR in slice coverage (9 over 6), SNR (mean 20% improvement), and provided cardiac phase resolution. CRIMP and SMS-SR sequences provided comparable MBF values (rest systolic CRIMP = 0.58 ± 0.07, SMS-SR = 0.61 ± 0.16).

CONCLUSION

CRIMP allows for whole-heart, cardiac-phase-resolved myocardial perfusion images without ECG-gating or breath-holding. The sequence can provide MBF if an accurate arterial input function is obtained separately.

A variable flip angle golden-angle-ordered 3D stack-of-radial MRI technique for simultaneous proton resonant frequency shift and T1 -based thermometry

PURPOSE

To develop and evaluate a variable-flip-angle golden-angle-ordered 3D stack-of-radial MRI technique for simultaneous proton resonance frequency shift (PRF) and T₁ -based thermometry in aqueous and adipose tissues, respectively.

METHODS

The proposed technique acquires multiecho radial k-space data in segments with alternating flip angles to measure 3D temperature maps dynamically on the basis of PRF and T₁ . A sliding-window k-space weighted image contrast filter is used to increase temporal resolution. PRF is measured in aqueous tissues and T₁ in adipose tissues using fat/water masks. The accuracy for T₁ quantification was evaluated in a reference T₁ /T₂ phantom. In vivo nonheating experiments were conducted in healthy subjects to evaluate the stability of PRF and T₁ in the brain, prostate, and breast. The proposed technique was used to monitor high-intensity focused ultrasound (HIFU) ablation in ex vivo porcine fat/muscle tissues and compared to temperature probe readings.

RESULTS

The proposed technique achieved 3D coverage with 1.1-mm to 1.3-mm in-plane resolution and 2-s to 5-s temporal resolution. During 20 to 30 min of nonheating in vivo scans, the temporal coefficient of variation for T₁ was <5% in the brain, prostate, and breast fatty tissues, while the standard deviation of relative PRF temperature change was within 3°C in aqueous tissues. During ex vivo HIFU ablation, the temperatures measured by PRF and T₁ were consistent with temperature probe readings, with an absolute mean difference within 2°C.

CONCLUSION

The proposed technique achieves simultaneous PRF and T₁ -based dynamic 3D MR temperature mapping in aqueous and adipose tissues. It may be used to improve MRI-guided thermal procedures.

A GRAPPA algorithm for arbitrary 2D/3D non-Cartesian sampling trajectories with rapid calibration

PURPOSE

GRAPPA is a popular reconstruction method for Cartesian parallel imaging, but is not easily extended to non-Cartesian sampling. We introduce a general and practical GRAPPA algorithm for arbitrary non-Cartesian imaging.

METHODS

We formulate a general GRAPPA reconstruction by associating a unique kernel with each unsampled k-space location with a distinct constellation, that is, local sampling pattern. We calibrate these generalized kernels using the Fourier transform phase shift property applied to fully gridded or separately acquired Cartesian Autocalibration signal (ACS) data. To handle the resulting large number of different kernels, we introduce a fast calibration algorithm based on nonuniform FFT (NUFFT) and adoption of circulant ACS boundary conditions. We applied our method to retrospectively under-sampled rotated stack-of-stars/spirals in vivo datasets, and to a prospectively under-sampled rotated stack-of-spirals functional MRI acquisition with a finger-tapping task.

RESULTS

We reconstructed all datasets without performing any trajectory-specific manual adaptation of the method. For the retrospectively under-sampled experiments, our method achieved image quality (i.e., error and g-factor maps) comparable to conjugate gradient SENSE (cg-SENSE) and SPIRiT. Functional activation maps obtained from our method were in good agreement with those obtained using cg-SENSE, but required a shorter total reconstruction time (for the whole time-series): 3 minutes (proposed) vs 15 minutes (cg-SENSE).

CONCLUSIONS

This paper introduces a general 3D non-Cartesian GRAPPA that is fast enough for practical use on today's computers. It is a direct generalization of original GRAPPA to non-Cartesian scenarios. The method should be particularly useful in dynamic imaging where a large number of frames are reconstructed from a single set of ACS data.

Simultaneous proton resonance frequency shift thermometry and T1 measurements using a single reference variable flip angle T1 method

PURPOSE

Implement simultaneous proton resonance frequency (PRF) shift and T₁ measurements with equivalent temporal resolution using a single reference variable flip angle method. This novel method allows for simultaneous thermometry in both aqueous and fatty tissue.

METHODS

This method acquires a single reference image at the lower flip angle and all dynamic images at the higher angle. T₁ is calculated using a single reference variable flip angle method, which accounts for the reference image temperature remaining constant. Monte Carlo simulations determined the optimal dynamic flip angle for combined PRF and T₁ measurements. This method was evaluated in MR-guided focused ultrasound heating experiments using a gelatin phantom and human cadaver breasts. In vivo measurement precision was demonstrated in healthy female volunteers under nonheating conditions.

RESULTS

Temperature rise during MR-guided focused ultrasound heating was measured in aqueous tissue with both PRF and T₁ . Both measures show good qualitative agreement in both space and time in aqueous tissue. The T₁ change due to temperature increase was measured in fat, demonstrating the expected temporal response. The dynamic flip angle that produces optimal SNR for PRF measurements is lower than the optimal angle for T₁ measurements, necessitating the selection of a compromise angle.

CONCLUSION

The single reference variable flip angle method provides a reliable way to simultaneously measure PRF temperature and T₁ change and overcomes PRF's inability to simultaneously monitor temperature in aqueous and adipose tissues. Future work will calibrate T₁ change to temperature, enabling real-time temperature in fat and increasing patient safety and treatment efficacy during thermal interventional treatments.

Frequency-modulated SSFP with radial sampling and subspace reconstruction: A time-efficient alternative to phase-cycled bSSFP

PURPOSE

A novel subspace-based reconstruction method for frequency-modulated balanced steady-state free precession (fmSSFP) MRI is presented. In this work, suitable data acquisition schemes, subspace sizes, and efficiencies for banding removal are investigated.

THEORY AND METHODS

By combining a fmSSFP MRI sequence with a 3D stack-of-stars trajectory, scan efficiency is maximized as spectral information is obtained without intermediate preparation phases. A memory-efficient reconstruction routine is implemented by introducing the low-frequency Fourier transform as a subspace which allows for the formulation of a convex reconstruction problem. The removal of banding artifacts is investigated by comparing the proposed acquisition and reconstruction technique to phase-cycled bSSFP MRI. Aliasing properties of different undersampling schemes are analyzed and water/fat separation is demonstrated by reweighting the reconstructed subspace coefficients to generate virtual spectral responses in a post-processing step.

RESULTS

A simple root-of-sum-of-squares combination of the reconstructed subspace coefficients yields high-SNR images with the characteristic bSSFP contrast but without banding artifacts. Compared to Golden-Angle trajectories, turn-based sampling schemes were superior in minimizing aliasing across reconstructed subspace coefficients. Water/fat separated images of the human knee were obtained by reweighting subspace coefficients.

CONCLUSIONS

The novel subspace-based fmSSFP MRI technique emerges as a time-efficient alternative to phase-cycled bSFFP. The method does not need intermediate preparation phases, offers high SNR and avoids banding artifacts. Reweighting of the reconstructed subspace coefficients allows for generating virtual spectral responses with applications to water/fat separation.

Effect of k-space-weighted image contrast and ultrasound focus size on the accuracy of proton resonance frequency thermometry

PURPOSE

To construct a predictive model that describes how the duration and symmetry of a k-space-weighted image contrast (KWIC) window affects the temporal resolution of differently sized ultrasound foci when using a pseudo-golden angle stack-of-stars acquisition.

METHODS

We performed a modulation analysis of proton resonance frequency temperature measurements to create the temporal modulation transfer function for KWIC windows of different symmetry and temporal duration. We reconstructed simulated ultrasound heating trajectories and stack-of-stars k-space data as well as experimental phantom data using the same trajectories. Images were reconstructed using symmetric and asymmetric KWIC windows of 3 different temporal durations. Simulated results were compared against temporal modulation transfer function predictions, experimental results, and the original simulated temperatures.

RESULTS

The temporal modulation transfer function shows that temporal resolution with KWIC reconstructions depend on the object size. The KWIC window duration affected SNR and severity of undersampling artifacts. Accuracy and response delay improved as the KWIC window duration decreased or the size of the heated region within the KWIC plane increased. Precision worsened as the window duration decreased. Using a symmetric window eliminated the response delay to heated region size but introduced a large reconstruction delay.

CONCLUSION

The accuracy and precision of proton resonance frequency temperature measurements from a stack-of-stars acquisition using a sliding KWIC window reconstruction are dependent on the size of the KWIC window and the size and shape of the heated region. The temporal modulation transfer function of KWIC reconstructions for any object size can predict the temporal response to changes in signal being acquired, such as temperature and contrast enhancement.

The Influence of Radial Undersampling Schemes on Compressed Sensing in Cardiac DTI

Diffusion tensor imaging (DTI) is known to suffer from long acquisition time, which greatly limits its practical and clinical use. Undersampling of k-space data provides an effective way to reduce the amount of data to acquire while maintaining image quality. Radial undersampling is one of the most popular non-Cartesian k-space sampling schemes, since it has relatively lower sensitivity to motion than Cartesian trajectories, and artifacts from linear reconstruction are more noise-like. Therefore, radial imaging is a promising strategy of undersampling to accelerate acquisitions. The purpose of this study is to investigate various radial sampling schemes as well as reconstructions using compressed sensing (CS). In particular, we propose two randomly perturbed radial undersampling schemes: golden-angle and random angle. The proposed methods are compared with existing radial undersampling methods, including uniformity-angle, randomly perturbed uniformity-angle, golden-angle, and random angle. The results on both simulated and real human cardiac diffusion weighted (DW) images show that, for the same amount of k-space data, randomly sampling around a random radial line results in better reconstruction quality for DTI indices, such as fractional anisotropy (FA), mean diffusivities (MD), and that the randomly perturbed golden-angle undersampling yields the best results for cardiac CS-DTI image reconstruction.

Frequency-modulated bSSFP for phase-sensitive separation of water and fat

Our study proposes the use of a frequency-modulated acquisition which suppresses banding artefacts in combination with a phase-sensitive water-fat separation algorithm. The performance of the phase-sensitive separation for standard bSSFP, complex sum combination thereof, and frequency-modulated bSSFP were compared in in vivo measurements of the upper and lower legs at 1.5 and 3 T. It is shown, that the standard acquisition suffered from banding artefacts and major swaps between tissues. The dual-acquisition bSSFP could alleviate banding artefacts and only minor swaps occurred, but it comes at the expense of a doubled acquisition. In the frequency-modulated acquisitions all banding artefacts and the associated phase jumps were eliminated and no swaps between tissues occurred. It therefore provides a means to robustly separate water and fat, in one single radial bSSFP scan, using the phase-sensitive approach, even in the presence of high field inhomogeneities.

Accelerated noncontrast-enhanced 4-dimensional intracranial MR angiography using golden-angle stack-of-stars trajectory and compressed sensing with magnitude subtraction

PURPOSE

To evaluate the feasibility and performance of compressed sensing (CS) with magnitude subtraction regularization in accelerating non-contrast-enhanced dynamic intracranial MR angiography (NCE-dMRA).

METHODS

A CS algorithm was introduced in NCE-dMRA by exploiting the sparsity of the magnitude difference of the control and label images. The NCE-dMRA data were acquired using golden-angle stack-of-stars trajectory on six healthy volunteers and one patient with arteriovenous fistula. Images were reconstructed using (i) the proposed magnitude-subtraction CS (MS-CS); (ii) complex-subtraction CS; (iii) independent CS; and (iv) view-sharing with k-space weighted image contrast (KWIC). The dMRA image quality was compared across the four reconstruction strategies. The proposed MS-CS method was further compared with KWIC for temporal fidelity of depicting dynamic flow.

RESULTS

The proposed MS-CS method was able to reconstruct NCE-dMRA images with detailed vascular structures and clean background. It provided better subjective image quality than the other two CS strategies (P < 0.05). Compared with KWIC, MS-CS showed similar image quality, but reduced temporal blurring in delineating the fine distal arteries.

CONCLUSIONS

The MS-CS method is a promising CS technique for accelerating NCE-dMRA acquisition without compromising image quality and temporal fidelity. Magn Reson Med 79:867-878, 2018. © 2017 International Society for Magnetic Resonance in Medicine.