Influence of deformable image registration on 4D dose simulation for extracranial SBRT: A multi-registration framework study

Stability analysis of patient-specific 4DCT- and 4DCBCT-based correspondence models

BACKGROUND

Surrogate-based motion compensation in stereotactic body radiation therapy (SBRT) strongly relies on a constant relationship between an external breathing signal and the internal tumor motion over the course of treatment, that is, a stable patient-specific correspondence model.

PURPOSE

This study aims to develop methods for analyzing the stability of correspondence models by integrating planning 4DCT and pretreatment 4D cone-beam computed tomography (4DCBCT) data and assessing the relation to patient-specific clinical parameters.

METHODS

For correspondence modeling, a regression-based approach is applied, correlating patient-specific internal motion (vector fields computed by deformable image registration) and external breathing signals (recorded by Varian's RPM and RGSC system). To analyze correspondence model stability, two complementary methods are proposed. (1) Target volume-based analysis: 4DCBCT-based correspondence models predict clinical target volumes (GTV and internal target volume [ITV]) within the planning 4DCT, which are evaluated by overlap and distance measures (Dice similarity coefficient [DSC]/average symmetric surface distance [ASSD]). (2) System matrix-based analysis: 4DCBCT-based regression models are compared to 4DCT-based models using mean squared difference (MSD) and principal component analysis of the system matrices. Stability analysis results are correlated with clinical parameters. Both methods are applied to a dataset of 214 pretreatment 4DCBCT scans (Varian TrueBeam) from a cohort of 46 lung tumor patients treated with ITV-based SBRT (planning 4DCTs acquired with Siemens AS Open and SOMATOM go.OPEN Pro CT scanners).

RESULTS

Consistent results across the two complementary analysis approaches (Spearman correlation coefficient of0.6 / 0.7 $0.6/ 0.7$ between system matrix-based MSD and GTV-based DSC/ASSD) were observed. Analysis showed that stability was not predominant, with 114/214 fraction-wise models not surpassing a threshold ofD S C > 0.7 $DSC > 0.7$ for the GTV, and only 14/46 patients demonstrating aD S C > 0.7 $DSC > 0.7$ in all fractions. Model stability did not degrade over the course of treatment. The mean GTV-based DSC is0.59 ± 0.26 $0.59\pm 0.26$ (mean ASSD of2.83 ± 3.37 $2.83\pm 3.37$ ) and the respective ITV-based DSC is0.69 ± 0.20 $0.69\pm 0.20$ (mean ASSD of2.35 ± 1.81 $2.35\pm 1.81$ ). The clinical parameters showed a strong correlation between smaller tumor motion ranges and increased stability.

CONCLUSIONS

The proposed methods identify patients with unstable correspondence models prior to each treatment fraction, serving as direct indicators for the necessity of replanning and adaptive treatment approaches to account for internal-external motion variations throughout the course of treatment.

Cardiorespiratory motion characteristics and their dosimetric impact on cardiac stereotactic body radiotherapy

BACKGROUND

Cardiac stereotactic body radiotherapy (CSBRT) is an emerging and promising noninvasive technique for treating refractory arrhythmias utilizing highly precise, single or limited-fraction high-dose irradiations. This method promises to revolutionize the treatment of cardiac conditions by delivering targeted therapy with minimal exposure to surrounding healthy tissues. However, the dynamic nature of cardiorespiratory motion poses significant challenges to the precise delivery of dose in CSBRT, introducing potential variabilities that can impact treatment efficacy. The complexities of the influence of cardiorespiratory motion on dose distribution are compounded by interplay and blurring effects, introducing additional layers of dose uncertainty. These effects, critical to the understanding and improvement of the accuracy of CSBRT, remain unexplored, presenting a gap in current clinical literature.

PURPOSE

To investigate the cardiorespiratory motion characteristics in arrhythmia patients and the dosimetric impact of interplay and blurring effects induced by cardiorespiratory motion on CSBRT plan quality.

METHODS

The position and volume variations in the substrate target and cardiac substructures were evaluated in 12 arrhythmia patients using displacement maximum (DMX) and volume metrics. Moreover, a four-dimensional (4D) dose reconstruction approach was employed to examine the dose uncertainty of the cardiorespiratory motion.

RESULTS

Cardiac pulsation induced lower DMX than respiratory motion but increased the coefficient of variation and relative range in cardiac substructure volumes. The mean DMX of the substrate target was 0.52 cm (range: 0.26-0.80 cm) for cardiac pulsation and 0.82 cm (range: 0.32-2.05 cm) for respiratory motion. The mean DMX of the cardiac structure ranged from 0.15 to 1.56 cm during cardiac pulsation and from 0.35 to 1.89 cm during respiratory motion. Cardiac pulsation resulted in an average deviation of -0.73% (range: -4.01%-4.47%) in V25 between the 3D and 4D doses. The mean deviations in the homogeneity index (HI) and gradient index (GI) were 1.70% (range: -3.10%-4.36%) and 0.03 (range: -0.14-0.11), respectively. For cardiac substructures, the deviations in D50 due to cardiac pulsation ranged from -1.88% to 1.44%, whereas the deviations in Dmax ranged from -2.96% to 0.88% of the prescription dose. By contrast, the respiratory motion led to a mean deviation of -1.50% (range: -10.73%-4.23%) in V25. The mean deviations in HI and GI due to respiratory motion were 4.43% (range: -3.89%-13.98%) and 0.18 (range: -0.01-0.47) (p < 0.05), respectively. Furthermore, the deviations in D50 and Dmax in cardiac substructures for the respiratory motion ranged from -0.28% to 4.24% and -4.12% to 1.16%, respectively.

CONCLUSIONS

Cardiorespiratory motion characteristics vary among patients, with the respiratory motion being more significant. The intricate cardiorespiratory motion characteristics and CSBRT plan complexity can induce substantial dose uncertainty. Therefore, assessing individual motion characteristics and 4D dose reconstruction techniques is critical for implementing CSBRT without compromising efficacy and safety.

Deep learning-based conditional inpainting for restoration of artifact-affected 4D CT images

BACKGROUND

4D CT imaging is an essential component of radiotherapy of thoracic and abdominal tumors. 4D CT images are, however, often affected by artifacts that compromise treatment planning quality and image information reliability.

PURPOSE

In this work, deep learning (DL)-based conditional inpainting is proposed to restore anatomically correct image information of artifact-affected areas.

METHODS

The restoration approach consists of a two-stage process: DL-based detection of common interpolation (INT) and double structure (DS) artifacts, followed by conditional inpainting applied to the artifact areas. In this context, conditional refers to a guidance of the inpainting process by patient-specific image data to ensure anatomically reliable results. The study is based on 65 in-house 4D CT images of lung cancer patients (48 with only slight artifacts, 17 with pronounced artifacts) and two publicly available 4D CT data sets that serve as independent external test sets.

RESULTS

Automated artifact detection revealed a ROC-AUC of 0.99 for INT and of 0.97 for DS artifacts (in-house data). The proposed inpainting method decreased the average root mean squared error (RMSE) by 52 % (INT) and 59 % (DS) for the in-house data. For the external test data sets, the RMSE improvement is similar (50 % and 59 %, respectively). Applied to 4D CT data with pronounced artifacts (not part of the training set), 72 % of the detectable artifacts were removed.

CONCLUSIONS

The results highlight the potential of DL-based inpainting for restoration of artifact-affected 4D CT data. Compared to recent 4D CT inpainting and restoration approaches, the proposed methodology illustrates the advantages of exploiting patient-specific prior image information.

Review and recommendations on deformable image registration uncertainties for radiotherapy applications

Deformable image registration (DIR) is a versatile tool used in many applications in radiotherapy (RT). DIR algorithms have been implemented in many commercial treatment planning systems providing accessible and easy-to-use solutions. However, the geometric uncertainty of DIR can be large and difficult to quantify, resulting in barriers to clinical practice. Currently, there is no agreement in the RT community on how to quantify these uncertainties and determine thresholds that distinguish a good DIR result from a poor one. This review summarises the current literature on sources of DIR uncertainties and their impact on RT applications. Recommendations are provided on how to handle these uncertainties for patient-specific use, commissioning, and research. Recommendations are also provided for developers and vendors to help users to understand DIR uncertainties and make the application of DIR in RT safer and more reliable.

Accuracy and consistency of intensity-based deformable image registration in 4DCT for tumor motion estimation in liver radiotherapy planning

We investigate the accuracy of intensity-based deformable image registration (DIR) for tumor localization in liver stereotactic body radiotherapy (SBRT). We included 4DCT scans to capture the breathing motion of eight patients receiving SBRT for liver metastases within a retrospective clinical study. Each patient had three fiducial markers implanted. The liver and the tumor were delineated in the mid-ventilation phase, and their positions in the other phases were estimated with deformable image registration. We tested referenced and sequential registrations strategies. The fiducial markers were the gold standard to evaluate registration accuracy. The registration errors related to measured versus estimated fiducial markers showed a mean value less than 1.6mm. The positions of some fiducial markers appeared not stable on the 4DCT throughout the respiratory phases. Markers’ center of mass tends to be a more reliable measurement. Distance errors of tumor location based on registration versus markers center of mass were less than 2mm. There were no statistically significant differences between the reference and the sequential registration, i.e., consistency and errors were comparable to resolution errors. We demonstrated that intensity-based DIR is accurate up to resolution level for locating the tumor in the liver during breathing motion.

Improving deformable image registration with point metric and masking technique for postoperative breast cancer radiotherapy

Background

Deformable image registration (DIR) is increasingly used for target volume definition in radiotherapy. However, this method is challenging for postoperative breast cancer patients due to the large deformations and non-correspondence caused by tumor resection and clip insertion. In this study, an improved B-splines based DIR method was developed to address this issue for higher registration accuracy.

Methods

The conventional B-splines based DIR method was improved with the introduction of point metric and masking technique. The point metric minimizes the distance between 2 point sets with known correspondence for regularization of intensity-based B-splines registration. The masking technique reduces the influence of non-corresponding regions in breast computed tomography (CT) images. Two sets of CT images before and after breast surgery were used for image registration. One set was the diagnostic CT image acquired before surgery, and another set was the planning CT image acquired after surgery for breast cancer radiotherapy. A total of 26 sets of CT images from 13 patients were collected retrospectively for the test. The improved DIR method's registration accuracy was evaluated by target registration error (TRE), the Jacobian determinant, and visual assessment.

Results

For soft tissue, the difference in the median TRE between the improved DIR method and the conventional DIR method was statistically significant (2.27 vs. 5.88, P<0.05). The Jacobian determinant of the deformation field was positive for all patients. For visual assessment, the improved DIR method with point metric achieved better matching for soft tissue.

Conclusions

The improved DIR method's registration accuracy was higher than the conventional DIR method based on the preliminary results. With point metric and masking technique, the influence of large deformations and non-correspondence on registration between pre- and post-operative CT images can be effectively reduced. Therefore, this method provides a feasible way for target volume definition in postoperative breast cancer radiotherapy treatment planning.

Technological quality requirements for stereotactic radiotherapy : Expert review group consensus from the DGMP Working Group for Physics and Technology in Stereotactic Radiotherapy

This review details and discusses the technological quality requirements to ensure the desired quality for stereotactic radiotherapy using photon external beam radiotherapy as defined by the DEGRO Working Group Radiosurgery and Stereotactic Radiotherapy and the DGMP Working Group for Physics and Technology in Stereotactic Radiotherapy. The covered aspects of this review are 1) imaging for target volume definition, 2) patient positioning and target volume localization, 3) motion management, 4) collimation of the irradiation and beam directions, 5) dose calculation, 6) treatment unit accuracy, and 7) dedicated quality assurance measures. For each part, an expert review for current state-of-the-art techniques and their particular technological quality requirement to reach the necessary accuracy for stereotactic radiotherapy divided into intracranial stereotactic radiosurgery in one single fraction (SRS), intracranial fractionated stereotactic radiotherapy (FSRT), and extracranial stereotactic body radiotherapy (SBRT) is presented. All recommendations and suggestions for all mentioned aspects of stereotactic radiotherapy are formulated and related uncertainties and potential sources of error discussed. Additionally, further research and development needs in terms of insufficient data and unsolved problems for stereotactic radiotherapy are identified, which will serve as a basis for the future assignments of the DGMP Working Group for Physics and Technology in Stereotactic Radiotherapy. The review was group peer-reviewed, and consensus was obtained through multiple working group meetings.

Investigation of 4D dose in volumetric modulated arc therapy-based stereotactic body radiation therapy: does fractional dose or number of arcs matter?

The aim of this study was to assess the impact of fractional dose and the number of arcs on interplay effects when volumetric modulated arc therapy (VMAT) is used to treat lung tumors with large respiratory motions. A three (fractional dose of 4, 7.5 or 12.5 Gy) by two (number of arcs, one or two) VMAT plan was created for 10 lung cancer cases. The median 3D tumor motion was 17.9 mm (range: 8.2-27.2 mm). Ten phase-specific subplans were generated by calculating the dose on each respiratory phase computed tomography (CT) scan using temporally assigned VMAT arcs. We performed temporal assignment of VMAT arcs using respiratory information obtained from infrared markers placed on the abdomens of the patients during CT simulations. Each phase-specific dose distribution was deformed onto exhale phase CT scans using contour-based deformable image registration, and a 4D plan was created by dose accumulation. The gross tumor volume dose of each 4D plan (4D GTV dose) was compared with the internal target volume dose of the original plan (3D ITV dose). The near-minimum 4D GTV dose (D99%) was higher than the near-minimum 3D internal target volume (ITV) dose, whereas the near-maximum 4D GTV dose (D1%) was lower than the near-maximum 3D ITV dose. However, the difference was negligible, and thus the 4D GTV dose corresponded well with the 3D ITV dose, regardless of the fractional dose and number of arcs. Therefore, interplay effects were negligible in VMAT-based stereotactic body radiation therapy for lung tumors with large respiratory motions.

Evaluation of Intensity- and Contour-Based Deformable Image Registration Accuracy in Pancreatic Cancer Patients

We aimed to clarify the accuracy of rigid image registration and deformable image registration (DIR) in carbon-ion radiotherapy (CIRT) for pancreatic cancer. Six patients with pancreatic cancer who were treated with passive irradiation CIRT were enrolled. Three registration patterns were evaluated: treatment planning computed tomography images (TPCT) to CT images acquired in the treatment room (IRCT) in the supine position, TPCT to IRCT in the prone position, and TPCT in the supine position to the prone position. After warping the contours of the original CT images to the destination CT images using deformation matrices from the registration, the warped delineated contours on the destination CT images were compared with the original ones using mean displacement to agreement (MDA). Four contours (clinical target volume (CTV), gross tumor volume (GTV), stomach, duodenum) and four registration algorithms (rigid image registration [RIR], intensity-based DIR [iDIR], contour-based DIR [cDIR], and a hybrid iDIR-cDIR ([hDIR]) were evaluated. The means ± standard deviation of the MDAs of all contours for RIR, iDIR, cDIR, and hDIR were 3.40 ± 3.30, 2.2 1± 2.48, 1.46 ± 1.49, and 1.46 ± 1.37 mm, respectively. There were significant differences between RIR and iDIR, and between RIR/iDIR and cDIR/hDIR. For the pancreatic cancer patient images, cDIR and hDIR had better accuracy than RIR and iDIR.

Validation of a model for physical dose variations in irregularly moving targets treated with carbon ion beams

PURPOSE

In particle therapy, conventional treatment planning systems rely on an imaging representation of the irradiated region to compute the dose. For irregular breathing, when an imaging dataset describing the actual motion is not available, a different approach for dose estimation is needed. To this aim, we validate a method for the estimation of physical dose variations in gated carbon ion treatments, providing also a demonstration of the feasibility of physical dose metrics to assess the method performance. Finally, we describe a sample use case, in which this method is used to assess plan robustness with respect to undetected irregular tumor motion.

METHODS

The method entails the definition of a patient- and beam-specific water equivalent depth (WED) space, the simulation of motion as a translation equal to tumor displacement, and the reconstruction of the altered dose. We validated the approach using four-dimensional computed tomographies (4DCTs) and clinical plans in 12 patients, treated with respiratory gated carbon ion beams at the National Centre for Oncological Hadrontherapy (Pavia, Italy). Using the end-exhale CT and dose distribution as a reference, the physical dose delivered at the end-inhale tumor position was estimated and compared to the ground-truth dose recalculation on the end-inhale CT. Biologically effective and physical dose variations between the plan and the recalculation were compared as well. As a use case, we evaluated dose changes caused by simulated irregular tumor motion, that is, linear and nonlinear baseline shifts and/or amplitude variations with hysteresis.

RESULTS

The ratio between biologically effective and physical equivalent uniform dose (EUD) variations due to end-exhale to end-inhale motion was less than one for 96% of investigated structures. In the validation study, we found a median error corresponding to a 14% EUD overestimation for the tumor and 4% EUD underestimation for a subgroup of organs at risk, together with a high EUD variation due to motion [median 352% EUD variation between end-exhale and end-inhale doses in the planning tumor volume (PTV)]. Considering relevant dose-volume histogram (DVH) metrics, the median difference between estimated and ground truth doses was ≤ 4%. Gamma analysis between estimated and recalculated dose distributions resulted in a pass rate > 80% for 83% of the target volumes. For the two patients selected for the sample use case, a patient-specific assessment of the method performance was performed on the 4DCT and it was possible to relate EUD variations of both tumor and organs at risk to the simulated target motion.

CONCLUSIONS

The physical dose distribution was found to be more sensitive to motion with respect to the biologically effective one, suggesting the suitability of the physical dose metrics for the WED-space method validation. We showed that the method can compensate for intra-fractional tumor motion with proper accuracy in the selected patient group, although its use is recommended when limited deformations are expected. In conclusion, the WED-space method can provide simulations of dose alteration due to irregular breathing when imaging data are lacking, and, once integrated with relative biological effectiveness (RBE) modeling, it would be useful in evaluating the robustness of carbon ion treatment plans.

Deformable image registration for radiation therapy: principle, methods, applications and evaluation

Background: Deformable image registration (DIR) is increasingly used in the field of radiation therapy (RT) to account for anatomical deformations. The aims of this paper are to describe the main applications of DIR in RT and discuss current DIR evaluation methods. Methods: Articles on DIR published from January 2000 to October 2018 were extracted from PubMed and Science Direct. Our search was restricted to articles that report data obtained from humans, were written in English, and address DIR methods for RT. A total of 207 articles were selected from among 2506 identified in the search process. Results: At planning, DIR is used for organ delineation using atlas-based segmentation, deformation-based planning target volume definition, functional planning and magnetic resonance imaging-based dose calculation. In image-guided RT, DIR is used for contour propagation and dose calculation on per-treatment imaging. DIR is also used to determine the accumulated dose from fraction to fraction in external beam RT and brachytherapy, both for dose reporting and adaptive RT. In the case of re-irradiation, DIR can be used to estimate the cumulated dose of the two irradiations. Finally, DIR can be used to predict toxicity in voxel-wise population analysis. However, the evaluation of DIR remains an open issue, especially when dealing with complex cases such as the disappearance of matter. To quantify DIR uncertainties, most evaluation methods are limited to geometry-based metrics. Software companies have now integrated DIR tools into treatment planning systems for clinical use, such as contour propagation and fraction dose accumulation. Conclusions: DIR is increasingly important in RT applications, from planning to toxicity prediction. DIR is routinely used to reduce the workload of contour propagation. However, its use for complex dosimetric applications must be carefully evaluated by combining quantitative and qualitative analyses.