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

Accuracy of contour propagation from planning computed tomography to iterative cone-beam computed tomography using a deformable image registration algorithm for assisting head and neck radiotherapy

BACKGROUND

The deteriorated image quality of cone-beam computed tomography (CBCT) reduces the accuracy of contour propagation. We investigated the accuracy of contour propagation from planning CT (pCT) to iterative CBCT (iCBCT) using deformable image registration and compared it with that of replanning CT (reCT) and Feldkamp-Davis-Kress algorithm-based CBCT (FDK-CBCT). No report exists regarding iCBCT improving the accuracy of this technique for images of the head and neck region.

METHODS

We included 29 patients who underwent radiotherapy for head and neck cancer. ReCT and CBCT were performed on the same day. The gross tumor volume (GTV) and organs at risk, including the brain stem, spinal cord, mandible, parotid glands, submandibular glands, and larynx, were manually contoured by radiation oncologists on pCT and reCT images. Contour propagation was performed using MIM software. Manually delineated contours on reCT images and deformably generated contours on reCT, FDK-CBCT, and iCBCT images were compared to determine the accuracy of contour propagation using the Dice similarity coefficient (DSC), mean distance to agreement (MDA), and Hausdorff distance (HD).

RESULTS

The mean DSC values for all contoured organs were 0.84 across reCT, FDK-CBCT and iCBCT. A mean DSC value of >0.8 was observed for all organs, except for the larynx and the GTV. The MDA was <1.5 mm for all organs and images, whereas the HD value showed a variation of >3 mm.

CONCLUSION

The results demonstrated no statistically significant difference in contour propagation from pCT to iCBCT compared to pCT to reCT.

Accuracy Evaluation of Dose Warping Using Deformable Image Registration in Carbon Ion Therapy

Purpose

The study's purpose was to use a simple geometry phantom to validate the deformable image registration (DIR) accuracy and dose warping accuracy in carbon ion radiotherapy (CIRT) and to provide an index for dosimetry in CIRT.

Materials and methods

We used geometric and anatomical phantoms provided by AAPM TG-132. The DIRs of 3 different settings were performed between reference and translational images for each phantom. CIRT and photon therapy (3-dimensional conformal radiotherapy and volumetric modulated radiotherapy) treatment plans were transformed by the use of deformation vector fields calculated from the DIR of each setting. The dose distribution calculated on the basis of rigid registration between images was used as the ground truth and compared with the warped dose determined by DIR to evaluate the error.

Results

The photon therapy treatment plans showed a dose warping error of <2% for a DIR error of <2 mm, whereas CIRT showed a dose warping error >10% for the same DIR accuracy. From this, even with similar DIR accuracy, the errors tended to be larger for CIRT than for photon therapy dose distributions. Due to the steepness of the CIRT dose gradient, the dose difference increased by about 5% for a DIR error of 5 mm, which was larger than the 3% dose difference generated for a 5 mm DIR error in photon therapy.

Conclusion

We evaluated the relationship between DIR accuracy and dose warping accuracy in the CIRT dose distribution. Due to the steepness of the dose gradient in CIRT, we concluded that dose warping based on DIR accuracy should be required to be sufficiently higher than that in photon therapy.

Interpreting the patterns of local failure following postoperative volumetric-modulated arctherapy in oral cavity and oropharynx cancers: Impact of the different methods of analysis

PURPOSE

Intensity-modulated radiation therapy or volumetric-modulated arctherapy is nowadays the recommended radiation technique for the treatment of head and neck cancers. However, by providing a significant dose gradient between target volumes and organs at risk, there is a risk of target missing and thus recurrence in case of inadequate delineation. It is therefore necessary to determine the origin of these recurrences to improve clinical practice. Over the past years, different methods have been described for the analysis of recurrences. Using the patterns of failure of patients with oral cavity and oropharynx carcinoma, treated with postoperative volumetric-modulated arctherapy in our institution, the purpose of this work was to analyse the sites of local recurrences and to evaluate the disparity in the classification of recurrences when different methods were used.

MATERIAL AND METHODS

Between 2011 and 2019, 167 patients who underwent postoperative volumetric-modulated arctherapy for oral cavity or oropharyngeal cancers were included (60 and 40 % respectively). Two or three dose levels were prescribed (54Gy, 59.4/60Gy±66Gy). Local recurrence occurred in 17 patients (10.2 %). We assessed the patterns of local recurrences according to four methods: 1/ volume-based method using the volume overlap between the recurrence volume and initial target volumes; 2/ volume-based method of overlap between the recurrence volume and the 95 % treatment isodose; 3/ point-based method using the position of the barycentre of the recurrence volume; 4/ combined centroid method classifying recurrences according to both the initial target volumes and dose distribution. Each case was reviewed to make a clinical judgment on these classifications and assessed them as "appropriate", "possible", or "inappropriate".

RESULTS

For the volume-based method using overlap between the recurrence volume and the initial clinical target volume, this classification was clinically judged as inappropriate in 11 out of 17 cases (65 %). For the volume-based method using overlap between the recurrence volume and the 95 % prescribed isodose, this classification was clinically judged as appropriate in 15 out of 17 cases (88 %). For the point-based method, this classification was clinically judged as appropriate in 14 out of 17 cases (82 %). Thirteen out of 17 local recurrences had the same classification between this point-based method and the volume-based method of overlap between the recurrence volume and the 95 % prescribed isodose. For the combined centroid method, among 17 local recurrences nine were classified as type A, two as type B, two as type C, three as type D and one as type E. This classification was clinically judged as appropriate in 15 out of 17 cases (88 %). Only five out of 17 of the local recurrences were classified the same way according to the four different methods (29 %).

CONCLUSION

Recurrences that are "marginal" or "outfield" represent a major challenge for intensity-modulated radiation therapy/volumetric-modulated arctherapy quality assurance and improvement of delineation recommendations. To date, there are no published methods that give complete satisfaction.

A unified generation-registration framework for improved MR-based CT synthesis in proton therapy

BACKGROUND

The use of magnetic resonance (MR) imaging for proton therapy treatment planning is gaining attention as a highly effective method for guidance. At the core of this approach is the generation of computed tomography (CT) images from MR scans. However, the critical issue in this process is accurately aligning the MR and CT images, a task that becomes particularly challenging in frequently moving body areas, such as the head-and-neck. Misalignments in these images can result in blurred synthetic CT (sCT) images, adversely affecting the precision and effectiveness of the treatment planning.

PURPOSE

This study introduces a novel network that cohesively unifies image generation and registration processes to enhance the quality and anatomical fidelity of sCTs derived from better-aligned MR images.

METHODS

The approach synergizes a generation network (G) with a deformable registration network (R), optimizing them jointly in MR-to-CT synthesis. This goal is achieved by alternately minimizing the discrepancies between the generated/registered CT images and their corresponding reference CT counterparts. The generation network employs a UNet architecture, while the registration network leverages an implicit neural representation (INR) of the displacement vector fields (DVFs). We validated this method on a dataset comprising 60 head-and-neck patients, reserving 12 cases for holdout testing.

RESULTS

Compared to the baseline Pix2Pix method with MAE 124.95 ± $\pm$ 30.74 HU, the proposed technique demonstrated 80.98 ± $\pm$ 7.55 HU. The unified translation-registration network produced sharper and more anatomically congruent outputs, showing superior efficacy in converting MR images to sCTs. Additionally, from a dosimetric perspective, the plan recalculated on the resulting sCTs resulted in a remarkably reduced discrepancy to the reference proton plans.

CONCLUSIONS

This study conclusively demonstrates that a holistic MR-based CT synthesis approach, integrating both image-to-image translation and deformable registration, significantly improves the precision and quality of sCT generation, particularly for the challenging body area with varied anatomic changes between corresponding MR and CT.

Deformable anthropomorphic pelvis phantom for dose accumulation verification

Objective.The validation of deformable image registration (DIR) for contour propagation is often done using contour-based metrics. Meanwhile, dose accumulation requires evaluation of voxel mapping accuracy, which might not be accurately represented by contour-based metrics. By fabricating a deformable anthropomorphic pelvis phantom, we aim to (1) quantify the voxel mapping accuracy for various deformation scenarios, in high- and low-contrast regions, and (2) identify any correlation between dice similarity coefficient (DSC), a commonly used contour-based metric, and the voxel mapping accuracy for each organ.Approach. Four organs, i.e. pelvic bone, prostate, bladder and rectum (PBR), were 3D printed using PLA and a Polyjet digital material, and assembled. The latter three were implanted with glass bead and CT markers within or on their surfaces. Four deformation scenarios were simulated by varying the bladder and rectum volumes. For each scenario, nine DIRs with different parameters were performed on RayStation v10B. The voxel mapping accuracy was quantified by finding the discrepancy between true and mapped marker positions, termed the target registration error (TRE). Pearson correlation test was done between the DSC and mean TRE for each organ.Main results. For the first time, we fabricated a deformable phantom purely from 3D printing, which successfully reproduced realistic anatomical deformations. Overall, the voxel mapping accuracy dropped with increasing deformation magnitude, but improved when more organs were used to guide the DIR or limit the registration region. DSC was found to be a good indicator of voxel mapping accuracy for prostate and rectum, but a comparatively poorer one for bladder. DSC > 0.85/0.90 was established as the threshold of mean TRE ⩽ 0.3 cm for rectum/prostate. For bladder, extra metrics in addition to DSC should be considered.Significance. This work presented a 3D printed phantom, which enabled quantification of voxel mapping accuracy and evaluation of correlation between DSC and voxel mapping accuracy.

A multi-center survey on adaptive radiation therapy for head and neck cancer in Malaysia

The survey is to assess the current state of adaptive radiation therapy (ART) for head and neck (H&N) cases among radiotherapy centers in Malaysia and to identify any implementation limitations. An online questionnaire was sent to all radiotherapy centers in Malaysia. The 24-question questionnaire consists of general information about the center, ART practices, and limitations faced in implementing ART. 28 out of 36 radiotherapy centers responded, resulting in an overall response rate of 78%. About 52% of the responding centers rescanned and replanned less than 5% of their H&N patients. The majority (88.9%) of the respondents reported the use Cone Beam Computed Tomography alone or in combination with other modalities to trigger the ART process. The main reasons cited for adopting ART were weight loss, changes in the immobilization fitting, and anatomical variation. The adaptation process typically occurred during week 3 or week 4 of treatment. More than half of the respondents require three days or more from re-simulation to starting a new treatment plan. Both target and organ at risk delineation on new planning CT relied heavily on manual delineation by physicians and physicists, respectively. All centers perform patient-specific quality assurance for their new adaptive plans. Two main limitations in implementing ART are "limited financial resources or equipment" and "limitation on technical knowledge". There is a need for a common consensus to standardize the practice of ART and address these limitations to improve the implementation of ART in Malaysia.

Practical and technical key challenges in head and neck adaptive radiotherapy: The GORTEC point of view

Anatomical variations occur during head and neck (H&N) radiotherapy (RT) treatment. These variations may result in underdosage to the target volume or overdosage to the organ at risk. Replanning during the treatment course can be triggered to overcome this issue. Due to technological, methodological and clinical evolutions, tools for adaptive RT (ART) are becoming increasingly sophisticated. The aim of this paper is to give an overview of the key steps of an H&N ART workflow and tools from the point of view of a group of French-speaking medical physicists and physicians (from GORTEC). Focuses are made on image registration, segmentation, estimation of the delivered dose of the day, workflow and quality assurance for an implementation of H&N offline and online ART. Practical recommendations are given to assist physicians and medical physicists in a clinical workflow.

Feasibility of delivered dose reconstruction for MR-guided SBRT of pancreatic tumors with fast, real-time 3D cine MRI

BACKGROUND AND PURPOSE

In MR-guided SBRT of pancreatic cancer, intrafraction motion is typically monitored with (interleaved) 2D cine MRI. However, tumor surroundings are often not fully captured in these images, and motion might be distorted by through-plane movement. In this study, the feasibility of highly accelerated 3D cine MRI to reconstruct the delivered dose during MR-guided SBRT was assessed.

MATERIALS AND METHODS

A 3D cine MRI sequence was developed for fast, time-resolved 4D imaging, featuring a low spatial resolution that allows for rapid volumetric imaging at 430 ms. The 3D cines were acquired during the entire beam-on time of 23 fractions of online adaptive MR-guided SBRT for pancreatic tumors on a 1.5 T MR-Linac. A 3D deformation vector field (DVF) was extracted for every cine dynamic using deformable image registration. Next, these DVFs were used to warp the partial dose delivered in the time interval between consecutive cine acquisitions. The warped dose plans were summed to obtain a total delivered dose. The delivered dose was also calculated under various motion correction strategies. Key DVH parameters of the GTV, duodenum, small bowel and stomach were extracted from the delivered dose and compared to the planned dose. The uncertainty of the calculated DVFs was determined with the inverse consistency error (ICE) in the high-dose regions.

RESULTS

The mean (SD) relative (ratio delivered/planned) D_99% of the GTV was 0.94 (0.06), and the mean (SD) relative D_0.5cc of the duodenum, small bowel, and stomach were respectively 0.98 (0.04), 1.00 (0.07), and 0.98 (0.06). In the fractions with the lowest delivered tumor coverage, it was found that significant lateral drifts had occurred. The DVFs used for dose warping had a low uncertainty with a mean (SD) ICE of 0.65 (0.07) mm.

CONCLUSION

We employed a fast, real-time 3D cine MRI sequence for dose reconstruction in the upper abdomen, and demonstrated that accurate DVFs, acquired directly from these images, can be used for dose warping. The reconstructed delivered dose showed only a modest degradation of tumor coverage, mostly attainable to baseline drifts. This emphasizes the need for motion monitoring and development of intrafraction treatment adaptation solutions, such as baseline drift corrections.

Dose evaluations of organs at risk and predictions of gastrointestinal toxicity after re-irradiation with stereotactic body radiation therapy for pancreatic cancer by deformable image registration

Purpose

Re-irradiation of locally recurrent pancreatic cancer may be an optimal choice as a local ablative therapy. However, dose constraints of organs at risk (OARs) predictive of severe toxicity remain unknown. Therefore, we aim to calculate and identify accumulated dose distributions of OARs correlating with severe adverse effects and determine possible dose constraints regarding re-irradiation.

Methods

Patients receiving two courses of stereotactic body radiation therapy (SBRT) for the same irradiated regions (the primary tumors) due to local recurrence were included. All doses of the first and second plans were recalculated to an equivalent dose of 2 Gy per fraction (EQD₂). Deformable image registration with the workflow "Dose Accumulation-Deformable" of the MIM^® System (version: 6.6.8) was performed for dose summations. Dose-volume parameters predictive of grade 2 or more toxicities were identified, and the receiver operating characteristic (ROC) curve was used to determine optimal thresholds of dose constraints.

Results

Forty patients were included in the analysis. Only the V ₁₀ of the stomach [hazard ratio (HR): 1.02 (95% CI:1.00-1.04), P = 0.035] and D _mean of the intestine [HR: 1.78 (95% CI: 1.00-3.18), P = 0.049] correlated with grade 2 or more gastrointestinal toxicity. Hence, the equation of probability of such toxicity was P = 1 1 + e - ( - 4.155 + 0.579 D mean of the intestine + 0.021 V 10  of the stomach ) Additionally, the area under the ROC curve and threshold of dose constraints of V ₁₀ of the stomach and D _mean of the intestine were 0.779 and 77.575 cc, 0.769 and 4.22 Gy₃ (α/β = 3), respectively. The area under the ROC curve of the equation was 0.821.

Conclusion

The V ₁₀ of the stomach and D _mean of the intestine may be vital parameters to predict grade 2 or more gastrointestinal toxicity, of which the threshold of dose constraints may be beneficial for the practice of re-irradiation of locally relapsed pancreatic cancer.

Dose accumulation for MR-guided adaptive radiotherapy: From practical considerations to state-of-the-art clinical implementation

MRI-linear accelerator (MR-linac) devices have been introduced into clinical practice in recent years and have enabled MR-guided adaptive radiation therapy (MRgART). However, by accounting for anatomical changes throughout radiation therapy (RT) and delivering different treatment plans at each fraction, adaptive radiation therapy (ART) highlights several challenges in terms of calculating the total delivered dose. Dose accumulation strategies-which typically involve deformable image registration between planning images, deformable dose mapping, and voxel-wise dose summation-can be employed for ART to estimate the delivered dose. In MRgART, plan adaptation on MRI instead of CT necessitates additional considerations in the dose accumulation process because MRI pixel values do not contain the quantitative information used for dose calculation. In this review, we discuss considerations for dose accumulation specific to MRgART and in relation to current MR-linac clinical workflows. We present a general dose accumulation framework for MRgART and discuss relevant quality assurance criteria. Finally, we highlight the clinical importance of dose accumulation in the ART era as well as the possible ways in which dose accumulation can transform clinical practice and improve our ability to deliver personalized RT.

Incorporating the synthetic CT image for improving the performance of deformable image registration between planning CT and cone-beam CT

Objective

To develop a contrast learning-based generative (CLG) model for the generation of high-quality synthetic computed tomography (sCT) from low-quality cone-beam CT (CBCT). The CLG model improves the performance of deformable image registration (DIR).

Methods

This study included 100 post-breast-conserving patients with the pCT images, CBCT images, and the target contours, which the physicians delineated. The CT images were generated from CBCT images via the proposed CLG model. We used the Sct images as the fixed images instead of the CBCT images to achieve the multi-modality image registration accurately. The deformation vector field is applied to propagate the target contour from the pCT to CBCT to realize the automatic target segmentation on CBCT images. We calculate the Dice similarity coefficient (DSC), 95 % Hausdorff distance (HD95), and average surface distance (ASD) between the prediction and reference segmentation to evaluate the proposed method.

Results

The DSC, HD95, and ASD of the target contours with the proposed method were 0.87 ± 0.04, 4.55 ± 2.18, and 1.41 ± 0.56, respectively. Compared with the traditional method without the synthetic CT assisted (0.86 ± 0.05, 5.17 ± 2.60, and 1.55 ± 0.72), the proposed method was outperformed, especially in the soft tissue target, such as the tumor bed region.

Conclusion

The CLG model proposed in this study can create the high-quality sCT from low-quality CBCT and improve the performance of DIR between the CBCT and the pCT. The target segmentation accuracy is better than using the traditional DIR.

Evaluation of DIR algorithm performance in real patients for radiotherapy treatments: A systematic review of operator-dependent strategies

PURPOSE

The performance of deformable medical image registration (DIR) algorithms has become a major concern.

METHODS

We aimed to obtain updated information on DIR algorithm performance quantification through a literature review of articles published between 2010 and 2022. We focused only on studies using operator-based methods to treat real patients. The PubMed, Google Scholar and Embase databases were searched following PRISMA guidelines.

RESULTS

One hundred and seven articles were identified. The mean number of patients and registrations per publication was 20 and 63, respectively. We found 23 different geometric metrics appearing at least twice, and the dosimetric impact of DIR was quantified in 32 articles. Forty-eight different at-risk organs were described, and target volumes were studied in 43 publications. Prostate, head-and-neck and thoracic locations represented more than ¾ of the studied locations. We summarized the type of DIR and the images used, and other key elements. Intra/interobserver variability, threshold values and the correlation between metrics were also discussed.

CONCLUSIONS

This literature review covers the past decade and should facilitate the implementation of DIR algorithms in clinical practice by providing practical and pertinent information to quantify their performance on real patients.

Comparison of an in-house hybrid DIR method to NiftyReg on CBCT and CT images for head and neck cancer

An in-house hybrid deformable image registration (DIR) method, which combines free-form deformation (FFD) and the viscous fluid registration method, is proposed. Its results on the planning computed tomography (CT) and the day 1 treatment cone-beam CT (CBCT) image from 68 head and neck cancer patients are compared with the results of NiftyReg, which uses B-spline FFD alone. Several similarity metrics, the target registration error (TRE) of annotated points, as well as the Dice similarity coefficient (DSC) and Hausdorff distance (HD) of the propagated organs at risk are employed to analyze their registration accuracy. According to quantitative analysis on mutual information, normalized cross-correlation, and the absolute pixel value differences, the results of the proposed DIR are more similar to the CBCT images than the NiftyReg results. Smaller TRE of the annotated points is observed in the proposed method, and the overall mean TRE for the proposed method and NiftyReg was 2.34 and 2.98 mm, respectively (p < 0.001). The mean DSC in the larynx, spinal cord, oral cavity, mandible, and parotid given by the proposed method ranged from 0.78 to 0.91, significantly higher than the NiftyReg results (ranging from 0.77 to 0.90), and the HD was significantly lower compared to NiftyReg. Furthermore, the proposed method did not suffer from unrealistic deformations as the NiftyReg did in the visual evaluation. Meanwhile, the execution time of the proposed method was much higher than NiftyReg (96.98 ± 11.88 s vs. 4.60 ± 0.49 s). In conclusion, the in-house hybrid method gave better accuracy and more stable performance than NiftyReg.

Self-channel-and-spatial-attention neural network for automated multi-organ segmentation on head and neck CT images

Accurate segmentation of organs at risk (OARs) is necessary for adaptive head and neck (H&N) cancer treatment planning, but manual delineation is tedious, slow, and inconsistent. A self-channel-and-spatial-attention neural network (SCSA-Net) is developed for H&N OAR segmentation on CT images. To simultaneously ease the training and improve the segmentation performance, the proposed SCSA-Net utilizes the self-attention ability of the network. Spatial and channel-wise attention learning mechanisms are both employed to adaptively force the network to emphasize the meaningful features and weaken the irrelevant features simultaneously. The proposed network was first evaluated on a public dataset, which includes 48 patients, then on a separate serial CT dataset, which contains ten patients who received weekly diagnostic fan-beam CT scans. On the second dataset, the accuracy of using SCSA-Net to track the parotid and submandibular gland volume changes during radiotherapy treatment was quantified. The Dice similarity coefficient (DSC), positive predictive value (PPV), sensitivity (SEN), average surface distance (ASD), and 95% maximum surface distance (95SD) were calculated on the brainstem, optic chiasm, optic nerves, mandible, parotid glands, and submandibular glands to evaluate the proposed SCSA-Net. The proposed SCSA-Net consistently outperforms the state-of-the-art methods on the public dataset. Specifically, compared with Res-Net and SE-Net, which is constructed from squeeze-and-excitation block equipped residual blocks, the DSC of the optic nerves and submandibular glands is improved by 0.06, 0.03 and 0.05, 0.04 by the SCSA-Net. Moreover, the proposed method achieves statistically significant improvements in terms of DSC on all and eight of nine OARs over Res-Net and SE-Net, respectively. The trained network was able to achieve good segmentation results on the serial dataset, but the results were further improved after fine-tuning of the model using the simulation CT images. For the parotids and submandibular glands, the volume changes of individual patients are highly consistent between the automated and manual segmentation (Pearson's correlation 0.97-0.99). The proposed SCSA-Net is computationally efficient to perform segmentation (sim 2 s/CT).

Computed Tomography/Magnetic Resonance Imaging (CT/MRI) Image Registration and Fusion Assessment for Accurate Glioblastoma Radiotherapy Treatment Planning

Background: In this study, computed tomography/magnetic resonance imaging (CT/MRI) image registration and fusion in the 3D conformal radiotherapy treatment planning of Glioblastoma brain tumor was investigated. Good CT/MRI image registration and fusion made a great impact on dose calculation and treatment planning accuracy. Indeed, the uncertainly associated with the registration and fusion methods must be well verified and communicated. Unfortunately, there is no standard procedure or mathematical formalism to perform this verification due to noise, distortion, and complicated anatomical situations. Objectives: This study aimed at assessing the effective contribution of MRI in Glioma radiotherapy treatment by improving the localization of target volumes and organs at risk (OARs). It is also a question to provide clinicians with some suitable metrics to evaluate the CT/MRI image registration and fusion results. Methods: Quantitative image registration and fusion evaluation were used in this study to compare Eclipse TPS tools and Elastix CT/MRI image registration fusion. Thus, Dice score coefficient (DSC), Jaccard similarity coefficient (JSC), and Hausdorff distance (HD) were found to be suitable metrics for the evaluation and comparison of the image registration and fusion methods of Eclipse TPS and Elastix. Results: The programmed tumor’s volumes (PTV) delineated on CT slices were approximately 1.38 times smaller than those delineated on CT/MRI fused images. Large differences were observed for the edema and the brainstem. It was also found that MRI considerably optimized the dose to be delivered to the optic nerve and brainstem. Conclusions: Image registration and fusion is a fundamental step for suitable and efficient Glioma treatment planning in 3D conformal radiotherapy that ensure accurate dose delivery and unnecessary OAR irradiation. MRI can provide accurate localization of targeted volumes leading to better irradiation control of Glioma tumor.

Mixed-beam approach for high-risk prostate cancer: Carbon-ion boost followed by photon intensity-modulated radiotherapy. Dosimetric and geometric evaluations (AIRC IG-14300)

BACKGROUND AND PURPOSE

The aim was to evaluate dosimetric uncertainties of a mixed beam approach for patients with high-risk prostate cancer (PCa). The treatment consists of a carbon ion radiotherapy (CIRT) boost followed by whole-pelvis intensity-modulated RT (IMRT).

MATERIALS AND METHODS

Patients were treated with a CIRT boost of 16.6 Gy/4 fractions followed by whole-pelvis IMRT of 50 Gy/25 fractions, with consequent long term androgen deprivation therapy. Deformable computed tomography image registration (DIR) was performed and corresponding doses were used for plan sum. A comparative IMRT photon plan was obtained as whole-pelvis IMRT of 50 Gy/25 fractions followed by a boost of 28 Gy/14 fractions. DIR performances were evaluated through structure-related and image characteristics parameters.

RESULTS

Until now, five patients out of ten total enrolled ended the treatment. Dosimetric parameters were lower in CIRT + IMRT than IMRT-only plans for all organs at risk (OARs) except femoral heads. Regarding DIR evaluation, femoral heads were the less deformed OAR. Penile bulb, bladder and anal canal showed intermediate deformation. Rectum was the most deformed. DIR algorithms were patient (P)-dependent, as performances were the highest for P3 and P4, intermediate for P2 and P5, and the lowest for P1.

CONCLUSIONS

CIRT allows better OARs sparing while increasing the efficacy due to the higher radio-biological effect of carbon ions. However, a mixed beam approach could introduce DIR problems in multi-centric treatments with different operative protocols. The development of this prospective trial will lead to more mature data concerning the clinical impact of implementing DIR procedures in dose accumulation applications for high-risk PCa treatments.

Free-to-use DIR solutions in radiotherapy: Benchmark against commercial platforms through a contour-propagation study

PURPOSE

A contour propagation study has been conducted to benchmark three algorithms for Deformable Image Registration (DIR) freely available online against well-established commercial solutions.

METHODS

ElastiX, BRAINS and Plastimach, available as modules in the open source platform 3DSlicer, were tested as the recent AAPM Task group 132 guidelines proposes. The overlap of the DIR-mapped ROIs in four computational anthropomorphic phantoms was measured. To avoid bias every algorithm was left to run without any human interaction nor particular registration strategy. The accuracy of the algorithms was measured using the Dice Similarity Coefficient (DSC) and Mean Distance to Conformity (MDC) metrics. The registration quality was compared to the recommended geometrical accuracy suggested by AAPM TG132 and to the results of a large population-based study performed with commercial DIR solutions.

RESULTS

The considered free-to-use DIR solutions generally meet acceptable accuracy and good overlap (DSC > 0.85). Mild failures (DSC < 0.75) were detected only for the smallest structures. In case of extremely severe deformations acceptable accuracy was not met (MDC > 3 mm). The morphing capability of the tested algorithms equals those of commercial systems when the user interaction is avoided. Underperformances were detected only in cases where a specific registration strategy is mandatory to obtain a satisfying match.

CONCLUSIONS

All of the considered algorithms show performances not inferior to previously published data and have the potential to be good candidates for use in the clinical routine. The results and conclusions only apply to the considered phantoms and should not be considered to be generally applicable and extendable to patient cases.

Quantifying the dose accumulation uncertainty after deformable image registration in head-and-neck radiotherapy

BACKGROUND AND PURPOSE

Deformable image registration (DIR) facilitated dose reconstruction and accumulation can be applied to assess delivered dose and verify the validity of the treatment plan during treatment. This retrospective study used in silico deformations based on clinically observed anatomical changes as ground truth to investigate the uncertainty of reconstructed and accumulated dose in head-and-neck radiotherapy (HNRT).

MATERIALS AND METHODS

A planning CT (pCT), cone beam CT (CBCT) from week one of treatment and three later CBCTs were selected for 12 HNRT patients. These images were used to generate in silico reference CBCTs and deformation vector fields (DVFs) as ground truth with B-spline DIR. Inverse consistency (IC) of voxels was assessed by determining their net displacement after successive application of the forward and backward DVF. The reconstructed dose based on demons DIR was compared to the ground truth to assess the structure-specific uncertainties of this DIR algorithm for inverse consistent and inverse inconsistent voxels.

RESULTS

Overall, 98.5% of voxels were inverse consistent with the 95% level of confidence range for dose reconstruction of a single fraction equal to [-2.3%; +2.1%], [-10.2%; +15.2%] and [-9.5%; +12.5%] relative to their planned dose for target structures, critical organs at risk (OARs) and non-critical OARs, respectively. Inverse inconsistent voxels generally showed a higher level of uncertainty.

CONCLUSION

The uncertainty in accumulated dose using DIR can be accurately quantified and incorporated in dose-volume histograms (DVHs). This method can be used to prospectively assess the adequacy of target coverage during treatment in an objective manner.

"Dose of the day" based on cone beam computed tomography and deformable image registration for lung cancer radiotherapy

Purpose

Adaptive radiotherapy (ART) has potential to reduce toxicity and facilitate safe dose escalation. Dose calculations with the planning CT deformed to cone beam CT (CBCT) have shown promise for estimating the “dose of the day”. The purpose of this study is to investigate the “dose of the day” calculation accuracy based on CBCT and deformable image registration (DIR) for lung cancer radiotherapy.

Methods

A total of 12 lung cancer patients were identified, for which daily CBCT imaging was performed for treatment positioning. A re‐planning CT (rCT) was acquired after 20 Gy for all patients. A virtual CT (vCT) was created by deforming initial planning CT (pCT) to the simulated CBCT that was generated from deforming CBCT to rCT acquired on the same day. Treatment beams from the initial plan were copied to the vCT and rCT for dose calculation. Dosimetric agreement between vCT‐based and rCT‐based accumulated doses was evaluated using the Bland‐Altman analysis.

Results

Mean differences in dose‐volume metrics between vCT and rCT were smaller than 1.5%, and most discrepancies fell within the range of ± 5% for the target volume, lung, esophagus, and heart. For spinal cord D_max, a large mean difference of −5.55% was observed, which was largely attributed to very limited CBCT image quality (e.g., truncation artifacts).

Conclusion

This study demonstrated a reasonable agreement in dose‐volume metrics between dose accumulation based on vCT and rCT, with the exception for cases with poor CBCT image quality. These findings suggest potential utility of vCT for providing a reasonable estimate of the “dose of the day”, and thus facilitating the process of ART for lung cancer.

Computed Tomography to Cone Beam Computed Tomography Deformable Image Registration for Contour Propagation Using Head and Neck, Patient-Based Computational Phantoms: A Multicenter Study

PURPOSE

To investigate the performance of various algorithms for deformable image registration (DIR) for propagating regions of interest (ROIs) using multiple commercial platforms, from computed tomography to cone beam computed tomography (CBCT) and megavoltage computed tomography.

METHODS AND MATERIALS

Fourteen institutions participated in the study using 5 commercial platforms: RayStation (RaySearch Laboratories, Stockholm, Sweden), MIM (Cleveland, OH), VelocityAI and SmartAdapt (Varian Medical Systems, Palo Alto, CA), and ABAS (Elekta AB, Stockholm, Sweden). Algorithms were tested on synthetic images generated with the ImSimQA (Oncology Systems Limited, Shrewsbury, UK) package by applying 2 specific deformation vector fields (DVF) to real head and neck patient datasets. On-board images from 3 systems were used: megavoltage computed tomography from Tomotherapy and 2 kinds of CBCT from a clinical linear accelerator. Image quality of the system was evaluated. The algorithms' accuracy was assessed by comparing the DIR-mapped ROIs returned by each center with those of the reference, using the Dice similarity coefficient and mean distance to conformity metrics. Statistical inference on the validation results was carried out to identify the prognostic factors of DIR performance.

RESULTS

Analyzing 840 DIR-mapped ROIs returned by the centers, it was demonstrated that DVF intensity and image quality were significant prognostic factors of DIR performance. The accuracy of the propagated contours was generally high, and acceptable DIR performance can be obtained with lower-dose CBCT image protocols.

CONCLUSIONS

The performance of the systems proved to be image quality specific, depending on the DVF type and only partially on the platforms. All systems proved to be robust against image artifacts and noise, except the demon-based software.

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.

Effect of Deformation Methods on the Accuracy of Deformable Image Registration From Kilovoltage CT to Tomotherapy Megavoltage CT

INTRODUCTION

The registration accuracy of megavoltage computed tomography images is limited by low image contrast when compared to that of kilovoltage computed tomography images. Such issues may degrade the deformable image registration accuracy. This study evaluates the deformable image registration from kilovoltage to megavoltage images when using different deformation methods and assessing nasopharyngeal carcinoma patient images.

METHODS

The kilovoltage and the megavoltage images from the first day and the 20th fractions of the treatment day of 12 patients with nasopharyngeal carcinoma were used to evaluate the deformable image registration application. The deformable image registration image procedures were classified into 3 groups, including kilovoltage to kilovoltage, megavoltage to megavoltage, and kilovoltage to megavoltage. Three deformable image registration methods were employed using the deformable image registration and adaptive radiotherapy software. The validation was compared by volume-based, intensity-based, and deformation field analyses.

RESULTS

The use of different deformation methods greatly affected the deformable image registration accuracy from kilovoltage to megavoltage. The asymmetric transformation with the demon method was significantly better than other methods and illustrated satisfactory value for adaptive applications. The deformable image registration accuracy from kilovoltage to megavoltage showed no significant difference from the kilovoltage to kilovoltage images when using the appropriate method of registration.

CONCLUSIONS

The choice of deformation method should be considered when applying the deformable image registration from kilovoltage to megavoltage images. The deformable image registration accuracy from kilovoltage to megavoltage revealed a good agreement in terms of intensity-based, volume-based, and deformation field analyses and showed clinically useful methods for nasopharyngeal carcinoma adaptive radiotherapy in tomotherapy applications.

A density assignment method for dose monitoring in head-and-neck radiotherapy

BACKGROUND AND PURPOSE

During head-and-neck (H&N) radiotherapy, the parotid glands (PGs) may be overdosed; thus, a tool is required to monitor the delivered dose. This study aimed to assess the dose accuracy of a patient-specific density assignment method (DAM) for dose calculation to monitor the dose to PGs during treatment.

PATIENTS AND METHODS

Forty patients with H&N cancer received an intensity modulated radiation therapy (IMRT), among whom 15 had weekly CTs. Dose distributions were calculated either on the CTs (CT_ref), on one-class CTs (1C-CT, water), or on three-class CTs (3C-CT, water-air-bone). The inter- and intra-patient DAM uncertainties were evaluated by the difference between doses calculated on CT_ref and 1C-CTs or 3C-CTs. PG mean dose (D_mean) and spinal cord maximum dose (D_2%) were considered. The cumulated dose to the PGs was estimated by the mean D_mean of the weekly CTs.

RESULTS

The mean (maximum) inter-patient DAM dose uncertainties for the PGs (in cGy) were 23 (75) using 1C-CTs and 12 (50) using 3C-CTs (p ≤ 0.001). For the spinal cord, these uncertainties were 118 (245) and 15 (67; p ≤ 0.001). The mean (maximum) DAM dose uncertainty between cumulated doses calculated on CTs and 3C-CTs was 7 cGy (45 cGy) for the PGs. Considering the difference between the planned and cumulated doses, 53% of the ipsilateral and 80% of the contralateral PGs were overdosed by +3.6 Gy (up to 8.2 Gy) and +1.9 Gy (up to 5.2 Gy), respectively.

CONCLUSION

The uncertainty of the three-class DAM appears to be clinically non-significant (<0.5 Gy) compared with the PG overdose (up to 8.2 Gy). This DAM could therefore be used to monitor PG doses and trigger replanning.

Anatomically plausible models and quality assurance criteria for online mono- and multi-modal medical image registration

Medical imaging is currently employed in the diagnosis, planning, delivery and response monitoring of cancer treatments. Due to physiological motion and/or treatment response, the shape and location of the pathology and organs-at-risk may change over time. Establishing their location within the acquired images is therefore paramount for an accurate treatment delivery and monitoring. A feasible solution for tracking anatomical changes during an image-guided cancer treatment is provided by image registration algorithms. Such methods are, however, often built upon elements originating from the computer vision/graphics domain. Since the original design of such elements did not take into consideration the material properties of particular biological tissues, the anatomical plausibility of the estimated deformations may not be guaranteed. In the current work we adapt two existing variational registration algorithms, namely Horn-Schunck and EVolution, to online soft tissue tracking. This is achieved by enforcing an incompressibility constraint on the estimated deformations during the registration process. The existing and the modified registration methods were comparatively tested against several quality assurance criteria on abdominal in vivo MR and CT data. These criteria included: the Dice similarity coefficient (DSC), the Jaccard index, the target registration error (TRE) and three additional criteria evaluating the anatomical plausibility of the estimated deformations. Results demonstrated that both the original and the modified registration methods have similar registration capabilities in high-contrast areas, with DSC and Jaccard index values predominantly in the 0.8-0.9 range and an average TRE of 1.6-2.0 mm. In contrast-devoid regions of the liver and kidneys, however, the three additional quality assurance criteria have indicated a considerable improvement of the anatomical plausibility of the deformations estimated by the incompressibility-constrained methods. Moreover, the proposed registration models maintain the potential of the original methods for online image-based guidance of cancer treatments.

Synthetic head and neck and phantom images for determining deformable image registration accuracy in magnetic resonance imaging

PURPOSE

Magnetic resonance imaging (MRI) provides noninvasive evaluation of patient's anatomy without using ionizing radiation. Due to this and the high soft-tissue contrast, MRI use has increased and has potential for use in longitudinal studies where changes in patients' anatomy or tumors at different time points are compared. Deformable image registration can be useful for these studies. Here, we describe two datasets that can be used to evaluate the registration accuracy of systems for MR images, as it cannot be assumed to be the same as that measured on CT images.

ACQUISITION AND VALIDATION METHODS

Two sets of images were created to test registration accuracy. (a) A porcine phantom was created by placing ten 0.35-mm gold markers into porcine meat. The porcine phantom was immobilized in a plastic container with movable dividers. T1-weighted, T2-weighted, and CT images were acquired with the porcine phantom compressed in four different ways. The markers were not visible on the MR images, due to the selected voxel size, so they did not interfere with the measured registration accuracy, while the markers were visible on the CT images and were used to identify the known deformation between positions. (b) Synthetic images were created using 28 head and neck squamous cell carcinoma patients who had MR scans pre-, mid-, and postradiotherapy treatment. An inter- and intrapatient variation model was created using these patient scans. Four synthetic pretreatment images were created using the interpatient model, and four synthetic post-treatment images were created for each of the pretreatment images using the intrapatient model.

DATA FORMAT AND USAGE NOTES

The T1-weighted, T2-weighted, and CT scans of the porcine phantom in the four positions are provided. Four T1-weighted synthetic pretreatment images each with four synthetic post-treatment images, and four T2-weighted synthetic pretreatment images each with four synthetic post-treatment images are provided. Additionally, the applied deformation vector fields to generate the synthetic post-treatment images are provided. The data are available on TCIA under the collection MRI-DIR.

POTENTIAL APPLICATIONS

The proposed database provides two sets of images (one acquired and one computer generated) for use in evaluating deformable image registration accuracy. T1- and T2-weighted images are available for each technique as the different image contrast in the two types of images may impact the registration accuracy.

The impact of robustness of deformable image registration on contour propagation and dose accumulation for head and neck adaptive radiotherapy

Deformable image registration (DIR) is the key process for contour propagation and dose accumulation in adaptive radiation therapy (ART). However, currently, ART suffers from a lack of understanding of “robustness” of the process involving the image contour based on DIR and subsequent dose variations caused by algorithm itself and the presetting parameters. The purpose of this research is to evaluate the DIR caused variations for contour propagation and dose accumulation during ART using the RayStation treatment planning system. Ten head and neck cancer patients were selected for retrospective studies. Contours were performed by a single radiation oncologist and new treatment plans were generated on the weekly CT scans for all patients. For each DIR process, four deformation vector fields (DVFs) were generated to propagate contours and accumulate weekly dose by the following algorithms: (a) ANACONDA with simple presetting parameters, (b) ANACONDA with detailed presetting parameters, (c) MORFEUS with simple presetting parameters, and (d) MORFEUS with detailed presetting parameters. The geometric evaluation considered DICE coefficient and Hausdorff distance. The dosimetric evaluation included D₉₅, D_max, D_mean, D_min, and Homogeneity Index. For geometric evaluation, the DICE coefficient variations of the GTV were found to be 0.78 ± 0.11, 0.96 ± 0.02, 0.64 ± 0.15, and 0.91 ± 0.03 for simple ANACONDA, detailed ANACONDA, simple MORFEUS, and detailed MORFEUS, respectively. For dosimetric evaluation, the corresponding Homogeneity Index variations were found to be 0.137 ± 0.115, 0.006 ± 0.032, 0.197 ± 0.096, and 0.006 ± 0.033, respectively. The coherent geometric and dosimetric variations also consisted in large organs and small organs. Overall, the results demonstrated that the contour propagation and dose accumulation in clinical ART were influenced by the DIR algorithm, and to a greater extent by the presetting parameters. A quality assurance procedure should be established for the proper use of a commercial DIR for adaptive radiation therapy.

Inter-patient image registration algorithms to disentangle regional dose bioeffects

Radiation therapy (RT) technological advances call for a comprehensive reconsideration of the definition of dose features leading to radiation induced morbidity (RIM). In this context, the voxel-based approach (VBA) to dose distribution analysis in RT offers a radically new philosophy to evaluate local dose response patterns, as an alternative to dose-volume-histograms for identifying dose sensitive regions of normal tissue. The VBA relies on mapping patient dose distributions into a single reference case anatomy which serves as anchor for local dosimetric evaluations. The inter-patient elastic image registrations (EIRs) of the planning CTs provide the deformation fields necessary for the actual warp of dose distributions. In this study we assessed the impact of EIR on the VBA results in thoracic patients by identifying two state-of-the-art EIR algorithms (Demons and B-Spline). Our analysis demonstrated that both the EIR algorithms may be successfully used to highlight subregions with dose differences associated with RIM that substantially overlap. Furthermore, the inclusion for the first time of covariates within a dosimetric statistical model that faces the multiple comparison problem expands the potential of VBA, thus paving the way to a reliable voxel-based analysis of RIM in datasets with strong correlation of the outcome with non-dosimetric variables.

Variation assessment of deformable registration in stereotactic radiosurgery

INTRODUCTION

The regular functions of CT-MRI registration include delineation of targets and organs-at-risk (OARs) in radiosurgery planning. The question of whether deformable image registration (DIR) could be applied to stereotactic radiosurgery (SRS) in its place remains a subject of debate.

METHODS

This study collected data regarding 16 patients who had undergone single-fraction SRS treatment. All lesions were located close to the brainstem. CT and MRI two image sets were registered by both rigid image registration (RIR) and DIR algorithms. The contours of the OARs were drawn individually on the rigid and deformable CT-MRI image sets by qualified radiation oncologists and dosimetrists. The evaluation metrics included volume overlapping (VO), Dice similarity coefficient (DSC), and dose. The modified demons deformable algorithm (VARIAN SmartAdapt) was used for evaluation in this study.

RESULTS

The mean range of VO for OARs was 0.84 ± 0.08, and DSC was 0.82 ± 0.07. The maximum average volume difference was at normal brain (17.18 ± 14.48 cm³) and the second highest was at brainstem (2.26 cm³ ± 1.18). Pearson correlation testing showed that all DIRs' OAR volumes were linearly and significantly correlated with RIRs' volume (0.679-0.992, two tailed, P << 0.001). The 100% dose was prescribed at gross tumor volume (GTV). The average maximum percent dose difference was observed in brainstem (26.54% ± 27.027), and the average mean dose difference has found at same organ (1.6% ± 1.66).

CONCLUSION

The change in image-registration method definitely produces dose variance, and is significantly more what depending on the target location. The volume size of OARs, however, was not statistical significantly correlated with dose variance.

Evaluation of Deformable Image Registration (DIR) Methods for Dose Accumulation in Nasopharyngeal Cancer Patients during Radiotherapy

Introduction

Deformable image registration (DIR) is used to modify structures according to anatomical changes for observing the dosimetric effect. In this study, megavoltage computed tomography (MVCT) images were used to generate cumulative doses for nasopharyngeal cancer (NPC) patients by various DIR methods. The performance of the multiple DIR methods was analysed, and the impact of dose accumulation was assessed.

Patients and methods

The study consisted of five NPC patients treated with a helical tomotherapy unit. The weekly MVCT images at the 1^st, 6^th, 11^th, 16^th, 21^st, 26^th, and 31^st fractions were used to assess the dose accumulation by the four DIR methods. The cumulative dose deviations from the initial treatment plan were analysed, and correlations of these variations with the anatomic changes and DIR methods were explored.

Results

The target dose received a slightly different result from the initial plan at the end of the treatment. The organ dose differences increased as the treatment progressed to 6.8% (range: 2.2 to 10.9%), 15.2% (range: -1.7 to 36.3%), and 6.4% (range: -1.6 to 13.2%) for the right parotid, the left parotid, and the spinal cord, respectively. The mean uncertainty values to estimate the accumulated doses for all the DIR methods were 0.21 ± 0.11 Gy (target dose), 1.99 ± 0.76 Gy (right parotid), 1.19 ± 0.24 Gy (left parotid), and 0.41 ± 0.04 Gy (spinal cord).

Conclusions

Accuracy of the DIR methods affects the estimation of dose accumulation on both the target dose and the organ dose. The DIR methods provide an adequate dose estimation technique for observation as a result of inter-fractional anatomic changes and are beneficial for adaptive treatment strategies.

Representing the dosimetric impact of deformable image registration errors

Deformable image registration (DIR) is emerging as a tool in radiation therapy for calculating the cumulative dose distribution across multiple fractions of treatment. Unfortunately, due to the variable nature of DIR algorithms and dependence of performance on image quality, registration errors can result in dose accumulation errors. In this study, landmarked images were used to characterize the DIR error throughout an image space and determine its impact on dosimetric analysis. Ten thoracic 4DCT images with 300 landmarks per image study matching the end-inspiration and end-expiration phases were obtained from 'dir-labs'. DIR was performed using commercial software MIM Maestro. The range of dose uncertainty (RDU) was calculated at each landmark pair as the maximum and minimum of the doses within a sphere around the landmark in the end-expiration phase. The radius of the sphere was defined by a measure of DIR error which included either the actual DIR error, mean DIR error per study, constant errors of 2 or 5 mm, inverse consistency error, transitivity error or the distance discordance metric (DDM). The RDUs were evaluated using the magnitude of dose uncertainty (MDU) and inclusion rate (IR) of actual error lying within the predicted RDU. The RDU was calculated for 300 landmark pairs on each 4DCT study for all measures of DIR error. The most representative RDU was determined using the actual DIR error with a MDU of 2.5 Gy and IR of 97%. Across all other measures of DIR error, the DDM was most predictive with a MDU of 2.5 Gy and IR of 86%, closest to the actual DIR error. The proposed method represents the range of dosimetric uncertainty of DIR error using either landmarks at specific voxels or measures of registration accuracy throughout the volume.

Voxel-based analysis unveils regional dose differences associated with radiation-induced morbidity in head and neck cancer patients

The risk of radiation-induced toxicity in patients treated for head and neck (HN) cancer with radiation therapy (RT) is traditionally estimated by condensing the 3D dose distribution into a monodimensional cumulative dose-volume histogram which disregards information on dose localization. We hypothesized that a voxel-based approach would identify correlations between radiation-induced morbidity and local dose release, thus providing a new insight into spatial signature of radiation sensitivity in composite regions like the HN district. This methodology was applied to a cohort of HN cancer patients treated with RT at risk of radiation-induced acute dysphagia (RIAD). We implemented an inter-patient elastic image registration framework that proved robust enough to match even the most elusive HN structures and to provide accurate dose warping. A voxel-based statistical analysis was then performed to test regional dosimetric differences between patients with and without RIAD. We identified a significantly higher dose delivered to RIAD patients in two voxel clusters in correspondence of the cricopharyngeus muscle and cervical esophagus. Our study goes beyond the well-established organ-based philosophy exploring the relationship between radiation-induced morbidity and local dose differences in the HN region. This approach is generally applicable to different HN toxicity endpoints and is not specific to RIAD.

Accuracy of eight deformable image registration (DIR) methods for tomotherapy megavoltage computed tomography (MVCT) images

Introduction

The application of deformable image registration (DIR) to megavoltage computed tomography (MVCT) images benefits adaptive radiotherapy. This study aims to quantify the accuracy of DIR for MVCT images when using different deformation methods assessed in a cubic phantom and nasopharyngeal carcinoma (NPC) patients.

Methods

In the control studies, the DIR accuracy in air‐tissue and tissue‐tissue interface areas was observed using twelve shapes of acrylic and tissue‐equivalent material inserted in the phantom. In the clinical studies, the 1st and 20th fraction MVCT images of seven NPC patients were used to evaluate application of DIR. The eight DIR methods used in the DIRART software varied in (i) transformation framework (asymmetric or symmetric), (ii) DIR registration algorithm (Demons or Optical Flow) and (iii) mapping direction (forward or backward). The accuracy of the methods was compared using an intensity‐based criterion (correlation coefficient, CC) and volume‐based criterion (Dice's similarity coefficient, DSC).

Results

The asymmetric transformation with Optical Flow showed the best performance for air‐tissue interface areas, with a mean CC and DSC of 0.97 ± 0.03 and 0.79 ± 0.11 respectively. The symmetric transformation with Optical Flow showed good agreement for tissue‐tissue interface areas with a CC of (0.99 ± 0.01) and DSC of (0.89 ± 0.03). The sequences of target domains were significantly different in tissue‐tissue interface areas.

Conclusions

The deformation method and interface area affected the accuracy of DIR. The validation techniques showed satisfactory volume matching of greater than 0.7 with DSC analysis. The methods can yield acceptable results for clinical applications.

Accuracy of deformable image registration on magnetic resonance images in digital and physical phantoms

PURPOSE

Accurate deformable image registration is necessary for longitudinal studies. The error associated with commercial systems has been evaluated using computed tomography (CT). Several in-house algorithms have been evaluated for use with magnetic resonance imaging (MRI), but there is still relatively little information about MRI deformable image registration. This work presents an evaluation of two deformable image registration systems, one commercial (Velocity) and one in-house (demons-based algorithm), with MRI using two different metrics to quantify the registration error.

METHODS

The registration error was analyzed with synthetic MR images. These images were generated from interpatient and intrapatient variation models trained on 28 patients. Four synthetic post-treatment images were generated for each of four synthetic pretreatment images, resulting in 16 image registrations for both the T1- and T2-weighted images. The synthetic post-treatment images were registered to their corresponding synthetic pretreatment image. The registration error was calculated between the known deformation vector field and the generated deformation vector field from the image registration system. The registration error was also analyzed using a porcine phantom with ten implanted 0.35-mm diameter gold markers. The markers were visible on CT but not MRI. CT, T1-weighted MR, and T2-weighted MR images were taken in four different positions. The markers were contoured on the CT images and rigidly registered to their corresponding MR images. The MR images were deformably registered and the distance between the projected marker location and true marker location was measured as the registration error.

RESULTS

The synthetic images were evaluated only on Velocity. Root mean square errors (RMSEs) of 0.76 mm in the left-right (LR) direction, 0.76 mm in the anteroposterior (AP) direction, and 0.69 mm in the superior-inferior (SI) direction were observed for the T1-weighted MR images. RMSEs of 1.1 mm in the LR direction, 0.75 mm in the AP direction, and 0.81 mm in the SI direction were observed for the T2-weighted MR images. The porcine phantom MR images, when evaluated with Velocity, had RMSEs of 1.8, 1.5, and 2.7 mm in the LR, AP, and SI directions for the T1-weighted images and 1.3, 1.2, and 1.6 mm in the LR, AP, and SI directions for the T2-weighted images. When the porcine phantom images were evaluated with the in-house demons-based algorithm, RMSEs were 1.2, 1.5, and 2.1 mm in the LR, AP, and SI directions for the T1-weighted images and 0.81, 1.1, and 1.1 mm in the LR, AP, and SI directions for the T2-weighted images.

CONCLUSIONS

The MRI registration error was low for both Velocity and the in-house demons-based algorithm according to both image evaluation methods, with all RMSEs below 3 mm. This implies that both image registration systems can be used for longitudinal studies using MRI.

Three-dimensional finite-element based deformable image registration for evaluation of pleural cavity irradiation during photodynamic therapy

PURPOSE

Photodynamic therapy (PDT) is used after surgical resection to treat the microscopic disease for malignant pleural mesothelioma and to increase survival rates. As accurate light delivery is imperative to PDT efficacy, the deformation of the pleural volume during the surgery is studied on its impact on the delivered light fluence. In this study, a three-dimensional finite element-based (3D FEM) deformable image registration is proposed to directly match the volume of lung to the volume of pleural cavity obtained during PDT to have accurate representation of the light fluence accumulated in the lung, heart and liver (organs-at-risk) during treatment.

METHODS

A wand, comprised of a modified endotrachial tube filled with Intralipid and an optical fiber inside the tube, is used to deliver the treatment light. The position of the treatment is tracked using an optical tracking system with an attachment comprised of nine reflective passive markers that are seen by an infrared camera-based navigation system. This information is used to obtain the surface contours of the plural cavity and the cumulative light fluence on every point of the cavity surface that is being treated. The lung, heart, and liver geometry are also reconstructed from a series of computed tomography (CT) scans of the organs acquired in the same patient before and after the surgery. The contours obtained with the optical tracking system and CTs are imported into COMSOL Multiphysics, where the 3D FEM-based deformable image registration is obtained. The delivered fluence values are assigned to the respective positions (x, y, and z) on the optical tracking contour. The optical tracking contour is considered as the reference, and the CT contours are used as the target, which will be deformed. The data from three patients formed the basis for this study.

RESULTS

The physical correspondence between the CT and optical tracking geometries, taken at different times, from different imaging devices was established using the 3D FEM-based image deformable registration. The volume of lung was matched to the volume of pleural cavity and the distribution of light fluence on the surface of the heart, liver and deformed lung volumes was obtained.

CONCLUSION

The method used is appropriate for analyzing problems over complicated domains, such as when the domain changes (as in a solid-state reaction with a moving boundary), when the desired precision varies over the entire domain, or when the solution lacks smoothness. Implementing this method in real-time for clinical applications and in situ monitoring of the under- or over- exposed regions to light during PDT can significantly improve the treatment for mesothelioma.

Evaluation of deformable image registration for contour propagation between CT and cone-beam CT images in adaptive head and neck radiotherapy

Deformable image registration (DIR) is a critical technic in adaptive radiotherapy (ART) to propagate contours between planning computerized tomography (CT) images and treatment CT/Cone-beam CT (CBCT) image to account for organ deformation for treatment re-planning. To validate the ability and accuracy of DIR algorithms in organ at risk (OAR) contours mapping, seven intensity-based DIR strategies are tested on the planning CT and weekly CBCT images from six Head & Neck cancer patients who underwent a 6 ∼ 7 weeks intensity-modulated radiation therapy (IMRT). Three similarity metrics, i.e. the Dice similarity coefficient (DSC), the percentage error (PE) and the Hausdorff distance (HD), are employed to measure the agreement between the propagated contours and the physician delineated ground truths. It is found that the performance of all the evaluated DIR algorithms declines as the treatment proceeds. No statistically significant performance difference is observed between different DIR algorithms (p> 0.05), except for the double force demons (DFD) which yields the worst result in terms of DSC and PE. For the metric HD, all the DIR algorithms behaved unsatisfactorily with no statistically significant performance difference (p= 0.273). These findings suggested that special care should be taken when utilizing the intensity-based DIR algorithms involved in this study to deform OAR contours between CT and CBCT, especially for those organs with low contrast.

Comprehensive evaluation of ten deformable image registration algorithms for contour propagation between CT and cone-beam CT images in adaptive head & neck radiotherapy

Deformable image registration (DIR) is a critical technic in adaptive radiotherapy (ART) for propagating contours between planning computerized tomography (CT) images and treatment CT/cone-beam CT (CBCT) images to account for organ deformation for treatment re-planning. To validate the ability and accuracy of DIR algorithms in organ at risk (OAR) contour mapping, ten intensity-based DIR strategies, which were classified into four categories—optical flow-based, demons-based, level-set-based and spline-based—were tested on planning CT and fractional CBCT images acquired from twenty-one head & neck (H&N) cancer patients who underwent 6~7-week intensity-modulated radiation therapy (IMRT). Three similarity metrics, i.e., the Dice similarity coefficient (DSC), the percentage error (PE) and the Hausdorff distance (HD), were employed to measure the agreement between the propagated contours and the physician-delineated ground truths of four OARs, including the vertebra (VTB), the vertebral foramen (VF), the parotid gland (PG) and the submandibular gland (SMG). It was found that the evaluated DIRs in this work did not necessarily outperform rigid registration. DIR performed better for bony structures than soft-tissue organs, and the DIR performance tended to vary for different ROIs with different degrees of deformation as the treatment proceeded. Generally, the optical flow-based DIR performed best, while the demons-based DIR usually ranked last except for a modified demons-based DISC used for CT-CBCT DIR. These experimental results suggest that the choice of a specific DIR algorithm depends on the image modality, anatomic site, magnitude of deformation and application. Therefore, careful examinations and modifications are required before accepting the auto-propagated contours, especially for automatic re-planning ART systems.

Construction of a biomechanical head and neck motion model as a guide to evaluation of deformable image registration

The use of deformable image registration methods in the context of adaptive radiotherapy leads to uncertainties in the simulation of the administered dose distributions during the treatment course. Evaluation of these methods is a prerequisite to decide if a plan adaptation will improve the individual treatment. Current approaches using manual references limit the validity of evaluation, especially for low-contrast regions. In particular, for the head and neck region, the highly flexible anatomy and low soft tissue contrast in control images pose a challenge to image registration and its evaluation. Biomechanical models promise to overcome this issue by providing anthropomorphic motion modelling of the patient. We introduce a novel biomechanical motion model for the generation and sampling of different postures of the head and neck anatomy. Motion propagation behaviour of the individual bones is defined by an underlying kinematic model. This model interconnects the bones by joints and thus is capable of providing a wide range of motion. Triggered by the motion of the individual bones, soft tissue deformation is described by an extended heterogeneous tissue model based on the chainmail approach. This extension, for the first time, allows the propagation of decaying rotations within soft tissue without the necessity for explicit tissue segmentation. Overall motion simulation and sampling of deformed CT scans including a basic noise model is achieved within 30 s. The proposed biomechanical motion model for the head and neck site generates displacement vector fields on a voxel basis, approximating arbitrary anthropomorphic postures of the patient. It was developed with the intention of providing input data for the evaluation of deformable image registration.

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.

A multiple-image-based method to evaluate the performance of deformable image registration in the pelvis

Deformable image registration (DIR) is essential for adaptive radiotherapy (RT) for tumor sites subject to motion, changes in tumor volume, as well as changes in patient normal anatomy due to weight loss. Several methods have been published to evaluate DIR-related uncertainties but they are not widely adopted. The aim of this study was, therefore, to evaluate intra-patient DIR for two highly deformable organs-the bladder and the rectum-in prostate cancer RT using a quantitative metric based on multiple image registration, the distance discordance metric (DDM). Voxel-by-voxel DIR uncertainties of the bladder and rectum were evaluated using DDM on weekly CT scans of 38 subjects previously treated with RT for prostate cancer (six scans/subject). The DDM was obtained from group-wise B-spline registration of each patient's collection of repeat CT scans. For each structure, registration uncertainties were derived from DDM-related metrics. In addition, five other quantitative measures, including inverse consistency error (ICE), transitivity error (TE), Dice similarity (DSC) and volume ratios between corresponding structures from pre- and post- registered images were computed and compared with the DDM. The DDM varied across subjects and structures; DDMmean of the bladder ranged from 2 to 13 mm and from 1 to 11 mm for the rectum. There was a high correlation between DDMmean of the bladder and the rectum (Pearson's correlation coefficient, R p  =  0.62). The correlation between DDMmean and the volume ratios post-DIR was stronger (R p  =  0.51; 0.68) than the correlation with the TE (bladder: R p  =  0.46; rectum: R p  =  0.47), or the ICE (bladder: R p  =  0.34; rectum: R p  =  0.37). There was a negative correlation between DSC and DDMmean of both the bladder (R p  =  -0.23) and the rectum (R p  =  -0.63). The DDM uncertainty metric indicated considerable DIR variability across subjects and structures. Our results show a stronger correlation with volume ratios and with the DSC using DDM compared to using ICE and TE. The DDM has the potential to quantitatively identify regions of large DIR uncertainties and consequently identify anatomical/scan outliers. The DDM can, thus, be applied to improve the adaptive RT process for tumor sites subject to motion.

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.

A Nomogram to predict parotid gland overdose in head and neck IMRT

Purposes

To generate a nomogram to predict parotid gland (PG) overdose and to quantify the dosimetric benefit of weekly replanning based on its findings, in the context of intensity-modulated radiotherapy (IMRT) for locally-advanced head and neck carcinoma (LAHNC).

Material and methods

Twenty LAHNC patients treated with radical IMRT underwent weekly computed tomography (CT) scans during IMRT. The cumulated PG dose was estimated by elastic registration. Early predictors of PG overdose (cumulated minus planned doses) were identified, enabling a nomogram to be generated from a linear regression model. Its performance was evaluated using a leave-one-out method. The benefit of weekly replanning was then estimated for the nomogram-identified PG overdose patients.

Results

Clinical target volume 70 (CTV70) and the mean PG dose calculated from the planning and first weekly CTs were early predictors of PG overdose, enabling a nomogram to be generated. A mean PG overdose of 2.5Gy was calculated for 16 patients, 14 identified by the nomogram. All patients with PG overdoses >1.5Gy were identified. Compared to the cumulated delivered dose, weekly replanning of these 14 targeted patients enabled a 3.3Gy decrease in the mean PG dose.

Conclusion

Based on the planning and first week CTs, our nomogram allowed the identification of all patients with PG overdoses >2.5Gy to be identified, who then benefitted from a final 4Gy decrease in mean PG overdose by means of weekly replanning.