A Novel method to generate on-board 4D MRI using prior 4D MRI and on-board kV projections from a conventional LINAC for target localization in liver SBRT

A review on 4D cone-beam CT (4D-CBCT) in radiation therapy: Technical advances and clinical applications

Cone-beam CT (CBCT) is the most commonly used onboard imaging technique for target localization in radiation therapy. Conventional 3D CBCT acquires x-ray cone-beam projections at multiple angles around the patient to reconstruct 3D images of the patient in the treatment room. However, despite its wide usage, 3D CBCT is limited in imaging disease sites affected by respiratory motions or other dynamic changes within the body, as it lacks time-resolved information. To overcome this limitation, 4D-CBCT was developed to incorporate a time dimension in the imaging to account for the patient's motion during the acquisitions. For example, respiration-correlated 4D-CBCT divides the breathing cycles into different phase bins and reconstructs 3D images for each phase bin, ultimately generating a complete set of 4D images. 4D-CBCT is valuable for localizing tumors in the thoracic and abdominal regions where the localization accuracy is affected by respiratory motions. This is especially important for hypofractionated stereotactic body radiation therapy (SBRT), which delivers much higher fractional doses in fewer fractions than conventional fractionated treatments. Nonetheless, 4D-CBCT does face certain limitations, including long scanning times, high imaging doses, and compromised image quality due to the necessity of acquiring sufficient x-ray projections for each respiratory phase. In order to address these challenges, numerous methods have been developed to achieve fast, low-dose, and high-quality 4D-CBCT. This paper aims to review the technical developments surrounding 4D-CBCT comprehensively. It will explore conventional algorithms and recent deep learning-based approaches, delving into their capabilities and limitations. Additionally, the paper will discuss the potential clinical applications of 4D-CBCT and outline a future roadmap, highlighting areas for further research and development. Through this exploration, the readers will better understand 4D-CBCT's capabilities and potential to enhance radiation therapy.

Evaluation of Multisource Adaptive MRI Fusion for Gross Tumor Volume Delineation of Hepatocellular Carcinoma

Purpose

Tumor delineation plays a critical role in radiotherapy for hepatocellular carcinoma (HCC) patients. The incorporation of MRI might improve the ability to correctly identify tumor boundaries and delineation consistency. In this study, we evaluated a novel Multisource Adaptive MRI Fusion (MAMF) method in HCC patients for tumor delineation.

Methods

Ten patients with HCC were included in this study retrospectively. Contrast-enhanced T1-weighted MRI at portal-venous phase (T1W_PP), contrast-enhanced T1-weighted MRI at 19-min delayed phase (T1W_DP), T2-weighted (T2W), and diffusion-weighted MRI (DWI) were acquired on a 3T MRI scanner and imported to in-house-developed MAMF software to generate synthetic MR fusion images. The original multi-contrast MR image sets were registered to planning CT by deformable image registration (DIR) using MIM. Four observers independently delineated gross tumor volumes (GTVs) on the planning CT, four original MR image sets, and the fused MRI for all patients. Tumor contrast-to-noise ratio (CNR) and Dice similarity coefficient (DSC) of the GTVs between each observer and a reference observer were measured on the six image sets. Inter-observer and inter-patient mean, SD, and coefficient of variation (CV) of the DSC were evaluated.

Results

Fused MRI showed the highest tumor CNR compared to planning CT and original MR sets in the ten patients. The mean ± SD tumor CNR was 0.72 ± 0.73, 3.66 ± 2.96, 4.13 ± 3.98, 4.10 ± 3.17, 5.25 ± 2.44, and 9.82 ± 4.19 for CT, T1W_PP, T2W, DWI, T1W_DP, and fused MRI, respectively. Fused MRI has the minimum inter-observer and inter-patient variations as compared to original MR sets and planning CT sets. GTV delineation inter-observer mean DSC across the ten patients was 0.81 ± 0.09, 0.85 ± 0.08, 0.88 ± 0.04, 0.89 ± 0.08, 0.90 ± 0.04, and 0.95 ± 0.02 for planning CT, T1W_PP, T2W, DWI, T1W_DP, and fused MRI, respectively. The patient mean inter-observer CV of DSC was 3.3%, 3.2%, 1.7%, 2.6%, 1.5%, and 0.9% for planning CT, T1W_PP, T2W, DWI, T1W_DP, and fused MRI, respectively.

Conclusion

The results demonstrated that the fused MRI generated using the MAMF method can enhance tumor CNR and improve inter-observer consistency of GTV delineation in HCC as compared to planning CT and four commonly used MR image sets (T1W_PP, T1W_DP, T2W, and DWI). The MAMF method holds great promise in MRI applications in HCC radiotherapy treatment planning.

A Dual-supervised Deformation Estimation Model (DDEM) for constructing ultra-quality 4D-MRI based on a commercial low-quality 4D-MRI for liver cancer radiation therapy

BACKGROUND

Most available 4D-MRI techniques are limited by insufficient image quality and long acquisition times or require specially designed sequences or hardware that are not available in the clinic. These limitations have greatly hindered the clinical implementation of 4D-MRI.

PURPOSE

This study aims to develop a fast ultra-quality (UQ) 4D-MRI reconstruction method using a commercially available 4D-MRI sequence and dual-supervised deformation estimation model (DDEM).

METHODS

Thirty-nine patients receiving radiotherapy for liver tumors were included. Each patient was scanned using a TWIST-VIBE MRI sequence to acquire 4D-MR images. They also received 3D T1-/T2-weighted MRI scans as prior images and UQ 4D-MRI at any instant was considered a deformation of them. A DDEM was developed to obtain a 4D deformable vector field (DVF) from 4D-MRI data, and the prior images were deformed using this 4D-DVF to generate UQ 4D-MR images. The registration accuracies of the DDEM, VoxelMorph (normalized cross-correlation (NCC) supervised), VoxelMorph (end-to-end point error (EPE) supervised), and the parametric total variation (pTV) algorithm were compared. Tumor motion on UQ 4D-MRI was evaluated quantitatively using region-of-interest (ROI) tracking errors, while image quality was evaluated using the contrast-to-noise ratio (CNR), lung-liver edge sharpness, and perceptual blur metric (PBM).

RESULTS

The registration accuracy of the DDEM was significantly better than those of VoxelMorph (NCC supervised), VoxelMorph (EPE supervised) and the pTV algorithm (all, p < 0.001), with an inference time of 69.3 ± 5.9 ms. UQ 4D-MRI yielded ROI tracking errors of 0.79 ± 0.65, 0.50 ± 0.55, and 0.51 ± 0.58 mm in the superior-inferior, anterior-posterior, and mid-lateral directions, respectively. From the original 4D-MRI to UQ 4D-MRI, the CNR increased from 7.25 ± 4.89 to 18.86 ± 15.81; the lung-liver edge full-width-at-half-maximum decreased from 8.22 ± 3.17 to 3.65 ± 1.66 mm in the in-plane direction and from 8.79 ± 2.78 to 5.04 ± 1.67 mm in the cross-plane direction, and the PBM decreased from 0.68 ± 0.07 to 0.38 ± 0.01.

CONCLUSION

This novel DDEM method successfully generated UQ 4D-MR images based on a commercial 4D-MRI sequence. It shows great promise for improving liver tumor motion management during radiation therapy. This article is protected by copyright. All rights reserved.

Probabilistic 4D predictive model from in-room surrogates using conditional generative networks for image-guided radiotherapy

Shape and location organ variability induced by respiration constitutes one of the main challenges during dose delivery in radiotherapy. Providing up-to-date volumetric information during treatment can improve tumor tracking, thereby increasing treatment efficiency and reducing damage to healthy tissue. We propose a novel probabilistic model to address the problem of volumetric estimation with scalable predictive horizon from image-based surrogates during radiotherapy treatments, thus enabling out-of-plane tracking of targets. This problem is formulated as a conditional learning task, where the predictive variables are the 2D surrogate images and a pre-operative static 3D volume. The model learns a distribution of realistic motion fields over a population dataset. Simultaneously, a seq-2-seq inspired temporal mechanism acts over the surrogate images yielding extrapolated-in-time representations. The phase-specific motion distributions are associated with the predicted temporal representations, allowing the recovery of dense organ deformation in multiple times. Due to its generative nature, this model enables uncertainty estimations by sampling the latent space multiple times. Furthermore, it can be readily personalized to a new subject via fine-tuning, and does not require inter-subject correspondences. The proposed model was evaluated on free-breathing 4D MRI and ultrasound datasets from 25 healthy volunteers, as well as on 11 cancer patients. A navigator-based data augmentation strategy was used during the slice reordering process to increase model robustness against inter-cycle variability. The patient data was used as a hold-out test set. Our approach yields volumetric prediction from image surrogates with a mean error of 1.67 ± 1.68 mm and 2.17 ± 0.82 mm in unseen cases of the patient MRI and US datasets, respectively. Moreover, model personalization yields a mean landmark error of 1.4 ± 1.1 mm compared to ground truth annotations in the volunteer MRI dataset, with statistically significant improvements over state-of-the-art.

A Monte Carlo study to investigate the feasibility of an on-board SPECT/spectral-CT/CBCT imager for medical linear accelerator

PURPOSE

The on-board flat-panel cone-beam computed tomography (CBCT) lacks molecular/functional information for current online image-guided radiation therapy (IGRT). It might not be adequate for adaptive radiation therapy (ART), particularly for biologically guided tumor delineation and targeting which might be shifted and/or distorted during the course of RT. A linear accelerator (Linac) gantry-mounted on-board imager (OBI) was proposed using a single photon counting detector (PCD) panel to achieve single photon emission computed tomography (SPECT), energy-resolved spectral CT, and conventional CBCT triple on-board imaging, which might facilitate online ART with an addition of volumetric molecular/functional imaging information.

METHODS

The system was designed and evaluated in the GATE Monte Carlo platform. The OBI system including a kV-beam source and a pixelated cadmium zinc telluride (CZT) detector panel mounted on a medical Linac orthogonally to the MV beam direction was designed to obtain online CBCT, spectral CT, and SPECT tri-modal imaging of patients in the treatment room. The spatial resolutions of the OBI system were determined by imaging simulated phantoms. The CBCT imaging was evaluated by a simulated contrast phantom. A PMMA phantom containing gadolinium was imaged to demonstrate quantitative imaging of spectral-CT/CBCT of the system. The capability of tri-modal imaging of the OBI was demonstrated using three different spectral CT imaging methods to differentiate gadolinium, gold, calcium within simulated PMMA and the SPECT to image radioactive ^99m Tc distribution. The dual-isotope SPECT imaging of the system was also evaluated by imaging a phantom containing ^99m Tc and ¹²³ I. The radiotherapy-related parameters of iodine contrast fraction and virtual non-contrast (VNC) tissue electron density in the Kidney1 inserts of a simulated phantom were decomposed using the Bayesian eigentissue decomposition method for contrast-enhanced CBCT/spectral-CT of the OBI in a single scan.

RESULTS

The spatial resolutions of CBCT and SPECT of the OBI were determined to be 15.1 lp/cm at 10% MTF and 4.8-12 mm for radii of rotation of 10-40 cm, respectively. In CBCT image of the contrast phantom, most of the soft-tissue inserts were visible with sufficient spatial structure details. As compared to the CBCT image of gadolinium, the spectral CT image provided higher image contrasts. Calcium, gadolinium, and gold were separated well by using the spectral CT material imaging methods. The reconstructed distribution of ^99m Tc agreed with the spatial position within the phantom. The two isotopes were separated from each other in dual-isotope SPECT imaging of the OBI. The iodine fractions and the VNC electron densities were estimated in the iodine-enhanced Kidney1 tissue inserts with reasonable RMS errors. The main procedures of the tri-modal imaging guided online ART workflow were presented with new functional features included.

CONCLUSIONS

Using a single photon counting CZT detector panel, an on-board SPECT, spectral CT, and CBCT tri-modal imaging could be realized in Linacs. With the added online molecular/functional imaging obtained from the new OBI for the online ART proposed, the accuracy of radiation treatment delivery could be further improved.

A modern review of the uncertainties in volumetric imaging of respiratory-induced target motion in lung radiotherapy

Radiotherapy has become a critical component for the treatment of all stages and types of lung cancer, often times being the primary gateway to a cure. However, given that radiation can cause harmful side effects depending on how much surrounding healthy tissue is exposed, treatment of the lung can be particularly challenging due to the presence of moving targets. Careful implementation of every step in the radiotherapy process is absolutely integral for attaining optimal clinical outcomes. With the advent and now widespread use of stereotactic body radiation therapy (SBRT), where extremely large doses are delivered, accurate, and precise dose targeting is especially vital to achieve an optimal risk to benefit ratio. This has largely become possible due to the rapid development of image-guided technology. Although imaging is critical to the success of radiotherapy, it can often be plagued with uncertainties due to respiratory-induced target motion. There has and continues to be an immense research effort aimed at acknowledging and addressing these uncertainties to further our abilities to more precisely target radiation treatment. Thus, the goal of this article is to provide a detailed review of the prevailing uncertainties that remain to be investigated across the different imaging modalities, as well as to highlight the more modern solutions to imaging motion and their role in addressing the current challenges.

Continuous generation of volumetric images during stereotactic body radiation therapy using periodic kV imaging and an external respiratory surrogate

We present a technique for continuous generation of volumetric images during SBRT using periodic kV imaging and an external respiratory surrogate signal to drive a patient-specific PCA motion model. Using the on-board imager, kV radiographs are acquired every 3 s and used to fit the parameters of a motion model so that it matches observed changes in internal patient anatomy. A multi-dimensional correlation model is established between the motion model parameters and the external surrogate position and velocity, enabling volumetric image reconstruction between kV imaging time points. Performance of the algorithm was evaluated using 10 realistic eXtended CArdiac-Torso (XCAT) digital phantoms including 3D anatomical respiratory deformation programmed with 3D tumor positions measured with orthogonal kV imaging of implanted fiducial gold markers. The clinically measured ground truth 3D tumor positions provided a dataset with realistic breathing irregularities, and the combination of periodic on-board kV imaging with recorded external respiratory surrogate signal was used for correlation modeling to account for any changes in internal-external correlation. The three-dimensional tumor positions are reconstructed with an average root mean square error (RMSE) of 1.47 mm, and an average 95th percentile 3D positional error of 2.80 mm compared with the clinically measured ground truth 3D tumor positions. This technique enables continuous 3D anatomical image generation based on periodic kV imaging of internal anatomy without the additional dose of continuous kV imaging. The 3D anatomical images produced using this method can be used for treatment verification and delivered dose computation in the presence of irregular respiratory motion.