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

Ultra-hypofractionated prostate cancer radiotherapy: Dosimetric impact of real-time intrafraction prostate motion and daily anatomical changes

PURPOSE

To retrospectively assess the differences between planned and delivered dose during ultra-hypofractionated (UHF) prostate cancer treatments, by evaluating the dosimetric impact of daily anatomical variations alone, and in combination with prostate intrafraction motion.

METHODS

Prostate intrafraction motion was recorded with a transperineal ultrasound probe in 15 patients treated by UHF radiotherapy (36.25 Gy/5 fractions). The dosimetric objective was to cover 99 % of the clinical target volume with the 100 % prescription isodose line. After treatment, planning CT (pCT) images were deformably registered onto daily Cone Beam CT to generate pseudo-CT for dose accumulation (accumulated CT, aCT). The interplay effect was accounted by synchronizing prostatic shifts and beam geometry. Finally, the shifted dose maps were accumulated (moved-accumulated CT, maCT).

RESULTS

No significant change in daily CTV volumes was observed. Conversely, CTV V100% was 98.2 ± 0.8 % and 94.7 ± 2.6 % on aCT and maCT, respectively, compared with 99.5 ± 0.2 % on pCT (p < 0.0001). Bladder volume was smaller than planned in 76 % of fractions and D5cc was 33.8 ± 3.2 Gy and 34.4 ± 3.4 Gy on aCT (p = 0.02) and maCT (p = 0.01) compared with the pCT (36.0 ± 1.1 Gy). The rectum was smaller than planned in 50.3 % of fractions, but the dosimetric differences were not statistically significant, except for D1cc, found smaller on the maCT (33.2 ± 3.2 Gy, p = 0.02) compared with the pCT (35.3 ± 0.7 Gy).

CONCLUSIONS

Anatomical variations and prostate movements had more important dosimetric impact than anatomical variations alone, although, in some cases, the two phenomena compensated. Therefore, an efficient IGRT protocol is required for treatment implementation to reduce setup errors and control intrafraction motion.

Camouflaged Nanoreactors Mediated Radiotherapy-Adjuvant Chemodynamic Synergistic Therapy

Chemodynamic therapy based on the Fenton-like catalysis ability of Fe3O4 has the advantages of no involvement of chemical drugs and minimal adverse effects as well as the limitation of depletable efficacy. Radiotherapy based on high-energy radiation offers the convenience of treatment and cost-effectiveness but lacks precision and cellular adaptation of tumor cells. Approaching such dilemmas from a nanoscale materials perspective, we aim to bridge the weaknesses of both treatment methods by combining the principles of two therapeutics reciprocally. We have designed a camouflaged Fe3O4@HfO2 composite nanoreactor (FHCM), which combines a chemodynamic therapeutic agent Fe3O4 and a radiosensitizer HfO2 that both has passed clinical trials and was inspired by a cell membrane biomimetic technique. FHCM is employed as conceived radiotherapy-adjuvant chemodynamic synergistic therapy of malignant tumors, which has undergone dual scrutiny from both the physical and biological aspects. Experimental results obtained at different levels, including theory, material characterizations, and in vitro and in vivo verifications, suggest that FHCM effectively impaired tumor cells through physical and molecular biological mechanisms involving a HfO2-Fe3O4 photoelectron-electron transfer chain and DNA damage-ferroptosis-immunity chain. It is worth noting that compared to single therapies such as only chemodynamic therapy or radiotherapy, FHCM-mediated radiotherapy-adjuvant chemodynamic synergistic therapy exhibits stronger tumor inhibition efficacy. It significantly addresses the inherent limitations of chemodynamic therapy and radiotherapy and underscores the feasibility and importance of using existing clinical weapons, such as radiotherapy, as auxiliary strategies to overcome certain flaws of emerging antitumor therapeutics like chemodynamic therapy.

Deep learning-based lung image registration: A review

Lung image registration can effectively describe the relative motion of lung tissues, thereby helping to solve series problems in clinical applications. Since the lungs are soft and fairly passive organs, they are influenced by respiration and heartbeat, resulting in discontinuity of lung motion and large deformation of anatomic features. This poses great challenges for accurate registration of lung image and its applications. The recent application of deep learning (DL) methods in the field of medical image registration has brought promising results. However, a versatile registration framework has not yet emerged due to diverse challenges of registration for different regions of interest (ROI). DL-based image registration methods used for other ROI cannot achieve satisfactory results in lungs. In addition, there are few review articles available on DL-based lung image registration. In this review, the development of conventional methods for lung image registration is briefly described and a more comprehensive survey of DL-based methods for lung image registration is illustrated. The DL-based methods are classified according to different supervision types, including fully-supervised, weakly-supervised and unsupervised. The contributions of researchers in addressing various challenges are described, as well as the limitations of these approaches. This review also presents a comprehensive statistical analysis of the cited papers in terms of evaluation metrics and loss functions. In addition, publicly available datasets for lung image registration are also summarized. Finally, the remaining challenges and potential trends in DL-based lung image registration are discussed.

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

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

Dosimetric impact of deformable image registration using radiophotoluminescent glass dosimeters with a physical geometric phantom

PURPOSE

To study the dosimetry impact of deformable image registration (DIR) using radiophotoluminescent glass dosimeter (RPLD) and custom developed phantom with various inserts.

METHODS

The phantom was developed to facilitate simultaneous evaluation of geometric and dosimetric accuracy of DIR. Four computed tomography (CT) images of the phantom were acquired with four different configurations. Four volumetric modulated arc therapy (VMAT) plans were computed for different phantom. Two different patterns were applied to combination of four phantom configurations. RPLD dose measurement was combined between corresponding two phantom configurations. DIR-based dose accumulation was calculated between corresponding two CT images with two commercial DIR software and various DIR parameter settings, and an open source software. Accumulated dose calculated using DIR was then compared with measured dose using RPLD.

RESULTS

The mean ± standard deviation (SD) of dose difference was 2.71 ± 0.23% (range, 2.22%-3.01%) for tumor-proxy and 3.74 ± 0.79% (range, 1.56%-4.83%) for rectum-proxy. The mean ± SD of target registration error (TRE) was 1.66 ± 1.36 mm (range, 0.03-4.43 mm) for tumor-proxy and 6.87 ± 5.49 mm (range, 0.54-17.47 mm) for rectum-proxy. These results suggested that DIR accuracy had wide range among DIR parameter setting.

CONCLUSIONS

The dose difference observed in our study was 3% for tumor-proxy and within 5% for rectum-proxy. The custom developed physical phantom with inserts showed potential for accurate evaluation of DIR-based dose accumulation. The prospect of simultaneous evaluation of geometric and dosimetric DIR accuracy in a single phantom may be useful for validation of DIR for clinical use.

Dynamic Characteristics and Predictive Capability of Tumor Voxel Dose-Response Assessed Using 18F-FDG PET/CT Imaging Feedback

Purpose

Tumor voxel dose–response matrix (DRM) can be quantified using feedback from serial FDG-PET/CT imaging acquired during radiotherapy. This study investigated the dynamic characteristics and the predictive capability of DRM.

Methods

FDG-PET/CT images were acquired before and weekly during standard chemoradiotherapy with the treatment dose 2 Gy × 35 from 31 head and neck cancer patients. For each patient, deformable image registration was performed between the pretreatment/baseline PET/CT image and each weekly PET/CT image. Tumor voxel DRM was derived using linear regression on the logarithm of the weekly standard uptake value (SUV) ratios for each tumor voxel, such as SUV measured at a dose level normalized to the baseline SUV₀. The dynamic characteristics were evaluated by comparing the DRM_i estimated using a single feedback image acquired at the ith treatment week (i = 1, 2, 3, or 4) to the DRM estimated using the last feedback image for each patient. The predictive capability of the DRM estimated using 1 or 2 feedback images was evaluated using the receiver operating characteristic test with respect to the treatment outcome of tumor local–regional control or failure.

Results

The mean ± SD of tumor voxel SUV measured at the pretreatment and the 1st, 2nd, 3rd, 4th, and last treatment weeks was 6.76 ± 3.69, 5.72 ± 3.43, 3.85 ± 2.22, 3.27 ± 2.25, 2.5 ± 1.79, and 2.23 ± 1.27, respectively. The deviations between the DRM_i estimated using the single feedback image obtained at the ith week and the last feedback image were 0.86 ± 4.87, −0.06 ± 0.3, −0.09 ± 0.17, and −0.09 ± 0.12 for DRM₁, DRM₂, DRM₃, and DRM₄, respectively. The predictive capability of DRM₃ and DRM₄ was significant (p < 0.001). The area under the curve (AUC) was increased with the increase in treatment dose level. The DRMs constructed using the single feedback image achieved an AUC of 0.86~1. The AUC was slightly improved to 0.94~1 for the DRMs estimated using 2 feedback images.

Conclusion

Tumor voxel metabolic activity measured using FDG-PET/CT fluctuated noticeably during the first 2 treatment weeks and obtained a stabilized reduction rate thereafter. Tumor voxel DRM constructed using a single FDG-PET/CT feedback image after the 2nd treatment week (>20 Gy) has a good predictive capability. The predictive capability improved continuously using a later feedback image and marginally improved when two feedback images were applied.

Evaluation of DIR schemes on tumor/organ with progressive shrinkage: accuracy of tumor/organ internal tissue tracking during the radiation treatment

PURPOSE

Accuracy of intratumoral treatment dose accumulation and response assessment highly depends on the accuracy of a DIR method. However, achievable accuracy of the existing DIR methods for tumor/organ with large and progressive shrinkage during the radiotherapy course have not been explored. This study aimed to use a bio-tissue phantom to quantify the achievable accuracy of different DIR schemes.

MATERIALS /METHODS

A fresh porcine liver was used for phantom material. Sixty gold markers were implanted on the surface and inside of the liver. To simulate the progressive radiation-induced tumor/organ shrinkage, the phantom was heated using a microwave oven incrementally from 30s to 200s in 8 phases. For each phase, the phantom was scanned by CT. Two extra image sets were generated from the original images: 1) the image set with overriding the high-density gold markers (feature image); 2) the image set with overriding the entire phantom to the mean soft tissue intensity (featureless image). Ten DIR schemes were evaluated to mimic clinical treatment situations of tumor/critical organ with respect to their surface and internal condition of image features, availability of intermediate feedback images and DIR methods. The internal marker's positions were utilized to evaluate DIR accuracy quantified by target registration error (TRE).

RESULTS

Volume reduction was about 20% ∼ 40% of the initial volume after 90s ∼ 200s of the heating. Without image features on the surface and inside of the phantom, the hybrid-DIR (image-based DIR followed by biomechanical model-based refinement) with the surface constraint achieved the registration TRE from 2.6 ± 1.2mm to 5.3 ± 2.6mm proportional to the %volume shrinkage. Meanwhile, the hybrid-DIR with the surface-marker constraint achieved the TRE from 2.4 ± 1.2mm to 2.6 ± 1.0mm. If both the surface and internal image features would be viable on the feedback images, the achievable accuracy could be minimal with the TRE from 1.6±0.9mm to 1.9 ± 1.2mm.

CONCLUSIONS

Standard DIR methods cannot guarantee intratumoral tissue registration accuracy for tumor/organ with large progressive shrinkage. Achievable accuracy with using the hybrid DIR method is highly dependent on the surface registration accuracy. If the surface registration mean TRE can be controlled within 2mm, the mean TRE of internal tissue can be controlled within 3mm.

Trajectory log analysis and cone-beam CT-based daily dose calculation to investigate the dosimetric accuracy of intensity-modulated radiotherapy for gynecologic cancer

This study evaluated unexpected dosimetric errors caused by machine control accuracy, patient setup errors, and patient weight changes/internal organ deformations. Trajectory log files for 13 gynecologic plans with seven‐ or nine‐beam dynamic multileaf collimator (MLC) intensity‐modulated radiation therapy (IMRT), and differences between expected and actual MLC positions and MUs were evaluated. Effects of patient setup errors on dosimetry were estimated by in‐house software. To simulate residual patient setup errors after image‐guided patient repositioning, planned dose distributions were recalculated (blurred dose) after the positions were randomly moved in three dimensions 0–2 mm (translation) and 0°–2° (rotation) 28 times per patient. Differences between planned and blurred doses in the clinical target volume (CTV) D_98% and D_2% were evaluated. Daily delivered doses were calculated from cone‐beam computed tomography by the Hounsfield unit‐to‐density conversion method. Fractional and accumulated dose differences between original plans and actual delivery were evaluated by CTV D_98% and D_2%. The significance of accumulated doses was tested by the paired t test. Trajectory log file analysis showed that MLC positional errors were −0.01 ± 0.02 mm and MU delivery errors were 0.10 ± 0.10 MU. Differences in CTV D_98% and D_2% were <0.5% for simulated patient setup errors. Differences in CTV D_98% and D_2% were 2.4% or less between the fractional planned and delivered doses, but were 1.7% or less for the accumulated dose. Dosimetric errors were primarily caused by patient weight changes and internal organ deformation in gynecologic radiation therapy.

Effect of uncertainties in quantitative 18 F-FDG PET/CT imaging feedback for intratumoral dose-response assessment and dose painting by number

PURPOSE

Intratumoral dose response can be detected using serial fluoro-2-deoxyglucose-(FDG) positron emission tomography (PET)/computed tomography (CT) imaging feedback during treatment and used to guide adaptive dose painting by number (DPbN). However, to reliably implement this technique, the effect of uncertainties in quantitative PET/CT imaging feedback on tumor voxel dose-response assessment and DPbN needs to be determined and reduced.

METHODS

Three major uncertainties, induced by (a) PET imaging partial volume effect (PVE) and (b) tumor deformable image registration (DIR), and (c) variation of the time interval between FDG injection and PET image acquisition (TI), were determined using serial FDG-PET/CT images acquired during chemoradiotherapy of 18 head and neck cancer patients. PET imaging PVE was simulated using the discrepancy between with and without iterative deconvolution-based PVE corrections. Effect of tumor DIR uncertainty was simulated using the discrepancy between two DIR algorithms, including one with and one without soft-tissue mechanical correction for the voxel displacement. The effect of TI variation was simulated using linear interpolation on the dual-point PET/CT images. Tumor voxel pretreatment metabolic activity (SUV₀ ) and dose-response matrix (DRM) discrepancies induced by each of the three uncertainties were quantified, respectively. Adverse effects of tumor voxel SUV₀ and DRM discrepancies on tumor control probability (TCP) in DPbN were assessed.

RESULTS

Partial volume effect and TI variations of 10 mins induced a mean ± standard deviation (SD) of tumor voxel SUV₀ discrepancies to be -0.7% ± 9.2% and 0% ± 4.8%, respectively. Tumor voxel DRM discrepancies induced by PVE, tumor DIR discrepancy, and TI variations were 0.6% ± 8.9%, 1.7% ± 9.1%, and 0% ± 7%, respectively. Partial volume effect induced SUV₀ and DRM discrepancies correlated significantly with the tumor shape and FDG uptake heterogeneity. Tumor DIR uncertainty-induced DRM discrepancy correlated significantly with the tumor volume and shrinkage during treatment. Among the three uncertainties, PVE dominated the adverse effects on the TCP, with a mean ± SD of TCP reduction to be 12.7% ± 9.8% for all tumors if no compensation was applied for.

CONCLUSIONS

Effect of uncertainties in quantitative FDG-PET/CT imaging feedback on intratumoral dose-response quantification was not negligible. These uncertainties primarily caused by PVE and tumor DIR were highly dependent on individual tumor shape, volume, shrinkage during treatment, and pretreatment SUV heterogeneity, which can be managed individually. The adverse effects of these uncertainties could be minimized by using proper PVE corrections and DIR methods and compensated for in the clinical implementation of DPbN.

Influence of different urinary bladder filling levels and controlling regions of interest selection on deformable image registration algorithms

PURPOSE

Evaluation of Raystation ANAtomically CONstrained Deformation Algorithm (ANACONDA) performance to different urinary bladder filling levels in male pelvis anatomic site varying the controlling Regions Of Interest (ROIs).

METHODS

Different image datasets were obtained with ImSimQA (Oncology System Limited, Shrewsbury, UK) to evaluate ANACONDA performances (RaySearch Laboratories, Stockholm, Sweden). Deformation vector fields were applied to a synthetic man pelvis and a real patient computed tomography (CT) dataset (reference CTs) resulting in deformed CTs (target CTs) with various bladder filling levels. Different deformable image registrations (DIRs) were generated between each target CTs and reference CTs varying the controlling ROIs subset. Deformed ROIs were mapped from target CT to reference CT and then compared to reference ROIs. Evaluation was performed by Dice Similarity Coefficient (DSC), Correlation Coefficient (CC), Mean Distance to Agreement (MDA), maximum Distance to Agreement (maxDA) and with the introduction of global DSC (global_DSC) and global CC (global_CC) parameters.

RESULTS

In both synthetic and real patient CT cases, DSC scored less than 0.75 and MDA greater than 3 mm when no ROIs or only bladder were exploited as controlling ROI. DSC and CC increased by increasing the number of controlling ROIs selected whereas, an opposite behavior was observed for MDA and maxDA.

CONCLUSIONS

ANACONDA performances can be influenced by bladder filling fluctuation if no controlling ROIs are selected. Global_DSC and global_CC are useful parameters to quantitatively compare DIR algorithms. DIR performances improve by increasing the number of controlling ROIs selected, reaching a saturation level after a defined ROIs subset selection.

Modeling Complex Deformations of the Sigmoid Colon Between External Beam Radiation Therapy and Brachytherapy Images of Cervical Cancer

PURPOSE

In this study, we investigated registration methods for estimating the large interfractional sigmoid deformations that occur between external beam radiation therapy (EBRT) and brachytherapy (BT) for cervical cancer.

METHODS AND MATERIALS

Sixty-three patients were retrospectively analyzed. The sigmoid colon was delineated on 2 computed tomography images acquired during EBRT (without applicator) and BT (with applicator) for each patient. Five registration approaches were compared to propagate the contour of the sigmoid from BT to EBRT anatomies: rigid registration, commercial hybrid (ANAtomically CONstrained Deformation Algorithm), controlling ROI surface projection of RayStation, and the classical and constrained symmetrical thin-plate spline robust point matching (sTPS-RPM) methods. Deformation of the sigmoid due to insertion of the BT applicator was reported. Registration performance was compared by using the Dice similarity coefficient (DSC), distance to agreement, and Hausdorff distance. The 2 sTPS-RPM methods were compared by using surface triangle quality criteria between deformed surfaces. Using the deformable approaches, the BT dose of the sigmoid was deformed toward the EBRT anatomy. The displacement and discrepancy between the deformable methods to propagate the planned D1cm³ and D2cm³ of the sigmoid from BT to EBRT anatomies were reported for 55 patients.

RESULTS

Large and complex deformations of the sigmoid were observed for each patient. Rigid registration resulted in poor sigmoid alignment with a mean DSC of 0.26. Using the contour to drive the deformation, ANAtomically CONstrained Deformation Algorithm was able to slightly improve the alignment of the sigmoid with a mean DSC of 0.57. Using only the sigmoid surface as controlling ROI, the mean DSC was improved to 0.79. The classical and constrained sTPS-RPM methods provided mean DSCs of 0.95 and 0.96, respectively, with an average inverse consistency error <1 mm. The constrained sTPS-RPM provided more realistic deformations and better surface topology of the deformed sigmoids. The planned mean (range) D1cm³ and D2cm³ of the sigmoid were 13.4 Gy (1-24.1) and 12.2 Gy (1-21.5) on the BT anatomy, respectively. Using the constrained sTPS-RPM to deform the sigmoid from BT to EBRT anatomies, these hotspots had a mean (range) displacement of 27.1 mm (6.8-81).

CONCLUSIONS

Large deformations of the sigmoid were observed between the EBRT and BT anatomies, suggesting that the D1cm³ and D2cm³ of the sigmoid would unlikely to be at the same position throughout treatment. The proposed constrained sTPS-RPM seems to be the preferred approach to manage the large deformation due to BT applicator insertion. Such an approach could be used to map the EBRT dose to the BT anatomy for personalized BT planning optimization.

Differences between planned and delivered dose for head and neck cancer, and their consequences for normal tissue complication probability and treatment adaptation

BACKGROUND AND PURPOSE

Anatomical changes induce differences between planned and delivered dose. Adaptive radiotherapy (ART) may reduce these differences but the optimal implementation is insufficiently clear. The aims of this study were to quantify the difference between planned and delivered dose in HNC patients, assess the consequential difference in normal tissue complication probability (ΔNTCP) and to explore the value of ΔNTCP as an objective selection strategy for ART.

MATERIALS AND METHODS

For 52 patients, daily doses were accumulated to estimate the delivered dose. The difference from planned dose was analyzed for CTVs and 9 organs-at-risk (OAR). ΔNTCP was calculated for xerostomia, dysphagia, parotid gland dysfunction and tube feeding dependency at 6 months. ART was deemed necessary if ΔNTCP was >5%. The positive predicted value (PPV) was calculated for identification of ART-patients by clinical judgement, and ΔNTCP at fraction 10 and 15.

RESULTS

ΔNTCP >5% was seen five times for dysphagia and twice for the other toxicities. Only 5/9 patients with any ΔNTCP >5% clinically received ART, although ART had been done for 13/52 patients (PPV: 0.38). PPV was 0.86 and 0.75 for accumulated dose at fraction 10 and 15, respectively, using a ΔNTCP cut-off for the allocation of ART of 5%. Using other ΔNTCP cut-offs did not substantially improve PPV. With this cut-off the negative predictive value was 0.93 for ΔNTCP method of fraction 10 and fraction 15, and 0.90 for clinical judgement.

CONCLUSION

To identify patients accurately for ART, NTCP calculations based on the dose differences between planned and delivered dose at fraction 10 are superior to clinical judgement.

Deformable image registration for radiation therapy: principle, methods, applications and evaluation

Background: Deformable image registration (DIR) is increasingly used in the field of radiation therapy (RT) to account for anatomical deformations. The aims of this paper are to describe the main applications of DIR in RT and discuss current DIR evaluation methods. Methods: Articles on DIR published from January 2000 to October 2018 were extracted from PubMed and Science Direct. Our search was restricted to articles that report data obtained from humans, were written in English, and address DIR methods for RT. A total of 207 articles were selected from among 2506 identified in the search process. Results: At planning, DIR is used for organ delineation using atlas-based segmentation, deformation-based planning target volume definition, functional planning and magnetic resonance imaging-based dose calculation. In image-guided RT, DIR is used for contour propagation and dose calculation on per-treatment imaging. DIR is also used to determine the accumulated dose from fraction to fraction in external beam RT and brachytherapy, both for dose reporting and adaptive RT. In the case of re-irradiation, DIR can be used to estimate the cumulated dose of the two irradiations. Finally, DIR can be used to predict toxicity in voxel-wise population analysis. However, the evaluation of DIR remains an open issue, especially when dealing with complex cases such as the disappearance of matter. To quantify DIR uncertainties, most evaluation methods are limited to geometry-based metrics. Software companies have now integrated DIR tools into treatment planning systems for clinical use, such as contour propagation and fraction dose accumulation. Conclusions: DIR is increasingly important in RT applications, from planning to toxicity prediction. DIR is routinely used to reduce the workload of contour propagation. However, its use for complex dosimetric applications must be carefully evaluated by combining quantitative and qualitative analyses.

Head and Neck Cancer Adaptive Radiation Therapy (ART): Conceptual Considerations for the Informed Clinician

For nearly 2 decades, adaptive radiation therapy (ART) has been proposed as a method to account for changes in head and neck tumor and normal tissue to enhance therapeutic ratios. While technical advances in imaging, planning and delivery have allowed greater capacity for ART delivery, and a series of dosimetric explorations have consistently shown capacity for improvement, there remains a paucity of clinical trials demonstrating the utility of ART. Furthermore, while ad hoc implementation of head and neck ART is reported, systematic full-scale head and neck ART remains an as yet unreached reality. To some degree, this lack of scalability may be related to not only the complexity of ART, but also variability in the nomenclature and descriptions of what is encompassed by ART. Consequently, we present an overview of the history, current status, and recommendations for the future of ART, with an eye toward improving the clarity and description of head and neck ART for interested clinicians, noting practical considerations for implementation of an ART program or clinical trial. Process level considerations for ART are noted, reminding the reader that, paraphrasing the writer Elbert Hubbard, "Art is not a thing, it is a way."

The evaluation of a hybrid biomechanical deformable registration method on a multistage physical phantom with reproducible deformation

BACKGROUND

Advanced clinical applications, such as dose accumulation and adaptive radiation therapy, require deformable image registration (DIR) algorithms capable of voxel-wise accurate mapping of treatment dose or functional imaging. By utilizing a multistage deformable phantom, the authors investigated scenarios where biomechanical refinement method (BM-DIR) may be better than the pure image intensity based deformable registration (IM-DIR).

METHODS

The authors developed a biomechanical-model based DIR refinement method (BM-DIR) to refine the deformable vector field (DVF) from any initial intensity-based DIR (IM-DIR). The BM-DIR method was quantitatively evaluated on a novel phantom capable of ten reproducible gradually-increasing deformation stages using the urethra tube as a surrogate. The internal DIR accuracy was inspected in term of the Dice similarity coefficient (DSC), Hausdorff and mean surface distance as defined in of the urethra structure inside the phantom and compared with that of the initial IM-DIR under various stages of deformation. Voxel-wise deformation vector discrepancy and Jacobian regularity were also inspected to evaluate the output DVFs. In addition to phantom, two pairs of Head&Neck patient MR images with expert-defined landmarks inside parotids were utilized to evaluate the BM-DIR accuracy with target registration error (TRE).

RESULTS

The DSC and surface distance measures of the inner urethra tube indicated the BM-DIR method can improve the internal DVF accuracy on masked MR images for the phases of a large degree of deformation. The smoother Jacobian distribution from the BM-DIR suggests more physically-plausible internal deformation. For H&N cancer patients, the BM-DIR improved the TRE from 0.339 cm to 0.210 cm for the landmarks inside parotid on the masked MR images.

CONCLUSIONS

We have quantitatively demonstrated on a multi-stage physical phantom and limited patient data that the proposed BM-DIR can improve the accuracy inside solid organs with large deformation where distinctive image features are absent.

A clinical 3D/4D CBCT-based treatment dose monitoring system

To monitor delivered dose and trigger plan adaptation when deviation becomes unacceptable, a clinical treatment dose (Tx‐Dose) reconstruction system based on three‐dimensional (3D)/four‐dimensional (4D)‐cone beam computed tomograpy (CBCT) images was developed and evaluated on various treatment sites, particularly for lung cancer patient treated by stereotactic body radiation therapy (SBRT). This system integrates with our treatment planning system (TPS), Linacs recording and verification system (R&V), and CBCT imaging system, consisting of three modules: Treatment Schedule Monitoring module (TSM), pseudo‐CT Generating module (PCG), and Treatment Dose Reconstruction/evaluation module (TDR). TSM watches the treatment progress in the R&V system and triggers the PCG module when new CBCT is available. PCG retrieves the CBCTs and performs planning CT to CBCT deformable registration (DIR) to generate pseudo‐CT. The 4D‐CBCT images are taken for target localization and correction in lung cancer patient before treatment. To take full advantage of the valuable information carried by 4D‐CBCT, a novel phase‐matching DIR scheme was developed to generate 4D pseudo‐CT images for 4D dose reconstruction. Finally, TDR module creates TPS scripts to automate Tx‐Dose calculation on the pseudo‐CT images. Both initial quantitative commissioning and patient‐specific qualitative quality assurance of the DIR tool were utilized to ensure the DIR quality. The treatment doses of ten patients (six SBRT‐lung, two head and neck (HN), one breast and one prostate cancer patients) were retrospectively constructed and evaluated. The target registration error (mean ± STD: 1.05 ± 1.13 mm) of the DIR tool is comparable to the interobserver uncertainty (0.88 ± 1.31 mm) evaluated by a publically available lung‐landmarks dataset. For lung SBRT patients, the D₉₉ of the final cumulative Tx‐Dose of GTV is 93.8 ± 5.5% (83.7–100.1%) of the originally planned D₉₉. CTV D₉₉ decreases by 3% and mean ipsilateral parotid dose increases by 11.5% for one of the two HN patients. In conclusion, we have demonstrated the feasibility and effectiveness of a treatment dose verification system in our clinical setting.

Practical quantification of image registration accuracy following the AAPM TG-132 report framework

The AAPM TG 132 Report enumerates important steps for validation of the medical image registration process. While the Report outlines the general goals and criteria for the tests, specific implementation may be obscure to the wider clinical audience. We endeavored to provide a detailed step‐by‐step description of the quantitative tests’ execution, applied as an example to a commercial software package (Mirada Medical, Oxford, UK), while striving for simplicity and utilization of readily available software. We demonstrated how the rigid registration data could be easily extracted from the DICOM registration object and used, following some simple matrix math, to quantify accuracy of rigid translations and rotations. The options for validating deformable image registration (DIR) were enumerated, and it was shown that the most practically viable ones are comparison of propagated internal landmark points on the published datasets, or of segmented contours that can be generated locally. The multimodal rigid registration in our example did not always result in the desired registration error below ½ voxel size, but was considered acceptable with the maximum errors under 1.3 mm and 1°. The DIR target registration errors in the thorax based on internal landmarks were far in excess of the Report recommendations of 2 mm average and 5 mm maximum. On the other hand, evaluation of the DIR major organs’ contours propagation demonstrated good agreement for lung and abdomen (Dice Similarity Coefficients, DSC, averaged over all cases and structures of 0.92 ± 0.05 and 0.91 ± 0.06, respectively), and fair agreement for Head and Neck (average DSC = 0.73 ± 0.14). The average for head and neck is reduced by small volume structures such as pharyngeal constrictor muscles. Even these relatively simple tests show that commercial registration algorithms cannot be automatically assumed sufficiently accurate for all applications. Formalized task‐specific accuracy quantification should be expected from the vendors.