Search

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., (© 2023 Institute of Physics and Engineering in Medicine.)

This study aims to develop a silent, fast and 3D method for T1 and proton density (PD) mapping, while generating time series of T1-weighted (T1w) images with bias-field correction. Undersampled T1w images at different effective inversion times (TIs) were acquired using the inversion recovery prepared RUFIS sequence with an interleaved k-space trajectory. Unaliased images were reconstructed by constraining the signal evolution to a temporal subspace which was learned from the signal model. Parameter maps were obtained by fitting the data to the signal model, and bias-field correction was conducted on T1w images. Accuracy and repeatability of the method was accessed in repeated experiments with phantom and volunteers. For the phantom study, T1 values obtained by the proposed method were highly consistent with values from the gold standard method, R 2 = 0.9976. Coefficients of variation (CVs) ranged from 0.09% to 0.83%. For the volunteer study, T1 values from gray and white matter regions were consistent with literature values, and peaks of gray and white matter can be clearly delineated on whole-brain T1 histograms. CVs ranged from 0.01% to 2.30%. The acoustic noise measured at the scanner isocenter was 2.6 dBA higher compared to the in-bore background. Rapid and with low acoustic noise, the proposed method is shown to produce accurate T1 and PD maps with high repeatability by reconstructing sparsely sampled T1w images at different TIs using temporal subspace. Our approach can greatly enhance patient comfort during examination and therefore increase the acceptance of the procedure.

Respiratory motion causes misalignments between positron emission tomography (PET) and magnetic resonance (MR)-derived attenuation maps (µ-maps) in addition to artifacts on both PET and MR images in simultaneous PET/MRI for organs such as liver that can experience motion of several centimeters. To address this problem, we developed an efficient MR-based attenuation correction (MRAC) method to generate phase-matched µ-maps for quiescent period PET (PET Q ) in abdominal PET/MRI. MRAC data was acquired with CIRcular Cartesian UnderSampling (CIRCUS) sampling during 100 s in free-breathing as an accelerated data acquisition strategy for phase-matched MRAC (MRAC PM-CIRCUS ). For comparison, MRAC data with raster (Default) k-space sampling was also acquired during 100 s in free-breathing (MRAC PM-Default ), and used to evaluate MRAC PM-CIRCUS as well as un-matched MRAC (MRAC UM ) that was un-gated. We purposefully oversampled the MRAC PM data to ensure we had enough information to capture all respiratory phases to make this comparison as robust as possible. The proposed MRAC PM-CIRCUS was evaluated in 17 patients with 68 Ga-DOTA-TOC PET/MRI exams, suspected of having neuroendocrine tumors or liver metastases. Effects of CIRCUS sampling for accelerating a data acquisition were evaluated by simulating the data acquisition time retrospectively in increments of 5 s. Effects of MRAC PM-CIRCUS on PET Q were evaluated using uptake differences in the liver lesions (n  =  35), compared to PET Q with MRAC PM-Default and MRAC UM . A Wilcoxon signed-rank test was performed to compare lesion uptakes between the MRAC methods. MRAC PM-CIRCUS showed higher image quality compared to MRAC PM-Default for the same acquisition times, demonstrating that a data acquisition time of 30 s was reasonable to achieve phase-matched µ-maps. Lesion update differences between MRAC PM-CIRCUS (30 s) versus MRAC PM-Default (reference, 100 s) were 0.1%  ±  1.4% (range of  -2.7% to 3.2%) and not significant (P  >  .05); while, the differences between MRAC UM versus MRAC PM-Default were 0.6%  ±  11.4% with a large variation (range of  -37% to 20%) and significant (P  <  .05). In conclusion, we demonstrated that a data acquisition of 30 s achieved phase-matched µ-maps when using specialized CIRCUS data sampling and phase-matched µ-maps improved PET Q quantification significantly.

Accurate and robust attenuation correction remains challenging in hybrid PET/MR particularly for torsos because it is difficult to segment bones, lungs and internal air in MR images. Additionally, MR suffers from susceptibility artifacts when a metallic implant is present. Recently, joint estimation (JE) of activity and attenuation based on PET data, also known as maximum likelihood reconstruction of activity and attenuation, has gained considerable interest because of (1) its promise to address the challenges in MR-based attenuation correction (MRAC), and (2) recent advances in time-of-flight (TOF) technology, which is known to be the key to the success of JE. In this paper, we implement a JE algorithm using an MR-based prior and evaluate the algorithm using whole-body PET/MR patient data, for both FDG and non-FDG tracers, acquired from GE SIGNA PET/MR scanners with TOF capability. The weight of the MR-based prior is spatially modulated, based on MR signal strength, to control the balance between MRAC and JE. Large prior weights are used in strong MR signal regions such as soft tissue and fat (i.e. MR tissue classification with a high degree of certainty) and small weights are used in low MR signal regions (i.e. MR tissue classification with a low degree of certainty). The MR-based prior is pragmatic in the sense that it is convex and does not require training or population statistics while exploiting synergies between MRAC and JE. We demonstrate the JE algorithm has the potential to improve the robustness and accuracy of MRAC by recovering the attenuation of metallic implants, internal air and some bones and by better delineating lung boundaries, not only for FDG but also for more specific non-FDG tracers such as 68 Ga-DOTATOC and 18 F-Fluoride.