Surface-constrained nonrigid registration for dose monitoring in prostate cancer radiotherapy

Contour-guided deep learning based deformable image registration for dose monitoring during CBCT-guided radiotherapy of prostate cancer

PURPOSE

To evaluate deep learning (DL)-based deformable image registration (DIR) for dose accumulation during radiotherapy of prostate cancer patients.

METHODS AND MATERIALS

Data including 341 CBCTs (209 daily, 132 weekly) and 23 planning CTs from 23 patients was retrospectively analyzed. Anatomical deformation during treatment was estimated using free-form deformation (FFD) method from Elastix and DL-based VoxelMorph approaches. The VoxelMorph method was investigated using anatomical scans (VMorph_Sc) or label images (VMorph_Msk), or the combination of both (VMorph_Sc_Msk). Accumulated doses were compared with the planning dose.

RESULTS

The DSC ranges, averaged for prostate, rectum and bladder, were 0.60-0.71, 0.67-0.79, 0.93-0.98, and 0.89-0.96 for the FFD, VMorph_Sc, VMorph_Msk, and VMorph_Sc_Msk methods, respectively. When including both anatomical and label images, VoxelMorph estimated more complex deformations resulting in heterogeneous determinant of Jacobian and higher percentage of deformation vector field (DVF) folding (up to a mean value of 1.90% in the prostate). Large differences were observed between DL-based methods regarding estimation of the accumulated dose, showing systematic overdosage and underdosage of the bladder and rectum, respectively. The difference between planned mean dose and accumulated mean dose with VMorph_Sc_Msk reached a median value of +6.3 Gy for the bladder and -5.1 Gy for the rectum.

CONCLUSION

The estimation of the deformations using DL-based approach is feasible for male pelvic anatomy but requires the inclusion of anatomical contours to improve organ correspondence. High variability in the estimation of the accumulated dose depending on the deformable strategy suggests further investigation of DL-based techniques before clinical deployment.

Evaluation of an Open Source Registration Package for Automatic Contour Propagation in Online Adaptive Intensity-Modulated Proton Therapy of Prostate Cancer

Objective: Our goal was to investigate the performance of an open source deformable image registration package, elastix, for fast and robust contour propagation in the context of online-adaptive intensity-modulated proton therapy (IMPT) for prostate cancer.

Methods: A planning and 7–10 repeat CT scans were available of 18 prostate cancer patients. Automatic contour propagation of repeat CT scans was performed using elastix and compared with manual delineations in terms of geometric accuracy and runtime. Dosimetric accuracy was quantified by generating IMPT plans using the propagated contours expanded with a 2 mm (prostate) and 3.5 mm margin (seminal vesicles and lymph nodes) and calculating dosimetric coverage based on the manual delineation. A coverage of V_95% ≥ 98% (at least 98% of the target volumes receive at least 95% of the prescribed dose) was considered clinically acceptable.

Results: Contour propagation runtime varied between 3 and 30 s for different registration settings. For the fastest setting, 83 in 93 (89.2%), 73 in 93 (78.5%), and 91 in 93 (97.9%) registrations yielded clinically acceptable dosimetric coverage of the prostate, seminal vesicles, and lymph nodes, respectively. For the prostate, seminal vesicles, and lymph nodes the Dice Similarity Coefficient (DSC) was 0.87 ± 0.05, 0.63 ± 0.18, and 0.89 ± 0.03 and the mean surface distance (MSD) was 1.4 ± 0.5 mm, 2.0 ± 1.2 mm, and 1.5 ± 0.4 mm, respectively.

Conclusion: With a dosimetric success rate of 78.5–97.9%, this software may facilitate online adaptive IMPT of prostate cancer using a fast, free and open implementation.

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.

Technical Note: The impact of deformable image registration methods on dose warping

PURPOSE

The purpose of this study was to investigate the clinical-relevant discrepancy between doses warped by pure image based deformable image registration (IM-DIR) and by biomechanical model based DIR (BM-DIR) on intensity-homogeneous organs.

METHODS AND MATERIALS

Ten patients (5Head&Neck, 5Prostate) were included. A research DIR tool (ADMRIE_v1.12) was utilized for IM-DIR. After IM-DIR, BM-DIR was carried out for organs (parotids, bladder, and rectum) which often encompass sharp dose gradient. Briefly, high-quality tetrahedron meshes were generated and deformable vector fields (DVF) from IM-DIR were interpolated to the surface nodes of the volume meshes as boundary condition. Then, a FEM solver (ABAQUS_v6.14) was used to simulate the displacement of internal nodes, which were then interpolated to image-voxel grids to get the more physically plausible DVF. Both geometrical and subsequent dose warping discrepancies were quantified between the two DIR methods. Target registration discrepancy(TRD) was evaluated to show the geometry difference. The re-calculated doses on second CT were warped to the pre-treatment CT via two DIR. Clinical-relevant dose parameters and γ passing rate were compared between two types of warped dose. The correlation was evaluated between parotid shrinkage and TRD/dose discrepancy.

RESULT

The parotid shrunk to 75.7% ± 9% of its pre-treatment volume and the percentage of volume with TRD>1.5 mm) was 6.5% ± 4.7%. The normalized mean-dose difference (NMDD) of IM-DIR and BM-DIR was -0.8% ± 1.5%, with range (-4.7% to 1.5%). 2 mm/2% passing rate was 99.0% ± 1.4%. A moderate correlation was found between parotid shrinkage and TRD and NMDD. The bladder had a NMDD of -9.9% ± 9.7%, with BM-DIR warped dose systematically higher. Only minor deviation was observed for rectum NMDD (0.5% ± 1.1%).

CONCLUSION

Impact of DIR method on treatment dose warping is patient and organ-specific. Generally, intensity-homogeneous organs, which undergo larger deformation/shrinkage during treatment and encompass sharp dose gradient, will have greater dose warping uncertainty. For these organs, BM-DIR could be beneficial to the evaluation of DIR/dose-warping uncertainty.

Implementing Radiation Dose-Volume Liver Response in Biomechanical Deformable Image Registration

PURPOSE

Understanding anatomic and functional changes in the liver resulting from radiation therapy is fundamental to the improvement of normal tissue complication probability models needed to advance personalized medicine. The ability to link pretreatment and posttreatment imaging is often compromised by significant dose-dependent volumetric changes within the liver that are currently not accounted for in deformable image registration (DIR) techniques. This study investigated using delivered dose, in combination with other patient factors, to biomechanically model longitudinal changes in liver anatomy for follow-up care and re-treatment planning.

METHODS AND MATERIALS

Population models describing the relationship between dose and hepatic volume response were produced using retrospective data from 33 patients treated with focal radiation therapy. A DIR technique was improved by implementing additional boundary conditions associated with the dose-volume response in series with a previously developed biomechanical DIR algorithm. Evaluation of this DIR technique was performed on computed tomography imaging from 7 patients by comparing the model-predicted volumetric change within the liver with the observed change, tracking vessel bifurcations within the liver through the deformation process, and then determining target registration error between the predicted and identified posttreatment bifurcation points.

RESULTS

Evaluation of the proposed DIR technique showed that all lobes were volumetrically deformed to within the respective contour variability of each lobe. The average target registration error achieved was 7.3 mm (2.8 mm left-right and anterior-posterior and 5.1 mm superior-inferior), with the superior-inferior component within the average limiting slice thickness (6.0 mm). This represented a significant improvement (P<.01, Wilcoxon test) over the application of the previously published biomechanical DIR algorithm (10.9 mm).

CONCLUSIONS

This study demonstrates the feasibility of implementing dose-driven volumetric response in deformable registration, enabling improved accuracy of modeling liver anatomy changes, which could allow for improved dose accumulation, particularly for patients who require additional liver radiation therapy.

Roles of Deformable Image Registration in adaptive RT: From contour propagation to dose monitoring

Adaptive radiation therapy (ART) is based on the optimization of the treatment plan during the treatment delivery to compensate for anatomical deformations. Deformable Image Registration (DIR) then constitutes a key step in order to analyze the huge amount of daily or weekly images to provide clinically usefull information. Two main applications of DIR have been developped in ART: delineation propagation and dose accumulation. If delineation propagation is well validated and transfered in the clinic, some challenges remain to address for dose accumulation. In this paper, we review the recent developments of DIR in ART, particularly in prostate and head-and-neck (H&N), with a focus on their evaluation.

Optimal adaptive IMRT strategy to spare the parotid glands in oropharyngeal cancer

PURPOSE

In oropharyngeal cancer adaptive radiation therapy (ART), this study aimed to quantify the dosimetric benefit of numerous replanning strategies, defined by various numbers and timings of replannings, with regard to parotid gland (PG) sparing.

MATERIAL AND METHODS

Thirteen oropharyngeal cancer patients had one planning and then six weekly CT scans during the seven weeks of IMRT. Weekly doses were recalculated without replanning or with replanning to spare the PG. Sixty-three ART scenarios were simulated by considering all the combinations of numbers and timings of replanning. The PG cumulated doses corresponding to "standard" IMRT and ART scenarios were estimated and compared, either by calculating the average of weekly doses or using deformable image registration (DIR).

RESULTS

Considering average weekly doses, the mean PG overdose using standard IMRT, compared to the planned dose, was 4.1Gy. The mean dosimetric benefit of 6 replannings was 3.3Gy. Replanning at weeks 1, 1-5, 1-2-5, 1-2-4-5 and 1-2-4-5-6 produced the lowest PG mean doses, 94% of the maximum benefit being obtained with 3 replannings. The percentage of patients who had a benefit superior to 5Gy for the contralateral PG was 31% for the three-replannings strategy. The same conclusions were found using DIR.

CONCLUSION

Early replannings proved the most beneficial for PG sparing, three replannings (weeks 1-2-5), representing an attractive combination for ART in oropharyngeal cancer.

Phase correlation applied to the 3D registration of CT and CBCT image volumes

PURPOSE

In this study, a 3D phase correlation algorithm was investigated to test feasibility for use in determining the anatomical changes that occur throughout a patient's radiotherapy treatment. The algorithm determines the transformations between two image volumes through analysis in the Fourier domain and has not previously been used in radiotherapy for 3D registration of CT and CBCT volumes.

METHODS

Various known transformations were applied to a patient's prostate CT image volume to create 12 different test cases. The mean absolute error and standard deviation were determined by evaluating the difference between the known contours and those calculated from the registration process on a point-by-point basis. Similar evaluations were performed on images with increasing levels of noise added. The improvement in structure overlap offered by the algorithm in registering clinical CBCT to CT images was evaluated using the Dice Similarity Coefficient (DSC).

RESULTS

A mean error of 2.35 (σ = 1.54) mm was calculated for the 12 deformations applied. When increasing levels of noise were introduced to the images, the mean errors were observed to rise up to a maximum increase of 1.77 mm. For CBCT to CT registration, maximum improvements in the DSC of 0.09 and 0.46 were observed for the bladder and rectum, respectively.

CONCLUSIONS

The Fourier-based 3D phase correlation registration algorithm investigated displayed promising results in CT to CT and CT to CBCT registration, offers potential in terms of efficiency and robustness to noise, and is suitable for use in radiotherapy for monitoring patient anatomy throughout treatment.

Evaluation of deformable image registration methods for dose monitoring in head and neck radiotherapy

In the context of head and neck cancer (HNC) adaptive radiation therapy (ART), the two purposes of the study were to compare the performance of multiple deformable image registration (DIR) methods and to quantify their impact for dose accumulation, in healthy structures. Fifteen HNC patients had a planning computed tomography (CT0) and weekly CTs during the 7 weeks of intensity-modulated radiation therapy (IMRT). Ten DIR approaches using different registration methods (demons or B-spline free form deformation (FFD)), preprocessing, and similarity metrics were tested. Two observers identified 14 landmarks (LM) on each CT-scan to compute LM registration error. The cumulated doses estimated by each method were compared. The two most effective DIR methods were the demons and the FFD, with both the mutual information (MI) metric and the filtered CTs. The corresponding LM registration accuracy (precision) was 2.44 mm (1.30 mm) and 2.54 mm (1.33 mm), respectively. The corresponding LM estimated cumulated dose accuracy (dose precision) was 0.85 Gy (0.93 Gy) and 0.88 Gy (0.95 Gy), respectively. The mean uncertainty (difference between maximal and minimal dose considering all the 10 methods) to estimate the cumulated mean dose to the parotid gland (PG) was 4.03 Gy (SD = 2.27 Gy, range: 1.06-8.91 Gy).

Impact of head and neck cancer adaptive radiotherapy to spare the parotid glands and decrease the risk of xerostomia

BACKGROUND

Large anatomical variations occur during the course of intensity-modulated radiation therapy (IMRT) for locally advanced head and neck cancer (LAHNC). The risks are therefore a parotid glands (PG) overdose and a xerostomia increase. The purposes of the study were to estimate: - the PG overdose and the xerostomia risk increase during a "standard" IMRT (IMRTstd); - the benefits of an adaptive IMRT (ART) with weekly replanning to spare the PGs and limit the risk of xerostomia.

MATERIAL AND METHODS

Fifteen patients received radical IMRT (70 Gy) for LAHNC. Weekly CTs were used to estimate the dose distributions delivered during the treatment, corresponding either to the initial planning (IMRTstd) or to weekly replanning (ART). PGs dose were recalculated at the fraction, from the weekly CTs. PG cumulated doses were then estimated using deformable image registration. The following PG doses were compared: pre-treatment planned dose, per-treatment IMRTstd and ART. The corresponding estimated risks of xerostomia were also compared. Correlations between anatomical markers and dose differences were searched.

RESULTS

Compared to the initial planning, a PG overdose was observed during IMRTstd for 59% of the PGs, with an average increase of 3.7 Gy (10.0 Gy maximum) for the mean dose, and of 8.2% (23.9% maximum) for the risk of xerostomia. Compared to the initial planning, weekly replanning reduced the PG mean dose for all the patients (p<0.05). In the overirradiated PG group, weekly replanning reduced the mean dose by 5.1 Gy (12.2 Gy maximum) and the absolute risk of xerostomia by 11% (p<0.01) (30% maximum). The PG overdose and the dosimetric benefit of replanning increased with the tumor shrinkage and the neck thickness reduction (p<0.001).

CONCLUSION

During the course of LAHNC IMRT, around 60% of the PGs are overdosed of 4 Gy. Weekly replanning decreased the PG mean dose by 5 Gy, and therefore by 11% the xerostomia risk.