Review of Real-Time 3-Dimensional Image Guided Radiation Therapy on Standard-Equipped Cancer Radiation Therapy Systems: Are We at the Tipping Point for the Era of Real-Time Radiation Therapy?

Reduction of Postradiation Therapy Urinary Toxicity Via Intrafractional Megavoltage-Kilovoltage Prostate Location Monitoring

PURPOSE

We hypothesized that an in-house developed system using megavoltage and kilovoltage image guidance (MKIG) to ensure correct prostate positioning during stereotactic body radiation therapy (SBRT) could potentially avoid unwanted doses to nontarget tissues, leading to reduced toxicities.

METHODS AND MATERIALS

We built a 3-dimensional MKIG platform that accurately tracks prostate implanted fiducials in real time and clinically translated the system to replace a commercial approach, intrafraction motion review (IMR), which only tracks fiducials in the 2-dimensional kilovoltage views. From 2017 to 2019, 150 patients with prostate cancer were treated with SBRT and monitored using MKIG. The motion trace of the fiducials alerts therapists to interrupt and reposition the prostate when displacement exceeds a 1.5 mm threshold. A comparison cohort of 121 patients was treated with the same dose regimen and treatment technique but managed by IMR. Statistics of intrafractional patient shifts and delivery time were collected to evaluate the workflow efficacy. The incidence of grade ≥2 urinary toxicities was analyzed to assess clinical complications. The median follow-up time was 3.7 years (0.2-8.2 years).

RESULTS

MKIG treatments had more treatment shifts (1.09 vs 0.28) and a longer average delivery time per fraction (579 ± 205 seconds vs 357 ± 117 seconds) than IMR treatments. Three-quarters (75%) of shifts resulting from MKIG were ≤3 mm, versus 51% in IMR, indicating that MKIG detected and corrected smaller deviations. The incidence of grade ≥2 urinary toxicity was lower in the MKIG than the IMR cohort: 10.7% versus 19.8% (P = .047). On multivariate analysis of late urinary toxicity, only high (>7) preradiation therapy international prostate symptom score (P < .043) and the use of MKIG were selected (P < .029).

CONCLUSIONS

Automated and quantitative MKIG introduced minimal workflow impact and was superior to IMR in localizing the prostate during SBRT, which correlated with a clinically significant reduction in late urinary toxicity. Further clinical testing using randomized trials will be required to validate the impact on outcomes.

Full-dose whole-body PET synthesis from low-dose PET using high-efficiency denoising diffusion probabilistic model: PET consistency model

PURPOSE

Positron Emission Tomography (PET) has been a commonly used imaging modality in broad clinical applications. One of the most important tradeoffs in PET imaging is between image quality and radiation dose: high image quality comes with high radiation exposure. Improving image quality is desirable for all clinical applications while minimizing radiation exposure is needed to reduce risk to patients.

METHODS

We introduce PET Consistency Model (PET-CM), an efficient diffusion-based method for generating high-quality full-dose PET images from low-dose PET images. It employs a two-step process, adding Gaussian noise to full-dose PET images in the forward diffusion, and then denoising them using a PET Shifted-window Vision Transformer (PET-VIT) network in the reverse diffusion. The PET-VIT network learns a consistency function that enables direct denoising of Gaussian noise into clean full-dose PET images. PET-CM achieves state-of-the-art image quality while requiring significantly less computation time than other methods. Evaluation with normalized mean absolute error (NMAE), peak signal-to-noise ratio (PSNR), multi-scale structure similarity index (SSIM), normalized cross-correlation (NCC), and clinical evaluation including Human Ranking Score (HRS) and Standardized Uptake Value (SUV) Error analysis shows its superiority in synthesizing full-dose PET images from low-dose inputs.

RESULTS

In experiments comparing eighth-dose to full-dose images, PET-CM demonstrated impressive performance with NMAE of 1.278 ± 0.122%, PSNR of 33.783 ± 0.824 dB, SSIM of 0.964 ± 0.009, NCC of 0.968 ± 0.011, HRS of 4.543, and SUV Error of 0.255 ± 0.318%, with an average generation time of 62 s per patient. This is a significant improvement compared to the state-of-the-art diffusion-based model with PET-CM reaching this result 12× faster. Similarly, in the quarter-dose to full-dose image experiments, PET-CM delivered competitive outcomes, achieving an NMAE of 0.973 ± 0.066%, PSNR of 36.172 ± 0.801 dB, SSIM of 0.984 ± 0.004, NCC of 0.990 ± 0.005, HRS of 4.428, and SUV Error of 0.151 ± 0.192% using the same generation process, which underlining its high quantitative and clinical precision in both denoising scenario.

CONCLUSIONS

We propose PET-CM, the first efficient diffusion-model-based method, for estimating full-dose PET images from low-dose images. PET-CM provides comparable quality to the state-of-the-art diffusion model with higher efficiency. By utilizing this approach, it becomes possible to maintain high-quality PET images suitable for clinical use while mitigating the risks associated with radiation. The code is availble at https://github.com/shaoyanpan/Full-dose-Whole-body-PET-Synthesis-from-Low-dose-PET-Using-Consistency-Model.

CBCT-DRRs superior to CT-DRRs for target-tracking applications for pancreatic SBRT

Objective.In current radiograph-based intra-fraction markerless target-tracking, digitally reconstructed radiographs (DRRs) from planning CTs (CT-DRRs) are often used to train deep learning models that extract information from the intra-fraction radiographs acquired during treatment. Traditional DRR algorithms were designed for patient alignment (i.e.bone matching) and may not replicate the radiographic image quality of intra-fraction radiographs at treatment. Hypothetically, generating DRRs from pre-treatment Cone-Beam CTs (CBCT-DRRs) with DRR algorithms incorporating physical modelling of on-board-imagers (OBIs) could improve the similarity between intra-fraction radiographs and DRRs by eliminating inter-fraction variation and reducing image-quality mismatches between radiographs and DRRs. In this study, we test the two hypotheses that intra-fraction radiographs are more similar to CBCT-DRRs than CT-DRRs, and that intra-fraction radiographs are more similar to DRRs from algorithms incorporating physical models of OBI components than DRRs from algorithms omitting these models.Approach.DRRs were generated from CBCT and CT image sets collected from 20 patients undergoing pancreas stereotactic body radiotherapy. CBCT-DRRs and CT-DRRs were generated replicating the treatment position of patients and the OBI geometry during intra-fraction radiograph acquisition. To investigate whether the modelling of physical OBI components influenced radiograph-DRR similarity, four DRR algorithms were applied for the generation of CBCT-DRRs and CT-DRRs, incorporating and omitting different combinations of OBI component models. The four DRR algorithms were: a traditional DRR algorithm, a DRR algorithm with source-spectrum modelling, a DRR algorithm with source-spectrum and detector modelling, and a DRR algorithm with source-spectrum, detector and patient material modelling. Similarity between radiographs and matched DRRs was quantified using Pearson's correlation and Czekanowski's index, calculated on a per-image basis. Distributions of correlations and indexes were compared to test each of the hypotheses. Distribution differences were determined to be statistically significant when Wilcoxon's signed rank test and the Kolmogorov-Smirnov two sample test returnedp≤ 0.05 for both tests.Main results.Intra-fraction radiographs were more similar to CBCT-DRRs than CT-DRRs for both metrics across all algorithms, with allp≤ 0.007. Source-spectrum modelling improved radiograph-DRR similarity for both metrics, with allp< 10-6. OBI detector modelling and patient material modelling did not influence radiograph-DRR similarity for either metric.Significance.Generating DRRs from pre-treatment CBCT-DRRs is feasible, and incorporating CBCT-DRRs into markerless target-tracking methods may promote improved target-tracking accuracies. Incorporating source-spectrum modelling into a treatment planning system's DRR algorithms may reinforce the safe treatment of cancer patients by aiding in patient alignment.

Biomedical advances and future prospects of high-precision three-dimensional radiotherapy and four-dimensional radiotherapy

Biomedical advances of external-beam radiotherapy (EBRT) with improvements in physical accuracy are reviewed. High-precision (±1 mm) three-dimensional radiotherapy (3DRT) can utilize respective therapeutic open doors in the tumor control probability curve and in the normal tissue complication probability curve instead of the one single therapeutic window in two-dimensional EBRT. High-precision 3DRT achieved higher tumor control and probable survival rates for patients with small peripheral lung and liver cancers. Four-dimensional radiotherapy (4DRT), which can reduce uncertainties in 3DRT due to organ motion by real-time (every 0.1-1 s) tumor-tracking and immediate (0.1-1 s) irradiation, have achieved reduced adverse effects for prostate and pancreatic tumors near the digestive tract and with similar or better tumor control. Particle beam therapy improved tumor control and probable survival for patients with large liver tumors. The clinical outcomes of locally advanced or multiple tumors located near serial-type organs can theoretically be improved further by integrating the 4DRT concept with particle beams.

Patient specific prior cross attention for kV decomposition in paraspinal motion tracking

BACKGROUND

X-ray image quality is critical for accurate intrafraction motion tracking in radiation therapy.

PURPOSE

This study aims to develop a deep-learning algorithm to improve kV image contrast by decomposing the image into bony and soft tissue components. In particular, we designed a priori attention mechanism in the neural network framework for optimal decomposition. We show that a patient-specific prior cross-attention (PCAT) mechanism can boost the performance of kV image decomposition. We demonstrate its use in paraspinal SBRT motion tracking with online kV imaging.

METHODS

Online 2D kV projections were acquired during paraspinal SBRT for patient motion monitoring. The patient-specific prior images were generated by randomly shifting and rotating spine-only DRR created from the setup CBCT, simulating potential motions. The latent features of the prior images were incorporated into the PCAT using multi-head cross attention. The neural network aimed to learn to selectively amplify the transmission of the projection image features that correlate with features of the priori. The PCAT network structure consisted of (1) a dual-branch generator that separates the spine and soft tissue component of the kV projection image and (2) a dual-function discriminator (DFD) that provides the realness score of the predicted images. For supervision, we used a loss combining mean absolute error loss, discriminator loss, perceptual loss, total variation, and mean squared error loss for soft tissues. The proposed PCAT approach was benchmarked against previous work using the ResNet generative adversarial network (ResNetGAN) without prior information.

RESULTS

The trained PCAT had improved performance in effectively retaining and preserving the spine structure and texture information while suppressing the soft tissues from the kV projection images. The decomposed spine-only x-ray images had the submillimeter matching accuracy at all beam angles. The decomposed spine-only x-ray significantly reduced the maximum errors to 0.44 mm (<2 pixels) in comparison to 0.92 mm (∼4 pixels) of ResNetGAN. The PCAT decomposed spine images also had higher PSNR and SSIM (p-value < 0.001).

CONCLUSION

The PCAT selectively learned the important latent features by incorporating the patient-specific prior knowledge into the deep learning algorithm, significantly improving the robustness of the kV projection image decomposition, and leading to improved motion tracking accuracy in paraspinal SBRT.

Identification of reproducible radiomic features from on-board volumetric images: A multi-institutional phantom study

BACKGROUND

Radiomics analysis using on-board volumetric images has attracted research attention as a method for predicting prognosis during treatment; however, the lack of standardization is still one of the main concerns.

PURPOSE

This study investigated the factors that influence the reproducibility of radiomic features extracted from on-board volumetric images using an anthropomorphic radiomics phantom. Furthermore, a phantom experiment was conducted with different treatment machines from multiple institutions as external validation to identify reproducible radiomic features.

METHODS

The phantom was designed to be 35 × 20 × 20 cm with eight types of heterogeneous spheres (⌀ = 1, 2, and 3 cm). On-board volumetric images were acquired using 15 treatment machines from eight institutions. Of these, kilovoltage cone-beam computed tomography (kV-CBCT) image data acquired from four treatment machines at one institution were used as an internal evaluation dataset to explore the reproducibility of radiomic features. The remaining image data, including kV-CBCT, megavoltage-CBCT (MV-CBCT), and megavoltage computed tomography (MV-CT) provided by seven different institutions (11 treatment machines), were used as an external validation dataset. A total of 1,302 radiomic features, including 18 first-order, 75 texture, 465 (i.e., 93 × 5) Laplacian of Gaussian (LoG) filter-based, and 744 (i.e., 93 × 8) wavelet filter-based features, were extracted within the spheres. The intraclass correlation coefficient (ICC) was calculated to explore feature repeatability and reproducibility using an internal evaluation dataset. Subsequently, the coefficient of variation (COV) was calculated to validate the feature variability of external institutions. An absolute ICC exceeding 0.85 or COV under 5% was considered indicative of a highly reproducible feature.

RESULTS

For internal evaluation, ICC analysis showed that the median percentage of radiomic features with high repeatability was 95.2%. The ICC analysis indicated that the median percentages of highly reproducible features for inter-tube current, reconstruction algorithm, and treatment machine were decreased by 20.8%, 29.2%, and 33.3%, respectively. For external validation, the COV analysis showed that the median percentage of reproducible features was 31.5%. A total of 16 features, including nine LoG filter-based and seven wavelet filter-based features, were indicated as highly reproducible features. The gray-level run-length matrix (GLRLM) was classified as containing the most frequent features (N = 8), followed by the gray-level dependence matrix (N = 7) and gray-level co-occurrence matrix (N = 1) features.

CONCLUSIONS

We developed the standard phantom for radiomics analysis of kV-CBCT, MV-CBCT, and MV-CT images. With this phantom, we revealed that the differences in the treatment machine and image reconstruction algorithm reduce the reproducibility of radiomic features from on-board volumetric images. Specifically, the most reproducible features for external validation were LoG or wavelet filter-based GLRLM features. However, the acceptability of the identified features should be examined in advance at each institution before applying the findings to prognosis prediction.

Real-time motion monitoring using orthogonal cine MRI during MR-guided adaptive radiation therapy for abdominal tumors on 1.5T MR-Linac

BACKGROUND

Real-time motion monitoring (RTMM) is necessary for accurate motion management of intrafraction motions during radiation therapy (RT).

PURPOSE

Building upon a previous study, this work develops and tests an improved RTMM technique based on real-time orthogonal cine magnetic resonance imaging (MRI) acquired during magnetic resonance-guided adaptive RT (MRgART) for abdominal tumors on MR-Linac.

METHODS

A motion monitoring research package (MMRP) was developed and tested for RTMM based on template rigid registration between beam-on real-time orthogonal cine MRI and pre-beam daily reference 3D-MRI (baseline). The MRI data acquired under free-breathing during the routine MRgART on a 1.5T MR-Linac for 18 patients with abdominal malignancies of 8 liver, 4 adrenal glands (renal fossa), and 6 pancreas cases were used to evaluate the MMRP package. For each patient, a 3D mid-position image derived from an in-house daily 4D-MRI was used to define a target mask or a surrogate sub-region encompassing the target. Additionally, an exploratory case reviewed for an MRI dataset of a healthy volunteer acquired under both free-breathing and deep inspiration breath-hold (DIBH) was used to test how effectively the RTMM using the MMRP can address through-plane motion (TPM). For all cases, the 2D T2/T1-weighted cine MRIs were captured with a temporal resolution of 200 ms interleaved between coronal and sagittal orientations. Manually delineated contours on the cine frames were used as the ground-truth motion. Common visible vessels and segments of target boundaries in proximity to the target were used as anatomical landmarks for reproducible delineations on both the 3D and the cine MRI images. Standard deviation of the error (SDE) between the ground-truth and the measured target motion from the MMRP package were analyzed to evaluate the RTMM accuracy. The maximum target motion (MTM) was measured on the 4D-MRI for all cases during free-breathing.

RESULTS

The mean (range) centroid motions for the 13 abdominal tumor cases were 7.69 (4.71-11.15), 1.73 (0.81-3.05), and 2.71 (1.45-3.93) mm with an overall accuracy of <2 mm in the superior-inferior (SI), the left-right (LR), and the anterior-posterior (AP) directions, respectively. The mean (range) of the MTM from the 4D-MRI was 7.38 (2-11) mm in the SI direction, smaller than the monitored motion of centroid, demonstrating the importance of the real-time motion capture. For the remaining patient cases, the ground-truth delineation was challenging under free-breathing due to the target deformation and the large TPM in the AP direction, the implant-induced image artifacts, and/or the suboptimal image plane selection. These cases were evaluated based on visual assessment. For the healthy volunteer, the TPM of the target was significant under free-breathing which degraded the RTMM accuracy. RTMM accuracy of <2 mm was achieved under DIBH, indicating DIBH is an effective method to address large TPM.

CONCLUSIONS

We have successfully developed and tested the use of a template-based registration method for an accurate RTMM of abdominal targets during MRgART on a 1.5T MR-Linac without using injected contrast agents or radio-opaque implants. DIBH may be used to effectively reduce or eliminate TPM of abdominal targets during RTMM.

Markerless motion tracking with simultaneous MV and kV imaging in spine SBRT treatment-a feasibility study

Objective. Motion tracking with simultaneous MV-kV imaging has distinct advantages over single kV systems. This research is a feasibility study of utilizing this technique for spine stereotactic body radiotherapy (SBRT) through phantom and patient studies.Approach. A clinical spine SBRT plan was developed using 6xFFF beams and nine sliding-window IMRT fields. The plan was delivered to a chest phantom on a linear accelerator. Simultaneous MV-kV image pairs were acquired during beam delivery. KV images were triggered at predefined intervals, and synthetic MV images showing enlarged MLC apertures were created by combining multiple raw MV frames with corrections for scattering and intensity variation. Digitally reconstructed radiograph (DRR) templates were generated using high-resolution CBCT reconstructions (isotropic voxel size (0.243 mm)3) as the reference for 2D-2D matching. 3D shifts were calculated from triangulation of kV-to-DRR and MV-to-DRR registrations. To evaluate tracking accuracy, detected shifts were compared to known phantom shifts as introduced before treatment. The patient study included a T-spine patient and an L-spine patient. Patient datasets were retrospectively analyzed to demonstrate the performance in clinical settings.Main results. The treatment plan was delivered to the phantom in five scenarios: no shift, 2 mm shift in one of the longitudinal, lateral and vertical directions, and 2 mm shift in all the three directions. The calculated 3D shifts agreed well with the actual couch shifts, and overall, the uncertainty of 3D detection is estimated to be 0.3 mm. The patient study revealed that with clinical patient image quality, the calculated 3D motion agreed with the post-treatment cone beam CT. It is feasible to automate both kV-to-DRR and MV-to-DRR registrations using a mutual information-based method, and the difference from manual registration is generally less than 0.3 mm.Significance. The MV-kV imaging-based markerless motion tracking technique was validated through a feasibility study. It is a step forward toward effective motion tracking and accurate delivery for spinal SBRT.

Explainability and controllability of patient-specific deep learning with attention-based augmentation for markerless image-guided radiotherapy

BACKGROUND

We reported the concept of patient-specific deep learning (DL) for real-time markerless tumor segmentation in image-guided radiotherapy (IGRT). The method was aimed to control the attention of convolutional neural networks (CNNs) by artificial differences in co-occurrence probability (CoOCP) in training datasets, that is, focusing CNN attention on soft tissues while ignoring bones. However, the effectiveness of this attention-based data augmentation has not been confirmed by explainable techniques. Furthermore, compared to reasonable ground truths, the feasibility of tumor segmentation in clinical kilovolt (kV) X-ray fluoroscopic (XF) images has not been confirmed.

PURPOSE

The first aim of this paper was to present evidence that the proposed method provides an explanation and control of DL behavior. The second purpose was to validate the real-time lung tumor segmentation in clinical kV XF images for IGRT.

METHODS

This retrospective study included 10 patients with lung cancer. Patient-specific and XF angle-specific image pairs comprising digitally reconstructed radiographs (DRRs) and projected-clinical-target-volume (pCTV) images were calculated from four-dimensional computer tomographic data and treatment planning information. The training datasets were primarily augmented by random overlay (RO) and noise injection (NI): RO aims to differentiate positional CoOCP in soft tissues and bones, and NI aims to make a difference in the frequency of occurrence of local and global image features. The CNNs for each patient-and-angle were automatically optimized in the DL training stage to transform the training DRRs into pCTV images. In the inference stage, the trained CNNs transformed the test XF images into pCTV images, thus identifying target positions and shapes.

RESULTS

The visual analysis of DL attention heatmaps for a test image demonstrated that our method focused CNN attention on soft tissue and global image features rather than bones and local features. The processing time for each patient-and-angle-specific dataset in the training stage was ∼30 min, whereas that in the inference stage was 8 ms/frame. The estimated three-dimensional 95 percentile tracking error, Jaccard index, and Hausdorff distance for 10 patients were 1.3-3.9 mm, 0.85-0.94, and 0.6-4.9 mm, respectively.

CONCLUSIONS

The proposed attention-based data augmentation with both RO and NI made the CNN behavior more explainable and more controllable. The results obtained demonstrated the feasibility of real-time markerless lung tumor segmentation in kV XF images for IGRT.

Dosimetric impact of rotational set-up errors in high-risk prostate cancer

Introduction: Cone-beam computed tomography (CBCT) provides an excellent solution to quantitative assessment and correction of patient set-up errors during radiotherapy. However, most linear accelerators are equipped with conventional therapy tables that can be moved in three translational directions and perform only yaw rotation. Uncorrected roll and pitch result in rotational set-up errors, particularly when the distance from the isocenter to the target border is large. The aim of this study was to investigate the impact of rotational errors on the dose delivered to the clinical target volume (CTV), the planning target volume (PTV) and organs at risk (OAR).

Material and methods: 30 patients with prostate cancer treated with VMAT technique had daily CBCT scans (840 CBCTs in total) prior to treatment delivery. The rotational errors remaining after on-line correction were retrospectively analysed. The sum plans simulating the dose distribution during the treatment course were calculated for selected patients with significant rotational errors.

Results: The dose delivered to the prostate bed CTV reported in the sum plan was not lower than in the original plan for all selected patients. For four patients from the selected group, the D98% for prostate bed PTV was less than 95%. The V47.88Gy for pelvic lymph nodes PTV was less than 98% for two of the selected patients.

Conclusions: The analysis of the dosimetric parameters showed that the impact of uncorrected rotations is not clinically significant in terms of the dose delivered to OAR and the dose coverage of CTV. However, the PTV dose coverage is correlated with distance away from the isocenter and is smaller than planned.

Chemical Overview of Gel Dosimetry Systems: A Comprehensive Review

Advances in radiotherapy technology during the last 25 years have significantly improved both dose conformation to tumors and the preservation of healthy tissues, achieving almost real-time feedback by means of high-precision treatments and theranostics. Owing to this, developing high-performance systems capable of coping with the challenging requirements of modern ionizing radiation is a key issue to overcome the limitations of traditional dosimeters. In this regard, a deep understanding of the physicochemical basis of gel dosimetry, as one of the most promising tools for the evaluation of 3D high-spatial-resolution dose distributions, represents the starting point for developing new and innovative systems. This review aims to contribute thorough descriptions of the chemical processes and interactions that condition gel dosimetry outputs, often phenomenologically addressed, and particularly formulations reported since 2017.

Using 4D dose accumulation to calculate organ-at-risk dose deviations from motion-synchronized liver and lung tomotherapy treatments

Tracking systems such as Radixact Synchrony change the planned delivery of radiation during treatment to follow the target. This is typically achieved without considering the location changes of organs at risk (OARs). The goal of this work was to develop a novel 4D dose accumulation framework to quantify OAR dose deviations due to the motion and tracked treatment. The framework obtains deformation information and the target motion pattern from a four-dimensional computed tomography dataset. The helical tomotherapy treatment plan is split into 10 plans and motion correction is applied separately to the jaw pattern and multi-leaf collimator (MLC) sinogram for each phase based on the location of the target in each phase. Deformable image registration (DIR) is calculated from each phase to the references phase using a commercial algorithm, and doses are accumulated according to the DIR. The effect of motion synchronization on OAR dose was analyzed for five lung and five liver subjects by comparing planned versus synchrony-accumulated dose. The motion was compensated by an average of 1.6 cm of jaw sway and by an average of 5.7% of leaf openings modified, indicating that most of the motion compensation was from jaw sway and not MLC changes. OAR dose deviations as large as 19 Gy were observed, and for all 10 cases, dose deviations greater than 7 Gy were observed. Target dose remained relatively constant (D95% within 3 Gy), confirming that motion-synchronization achieved the goal of maintaining target dose. Dose deviations provided by the framework can be leveraged during the treatment planning process by identifying cases where OAR doses may change significantly from their planned values with respect to the critical constraints. The framework is specific to synchronized helical tomotherapy treatments, but the OAR dose deviations apply to any real-time tracking technique that does not consider location changes of OARs.

An abdominal phantom with anthropomorphic organ motion and multimodal imaging contrast for MR-guided radiotherapy

Purpose. Improvements in image-guided radiotherapy (IGRT) enable accurate and precise treatment of moving tumors in the abdomen while simultaneously sparing healthy tissue. However, the lack of validation tools for newly developed MR-guided radiotherapy hybrid devices such as the MR-Linac is an open issue. This study presents a custom developed abdominal phantom with respiratory organ motion and multimodal imaging contrast to perform end-to-end tests for IGRT treatment planning scenarios.Methods. The abdominal phantom contains deformable and anatomically shaped liver and kidney models made of Ni-DTPA and KCl-doped agarose mixtures that can be reproducibly positioned within the phantom. Organ models are wrapped in foil to avoid ion exchange with the surrounding agarose and to provide stable T1 and T2 relaxation times as well as HU numbers. Breathing motion is realized by a diaphragm connected to an actuator that is hydraulically controlled via a programmable logic controller. With this system, artificial and patient-specific breathing patterns can be carried out. In 1.5 T magnetic resonance imaging (MRI), diaphragm, liver and kidney motion was measured and compared to the breathing motion of a healthy male volunteer for different breathing amplitudes including shallow, normal and deep breathing.Results. The constructed abdominal phantom demonstrated organ-equivalent intensity values in CT as well as in MRI. T1-weighted (T1w) and T2-weighted (T2w) relaxation times for 1.5 T and CT numbers were 552.9 ms, 48.2 ms and 48.8 HU (liver) as well as 950.42 ms, 79 ms and 28.2 HU (kidney), respectively. These values were stable for more than six months. Extracted breathing motion from a healthy volunteer revealed a liver to diaphragm motion ratio (LDMR) of 64.4% and a kidney to diaphragm motion ratio (KDMR) of 30.7%. Well-comparable values were obtained for the phantom (LDMR: 65.5%, KDMR: 27.5%).Conclusions. The abdominal phantom demonstrated anthropomorphic T1 and T2 relaxation times as well as HU numbers and physiological motion pattern in MRI and CT. This allows for wide use in the validation of IGRT including MRgRT.

Artificial intelligence in image-guided radiotherapy: a review of treatment target localization

Modern conformal beam delivery techniques require image-guidance to ensure the prescribed dose to be delivered as planned. Recent advances in artificial intelligence (AI) have greatly augmented our ability to accurately localize the treatment target while sparing the normal tissues. In this paper, we review the applications of AI-based algorithms in image-guided radiotherapy (IGRT), and discuss the indications of these applications to the future of clinical practice of radiotherapy. The benefits, limitations and some important trends in research and development of the AI-based IGRT techniques are also discussed. AI-based IGRT techniques have the potential to monitor tumor motion, reduce treatment uncertainty and improve treatment precision. Particularly, these techniques also allow more healthy tissue to be spared while keeping tumor coverage the same or even better.

Six degrees of freedom dynamic motion-including dose reconstruction in a commercial treatment planning system

PURPOSE

Intrafractional motion during radiotherapy delivery can deteriorate the delivered dose. Dynamic rotational motion of up to 38 degrees has been reported during prostate cancer radiotherapy, but methods to determine the dosimetric consequences of such rotations are lacking. Here, we create and experimentally validate a dose reconstruction method that accounts for dynamic rotations and translations in a commercial treatment planning system (TPS). Interplay effects are quantified by comparing dose reconstructions with dynamic and constant rotations.

METHODS

The dose reconstruction accumulates the dose in points of interest while the points are moved in six degrees of freedom (6DoF) in a precalculated time-resolved four-dimensional (4D) dose matrix to emulate dynamic motion in a patient. The required 4D dose matrix was generated by splitting the original treatment plan into multiple sub-beams, each representing 0.4 s dose delivery, and recalculating the dose of the split plan in the TPS (Eclipse). The dose accumulation was performed via TPS scripting by querying the dose of each sub-beam in dynamically moving points, allowing dose reconstruction with any dynamic motion. The dose reconstruction was validated with film dosimetry for two prostate dual arc VMAT plans with intra-prostatic lesion boosts. The plans were delivered to a pelvis phantom with internal dynamic rotational motion of a film stack (21 films with 2.5 mm separation). Each plan was delivered without motion and with three prostate motion traces. Motion-including dose reconstruction was performed for each motion experiment using the actual dynamic rotation as well as a constant rotation equal to the mean rotation during the experiment. For each experiment, the 3%/2 mm γ failure rate of the TPS dose reconstruction was calculated with the film measurement being the reference. For each motion experiment, the motion-induced 3%/2 mm γ failure rate was calculated using the static delivery as the reference and compared between film measurements and TPS dose reconstruction. DVH metrics for RT structures fully contained in the film volume were also compared between film and TPS.

RESULTS

The mean γ failure rate of the TPS dose reconstructions when compared to film doses was 0.8% (two static experiments) and 1.7% (six dynamic experiments). The mean (range) of the motion-induced γ failure rate in film measurements was 35.4% (21.3-59.2%). The TPS dose reconstruction agreed with these experimental γ failure rates with root-mean-square errors of 2.1% (dynamic rotation dose reconstruction) and 17.1% (dose reconstruction assuming constant rotation). By DVH metrics, the mean (range) difference between dose reconstructions with dynamic and constant rotation was 4.3% (-0.3-10.6%) (urethra D 2 % ), -0.6% (-5.6%-2.5%) (urethra D 99 % ), 1.1% (-7.1-7.7%) (GTV D 2 % ), -1.4% (-17.4-7.1%) (GTV D 95 % ), -1.2% (-17.1-5.7%) (GTV D 99 % ), and -0.1% (-3.2-7.6%) (GTV mean dose). Dose reconstructions with dynamic motion revealed large interplay effects (cold and hot spots).

CONCLUSIONS

A method to perform dose reconstructions for dynamic 6DoF motion in a TPS was developed and experimentally validated. It revealed large differences in dose distribution between dynamic and constant rotations not identifiable through dose reconstructions with constant rotation.

Comprehensive characterization of ExacTrac stereoscopic image guidance system using Monte Carlo and Spektr simulations

The purpose of this work is to develop accurate computational methods to comprehensively characterize and model the clinical ExacTrac imaging system, which is used as an image guidance system for stereotactic treatment applications. The Spektr toolkit was utilized to simulate the spectral and imaging characterization of the system. Since Spektr only simulates the primary beam (ignoring scatter), a full model of ExacTrac was also developed in Monte Carlo (MC) to characterize the imaging system. To ensure proper performance of both simulation models, Spektr and MC data were compared to the measured spectral and half value layers (HVLs) values. To validate the simulation results, x-ray spectra of the ExacTrac system were measured for various tube potentials using a CdTe spectrometer with multiple added narrow collimators. The raw spectra were calibrated using a ⁵⁷Co source and corrected for the escape peaks and detector efficiency. HVLs in mm of Al for various energies were measured using a calibrated RaySafe detector. Spektr and MC HVLs were calculated and compared to the measured values. The patient surface dose was calculated for different clinical imaging protocols from the measured air kerma and HVL values following the TG-61 methodology. The x-ray focal spot was measured by slanted edge technique using gafchromic films. ExacTrac imaging system beam profiles were simulated for various energies by MC simulation and the results were benchmarked by experimentally acquired beam profiles using gafchromic films. The effect of 6D IGRT treatment couch on beam hardening, dynamic range of the flat panel detector and scatter effect were determined using both Spektr simulation and experimental measurements. The measured and simulated spectra (of both MC and Spektr) for various kVps were compared and agreed within acceptable error. As another validation, the measured HVLs agreed with the Spektr and MC simulated HVLs on average within 1.0% for all kVps. The maximum and minimum patient surface doses were found to be 1.06 mGy for shoulder (high) and 0.051 mGy for cranial (low) imaging protocols, respectively. The MC simulated beam profiles were well matched with experimental results and replicated the penumbral slopes, the heel effect, and out-of-field regions. Dynamic range of detector (in terms of air kerma at detector surface) was found to be in the range of [6.1 × 10^-6, 5.3 × 10^-3] mGy. Accurate MC and Spektr models of the ExacTrac image guidance system were successfully developed and benchmarked via experimental validation. While patient surface dose for available imaging protocols were reported in this study, the established MC model may be used to obtain 3D imaging dose distribution for real patient geometries.

Quantification and correction of the scattered X-rays from a megavoltage photon beam to a linac-mounted kilovoltage imaging subsystem

Objective

To quantify and correct megavoltage (MV) scattered X-rays (MV-scatter) on an image acquired using a linac-mounted kilovoltage (kV) imaging subsystem.

Methods and materials

A linac-mounted flat-panel detector (FPD) was used to acquire an image containing MV-scatter by activating the FPD only during MV beam irradiation. 6-, 10-, and 15 MV with a flattening-filter (FF; 6X-FF, 10X-FF, 15X-FF), and 6- and 10 MV without an FF (6X-FFF, 10X-FFF) were used. The maps were acquired by changing one of the irradiation parameters while the others remained fixed. The mean pixel values of the MV-scatter were normalized to the 6X-FF reference condition (MV-scatter value). An MV-scatter database was constructed using these values. An MV-scatter correction experiment with one full arc image acquisition and two square field sizes (FSs) was conducted. Measurement- and estimation-based corrections were performed using the database. The image contrast was calculated at each angle.

Results

The MV-scatter increased with a larger FS and dose rate. The MV-scatter value factor varied substantially depending on the FPD position or collimator rotation. The median relative error ranges of the contrast for the image without, and with the measurement- and estimation-based correction were -10.9 to -2.9, and -1.5 to 4.8 and -7.4 to 2.6, respectively, for an FS of 10.0 × 10.0 cm².

Conclusions

The MV-scatter was strongly dependent on the FS, dose rate, and FPD position. The MV-scatter correction improved the image contrast.

Advances in knowledge

The MV-scatters on the TrueBeam linac kV imaging subsystem were quantified with various MV beam parameters, and strongly depended on the fieldsize, dose rate, and flat panel detector position. The MV-scatter correction using the constructed database improved the image quality.

Technical Note: Synthetic treatment beam imaging for motion monitoring during spine SBRT treatments - a phantom study

PURPOSE

One of the biggest challenges in applying megavoltage (MV) treatment beam imaging for monitoring spine motion in stereotactic body radiotherapy (SBRT) is the small beam apertures in the images due to strong beam modulations in IMRT planning. The purpose of this study is to investigate the feasibility of a markerless motion tracking method in spine SBRT delivery using a novel enhanced synthetic treatment beam (ESTB) imaging technique.

METHODS

Three clinical spine SBRT plans using 6XFFF beams and sliding window IMRT technique were transferred to a thorax phantom and delivered by a TrueBeam machine. Before delivery, the phantom was aligned to the plan isocenter using CBCT setup and verified with a second CBCT, and then, 2 mm shifts were introduced in both the craniocaudal (CC) and the left-right (LR) directions with the couch. During beam delivery, MV images were continuously taken with an electronic portal imaging device (EPID) and automatically grabbed by Varian iTools Capture software with a frame rate of 11.6 Hz. After preprocessing for scatter correction and beam intensity compensation, every 50 frames of MV images were combined to generate a series of ESTB images for each beam. The ESTB images were then registered to the projections of the verification CBCT at the matched beam angles to detect the 2 mm shifts.

RESULTS

Compared to snapshot MV images, the ESTB images had significantly enlarged fields of view (FOVs) and improved image quality. Based on two-dimensional (2D) rigid registration, the ESTB image to CBCT projection matching showed submillimeter accuracy in detecting motion. Specifically, the root mean square errors in detecting the LR/CC shifts were 0.35/0.28, 0.32/0.35, 0.63/0.44, 0.55/0.51, and 0.69/0.42 mm at gantry angles 180, 160, 140, 120, and 100, respectively.

CONCLUSION

Our results in the phantom study suggest that ESTB images from a sliding window IMRT plan can be used to detect spine motion, with submillimeter precision in the 2D plane perpendicular to the beam.

ACR-ASTRO Practice Parameter for Image-guided Radiation Therapy (IGRT)

AIM/OBJECTIVES/BACKGROUND

The American College of Radiology (ACR) and the American Society for Radiation Oncology (ASTRO) have jointly developed the following practice parameter for image-guided radiation therapy (IGRT). IGRT is radiation therapy that employs imaging to maximize accuracy and precision throughout the entire process of treatment delivery with the goal of optimizing accuracy and reliability of radiation therapy to the target, while minimizing dose to normal tissues.

METHODS

The ACR-ASTRO Practice Parameter for IGRT was revised according to the process described on the ACR website ("The Process for Developing ACR Practice Parameters and Technical Standards," www.acr.org/ClinicalResources/Practice-Parametersand-Technical-Standards) by the Committee on Practice Parameters of the ACR Commission on Radiation Oncology in collaboration with the ASTRO. Both societies then reviewed and approved the document.

RESULTS

This practice parameter is developed to serve as a tool in the appropriate application of IGRT in the care of patients with conditions where radiation therapy is indicated. It addresses clinical implementation of IGRT including personnel qualifications, quality assurance standards, indications, and suggested documentation.

CONCLUSIONS

This practice parameter is a tool to guide clinical use of IGRT and does not make recommendations on site-specific IGRT directives. It focuses on the best practices and principles to consider when using IGRT effectively, especially with the significant increase in imaging data that is now available with IGRT. The clinical benefit and medical necessity of the imaging modality and frequency of IGRT should be assessed for each patient.

Image-guided Radiotherapy to Manage Respiratory Motion: Lung and Liver

Organ motion as a result of respiratory and cardiac motion poses significant challenges for the accurate delivery of radiotherapy to both the thorax and the upper abdomen. Modern imaging techniques during radiotherapy simulation and delivery now permit better quantification of organ motion, which in turn reduces tumour and organ at risk position uncertainty. These imaging advances, coupled with respiratory correlated radiotherapy delivery techniques, have led to the development of a range of approaches to manage respiratory motion. This review summarises the key strategies of image-guided respiratory motion management with a focus on lung and liver radiotherapy.

Improvement of kilovoltage intrafraction monitoring accuracy through gantry angles selection

Kilovoltage intrafraction monitoring (KIM) is a method allowing to precisely infer the tumour trajectory based on cone beam computed tomography (CBCT) 2D-projections. However, its accuracy is deteriorated in the case of highly mobile tumours involving hysteresis. A first adaptation of KIM consisting of a prior amplitude based binning step has been developed in order to minimize the errors of the original model (phase-KIM). In this work, we propose enhanced methods (KIM_{sub-arc optim}and phase-KIM_{sub-arc optim}) to improve the accuracy of KIM and phase-KIM which relies on the selection of the optimal starting CBCT gantry angle. Aiming at demonstrating the interest of our approach, we carried out a simulation study and an experimental study: we compared the accuracy of the conventional versus sub-arc optim methods on simulated realistic tumour motions with amplitudes ranging from 5 to 30 mm in 1 mm increments. The same approach was performed using a lung dynamic phantom generating a 30 mm amplitude sinusoidal motion. The results show that for in-silico simulated motions of 10, 20 and 30 mm amplitude, the three-dimensional root mean square error (3D-RMSE) can be reduced by 0.67 mm, 0.91 mm, 0.94 mm and 0.18 mm, 0.25 mm, 0.28 mm using KIM_{sub-arc optim}and phase-KIM_{sub-arc optim}respectively. Considering all in-silico simulated trajectories, the percentage of errors larger than 1 mm decreases from 21.9% down to 1.6% for KIM (p < 0.001) and from 6.6% down to 1.2% for phase-KIM (p < 0.001). Experimentally, the 3D-RMSE is lowered by 0.5732 mm for KIM and by 0.1 mm for phase-KIM. The percentage of errors larger than 1 mm falls from 39.7% down to 18.5% for KIM and from 23.2% down to 11.1% for phase-KIM. In conclusion, our method efficiently anticipates CBCT gantry angles associated with a significantly better accuracy by using KIM and phase-KIM.

Boosting radiotherapy dose calculation accuracy with deep learning

In radiotherapy, a trade‐off exists between computational workload/speed and dose calculation accuracy. Calculation methods like pencil‐beam convolution can be much faster than Monte‐Carlo methods, but less accurate. The dose difference, mostly caused by inhomogeneities and electronic disequilibrium, is highly correlated with the dose distribution and the underlying anatomical tissue density. We hypothesize that a conversion scheme can be established to boost low‐accuracy doses to high‐accuracy, using intensity information obtained from computed tomography (CT) images. A deep learning‐driven framework was developed to test the hypothesis by converting between two commercially available dose calculation methods: Anisotropic analytic algorithm (AAA) and Acuros XB (AXB). A hierarchically dense U‐Net model was developed to boost the accuracy of AAA dose toward the AXB level. The network contained multiple layers of varying feature sizes to learn their dose differences, in relationship to CT, both locally and globally. Anisotropic analytic algorithm and AXB doses were calculated in pairs for 120 lung radiotherapy plans covering various treatment techniques, beam energies, tumor locations, and dose levels. For each case, the CT and the AAA dose were used as the input and the AXB dose as the “ground‐truth” output, to train and test the model. The mean squared errors (MSEs) and gamma passing rates (2 mm/2% & 1 mm/1%) were calculated between the boosted AAA doses and the “ground‐truth” AXB doses. The boosted AAA doses demonstrated substantially improved match to the “ground‐truth” AXB doses, with average (± s.d.) gamma passing rate (1 mm/1%) 97.6% (±2.4%) compared to 87.8% (±9.0%) of the original AAA doses. The corresponding average MSE was 0.11(±0.05) vs 0.31(±0.21). Deep learning is able to capture the differences between dose calculation algorithms to boost the low‐accuracy algorithms. By combining a less accurate dose calculation algorithm with a trained deep learning model, dose calculation can potentially achieve both high accuracy and efficiency.

Real-Time Image Guided Ablative Prostate Cancer Radiation Therapy: Results From the TROG 15.01 SPARK Trial

PURPOSE

Kilovoltage intrafraction monitoring (KIM) is a novel software platform implemented on standard radiation therapy systems and enabling real-time image guided radiation therapy (IGRT). In a multi-institutional prospective trial, we investigated whether real-time IGRT improved the accuracy of the dose patients with prostate cancer received during radiation therapy.

METHODS AND MATERIALS

Forty-eight patients with prostate cancer were treated with KIM-guided SABR with 36.25 Gy in 5 fractions. During KIM-guided treatment, the prostate motion was corrected for by either beam gating with couch shifts or multileaf collimator tracking. A dose reconstruction method was used to evaluate the dose delivered to the target and organs at risk with and without real-time IGRT. Primary outcome was the effect of real-time IGRT on dose distributions. Secondary outcomes included patient-reported outcomes and toxicity.

RESULTS

Motion correction occurred in ≥1 treatment for 88% of patients (42 of 48) and 51% of treatments (121 of 235). With real-time IGRT, no treatments had prostate clinical target volume (CTV) D98% dose 5% less than planned. Without real-time IGRT, 13 treatments (5.5%) had prostate CTV D98% doses 5% less than planned. The prostate CTV D98% dose with real-time IGRT was closer to the plan by an average of 1.0% (range, -2.8% to 20.3%). Patient outcomes showed no change in the 12-month patient-reported outcomes compared with baseline and no grade ≥3 genitourinary or gastrointestinal toxicities.

CONCLUSIONS

Real-time IGRT is clinically effective for prostate cancer SABR.

Technical Note: 3D localization of lung tumors on cone beam CT projections via a convolutional recurrent neural network

PURPOSE

To design a convolutional recurrent neural network (CRNN) that calculates three-dimensional (3D) positions of lung tumors from continuously acquired cone beam computed tomography (CBCT) projections, and facilitates the sorting and reconstruction of 4D-CBCT images.

METHOD

Under an IRB-approved clinical lung protocol, kilovoltage (kV) projections of the setup CBCT were collected in free-breathing. Concurrently, an electromagnetic signal-guided system recorded motion traces of three transponders implanted in or near the tumor. Convolutional recurrent neural network was designed to utilize a convolutional neural network (CNN) for extracting relevant features of the kV projections around the tumor, followed by a recurrent neural network for analyzing the temporal patterns of the moving features. Convolutional recurrent neural network was trained on the simultaneously collected kV projections and motion traces, subsequently utilized to calculate motion traces solely based on the continuous feed of kV projections. To enhance performance, CRNN was also facilitated by frequent calibrations (e.g., at 10° gantry rotation intervals) derived from cross-correlation-based registrations between kV projections and templates created from the planning 4DCT. Convolutional recurrent neural network was validated on a leave-one-out strategy using data from 11 lung patients, including 5500 kV images. The root-mean-square error between the CRNN and motion traces was calculated to evaluate the localization accuracy.

RESULT

Three-dimensional displacement around the simulation position shown in the Calypso traces was 3.4 ± 1.7 mm. Using motion traces as ground truth, the 3D localization error of CRNN with calibrations was 1.3 ± 1.4 mm. CRNN had a success rate of 86 ± 8% in determining whether the motion was within a 3D displacement window of 2 mm. The latency was 20 ms when CRNN ran on a high-performance computer cluster.

CONCLUSIONS

CRNN is able to provide accurate localization of lung tumors with aid from frequent recalibrations using the conventional cross-correlation-based registration approach, and has the potential to remove reliance on the implanted fiducials.

A Review on Curability of Cancers: More Efforts for Novel Therapeutic Options Are Needed

Cancer remains a major cause of death globally. Given its relapsing and fatal features, curing cancer seems to be something hardly possible for the majority of patients. In view of the development in cancer therapies, this article summarizes currently available cancer therapeutics and cure potential by cancer type and stage at diagnosis, based on literature and database reviews. Currently common cancer therapeutics include surgery, chemotherapy, radiotherapy, targeted therapy, and immunotherapy. However, treatment with curative intent by these methods are mainly eligible for patients with localized disease or treatment-sensitive cancers and therefore their contributions to cancer curability are relatively limited. The prognosis for cancer patients varies among different cancer types with a five-year relative survival rate (RSR) of more than 80% in thyroid cancer, melanoma, breast cancer, and Hodgkin's lymphoma. The most dismal prognosis is observed in patients with small-cell lung cancer, pancreatic cancer, hepatocellular carcinoma, oesophagal cancer, acute myeloid leukemia, non-small cell lung cancer, and gastric cancer with a five-year RSR ranging between 7% and 28%. The current review is intended to provide a general view about how much we have achieved in curing cancer as regards to different therapies and cancer types. Finally, we propose a small molecule dual-targeting broad-spectrum anticancer strategy called OncoCiDia, in combination with emerging highly sensitive liquid biopsy, with theoretical curative potential for the management of solid malignancies, especially at the micro-cancer stage.

Stereotactic body radiation therapy with optional focal lesion ablative microboost in prostate cancer: Topical review and multicenter consensus

Stereotactic body radiotherapy (SBRT) for prostate cancer (PCa) is gaining interest by the recent publication of the first phase III trials on prostate SBRT and the promising results of many other phase II trials. Before long term results became available, the major concern for implementing SBRT in PCa in daily clinical practice was the potential risk of late genitourinary (GU) and gastrointestinal (GI) toxicity. A number of recently published trials, including late outcome and toxicity data, contributed to the growing evidence for implementation of SBRT for PCa in daily clinical practice. However, there exists substantial variability in delivering SBRT for PCa. The aim of this topical review is to present a number of prospective trials and retrospective analyses of SBRT in the treatment of PCa. We focus on the treatment strategies and techniques used in these trials. In addition, recent literature on a simultaneous integrated boost to the tumor lesion, which could create an additional value in the SBRT treatment of PCa, was described. Furthermore, we discuss the multicenter consensus of the FLAME consortium on SBRT for PCa with a focal boost to the macroscopic intraprostatic tumor nodule(s).

Real-time intrafraction motion monitoring in external beam radiotherapy

Radiotherapy (RT) aims to deliver a spatially conformal dose of radiation to tumours while maximizing the dose sparing to healthy tissues. However, the internal patient anatomy is constantly moving due to respiratory, cardiac, gastrointestinal and urinary activity. The long term goal of the RT community to 'see what we treat, as we treat' and to act on this information instantaneously has resulted in rapid technological innovation. Specialized treatment machines, such as robotic or gimbal-steered linear accelerators (linac) with in-room imaging suites, have been developed specifically for real-time treatment adaptation. Additional equipment, such as stereoscopic kilovoltage (kV) imaging, ultrasound transducers and electromagnetic transponders, has been developed for intrafraction motion monitoring on conventional linacs. Magnetic resonance imaging (MRI) has been integrated with cobalt treatment units and more recently with linacs. In addition to hardware innovation, software development has played a substantial role in the development of motion monitoring methods based on respiratory motion surrogates and planar kV or Megavoltage (MV) imaging that is available on standard equipped linacs. In this paper, we review and compare the different intrafraction motion monitoring methods proposed in the literature and demonstrated in real-time on clinical data as well as their possible future developments. We then discuss general considerations on validation and quality assurance for clinical implementation. Besides photon RT, particle therapy is increasingly used to treat moving targets. However, transferring motion monitoring technologies from linacs to particle beam lines presents substantial challenges. Lessons learned from the implementation of real-time intrafraction monitoring for photon RT will be used as a basis to discuss the implementation of these methods for particle RT.

Comparison of treatment plan quality among MRI-based IMRT with a linac, MRI-based IMRT with tri-Co-60 sources, and VMAT for spine SABR

PURPOSE

This study compares the plan quality of magnetic-resonance image (MRI)-based intensity modulated radiation therapy (IMRT) using a linac (MR-linac-IMRT), MRI-based IMRT using tri-Co-60 sources (MR-Co-60-IMRT), and volumetric modulated arc therapy (VMAT) for spine stereotactic ablative radiotherapy (SABR).

METHODS

Twenty patients with thoracic spine metastasis were retrospectively selected for this study. For each patient, the MR-linac-IMRT, MR-Co-60-IMRT, and VMAT plans were generated using an identical CT image set and structures, except for the spinal cord and spinal cord planning organ-at-risk volume (PRV). Those two structures were contoured based on CT image sets for VMAT planning while those were contoured based on MR image sets for MR-linac-IMRT and MR-Co-60-IMRT planning. The initial prescription doses were 18 Gy in a single fraction for every plan in this study. If the tolerance level of the spinal cord was not met, the prescription doses were reduced to meet the tolerance level of the spinal cord. Dose-volumetric parameters of each plan were analyzed.

RESULTS

The average spinal cord volumes contoured based on the CT and MR images were 3.8±1.6 cm3 and 1.1±1.0 cm3, respectively (p<0.001). For four patients, the prescription doses of VMAT plans were reduced to 16 Gy to satisfy the spinal cord tolerance level. For thirteen patients, the prescription doses of MR-Co-60-IMRT plans were reduced to be less than 16 Gy to meet the spinal cord tolerance level. However, for every MR-linac-IMRT plan, the initial prescription doses of 18 Gy could be delivered to the target volume while satisfying the spinal cord tolerance. The average values of D10%, V10Gy, and V14Gy of the spinal cord PRV consistently indicated that the doses to the spinal cord PRV in the MR-linac-IMRT plans were the lowest among three types of plans in this study (all with p≤0.003).

CONCLUSION

MR-linac-IMRT appears promising for spine SABR.

Incorporating imaging information from deep neural network layers into image guided radiation therapy (IGRT)

BACKGROUND AND PURPOSE

To investigate a novel markerless prostate localization strategy using a pre-trained deep learning model to interpret routine projection kilovoltage (kV) X-ray images in image-guided radiation therapy (IGRT).

MATERIALS AND METHODS

We developed a personalized region-based convolutional neural network to localize the prostate treatment target without implanted fiducials. To train the deep neural network (DNN), we used the patient's planning computed tomography (pCT) images with pre-delineated prostate target to generate a large amount of synthetic kV projection X-ray images in the geometry of onboard imager (OBI) system. The DNN model was evaluated by retrospectively studying 10 patients who underwent prostate IGRT. Three out of the ten patients who had implanted fiducials and the fiducials' positions in the OBI images acquired for treatment setup were examined to show the potential of the proposed method for prostate IGRT. Statistical analysis using Lin's concordance correlation coefficient was calculated to assess the results along with the difference between the digitally reconstructed radiographs (DRR) derived and DNN predicted locations of the prostate.

RESULTS

Differences between the predicted target positions using DNN and their actual positions are (mean ± standard deviation) 1.58 ± 0.43 mm, 1.64 ± 0.43 mm, and 1.67 ± 0.36 mm in anterior-posterior, lateral, and oblique directions, respectively. Prostate position identified on the OBI kV images is also found to be consistent with that derived from the implanted fiducials.

CONCLUSIONS

Highly accurate, markerless prostate localization based on deep learning is achievable. The proposed method is useful for daily patient positioning and real-time target tracking during prostate radiotherapy.

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.

Developing a MLC modifier program to improve fiducial detection for MV/kV imaging during hypofractionated prostate volumetric modulated arc therapy

Purpose

To develop an Eclipse plug‐in (MLC_MODIFIER) that automatically modifies control points to expose fiducials obscured by MLC during VMAT, thereby facilitating tracking using periodic MV/kV imaging.

Method

Three‐dimensional fiducial tracking was performed during VMAT by pairing short‐arc (3°) MV digital tomosynthesis (DTS) images to triggered kV images. To evaluate MLC_MODIFIER efficacy, two cohorts of patients were considered. For first 12 patients, plans were manually edited to expose one fiducial marker. Next for 15 patients, plans were modified using MLC_MODIFIER script. MLC_MODIFIER evaluated MLC apertures at appropriate angles for marker visibility. Angles subtended by control points were compressed and low‐dose “imaging” control points were inserted and exposed one marker with 1 cm margin. Patient's images were retrospectively reviewed to determine rate of MV registration failures. Failure categories were poor DTS image quality, MLC blockage of fiducials, or unknown reasons. Dosimetric differences in rectum, bladder, and urethra D1 cc, PTV maximum dose, and PTV dose homogeneity (PTV HI) were evaluated. Statistical significance was evaluated using Fisher's exact and Student's t test.

Result

Overall MV registration failures, failures due to poor image quality, MLC blockage, and unknown reasons were 33% versus 8.9% (P < 0.0001), 8% versus 6.4% (P < 0.05), 13.6% versus 0.1% (P < 0.0001), and 7.6% versus 2.4% (P < 0.0001) for manually edited and MLC_MODIFIER plans, respectively. PTV maximum and HI increased on average from unmodified plans by 2.1% and 0.3% (P < 0.004) and 22.0% and 3.3% (P < 0.004) for manually edited and MLC_MODIFIED plans, respectively. Changes in bladder, rectum, and urethra D1CC were similar for each method and less than 0.7%.

Conclusion

Increasing fiducial visibility via an automated process comprised of angular compression of control points and insertion of additional “imaging” control points is feasible. Degradation of plan quality is minimal. Fiducial detection and registration success rates are significantly improved compared to manually edited apertures.

See, Think, and Act: Real-Time Adaptive Radiotherapy

The world is embracing the information age, with real-time data at hand to assist with many decisions. Similarly, in cancer radiotherapy we are inexorably moving toward using information in a smarter and faster fashion, to usher in the age of real-time adaptive radiotherapy. The three critical steps of real-time adaptive radiotherapy, aligned with driverless vehicle technology are a continuous see, think, and act loop. See: use imaging systems to probe the patient anatomy or physiology as it evolves with time. Think: use current and prior information to optimize the treatment using the available adaptive degrees of freedom. Act: deliver the real-time adapted treatment. This paper expands upon these three critical steps for real-time adaptive radiotherapy, provides a historical context, reviews the clinical rationale, and gives a future outlook for real-time adaptive radiotherapy.