RACER-GRASP: Respiratory-weighted, aortic contrast enhancement-guided and coil-unstreaking golden-angle radial sparse MRI

High-resolution free-breathing hepatobiliary phase MRI of the liver using XD-GRASP

PURPOSE

To evaluate the performance of high-resolution free-breathing (FB) hepatobiliary phase imaging of the liver using the eXtra-Dimension Golden-angle RAdial Sparse Parallel (XD-GRASP) MRI technique.

METHODS

Fifty-eight clinical patients (41 males, mean age = 52.9 ± 12.9) with liver lesions who underwent dynamic contrast-enhanced MRI with a liver-specific contrast agent were prospectively recruited for this study. Both breath-hold volumetric interpolated examination (BH-VIBE) imaging and FB imaging were performed during the hepatobiliary phase. FB images were acquired using a stack-of-stars golden-angle radial sequence and were reconstructed using the XD-GRASP method. Two experienced radiologists blinded to acquisition schemes independently scored the overall image quality, liver edge sharpness, hepatic vessel clarity, conspicuity of lesion, and overall artifact level of each image. The non-parametric paired two-tailed Wilcoxon signed-rank test was used for statistical analysis.

RESULTS

Compared to BH-VIBE images, XD-GRASP images received significantly higher scores (P < 0.05) for the liver edge sharpness (4.83 ± 0.45 vs 4.29 ± 0.46), the hepatic vessel clarity (4.64 ± 0.67 vs 4.15 ± 0.56) and the conspicuity of lesion (4.75 ± 0.53 vs 4.31 ± 0.50). There were no significant differences (P > 0.05) between BH-VIBE and XD-GRASP images for the overall image quality (4.61 ± 0.50 vs 4.74 ± 0.47) and the overall artifact level (4.13 ± 0.44 vs 4.05 ± 0.61).

CONCLUSION

Compared to conventional BH-VIBE MRI, FB radial acquisition combined with XD-GRASP reconstruction facilitates higher spatial resolution imaging of the liver during the hepatobiliary phase. This enhancement can significantly improve the visualization and evaluation of the liver.

Review of oxygen-enhanced lung mri: Pulse sequences for image acquisition and T1 measurement

Oxygen-enhanced MR imaging (OE-MRI) is a special proton imaging technique that can be performed without modifying the scanner hardware. Many fundamental studies have been conducted following the initial reporting of this technique in 1996, illustrating the high potential for its clinical application. This review aims to summarise and analyse current pulse sequences and T1 measurement methods for OE-MRI, including fundamental theories, existing pulse sequences applied to OE-MRI acquisition and T1 mapping. Wash-in and wash-out time identify lung function and are sensitive to ventilation; thus, dynamic OE-MRI is also discussed in this review. We compare OE-MRI with the primary competitive technique, hyperpolarised gas MRI. Finally, an overview of lower-field applications of OE-MRI is highlighted, as relatively recent publications demonstrated positive results. Lower-field OE-MRI, which is lower than 1.5 T, could be an alternative modality for detecting lung diseases. This educational review is aimed at researchers who want a quick summary of the steps needed to perform pulmonary OE-MRI with a particular focus on sequence design, settings, and quantification methods.

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.

4D Golden-Angle Radial MRI at Subsecond Temporal Resolution

Intraframe motion blurring, as a major challenge in free-breathing dynamic MRI, can be reduced if high temporal resolution can be achieved. To address this challenge, this work proposes a highly accelerated 4D (3D + time) dynamic MRI framework with subsecond temporal resolution that does not require explicit motion compensation. The method combines standard stack-of-stars golden-angle radial sampling and tailored GRASP-Pro (Golden-angle RAdial Sparse Parallel imaging with imProved performance) reconstruction. Specifically, 4D dynamic MRI acquisition is performed continuously without motion gating or sorting. The k-space centers in stack-of-stars radial data are organized to guide estimation of a temporal basis, with which GRASP-Pro reconstruction is employed to enforce joint low-rank subspace and sparsity constraints. This new basis estimation strategy is the new feature proposed for subspace-based reconstruction in this work to achieve high temporal resolution (e.g., subsecond/3D volume). It does not require sequence modification to acquire additional navigation data, it is compatible with commercially available stack-of-stars sequences, and it does not need an intermediate reconstruction step. The proposed 4D dynamic MRI approach was tested in abdominal motion phantom, free-breathing abdominal MRI, and dynamic contrast-enhanced MRI (DCE-MRI). Our results have shown that GRASP-Pro reconstruction with the new basis estimation strategy enables highly-accelerated 4D dynamic imaging at subsecond temporal resolution (with five spokes or less for each dynamic frame per image slice) for both free-breathing non-DCE-MRI and DCE-MRI. In the abdominal phantom, better image quality with lower root mean square error and higher structural similarity index was achieved using GRASP-Pro compared with standard GRASP. With the ability to acquire each 3D image in less than 1 s, intraframe respiratory blurring can be intrinsically reduced for body applications with our approach, which eliminates the need for explicit motion detection and motion compensation.

Cancellation of streak artifacts in radial abdominal imaging using interference null space projection

PURPOSE

In radial abdominal imaging, it has been commonly observed that signal from the arms cause streaks due to system imperfections. We previously introduced a streak removal technique (B-STAR), which is inherently spatially variant and limited to work in image space. In this work, we propose a spatially invariant streak cancellation technique (CACTUS), which can be applied in either image space or k-space and is compatible with iterative reconstructions.

THEORY AND METHODS

Streak sources are typically spatially localized and can be represented using a low-dimensional subspace. CACTUS identifies the streak subspace by leveraging the spatial redundancy of receiver coils and projects the data onto the streak null space to eliminate the streaks. When applied in k-space, CACTUS can be combined with iterative reconstructions. CACTUS was tested in phantoms and in vivo abdominal imaging using a radial turbo spin-echo pulse sequence.

RESULTS

In phantoms, CACTUS improved T2 estimation in comparison to previous de-streaking methods. In vivo experiments showed that CACTUS reduced streaks and yielded T2 estimation, in regions affected by streaks, closer to a streak-free reference. Evaluation using a clinical abdominal dataset (n = 20) showed that CACTUS is comparable to B-STAR and yields significantly better signal preservation and streak cancellation than coil removal and suppression methods.

CONCLUSION

CACTUS provides superior signal preservation and streak reduction performance compared to coil removal and suppression methods. As a clear advantage over B-STAR, CACTUS can be integrated with iterative reconstruction methods. In abdominal T2 mapping, CACTUS improves the accuracy of parameter estimation in areas affected by streaks.

Repeatability and robustness of MP-GRASP T1 mapping

PURPOSE

To demonstrate the repeatability of fast 3D T1 mapping using Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI and its robustness to variation of imaging parameters including flip angle and spatial resolution in phantoms and the brain.

THEORY AND METHODS

Multiple imaging experiments were performed to (1) assess the robustness of MP-GRASP T1 mapping to B1 inhomogeneity using a single tube phantom filled with uniform MnCl2 liquid; (2) compare the repeatability of T1 mapping between MP-GRASP and inversion recovery-based spin-echo (IR-SE; over 12 scans), using a commercial T1MES phantom; (3) evaluate the longitudinal variation of T1 estimation using MP-GRASP with varying imaging parameters, including spatial resolution, flip angle, TR/TE, and acceleration rate, using the T1MES phantom (106 scans performed over a period of 12 months); and (4) evaluate the variation of T1 estimation using MP-GRASP with varying imaging parameters in the brain (24 scans in a single visit). In addition, the accuracy of MP-GRASP T1 mapping was also validated against IR-SE by performing linear correlation and calculating the Lin's concordance correlation coefficient (CCC).

RESULTS

MP-GRASP demonstrates good robustness to B1 inhomogeneity, with intra-slice variability below 1% in the single tube phantom experiment. The longitudinal variability is good both in the phantom (below 2.5%) and in the brain (below 2%) with varying imaging parameters. The T1 values estimated from MP-GRASP are accurate compared to that from the IR-SE imaging (R2 = 0.997, Lin's CCC = 0.996).

CONCLUSION

MP-GRASP shows excellent repeatability of T1 estimation over time, and it is also robust to variation of different imaging parameters evaluated in this study.

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.

Dynamic pulmonary MRI using motion-state weighted motion-compensation (MostMoCo) reconstruction with ultrashort TE: A structural and functional study

PURPOSE

To improve the quality of structural images and the quantification of ventilation in free-breathing dynamic pulmonary MRI.

METHODS

A 3D radial ultrashort TE (UTE) sequence with superior-inferior navigators was used to acquire pulmonary data during free breathing. All acquired data were binned into different motion states according to the respiratory signal extracted from superior-inferior navigators. Motion-resolved images were reconstructed using eXtra-Dimensional (XD) UTE reconstruction. The initial motion fields were generated by registering images at each motion state to other motion states in motion-resolved images. A motion-state weighted motion-compensation (MostMoCo) reconstruction algorithm was proposed to reconstruct the dynamic UTE images. This technique, termed as MostMoCo-UTE, was compared with XD-UTE and iterative motion-compensation (iMoCo) on a porcine lung and 10 subjects.

RESULTS

MostMoCo reconstruction provides higher peak SNR (37.0 vs. 35.4 and 34.2) and structural similarity (0.964 vs. 0.931 and 0.947) compared to XD-UTE and iMoCo in the porcine lung experiment. Higher apparent SNR and contrast-to-noise ratio are achieved using MostMoCo in the human experiment. MostMoCo reconstruction better preserves the temporal variations of signal intensity of parenchyma compared to iMoCo, shows reduced random noise and improved sharpness of anatomical structures compared to XD-UTE. In the porcine lung experiment, the quantification of ventilation using MostMoCo images is more accurate than that using XD-UTE and iMoCo images.

CONCLUSION

The proposed MostMoCo-UTE provides improved quality of structural images and quantification of ventilation for free-breathing pulmonary MRI. It has the potential for the detection of structural and functional disorders of the lung in clinical settings.

Revisiting DCE-MRI: Classification of Prostate Tissue Using Descriptive Signal Enhancement Features Derived From DCE-MRI Acquisition With High Spatiotemporal Resolution

METHODS

A retrospective study (from January 2016 to July 2019) including 75 subjects (mean, 65 years; 46-80 years) with 2.5-second temporal resolution DCE-MRI and PIRADS 4 or 5 lesions was performed. Fifty-four subjects had biopsy-proven prostate cancer (Gleason 6, 15; Gleason 7, 20; Gleason 8, 13; Gleason 9, 6), whereas 21 subjects had negative MRI/ultrasound fusion-guided biopsies. Voxel-wise analysis of contrast signal enhancement was performed for all time points using custom-developed software, including automatic arterial input function detection. Seven descriptive parameter maps were calculated: normalized maximum signal intensity, time to start, time to maximum, time-to-maximum slope, and maximum slope with normalization on maximum signal and the arterial input function (SMN1, SMN2). The parameters were compared with ADC using multiparametric machine-learning models to determine classification accuracy. A Wilcoxon test was used for the hypothesis test and the Spearman coefficient for correlation.

RESULTS

There were significant differences (P < 0.05) for all 7 DCE-derived parameters between the normal peripheral zone versus PIRADS 4 or 5 lesions and the biopsy-positive versus biopsy-negative lesions. Multiparametric analysis showed better performance when combining ADC + DCE as input (accuracy/sensitivity/specificity, 97%/93%/100%) relative to ADC alone (accuracy/sensitivity/specificity, 94%/95%/95%) and to DCE alone (accuracy/sensitivity/specificity, 78%/79%/77%) in differentiating the normal peripheral zone from PIRADS lesions, biopsy-positive versus biopsy-negative lesions (accuracy/sensitivity/specificity, 68%/33%/81%), and Gleason 6 versus ≥7 prostate cancer (accuracy/sensitivity/specificity, 69%/60%/72%).

CONCLUSIONS

Descriptive perfusion characteristics derived from high-resolution DCE-MRI using model-free computations show significant differences between normal and cancerous tissue but do not reach the accuracy achieved with solely ADC-based classification. Combining ADC with DCE-based input features improved classification accuracy for PIRADS lesions, discrimination of biopsy-positive versus biopsy-negative lesions, and differentiation between Gleason 6 versus Gleason ≥7 lesions.

Magnetization-prepared GRASP MRI for rapid 3D T1 mapping and fat/water-separated T1 mapping

PURPOSE

This study aimed to (i) develop Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI using a stack-of-stars trajectory for rapid free-breathing T1 mapping and (ii) extend MP-GRASP to multi-echo acquisition (MP-Dixon-GRASP) for fat/water-separated (water-specific) T1 mapping.

METHODS

An adiabatic non-selective 180° inversion-recovery pulse was added to a gradient-echo-based golden-angle stack-of-stars sequence for magnetization-prepared 3D single-echo or 3D multi-echo acquisition. In combination with subspace-based GRASP-Pro reconstruction, the sequence allows for standard T1 mapping (MP-GRASP) or fat/water-separated T1 mapping (MP-Dixon-GRASP), respectively. The accuracy of T1 mapping using MP-GRASP was evaluated in a phantom and volunteers (brain and liver) against clinically accepted reference methods. The repeatability of T1 estimation was also assessed in the phantom and volunteers. The performance of MP-Dixon-GRASP for water-specific T1 mapping was evaluated in a fat/water phantom and volunteers (brain and liver).

RESULTS

ROI-based mean T1 values are correlated between the references and MP-GRASP in the phantom (R2 = 1.0), brain (R2 = 0.96), and liver (R2 = 0.73). MP-GRASP achieved good repeatability of T1 estimation in the phantom (R2 = 1.0), brain (R2 = 0.99), and liver (R2 = 0.82). Water-specific T1 is different from in-phase and out-of-phase composite T1 (composite T1 when fat and water signal are mixed in phase or out of phase) both in the phantom and volunteers.

CONCLUSION

This work demonstrated the initial performance of MP-GRASP and MP-Dixon-GRASP MRI for rapid 3D T1 mapping and 3D fat/water-separated T1 mapping in the brain (without motion) and in the liver (during free breathing). With fat/water-separated T1 estimation, MP-Dixon-GRASP could be potentially useful for imaging patients with fatty-liver diseases.

Intraindividual Comparison of Compressed Sensing-Accelerated Cartesian and Radial Arterial Phase Imaging of the Liver in an Experimental Tumor Model

OBJECTIVES

The aim of this study was to intraindividually compare the performance of 2 compressed sensing (CS)-accelerated magnetic resonance imaging (MRI) sequences, 1 featuring Cartesian (compressed sensing volumetric interpolated breath-hold examination [CS-VIBE]) and the other radial (golden-angle radial sparse parallel [GRASP]) k-space sampling in continuous dynamic imaging during hepatic vascular phases, using extracellular and hepatocyte-specific contrast agents.

MATERIALS AND METHODS

Seven New Zealand white rabbits, with induced VX2 liver tumors (median number of lesions, 2 ± 0.83; range, 1-3), received 2 continuously acquired T1-weighted prototype CS-accelerated MRI sequences (CS-VIBE and GRASP) with high spatial (0.8 × 0.8 × 1.5 mm) and temporal resolution (3.5 seconds) in randomized order on 2 separate days using a 1.5-T scanner. In all animals, imaging was performed using first gadobutrol at a dose of 0.1 mmol/kg and, then 45 minutes later, gadoxetic acid at a dose of 0.025 mmol/kg.The following qualitative parameters were assessed using 3- and 5-point Likert scales (3 and 5 being the highest scores respectively): image quality (IQ), arterial and venous vessel delineation, tumor enhancement, motion artifacts, and sequence-specific artifacts. Furthermore, the following quantitative parameters were obtained: relative peak signal enhancement, time to peak, mean transit time, and plasma flow ratios. Paired sampled t tests and Wilcoxon signed rank tests were used for intraindividual comparison. Image analysis was performed by 2 radiologists.

RESULTS

Six of 7 animals underwent the full imaging protocol and obtained data were analyzed statistically. Overall IQ was rated moderate to excellent, not differing significantly between the 2 sequences.Gadobutrol-enhanced CS-VIBE examinations revealed the highest mean Likert scale values in terms of vessel delineation and tumor enhancement (arterial 4.4 [4-5], venous 4.3 [3-5], and tumor 2.9 [2-3]). Significantly, more sequence-specific artifacts were seen in GRASP examinations (P = 0.008-0.031). However, these artifacts did not impair IQ. Excellent Likert scale ratings were found for motion artifacts in both sequences. In both sequences, a maximum of 4 hepatic arterial dominant phases were obtained. Regarding the relative peak signal enhancement, CS-VIBE and GRASP showed similar results. The relative peak signal enhancement values did not differ significantly between the 2 sequences in the aorta, the hepatic artery, or the inferior vena cava (P = 0.063-0.536). However, significantly higher values were noted for CS-VIBE in gadoxetic acid-enhanced examinations in the portal vein (P = 0.031) and regarding the tumor enhancement (P = 0.005). Time to peak and mean transit time or plasma flow ratios did not differ significantly between the sequences.

CONCLUSIONS

Both CS-VIBE and GRASP provide excellent results in dynamic liver MRI using extracellular and hepatocyte-specific contrast agents, in terms of IQ, peak signal intensity, and presence of artifacts.

Rapid 4D-MRI reconstruction using a deep radial convolutional neural network: Dracula

BACKGROUND AND PURPOSE

4D and midposition MRI could inform plan adaptation in lung and abdominal MR-guided radiotherapy. We present deep learning-based solutions to overcome long 4D-MRI reconstruction times while maintaining high image quality and short scan times.

METHODS

Two 3D U-net deep convolutional neural networks were trained to accelerate the 4D joint MoCo-HDTV reconstruction. For the first network, gridded and joint MoCo-HDTV-reconstructed 4D-MRI were used as input and target data, respectively, whereas the second network was trained to directly calculate the midposition image. For both networks, input and target data had dimensions of 256 × 256 voxels (2D) and 16 respiratory phases. Deep learning-based MRI were verified against joint MoCo-HDTV-reconstructed MRI using the structural similarity index (SSIM) and the naturalness image quality evaluator (NIQE). Moreover, two experienced observers contoured the gross tumour volume and scored the images in a blinded study.

RESULTS

For 12 subjects, previously unseen by the networks, high-quality 4D and midposition MRI (1.25 × 1.25 × 3.3 mm3) were each reconstructed from gridded images in only 28 seconds per subject. Excellent agreement was found between deep-learning-based and joint MoCo-HDTV-reconstructed MRI (average SSIM ≥ 0.96, NIQE scores 7.94 and 5.66). Deep-learning-based 4D-MRI were clinically acceptable for target and organ-at-risk delineation. Tumour positions agreed within 0.7 mm on midposition images.

CONCLUSION

Our results suggest that the joint MoCo-HDTV and midposition algorithms can each be approximated by a deep convolutional neural network. This rapid reconstruction of 4D and midposition MRI facilitates online treatment adaptation in thoracic or abdominal MR-guided radiotherapy.

Aliasing artifact reduction in spiral real-time MRI

PURPOSE

To mitigate a common artifact in spiral real-time MRI, caused by aliasing of signal outside the desired FOV. This artifact frequently occurs in midsagittal speech real-time MRI.

METHODS

Simulations were performed to determine the likely origin of the artifact. Two methods to mitigate the artifact are proposed. The first approach, denoted as "large FOV" (LF), keeps an FOV that is large enough to include the artifact signal source during reconstruction. The second approach, denoted as "estimation-subtraction" (ES), estimates the artifact signal source before subtracting a synthetic signal representing that source in multicoil k-space raw data. Twenty-five midsagittal speech-production real-time MRI data sets were used to evaluate both of the proposed methods. Reconstructions without and with corrections were evaluated by two expert readers using a 5-level Likert scale assessing artifact severity. Reconstruction time was also compared.

RESULTS

The origin of the artifact was found to be a combination of gradient nonlinearity and imperfect anti-aliasing in spiral sampling. The LF and ES methods were both able to substantially reduce the artifact, with an averaged qualitative score improvement of 1.25 and 1.35 Likert levels for LF correction and ES correction, respectively. Average reconstruction time without correction, with LF correction, and with ES correction were 160.69 ± 1.56, 526.43 ± 5.17, and 171.47 ± 1.71 ms/frame.

CONCLUSION

Both proposed methods were able to reduce the spiral aliasing artifacts, with the ES-reduction method being more effective and more time efficient.

Free-breathing radial imaging using a pilot-tone radiofrequency transmitter for detection of respiratory motion

PURPOSE

To describe an approach for detection of respiratory signals using a transmitted radiofrequency (RF) reference signal called Pilot-Tone (PT) and to use the PT signal for creation of motion-resolved images based on 3D stack-of-stars imaging under free-breathing conditions.

METHODS

This work explores the use of a reference RF signal generated by a small RF transmitter, placed outside the MR bore. The reference signal is received in parallel to the MR signal during each readout. Because the received PT amplitude is modulated by the subject's breathing pattern, a respiratory signal can be obtained by detecting the strength of the received PT signal over time. The breathing-induced PT signal modulation can then be used for reconstructing motion-resolved images from free-breathing scans. The PT approach was tested in volunteers using a radial stack-of-stars 3D gradient echo (GRE) sequence with golden-angle acquisition.

RESULTS

Respiratory signals derived from the proposed PT method were compared to signals from a respiratory cushion sensor and k-space-center-based self-navigation under different breathing conditions. Moreover, the accuracy was assessed using a modified acquisition scheme replacing the golden-angle scheme by a zero-angle acquisition. Incorporating the PT signal into eXtra-Dimensional (XD) motion-resolved reconstruction led to improved image quality and clearer anatomical depiction of the lung and liver compared to k-space-center signal and motion-averaged reconstruction, when binned into 6, 8, and 10 motion states.

CONCLUSION

PT is a novel concept for tracking respiratory motion. Its small dimension (8 cm), high sampling rate, and minimal interaction with the imaging scan offers great potential for resolving respiratory motion.

Noninvasive imaging assessment of portal hypertension

Portal hypertension (PH) is a spectrum of complications of chronic liver disease (CLD) and cirrhosis, with manifestations including ascites, gastroesophageal varices, splenomegaly, hypersplenism, hepatic hydrothorax, hepatorenal syndrome, hepatopulmonary syndrome and portopulmonary hypertension. PH can vary in severity and is diagnosed via invasive hepatic venous pressure gradient measurement (HVPG), which is considered the reference standard. Accurate diagnosis of PH and assessment of severity are highly relevant as patients with clinically significant portal hypertension (CSPH) are at higher risk for developing acute variceal bleeding and mortality. In this review, we discuss current and upcoming noninvasive imaging methods for diagnosis and assessment of severity of PH.

Delayed contrast dynamics as marker of regional impairment in pulmonary fibrosis using 5D MRI - a pilot study

OBJECTIVE

To analyse delayed contrast dynamics of fibrotic lesions in interstitial lung disease (ILD) using five dimensional (5D) MRI and to correlate contrast dynamics with disease severity.

METHODS

20 patients (mean age: 71 years; M:F, 13:7), with chronic fibrosing ILD: n = 12 idiopathic pulmonary fibrosis (IPF) and n = 8 non-IPF, underwent thin-section multislice CT as part of the standard diagnostic workup and additionally MRI of the lung. 2 min after contrast injection, a radial gradient echo sequence with golden-angle spacing was acquired during 5 min of free-breathing, followed by 5D image reconstruction. Disease was categorized as severe or non-severe according to CT morphological regional severity. For each patient, 10 lesions were analysed.

RESULTS

IPF lesions showed later peak enhancement compared to non-IPF (severe: p = 0.01, non-severe: p = 0.003). Severe lesions showed later peak enhancement compared to non-severe lesions, in non-IPF (p = 0.04), but not in IPF (p = 0.35). There was a tendency towards higher accumulation and washout rates in IPF compared to non-IPF in non-severe disease. Severe lesions had lower washout rate than non-severe ones in both IPF (p = 0.003) and non-IPF (p = 0.005). Continuous contrast agent accumulation, without washout, was found only in IPF lesions.

CONCLUSIONS

Contrast agent dynamics are influenced by type and severity of pulmonary fibrosis, which might enable a more thorough characterisation of disease burden. The regional impairment is of particular interest in the context of antifibrotic treatments and was characterised using a non-invasive, non-irradiating, free-breathing method.

ADVANCES IN KNOWLEDGE

Delayed contrast enhancement patterns allow the assessment of regional lung impairment which could represent different disease stages or phenotypes in ILD.

MRSIGMA: Magnetic Resonance SIGnature MAtching for real-time volumetric imaging

PURPOSE

To propose a real-time 3D MRI technique called MR SIGnature MAtching (MRSIGMA) for high-resolution volumetric imaging and motion tracking with very low imaging latency.

METHODS

MRSIGMA consists of two steps: (1) offline learning of a database of possible 3D motion states and corresponding motion signature ranges and (2) online matching of new motion signatures acquired in real time with prelearned motion states. Specifically, the offline learning step (non-real-time) reconstructs motion-resolved 4D images representing different motion states and assigns a unique motion range to each state. The online matching step (real-time) acquires motion signatures only and selects one of the prelearned 3D motion states for each newly acquired signature, which generates 3D images efficiently in real time. The MRSIGMA technique was evaluated on 15 golden-angle stack-of-stars liver data sets, and the performance of respiratory motion tracking with the online-generated real-time 3D MRI was compared with the corresponding 2D projections acquired in real time.

RESULTS

The total latency of generating each 3D image during online matching was about 300 ms, including acquisition of the motion signature data (~138 ms) and corresponding matching process (~150 ms). Linear correlation assessment suggested excellent correlation (R²  = 0.948) between motion displacement measured from the online-generated real-time 3D images and the 2D real-time projections.

CONCLUSION

This proof-of-concept study demonstrates the feasibility of MRSIGMA for high-resolution real-time volumetric imaging, which shifts the acquisition and reconstruction burden to an offline learning step and leaves fast online matching for online imaging with very low imaging latency. The MRSIGMA technique can potentially be used for real-time motion tracking in MRI-guided radiation therapy.

Improving dynamic contrast-enhanced MRI of the lung using motion-weighted sparse reconstruction: Initial experiences in patients

PURPOSE

The purpose of this study was to evaluate the performance of motion-weighted Golden-angle RAdial Sparse Parallel MRI (motion-weighted GRASP) for free-breathing dynamic contrast-enhanced MRI (DCE-MRI) of the lung.

METHODS

Motion-weighted GRASP incorporates a soft-gating motion compensation algorithm into standard GRASP reconstruction, so that motion-corrupted motion k-space (e.g., k-space acquired in inspiratory phases) contributes less to the final reconstructed images. Lung MR data from 20 patients (mean age = 57.9 ± 13.5) with known pulmonary lesions were retrospectively collected for this study. Each subject underwent a free-breathing DCE-MR scan using a fat-statured T₁-weighted stack-of-stars golden-angle radial sequence and a post-contrast breath-hold MR scan using a Cartesian volumetric-interpolated imaging sequence (BH-VIBE). Each radial dataset was reconstructed using GRASP without motion compensation and motion-weighted GRASP. All MR images were visually evaluated by two experienced radiologists blinded to reconstruction and acquisition schemes independently. In addition, the influence of motion-weighted reconstruction on dynamic contrast-enhancement patterns was also investigated.

RESULTS

For image quality assessment, motion-weighted GRASP received significantly higher visual scores than GRASP (P < 0.05) for overall image quality (3.68 vs. 3.39), lesion conspicuity (3.54 vs. 3.18) and overall artifact level (3.53 vs. 3.15). There was no significant difference (P > 0.05) between the breath-hold BH-VIBE and motion-weighted GRASP images. For assessment of temporal fidelity, motion-weighted GRASP maintained a good agreement with respect to GRASP.

CONCLUSION

Motion-weighted GRASP achieved better reconstruction performance in free-breathing DCE-MRI of the lung compared to standard GRASP, and it may enable improved assessment of pulmonary lesions.

SANTIS: Sampling-Augmented Neural neTwork with Incoherent Structure for MR image reconstruction

PURPOSE

To develop and evaluate a novel deep learning-based reconstruction framework called SANTIS (Sampling-Augmented Neural neTwork with Incoherent Structure) for efficient MR image reconstruction with improved robustness against sampling pattern discrepancy.

METHODS

With a combination of data cycle-consistent adversarial network, end-to-end convolutional neural network mapping, and data fidelity enforcement for reconstructing undersampled MR data, SANTIS additionally utilizes a sampling-augmented training strategy by extensively varying undersampling patterns during training, so that the network is capable of learning various aliasing structures and thereby removing undersampling artifacts more effectively and robustly. The performance of SANTIS was demonstrated for accelerated knee imaging and liver imaging using a Cartesian trajectory and a golden-angle radial trajectory, respectively. Quantitative metrics were used to assess its performance against different references. The feasibility of SANTIS in reconstructing dynamic contrast-enhanced images was also demonstrated using transfer learning.

RESULTS

Compared to conventional reconstruction that exploits image sparsity, SANTIS achieved consistently improved reconstruction performance (lower errors and greater image sharpness). Compared to standard learning-based methods without sampling augmentation (e.g., training with a fixed undersampling pattern), SANTIS provides comparable reconstruction performance, but significantly improved robustness, against sampling pattern discrepancy. SANTIS also achieved encouraging results for reconstructing liver images acquired at different contrast phases.

CONCLUSION

By extensively varying undersampling patterns, the sampling-augmented training strategy in SANTIS can remove undersampling artifacts more robustly. The novel concept behind SANTIS can particularly be useful for improving the robustness of deep learning-based image reconstruction against discrepancy between training and inference, an important, but currently less explored, topic.

GRASP-Pro: imProving GRASP DCE-MRI through self-calibrating subspace-modeling and contrast phase automation

PURPOSE

To propose a highly accelerated, high-resolution dynamic contrast-enhanced MRI (DCE-MRI) technique called GRASP-Pro (golden-angle radial sparse parallel imaging with imProved performance) through a joint sparsity and self-calibrating subspace constraint with automated selection of contrast phases.

METHODS

GRASP-Pro reconstruction enforces a combination of an explicit low-rank subspace-constraint and a temporal sparsity constraint. The temporal basis used to construct the subspace is learned from an intermediate reconstruction step using the low-resolution portion of radial k-space, which eliminates the need for generating the basis using auxiliary data or a physical signal model. A convolutional neural network was trained to generate the contrast enhancement curve in the artery, from which clinically relevant contrast phases are automatically selected for evaluation. The performance of GRASP-Pro was demonstrated for high spatiotemporal resolution DCE-MRI of the prostate and was compared against standard GRASP in terms of overall image quality, image sharpness, and residual streaks and/or noise level.

RESULTS

Compared to GRASP, GRASP-Pro reconstructed dynamic images with enhanced sharpness, less residual streaks and/or noise, and finer delineation of the prostate without prolonging reconstruction time. The image quality improvement reached statistical significance (P < 0.05) in all the assessment categories. The neural network successfully generated contrast enhancement curves in the artery, and corresponding peak enhancement indexes correlated well with that from the manual selection.

CONCLUSION

GRASP-Pro is a promising method for rapid and continuous DCE-MRI. It enables superior reconstruction performance over standard GRASP and allows reliable generation of artery enhancement curve to guide the selection of desired contrast phases for improving the efficiency of GRASP MRI workflow.

Radial streak artifact reduction using phased array beamforming

PURPOSE

A new method for streak artifact reduction in radial MRI based on phased array filtering.

THEORY

Radial imaging in applications that require large fields-of-view can be susceptible to streaking artifacts due to gradient nonlinearities. Coil removal methods prune the coils contributing the most to streaking artifacts at the expense of signal loss. Phased array beamforming is a form of spatial filtering used to suppress unwanted signals. The proposed method uses interference covariance generated from the streaking artifact samples which are manually extracted with phased array beamforming to suppress streaking in the images.

METHODS

The performance of the proposed method was evaluated on abdomen radial fast spin echo images acquired on a 1.5T Siemens scanner and compared with previously proposed methods.

RESULTS

Our results demonstrate that the proposed method can effectively suppress streaking artifacts without any noticeable loss in signal levels. Coil removal methods can suppress streaks as well but they may incur significant signal loss due to coil pruning. Quantitative metrics also demonstrate the superiority of the proposed method over earlier methods.

CONCLUSION

The use of interference covariance with phased array beamforming can help reduce streaking artifacts.

Free-breathing cine imaging with motion-corrected reconstruction at 3T using SPiral Acquisition with Respiratory correction and Cardiac Self-gating (SPARCS)

PURPOSE

To develop a continuous-acquisition cardiac self-gated spiral pulse sequence and a respiratory motion-compensated reconstruction strategy for free-breathing cine imaging.

METHODS

Cine data were acquired continuously on a 3T scanner for 8 seconds per slice without ECG gating or breath-holding, using a golden-angle gradient echo spiral pulse sequence. Cardiac motion information was extracted by applying principal component analysis on the gridded 8 × 8 k-space center data. Respiratory motion was corrected by rigid registration on each heartbeat. Images were reconstructed using a low-rank and sparse (L+S) technique. This strategy was evaluated in 37 healthy subjects and 8 subjects undergoing clinical cardiac MR studies. Image quality was scored (1-5 scale) in a blinded fashion by 2 experienced cardiologists. In 13 subjects with whole-heart coverage, left ventricular ejection fraction (LVEF) from SPiral Acquisition with Respiratory correction and Cardiac Self-gating (SPARCS) was compared to that from a standard ECG-gated breath-hold balanced steady-state free precession (bSSFP) cine sequence.

RESULTS

The self-gated signal was successfully extracted in all cases and demonstrated close agreement with the acquired ECG signal (mean bias, -0.22 ms). The mean image score across all subjects was 4.0 for reconstruction using the L+S model. There was good agreement between the LVEF derived from SPARCS and the gold-standard bSSFP technique.

CONCLUSION

SPARCS successfully images cardiac function without the need for ECG gating or breath-holding. With an 8-second data acquisition per slice, whole-heart cine images with clinically acceptable spatial and temporal resolution and image quality can be acquired in <90 seconds of free-breathing acquisition.

Pancreas deformation in the presence of tumors using feature tracking from free-breathing XD-GRASP MRI

BACKGROUND

Quantifying the biomechanical properties of pancreatic tumors could potentially help with assessment of tumor aggressiveness, prognosis, and prediction of therapy response.

PURPOSE

To quantify respiratory-induced deformation in the pancreas and pancreatic lesions using XD-GRASP (eXtra-Dimensional Golden-angle RAdial Sparse Parallel), MRI.

STUDY TYPE

Retrospective study where patients undergoing clinically indicated abdominal MRI which included free-breathing radial T₁ -weighted (T₁ W) imaging were studied.

SUBJECTS

Thirty-two patients (12 male and 20 female) including nine with pancreatic lesions constituted our study cohort.

FIELD STRENGTH/SEQUENCE

3.0 T with T₁ WI contrast-enhanced gradient echo radial free-breathing acquisition.

ASSESSMENT

Using the XD-GRASP imaging technique, the acquired free-breathing radial data were sorted and binned into 10 consecutive respiratory motion states that were jointly reconstructed. 3D deformation fields along the respiratory dimension were computed using an optical flow method and were analyzed in the pancreas.

STATISTICAL TESTS

The Wilcoxon signed-rank test was used to assess the difference in average displacement across pancreatic regions, while the Wilcoxon rank-sum test was used for displacement differences between patients with and without tumors. The interclass correlation coefficient (ICC) was computed to assess consistency between observers for each image quality measure.

RESULTS

There was a significantly larger displacement in the pancreatic tail compared with the head (8.2 ± 3.7 mm > 5.8 ± 2.4 mm; P < 0.001) and body regions (8.2 ± 3.7 mm > 6.6 ± 2.9 mm; P < 0.001). Furthermore, there was reduced normalized average displacement in patients with pancreatic lesions compared with subjects without lesions (0.33 ± 0.1 < 0.69 ± 0.26, P < 0.001 for the head; 0.30 ± 0.1 < 0.84 ± 0.31, P < 0.001 for the body; and 0.44 ± 0.31 < 1.08 ± 0.53, P < 0.001 for the tail, respectively).

DATA CONCLUSION

Free-breathing respiratory motion-sorted XD-GRASP MRI has the potential to noninvasively characterize the biomechanical properties of the pancreas by quantifying breathing-induced mechanical displacement.

LEVEL OF EVIDENCE

4 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;50:1633-1640.

Feasibility of multiple-view myocardial perfusion MRI using radial simultaneous multi-slice acquisitions

Purpose

Dynamic contrast enhanced MRI of the heart typically acquires 2–4 short-axis (SA) slices to detect and characterize coronary artery disease. This acquisition scheme is limited by incomplete coverage of the left ventricle. We studied the feasibility of using radial simultaneous multi-slice (SMS) technique to achieve SA, 2-chamber and/or 4-chamber long-axis (2CH LA and/or 4CH LA) coverage with and without electrocardiography (ECG) gating using a motion-robust reconstruction framework.

Methods

12 subjects were scanned at rest and/or stress, free breathing, with or without ECG gating. Multiple sets of radial SMS k-space were acquired within each cardiac cycle, and each SMS set sampled 3 parallel slices that were either SA, 2CH LA, or 4CH LA slices. The radial data was interpolated onto Cartesian space using an SMS GRAPPA operator gridding method. Self-gating and respiratory states binning of the data were done. The binning information as well as a pixel tracking spatiotemporal constrained reconstruction method were applied to obtain motion-robust image reconstructions. Reconstructions with and without the pixel tracking method were compared for signal-to-noise ratio and contrast-to-noise ratio.

Results

Full coverage of the heart (at least 3 SA and 3 LA slices) during the first pass of contrast at every heartbeat was achieved by using the radial SMS acquisition. The proposed pixel tracking reconstruction improves the average SNR and CNR by 21% and 30% respectively, and reduces temporal blurring for both gated and ungated acquisitions.

Conclusion

Acquiring simultaneous multi-slice SA, 2CH LA and/or 4CH LA myocardial perfusion images in every heartbeat is feasible in both gated and ungated acquisitions. This can add confidence when detecting and characterizing coronary artery disease by revealing ischemia in different views, and by providing apical coverage that is improved relative to SA slices alone. The proposed pixel tracking framework improves the reconstruction while adding little computational cost.

Rapid dynamic contrast-enhanced MRI for small animals at 7T using 3D ultra-short echo time and golden-angle radial sparse parallel MRI

PURPOSE

To develop a rapid dynamic contrast-enhanced MRI method with high spatial and temporal resolution for small-animal imaging at 7 Tesla.

METHODS

An ultra-short echo time (UTE) pulse sequence using a 3D golden-angle radial sampling was implemented to achieve isotropic spatial resolution with flexible temporal resolution. Continuously acquired radial spokes were grouped into subsets for image reconstruction using a multicoil compressed sensing approach (Golden-angle RAdial Sparse Parallel; GRASP). The proposed 3D-UTE-GRASP method with high temporal and spatial resolutions was tested using 7 mice with GL261 intracranial glioma models.

RESULTS

Iterative reconstruction with different temporal resolutions and regularization factors λ showed that, in all cases, the cost function decreased to less than 2.5% of its starting value within 20 iterations. The difference between the time-intensity curves of 3D-UTE-GRASP and nonuniform fast Fourier transform (NUFFT) images was minimal when λ was 1% of the maximum signal intensity of the initial NUFFT images. The 3D isotropic images were used to generate pharmacokinetic parameter maps to show the detailed images of the tumor characteristics in 3D and also to show longitudinal changes during tumor growth.

CONCLUSION

This feasibility study demonstrated that the proposed 3D-UTE-GRASP method can be used for effective measurement of the 3D spatial heterogeneity of tumor pharmacokinetic parameters.

Simultaneous Evaluation of Lung Anatomy and Ventilation Using 4D Respiratory-Motion-Resolved Ultrashort Echo Time Sparse MRI

BACKGROUND

Computed tomography (CT) and spirometry are the current standard methods for assessing lung anatomy and pulmonary ventilation, respectively. However, CT provides limited ventilation information and spirometry only provides global measures of lung ventilation. Thus, a method that can enable simultaneous examination of lung anatomy and ventilation is of clinical interest.

PURPOSE

To develop and test a 4D respiratory-resolved sparse lung MRI (XD-UTE: eXtra-Dimensional Ultrashort TE imaging) approach for simultaneous evaluation of lung anatomy and pulmonary ventilation.

STUDY TYPE

Prospective.

POPULATION

In all, 23 subjects (11 volunteers and 12 patients, mean age = 63.6 ± 8.4).

FIELD STRENGTH/SEQUENCE

3T MR; a prototype 3D golden-angle radial UTE sequence, a Cartesian breath-hold volumetric-interpolated examination (BH-VIBE) sequence.

ASSESSMENT

All subjects were scanned using the 3D golden-angle radial UTE sequence during normal breathing. Ten subjects underwent an additional scan during alternating normal and deep breathing. Respiratory-motion-resolved sparse reconstruction was performed for all the acquired data to generate dynamic normal-breathing or deep-breathing image series. For comparison, BH-VIBE was performed in 12 subjects. Lung images were visually scored by three experienced chest radiologists and were analyzed by two observers who segmented the left and right lung to derive ventilation parameters in comparison with spirometry.

STATISTICAL TESTS

Nonparametric paired two-tailed Wilcoxon signed-rank test; intraclass correlation coefficient, Pearson correlation coefficient.

RESULTS

XD-UTE achieved significantly improved image quality compared both with Cartesian BH-VIBE and radial reconstruction without motion compensation (P < 0.05). The global ventilation parameters (a sum of the left and right lung measures) were in good correlation with spirometry in the same subjects (correlation coefficient = 0.724). There were excellent correlations between the results obtained by two observers (intraclass correlation coefficient ranged from 0.8855-0.9995).

DATA CONCLUSION

Simultaneous evaluation of lung anatomy and ventilation using XD-UTE is demonstrated, which have shown good potential for improved diagnosis and management of patients with heterogeneous lung diseases.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019;49:411-422.

Improved 4D cardiac functional assessment for pediatric patients using motion-weighted image reconstruction

OBJECTIVE

Our aim was to develop and evaluate a motion-weighted reconstruction technique for improved cardiac function assessment in 4D magnetic resonance imaging (MRI).

MATERIALS AND METHODS

A flat-topped, two-sided Gaussian kernel was used to weigh k-space data in each target cardiac phase and adjacent two temporal phases during the proposed phase-by-phase reconstruction algorithm. The proposed method (Strategy 3) was used to reconstruct 18 cardiac phases based on data acquired using a previously proposed technique [4D multiphase steady-state imaging with contrast enhancement (MUSIC) technique and its self-gated extension using rotating Cartesian k-space (ROCK-MUSIC) from 12 pediatric patients. As a comparison, the same data set was reconstructed into nine phases using a phase-by-phase method (Strategy 1), 18 phases using view sharing (Strategy 4), and 18 phases using a temporal regularized method (Strategy 2). Regional image sharpness and left ventricle volumetric measurements were used to compare the four reconstructions quantitatively.

RESULTS

Strategies 1 and 4 generated significantly sharper images of static structures (P ≤ 0.018) than Strategies 2 and 3 but significantly more blurry (P ≤ 0.021) images of the heart. Left ventricular volumetric measurements from the nine-phase reconstruction (Strategy 1) correlated moderately (r < 0.8) with the 2D cine, whereas the remaining three techniques had a higher correlation (r > 0.9). The computational burden of Strategy 2 was six times that of Strategy 3.

CONCLUSION

The proposed method of motion-weighted reconstruction improves temporal resolution in 4D cardiac imaging with a clinically practical workflow.

Free-breathing dynamic contrast-enhanced MRI for assessment of pulmonary lesions using golden-angle radial sparse parallel imaging

BACKGROUND

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) has been shown to be a promising technique for assessing lung lesions. However, DCE-MRI often suffers from motion artifacts and insufficient imaging speed. Therefore, highly accelerated free-breathing DCE-MRI is of clinical interest for lung exams.

PURPOSE

To test the performance of rapid free-breathing DCE-MRI for simultaneous qualitative and quantitative assessment of pulmonary lesions using Golden-angle RAdial Sparse Parallel (GRASP) imaging.

STUDY TYPE

Prospective.

POPULATION

Twenty-six patients (17 males, mean age = 55.1 ± 14.4) with known pulmonary lesions.

FIELD STRENGTH/SEQUENCE

3T MR scanner; a prototype fat-saturated, T₁ -weighted stack-of-stars golden-angle radial sequence for data acquisition and a Cartesian breath-hold volumetric-interpolated examination (BH-VIBE) sequence for comparison.

ASSESSMENT

After a dual-mode GRASP reconstruction, one with 3-second temporal resolution (3s-GRASP) and the other with 15-second temporal resolution (15s-GRASP), all GRASP and BH-VIBE images were pooled together for blind assessment by two experienced radiologists, who independently scored the overall image quality, lesion delineation, overall artifact level, and diagnostic confidence of each case. Perfusion analysis was performed for the 3s-GRASP images using a Tofts model to generate the volume transfer coefficient (K^trans ) and interstitial volume (V_e ).

STATISTICAL TESTS

Nonparametric paired two-tailed Wilcoxon signed-rank test; Cohen's kappa; unpaired Student's t-test.

RESULTS

15s-GRASP achieved comparable image quality with conventional BH-VIBE (P > 0.05), except for the higher overall artifact level in the precontrast phase (P = 0.018). The K^trans and V_e in inflammation were higher than those in malignant lesions (K^trans : 0.78 ± 0.52 min^-1 vs. 0.37 ± 0.22 min^-1 , P = 0.020; V_e : 0.36 ± 0.16 vs. 0.26 ± 0.1, P = 0.177). Also, the K^trans and V_e in malignant lesions were also higher than those in benign lesions (K^trans : 0.37 ± 0.22 min^-1 vs. 0.04 ± 0.04 min^-1 , P = 0.001; V_e : 0.26 ± 0.12 vs. 0.10 ± 0.00, P = 0.063).

DATA CONCLUSION

This feasibility study demonstrated the performance of high spatiotemporal resolution free-breathing DCE-MRI of the lung using GRASP for qualitative and quantitative assessment of pulmonary lesions.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2018;48:459-468.