Validation of a deformable image registration produced by a commercial treatment planning system in head and neck

Experimental validation of a 4D dynamic dose calculation model for proton pencil beam scanning without spot time stamp considering free-breathing motion

PURPOSE

We developed a 4-dimensional dynamic dose (4DDD) calculation model for proton pencil beam scanning (PBS). This model incorporates the spill start time for all energies and uses the remaining irradiated spot time model instead of irradiated spot time logs. This study aimed to validate the calculation accuracy of a log file-based 4DDD model by comparing it with dose measurements performed under free-breathing conditions, thereby serving as an alternative approach to the conventional log file-based system.

METHODS

Three cubic verification plans were created using a heterogeneous block phantom; these plans were created using 10 phase 4D-CT datasets of the phantom. The CIRS dynamic platform was used to simulate motion with amplitudes of 2.5, 3.75, and 5.0 mm. These plans consisted of eight- and two-layered rescanning techniques. The lateral profiles were measured using a 2D ionization chamber array (2D-array) and EBT3 Gafchromic films at four starting phases, including three sinusoidal curves (periods of 3, 4, and 6 s) and a representative patient curve during actual treatment. 4DDDs were calculated using in-house scripting that assigned a time stamp to each spot and performed dose accumulation using deformable image registration. Furthermore, to evaluate the impact of parameter selection on our 4DDD model calculations, simulations were performed assuming a ±10% change in irradiation time stamp (0.8 ± 0.08 s) and spot scan speed. We evaluated the 2D gamma index and the absolute point doses between the calculated values and the measurements.

RESULTS

The 2D-array measurements revealed that the gamma scores for the static plans (no motion) and 4DDD plans exceeded 97.5% and 93.9% at 3%/3 mm, respectively. The average gamma score of the 4DDD plans was at least 96.1%. When using EBT3 films, the gamma scores of the 4DDD model exceeded 92.4% and 98.7% at 2%/2 mm and 3%/3 mm, respectively. Regarding the 4DDD point dose differences, more than 95% of the dose regions exhibited discrepancies within ±5.0% for 97.7% of the total points across all plans. The spot time assignment accuracy of our 4DDD model was acceptable even with ±10% sensitivity. However, the accuracy of the scan speed, when varied within ±10% sensitivity, was not acceptable (minimum gamma scores of 82.6% and maximum dose difference of 12.9%).

CONCLUSIONS

Our 4DDD calculations under free-breathing conditions using amplitudes of less than 5.0 mm were in good agreement with the measurements regardless of the starting phases, breathing curve patterns (between 3 and 6 s periods), and varying numbers of layered rescanning. The proposed system allows us to evaluate actual irradiated doses in various breathing periods, amplitudes, and starting phases, even on PBS machines without the ability to record spot logs.

Cone-beam computed-tomography-based delta-radiomic analysis for investigating prognostic power for esophageal squamous cell cancer patients undergoing concurrent chemoradiotherapy

PURPOSE

To investigate the prognostic power of cone-beam computed-tomography (CBCT)-based delta-radiomics in esophageal squamous cell cancer (ESCC) patients treated with concurrent chemoradiotherapy (CCRT).

METHODS

We collected data from 26 ESCC patients treated with CCRT. CBCT images acquired at five time points (1st-5th week) per patient during CCRT were used in this study. Radiomic features were extracted from the five CBCT images on the gross tumor volumes. Then, 17 delta-radiomic feature sets derived from five types of calculations were obtained for all the cases. Leave-one-out cross-validation was applied to investigate the prognostic power of CBCT-based delta-radiomic features. Feature selection and construction of a prediction model using Coxnet were performed using training samples. Then, the test sample was classified into high or low risk in each cross-validation fold. Survival analysis for the two groups were performed to evaluate the prognostic power of the extracted CBCT-based delta-radiomic features.

RESULTS

Four delta-radiomic feature sets indicated significant differences between the high- and low-risk groups (p < 0.05). The highest C-index in the 17 delta-radiomic feature sets was 0.821 (95 % confidence interval, 0.735-0.907). That feature set had p-value of the log-rank test and hazard ratio of 0.003 and 4.940 (95 % confidence interval, 1.391-17.544), respectively.

CONCLUSIONS

We investigated the potential of using CBCT-based delta-radiomics for prognosis of ESCC patients treated with CCRT. It was demonstrated that delta-radiomic feature sets based on the absolute value of relative difference obtained from the early to the middle treatment stages have high prognostic power for ESCC.

Assessment of residual geometrical errors of clinical target volumes and their impact on dose accumulation for head and neck radiotherapy

PURPOSE

To assess the residual geometrical errors (dr) and their impact on the clinical target volumes (CTV) dose coverage for head and neck cancer (HNC) proton therapy patients.

METHODS

We analysed 28 HNC patients treated with 70 Gy (RBE) and 54.25 Gy (RBE) to the therapeutic CTV70 and prophylactic CTV54.25, respectively. Daily cone beam CTs were converted to high quality synthetic CTs (sCTs). The CTVs from the nominal CT were propagated to the corresponding sCTs using a hybrid deformable image registration (propagated CTVs) in RayStation 11B. For 11 patients, all propagated CTVs were reviewed by our HNC radiation oncologist (physician corrected CTVs). The residual geometrical error dr was quantified as a function of the daily CTVs volume overlap with the nominal plan CTV. The errors dr(propagated CTVs) and dr(physician corrected CTVs) and the difference in dice similarity coefficients (ΔDSC) were determined. Using clinical plans, dose coverage and the tumor control probability (TCP) for the nominal, accumulated and voxel-wise minimum scenarios were determined.

RESULTS

The difference in the residual geometrical error dr (propagated CTVs - physician corrected CTVs) and mean DSC (|ΔDSC|mean) were minor: Δdr(CTV70) = 0.16 mm, Δdr(CTV54.25) = 0.26 mm, |ΔDSC|mean < 0.9%. For all 28 patients, dr(CTV70) = 1.91 mm and dr(CTV54.25) = 1.90 mm. However, CTV54.25 above and below the cricoid cartilage differed substantially (1.00 mm c.f. 3.93 mm). The CTV54.25 coverage below the cricoid was then almost always lower, although the TCP of the accumulated dose was higher than the TCP of the voxel-wise minimum dose.

CONCLUSIONS

Setup uncertainty setting of 2 mm is possible. The feasibility of using propagated CTVs for error determination is demonstrated.

The geometric and dosimetric accuracy of kilovoltage cone beam computed tomography images for adaptive treatment: a systematic review

Objectives

To provide an overview and meta-analysis of different techniques adopted to accomplish kVCBCT for dose calculation and automated segmentation.

Methods

A systematic review and meta-analysis were performed on eligible studies demonstrating kVCBCT-based dose calculation and automated contouring of different tumor features. Meta-analysis of the performance was accomplished on the reported γ analysis and dice similarity coefficient (DSC) score of both collected results as three subgroups (head and neck, chest, and abdomen).

Results

After the literature scrutinization (n = 1008), 52 papers were recognized for the systematic review. Nine studies of dosimtric studies and eleven studies of geometric analysis were suitable for inclusion in meta-analysis. Using kVCBCT for treatment replanning depends on a method used. Deformable Image Registration (DIR) methods yielded small dosimetric error (≤2%), γ pass rate (≥90%) and DSC (≥0.8). Hounsfield Unit (HU) override and calibration curve-based methods also achieved satisfactory yielded small dosimetric error (≤2%) and γ pass rate ((≥90%), but they are prone to error due to their sensitivity to a vendor-specific variation in kVCBCT image quality.

Conclusions

Large cohorts of patients ought to be undertaken to validate methods achieving low levels of dosimetric and geometric errors. Quality guidelines should be established when reporting on kVCBCT, which include agreed metrics for reporting on the quality of corrected kVCBCT and defines protocols of new site-specific standardized imaging used when obtaining kVCBCT images for adaptive radiotherapy.

Advances in knowledge

This review gives useful knowledge about methods making kVCBCT feasible for kVCBCT-based adaptive radiotherapy, simplifying patient pathway and reducing concomitant imaging dose to the patient.

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.

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.

Estimating the accumulative dose uncertainty for intracavitary and interstitial brachytherapy

BACKGROUND

Image-guided adaptive brachytherapy shows the ability to deliver high doses to tumors while sparing normal tissues. However, interfraction dose delivery introduces uncertainties to high dose estimation, which relates to normal tissue toxicity. The purpose of this study was to investigate the high-dose regions of two applicator approaches in brachytherapy.

METHOD

For 32 cervical cancer patients, the CT images from each fraction were wrapped to a reference image, and the displacement vector field (DVF) was calculated with a hybrid intensity-based deformable registration algorithm. The fractional dose was then accumulated to calculate the position and the overlap of high dose (D2cc) during multiple fractions.

RESULT

The overall Dice similarity coefficient (DSC) of the deformation algorithm for the bladder and the rectum was (0.97 and 0.91). No significant difference was observed between the two applicators. However, the location of the intracavitary brachytherapy (ICBT) high-dose region was relatively concentrated. The overlap volume of bladder and rectum D2cc was 0.42 and 0.71, respectively, which was higher than that of interstitial brachytherapy (ISBT) (0.26 and 0.31). The cumulative dose was overestimated in ISBT cases when using the GEC-recommended method. The ratio of bladder and rectum D2cc to the GEC method was 0.99 and 1, respectively, which was higher than that of the ISBT method (0.96 and 0.94).

CONCLUSION

High-dose regions for brachytherapy based on different applicator types were different. The 3D-printed ICBT has better high-dose region consistency than freehand ISBT and hence is more predictable.

Deformable image registration to assist clinical decision for radiotherapy treatment adaptation for head and neck cancer patients

Head and neck (H&N) cancer patients often present anatomical and geometrical changes in tumors and organs at risk (OARs) during radiotherapy treatment. These changes may result in the need to adapt the existing treatment planning, using an expert's subjective opinion, for offline adaptive radiotherapy and a new treatment planning before each treatment, for online adaptive radiotherapy. In the present study, a fast methodology is proposed to assist in planning adaptation clinical decision using tumor and parotid glands percentage volume changes during treatment. The proposed approach was applied to 40 Η&Ν cases, with one planning Computed Tomography (pCT) image and CBCT scans for 6 weeks of treatment per case. Deformable registration was used for each patient's pCT image alignment to its weekly CBCT. The calculated transformations were used to align each patient's anatomical structures to the weekly anatomy. Clinical target volume (CTV) and parotid gland volume percentage changes were calculated in each case. The accuracy of the achieved image alignment was validated qualitatively and quantitatively. Furthermore, statistical analysis was performed to test if there is a statistically significant correlation between CTV and parotid glands volume percentage changes. Average MDA for CTV and parotid glands between corresponding structures defined by an expert in CBCTs and automatically calculated through registration was 1.4 ± 0.1 mm and 1.5 ± 0.1 mm, respectively. The mean registration time of the first CBCT image registration for 40 cases was lower than 3.4 min. Five patients show more than 20% tumor volume change. Six patients show more than 30% parotid glands volume change. Ten out of 40 patients proposed for planning adaptation. All the statistical tests performed showed no correlation between CTV/parotid glands percentage volume changes. The aim to assist in clinical decision making on a fast and automatic way was achieved using the proposed methodology, thereby reducing workload in clinical practice.

Benchmarking of Deformable Image Registration for Multiple Anatomic Sites Using Digital Data Sets With Ground-Truth Deformation Vector Fields

PURPOSE

This study aimed to evaluate the accuracy of deformable image registration (DIR) algorithms using data sets with different levels of ground-truth deformation vector fields (DVFs) and to investigate the correlation between DVF errors and contour-based metrics.

METHODS AND MATERIALS

Nine pairs of digital data sets were generated through contour-controlled deformations based on 3 anonymized patients' CTs (head and neck, thorax/abdomen, and pelvis) with low, medium, and high deformation intensity for each site using the ImSimQA software. Image pairs and their associated contours were imported to MIM-Maestro, Raystation, and Velocity systems, followed by DIR and contour propagation. The system-generated DVF and propagated contours were compared with the ground-truth data. The correlation between DVF errors and contour-based metrics was evaluated using the Pearson correlation coefficient (r), while their correlation with volumes were calculated using Spearman correlation coefficient (rho).

RESULTS

The DVF errors increased with increasing deformation intensity. All DIR algorithms performed well for esophagus, trachea, left femoral, right femoral, and urethral (mean and maximum DVF errors <2.50 mm and <4.27 mm, respectively; Dice similarity coefficient: 0.93-0.99). Brain, liver, left lung, and bladder showed large DVF errors for all 3 systems (d_max: 2.8-91.90 mm). The minimum and maximum DVF errors, conformity index, and Dice similarity coefficient were correlated with volumes (|rho|: 0.41-0.64), especially for very large or small structures (|rho|: 0.64-0.80). Only mean distance to agreement of Raystation and Velocity correlated with some indices of DVF errors (r: 0.70-0.78).

CONCLUSIONS

Most contour-based metrics had no correlation with DVF errors. For adaptive radiation therapy, well-performed contour propagation does not directly indicate accurate dose deformation and summation/accumulation within each contour (determined by DVF accuracy). Tolerance values for DVF errors should vary as the acceptable accuracy for overall adaptive radiation therapy depends on anatomic site, deformation intensity, organ size, and so forth. This study provides benchmark tables for evaluating DIR accuracy in various clinical scenarios.

Use of a Head-Tilting Baseplate During Tomotherapy to Shorten the Irradiation Time and Protect the Hippocampus and Lens in Hippocampal Sparing-Whole Brain Radiotherapy

Purpose:

The aim of this study is to comparatively examine the possibility of reducing the exposure dose to organs at risk, such as the hippocampus and lens, and improving the dose distribution of the planned target volume with and without the use of a head-tilting base plate in hippocampal-sparing whole-brain radiotherapy using tomotherapy.

Methods:

Five paired images of planned head computed tomography without and with tilt were analyzed. The hippocampus and planning target volume were contoured according to the RTOG 0933 contouring atlas protocol. The hippocampal zone to be avoided was delineated using a 5-mm margin. The prescribed radiation dose was 30 Gy in 10 fractions. The absorbed dose to planning target volume dose, absorbed dose to the organ at risk, and irradiation time were evaluated. The paired t-test was used to analyze the differences between hippocampal-sparing whole-brain radiotherapy with head tilts and without head tilts.

Results:

Hippocampal-sparing whole-brain radiotherapy with tilt was not superior in planning target volume doses using the homogeneity index than that without tilt; however, it showed better values, and for D_mean and D_2%, the values were closer to 30 Gy. Regarding the hippocampus, dose reduction with tilt was significantly greater at D_max, D_mean, and D_min, whereas regarding the lens, it was significantly greater at D_max and D_min. The irradiation time was also predominantly shorter.

Conclusion:

In our study, a tilted hippocampal-sparing whole-brain radiotherapy reduced the irradiation time by >10%. Therefore, our study indicated that hippocampal-sparing whole-brain radiotherapy with tomotherapy should be performed with a tilt. The head-tilting technique might be useful during hippocampal-sparing whole-brain radiotherapy. This method could decrease the radiation exposure time, while sparing healthy organs, including the hippocampus and lens.

Head-and-neck organs-at-risk auto-delineation using dual pyramid networks for CBCT-guided adaptive radiotherapy

Organ-at-risk (OAR) delineation is a key step for cone-beam CT (CBCT) based adaptive radiotherapy planning that can be a time-consuming, labor-intensive, and subject-to-variability process. We aim to develop a fully automated approach aided by synthetic MRI for rapid and accurate CBCT multi-organ contouring in head-and-neck (HN) cancer patients. MRI has superb soft-tissue contrasts, while CBCT offers bony-structure contrasts. Using the complementary information provided by MRI and CBCT is expected to enable accurate multi-organ segmentation in HN cancer patients. In our proposed method, MR images are firstly synthesized using a pre-trained cycle-consistent generative adversarial network given CBCT. The features of CBCT and synthetic MRI (sMRI) are then extracted using dual pyramid networks for final delineation of organs. CBCT images and their corresponding manual contours were used as pairs to train and test the proposed model. Quantitative metrics including Dice similarity coefficient (DSC), Hausdorff distance 95% (HD95), mean surface distance, and residual mean square distance (RMS) were used to evaluate the proposed method. The proposed method was evaluated on a cohort of 65 HN cancer patients. CBCT images were collected from those patients who received proton therapy. Overall, DSC values of 0.87 ± 0.03, 0.79 ± 0.10/0.79 ± 0.11, 0.89 ± 0.08/0.89 ± 0.07, 0.90 ± 0.08, 0.75 ± 0.06/0.77 ± 0.06, 0.86 ± 0.13, 0.66 ± 0.14, 0.78 ± 0.05/0.77 ± 0.04, 0.96 ± 0.04, 0.89 ± 0.04/0.89 ± 0.04, 0.83 ± 0.02, and 0.84 ± 0.07 for commonly used OARs for treatment planning including brain stem, left/right cochlea, left/right eye, larynx, left/right lens, mandible, optic chiasm, left/right optic nerve, oral cavity, left/right parotid, pharynx, and spinal cord, respectively, were achieved. This study provides a rapid and accurate OAR auto-delineation approach, which can be used for adaptive radiation therapy.

Comparison of CBCT conversion methods for dose calculation in the head and neck region

The purpose of this study was to compare different methods of CBCT conversion respect to dose calculation accuracy. Twelve head and neck cancer patients treated with VMAT using simultaneous integrated boost technique were selected for the study. For each patient a planning CT (pCT), a control. CT acquired in the fourth week of treatment and a CBCT scan acquired on the closest day with the control CT were used. In order to re-calculate dose directly on CBCT image sets, a population based approach (CBCTPop) and a Histogram Matching (HM) approach based on rigid (CBCTHM-R) and deformable registration (CBCTHM-D) were used. Additionally, virtual CTs (vCTs) were generated using two deformable image registration algorithms (CTELX and CTANC) of the planning CT to the CBCT by using two different deformable image registration (DIR) algorithms. The corresponding control CTs were selected as ground truth and dose distributions on CBCT were analyzed using 3D global gamma index analysis applying a threshold of 10% with respect to the prescribed dose. Using the 2%/2mm gamma criterion, the results were 89.9%(±8.3%), 94.1%(±5.0%), 94.3%(±5.7%), 96.1%(±3.9%), 93.4%(±6.3%) for the CBCTPop, CBCTHM-R, CBCTHM-D, CTELX and CTANC, respectively. On average, the HM and DIR techniques showed a higher accuracy compared to the population based approach, but Kruskal-Wallis test did not show significant difference among the investigated dose calculation techniques assuming p<0.05. More sophisticated CBCT dose calculation methods seem to improve the dose calculation accuracy, but statistical significance remains to be demonstrated.

Noncoplanar Radiation using Tomotherapy: A Phantom Study

Background:

There are very few studies on noncoplanar radiation in tomotherapy because deformable image registration is not implemented in the TomoTherapy Planning Station, a treatment planning device used in tomotherapy. This study examined whether noncoplanar radiation can be performed on the head using a tilt-type head and neck fixture and deformable image registration.

Methods:

Planning target volume spheres with diameters of 2, 3, and 4 cm were set on a head phantom, and computed tomography images were taken at 0° and 40° using a tilt-type head and neck fixture. Irradiation plans were created in the Tomotherapy Planning Station. Noncoplanar radiation was simulated, and the dose volume was evaluated by adding the 0° dose distribution and 40° dose distribution using the deformable image registration of the RayStation treatment planning system.

Results:

The ratio of the phantom volume to the irradiation dose for 20% to 30% of the planning target volume in noncoplanar radiation was smaller than that for 40% to 90% of the planning target volume in single-section irradiation at 0° or 40°.

Conclusions:

Noncoplanar radiation on the head region using tomotherapy was possible by using a tilt-type head and neck fixture, and the dose distribution could be evaluated using deformable image registration. This method helps reduce the dose of the organ-at-risk region located slightly away from the planning target volume.

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.

Performance of a deformable image registration algorithm for CT and cone beam CT using physical multi-density geometric and digital anatomic phantoms

PURPOSE

To study the accuracy of deformable registration algorithm for CT and cone beam CT (CBCT) using a combination of physical and digital phantoms.

MATERIALS AND METHODS

The physical phantoms consisted of objects over a range of electron densities, shape and sizes. The system was tested for simple and complex scenarios including performance in the presence of metallic artefacts. Clinically present deformations were simulated using a set of five geometric and anatomic virtual phantoms.

RESULTS

The system could not account for large changes in size, shape and Hounsfield units. Deformations of low intensity structures and small objects were highly inaccurate, and errors were prominent for volume reduction scenario than volume growth. The presence of artefacts did alter the performance of the algorithm. Objects of low density and that close to artefacts were affected the most. Overall, deformations to CBCT were poor. In virtual phantoms, the system could not handle gas pockets and deformation errors in inverse direction were higher than that in forward direction.

CONCLUSION

The algorithm was tested for several non-clinical and clinical scenarios. The performance was acceptable for realistic and clinically present deformations. However, it is necessary to tread cautiously for structures with small volumes and large reductions in volume.

Proton pencil beam scanning treatment with feedback based voluntary moderate breath hold

Introduction The aim of this article is to introduce a novel protocol for proton pencil beam scanning treatment with moderate deep inspiration breath hold (mDIBH) and report on our clinical implementation results. Methods Three computed tomography (CT) scannings to build the patient's anatomy model were performed during the patient's voluntary mDIBH. All 3 CT scans were used in the optimization during the treatment planning process. Both orthogonal kV imaging and cone-beam computed tomography (CBCT) were implemented for patient alignment with BH prior to the treatment. The BH CBCT images were analyzed for BH reproducibility and the virtual total dose (VTD) retrospectively. To find the VTD, a series of deformable image registrations (DIR) were performed between CBCT and pCT. The effect of the variation of lung density on the dose distribution was also analyzed in the study. Results The values of the mean, standard deviation, maximum, and minimum of the tumor location difference between the CBCT and pCT were 1.9, 1.6, 4.7, and 0.0 mm, respectively. The percentage difference in D99% of CTVs between VTD and the nominal plan was within 1.5%. Conclusions The feedback-based voluntary moderate BH proton PBS treatment was successfully performed in our clinic. This study shows that there is a potential to implement the BH treatment widely in proton centers.

Quantifying the accuracy of deformable image registration for cone-beam computed tomography with a physical phantom

PURPOSE

Kilo-voltage cone-beam computed tomography (CBCT) is widely used for patient alignment, contour propagation, and adaptive treatment planning in radiation therapy. In this study, we evaluated the accuracy of deformable image registration (DIR) for CBCT under various imaging protocols with different noise and patient dose levels.

METHODS

A physical phantom previously developed to facilitate end-to-end testing of the DIR accuracy was used with Varian Velocity v4.0 software to evaluate the performance of image registration from CT to CT, CBCT to CT, and CBCT to CBCT. The phantom is acrylic and includes several inserts that simulate different tissue shapes and properties. Deformations and anatomic changes were simulated by changing the rotations of both the phantom and the inserts. CT images (from a head and neck protocol) and CBCT images (from pelvis, head and "Image Gently" protocols) were obtained with different image noise and dose levels. Large inserts were filled with Mobil DTE oil to simulate soft tissue, and small inserts were filled with bone materials. All inserts were contoured before the DIR process to provide a ground truth contour size and shape for comparison. After the DIR process, all deformed contours were compared with the originals using Dice similarity coefficient (DSC) and mean distance to agreement (MDA). Both large and small volume of interests (VOIs) for DIR volume selection were tested by simulating a DIR process that included whole patient image volume and clinical target volumes (CTV) only (for CTVs propagation).

RESULTS

For cross-modality DIR registration (CT to CBCT), the DSC were >0.8 and the MDA were <3 mm for CBCT pelvis, and CBCT head protocols. For CBCT to CBCT and CT to CT, the DIR accuracy was improved relative to the cross-modality tests. For smaller VOIs, the DSC were >0.8 and MDA <2 mm for all modalities.

CONCLUSIONS

The accuracy of DIR depends on the quality of the CBCT image at different dose and noise levels.

The GIRAFE phase II trial on MVCT-based "volumes of the day" and "dose of the day" addresses when and how to implement adaptive radiotherapy for locally advanced head and neck cancer

During exclusive curative radiotherapy for head and neck tumors, the patient's organs at risk (OAR) and target volumes frequently change size and shape, leading to a risk of higher toxicity and lower control than expected on planned dosimetry. Adaptive radiotherapy is often necessary but 1) tools are needed to define the optimal time for replanning, and 2) the subsequent workflow is time-consuming. We designed a prospective study to evaluate 1) the validity of automatically deformed contours on the daily MVCT, in order to safely use the "dose-of the day" tool to check daily if replanning is necessary; 2) the automatically deformed contours on the replanning CT and the time gained in the replanning workflow. Forty-eight patients with T3-T4 and/or involved node >2 cm head and neck squamous cell carcinomas, planned for curative radiotherapy without surgery, will be enrolled. They will undergo treatment with helical IMRT including daily repositioning MVCTs. The contours proposed will be compared weekly on intermediate planning CTs (iCTs) on weeks 3, 4, 5 and 6. On these iCTs both manual recontouring and automated deformable registration of the initial contours will be compared with the contours automatically defined on the MVCT. The primary objective is to evaluate the Dice similarity coefficient (DSC) of the volumes of each parotid gland. The secondary objectives will evaluate, for target volumes and all OARs: the DSC, the mean distance to agreement, and the average surface-to-surface distance. Time between the automatic and the manual recontouring workflows will be compared.

Motion effects in proton treatments of hepatocellular carcinoma-4D robustly optimised pencil beam scanning plans versus double scattering plans

Pencil beam scanning (PBS) proton therapy enables better dose conformality for complex anatomical geometries than passive proton scattering techniques, but is more susceptible to organ motion. This becomes an issue when treating moving tumours in the thorax or abdomen. Novel four-dimensional treatment planning approaches have been developed to increase the robustness of PBS plans against motion. However, their efficacy still needs to be examined by means of 4D dynamically accumulated dose (4DDD) analyses. This study investigates the potential use of 4D robust optimisation to maintain sufficient target coverage in the presence of organ motion, while sparing surrounding healthy tissue, for hepatocellular carcinoma (HCC). The liver is particularly suited to study motion interplay effects since the treatment region exhibits smaller density gradients and more homogeneous tissue than targets in the thorax, making it less prone to range errors. A facility-specific beam time model, developed and experimentally validated previously, was used for the clinical evaluation. 4DDD analyses of eleven target volumes did not show a significant improvement of the target coverage using 4D robust optimisation, but a reduction of the dose to close-by organs at risk. Interplay effects were averaged out for the applied fractionation scheme of 15 fractions. Contrary to PBS, passive double scattering (DS) plans yielded homogeneous 4DDD dose distributions in a single fraction. But, in some cases, they exceeded organ at risk dose limits, which were only satisfied in PBS. The average normal liver dose could be decreased by almost 6% compared to non-robustly optimised PBS plans and by 16% compared to DS plans when implementing 4D robust optimisation. Except for some very small tumours with large motion amplitudes, 4D robustly optimised PBS plans were found to be clinically acceptable even without supplementary motion mitigation techniques.

Registration methods in radiotherapy

Purpose

The aim of this study is to present a short and comprehensive review of the methods of medical image registration, their conditions and applications in radiotherapy. A particular focus was placed on the methods of deformable image registration.

Methods

To structure and deepen the knowledge on medical image registration in radiotherapy, a medical literature analysis was made using the Google Scholar browser and the medical database of the PubMed library.

Results

Chronological review of image registration methods in radiotherapy based on 34 selected articles. A particular attention was given to show: (i) potential regions of the application of different methods of registration, (ii) mathematical basis of the deformable methods and (iii) the methods of quality control for the registration process.

Conclusions

The primary aim of the medical image registration process is to connect the contents of images. What we want to achieve is a complementary or extended knowledge that can be used for more precise localisation of pathogenic lesions and continuous improvement of patient treatment. Therefore, the choice of imaging mode is dependent on the type of clinical study. It is impossible to visualise all anatomical details or functional changes using a single modality machine. Therefore, fusion of various modality images is of great clinical relevance. A natural problem in analysing the fusion of medical images is geographical errors related to displacement. The registered images are performed not at the same time and, very often, at different respiratory phases.

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.

Real-time energy/mass transfer mapping for online 4D dose reconstruction

In this work we describe an ultra-fast, low-latency implementation of the energy/mass transfer (EMT) mapping method to accumulate dose on deforming geometries such as lung using the central processing unit (CPU). It enables the computation of the actually delivered dose for intensity-modulated radiation therapy on 4D image data in real-time at 25 Hz. In order to accumulate the delivered dose onto a reference phase a pre-calculated deformable vector field is used. The aim of this study is to present an online dose accumulation technique that can be carried out in less than 40 ms to accommodate the machine log update rate of our research linac. Three speed optimisation strategies for the CPU are discussed: single-core optimisation, parallelisation for multiple cores and vectorisation. The single-core implementation accumulates dose in about 1.1 s on a typical high resolution grid for a lung stereotactic body radiation therapy case. Adding parallelisation decreased the runtime to about 50 ms while adding vectorisation satisfied our real-time constraint by further reducing the dose accumulation time to 15 ms without compromising on resolution or accuracy. The presented method allows real-time dose accumulation on deforming patient geometries and has the potential to enable online dose evaluation and re-planning scenarios.

Evaluation of rectum and bladder dose accumulation from external beam radiotherapy and brachytherapy for cervical cancer using two different deformable image registration techniques

We evaluated dose-volume histogram (DVH) parameters based on deformable image registration (DIR) between brachytherapy (BT) and external beam radiotherapy (EBRT) that included a center-shielded (CS) plan. Eleven cervical cancer patients were treated with BT, and their pelvic and CS EBRT were studied. Planning CT images for EBRT and BT (except for the first BT, used as the reference image) were deformed with DIR to reference image. We used two DIR parameter settings: intensity-based and hybrid. Mean Dice similarity coefficients (DSCs) comparing EBRT with the reference for the uterus, rectum and bladder were 0.81, 0.77 and 0.83, respectively, for hybrid DIR and 0.47, 0.37 and 0.42, respectively, for intensity-based DIR (P < 0.05). D1 cm3 for hybrid DIR, intensity-based DIR and DVH addition were 75.1, 81.2 and 78.2 Gy, respectively, for the rectum, whereas they were 93.5, 92.3 and 94.3 Gy, respectively, for the bladder. D2 cm3 for hybrid DIR, intensity-based DIR and DVH addition were 70.1, 74.0 and 71.4 Gy, respectively, for the rectum, whereas they were 85.4, 82.8 and 85.4 Gy, respectively, for the bladder. Overall, hybrid DIR obtained higher DSCs than intensity-based DIR, and there were moderate differences in DVH parameters between the two DIR methods, although the results varied among patients. DIR is only experimental, and extra care should be taken when comparing DIR-based dose values with dose-effect curves established using DVH addition. Also, a true evaluation of DIR-based dose accumulation would require ground truth data (e.g. measurement with physical phantom).

Evaluation of the performance of deformable image registration between planning CT and CBCT images for the pelvic region: comparison between hybrid and intensity-based DIR

This study aimed to evaluate the performance of the hybrid deformable image registration (DIR) method in comparison with intensity-based DIR for pelvic cone-beam computed tomography (CBCT) images, using intensity and anatomical information. Ten prostate cancer patients treated with intensity-modulated radiation therapy (IMRT) were studied. Nine or ten CBCT scans were performed for each patient. First, rigid registration was performed between the planning CT and all CBCT images using gold fiducial markers, and then DIR was performed. The Dice similarity coefficient (DSC) and center of mass (COM) displacement were used to evaluate the quantitative DIR accuracy. The average DSCs for intensity-based DIR for the prostate, rectum, bladder, and seminal vesicles were 0.84 ± 0.05, 0.75 ± 0.05, 0.69 ± 0.07 and 0.65 ± 0.11, respectively, whereas those values for hybrid DIR were 0.98 ± 0.00, 0.97 ± 0.01, 0.98 ± 0.00 and 0.94 ± 0.03, respectively (P < 0.05). The average COM displacements for intensity-based DIR for the prostate, rectum, bladder, and seminal vesicles were 2.0 ± 1.5, 3.7 ± 1.4, 7.8 ± 2.2 and 3.6 ± 1.2 mm, whereas those values for hybrid DIR were 0.1 ± 0.0, 0.3 ± 0.2, 0.2 ± 0.1 and 0.6 ± 0.6 mm, respectively (P < 0.05). These results showed that the DSC for hybrid DIR had a higher DSC value and smaller COM displacement for all structures and all patients, compared with intensity-based DIR. Thus, the accumulative dose based on hybrid DIR might be trusted as a high-precision dose estimation method that takes into account organ movement during treatment radiotherapy.

Dose tracking assessment for image-guided radiotherapy of the prostate bed and the impact on clinical workflow

Background

The cumulative dose was compared with the planned dose among fourteen patients undergoing image-guided, intensity-modulated radiotherapy of the prostate bed. Moreover, we investigated the feasibility of adding dose tracking to the routine workflow for radiotherapy.

Methods

Daily cone beam computed tomography was conducted for image-guided radiotherapy, and weekly cumulative delivered doses were calculated for dose tracking. Deformable image registration was applied to map weekly dose distributions to the original treatment plan and to create a cumulative dose distribution. The dose–volume histogram (DVH) cut-off points for the rectum and bladder and the planning target volume (PTV), were used to compare the planned and cumulative delivered doses. The additional time required by the departmental staff to complete these duties was recorded.

Results

The PTV coverage of the delivered treatment did not satisfy the expected goal for three patients (V98% >98%). In another three patients, the DVH cut-off point for the bladder was higher than the limits, while for the rectum, treatment was as expected in all cases (two patients failed both their bladder constraints and the PTV coverage). Overall, four patients did not satisfy one or more criteria at the end of their treatment.

Conclusions

A well-defined strategy for dose tracking assessment is feasible, would have minimal impact on the workload of a radiotherapy department, and may offer objective information to support radiation oncologists in making decisions about adaptive procedures.

Evaluation of deformable image registration between external beam radiotherapy and HDR brachytherapy for cervical cancer with a 3D-printed deformable pelvis phantom

PURPOSE

In this study, we developed a 3D-printed deformable pelvis phantom for evaluating spatial DIR accuracy. We then evaluated the spatial DIR accuracies of various DIR settings for cervical cancer.

METHODS

A deformable female pelvis phantom was created based on patient CT data using 3D printing. To create the deformable uterus phantom, we first 3D printed both a model of uterus and a model of the internal cavities of the vagina and uterus. We then made a mold using the 3D printed uterus phantom. Finally, urethane was poured into the mold with the model of the internal cavities in place, creating the deformable uterus phantom with a cavity into which an applicator could be inserted. To create the deformable bladder phantom, we first 3D printed models of the bladder and of the same bladder scaled down by 2 mm. We then made a mold using the larger bladder model. Finally, silicone was poured into the mold with the smaller bladder model in place to create the deformable bladder phantom with a wall thickness of 2 mm. To emulate the anatomical bladder, water was poured into the created bladder. We acquired phantom image without applicator for EBRT. Then, we inserted the applicator into the phantom to simulate BT. In this situation, we scanned the phantom again to obtain the phantom image for BT. We performed DIR using the two phantom images in two cases: Case A, with full bladder (170 ml) in both EBRT and BT images; and Case B with full bladder in the BT image and half-full bladder (100 ml) in the EBRT image. DIR was evaluated using Dice similarity coefficients (DSCs) and 31 landmarks for the uterus and 25 landmarks for the bladder. A hybrid intensity and structure DIR algorithm implemented in RayStation with four DIR settings was evaluated.

RESULTS

On visual inspection, reasonable agreement in shape of the uterus between the phantom and patient CT images was observed for both EBRT and BT, although some regional disagreements in shape of the bladder and rectum were apparent. The created phantom could reproduce the actual patient's uterus deformation by the applicator. For both Case A and B, large variation was seen in landmark error among the four DIR parameters. In addition, although DSCs were comparable, moderate differences in landmark error existed between the two different DIR parameters selected from the four DIR parameters (i.e., DSC = 0.96, landmark error = 13.2 ± 5.7 mm vs. DSC = 0.97, landmark error = 9.7 ± 4.0 mm). This result suggests that landmark error evaluation might thus be more effective than DSC for evaluating DIR accuracy.

CONCLUSIONS

Our developed phantom enabled the evaluation of spatial DIR accuracy for the female pelvic region for the first time. Although the DSCs are high, the spatial errors can still be significant and our developed phantom facilitates their quantification. Our results showed that optimization is needed to identify suitable DIR settings. For determining suitable DIR settings, our method of evaluating spatial DIR accuracy using the 3D-printed phantom may prove helpful.

Real-time 4D dose reconstruction for tracked dynamic MLC deliveries for lung SBRT

PURPOSE

This study provides a proof of concept for real-time 4D dose reconstruction for lung stereotactic body radiation therapy (SBRT) with multileaf collimator (MLC) tracking and assesses the impact of tumor tracking on the size of target margins.

METHODS

The authors have implemented real-time 4D dose reconstruction by connecting their tracking and delivery software to an Agility MLC at an Elekta Synergy linac and to their in-house treatment planning software (TPS). Actual MLC apertures and (simulated) target positions are reported to the TPS every 40 ms. The dose is calculated in real-time from 4DCT data directly after each reported aperture by utilization of precalculated dose-influence data based on a Monte Carlo algorithm. The dose is accumulated onto the peak-exhale (reference) phase using energy-mass transfer mapping. To investigate the impact of a potentially reducible safety margin, the authors have created and delivered treatment plans designed for a conventional internal target volume (ITV) + 5 mm, a midventilation approach, and three tracking scenarios for four lung SBRT patients. For the tracking plans, a moving target volume (MTV) was established by delineating the gross target volume (GTV) on every 4DCT phase. These were rigidly aligned to the reference phase, resulting in a unified maximum GTV to which a 1, 3, or 5 mm isotropic margin was added. All scenarios were planned for 9-beam step-and-shoot IMRT to meet the criteria of RTOG 1021 (3 × 18 Gy). The GTV 3D center-of-volume shift varied from 6 to 14 mm.

RESULTS

Real-time dose reconstruction at 25 Hz could be realized on a single workstation due to the highly efficient implementation of dose calculation and dose accumulation. Decreased PTV margins resulted in inadequate target coverage during untracked deliveries for patients with substantial tumor motion. MLC tracking could ensure the GTV target dose for these patients. Organ-at-risk (OAR) doses were consistently reduced by decreased PTV margins. The tracked MTV + 1 mm deliveries resulted in the following OAR dose reductions: lung V20 up to 3.5%, spinal cord D2 up to 0.9 Gy/Fx, and proximal airways D2 up to 1.4 Gy/Fx.

CONCLUSIONS

The authors could show that for patient data at clinical resolution and realistic motion conditions, the delivered dose could be reconstructed in 4D for the whole lung volume in real-time. The dose distributions show that reduced margins yield lower doses to healthy tissue, whilst target dose can be maintained using dynamic MLC tracking.

A support vector machine tool for adaptive tomotherapy treatments: Prediction of head and neck patients criticalities

PURPOSE

Adaptive radiation therapy (ART) is an advanced field of radiation oncology. Image-guided radiation therapy (IGRT) methods can support daily setup and assess anatomical variations during therapy, which could prevent incorrect dose distribution and unexpected toxicities. A re-planning to correct these anatomical variations should be done daily/weekly, but to be applicable to a large number of patients, still require time consumption and resources. Using unsupervised machine learning on retrospective data, we have developed a predictive network, to identify patients that would benefit of a re-planning.

METHODS

1200 MVCT of 40 head and neck (H&N) cases were re-contoured, automatically, using deformable hybrid registration and structures mapping. Deformable algorithm and MATLAB(®) homemade machine learning process, developed, allow prediction of criticalities for Tomotherapy treatments.

RESULTS

Using retrospective analysis of H&N treatments, we have investigated and predicted tumor shrinkage and organ at risk (OAR) deformations. Support vector machine (SVM) and cluster analysis have identified cases or treatment sessions with potential criticalities, based on dose and volume discrepancies between fractions. During 1st weeks of treatment, 84% of patients shown an output comparable to average standard radiation treatment behavior. Starting from the 4th week, significant morpho-dosimetric changes affect 77% of patients, suggesting need for re-planning. The comparison of treatment delivered and ART simulation was carried out with receiver operating characteristic (ROC) curves, showing monotonous increase of ROC area.

CONCLUSIONS

Warping methods, supported by daily image analysis and predictive tools, can improve personalization and monitoring of each treatment, thereby minimizing anatomic and dosimetric divergences from initial constraints.