A method to evaluate dose errors introduced by dose mapping processes for mass conserving deformations

Evaluating the accumulated dose distribution of organs at risk in combined radiotherapy for cervical carcinoma based on deformable image registration

OBJECTIVE

To evaluate the feasibility and value of deformable image registration (DIR) in calculating the cumulative doses of organs at risk (OARs) in the combined radiotherapy of cervical cancer.

PATIENTS AND METHODS

Thirty cervical cancer patients treated with external beam radiotherapy (EBRT) combined with intracavitary brachytherapy (ICBT) were reviewed. The simulation CT images of EBRT and ICBT were imported into Varian Velocity 4.1 for the DIR-based dose accumulation. Cumulative dose-volume parameters of D2cc for rectum and bladder were compared between the direct addition (DA) and DIR methods. The quantitative parameters were measured to evaluate the accuracy of DIR.

RESULTS

The three-dimensional cumulative dose distribution of the tumor and OARs were graphically well illustrated by composite isodose lines. In combined EBRT and ICBT, the mean cumulative bladder D2cc calculated by DIR and DA was 86.13 Gy and 86.27 Gy, respectively. The mean cumulative rectal D2cc calculated by DIR and DA was 72.97 Gy and 73.90 Gy, respectively. No significant differences were noted between these two methods (p > 0.05). As to the parameters used to evaluate the DIR accuracy, the mean DSC, Jacobian, MDA (mm) and Hausdorff distance (mm) were 0.79, 1.0, 3.84, and 22.01 respectively for the bladder and 0.53, 1.2, 7.31, and 29.58 respectively for the rectum. In this study, the DSC seemed to be slightly lower compared with previous studies.

CONCLUSION

Dose accumulation based on DIR might be an alternative method to illustrate and evaluate the cumulative doses of the OARs in combined radiotherapy for cervical cancer. However, DIR should be used with caution before overcoming the relevant limitations.

MIRSIG position paper: the use of image registration and fusion algorithms in radiotherapy

The report of the American Association of Physicists in Medicine (AAPM) Task Group No. 132 published in 2017 reviewed rigid image registration and deformable image registration (DIR) approaches and solutions to provide recommendations for quality assurance and quality control of clinical image registration and fusion techniques in radiotherapy. However, that report did not include the use of DIR for advanced applications such as dose warping or warping of other matrices of interest. Considering that DIR warping tools are now readily available, discussions were hosted by the Medical Image Registration Special Interest Group (MIRSIG) of the Australasian College of Physical Scientists & Engineers in Medicine in 2018 to form a consensus on best practice guidelines. This position statement authored by MIRSIG endorses the recommendations of the report of AAPM task group 132 and expands on the best practice advice from the ‘Deforming to Best Practice’ MIRSIG publication to provide guidelines on the use of DIR for advanced applications.

Online daily adaptive proton therapy

It is recognized that the use of a single plan calculated on an image acquired some time before the treatment is generally insufficient to accurately represent the daily dose to the target and to the organs at risk. This is particularly true for protons, due to the physical finite range. Although this characteristic enables the generation of steep dose gradients, which is essential for highly conformal radiotherapy, it also tightens the dependency of the delivered dose to the range accuracy. In particular, the use of an outdated patient anatomy is one of the most significant sources of range inaccuracy, thus affecting the quality of the planned dose distribution. A plan should be ideally adapted as soon as anatomical variations occur, ideally online. In this review, we describe in detail the different steps of the adaptive workflow and discuss the challenges and corresponding state-of-the art developments in particular for an online adaptive strategy.

Deformable Registration for Dose Accumulation

As deformable image registration makes its way into the clinical routine, the summation of doses from fractionated treatment regimens to evaluate cumulative doses to targets and healthy tissues is also becoming a frequently utilized tool in the context of image-guided adaptive radiotherapy. Accounting for daily geometric changes using deformable image registration and dose accumulation potentially enables a better understanding of dose-volume-effect relationships, with the goal of translation of this knowledge to personalization of treatment, to further enhance treatment outcomes. Treatment adaptation involving image deformation requires patient-specific quality assurance of the image registration and dose accumulation processes, to ensure that uncertainties in the 3D dose distributions are identified and appreciated from a clinical relevance perspective. While much research has been devoted to identifying and managing the uncertainties associated with deformable image registration and dose accumulation approaches, there are still many unanswered questions. Here, we provide a review of current deformable image registration and dose accumulation techniques, and related clinical application. We also discuss salient issues that need to be deliberated when applying deformable algorithms for dose mapping and accumulation in the context of adaptive radiotherapy and response assessment.

4DMRI-based investigation on the interplay effect for pencil beam scanning proton therapy of pancreatic cancer patients

BACKGROUND

Time-resolved volumetric magnetic resonance imaging (4DMRI) offers the potential to analyze 3D motion with high soft-tissue contrast without additional imaging dose. We use 4DMRI to investigate the interplay effect for pencil beam scanning (PBS) proton therapy of pancreatic cancer and to quantify the dependency of residual interplay effects on the number of treatment fractions.

METHODS

Based on repeated 4DMRI datasets for nine pancreatic cancer patients, synthetic 4DCTs were generated by warping static 3DCTs with 4DMRI deformation vector fields. 4D dose calculations for scanned proton therapy were performed to quantify the interplay effect by CTV coverage (v95) and dose homogeneity (d5/d95) for incrementally up to 28 fractions. The interplay effect was further correlated to CTV motion characteristics. For quality assurance, volume and mass conservation were evaluated by Jacobian determinants and volume-density comparisons.

RESULTS

For the underlying patient cohort with CTV motion amplitudes < 15 mm, we observed significant correlations between CTV motion amplitudes and both the length of breathing cycles and the interplay effect. For individual fractions, tumor underdosage down to v95 = 70% was observed with pronounced dose heterogeneity (d5/d95 = 1.3). For full × 28 fractionated treatments, we observed a mitigation of the interplay effect with increasing fraction numbers. On average, after seven fractions, a CTV coverage with 95-107% of the prescribed dose was reached with sufficient dose homogeneity. For organs at risk, no significant differences were found between the static and accumulated dose plans for 28 fractions.

CONCLUSION

Intrafractional organ motion exhibits a large interplay effect for PBS proton therapy of pancreatic cancer. The interplay effect correlates with CTV motion, but can be mitigated efficiently by fractionation, mainly due to different breathing starting phases in fractionated treatments. For hypofractionated treatments, a further restriction of motion may be required. Repeated 4DMRI measurements are a viable tool for pre- and post-treatment evaluations of the interplay effect.

Validation of a dose warping algorithm using clinically realistic scenarios

Objective:

Dose warping following deformable image registration (DIR) has been proposed for interfractional dose accumulation. Robust evaluation workflows are vital to clinically implement such procedures. This study demonstrates such a workflow and quantifies the accuracy of a commercial DIR algorithm for this purpose under clinically realistic scenarios.

Methods:

12 head and neck (H&N) patient data sets were used for this retrospective study. For each case, four clinically relevant anatomical changes have been manually generated. Dose distributions were then calculated on each artificially deformed image and warped back to the original anatomy following DIR by a commercial algorithm. Spatial registration was evaluated by quantitative comparison of the original and warped structure sets, using conformity index and mean distance to conformity (MDC) metrics. Dosimetric evaluation was performed by quantitative comparison of the dose–volume histograms generated for the calculated and warped dose distributions, which should be identical for the ideal “perfect” registration of mass-conserving deformations.

Results:

Spatial registration of the artificially deformed image back to the planning CT was accurate (MDC range of 1–2 voxels or 1.2–2.4 mm). Dosimetric discrepancies introduced by the DIR were low (0.02 ± 0.03 Gy per fraction in clinically relevant dose metrics) with no statistically significant difference found (Wilcoxon test, 0.6 ≥ p ≥ 0.2).

Conclusion:

The reliability of CT-to-CT DIR-based dose warping and image registration was demonstrated for a commercial algorithm with H&N patient data.

Advances in knowledge:

This study demonstrates a workflow for validation of dose warping following DIR that could assist physicists and physicians in quantifying the uncertainties associated with dose accumulation in clinical scenarios.

Assessing cumulative dose distributions in combined radiotherapy for cervical cancer using deformable image registration with pre-imaging preparations

Background

The purpose of the study was to evaluate the feasibility of deformable image registration (DIR) in assessing cumulative dose distributions of the combination of external beam radiotherapy (EBRT) and fractionated intracavitary brachytherapy (ICBT) for cervical cancer.

Materials and methods

Three-dimensional image data sets of five consecutive patients were used. The treatment plan consisted of whole pelvic EBRT (total dose: 45 Gy in 25 fractions) combined with computed tomography (CT)-based high-dose rate ICBT (≥24 Gy in 4 fractions to the high risk clinical target volume (HR-CTV)). Organs at risk and HR-CTV were contoured on each CT images and dose-volume parameters were acquired. Pre-imaging preparations were performed prior to each ICBT to minimize the uncertainty of the organ position. Physical doses of each treatment were converted to biologically equivalent doses in 2 Gy daily fractions by the linear quadratic model. Three-dimensional dose distributions of each treatment were accumulated on CT images of the first ICBT using DIR with commercially available image registration software (MIM Maestro®). To compare with DIR, 3D dose distributions were fused by rigid registration based on bony structure matching. To evaluate the accuracy of DIR, the Dice similarity coefficient (DSC) was measured between deformed contours and initial contours.

Results

The cumulative dose distributions were successfully illustrated on the CT images using DIR. Mean DSCs of the HR-CTV, rectum, and bladder were 0.46, 0.62 and 0.69, respectively, with rigid registration; and 0.78, 0.76, and 0.87, respectively, with DIR (p <0.05). The mean DSCs derived from our DIR procedure were comparable to those of previous reports describing the quality of DIR algorithms in the pelvic region. DVH parameters derived from the 2 methods showed no significant difference.

Conclusions

Our results suggest that DIR-based dose accumulation may be acceptable for assessing cumulative dose distributions to assess doses to the tumor and organs at risk in combined radiotherapy for cervical cancer under pre-imaging preparations.

The effect on dose accumulation accuracy of inverse-consistency and transitivity error reduced deformation maps

It has previously been shown that deformable image registrations (DIRs) often result in deformation maps that are neither inverse-consistent nor transitive, and that the dose accumulation based on these deformation maps can be inconsistent if different image pathways are used for dose accumulation. A method presented to reduce inverse consistency and transitivity errors has been shown to result in more consistent dose accumulation, regardless of the image pathway selected for dose accumulation. The present study investigates the effect on the dose accumulation accuracy of deformation maps processed to reduce inverse consistency and transitivity errors. A set of lung 4DCT phases were analysed, consisting of four images on which a dose grid was created. Dose to 75 corresponding anatomical locations was manually tracked. Dose accumulation was performed between all image sets with Demons derived deformation maps as well as deformation maps processed to reduce inverse consistency and transitivity errors. The ground truth accumulated dose was then compared with the accumulated dose derived from DIR. Two dose accumulation image pathways were considered. The post-processing method to reduce inverse consistency and transitivity errors had minimal effect on the dose accumulation accuracy. There was a statistically significant improvement in dose accumulation accuracy for one pathway, but for the other pathway there was no statistically significant difference. A post-processing technique to reduce inverse consistency and transitivity errors has a positive, yet minimal effect on the dose accumulation accuracy. Thus the post-processing technique improves consistency of dose accumulation with minimal effect on dose accumulation accuracy.

Dose mapping sensitivity to deformable registration uncertainties in fractionated radiotherapy - applied to prostate proton treatments

BACKGROUND

Calculation of accumulated dose in fractionated radiotherapy based on spatial mapping of the dose points generally requires deformable image registration (DIR). The accuracy of the accumulated dose thus depends heavily on the DIR quality. This motivates investigations of how the registration uncertainty influences dose planning objectives and treatment outcome predictions.A framework was developed where the dose mapping can be associated with a variable known uncertainty to simulate the DIR uncertainties in a clinical workflow. The framework enabled us to study the dependence of dose planning metrics, and the predicted treatment outcome, on the DIR uncertainty. The additional planning margin needed to compensate for the dose mapping uncertainties can also be determined. We applied the simulation framework to a hypofractionated proton treatment of the prostate using two different scanning beam spot sizes to also study the dose mapping sensitivity to penumbra widths.

RESULTS

The planning parameter most sensitive to the DIR uncertainty was found to be the target D95. We found that the registration mean absolute error needs to be ≤0.20 cm to obtain an uncertainty better than 3% of the calculated D95 for intermediate sized penumbras. Use of larger margins in constructing PTV from CTV relaxed the registration uncertainty requirements to the cost of increased dose burdens to the surrounding organs at risk.

CONCLUSIONS

The DIR uncertainty requirements should be considered in an adaptive radiotherapy workflow since this uncertainty can have significant impact on the accumulated dose. The simulation framework enabled quantification of the accuracy requirement for DIR algorithms to provide satisfactory clinical accuracy in the accumulated dose.

A reference dataset for deformable image registration spatial accuracy evaluation using the COPDgene study archive

Landmark point-pairs provide a strategy to assess deformable image registration (DIR) accuracy in terms of the spatial registration of the underlying anatomy depicted in medical images. In this study, we propose to augment a publicly available database (www.dir-lab.com) of medical images with large sets of manually identified anatomic feature pairs between breath-hold computed tomography (BH-CT) images for DIR spatial accuracy evaluation. Ten BH-CT image pairs were randomly selected from the COPDgene study cases. Each patient had received CT imaging of the entire thorax in the supine position at one-fourth dose normal expiration and maximum effort full dose inspiration. Using dedicated in-house software, an imaging expert manually identified large sets of anatomic feature pairs between images. Estimates of inter- and intra-observer spatial variation in feature localization were determined by repeat measurements of multiple observers over subsets of randomly selected features. 7298 anatomic landmark features were manually paired between the 10 sets of images. Quantity of feature pairs per case ranged from 447 to 1172. Average 3D Euclidean landmark displacements varied substantially among cases, ranging from 12.29 (SD: 6.39) to 30.90 (SD: 14.05) mm. Repeat registration of uniformly sampled subsets of 150 landmarks for each case yielded estimates of observer localization error, which ranged in average from 0.58 (SD: 0.87) to 1.06 (SD: 2.38) mm for each case. The additions to the online web database (www.dir-lab.com) described in this work will broaden the applicability of the reference data, providing a freely available common dataset for targeted critical evaluation of DIR spatial accuracy performance in multiple clinical settings. Estimates of observer variance in feature localization suggest consistent spatial accuracy for all observers across both four-dimensional CT and COPDgene patient cohorts.