Scale invariant feature transform in adaptive radiation therapy: a tool for deformable image registration assessment and re-planning indication

Structure-specific rigid dose accumulation dosimetric analysis of ablative stereotactic MRI-guided adaptive radiation therapy in ultracentral lung lesions

BACKGROUND

Definitive local therapy with stereotactic ablative radiation therapy (SABR) for ultracentral lung lesions is associated with a high risk of toxicity, including treatment related death. Stereotactic MR-guided adaptive radiation therapy (SMART) can overcome many of the challenges associated with SABR treatment of ultracentral lesions.

METHODS

We retrospectively identified 14 consecutive patients who received SMART to ultracentral lung lesions from 10/2019 to 01/2021. Patients had a median distance from the proximal bronchial tree (PBT) of 0.38 cm. Tumors were most often lung primary (64.3%) and HILUS group A (85.7%). A structure-specific rigid registration approach was used for cumulative dose analysis. Kaplan-Meier log-rank analysis was used for clinical outcome data and the Wilcoxon Signed Rank test was used for dosimetric data.

RESULTS

Here we show that SMART dosimetric improvements in favor of delivered plans over predicted non-adapted plans for PBT, with improvements in proximal bronchial tree DMax of 5.7 Gy (p = 0.002) and gross tumor 100% prescription coverage of 7.3% (p = 0.002). The mean estimated follow-up is 17.2 months and 2-year local control and local failure free survival rates are 92.9% and 85.7%, respectively. There are no grade ≥ 3 toxicities.

CONCLUSIONS

SMART has dosimetric advantages and excellent clinical outcomes for ultracentral lung tumors. Daily plan adaptation reliably improves target coverage while simultaneously reducing doses to the proximal airways. These results further characterize the therapeutic window improvements for SMART. Structure-specific rigid dose accumulation dosimetric analysis provides insights that elucidate the dosimetric advantages of SMART more so than per fractional analysis alone.

Detecting and quantifying spatial misalignment between longitudinal kilovoltage computed tomography (kVCT) scans of the head and neck by using convolutional neural networks (CNNs)

BACKGROUND

Adaptive radiotherapy (ART) aims to address anatomical modifications appearing during the treatment of patients by modifying the planning treatment according to the daily positioning image. Clinical implementation of ART relies on the quality of the deformable image registration (DIR) algorithms included in the ART workflow. To translate ART into clinical practice, automatic DIR assessment is needed.

OBJECTIVE

This article aims to estimate spatial misalignment between two head and neck kilovoltage computed tomography (kVCT) images by using two convolutional neural networks (CNNs).

METHODS

The first CNN quantifies misalignments between 0 mm and 15 mm and the second CNN detects and classifies misalignments into two classes (poor alignment and good alignment). Both networks take pairs of patches of 33x33x33 mm3 as inputs and use only the image intensity information. The training dataset was built by deforming kVCT images with basis splines (B-splines) to simulate DIR error maps. The test dataset was built using 2500 landmarks, consisting of hard and soft landmark tissues annotated by 6 clinicians at 10 locations.

RESULTS

The quantification CNN reaches a mean error of 1.26 mm (± 1.75 mm) on the landmark set which, depending on the location, has annotation errors between 1 mm and 2 mm. The errors obtained for the quantification network fit the computed interoperator error. The classification network achieves an overall accuracy of 79.32%, and although the classification network overdetects poor alignments, it performs well (i.e., it achieves a rate of 90.4%) in detecting poor alignments when given one.

CONCLUSION

The performances of the networks indicate the feasibility of using CNNs for an agnostic and generic approach to misalignment quantification and detection.

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.

Adaptive radiation therapy: When, how and what are the benefits that literature provides?

PURPOSE

To identify from the current literature when is the right time to replan and to assign thresholds for the optimum process of replanning. Nowadays, adaptive radiotherapy (ART) for head and neck cancer plays an exceptional role consisting of an evaluation procedure of the prominent anatomical and dosimetric variations. By performing complex radiotherapy methods, the credibility of the therapeutic result is crucial. Image guided radiotherapy (IGRT) was developed to ensure locoregional control and thus changes that might occur during radiotherapy be dealt with.

MATERIALS AND METHODS

An electronic research of articles published in PubMed/MEDLINE and Science Direct databases from January 2004 to October 2020 was performed. Among a total of 127 studies assessed for eligibility, 85 articles were ultimately retained for the review.

RESULTS

The most noticeable changes have been reported in the middle fraction of the treatment. Therefore, the suggested optimal time to replan is between the third and the fourth week. Anatomical deviations>1cm in the external contour, average weight loss>10%, violation in the dose coverage of the targets>5%, and violation in the dose of the peripherals were some of the thresholds that are currently used, and which lead to replanning.

CONCLUSION

ART may decrease toxicity and improve local-control. Whether it is beneficial or not, depends ultimately on each patient. However, more investigation of the changes should be performed in future prospective studies to obtain more accurate results.

Detection of vessel bifurcations in CT scans for automatic objective assessment of deformable image registration accuracy

PURPOSE

Objective assessment of deformable image registration (DIR) accuracy often relies on the identification of anatomical landmarks in image pairs, a manual process known to be extremely time-expensive. The goal of this study is to propose a method to automatically detect vessel bifurcations in images and assess their use for the computation of target registration errors (TREs).

MATERIALS AND METHODS

Three image datasets were retrospectively analyzed. The first dataset included 10 pairs of inhale/exhale phases from lung 4DCTs and full inhale and exhale breath-hold CT scans from 10 patients presenting with chronic obstructive pulmonary disease, with 300 corresponding landmarks available for each case (DIR-Lab). The second dataset included six pairs of inhale/exhale phases from lung 4DCTs (POPI dataset), with 100 pairs of landmarks for each case. The third dataset included 28 pairs of pre/post-radiotherapy liver contrast-enhanced CT scans, each with five manually picked vessel bifurcation correspondences. For all images, the vasculature was autosegmented by computing and thresholding a vesselness image. Images of the vasculature centerline were computed, and bifurcations were detected based on centerline voxel neighbors' count. The vasculature segmentations were independently registered using a Demons algorithm between representations of their surface with distance maps. Detected bifurcations were considered as corresponding when distant by less than 5 mm after vasculature DIR. The selected pairs of bifurcations were used to calculate TRE after registration of the images considering three algorithms: rigid registration, Anaconda, and a Demons algorithm. For comparison with the ground truth, TRE values calculated using the automatically detected correspondences were interpolated in the whole organs to generate TRE maps. The performance of the method in automatically calculating TRE after image registration was quantified by measuring the correlation with the TRE obtained when using the ground truth landmarks.

RESULTS

The median Pearson correlation coefficients between ground truth TRE and corresponding values in the generated TRE maps were r = 0.81 and r = 0.67 for the lung and liver cases, respectively. The correlation coefficients between mean TRE for each case were r = 0.99 and r = 0.64 for the lung and liver cases, respectively.

CONCLUSION

For lungs or liver CT scans DIR, a strong correlation was obtained between TRE calculated using manually picked or landmarks automatically detected with the proposed method. This tool should be particularly useful in studies requiring assessing the reliability of a high number of DIRs.

An approach for estimating dosimetric uncertainties in deformable dose accumulation in pencil beam scanning proton therapy for lung cancer

Deformable image registration (DIR) is an important component for dose accumulation and associated clinical outcome evaluation in radiotherapy. However, the resulting deformation vector field (DVF) is subject to unavoidable discrepancies when different algorithms are applied, leading to dosimetric uncertainties of the accumulated dose. We propose here an approach for proton therapy to estimate dosimetric uncertainties as a consequence of modeled or estimated DVF uncertainties. A patient-specific DVF uncertainty model was built on the first treatment fraction, by correlating the magnitude differences of five DIR results at each voxel to the magnitude of any single reference DIR. In the following fractions, only the reference DIR needs to be applied, and DVF geometric uncertainties were estimated by this model. The associated dosimetric uncertainties were then derived by considering the estimated geometric DVF uncertainty, the dose gradient of fractional recalculated dose distribution and the direction factor from the applied reference DIR of this fraction. This estimated dose uncertainty was respectively compared to the reference dose uncertainty when different DIRs were applied individually for each dose warping. This approach was validated on seven NSCLC patients, each with nine repeated CTs. The proposed model-based method is able to achieve dose uncertainty distribution on a conservative voxel-to-voxel comparison within ±5% of the prescribed dose to the 'reference' dosimetric uncertainty, for 77% of the voxels in the body and 66%-98% of voxels in investigated structures. We propose a method to estimate DIR induced uncertainties in dose accumulation for proton therapy of lung tumor treatments.

An MRI framework for respiratory motion modelling validation

Introduction

Respiratory motion models establish a correspondence between respiratory‐correlated (RC) 4‐dimensional (4D) imaging and respiratory surrogates, to estimate time‐resolved (TR) 3D breathing motion. To evaluate the performance of motion models on real patient data, a validation framework based on magnetic resonance imaging (MRI) is proposed, entailing the use of RC 4DMRI to build the model, and on both (i) TR 2D cine‐MRI and (ii) additional 4DMRI data for testing intra‐/inter‐fraction breathing motion variability.

Methods

Repeated MRI data were acquired in 7 patients with abdominal lesions. The considered model relied on deformable image registration (DIR) for building the model and compensating for inter‐fraction baseline variations. Both 2D and 3D validation were performed, by comparing model estimations with the ground truth 2D cine‐MRI and 4DMRI respiratory phases, respectively.

Results

The median DIR error was comparable to the voxel size (1.33 × 1.33 × 5 mm³), with higher values in the presence of large inter‐fraction motion (median value: 2.97 mm). In the 2D validation, the median estimation error on anatomical landmarks’ position resulted below 4 mm in every scenario, whereas in the 3D validation it was 1.33 mm and 4.21 mm when testing intra‐ and inter‐fraction motion, respectively. The range of motion described in the cine‐MRI was comparable to the motion of the building 4DMRI, being always above the estimation error. Overall, the model performance was dependent on DIR error, presenting reduced accuracy when inter‐fraction baseline variations occurred.

Conclusions

Results suggest the potential of the proposed framework in evaluating global motion models for organ motion management in MRI‐guided radiotherapy.

RegQCNET: Deep quality control for image-to-template brain MRI affine registration

Affine registration of one or several brain image(s) onto a common reference space is a necessary prerequisite for many image processing tasks, such as brain segmentation or functional analysis. Manual assessment of registration quality is a tedious and time-consuming task, especially in studies comprising a large amount of data. Automated and reliable quality control (QC) becomes mandatory. Moreover, the computation time of the QC must be also compatible with the processing of massive datasets. Therefore, automated deep neural network approaches have emerged as a method of choice to automatically assess registration quality. In the current study, a compact 3D convolutional neural network, referred to as RegQCNET, is introduced to quantitatively predict the amplitude of an affine registration mismatch between a registered image and a reference template. This quantitative estimation of registration error is expressed using the metric unit system. Therefore, a meaningful task-specific threshold can be manually or automatically defined in order to distinguish between usable and non-usable images. The robustness of the proposed RegQCNET is first analyzed on lifespan brain images undergoing various simulated spatial transformations and intensity variations between training and testing. Secondly, the potential of RegQCNET to classify images as usable or non-usable is evaluated using both manual and automatic thresholds. During our experiments, automatic thresholds are estimated using several computer-assisted classification models (logistic regression, support vector machine, Naive Bayes and random forest) through cross-validation. To this end we use an expert's visual QC estimated on a lifespan cohort of 3953 brains. Finally, the RegQCNET accuracy is compared to usual image features such as image correlation coefficient and mutual information. The results show that the proposed deep learning QC is robust, fast and accurate at estimating affine registration error in the processing pipeline.

Investigating the use of virtual 4DCT from 4DMRI in gated carbon ion radiation therapy of abdominal tumors

PURPOSE

To generate virtual 4DCT from 4DMRI with field of view (FOV) extended to the entire involved patient anatomy, in order to evaluate its use in carbon ion radiation therapy (CIRT) of the abdominal site in a clinical scenario.

MATERIALS AND METHODS

The virtual 4DCT was generated by deforming a reference CT in order to (1) match the anatomy depicted in the 4DMRI within its FOV, by calculating deformation fields with deformable image registration to describe inter-fractional and breathing motion, and (2) obtain physically plausible deformation outside of the 4DMRI FOV, by propagating and modulating the previously obtained deformation fields. The implemented method was validated on a digital anthropomorphic phantom, for which a ground truth (GT) 4DCT was available. A CIRT treatment plan was optimized at the end-exhale reference CT and the RBE-weighted dose distribution was recalculated on both the virtual and GT 4DCTs. The method estimation error was quantified by comparing the virtual and GT 4DCTs and the corresponding recomputed doses. The method was then evaluated on 8 patients with pancreas or liver tumors treated with CIRT using respiratory gating at end-exhale. The clinical treatment plans adopted at the National Center for Oncological Hadrontherapy (CNAO, Pavia, Italy) were considered and the dose distribution was recomputed on all respiratory phases of the planning and virtual 4DCTs. By comparing the two datasets and the corresponding dose distributions, the geometrical and dosimetric impact of organ motion was assessed.

RESULTS

For the phantom, the error outside of the 4DMRI FOV was up to 4.5mm, but it remained sub-millimetric in correspondence to the target within the 4DMRI FOV. Although the impact of motion on the target D_95% resulted in variations ranging from 22% to 90% between the planned dose and the doses recomputed on the GT 4DCT phases, the corresponding estimation error was ≤2.2%. In the patient cases, the variation of the baseline tumor position between the planning and the virtual end-exhale CTs presented a median (interquartile range) value of 6.0 (4.9) mm. For baseline variations larger than 5mm, the tumor D_95% variation between the plan and the dose recomputed on the end-exhale virtual CT resulted larger than 10%. Median variations higher than 10% in the target D_95% and gastro-intestinal OARs D_2% were quantified at the end-inhale, whereas close to the end-exhale phase, limited variations of relevant dose metrics were found for both tumor and OARs.

CONCLUSIONS

The negligible impact of the geometrical inaccuracy in the estimated anatomy outside of the 4DMRI FOV on the overall dosimetric accuracy suggests the feasibility of virtual 4DCT with extended FOV in CIRT of the abdominal site. In the analyzed patient group, inter-fractional variations such as baseline variation and breathing variability were quantified, demonstrating the method capability to support treatment planning in gated CIRT of the abdominal site.

Virtual 4DCT from 4DMRI for the management of respiratory motion in carbon ion therapy of abdominal tumors

PURPOSE

To evaluate a method for generating virtual four-dimensional computed tomography (4DCT) from four-dimensional magnetic resonance imaging (4DMRI) data in carbon ion radiotherapy with pencil beam scanning for abdominal tumors.

METHODS

Deformable image registration is used to: (a) register each respiratory phase of the 4DMRI to the end-exhale MRI; (b) register the reference end-exhale CT to the end-exhale MRI volume; (c) generate the virtual 4DCT by warping the registered CT according to the obtained deformation fields. A respiratory-gated carbon ion treatment plan is optimized on the planning 4DCT and the corresponding dose distribution is recalculated on the virtual 4DCT. The method was validated on a digital anthropomorphic phantom and tested on eight patients (18 acquisitions). For the phantom, a ground truth dataset was available to assess the method performances from the geometrical and dosimetric standpoints. For the patients, the virtual 4DCT was compared with the planning 4DCT.

RESULTS

In the phantom, the method exhibits a geometrical accuracy within the voxel size and Dose Volume Histograms deviations up to 3.3% for target V_95% (mean dose difference ≤ 0.2% of the prescription dose, gamma pass rate > 98%). For patients, the virtual and the planning 4DCTs show good agreement at end-exhale (3% median D_95% difference), whereas other respiratory phases exhibit moderate motion variability with consequent dose discrepancies, confirming the need for motion mitigation strategies during treatment.

CONCLUSIONS

The virtual 4DCT approach is feasible to evaluate treatment plan robustness against intra- and interfraction motion in carbon ion therapy delivered at the abdominal site.

Automatic large quantity landmark pairs detection in 4DCT lung images

PURPOSE

To automatically and precisely detect a large quantity of landmark pairs between two lung computed tomography (CT) images to support evaluation of deformable image registration (DIR). We expect that the generated landmark pairs will significantly augment the current lung CT benchmark datasets in both quantity and positional accuracy.

METHODS

A large number of landmark pairs were detected within the lung between the end-exhalation (EE) and end-inhalation (EI) phases of the lung four-dimensional computed tomography (4DCT) datasets. Thousands of landmarks were detected by applying the Harris-Stephens corner detection algorithm on the probability maps of the lung vasculature tree. A parametric image registration method (pTVreg) was used to establish initial landmark correspondence by registering the images at EE and EI phases. A multi-stream pseudo-siamese (MSPS) network was then developed to further improve the landmark pair positional accuracy by directly predicting three-dimensional (3D) shifts to optimally align the landmarks in EE to their counterparts in EI. Positional accuracies of the detected landmark pairs were evaluated using both digital phantoms and publicly available landmark pairs.

RESULTS

Dense sets of landmark pairs were detected for 10 4DCT lung datasets, with an average of 1886 landmark pairs per case. The mean and standard deviation of target registration error (TRE) were 0.47 ± 0.45 mm with 98% of landmark pairs having a TRE smaller than 2 mm for 10 digital phantom cases. Tests using 300 manually labeled landmark pairs in 10 lung 4DCT benchmark datasets (DIRLAB) produced TRE results of 0.73 ± 0.53 mm with 97% of landmark pairs having a TRE smaller than 2 mm.

CONCLUSION

A new method was developed to automatically and precisely detect a large quantity of landmark pairs between lung CT image pairs. The detected landmark pairs could be used as benchmark datasets for more accurate and informative quantitative evaluation of DIR algorithms.

Accuracy of low-dose proton CT image registration for pretreatment alignment verification in reference to planning proton CT

Purpose

Proton CT (pCT) has the ability to reduce inherent uncertainties in proton treatment by directly measuring the relative proton stopping power with respect to water, thereby avoiding the uncertain conversion of X‐ray CT Hounsfield unit to relative stopping power and the deleterious effect of X‐ ray CT artifacts. The purpose of this work was to further evaluate the potential of pCT for pretreatment positioning using experimental pCT data of a head phantom.

Methods

The performance of a 3D image registration algorithm was tested with pCT reconstructions of a pediatric head phantom. A planning pCT simulation scan of the phantom was obtained with 200 MeV protons and reconstructed with a 3D filtered back projection (FBP) algorithm followed by iterative reconstruction and a representative pretreatment pCT scan was reconstructed with FBP only to save reconstruction time. The pretreatment pCT scan was rigidly transformed by prescribing random errors with six degrees of freedom or deformed by the deformation field derived from a head and neck cancer patient to the pretreatment pCT reconstruction, respectively. After applying the rigid or deformable image registration algorithm to retrieve the original pCT image before transformation, the accuracy of the registration was assessed. To simulate very low‐dose imaging for patient setup, the proton CT images were reconstructed with 100%, 50%, 25%, and 12.5% of the total number of histories of the original planning pCT simulation scan, respectively.

Results

The residual errors in image registration were lower than 1 mm and 1° of magnitude regardless of the anatomic directions and imaging dose. The mean residual errors ranges found for rigid image registration were from −0.29 ± 0.09 to 0.51 ± 0.50 mm for translations and from −0.05 ± 0.13 to 0.08 ± 0.08 degrees for rotations. The percentages of sub‐millimetric errors found, for deformable image registration, were between 63.5% and 100%.

Conclusion

This experimental head phantom study demonstrated the potential of low‐dose pCT imaging for 3D image registration. Further work is needed to confirm the value pCT for pretreatment image‐guided proton therapy.

A clustering approach to 4D MRI retrospective sorting for the investigation of different surrogates

PURPOSE

In retrospective 4-Dimensional Magnetic Resonance Imaging (4D MRI) sorting, respiratory surrogate selection affects the image quality of reconstructed volumes. We propose a method for retrospective 4D MRI sorting based on clustering, which allowed us to compare the performance of single or multiple internal surrogates vs. a conventional external signal.

METHODS

A k-medoids clustering algorithm was exploited for sorting 2D MRI into 4D MRI, relying on (A) multiple or (B) single automatically tracked internal landmarks or (C) respiratory belt signal. 4D MRI reconstructions for seven liver cancer patients were compared to those of the state-of-the-art mutual information (MI) approach. Sorting artifacts were measured by the root mean square error (RMSE) between the diaphragm profile and a fitted second order curve. Diaphragm and tumor motions were evaluated.

RESULTS

The median RMSEs ranged 0.97-1.66 mm, 1.24-1.89 mm, 1.43-2.27 mm, 1.74-3.72 mm for the MI, (A), (B) and (C) methods, respectively. Significant differences (Friedman, α = 5%) were found between (C) and all other methods, and between (B) and MI approaches. The discrepancies between (A) and MI approaches ranged 1.1-6.2 mm and 0.7-5.3 mm respectively in diaphragm and tumor motions. Methods (A) and (B) showed similar ranges of motion.

CONCLUSION

With multiple internal points, our method yielded the description of a higher range of motion and similar image quality with respect to the MI approach. The single point method led to more artifacts, suggesting the superior suitability of multiple internal surrogates for retrospective 4D MRI sorting. Considering internal rather than external information favored superior performance.

The VAMPIRE challenge: A multi-institutional validation study of CT ventilation imaging

PURPOSE

CT ventilation imaging (CTVI) is being used to achieve functional avoidance lung cancer radiation therapy in three clinical trials (NCT02528942, NCT02308709, NCT02843568). To address the need for common CTVI validation tools, we have built the Ventilation And Medical Pulmonary Image Registration Evaluation (VAMPIRE) Dataset, and present the results of the first VAMPIRE Challenge to compare relative ventilation distributions between different CTVI algorithms and other established ventilation imaging modalities.

METHODS

The VAMPIRE Dataset includes 50 pairs of 4DCT scans and corresponding clinical or experimental ventilation scans, referred to as reference ventilation images (RefVIs). The dataset includes 25 humans imaged with Galligas 4DPET/CT, 21 humans imaged with DTPA-SPECT, and 4 sheep imaged with Xenon-CT. For the VAMPIRE Challenge, 16 subjects were allocated to a training group (with RefVI provided) and 34 subjects were allocated to a validation group (with RefVI blinded). Seven research groups downloaded the Challenge dataset and uploaded CTVIs based on deformable image registration (DIR) between the 4DCT inhale/exhale phases. Participants used DIR methods broadly classified into B-splines, Free-form, Diffeomorphisms, or Biomechanical modeling, with CT ventilation metrics based on the DIR evaluation of volume change, Hounsfield Unit change, or various hybrid approaches. All CTVIs were evaluated against the corresponding RefVI using the voxel-wise Spearman coefficient r S , and Dice similarity coefficients evaluated for low function lung ( DSC low ) and high function lung ( DSC high ).

RESULTS

A total of 37 unique combinations of DIR method and CT ventilation metric were either submitted by participants directly or derived from participant-submitted DIR motion fields using the in-house software, VESPIR. The r S and DSC results reveal a high degree of inter-algorithm and intersubject variability among the validation subjects, with algorithm rankings changing by up to ten positions depending on the choice of evaluation metric. The algorithm with the highest overall cross-modality correlations used a biomechanical model-based DIR with a hybrid ventilation metric, achieving a median (range) of 0.49 (0.27-0.73) for r S , 0.52 (0.36-0.67) for DSC low , and 0.45 (0.28-0.62) for DSC high . All other algorithms exhibited at least one negative r S value, and/or one DSC value less than 0.5.

CONCLUSIONS

The VAMPIRE Challenge results demonstrate that the cross-modality correlation between CTVIs and the RefVIs varies not only with the choice of CTVI algorithm but also with the choice of RefVI modality, imaging subject, and the evaluation metric used to compare relative ventilation distributions. This variability may arise from the fact that each of the different CTVI algorithms and RefVI modalities provides a distinct physiologic measurement. Ultimately this variability, coupled with the lack of a "gold standard," highlights the ongoing importance of further validation studies before CTVI can be widely translated from academic centers to the clinic. It is hoped that the information gleaned from the VAMPIRE Challenge can help inform future validation efforts.

Performance of a feature-based algorithm for 3D-3D registration of CT angiography to cone-beam CT for endovascular repair of complex abdominal aortic aneurysms

BACKGROUND

A crucial step in image fusion for intraoperative guidance during endovascular procedures is the registration of preoperative computed tomography angiography (CTA) with intraoperative Cone Beam CT (CBCT). Automatic tools for image registration facilitate the 3D image guidance workflow. However their performance is not always satisfactory. The aim of this study is to assess the accuracy of a new fully automatic, feature-based algorithm for 3D3D registration of CTA to CBCT.

METHODS

The feature-based algorithm was tested on clinical image datasets from 14 patients undergoing complex endovascular aortic repair. Deviations in Euclidian distances between vascular as well as bony landmarks were measured and compared to an intensity-based, normalized mutual information algorithm.

RESULTS

The results for the feature-based algorithm showed that the median 3D registration error between the anatomical landmarks of CBCT and CT images was less than 3 mm. The feature-based algorithm showed significantly better accuracy compared to the intensity-based algorithm (p < 0.001).

CONCLUSION

A feature-based algorithm for 3D image registration is presented.

Online daily assessment of dose change in head and neck radiotherapy without dose-recalculation

BACKGROUND

Head and neck cancers are commonly treated with radiation therapy, but due to possible volume changes, plan adaptation may be required during the course of treatment. Currently, plan adaptations consume significant clinical resources. Existing methods to evaluate the need for plan adaptation requires deformable image registration (DIR) to a new CT simulation or daily cone beam CT (CBCT) images and the recalculation of the dose distribution. In this study, we explore a tool to assist the decision for plan adaptation using a CBCT without re-computation of dose, allowing for rapid online assessment.

METHODS

This study involved 18 head and neck cancer patients treated with CBCT image guidance who had their treatment plan modified based on a new CT simulation (ReCT). Dose changes were estimated using different methods and compared to the current gold standard of using DIR between the planning CT scan (PCT) and ReCT with recomputed dose. The first and second methods used DIR between the PCT and daily CBCT with the planned dose or recalculated dose from the ReCT respectively, with the dose transferred to the CBCT using rigid registration. The necessity of plan adaptation was assessed by the change in dose to 95% of the planning target volume (D95) and mean dose to the parotids.

RESULTS

The treatment plans were adapted clinically for all 18 patients but only 7 actually needed an adaptation yielding 11 unnecessary adaptations. Applying a method using the daily CBCT with the planned dose distribution would have yielded only four unnecessary adaptations and no missed adaptations: a significant improvement from that done clinically.

CONCLUSION

Using the DIR between the planning CT and daily CBCT can flag cases for plan adaptation before every fraction while not requiring a new re-planning CT scan and dose recalculation.

A phantom study to create synthetic CT from orthogonal twodimensional cine MRI and evaluate the effect of irregular breathing

An exciting innovation in radiotherapy is the use of real-time MRI for treatment adaptation. This study proposes an in-silico framework for the generation of 3D synthetic CT (sCT) from orthogonal interleaved 2D cine MRI data to overcome the lack of electron density information in MR images. The method uses pre-treatment data to build a patient breathing motion model. This model is then driven by surrogates extracted from cine MR images during the treatment. The effect of irregular breathing on the motion model is also evaluated by simulating different motion components related to uncorrelated diaphragm, chest and tumor motion. 3D sCT were successfully created for each of the 512 cine MRI pairs in the digital phantom study. The analysis showed that the diaphragm position was a good surrogate to rescale the 3D breathing motion for the current regular breathing phase. However, respiratory and tumor motion were correlated in only 59% of the phases, resulting in tumor position uncertainties of up to 3mm. The inclusion of additional chest and tumor motion information improved the accuracy for irregular changes in breathing pattern and enhanced the tumor position uncertainties to less than 1mm. This study successfully demonstrated a proof-ofprinciple for a digital phantom dataset based on patient parameters, providing a way to create real-time 3D electron density volumes and enhancing the need to account for irregular breathing pattern.

Using gamma index to flag changes in anatomy during image-guided radiation therapy of head and neck cancer

During radiation therapy of head and neck cancer, the decision to consider replanning a treatment because of anatomical changes has significant resource implications. We developed an algorithm that compares cone‐beam computed tomography (CBCT) image pairs and provides an automatic alert as to when remedial action may be required. Retrospective CBCT data from ten head and neck cancer patients that were replanned during their treatment was used to train the algorithm on when to recommend a repeat CT simulation (re‐CT). An additional 20 patients (replanned and not replanned) were used to validate the predictive power of the algorithm. CBCT images were compared in 3D using the gamma index, combining Hounsfield Unit (HU) difference with distance‐to‐agreement (DTA), where the CBCT study acquired on the first fraction is used as the reference. We defined the match quality parameter (MQP_x) as a difference between the x^th percentiles of the failed‐pixel histograms calculated from the reference gamma comparison and subsequent comparisons, where the reference gamma comparison is taken from the first two CBCT images acquired during treatment. The decision to consider re‐CT was based on three consecutive MQP values being less than or equal to a threshold value, such that re‐CT recommendations were within ±3 fractions of the actual re‐CT order date for the training cases. Receiver‐operator characteristic analysis showed that the best trade‐off in sensitivity and specificity was achieved using gamma criteria of 3 mm DTA and 30 HU difference, and the 80^th percentile of the failed‐pixel histogram. A sensitivity of 82% and 100% was achieved in the training and validation cases, respectively, with a false positive rate of ~30%. We have demonstrated that gamma analysis of CBCT‐acquired anatomy can be used to flag patients for possible replanning in a manner consistent with local clinical practice guidelines.

A method to detect landmark pairs accurately between intra-patient volumetric medical images

PURPOSES

An image processing procedure was developed in this study to detect large quantity of landmark pairs accurately in pairs of volumetric medical images. The detected landmark pairs can be used to evaluate of deformable image registration (DIR) methods quantitatively.

METHODS

Landmark detection and pair matching were implemented in a Gaussian pyramid multi-resolution scheme. A 3D scale-invariant feature transform (SIFT) feature detection method and a 3D Harris-Laplacian corner detection method were employed to detect feature points, i.e., landmarks. A novel feature matching algorithm, Multi-Resolution Inverse-Consistent Guided Matching or MRICGM, was developed to allow accurate feature pairs matching. MRICGM performs feature matching using guidance by the feature pairs detected at the lower resolution stage and the higher confidence feature pairs already detected at the same resolution stage, while enforces inverse consistency.

RESULTS

The proposed feature detection and feature pair matching algorithms were optimized to process 3D CT and MRI images. They were successfully applied between the inter-phase abdomen 4DCT images of three patients, between the original and the re-scanned radiation therapy simulation CT images of two head-neck patients, and between inter-fractional treatment MRIs of two patients. The proposed procedure was able to successfully detect and match over 6300 feature pairs on average. The automatically detected landmark pairs were manually verified and the mismatched pairs were rejected. The automatic feature matching accuracy before manual error rejection was 99.4%. Performance of MRICGM was also evaluated using seven digital phantom datasets with known ground truth of tissue deformation. On average, 11855 feature pairs were detected per digital phantom dataset with TRE = 0.77 ± 0.72 mm.

CONCLUSION

A procedure was developed in this study to detect large number of landmark pairs accurately between two volumetric medical images. It allows a semi-automatic way to generate the ground truth landmark datasets that allow quantitatively evaluation of DIR algorithms for radiation therapy applications.

Evaluating the accuracy of 4D-CT ventilation imaging: First comparison with Technegas SPECT ventilation

PURPOSE

Computed tomography ventilation imaging (CTVI) is a highly accessible functional lung imaging modality that can unlock the potential for functional avoidance in lung cancer radiation therapy. Previous attempts to validate CTVI against clinical ventilation single-photon emission computed tomography (V-SPECT) have been hindered by radioaerosol clumping artifacts. This work builds on those studies by performing the first comparison of CTVI with ^99m Tc-carbon ('Technegas'), a clinical V-SPECT modality featuring smaller radioaerosol particles with less clumping.

METHODS

Eleven lung cancer radiotherapy patients with early stage (T1/T2N0) disease received treatment planning four-dimensional CT (4DCT) scans paired with Technegas V/Q-SPECT/CT. For each patient, we applied three different CTVI methods. Two of these used deformable image registration (DIR) to quantify breathing-induced lung density changes (CTVI_DIR-HU ), or breathing-induced lung volume changes (CTVI_DIR-Jac ) between the 4DCT exhale/inhale phases. A third method calculated the regional product of air-tissue densities (CTVI_HU ) and did not involve DIR. Corresponding CTVI and V-SPECT scans were compared using the Dice similarity coefficient (DSC) for functional defect and nondefect regions, as well as the Spearman's correlation r computed over the whole lung. The DIR target registration error (TRE) was quantified using both manual and computer-selected anatomic landmarks.

RESULTS

Interestingly, the overall best performing method (CTVI_HU ) did not involve DIR. For nondefect regions, the CTVI_HU , CTVI_DIR-HU , and CTVI_DIR-Jac methods achieved mean DSC values of 0.69, 0.68, and 0.54, respectively. For defect regions, the respective DSC values were moderate: 0.39, 0.33, and 0.44. The Spearman r-values were generally weak: 0.26 for CTVI_HU , 0.18 for CTVI_DIR-HU , and -0.02 for CTVI_DIR-Jac . The spatial accuracy of CTVI was not significantly correlated with TRE, however the DIR accuracy itself was poor with TRE > 3.6 mm on average, potentially indicative of poor quality 4DCT. Q-SPECT scans achieved good correlations with V-SPECT (mean r > 0.6), suggesting that the image quality of Technegas V-SPECT was not a limiting factor in this study.

CONCLUSIONS

We performed a validation of CTVI using clinically available 4DCT and Technegas V/Q-SPECT for 11 lung cancer patients. The results reinforce earlier findings that the spatial accuracy of CTVI exhibits significant interpatient and intermethod variability. We propose that the most likely factor affecting CTVI accuracy was poor image quality of clinical 4DCT.

Examination of a deformable motion model for respiratory movements and 4D dose calculations using different driving surrogates

PURPOSE

The aim of this study was to evaluate a surrogate-driven motion model based on four-dimensional computed tomography that is able to predict CT volumes corresponding to arbitrary respiratory phases. Furthermore, the comparison of three different driving surrogates is examined and the feasibility of using the model for 4D dose re-calculation will be discussed.

METHODS

The study is based on repeated 4DCTs of twenty patients treated for bronchial carcinoma and metastasis. The motion model was estimated from the planning 4DCT through deformable image registration. To predict a certain phase of a follow-up 4DCT, the model considers inter-fractional variations (baseline correction) and intra-fractional respiratory parameters (amplitude and phase) derived from surrogates. The estimated volumes resulting from the model were compared to ground-truth clinical 4DCTs using absolute HU differences in the lung region and landmarks localized using the Scale Invariant Feature Transform. Finally, the γ-index was used to evaluate the dosimetric effects of the intensity differences measured between the estimated and the ground-truth CT volumes.

RESULTS

The results show absolute HU differences between estimated and ground-truth images with median value (± standard deviation) of (61.3 ± 16.7) HU. Median 3D distances, measured on about 400 matching landmarks in each volume, were (2.9 ± 3.0) mm. 3D errors up to 28.2 mm were found for CT images with artifacts or reduced quality. Pass rates for all surrogate approaches were above 98.9% with a γ-criterion of 2%/2 mm.

CONCLUSION

The results depend mainly on the image quality of the initial 4DCT and the deformable image registration. All investigated surrogates can be used to estimate follow-up 4DCT phases, however, uncertainties decrease for volumetric approaches. Application of the model for 4D dose calculations is feasible.

Deformable registration of CT and cone-beam CT with local intensity matching

Cone-beam CT (CBCT) is a widely used intra-operative imaging modality in image-guided radiotherapy and surgery. A short scan followed by a filtered-backprojection is typically used for CBCT reconstruction. While data on the mid-plane (plane of source-detector rotation) is complete, off-mid-planes undergo different information deficiency and the computed reconstructions are approximate. This causes different reconstruction artifacts at off-mid-planes depending on slice locations, and therefore impedes accurate registration between CT and CBCT. In this paper, we propose a method to accurately register CT and CBCT by iteratively matching local CT and CBCT intensities. We correct CBCT intensities by matching local intensity histograms slice by slice in conjunction with intensity-based deformable registration. The correction-registration steps are repeated in an alternating way until the result image converges. We integrate the intensity matching into three different deformable registration methods, B-spline, demons, and optical flow that are widely used for CT-CBCT registration. All three registration methods were implemented on a graphics processing unit for efficient parallel computation. We tested the proposed methods on twenty five head and neck cancer cases and compared the performance with state-of-the-art registration methods. Normalized cross correlation (NCC), structural similarity index (SSIM), and target registration error (TRE) were computed to evaluate the registration performance. Our method produced overall NCC of 0.96, SSIM of 0.94, and TRE of 2.26 → 2.27 mm, outperforming existing methods by 9%, 12%, and 27%, respectively. Experimental results also show that our method performs consistently and is more accurate than existing algorithms, and also computationally efficient.

Investigating CT to CBCT image registration for head and neck proton therapy as a tool for daily dose recalculation

PURPOSE

Intensity modulated proton therapy (IMPT) of head and neck (H&N) cancer patients may be improved by plan adaptation. The decision to adapt the treatment plan based on a dose recalculation on the current anatomy requires a diagnostic quality computed tomography (CT) scan of the patient. As gantry-mounted cone beam CT (CBCT) scanners are currently being offered by vendors, they may offer daily or weekly updates of patient anatomy. CBCT image quality may not be sufficient for accurate proton dose calculation and it is likely necessary to perform CBCT CT number correction. In this work, the authors investigated deformable image registration (DIR) of the planning CT (pCT) to the CBCT to generate a virtual CT (vCT) to be used for proton dose recalculation.

METHODS

Datasets of six H&N cancer patients undergoing photon intensity modulated radiation therapy were used in this study to validate the vCT approach. Each dataset contained a CBCT acquired within 3 days of a replanning CT (rpCT), in addition to a pCT. The pCT and rpCT were delineated by a physician. A Morphons algorithm was employed in this work to perform DIR of the pCT to CBCT following a rigid registration of the two images. The contours from the pCT were deformed using the vector field resulting from DIR to yield a contoured vCT. The DIR accuracy was evaluated with a scale invariant feature transform (SIFT) algorithm comparing automatically identified matching features between vCT and CBCT. The rpCT was used as reference for evaluation of the vCT. The vCT and rpCT CT numbers were converted to stopping power ratio and the water equivalent thickness (WET) was calculated. IMPT dose distributions from treatment plans optimized on the pCT were recalculated with a Monte Carlo algorithm on the rpCT and vCT for comparison in terms of gamma index, dose volume histogram (DVH) statistics as well as proton range. The DIR generated contours on the vCT were compared to physician-drawn contours on the rpCT.

RESULTS

The DIR accuracy was better than 1.4 mm according to the SIFT evaluation. The mean WET differences between vCT (pCT) and rpCT were below 1 mm (2.6 mm). The amount of voxels passing 3%/3 mm gamma criteria were above 95% for the vCT vs rpCT. When using the rpCT contour set to derive DVH statistics from dose distributions calculated on the rpCT and vCT the differences, expressed in terms of 30 fractions of 2 Gy, were within [-4, 2 Gy] for parotid glands (D(mean)), spinal cord (D(2%)), brainstem (D(2%)), and CTV (D(95%)). When using DIR generated contours for the vCT, those differences ranged within [-8, 11 Gy].

CONCLUSIONS

In this work, the authors generated CBCT based stopping power distributions using DIR of the pCT to a CBCT scan. DIR accuracy was below 1.4 mm as evaluated by the SIFT algorithm. Dose distributions calculated on the vCT agreed well to those calculated on the rpCT when using gamma index evaluation as well as DVH statistics based on the same contours. The use of DIR generated contours introduced variability in DVH statistics.

Possibility Study of Scale Invariant Feature Transform (SIFT) Algorithm Application to Spine Magnetic Resonance Imaging

The purpose of this study is an application of scale invariant feature transform (SIFT) algorithm to stitch the cervical-thoracic-lumbar (C-T-L) spine magnetic resonance (MR) images to provide a view of the entire spine in a single image. All MR images were acquired with fast spin echo (FSE) pulse sequence using two MR scanners (1.5 T and 3.0 T). The stitching procedures for each part of spine MR image were performed and implemented on a graphic user interface (GUI) configuration. Moreover, the stitching process is performed in two categories; manual point-to-point (mPTP) selection that performed by user specified corresponding matching points, and automated point-to-point (aPTP) selection that performed by SIFT algorithm. The stitched images using SIFT algorithm showed fine registered results and quantitatively acquired values also indicated little errors compared with commercially mounted stitching algorithm in MRI systems. Our study presented a preliminary validation of the SIFT algorithm application to MRI spine images, and the results indicated that the proposed approach can be performed well for the improvement of diagnosis. We believe that our approach can be helpful for the clinical application and extension of other medical imaging modalities for image stitching.

Challenges and opportunities in image guided particle therapy

The application of biomedical imaging and image processing to radiation therapy with accelerated particles has unique challenges. The potential of particle therapy to precisely tailor the dose distribution around the target volume needs to account for the intrinsic sensitivity to uncertainties in dose deposition. These peculiar features motivate the use of image guided methods to consistently verify the accuracy in dose delivery. Dedicated imaging and image processing methods are required, from treatment planning to treatment verification phases, in order to reduce the effects of uncertainties. The scenario is also complicated by the lack of standardized layouts of treatment bunkers, which implies the relatively increased use of custom solutions. Conversely, imaging can be applied to verify the actual delivered dose, representing a valuable opportunity to validate specific protocols and visualize the efficacy of the intended treatment. In this contribution, challenges and opportunities in image guided particle therapy are overviewed, with a clear focus on research perspectives in biomedical imaging and image processing.

Comparing cone-beam CT intensity correction methods for dose recalculation in adaptive intensity-modulated photon and proton therapy for head and neck cancer

BACKGROUND

Adaptive intensity-modulated photon and proton radiotherapy (IMRT and IMPT) of head and neck (H&N) cancer requires frequent three-dimensional (3D) dose calculation. We compared two approaches for dose recalculation on the basis of intensity-corrected cone-beam (CB) x-ray computed tomography (CT) images.

MATERIAL AND METHODS

For nine H&N tumor patients, virtual CTs (vCT) were generated by deformable image registration of the planning CT (pCT) to the CBCT. The second intensity correction approach used population-based lookup tables for scaling CBCT intensities to the pCT HU range (CBCTLUT). IMRT and IMPT plans were generated with a commercial treatment planning system. Dose recalculations on vCT and CBCTLUT were analyzed using a (3%, 3 mm) gamma-index analysis and comparison of normal tissue and tumor dose/volume parameters. A replanning CT (rpCT) acquired within three days of the CBCT served as reference. Single field uniform dose (SFUD) proton plans were created and recalculated on vCT and CBCTLUT for proton range comparison.

RESULTS

Dose/volume parameters showed minor differences between rpCT, vCT and CBCTLUT in IMRT, but clinically relevant deviations between CBCTLUT and rpCT in the spinal cord for IMPT. Gamma-index pass-rates were found increased for vCT with respect to CBCTLUT in IMPT (by up to 21 percentage points) and IMRT (by up to 9 percentage points) for most cases. The SFUD-based proton range assessment showed improved agreement of vCT and rpCT, with 88-99% of the depth dose profiles in beam's eye view agreeing within 3 mm. For CBCTLUT, only 80-94% of the profiles fulfilled this criterion.

CONCLUSION

vCT and CBCTLUT are suitable options for dose recalculation in adaptive IMRT. In the scope of IMPT, the vCT approach is preferable.

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

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

Phantom based evaluation of CT to CBCT image registration for proton therapy dose recalculation

The ability to perform dose recalculation on the anatomy of the day is important in the context of adaptive proton therapy. The objective of this study was to investigate the use of deformable image registration (DIR) and cone beam CT (CBCT) imaging to generate the daily stopping power distribution of the patient. We investigated the deformation of the planning CT scan (pCT) onto daily CBCT images to generate a virtual CT (vCT) using a deformable phantom designed for the head and neck (H & N) region. The phantom was imaged at a planning CT scanner in planning configuration, yielding a pCT and in deformed, treatment day configuration, yielding a reference CT (refCT). The treatment day configuration was additionally scanned at a CBCT scanner. A Morphons DIR algorithm was used to generate a vCT. The accuracy of the vCT was evaluated by comparison to the refCT in terms of corresponding features as identified by an adaptive scale invariant feature transform (aSIFT) algorithm. Additionally, the vCT CT numbers were compared to those of the refCT using both profiles and regions of interest and the volumes and overlap (DICE coefficients) of various phantom structures were compared. The water equivalent thickness (WET) of the vCT, refCT and pCT were also compared to evaluate proton range differences. Proton dose distributions from the same initial fluence were calculated on the refCT, vCT and pCT and compared in terms of proton range. The method was tested on a clinical dataset using a replanning CT scan acquired close in time to a CBCT scan as reference using the WET evaluation. Results from the aSIFT investigation suggest a deformation accuracy of 2-3 mm. The use of the Morphon algorithm did not distort CT number intensity in uniform regions and WET differences between vCT and refCT were of the order of 2% of the proton range. This result was confirmed by proton dose calculations. The patient results were consistent with phantom observations. In conclusion, our phantom study suggests the vCT approach is adequate for proton dose recalculation on the basis of CBCT imaging.

Advances in 4D treatment planning for scanned particle beam therapy - report of dedicated workshops

We report on recent progress in the field of mobile tumor treatment with scanned particle beams, as discussed in the latest editions of the 4D treatment planning workshop. The workshop series started in 2009, with about 20 people from 4 research institutes involved, all actively working on particle therapy delivery and development. The first workshop resulted in a summary of recommendations for the treatment of mobile targets, along with a list of requirements to apply these guidelines clinically. The increased interest in the treatment of mobile tumors led to a continuously growing number of attendees: the 2012 edition counted more than 60 participants from 20 institutions and commercial vendors. The focus of research discussions among workshop participants progressively moved from 4D treatment planning to complete 4D treatments, aiming at effective and safe treatment delivery. Current research perspectives on 4D treatments include all critical aspects of time resolved delivery, such as in-room imaging, motion detection, beam application, and quality assurance techniques. This was motivated by the start of first clinical treatments of hepato cellular tumors with a scanned particle beam, relying on gating or abdominal compression for motion mitigation. Up to date research activities emphasize significant efforts in investigating advanced motion mitigation techniques, with a specific interest in the development of dedicated tools for experimental validation. Potential improvements will be made possible in the near future through 4D optimized treatment plans that require upgrades of the currently established therapy control systems for time resolved delivery. But since also these novel optimization techniques rely on the validity of the 4DCT, research focusing on alternative 4D imaging technique, such as MRI based 4DCT generation will continue.

Uncertainties encountered in implementation of adaptive planning with in vivo dosimeters

Under a previously approved institutional review board protocol for prostate cancer patients, implanted metal-oxide semiconductor field-effect transistor dosimeters (dose verification system, Sicel Technologies) were used for measurement of the in vivo delivered daily dose. This dosimetric information provided the ability to adapt the plan if the measured doses did not match the dose expected from the planning system. Due to the inherent uncertainty in the dosimeters, the decision to adapt the treatment plan was made only if readings differed by more than 7 % for three consecutive days. To validate this method, we acquired daily cone beam computed tomography images for five patients, and the dose delivered to the dosimeters was calculated by use of (1) an automated procedure (MIM Maestro, MIM Software) and (2) the treatment planning system (XIO, Elekta). 72 % of the doses calculated automatically fell within 1 % of the doses calculated in the planning system, and 99 % agreed within 2 %. When compared to the calculated dose, 53 % of the in vivo measurements fell within 3 % of the calculated dose, 80 % fell within 5 %, and 9.8 % were greater than 7 %, but never on three consecutive days. The measured doses agreed reasonably well with the calculated doses, supporting the decision to adapt the plan only if there were discrepancies of more than 7 % over three consecutive days. Even with the inherent uncertainty in the dosimeters, this adaptive planning method can detect delivery inaccuracies that would not otherwise be caught with the use of only daily image guidance or other dose calculation surrogates.

The frequency of re-planning and its variability dependent on the modification of the re-planning criteria and IGRT correction strategy in head and neck IMRT

Background

To analyse the frequency of re-planning and its variability dependent on the IGRT correction strategy and on the modification of the dosimetric criteria for re-planning for the spinal cord in head and neck IG-IMRT.

Methods

Daily kV-control-CTs of six head and neck patients (=175 CTs) were analysed. All volumes of interest were re-contoured using deformable image registration. Three IGRT correction strategies were simulated and the resulting dose distributions were computed for all fractions. Different sets of criteria with varying dose thresholds for re-planning were investigated. All sets of criteria ensure equivalent target coverage of both CTVs, but vary in the tolerance threshold of the spinal cord.

Results

The variations of the D95 and D2 in respect to the planned values ranged from -7% to +3% for both CTVs, and -2% to +6% for the spinal cord. Despite different correction vectors of the three IGRT strategies, the dosimetric differences were small. The number of fractions not requiring re-planning varied between 0% and 11% dependent on the applied IGRT correction strategy. In contrast, this number ranged between 32% and 70% dependent on the dosimetric thresholds, even though these thresholds were only gently modified.

Conclusions

The more precise the planned dose needs to be maintained over the treatment course, the more frequently re-planning is required. The influence of different IGRT correction strategies, even though geometrically notable, was found to be of only limited relevance for the re-planning frequency. In contrast, the definition and modification of thresholds for re-planning have a major impact on the re-planning frequency.

Efficient iris recognition based on optimal subfeature selection and weighted subregion fusion

In this paper, we propose three discriminative feature selection strategies and weighted subregion matching method to improve the performance of iris recognition system. Firstly, we introduce the process of feature extraction and representation based on scale invariant feature transformation (SIFT) in detail. Secondly, three strategies are described, which are orientation probability distribution function (OPDF) based strategy to delete some redundant feature keypoints, magnitude probability distribution function (MPDF) based strategy to reduce dimensionality of feature element, and compounded strategy combined OPDF and MPDF to further select optimal subfeature. Thirdly, to make matching more effective, this paper proposes a novel matching method based on weighted sub-region matching fusion. Particle swarm optimization is utilized to accelerate achieve different sub-region's weights and then weighted different subregions' matching scores to generate the final decision. The experimental results, on three public and renowned iris databases (CASIA-V3 Interval, Lamp, andMMU-V1), demonstrate that our proposed methods outperform some of the existing methods in terms of correct recognition rate, equal error rate, and computation complexity.

An optimised IGRT correction vector determined from a displacement vector field: a proof of principle of a decision-making aid for re-planning

BACKGROUND

To present a new method that determines an optimised IGRT couch correction vector from a displacement vector field (DVF). The DVF is computed by a deformable image registration (DIR) method. The proposed method can improve the quality of volume-of-interest (VOI) alignment in image guided radiation therapy (IGRT), and can serve as a decision-making aid for re-planning.

MATERIAL AND METHODS

The proposed method was demonstrated using the CT data sets of 11 head-and-neck cancer patients with daily kilovoltage control-CTs. A DVF was computed for each control-CT using a DIR method. The DVF was used for voxel tracking and re-contouring of the VOIs in the control-CTs. Then a rigid body transformation, which could be used as couch correction vector, was optimised. The aim of the optimisation process was to find a vector and rotations that map the deformed VOIs into a specified territory. This territory was defined by a margin extension of the VOIs at the time of the planning process. Within this extension, VOI motion and deformation was tolerated. The objective function in the optimisation process was the sum of all volume fractions outside the defined territories.

RESULTS

The proposed method was able to find a correction vector, which resulted in a coverage of the target volumes of at least 98% in 52.3% of all fractions. In contrast, a standard IGRT correction using a rigid registration method only fulfilled this criterion in 22.6% of all fractions. The optimisation process took an average of 1.5 minutes per fraction.

CONCLUSION

The knowledge of the deformation of the anatomy allows the determination of an optimised rigid correction vector using our method. The method ensures controlled mapping of the VOIs despite small deformations. If no optimised vector can be determined, re-planning should be considered. Thus, our method can also serve as a decision-making aid for re-planning.

The adaptive FEM elastic model for medical image registration

This paper proposes an adaptive mesh refinement strategy for the finite element method (FEM) based elastic registration model. The signature matrix for mesh refinement takes into account the regional intensity variance and the local deformation displacement. The regional intensity variance reflects detailed information for improving registration accuracy and the deformation displacement fine-tunes the mesh refinement for a more efficient algorithm. The gradient flows of two different similarity metrics, the sum of the squared difference and the spatially encoded mutual information for the mono-modal and multi-modal registrations, are used to derive external forces to drive the model to the equilibrium state. We compared our approach to three other models: (1) the conventional multi-resolution FEM registration algorithm; (2) the FEM elastic method that uses variation information for mesh refinement; and (3) the robust block matching based registration. Comparisons among different methods in a dataset with 20 CT image pairs upon artificial deformation demonstrate that our registration method achieved significant improvement in accuracies. Experimental results in another dataset of 40 real medical image pairs for both mono-modal and multi-modal registrations also show that our model outperforms the other three models in its accuracy.

Validation of Automatic Contour Propagation for 4D Treatment Planning Using Multiple Metrics

The aim of this work is to provide insights into multiple metrics clinical validation of deformable image registration and contour propagation methods in 4D lung radiotherapy planning. The following indices were analyzed and compared: Volume Difference (VD), Dice Similarity Coefficient (DSC), Positive Predictive Value (PPV) and Surface Distances (SD). The analysis was performed on three patient datasets, using as reference a ground-truth volume generated by means of Simultaneous Truth And Performance Level Estimation (STAPLE) algorithm from the outlines of five experts. Significant discrepancies in the quality assessment provided by the different metrics in all the examined cases were found. Metrics sensitivity was more evident in presence of image artifacts and particularly for tubular anatomical structures, such as esophagus or spinal cord. Volume Differences did not account for position and DSC exhibited criticalities due to its intrinsic symmetry ( i.e. over- and under-estimation of the reference contours cannot be discriminated) and dependency on the total volume of the structure. PPV analysis showed more robust performance, as each voxel concurs to the classification of the propagation, but was not able to detect inclusion of propagated and ground-truth volumes. Mesh distances could interpret the actual shape of the structures, but might report higher mismatches in case of large local differences in the contour surfaces. According to our study, the combination of VD and SD for the validation of contour propagation algorithms in 4D could provide the necessary failure detection accuracy.

A quantitative comparison of the performance of three deformable registration algorithms in radiotherapy

We present an evaluation of various non-rigid registration algorithms for the purpose of compensating interfractional motion of the target volume and organs at risk areas when acquiring CBCT image data prior to irradiation. Three different deformable registration (DR) methods were used: the Demons algorithm implemented in the iPlan Software (BrainLAB AG, Feldkirchen, Germany) and two custom-developed piecewise methods using either a Normalized Correlation or a Mutual Information metric (featurelet_NC and featurelet_MI). These methods were tested on data acquired using a novel purpose-built phantom for deformable registration and clinical CT/CBCT data of prostate and lung cancer patients. The Dice similarity coefficient (DSC) between manually drawn contours and the contours generated by a derived deformation field of the structures in question was compared to the result obtained with rigid registration (RR). For the phantom, the piecewise methods were slightly superior, the featurelet_NC for the intramodality and the featurelet_MI for the intermodality registrations. For the prostate cases in less than 50% of the images studied the DSC was improved over RR. Deformable registration methods improved the outcome over a rigid registration for lung cases and in the phantom study, but not in a significant way for the prostate study. A significantly superior deformation method could not be identified.