Simultaneous multislice diffusion-weighted imaging versus standard diffusion-weighted imaging in whole-body PET/MRI.

OBJECTIVE: To compare standard (STD-DWI) single-shot echo-planar imaging DWI and simultaneous multislice (SMS) DWI during whole-body positron emission tomography (PET)/MRI regarding acquisition time, image quality, and lesion detection. METHODS: Eighty-three adults (47 females, 57%), median age of 64 years (IQR 52-71), were prospectively enrolled from August 2018 to March 2020. Inclusion criteria were (a) abdominal or pelvic tumors and (b) PET/MRI referral from a clinician. Patients were excluded if whole-body acquisition of STD-DWI and SMS-DWI sequences was not completed. The evaluated sequences were axial STD-DWI at b-values 50-400-800 s/mm2 and the apparent diffusion coefficient (ADC), and axial SMS-DWI at b-values 50-300-800 s/mm2 and ADC, acquired with a 3-T PET/MRI scanner. Three radiologists rated each sequence's quality on a five-point scale. Lesion detection was quantified using the anatomic MRI sequences and PET as the reference standard. Regression models were constructed to quantify the association between all imaging outcomes/scores and sequence type. RESULTS: The median whole-body STD-DWI acquisition time was 14.8 min (IQR 14.1-16.0) versus 7.0 min (IQR 6.7-7.2) for whole-body SMS-DWI, p < 0.001. SMS-DWI image quality scores were higher than STD-DWI in the abdomen (OR 5.31, 95% CI 2.76-10.22, p < 0.001), but lower in the cervicothoracic junction (OR 0.21, 95% CI 0.10-0.43, p < 0.001). There was no significant difference in the chest, mediastinum, pelvis, and rectum. STD-DWI detected 276/352 (78%) lesions while SMS-DWI located 296/352 (84%, OR 1.46, 95% CI 1.02-2.07, p = 0.038). CONCLUSIONS: In cancer staging and restaging, SMS-DWI abbreviates acquisition while maintaining or improving the diagnostic yield in most anatomic regions. KEY POINTS: • Simultaneous multislice diffusion-weighted imaging enables faster whole-body image acquisition. • Simultaneous multislice diffusion-weighted imaging maintains or improves image quality when compared to single-shot echo-planar diffusion-weighted imaging in most anatomical regions. • Simultaneous multislice diffusion-weighted imaging leads to superior lesion detection.

MR-assisted PET respiratory motion correction using deep-learning based short-scan motion fields.

PURPOSE: We evaluated the impact of PET respiratory motion correction (MoCo) in a phantom and patients. Moreover, we proposed and examined a PET MoCo approach using motion vector fields (MVFs) from a deep-learning reconstructed short MRI scan. METHODS: The evaluation of PET MoCo was performed in a respiratory motion phantom study with varying lesion sizes and tumor to background ratios (TBRs) using a static scan as the ground truth. MRI-based MVFs were derived from either 2000 spokes (MoCo2000 , 5-6 min acquisition time) using a Fourier transform reconstruction or 200 spokes (MoCoP2P200 , 30-40 s acquisition time) using a deep-learning Phase2Phase (P2P) reconstruction and then incorporated into PET MoCo reconstruction. For six patients with hepatic lesions, the performance of PET MoCo was evaluated using quantitative metrics (SUVmax , SUVpeak , SUVmean , lesion volume) and a blinded radiological review on lesion conspicuity. RESULTS: MRI-assisted PET MoCo methods provided similar results to static scans across most lesions with varying TBRs in the phantom. Both MoCo2000 and MoCoP2P200 PET images had significantly higher SUVmax , SUVpeak , SUVmean and significantly lower lesion volume than non-motion-corrected (non-MoCo) PET images. There was no statistical difference between MoCo2000 and MoCoP2P200 PET images for SUVmax , SUVpeak , SUVmean or lesion volume. Both radiological reviewers found that MoCo2000 and MoCoP2P200 PET significantly improved lesion conspicuity. CONCLUSION: An MRI-assisted PET MoCo method was evaluated using the ground truth in a phantom study. In patients with hepatic lesions, PET MoCo images improved quantitative and qualitative metrics based on only 30-40 s of MRI motion modeling data.

Approximating anatomically-guided PET reconstruction in image space using a convolutional neural network.

In the last two decades, it has been shown that anatomically-guided PET reconstruction can lead to improved bias-noise characteristics in brain PET imaging. However, despite promising results in simulations and first studies, anatomically-guided PET reconstructions are not yet available for use in routine clinical because of several reasons. In light of this, we investigate whether the improvements of anatomically-guided PET reconstruction methods can be achieved entirely in the image domain with a convolutional neural network (CNN). An entirely image-based CNN post-reconstruction approach has the advantage that no access to PET raw data is needed and, moreover, that the prediction times of trained CNNs are extremely fast on state of the art GPUs which will substantially facilitate the evaluation, fine-tuning and application of anatomically-guided PET reconstruction in real-world clinical settings. In this work, we demonstrate that anatomically-guided PET reconstruction using the asymmetric Bowsher prior can be well-approximated by a purely shift-invariant convolutional neural network in image space allowing the generation of anatomically-guided PET images in almost real-time. We show that by applying dedicated data augmentation techniques in the training phase, in which 16 [18F]FDG and 10 [18F]PE2I data sets were used, lead to a CNN that is robust against the used PET tracer, the noise level of the input PET images and the input MRI contrast. A detailed analysis of our CNN in 36 [18F]FDG, 18 [18F]PE2I, and 7 [18F]FET test data sets demonstrates that the image quality of our trained CNN is very close to the one of the target reconstructions in terms of regional mean recovery and regional structural similarity.

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.

Cardiorespiratory motion-tracking via self-refocused rosette navigators.

PURPOSE: To develop a flexible method for tracking respiratory and cardiac motions throughout MR and PET-MR body examinations that requires no additional hardware and minimal sequence modification. METHODS: The incorporation of a contrast-neutral rosette navigator module following the RF excitation allows for robust cardiorespiratory motion tracking with minimal impact on the host sequence. Spatial encoding gradients are applied to the FID signal and the desired motion signals are extracted with a blind source separation technique. This approach is validated with an anthropomorphic, PET-MR-compatible motion phantom as well as in 13 human subjects. RESULTS: Both respiratory and cardiac motions were reliably extracted from the proposed rosette navigator in phantom and patient studies. In the phantom study, the MR-derived motion signals were additionally validated against the ground truth measurement of diaphragm displacement and left ventricle model triggering pulse. CONCLUSION: The proposed method yields accurate respiratory and cardiac motion-state tracking, requiring only a short (1.76 ms) additional navigator module, which is self-refocusing and imposes minimal constraints on sequence design.

Undersampling artifact reduction for free-breathing 3D stack-of-radial MRI based on a deep adversarial learning network.

PURPOSE: Stack-of-radial MRI allows free-breathing abdominal scans, however, it requires relatively long acquisition time. Undersampling reduces scan time but can cause streaking artifacts and degrade image quality. This study developed deep learning networks with adversarial loss and evaluated the performance of reducing streaking artifacts and preserving perceptual image sharpness. METHODS: A 3D generative adversarial network (GAN) was developed for reducing streaking artifacts in stack-of-radial abdominal scans. Training and validation datasets were self-gated to 5 respiratory states to reduce motion artifacts and to effectively augment the data. The network used a combination of three loss functions to constrain the anatomy and preserve image quality: adversarial loss, mean-squared-error loss and structural similarity index loss. The performance of the network was investigated for 3-5 times undersampled data from 2 institutions. The performance of the GAN for 5 times accelerated images was compared with a 3D U-Net and evaluated using quantitative NMSE, SSIM and region of interest (ROI) measurements as well as qualitative scores of radiologists. RESULTS: The 3D GAN showed similar NMSE (0.0657 vs. 0.0559, p = 0.5217) and significantly higher SSIM (0.841 vs. 0.798, p < 0.0001) compared to U-Net. ROI analysis showed GAN removed streaks in both the background air and the tissue and was not significantly different from the reference mean and variations. Radiologists' scores showed GAN had a significant improvement of 1.6 point (p = 0.004) on a 4-point scale in streaking score while no significant difference in sharpness score compared to the input. CONCLUSION: 3D GAN removes streaking artifacts and preserves perceptual image details.

Phase2Phase: Respiratory Motion-Resolved Reconstruction of Free-Breathing Magnetic Resonance Imaging Using Deep Learning Without a Ground Truth for Improved Liver Imaging.

OBJECTIVES: Respiratory binning of free-breathing magnetic resonance imaging data reduces motion blurring; however, it exacerbates noise and introduces severe artifacts due to undersampling. Deep neural networks can remove artifacts and noise but usually require high-quality ground truth images for training. This study aimed to develop a network that can be trained without this requirement. MATERIALS AND METHODS: This retrospective study was conducted on 33 participants enrolled between November 2016 and June 2019. Free-breathing magnetic resonance imaging was performed using a radial acquisition. Self-navigation was used to bin the k-space data into 10 respiratory phases. To simulate short acquisitions, subsets of radial spokes were used in reconstructing images with multicoil nonuniform fast Fourier transform (MCNUFFT), compressed sensing (CS), and 2 deep learning methods: UNet3DPhase and Phase2Phase (P2P). UNet3DPhase was trained using a high-quality ground truth, whereas P2P was trained using noisy images with streaking artifacts. Two radiologists blinded to the reconstruction methods independently reviewed the sharpness, contrast, and artifact-freeness of the end-expiration images reconstructed from data collected at 16% of the Nyquist sampling rate. The generalized estimating equation method was used for statistical comparison. Motion vector fields were derived to examine the respiratory motion range of 4-dimensional images reconstructed using different methods. RESULTS: A total of 15 healthy participants and 18 patients with hepatic malignancy (50 ± 15 years, 6 women) were enrolled. Both reviewers found that the UNet3DPhase and P2P images had higher contrast (P < 0.01) and fewer artifacts (P < 0.01) than the CS images. The UNet3DPhase and P2P images were reported to be sharper than the CS images by 1 reviewer (P < 0.01) but not by the other reviewer (P = 0.22, P = 0.18). UNet3DPhase and P2P were similar in sharpness and contrast, whereas UNet3DPhase had fewer artifacts (P < 0.01). The motion vector lengths for the MCNUFFT800 and P2P800 images were comparable (10.5 ± 4.2 mm and 9.9 ± 4.0 mm, respectively), whereas both were significantly larger than CS2000 (7.0 ± 3.9 mm; P < 0.0001) and UNnet3DPhase800 (6.9 ± 3.2; P < 0.0001) images. CONCLUSIONS: Without a ground truth, P2P can reconstruct sharp, artifact-free, and high-contrast respiratory motion-resolved images from highly undersampled data. Unlike the CS and UNet3DPhase methods, P2P did not artificially reduce the respiratory motion range.

Respiratory Motion Detection and Correction for MR Using the Pilot Tone: Applications for MR and Simultaneous PET/MR Examinations.

OBJECTIVES: The aim of this study was to develop a method for tracking respiratory motion throughout full MR or PET/MR studies that requires only minimal additional hardware and no modifications to the sequences. MATERIALS AND METHODS: Patient motion that is caused by respiration affects the quality of the signal of the individual radiofrequency receive coil elements. This effect can be detected as a modulation of a monofrequent signal that is emitted by a small portable transmitter placed inside the bore (Pilot Tone). The frequency is selected such that it is located outside of the frequency band of the actual MR readout experiment but well within the bandwidth of the radiofrequency receiver, that is, the oversampling area. Temporal variations of the detected signal indicate motion. After extraction of the signal from the raw data, principal component analysis was used to identify respiratory motion. The approach and potential applications during MR and PET/MR examinations that rely on a continuous respiratory signal were validated with an anthropomorphic, PET/MR-compatible motion phantom as well as in a volunteer study. RESULTS: Respiratory motion detection and correction were presented for MR and PET data in phantom and volunteer studies. The Pilot Tone successfully recovered the ground-truth respiratory signal provided by the phantom. CONCLUSIONS: The presented method provides reliable respiratory motion tracking during arbitrary imaging sequences throughout a full PET/MR study. All results can directly be transferred to MR-only applications as well.

Image reconstruction for interrupted-beam x-ray CT on diagnostic clinical scanners.

Low-dose x-ray CT is a major research area with high clinical impact. Compressed sensing using view-based sparse sampling and sparsity-promoting regularization has shown promise in simulations, but these methods can be difficult to implement on diagnostic clinical CT scanners since the x-ray beam cannot be switched on and off rapidly enough. An alternative to view-based sparse sampling is interrupted-beam sparse sampling. SparseCT is a recently-proposed interrupted-beam scheme that achieves sparse sampling by blocking a portion of the beam using a multislit collimator (MSC). The use of an MSC necessitates a number of modifications to the standard compressed sensing reconstruction pipeline. In particular, we find that SparseCT reconstruction is feasible within a model-based image reconstruction framework that incorporates data fidelity weighting to consider penumbra effects and source jittering to consider the effect of partial source obstruction. Here, we present these modifications and demonstrate their application in simulations and real-world prototype scans. In simulations compared to conventional low-dose acquisitions, SparseCT is able to achieve smaller normalized root-mean square differences and higher structural similarity measures on two reduction factors. In prototype experiments, we successfully apply our reconstruction modifications and maintain image resolution at quarter-dose reduction level. The SparseCT design requires only small hardware modifications to current diagnostic clinical scanners, opening up new possibilities for CT dose reduction.

Joint Reconstruction of Activity and Attenuation in Time-of-Flight PET: A Quantitative Analysis.

Methods for joint activity reconstruction and attenuation reconstruction of time-of-flight (TOF) PET data provide an effective solution to attenuation correction when no (or incomplete or inaccurate) information on attenuation is available. One of the main barriers limiting use of these methods in clinical practice is their lack of validation in a relatively large patient database. In this contribution, we aim to validate reconstruction performed with maximum-likelihood activity reconstruction and attenuation registration (MLRR) in a whole-body patient dataset. Furthermore, a partial validation (because the scale problem of the algorithm is avoided for now) of reconstruction performed with maximum-likelihood activity and attenuation (MLAA) is also provided. We present a quantitative comparison between these 2 methods of joint reconstruction and the current clinical gold standard, maximum-likelihood expectation maximization (MLEM) with CT-based attenuation correction. Methods: The whole-body TOF PET emission data of each patient dataset were processed as a whole to reconstruct an activity volume covering all the acquired bed positions, helping reduce the problem of a scale per bed position in MLAA to a global scale for the entire activity volume. Three reconstruction algorithms were used: MLEM, MLRR, and MLAA. A maximum-likelihood scaling of the single-scatter simulation estimate to the emission data was used for scatter correction. The reconstruction results for various regions of interest were then analyzed. Results: The joint reconstructions of the whole-body patient dataset provided better quantification than the gold standard in cases of PET and CT misalignment caused by patient or organ motion. Our quantitative analysis showed a difference of -4.2% ± 2.3% between MLRR and MLEM and a difference of -7.5% ± 4.6% between MLAA and MLEM, averaged over all regions of interest. Conclusion: Joint reconstruction of activity and attenuation provides a useful means to estimate tracer distribution when CT-based-attenuation images are subject to misalignment or are not available. With an accurate estimate of the scatter contribution in the emission measurements, the joint reconstructions of TOF PET data are within clinically acceptable accuracy.

Visual detection of regional brain hypometabolism in cognitively impaired patients is independent of positron emission tomography-magnetic resonance attenuation correction method.

Fluorodeoxyglucose (FDG) positron emission tomography-magnetic resonance (PET/MR) is useful for the evaluation of cognitively-impaired patients. This study aims to assess two different attenuation correction (AC) methods (Dixon-MR and atlas-based) versus index-standard computed tomography (CT) AC for the visual interpretation of regional hypometabolism in patients with cognitive impairment. Two board-certified nuclear medicine physicians blindly scored brain region FDG hypometabolism as normal versus hypometabolic using two-dimensional (2D) and 3D FDG PET/MR images generated by MIM software. Regions were quantitatively assessed as normal versus mildly, moderately, or severely hypometabolic. Hypometabolism scores obtained using the different methods of AC were compared, and interreader, as well as intra-reader agreement, was assessed. Regional hypometabolism versus normal metabolism was correctly classified in 16 patients on atlas-based and Dixon-based AC map PET reconstructions (vs. CT reference AC) for 94% (90%-96% confidence interval [CI]) and 93% (89%-96% CI) of scored regions, respectively. The averaged sensitivity/specificity for detection of any regional hypometabolism was 95%/94% (P = 0.669) and 90%/91% (P = 0.937) for atlas-based and Dixon-based AC maps. Interreader agreement for detection of regional hypometabolism was high, with similar outcome assessments when using atlas- and Dixon-corrected PET data in 93% (Κ =0.82) and 93% (Κ =0.84) of regions, respectively. Intrareader agreement for detection of regional hypometabolism was high, with concordant outcome assessments when using atlas- and Dixon-corrected data in 93%/92% (Κ =0.79) and 92/93% (Κ =0.78). Despite the quantitative advantages of atlas-based AC in brain PET/MR, routine clinical Dixon AC yields comparable visual ratings of regional hypometabolism in the evaluation of cognitively impaired patients undergoing brain PET/MR and is similar in performance to CT-based AC. Therefore, Dixon AC is acceptable for the routine clinical evaluation of dementia syndromes.

Plane-dependent ML scatter scaling: 3D extension of the 2D simulated single scatter (SSS) estimate.

Scatter correction is typically done using a simulation of the single scatter, which is then scaled to account for multiple scatters and other possible model mismatches. This scaling factor is determined by fitting the simulated scatter sinogram to the measured sinogram, using only counts measured along LORs that do not intersect the patient body, i.e. 'scatter-tails'. Extending previous work, we propose to scale the scatter with a plane dependent factor, which is determined as an additional unknown in the maximum likelihood (ML) reconstructions, using counts in the entire sinogram rather than only the 'scatter-tails'. The ML-scaled scatter estimates are validated using a Monte-Carlo simulation of a NEMA-like phantom, a phantom scan with typical contrast ratios of a 68Ga-PSMA scan, and 23 whole-body 18F-FDG patient scans. On average, we observe a 12.2% change in the total amount of tracer activity of the MLEM reconstructions of our whole-body patient database when the proposed ML scatter scales are used. Furthermore, reconstructions using the ML-scaled scatter estimates are found to eliminate the typical 'halo' artifacts that are often observed in the vicinity of high focal uptake regions.

Dosimetric evaluation of synthetic CT for magnetic resonance-only based radiotherapy planning of lung cancer.

BACKGROUND: Interest in MR-only treatment planning for radiation therapy is growing rapidly with the emergence of integrated MRI/linear accelerator technology. The purpose of this study was to evaluate the feasibility of using synthetic CT images generated from conventional Dixon-based MRI scans for radiation treatment planning of lung cancer. METHODS: Eleven patients who underwent whole-body PET/MR imaging following a PET/CT exam were randomly selected from an ongoing prospective IRB-approved study. Attenuation maps derived from the Dixon MR Images and atlas-based method was used to create CT data (synCT). Treatment planning for radiation treatment of lung cancer was optimized on the synCT and subsequently copied to the registered CT (planCT) for dose calculation. Planning target volumes (PTVs) with three sizes and four different locations in the lung were planned for irradiation. The dose-volume metrics comparison and 3D gamma analysis were performed to assess agreement between the synCT and CT calculated dose distributions. RESULTS: Mean differences between PTV doses on synCT and CT across all the plans were -0.1% ± 0.4%, 0.1% ± 0.5%, and 0.4% ± 0.5% for D95, D98 and D100, respectively. Difference in dose between the two datasets for organs at risk (OARs) had average differences of -0.14 ± 0.07 Gy, 0.0% ± 0.1%, and -0.1% ± 0.2% for maximum spinal cord, lung V20, and heart V40 respectively. In patient groups based on tumor size and location, no significant differences were observed in the PTV and OARs dose-volume metrics (p > 0.05), except for the maximum spinal-cord dose when the target volumes were located at the lung apex (p = 0.001). Gamma analysis revealed a pass rate of 99.3% ± 1.1% for 2%/2 mm (dose difference/distance to agreement) acceptance criteria in every plan. CONCLUSIONS: The synCT generated from Dixon-based MRI allows for dose calculation of comparable accuracy to the standard CT for lung cancer treatment planning. The dosimetric agreement between synCT and CT calculated doses warrants further development of a MR-only workflow for radiotherapy of lung cancer.