Evaluating a radiotherapy deep learning synthetic CT algorithm for PET-MR attenuation correction in the pelvis.

BACKGROUND: Positron emission tomography-magnetic resonance (PET-MR) attenuation correction is challenging because the MR signal does not represent tissue density and conventional MR sequences cannot image bone. A novel zero echo time (ZTE) MR sequence has been previously developed which generates signal from cortical bone with images acquired in 65 s. This has been combined with a deep learning model to generate a synthetic computed tomography (sCT) for MR-only radiotherapy. This study aimed to evaluate this algorithm for PET-MR attenuation correction in the pelvis. METHODS: Ten patients being treated with ano-rectal radiotherapy received a [Formula: see text]F-FDG-PET-MR in the radiotherapy position. Attenuation maps were generated from ZTE-based sCT (sCTAC) and the standard vendor-supplied MRAC. The radiotherapy planning CT scan was rigidly registered and cropped to generate a gold standard attenuation map (CTAC). PET images were reconstructed using each attenuation map and compared for standard uptake value (SUV) measurement, automatic thresholded gross tumour volume (GTV) delineation and GTV metabolic parameter measurement. The last was assessed for clinical equivalence to CTAC using two one-sided paired t tests with a significance level corrected for multiple testing of [Formula: see text]. Equivalence margins of [Formula: see text] were used. RESULTS: Mean whole-image SUV differences were -0.02% (sCTAC) compared to -3.0% (MRAC), with larger differences in the bone regions (-0.5% to -16.3%). There was no difference in thresholded GTVs, with Dice similarity coefficients [Formula: see text]. However, there were larger differences in GTV metabolic parameters. Mean differences to CTAC in [Formula: see text] were [Formula: see text] (± standard error, sCTAC) and [Formula: see text] (MRAC), and [Formula: see text] (sCTAC) and [Formula: see text] (MRAC) in [Formula: see text]. The sCTAC was statistically equivalent to CTAC within a [Formula: see text] equivalence margin for [Formula: see text] and [Formula: see text] ([Formula: see text] and [Formula: see text]), whereas the MRAC was not ([Formula: see text] and [Formula: see text]). CONCLUSION: Attenuation correction using this radiotherapy ZTE-based sCT algorithm was substantially more accurate than current MRAC methods with only a 40 s increase in MR acquisition time. This did not impact tumour delineation but did significantly improve the accuracy of whole-image and tumour SUV measurements, which were clinically equivalent to CTAC. This suggests PET images reconstructed with sCTAC would enable accurate quantitative PET images to be acquired on a PET-MR scanner.

Region of interest focused MRI to synthetic CT translation using regression and segmentation multi-task network.

Objective. In MR-only clinical workflow, replacing CT with MR image is of advantage for workflow efficiency and reduces radiation to the patient. An important step required to eliminate CT scan from the workflow is to generate the information provided by CT via an MR image. In this work, we aim to demonstrate a method to generate accurate synthetic CT (sCT) from an MR image to suit the radiation therapy (RT) treatment planning workflow. We show the feasibility of the method and make way for a broader clinical evaluation.Approach. We present a machine learning method for sCT generation from zero-echo-time (ZTE) MRI aimed at structural and quantitative accuracies of the image, with a particular focus on the accurate bone density value prediction. The misestimation of bone density in the radiation path could lead to unintended dose delivery to the target volume and results in suboptimal treatment outcome. We propose a loss function that favors a spatially sparse bone region in the image. We harness the ability of the multi-task network to produce correlated outputs as a framework to enable localization of region of interest (RoI) via segmentation, emphasize regression of values within RoI and still retain the overall accuracy via global regression. The network is optimized by a composite loss function that combines a dedicated loss from each task.Main results. We have included 54 brain patient images in this study and tested the sCT images against reference CT on a subset of 20 cases. A pilot dose evaluation was performed on 9 of the 20 test cases to demonstrate the viability of the generated sCT in RT planning. The average quantitative metrics produced by the proposed method over the test set were-(a) mean absolute error (MAE) of 70 ± 8.6 HU; (b) peak signal-to-noise ratio (PSNR) of 29.4 ± 2.8 dB; structural similarity metric (SSIM) of 0.95 ± 0.02; and (d) Dice coefficient of the body region of 0.984 ± 0.Significance. We demonstrate that the proposed method generates sCT images that resemble visual characteristics of a real CT image and has a quantitative accuracy that suits RT dose planning application. We compare the dose calculation from the proposed sCT and the real CT in a radiation therapy treatment planning setup and show that sCT based planning falls within 0.5% target dose error. The method presented here with an initial dose evaluation makes an encouraging precursor to a broader clinical evaluation of sCT based RT planning on different anatomical regions.

Comprehensive dose evaluation of a Deep Learning based synthetic Computed Tomography algorithm for pelvic Magnetic Resonance-only radiotherapy.

BACKGROUND AND PURPOSE: Magnetic Resonance (MR)-only radiotherapy enables the use of MR without the uncertainty of MR-Computed Tomography (CT) registration. This requires a synthetic CT (sCT) for dose calculations, which can be facilitated by a novel Zero Echo Time (ZTE) sequence where bones are visible and images are acquired in 65 seconds. This study evaluated the dose calculation accuracy for pelvic sites of a ZTE-based Deep Learning sCT algorithm developed by GE Healthcare. MATERIALS AND METHODS: ZTE and CT images were acquired in 56 pelvic radiotherapy patients in the radiotherapy position. A 2D U-net convolutional neural network was trained using pairs of deformably registered CT and ZTE images from 36 patients. In the remaining 20 patients the dosimetric accuracy of the sCT was assessed using cylindrical dummy Planning Target Volumes (PTVs) positioned at four different central axial locations, as well as the clinical treatment plans (for prostate (n = 10), rectum (n = 4) and anus (n = 6) cancers). The sCT was rigidly and deformably registered, the plan recalculated and the doses compared using mean differences and gamma analysis. RESULTS: Mean dose differences to the PTV D98% were ≤ 0.5% for all dummy PTVs and clinical plans (rigid registration). Mean gamma pass rates at 1%/1 mm were 98.0 ± 0.4% (rigid) and 100.0 ± 0.0% (deformable), 96.5 ± 0.8% and 99.8 ± 0.1%, and 95.4 ± 0.6% and 99.4 ± 0.4% for the clinical prostate, rectum and anus plans respectively. CONCLUSIONS: A ZTE-based sCT algorithm with high dose accuracy throughout the pelvis has been developed. This suggests the algorithm is sufficiently accurate for MR-only radiotherapy for all pelvic sites.

In-phase zero TE musculoskeletal imaging.

PURPOSE: To introduce a new method for in-phase zero TE (ipZTE) musculoskeletal MR imaging. METHODS: ZTE is a 3D radial imaging method, which is sensitive to chemical shift off-resonance signal interference, especially around fat-water tissue interfaces. The ipZTE method addresses this fat-water chemical shift artifact by acquiring each 3D radial spoke at least twice with varying readout gradient amplitude and hence varying effective sampling time. Using k-space-based chemical shift decomposition, the acquired data is then reconstructed into an in-phase ZTE image and an out-of-phase disturbance. RESULTS: The ipZTE method was tested for knee, pelvis, brain, and whole-body. The obtained images demonstrate exceptional soft-tissue uniformity free from out-of-phase disturbances apparent in the original ZTE images. The chemical shift decomposition was found to improve SNR at the cost of reduced image resolution. CONCLUSION: The ipZTE method can be used as an averaging mechanism to eliminate fat-water chemical shift artifacts and improve SNR. The method is expected to improve ZTE-based musculoskeletal imaging and pseudo CT conversion as required for PET/MR attenuation correction and MR-guided radiation therapy planning.

Zero TE-based pseudo-CT image conversion in the head and its application in PET/MR attenuation correction and MR-guided radiation therapy planning.

PURPOSE: To describe a method for converting Zero TE (ZTE) MR images into X-ray attenuation information in the form of pseudo-CT images and demonstrate its performance for (1) attenuation correction (AC) in PET/MR and (2) dose planning in MR-guided radiation therapy planning (RTP). METHODS: Proton density-weighted ZTE images were acquired as input for MR-based pseudo-CT conversion, providing (1) efficient capture of short-lived bone signals, (2) flat soft-tissue contrast, and (3) fast and robust 3D MR imaging. After bias correction and normalization, the images were segmented into bone, soft-tissue, and air by means of thresholding and morphological refinements. Fixed Hounsfield replacement values were assigned for air (-1000 HU) and soft-tissue (+42 HU), whereas continuous linear mapping was used for bone. RESULTS: The obtained ZTE-derived pseudo-CT images accurately resembled the true CT images (i.e., Dice coefficient for bone overlap of 0.73 ± 0.08 and mean absolute error of 123 ± 25 HU evaluated over the whole head, including errors from residual registration mismatches in the neck and mouth regions). The linear bone mapping accounted for bone density variations. Averaged across five patients, ZTE-based AC demonstrated a PET error of -0.04 ± 1.68% relative to CT-based AC. Similarly, for RTP assessed in eight patients, the absolute dose difference over the target volume was found to be 0.23 ± 0.42%. CONCLUSION: The described method enables MR to pseudo-CT image conversion for the head in an accurate, robust, and fast manner without relying on anatomical prior knowledge. Potential applications include PET/MR-AC, and MR-guided RTP.

Fourier domain image fusion for differential X-ray phase-contrast breast imaging.

X-Ray Phase-Contrast (XPC) imaging is a novel technology with a great potential for applications in clinical practice, with breast imaging being of special interest. This work introduces an intuitive methodology to combine and visualize relevant diagnostic features, present in the X-ray attenuation, phase shift and scattering information retrieved in XPC imaging, using a Fourier domain fusion algorithm. The method allows to present complementary information from the three acquired signals in one single image, minimizing the noise component and maintaining visual similarity to a conventional X-ray image, but with noticeable enhancement in diagnostic features, details and resolution. Radiologists experienced in mammography applied the image fusion method to XPC measurements of mastectomy samples and evaluated the feature content of each input and the fused image. This assessment validated that the combination of all the relevant diagnostic features, contained in the XPC images, was present in the fused image as well.

Low-dose, phase-contrast mammography with high signal-to-noise ratio.

Differential phase-contrast X-ray imaging using a Talbot-Lau interferometer has recently shown promising results for applications in medical imaging. However, reducing the applied radiation dose remains a major challenge. In this study, we consider the realization of a Talbot-Lau interferometer in a high Talbot order to increase the signal-to-noise ratio for low-dose applications. The quantitative performance of π and π/2 systems at high Talbot orders is analyzed through simulations, and the design energy and X-ray spectrum are optimized for mammography. It is found that operation even at very high Talbot orders is feasible and beneficial for image quality. As long as the X-ray spectrum is matched to the visibility spectrum, the SNR continuously increases with the Talbot order for π-systems. We find that the optimal X-ray spectra and design energies are almost independent of the Talbot order and that the overall imaging performance is robust against small variations in these parameters. Discontinuous spectra, such as that from molybdenum, are less robust because the characteristic lines may coincide with minima in the visibility spectra; however, they may offer slightly better performance. We verify this hypothesis by realizing a prototype system with a mean fringe visibility of above 40% at the seventh Talbot order. With this prototype, a proof-of-principle measurement of a freshly dissected breast at reasonable compression to 4 cm is conducted with a mean glandular dose of only 3 mGy but with a high SNR.

Fast one-dimensional wave-front propagation for x-ray differential phase-contrast imaging.

Numerical wave-optical simulations of X-ray differential phase-contrast imaging using grating interferometry require the oversampling of gratings and object structures in the range of few micrometers. Consequently, fields of view of few millimeters already use large amounts of a computer's main memory to store the propagating wave front, limiting the scope of the investigations to only small-scale problems. In this study, we apply an approximation to the Fresnel-Kirchhoff diffraction theory to overcome these restrictions by dividing the two-dimensional wave front up into 1D lines, which are processed separately. The approach enables simulations with samples of clinically relevant dimensions by significantly reducing the memory footprint and the execution time and, thus, allows the qualitative comparison of different setup configurations. We analyze advantages as well as limitations and present the simulation of a virtual mammography phantom of several centimeters of size.

A Fourier-domain algorithm for total-variation regularized phase retrieval in differential X-ray phase contrast imaging.

Phase retrieval in differential X-ray phase contrast imaging involves a one dimensional integration step. In the presence of noise, standard integration methods result in image blurring and streak artifacts. This work proposes a regularized integration method which takes the availability of two dimensional data as well as the integration-specific frequency-dependent noise amplification into account. In more detail, a Fourier-domain algorithm is developed comprising a frequency-dependent minimization of the total variation orthogonal to the direction of integration. For both simulated and experimental data, the novel method yielded strong artefact reduction without increased blurring superior to the results obtained by standard integration methods or regularization techniques in the image domain.