Use of image registration and fusion algorithms and techniques in radiotherapy: Report of the AAPM Radiation Therapy Committee Task Group No. 132

Impact of CT scan parameters on deformable image registration accuracy using deformable thorax phantom

The purpose of this study was to evaluate the deformable image registration (DIR) accuracy using various CT scan parameters with deformable thorax phantom. Our developed deformable thorax phantom (Dephan, Chiyoda Technol Corp, Tokyo, Japan) was used. The phantom consists of a base phantom, an inner phantom, and a motor-derived piston. The base phantom is an acrylic cylinder phantom with a diameter of 180 mm, which simulates the chest wall. The inner phantom consists of deformable, 20 mm thick disk-shaped sponges with 48 Lucite beads and 48 nylon cross-wires which simulate the vascular and bronchial bifurcations of the lung. Peak-exhale and peak-inhale images of the deformable phantom were acquired using a CT scanner (Aquilion LB, TOSHIBA). To evaluate the impact of CT scan parameters on DIR accuracy, we used the four tube voltages (80, 100, 120, and 135 kV) and six reconstruction algorithms (FC11, FC13, FC15, FC41, FC44, and FC52). Intensity-based DIR was performed between the two images using MIM Maestro (MIM software, Cleveland, USA). Fiducial markers (beads and cross-wires) based target registration error (TRE) was used for quantitative evaluation of DIR. In case with different tube voltages, the range of average TRE were 4.44-5.69 mm (reconstruction algorithm: FC13). In case with different reconstruction algorithms, the range of average TRE were 4.26-4.59 mm (tube voltage: 120 kV). The TRE were differed by up to 3.0 mm (3.96-6.96 mm) depending on the combination of tube voltage and reconstruction algorithm. Our result indicated that CT scan parameters had moderate impact of TRE, especially for reconstruction algorithms for the deformable thorax phantom.

Deep learning based automatic contour refinement for inaccurate auto-segmentation in MR-guided adaptive radiotherapy

Objective.Fast and accurate auto-segmentation is essential for magnetic resonance-guided adaptive radiation therapy (MRgART). Deep learning auto-segmentation (DLAS) is not always clinically acceptable, particularly for complex abdominal organs. We previously reported an automatic contour refinement (ACR) solution of using an active contour model (ACM) to partially correct the DLAS contours. This study aims to develop a DL-based ACR model to work in conjunction with ACM-ACR to further improve the contour accuracy.Approach.The DL-ACR model was trained and tested using bowel contours created by an in-house DLAS system from 160 MR sets (76 from MR-simulation and 84 from MR-Linac). The contours were classified into acceptable, minor-error and major-error groups using two approaches of contour quality classification (CQC), based on the AAPM TG-132 recommendation and an in-house classification model, respectively. For the major-error group, DL-ACR was applied subsequently after ACM-ACR to further refine the contours. For the minor-error group, contours were directly corrected by DL-ACR without applying an initial ACM-ACR. The ACR workflow was performed separately for the two CQC methods and was evaluated using contours from 25 image sets as independent testing data.Main results.The best ACR performance was observed in the MR-simulation testing set using CQC by TG-132: (1) for the major-error group, 44% (177/401) were improved to minor-error group and 5% (22/401) became acceptable by applying ACM-ACR; among these 177 contours that shifted from major-error to minor-error with ACM-ACR, DL-ACR further refined 49% (87/177) to acceptable; and overall, 36% (145/401) were improved to minor-error contours, and 30% (119/401) became acceptable after sequentially applying ACM-ACR and DL-ACR; (2) for the minor-error group, 43% (320/750) were improved to acceptable contours using DL-ACR.Significance.The obtained ACR workflow substantially improves the accuracy of DLAS bowel contours, minimizing the manual editing time and accelerating the segmentation process of MRgART.

The relationship between radiation dose and bevacizumab-related imaging abnormality in patients with brain tumors: A voxel-wise normal tissue complication probability (NTCP) analysis

PURPOSE

Bevacizumab-related imaging abnormality (BRIA), appearing as areas of restricted diffusion on magnetic resonance imaging (MRI) and representing atypical coagulative necrosis pathologically, has been observed in patients with brain tumors receiving radiotherapy and bevacizumab. We investigated the role of cumulative radiation dose in BRIA development in a voxel-wise analysis.

METHODS

Patients (n = 18) with BRIA were identified. All had high-grade gliomas or brain metastases treated with radiotherapy and bevacizumab. Areas of BRIA were segmented semi-automatically on diffusion-weighted MRI with apparent diffusion coefficient (ADC) images. To avoid confounding by possible tumor, hypoperfusion was confirmed with perfusion imaging. ADC images and radiation dose maps were co-registered to a high-resolution T1-weighted MRI and registration accuracy was verified. Voxel-wise normal tissue complication probability analyses were performed using a logistic model analyzing the relationship between cumulative voxel equivalent total dose in 2 Gy fractions (EQD2) and BRIA development at each voxel. Confidence intervals for regression model predictions were estimated with bootstrapping.

RESULTS

Among 18 patients, 39 brain tumors were treated. Patients received a median of 4.5 cycles of bevacizumab and 1-4 radiation courses prior to BRIA appearance. Most (64%) treated tumors overlapped with areas of BRIA. The median proportion of each BRIA region of interest volume overlapping with tumor was 98%. We found a dose-dependent association between cumulative voxel EQD2 and the relative probability of BRIA (β0 = -5.1, β1 = 0.03 Gy-1, γ = 1.3).

CONCLUSIONS

BRIA is likely a radiation dose-dependent phenomenon in patients with brain tumors receiving bevacizumab and radiotherapy. The combination of radiation effects and tumor microenvironmental factors in potentiating BRIA in this population should be further investigated.

Evaluating on-board kVCT- and MVCT-based dose calculation accuracy using a thorax phantom for helical tomotherapy treatments

Purpose.To evaluate the impact of CT number calibration and imaging parameter selection on dose calculation accuracy relative to the CT planning process in thoracic treatments for on-board helical CT imaging systems used in helical tomotherapy.Methods and Materials.Direct CT number calibrations were performed with appropriate protocols for each imaging system using an electron density phantom. Large volume and SBRT treatment plans were simulated and optimized for planning CT scans of an anthropomorphic thorax phantom and transferred to registered kVCT and MVCT scans of the phantom as appropriate. Relevant DVH metrics and dose-difference maps were used to evaluate and compare dose calculation accuracy relative to the planning CT based on a variation in imaging parameters applied for the on-board systems.Results.For helical kVCT scans of the thorax phantom, median differences in DVH parameters for the large volume treatment plan were less than ±1% with dose to the target volume either over- or underestimated depending on the imaging parameters utilized for CT number calibration and thorax phantom acquisition. For the lung SBRT plan calculated on helical kVCT scans, median dose differences were up to -2.7% with a more noticeable dependence on parameter selection. For MVCT scans, median dose differences for the large volume plan were within +2% with dose to the target overestimated regardless of the imaging protocol.Conclusion.Accurate dose calculations (median errors of <±1%) using a thorax phantom simulating realistic patient geometry and scatter conditions can be achieved with images acquired with a helical kVCT system on a helical tomotherapy unit. This accuracy is considerably improved relative to that achieved with the MV-based approach. In a clinical setting, careful consideration should be made when selecting appropriate kVCT imaging parameters for this process as dose calculation accuracy was observed to vary with both parameter selection and treatment type.

Applying Multi-Metric Deformable Image Registration for Dose Accumulation in Combined Cervical Cancer Radiotherapy

(1) Purpose: Challenges remain in dose accumulation for cervical cancer radiotherapy combined with external beam radiotherapy (EBRT) and brachytherapy (BT) as there are many large and complex organ deformations between different treatments. This study aims to improve deformable image registration (DIR) accuracy with the introduction of multi-metric objectives for dose accumulation of EBRT and BT. (2) Materials and methods: Twenty cervical cancer patients treated with EBRT (45-50 Gy/25 fractions) and high-dose-rate BT (≥20 Gy in 4 fractions) were included for DIR. The multi-metric DIR algorithm included an intensity-based metric, three contour-based metrics, and a penalty term. Nonrigid B-spine transformation was used to transform the planning CT images from EBRT to the first BT, with a six-level resolution registration strategy. To evaluate its performance, the multi-metric DIR was compared with a hybrid DIR provided by commercial software. The DIR accuracy was measured by the Dice similarity coefficient (DSC) and Hausdorff distance (HD) between deformed and reference organ contours. The accumulated maximum dose of 2 cc (D_2cc) of the bladder and rectum was calculated and compared to simply addition of D_2cc from EBRT and BT (ΔD_2cc). (3) Results: The mean DSC of all organ contours for the multi-metric DIR were significantly higher than those for the hybrid DIR (p ≤ 0.011). In total, 70% of patients had DSC > 0.8 using the multi-metric DIR, while 15% of patients had DSC > 0.8 using the commercial hybrid DIR. The mean ΔD_2cc of the bladder and rectum for the multi-metric DIR were 3.25 ± 2.29 and 3.54 ± 2.02 Gy_EQD2, respectively, whereas those for the hybrid DIR were 2.68 ± 2.56 and 2.32 ± 3.25 Gy_EQD2, respectively. The multi-metric DIR resulted in a much lower proportion of unrealistic D_2cc than the hybrid DIR (2.5% vs. 17.5%). (4) Conclusions: Compared with the commercial hybrid DIR, the introduced multi-metric DIR significantly improved the registration accuracy and resulted in a more reasonable accumulated dose distribution.

Dosimetric Impact of Intrafraction Prostate Motion and Interfraction Anatomical Changes in Dose-Escalated Linac-Based SBRT

The dosimetric impact of intrafraction prostate motion and interfraction anatomical changes and the effect of beam gating and motion correction were investigated in dose-escalated linac-based SBRT. Fifty-six gated fractions were delivered using a novel electromagnetic tracking device with a 2 mm threshold. Real-time prostate motion data were incorporated into the patient's original plan with an isocenter shift method. Delivered dose distributions were obtained by recalculating these motion-encoded plans on deformed CTs reflecting the patient's CBCT daily anatomy. Non-gated treatments were simulated using the prostate motion data assuming that no treatment interruptions have occurred. The mean relative dose differences between delivered and planned treatments were -3.0% [-18.5-2.8] for CTV D99% and -2.6% [-17.8-1.0] for PTV D95%. The median cumulative CTV coverage with 93% of the prescribed dose was satisfactory. Urethra sparing was slightly degraded, with the maximum dose increased by only 1.0% on average, and a mean reduction in the rectum and bladder doses was seen in almost all dose metrics. Intrafraction prostate motion marginally contributed in gated treatments, while in non-gated treatments, further deteriorations in the minimum target coverage and bladder dose metrics would have occurred on average. The implemented motion management strategy and the strict patient preparation regimen, along with other treatment optimization strategies, ensured no significant degradations of dose metrics in delivered treatments.

RMSim: controlled respiratory motion simulation on static patient scans

Objective.This work aims to generate realistic anatomical deformations from static patient scans. Specifically, we present a method to generate these deformations/augmentations via deep learning driven respiratory motion simulation that provides the ground truth for validating deformable image registration (DIR) algorithms and driving more accurate deep learning based DIR.Approach.We present a novel 3D Seq2Seq deep learning respiratory motion simulator (RMSim) that learns from 4D-CT images and predicts future breathing phases given a static CT image. The predicted respiratory patterns, represented by time-varying displacement vector fields (DVFs) at different breathing phases, are modulated through auxiliary inputs of 1D breathing traces so that a larger amplitude in the trace results in more significant predicted deformation. Stacked 3D-ConvLSTMs are used to capture the spatial-temporal respiration patterns. Training loss includes a smoothness loss in the DVF and mean-squared error between the predicted and ground truth phase images. A spatial transformer deforms the static CT with the predicted DVF to generate the predicted phase image. 10-phase 4D-CTs of 140 internal patients were used to train and test RMSim. The trained RMSim was then used to augment a public DIR challenge dataset for training VoxelMorph to show the effectiveness of RMSim-generated deformation augmentation.Main results.We validated our RMSim output with both private and public benchmark datasets (healthy and cancer patients). The structure similarity index measure (SSIM) for predicted breathing phases and ground truth 4D CT images was 0.92 ± 0.04, demonstrating RMSim's potential to generate realistic respiratory motion. Moreover, the landmark registration error in a public DIR dataset was improved from 8.12 ± 5.78 mm to 6.58mm ± 6.38 mm using RMSim-augmented training data.Significance.The proposed approach can be used for validating DIR algorithms as well as for patient-specific augmentations to improve deep learning DIR algorithms. The code, pretrained models, and augmented DIR validation datasets will be released athttps://github.com/nadeemlab/SeqX2Y.

Cardiac motion and its dosimetric impact during radioablation for refractory ventricular tachycardia

INTRODUCTION

Cardiac radioablation (CR) is a noninvasive treatment option for patients with refractory ventricular tachycardia (VT) during which high doses of radiation, typically 25 Gy, are delivered to myocardial scar. In this study, we investigate motion from cardiac cycle and evaluate the dosimetric impact in a cohort of patients treated with CR.

METHODS

This retrospective study included eight patients treated at our institution who had respiratory-correlated and ECG-gated 4DCT scans acquired within 2 weeks of CR. Deformable image registration was applied between maximum systole (SYS) and diastole (DIAS) CTs to assess cardiac motion. The average respiratory-correlated CT (AVG_resp ) was deformably registered to the average cardiac (AVG_cardiac ), SYS, and DIAS CTs, and contours were propagated using the deformation vector fields (DVFs). Finally, the original treatment plan was recalculated on the deformed AVG_resp CT for dosimetric assessment.

RESULTS

Motion magnitudes were measured as the mean (SD) value over the DVFs within each structure. Displacement during the cardiac cycle for all chambers was 1.4 (0.9) mm medially/laterally (ML), 1.6 (1.0) mm anteriorly/posteriorly (AP), and 3.0 (2.8) mm superiorly/inferiorly (SI). Displacement for the 12 distinct clinical target volumes (CTVs) was 1.7 (1.5) mm ML, 2.4 (1.1) mm AP, and 2.1 (1.5) SI. Displacements between the AVG_resp and AVG_cardiac scans were 4.2 (2.0) mm SI and 5.8 (1.4) mm total. Dose recalculations showed that cardiac motion may impact dosimetry, with dose to 95% of the CTV dropping from 27.0 (1.3) Gy on the AVG_resp to 20.5 (7.1) Gy as estimated on the AVG_cardiac .

CONCLUSIONS

Cardiac CTV motion in this patient cohort is on average below 3 mm, location-dependent, and when not accounted for in treatment planning may impact target coverage. Further study is needed to assess the impact of cardiac motion on clinical outcomes.

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.

Seq2Morph: A deep learning deformable image registration algorithm for longitudinal imaging studies and adaptive radiotherapy

PURPOSE

To simultaneously register all the longitudinal images acquired in a radiotherapy course for analyzing patients' anatomy changes for adaptive radiotherapy (ART).

METHODS

To address the unique needs of ART, we designed Seq2Morph, a novel deep learning-based deformable image registration (DIR) network. Seq2Morph was built upon VoxelMorph which is a general-purpose framework for learning-based image registration. The major upgrades are (1) expansion of inputs to all weekly cone-beam computed tomography (CBCTs) acquired for monitoring treatment responses throughout a radiotherapy course, for registration to their planning CT; (2) incorporation of 3D convolutional long short-term memory between the encoder and decoder of VoxelMorph, to parse the temporal patterns of anatomical changes; and (3) addition of bidirectional pathways to calculate and minimize inverse consistency errors (ICEs). Longitudinal image sets from 50 patients, including a planning CT and 6 weekly CBCTs per patient, were utilized for network training and cross-validation. The outputs were deformation vector fields for all the registration pairs. The loss function was composed of a normalized cross-correlation for image intensity similarity, a DICE for contour similarity, an ICE, and a deformation regularization term. For performance evaluation, DICE and Hausdorff distance (HD) for the manual versus predicted contours of tumor and esophagus on weekly basis were quantified and further compared with other state-of-the-art algorithms, including conventional VoxelMorph and large deformation diffeomorphic metric mapping (LDDMM).

RESULTS

Visualization of the hidden states of Seq2Morph revealed distinct spatiotemporal anatomy change patterns. Quantitatively, Seq2Morph performed similarly to LDDMM, but significantly outperformed VoxelMorph as measured by GTV DICE: (0.799±0.078, 0.798±0.081, and 0.773±0.078), and 50% HD (mm): (0.80±0.57, 0.88±0.66, and 0.95±0.60). The per-patient inference of Seq2Morph took 22 s, much less than LDDMM (∼30 min).

CONCLUSIONS

Seq2Morph can provide accurate and fast DIR for longitudinal image studies by exploiting spatial-temporal patterns. It closely matches the clinical workflow and has the potential to serve both online and offline ART.

The Impact of Interactive MRI-Based Radiologist Review on Radiotherapy Target Volume Delineation in Head and Neck Cancer

BACKGROUND AND PURPOSE

Peer review of head and neck cancer radiation therapy target volumes by radiologists was introduced in our center to optimize target volume delineation. Our aim was to assess the impact of MR imaging-based radiologist peer review of head and neck radiation therapy gross tumor and nodal volumes, through qualitative and quantitative analysis.

MATERIALS AND METHODS

Cases undergoing radical radiation therapy with a coregistered MR imaging, between April 2019 and March 2020, were reviewed. The frequency and nature of volume changes were documented, with major changes classified as per the guidance of The Royal College of Radiologists. Volumetric alignment was assessed using the Dice similarity coefficient, Jaccard index, and Hausdorff distance.

RESULTS

Fifty cases were reviewed between April 2019 and March 2020. The median age was 59 years (range, 29-83 years), and 72% were men. Seventy-six percent of gross tumor volumes and 41.5% of gross nodal volumes were altered, with 54.8% of gross tumor volume and 66.6% of gross nodal volume alterations classified as "major." Undercontouring of soft-tissue involvement and unidentified lymph nodes were predominant reasons for change. Radiologist review significantly altered the size of both the gross tumor volume (P  = .034) and clinical target tumor volume (P  = .003), but not gross nodal volume or clinical target nodal volume. The median conformity and surface distance metrics were the following: gross tumor volume Dice similarity coefficient = 0.93 (range, 0.82-0.96), Jaccard index = 0.87 (range, 0.7-0.94), Hausdorff distance = 7.45 mm (range, 5.6-11.7 mm); and gross nodular tumor volume Dice similarity coefficient = 0.95 (0.91-0.97), Jaccard index = 0.91 (0.83-0.95), and Hausdorff distance = 20.7 mm (range, 12.6-41.6). Conformity improved on gross tumor volume-to-clinical target tumor volume expansion (Dice similarity coefficient = 0.93 versus 0.95, P  = .003).

CONCLUSIONS

MR imaging-based radiologist review resulted in major changes to most radiotherapy target volumes and significant changes in volume size of both gross tumor volume and clinical target tumor volume, suggesting that this is a fundamental step in the radiotherapy workflow of patients with head and neck cancer.

Detailed analysis of failure patterns using deformable image registration in hypopharyngeal cancer patients treated with sequential boost intensity-modulated radiotherapy

INTRODUCTION

Sequential boost intensity-modulated radiotherapy (SQB-IMRT) uses two different planning CTs (pCTs) and treatment plans. SQB-IMRT is a form of adaptive radiotherapy that allows for responses to changes in the shape of the tumour and organs at risk (OAR). On the other hand, dose accumulation with the two plans can be difficult to evaluate. The purpose of this study was to analyse patterns of loco-regional failure using deformable image registration (DIR) in hypopharyngeal cancer patients treated with SQB-IMRT.

METHODS

Between 2013 and 2019, 102 patients with hypopharyngeal cancer were treated with definitive SQB-IMRT at our institution. Dose accumulation with the 1st and 2nd plans was performed, and the dose to the loco-regional recurrent tumour volume was calculated using the DIR workflow. Failure was classified as follows: (i) in-field (≥95% of the recurrent tumour volume received 95% of the prescribed dose); (ii) marginal (20-95%); or (iii) out-of-field (<20%).

RESULTS

After a median follow-up period of 25 months, loco-regional failure occurred in 34 patients. Dose-volume histogram analysis showed that all loco-regional failures occurred in the field within 95% of the prescribed dose, with no marginal or out-of-field recurrences observed.

CONCLUSION

The dosimetric analysis using DIR showed that all loco-regional failures were within the high-dose region. More aggressive treatment may be required for gross tumours.

A Review of the Metrics Used to Assess Auto-Contouring Systems in Radiotherapy

Auto-contouring could revolutionise future planning of radiotherapy treatment. The lack of consensus on how to assess and validate auto-contouring systems currently limits clinical use. This review formally quantifies the assessment metrics used in studies published during one calendar year and assesses the need for standardised practice. A PubMed literature search was undertaken for papers evaluating radiotherapy auto-contouring published during 2021. Papers were assessed for types of metric and the methodology used to generate ground-truth comparators. Our PubMed search identified 212 studies, of which 117 met the criteria for clinical review. Geometric assessment metrics were used in 116 of 117 studies (99.1%). This includes the Dice Similarity Coefficient used in 113 (96.6%) studies. Clinically relevant metrics, such as qualitative, dosimetric and time-saving metrics, were less frequently used in 22 (18.8%), 27 (23.1%) and 18 (15.4%) of 117 studies, respectively. There was heterogeneity within each category of metric. Over 90 different names for geometric measures were used. Methods for qualitative assessment were different in all but two papers. Variation existed in the methods used to generate radiotherapy plans for dosimetric assessment. Consideration of editing time was only given in 11 (9.4%) papers. A single manual contour as a ground-truth comparator was used in 65 (55.6%) studies. Only 31 (26.5%) studies compared auto-contours to usual inter- and/or intra-observer variation. In conclusion, significant variation exists in how research papers currently assess the accuracy of automatically generated contours. Geometric measures are the most popular, however their clinical utility is unknown. There is heterogeneity in the methods used to perform clinical assessment. Considering the different stages of system implementation may provide a framework to decide the most appropriate metrics. This analysis supports the need for a consensus on the clinical implementation of auto-contouring.

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

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

Evaluation of generalization ability for deep learning-based auto-segmentation accuracy in limited field of view CBCT of male pelvic region

PURPOSE

The aim of this study was to evaluate generalization ability of segmentation accuracy for limited FOV CBCT in the male pelvic region using a full-image CNN. Auto-segmentation accuracy was evaluated using various datasets with different intensity distributions and FOV sizes.

METHODS

A total of 171 CBCT datasets from patients with prostate cancer were enrolled. There were 151, 10, and 10 CBCT datasets acquired from Vero4DRT, TrueBeam STx, and Clinac-iX, respectively. The FOV for Vero4DRT, TrueBeam STx, and Clinac-iX was 20, 26, and 25 cm, respectively. The ROIs, including the bladder, prostate, rectum, and seminal vesicles, were manually delineated. The U² -Net CNN network architecture was used to train the segmentation model. A total of 131 limited FOV CBCT datasets from Vero4DRT were used for training (104 datasets) and validation (27 datasets); thereafter the rest were for testing. The training routine was set to save the best weight values when the DSC in the validation set was maximized. Segmentation accuracy was qualitatively and quantitatively evaluated between the ground truth and predicted ROIs in the different testing datasets.

RESULTS

The mean scores ± standard deviation of visual evaluation for bladder, prostate, rectum, and seminal vesicle in all treatment machines were 1.0 ± 0.7, 1.5 ± 0.6, 1.4 ± 0.6, and 2.1 ± 0.8 points, respectively. The median DSC values for all imaging devices were ≥0.94 for the bladder, 0.84-0.87 for the prostate and rectum, and 0.48-0.69 for the seminal vesicles. Although the DSC values for the bladder and seminal vesicles were significantly different among the three imaging devices, the DSC value of the bladder changed by less than 1% point. The median MSD values for all imaging devices were ≤1.2 mm for the bladder and 1.4-2.2 mm for the prostate, rectum, and seminal vesicles. The MSD values for the seminal vesicles were significantly different between the three imaging devices.

CONCLUSION

The proposed method is effective for testing datasets with different intensity distributions and FOV from training datasets.

Improving total body irradiation with a dedicated couch and 3D-printed patient-specific lung blocks: A feasibility study

Introduction

Total body irradiation (TBI) is an important component of the conditioning regimen in patients undergoing hematopoietic stem cell transplants. TBI is used in very few patients and therefore it is generally delivered with standard linear accelerators (LINACs) and not with dedicated devices. Severe pulmonary toxicity is the most common adverse effect after TBI, and patient-specific lead blocks are used to reduce mean lung dose. In this context, online treatment setup is crucial to achieve precise positioning of the lung blocks. Therefore, in this study we aim to report our experience at generating 3D-printed patient-specific lung blocks and coupling a dedicated couch (with an integrated onboard image device) with a modern LINAC for TBI treatment.

Material and methods

TBI was planned and delivered (2Gy/fraction given twice a day, over 3 days) to 15 patients. Online images, to be compared with planned digitally reconstructed radiographies, were acquired with the couch-dedicated Electronic Portal Imaging Device (EPID) panel and imported in the iView software using a homemade Graphical User Interface (GUI). In vivo dosimetry, using Metal-Oxide Field-Effect Transistors (MOSFETs), was used to assess the setup reproducibility in both supine and prone positions.

Results

3D printing of lung blocks was feasible for all planned patients using a stereolithography 3D printer with a build volume of 14.5×14.5×17.5 cm3. The number of required pre-TBI EPID-images generally decreases after the first fraction. In patient-specific quality assurance, the difference between measured and calculated dose was generally<2%. The MOSFET measurements reproducibility along each treatment and patient was 2.7%, in average.

Conclusion

The TBI technique was successfully implemented, demonstrating that our approach is feasible, flexible, and cost-effective. The use of 3D-printed patient-specific lung blocks have the potential to personalize TBI treatment and to refine the shape of the blocks before delivery, making them extremely versatile.

Machine-learning-based prediction of the effectiveness of the delivered dose by exhale-gated radiotherapy for locally advanced lung cancer: The additional value of geometric over dosimetric parameters alone

Purpose

This study aimed to assess interfraction stability of the delivered dose distribution by exhale-gated volumetric modulated arc therapy (VMAT) or intensity-modulated arc therapy (IMAT) for lung cancer and to determine dominant prognostic dosimetric and geometric factors.

Methods

Clinical target volume (CTV_Plan) from the planning CT was deformed to the exhale-gated daily CBCT scans to determine CTV_i, treated by the respective dose fraction. The equivalent uniform dose of the CTV_i was determined by the power law (gEUD_i) and cell survival model (EUD_iSF) as effectiveness measure for the delivered dose distribution. The following prognostic factors were analyzed: (I) minimum dose within the CTV_i (D_{min_i}), (II) Hausdorff distance (HDD_i) between CTV_i and CTV_Plan, (III) doses and deformations at the point in CTV_Plan at which the global minimum dose over all fractions per patient occurs (PD_{min_global_i}), and (IV) deformations at the point over all CTV_i margins per patient with the largest Hausdorff distance (HDPw_orst). Prognostic value and generalizability of the prognostic factors were examined using cross-validated random forest or multilayer perceptron neural network (MLP) classifiers. Dose accumulation was performed using back deformation of the dose distribution from CTV_i to CTV_Plan.

Results

Altogether, 218 dose fractions (10 patients) were evaluated. There was a significant interpatient heterogeneity between the distributions of the normalized gEUD_i values (p<0.0001, Kruskal-Wallis tests). Accumulated gEUD over all fractions per patient was 1.004-1.023 times of the prescribed dose. Accumulation led to tolerance of ~20% of fractions with gEUD_i <93% of the prescribed dose. Normalized D_min >60% was associated with predicted gEUD values above 95%. D_min had the highest importance for predicting the gEUD over all analyzed prognostic parameters by out-of-bag loss reduction using the random forest procedure. Cross-validated random forest classifier based on D_min as the sole input had the largest Pearson correlation coefficient (R=0.897) in comparison to classifiers using additional input variables. The neural network performed better than the random forest classifier, and the gEUD values predicted by the MLP classifier with D_min as the sole input were correlated with the gEUD values characterized by R=0.933 (95% CI, 0.913-0.948). The performance of the full MLP model with all geometric input parameters was slightly better (R=0.952) than that based on D_min (p=0.0034, Z-test).

Conclusion

Accumulated dose distributions over the treatment series were robust against interfraction CTV deformations using exhale gating and online image guidance. D_min was the most important parameter for gEUD prediction for a single fraction. All other parameters did not lead to a markedly improved generalizable prediction. Dosimetric information, especially location and value of D_min within the CTV _i , are vital information for image-guided radiation treatment.

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.

Optimization of mesh generation for geometric accuracy, robustness, and efficiency of biomechanical-model-based deformable image registration

BACKGROUND

Successful generation of biomechanical-model-based deformable image registration (BM-DIR) relies on user-defined parameters that dictate surface mesh quality. The trial-and-error process to determine the optimal parameters can be labor-intensive and hinder DIR efficiency and clinical workflow.

PURPOSE

To identify optimal parameters in surface mesh generation as boundary conditions for a BM-DIR in longitudinal liver and lung CT images to facilitate streamlined image registration processes.

METHODS

Contrast-enhanced CT images of 29 colorectal liver cancer patients and end-exhale four-dimensional CT images of 26 locally advanced non-small cell lung cancer patients were collected. Different combinations of parameters that determine the triangle mesh quality (voxel side length and triangle edge length) were investigated. The quality of DIRs generated using these parameters was evaluated with metrics for geometric accuracy, robustness, and efficiency. Metrics for geometric accuracy included target registration error (TRE) of internal vessel bifurcations, dice similar coefficient (DSC), mean distance to agreement (MDA), Hausdorff distance (HD) for organ contours, and number of vertices in the triangle mesh. American Association of Physicists in Medicine Task Group 132 was used to ensure parameters met TRE, DSC, MDA recommendations before the comparison among the parameters. Robustness was evaluated as the success rate of DIR generation, and efficiency was evaluated as the total time to generate boundary conditions and compute finite element analysis.

RESULTS

Voxel side length of 0.2 cm and triangle edge length of 3 were found to be the optimal parameters for both liver and lung, with success rate of 1.00 and 0.98 and average DIR computation time of 100 and 143 s, respectively. For this combination, the average TRE, DSC, MDA, and HD were 0.38-0.40, 0.96-0.97, 0.09-0.12, and 0.87-1.17 mm, respectively.

CONCLUSION

The optimal parameters were found for the analyzed patients. The decision-making process described in this study serves as a recommendation for BM-DIR algorithms to be used for liver and lung. These parameters can facilitate consistence in the evaluation of published studies and more widespread utilization of BM-DIR in clinical practice.

Accuracy and Feasibility of Synthetic CT for Lung Adaptive Radiotherapy: A Phantom Study

OBJECTIVES

The respiratory variations will lead to inconsistency between the actual delivery dose and the planning dose. How the minor interfractional amplitude changes affect the geometry and dose delivery accuracy remains to be investigated in the context of lung adaptive radiotherapy.

METHODS

Planning 4-dimensional-computed tomography and kV-cone beam computed tomography were scanned based on the Computerized Imaging Reference Systems phantom, which was employed to simulate the minor interfractional amplitude variations. The corresponding synthetic computed tomography for a particular motion pattern can be generated from Velocity program. Then a clinically meaningful synthetic computed tomography was analyzed through the geometrical and dosimetric assessment.

RESULTS

The image quality of synthetic computed tomography was improved obviously compared with cone beam computed tomography. Mean absolute error was minimized when no significant interfractional motion occurs and Velocity can be qualified for dealing with the regular breathing motion patterns. The mean percent hounsfield unit difference of the synthetic hounsfield unit values per organ relative to the planning 4-dimensional-computed tomography image was 22.3%. Under the same conditions, the mean percent hounsfield unit difference of the cone beam computed tomography hounsfield unit values per organ, relative to the planning 4-dimensional-computed tomography image was 83.9%. Overall, the accuracy of hounsfield unit in synthetic computed tomography was improved obviously and the variability of the synthetic image correlates with the planning 4-dimensional-computed tomography image variability. Meanwhile, the dose-volume histograms between planning 4-dimensional-computed tomography and synthetic computed tomography almost coincided each other, which indicates that Velocity program can qualify lung adaptive radiotherapy well when there were no interfractional respiratory variations. However, for cases with obvious interfractional amplitude change, the volume covered at least by 100% of the prescription dose was only 59.6% for that synthetic image.

CONCLUSION

The synthetic computed tomography images generated from Velocity were close to the real images in anatomy and dosimetry, which can make clinical lung adaptive radiotherapy possible based on the actual patient anatomy during treatment.

Prospective imaging comparison of anatomic delineation with rapid kV cone beam CT on a novel ring gantry radiotherapy device

INTRODUCTION

A kV imager coupled to a novel, ring-gantry radiotherapy system offers improved on-board kV-cone-beam computed tomography (CBCT) acquisition time (17-40 seconds) and image quality, which may improve CT radiotherapy image-guidance and enable online adaptive radiotherapy. We evaluated whether inter-observer contour variability over various anatomic structures was non-inferior using a novel ring gantry kV-CBCT (RG-CBCT) imager as compared to diagnostic-quality simulation CT (simCT).

MATERIALS/METHODS

Seven patients undergoing radiotherapy were imaged with the RG-CBCT system at breath hold (BH) and/or free breathing (FB) for various disease sites on a prospective imaging study. Anatomy was independently contoured by seven radiation oncologists on: 1. SimCT 2. Standard C-arm kV-CBCT (CA-CBCT), and 3. Novel RG-CBCT at FB and BH. Inter-observer contour variability was evaluated by computing simultaneous truth and performance level estimation (STAPLE) consensus contours, then computing average symmetric surface distance (ASSD) and Dice similarity coefficient (DSC) between individual raters and consensus contours for comparison across image types.

RESULTS

Across 7 patients, 18 organs-at-risk (OARs) were evaluated on 27 image sets. Both BH and FB RG-CBCT were non-inferior to simCT for inter-observer delineation variability across all OARs and patients by ASSD analysis (p < 0.001), whereas CA-CBCT was not (p = 0.923). RG-CBCT (FB and BH) also remained non-inferior for abdomen and breast subsites compared to simCT on ASSD analysis (p < 0.025). On DSC comparison, neither RG-CBCT nor CA-CBCT were non-inferior to simCT for all sites (p > 0.025).

CONCLUSIONS

Inter-observer ability to delineate OARs using novel RG-CBCT images was non-inferior to simCT by the ASSD criterion but not DSC criterion.

Prediction of adaptive strategies based on deformation vector field features for MR-guided adaptive radiotherapy of prostate cancer

BACKGROUND

The selection of adaptive strategies in MR-guided adaptive radiotherapy (MRgART) usually relies on subjective review of anatomical changes. However, this kind of review may lead to improper selection of adaptive strategy for some fractions.

PURPOSE

The purpose of this study was to develop prediction models based on deformation vector field (DVF) features for automatic and accurate strategy selection, using prostate cancer as an example.

METHODS

100 fractions of 20 prostate cancer patients were retrospectively selected in this study. Treatment plans using both adapt to position (ATP) strategy and adapt to shape (ATS) strategy were generated. Optimal adaptive strategy was determined according to dosimetric evaluation. DVFs of the deformable image registration (DIR) of daily MRI and CT simulation scans were extracted. The shape, first order statistics, and spatial features were extracted from the DVFs, subjected to further selection using the minimum redundancy maximum relevance (mRMR) method. The number of features (F_n ) was hyper-tuned using bootstrapping method, and then F_n indicating a peak area under the curve (AUC) value was used to construct three prediction models.

RESULTS

According to subjective review, the ATS strategy was adopted for all 100 fractions. However, the evaluation results showed that the ATP strategy could have met the clinical requirements for 23 (23%) fractions. The three prediction models showed high prediction performance, with the best performing model achieving an AUC value of 0.896, corresponding accuracy (ACC), sensitivity (SEN) and specificity (SPC) of 0.9, 0.958, and 0.667, respectively. The features used to construct prediction models included four features extracted from y direction of DVF (DVF_y ) and mask, one feature from z direction of DVF (DVF_z ). It indicated that the deformation along the anterior-posterior direction had a greater impact on determining the adaptive strategy than other directions.

CONCLUSIONS

DVF-feature-based models could accurately predict the adaptive strategy and avoid unnecessary selection of time-consuming ATS strategy, which consequently improves the efficiency of the MRgART process.

Interfractional dose accumulation for MR-guided liver SBRT: Variation among algorithms is highly patient- and fraction-dependent

BACKGROUND AND PURPOSE

Daily plan adaptations could take the dose delivered in previous fractions into account. Due to high dose delivered per fraction, low number of fractions, steep dose gradients, and large interfractional organ deformations, this might be particularly important for liver SBRT. This study investigates inter-algorithm variation of interfractional dose accumulation for MR-guided liver SBRT.

MATERIALS AND METHODS

We assessed 27 consecutive MR-guided liver SBRT treatments of 67.5 Gy in three (n = 15) or 50 Gy in five fractions (n = 12), both prescribed to the GTV. We calculated fraction doses on daily patient anatomy, warped these doses to the simulation MRI using seven different algorithms, and accumulated the warped doses. Thus, we obtained differences in planned doses and warped or accumulated doses for each algorithm. This enabled us to calculate the inter-algorithm variations in warped doses per fraction and in accumulated doses per treatment course.

RESULTS

The four intensity-based algorithms were more consistent with planned PTV dose than affine or contour-based algorithms. The mean (range) variation of the dose difference for PTV D_95% due to dose warping by these intensity-based algorithms was 10.4 percentage points (0.3 to 43.7) between fractions and 8.6 (0.3 to 24.9) between accumulated treatment doses. As seen by these ranges, the variation was very dependent on the patient and the fraction being analyzed. Nevertheless, no correlations between patient or plan characteristics on the one hand and inter-algorithm dose warping variation on the other hand was found.

CONCLUSION

Inter-algorithm dose accumulation variation is highly patient- and fraction-dependent for MR-guided liver SBRT. We advise against trusting a single algorithm for dose accumulation in liver SBRT.

Detailed dosimetric evaluation of inter-fraction and respiratory motion in lung stereotactic body radiation therapy based on daily 4D cone beam CT images

Objective. Periodic respiratory motion and inter-fraction variations are sources of geometric uncertainty in stereotactic body radiation therapy (SBRT) of pulmonary lesions. This study extensively evaluates and validates the separate and combined dosimetric effect of both factors using 4D-CT and daily 4D-cone beam CT (CBCT) dose accumulation scenarios.Approach. A first cohort of twenty early stage or metastatic disease lung cancer patients were retrospectively selected to evaluate each scenario. The planned-dose (3D_Ref) was optimized on a 3D mid-position CT. To estimate the dosimetric impact of respiratory motion (4D_Ref), inter-fractional variations (3D_Acc) and the combined effect of both factors (4D_Acc), three dose accumulation scenarios based on 4D-CT, daily mid-cone beam CT (CBCT) position and 4D-CBCT were implemented via CT-CT/CT-CBCT deformable image registration (DIR) techniques. Each scenario was compared to 3D_Ref.A separate cohort of ten lung SBRT patients was selected to validate DIR techniques. The distance discordance metric (DDM) was implemented per voxel and per patient for tumor and organs at risk (OARs), and the dosimetric impact for CT-CBCT DIR geometric errors was calculated.Main results.Median and interquartile range (IQR) of the dose difference per voxel were 0.05/2.69 Gy and -0.12/2.68 Gy for3DAcc-3DRefand4DAcc-3DRef.For4DRef-3DRefthe IQR was considerably smaller -0.15/0.78 Gy. These findings were confirmed by dose volume histogram parameters calculated in tumor and OARs. For CT-CT/CT-CBCT DIR validation, DDM (95th percentile) was highest for heart (6.26 mm)/spinal cord (8.00 mm), and below 3 mm for tumor and the rest of OARs. The dosimetric impact of CT-CBCT DIR errors was below 2 Gy for tumor and OARs.Significance. The dosimetric impact of inter-fraction variations were shown to dominate those of periodic respiration in SBRT for pulmonary lesions. Therefore, treatment evaluation and dose-effect studies would benefit more from dose accumulation focusing on day-to-day changes then those that focus on respiratory motion.

SWFT-Net: a deep learning framework for efficient fine-tuning spot weights towards adaptive proton therapy

Objective. One critical task for adaptive proton therapy is how to perform spot weight re-tuning and reoptimize plan, both of which are time-consuming and labor intensive. We proposed a deep learning framework (SWFT-Net) to speed up such a task, a starting point for us to move towards online adaptive proton therapy.Approach. For a H&N patient case, a reference intensity modulated proton therapy plan was generated. For data augmentation, spot weights were modified to generate three datasets (DS10, DS30, DS50), corresponding to different levels of weight adjustment. For each dataset, the samples were split into the training and testing groups at a ratio of 8:2 (6400 for training, 1706 for testing). To ease the difficulty of machine learning, the residuals of dose maps and spot weights (i.e. difference relative to a reference) were used as inputs and outputs, respectively. Quantitative analyses were performed in terms of normalized root mean square error (NRMSE) of spot weights, Gamma passing rate and dose difference within the PTV.Main results. The SWFT-Net is able to generate an adapted plan in less than a second with a NVIDIA GeForce RTX 3090 GPU. For the 1706 samples in the testing dataset, the NRMSE is 0.41% (DS10), 1.05% (DS30) and 2.04% (DS50), respectively. Cold/hot spots in the dose maps after adaptation are observed. The mean relative dose difference is 0.64% (DS10), 0.92% (DS30) and 0.88% (DS50), respectively. For all three datasets, the mean Gamma passing rate is consistently over 95% for both 1 mm/1% and 3 mm/3% settings.Significance. The proposed SWFT-Net is a promising tool to help realize adaptive proton therapy. It can be used as an alternative tool to other spot fine-tuning optimization algorithms, likely demonstrating superior performance in terms of speed, accuracy, robustness and minimum human interaction. This study lays down a foundation for us to move further incorporating other factors such as daily anatomical changes and propagated PTVs, and develop a truly online adaptive workflow in proton therapy.

A review of deep learning-based deformable medical image registration

The alignment of images through deformable image registration is vital to clinical applications (e.g., atlas creation, image fusion, and tumor targeting in image-guided navigation systems) and is still a challenging problem. Recent progress in the field of deep learning has significantly advanced the performance of medical image registration. In this review, we present a comprehensive survey on deep learning-based deformable medical image registration methods. These methods are classified into five categories: Deep Iterative Methods, Supervised Methods, Unsupervised Methods, Weakly Supervised Methods, and Latest Methods. A detailed review of each category is provided with discussions about contributions, tasks, and inadequacies. We also provide statistical analysis for the selected papers from the point of view of image modality, the region of interest (ROI), evaluation metrics, and method categories. In addition, we summarize 33 publicly available datasets that are used for benchmarking the registration algorithms. Finally, the remaining challenges, future directions, and potential trends are discussed in our review.

Phantom-based Quality Assurance of a Clinical Dose Accumulation Technique Used in an Online Adaptive Radiation Therapy Platform

Purpose

This study aimed to develop a routine quality assurance method for a dose accumulation technique provided by a radiation therapy platform for online treatment adaptation.

Methods and Materials

Two commonly used phantoms were selected for the dose accumulation QA: Electron density and anthropomorphic pelvis. On a computed tomography (CT) scan of the electron density phantom, 1 target (gross tumor volume [GTV]; insert at 6 o'clock), a subvolume within this target, and 7 organs at risk (OARs; other inserts) were contoured in the treatment planning system (TPS). Two adaptation sessions were performed in which the GTV was recontoured, first at 7 o'clock and then at 5 o'clock. The accumulated dose was exported from the TPS after delivery. Deformable vector fields were also exported to manually accumulate doses for comparison. For the pelvis phantom, synthetic Gaussian deformations were applied to the planning CT image to simulate organ changes. Two single-fraction adaptive plans were created based on the deformed planning CT and cone beam CT images acquired onboard the radiation therapy platform. A manual dose accumulation was performed after delivery using the exported deformable vector fields for comparison with the system-generated result.

Results

All plans were successfully delivered, and the accumulated dose was both manually calculated and derived from the TPS. For the electron density phantom, the average mean dose differences in the GTV, boost volume, and OARs 1 to 7 were 0.0%, -0.2%, 92.0%, 78.4%, 1.8%, 1.9%, 0.0%, 0.0%, and 2.3%, respectively, between the manually summed and platform-accumulated doses. The gamma passing rates for the 3-dimensional dose comparison between the manually generated and TPS-provided dose accumulations were >99% for both phantoms.

Conclusions

This study demonstrated agreement between manually obtained and TPS-generated accumulated doses in terms of both mean structure doses and local 3-dimensional dose distributions. Large disagreements were observed for OAR1 and OAR2 defined on the electron density phantom due to OARs having lower deformation priority over the target in addition to artificially large changes in position induced for these structures fraction-by-fraction. The tests applied in this study to a commercial platform provide a straightforward approach toward the development of routine quality assurance of dose accumulation in online adaptation.

Automatic breast lesion segmentation in phase preserved DCE-MRIs

We offer a framework for automatically and accurately segmenting breast lesions from Dynamic Contrast Enhanced (DCE) MRI in this paper. The framework is built using max flow and min cut problems in the continuous domain over phase preserved denoised images. Three stages are required to complete the proposed approach. First, post-contrast and pre-contrast images are subtracted, followed by image registrations that benefit to enhancing lesion areas. Second, a phase preserved denoising and pixel-wise adaptive Wiener filtering technique is used, followed by max flow and min cut problems in a continuous domain. A denoising mechanism clears the noise in the images by preserving useful and detailed features such as edges. Then, lesion detection is performed using continuous max flow. Finally, a morphological operation is used as a post-processing step to further delineate the obtained results. A series of qualitative and quantitative trials employing nine performance metrics on 21 cases with two different MR image resolutions were used to verify the effectiveness of the proposed method. Performance results demonstrate the quality of segmentation obtained from the proposed method.

Image-to-Graph Convolutional Network for 2D/3D Deformable Model Registration of Low-Contrast Organs

Organ shape reconstruction based on a single-projection image during treatment has wide clinical scope, e.g., in image-guided radiotherapy and surgical guidance. We propose an image-to-graph convolutional network that achieves deformable registration of a three-dimensional (3D) organ mesh for a low-contrast two-dimensional (2D) projection image. This framework enables simultaneous training of two types of transformation: from the 2D projection image to a displacement map, and from the sampled per-vertex feature to a 3D displacement that satisfies the geometrical constraint of the mesh structure. Assuming application to radiation therapy, the 2D/3D deformable registration performance is verified for multiple abdominal organs that have not been targeted to date, i.e., the liver, stomach, duodenum, and kidney, and for pancreatic cancer. The experimental results show shape prediction considering relationships among multiple organs can be used to predict respiratory motion and deformation from digitally reconstructed radiographs with clinically acceptable accuracy.

A clinically relevant online patient QA solution with daily CT scans and EPID-based in vivo dosimetry: a feasibility study on rectal cancer

Objective. Adaptive radiation therapy (ART) could protect organs at risk (OARs) while maintain high dose coverage to targets. However, there is still a lack of efficient online patient quality assurance (QA) methods, which is an obstacle to large-scale adoption of ART. We aim to develop a clinically relevant online patient QA solution for ART using daily CT scans and EPID-based in vivo dosimetry. Approach. Ten patients with rectal cancer at our center were included. Patients’ daily CT scans and portal images were collected to generate reconstructed 3D dose distributions. Contours of targets and OARs were recontoured on these daily CT scans by a clinician or an auto-segmentation algorithm, then dose-volume indices were calculated, and the percent deviation of these indices to their original plans were determined. This deviation was regarded as the metric for clinically relevant patient QA. The tolerance level was obtained using a 95% confidence interval of the QA metric distribution. These deviations could be further divided into anatomically relevant or delivery relevant indicators for error source analysis. Finally, our QA solution was validated on an additional six clinical patients. Main results. In rectal cancer, the 95% confidence intervals of the QA metric for PTV ΔD 95 (%) were [−3.11%, 2.35%], and for PTV ΔD 2 (%) were [−0.78%, 3.23%]. In validation, 68% for PTV ΔD 95 (%), and 79% for PTV ΔD 2 (%) of the 28 fractions are within tolerances of the QA metrics. one patient’s dosimetric impact of anatomical variations during treatment were observed through the source of error analysis. Significance. The online patient QA solution using daily CT scans and EPID-based in vivo dosimetry is clinically feasible. Source of error analysis has the potential for distinguishing sources of error and guiding ART for future treatments.

Integrating external beam and prostate seed implant dosimetry for intermediate and high-risk prostate cancer using biologically effective dose: Impact of image registration technique

PURPOSE

Combining external beam radiation therapy (EBRT) and prostate seed implant (PSI) is efficacious in treating intermediate- and high-risk prostate cancer at the cost of increased genitourinary toxicity. Accurate combined dosimetry remains elusive due to lack of registration between treatment plans and different biological effect. The current work proposes a method to convert physical dose to biological effective dose (BED) and spatially register the dose distributions for more accurate combined dosimetry.

METHODS AND MATERIALS

A PSI phantom was CT scanned with and without seeds under rigid and deformed transformations. The resulting CTs were registered using image-based rigid registration (RI), fiducial-based rigid registration (RF), or b-spline deformable image registration (DIR) to determine which was most accurate. Physical EBRT and PSI dose distributions from a sample of 91 previously-treated combined-modality prostate cancer patients were converted to BED and registered using RI, RF, and DIR. Forty-eight (48) previously-treated patients whose PSI occurred before EBRT were included as a "control" group due to inherent registration. Dose-volume histogram (DVH) parameters were compared for RI, RF, DIR, DICOM, and scalar addition of DVH parameters using ANOVA or independent Student's t tests (α = 0.05).

RESULTS

In the phantom study, DIR was the most accurate registration algorithm, especially in the case of deformation. In the patient study, dosimetry from RI was significantly different than the other registration algorithms, including the control group. Dosimetry from RF and DIR were not significantly different from the control group or each other.

CONCLUSIONS

Combined dosimetry with BED and image registration is feasible. Future work will utilize this method to correlate dosimetry with clinical outcomes.

[Evaluation of Accuracy of Deformable Image Registration with Newly Developed Treatment Planning Support Software for Thoracic Images]

This study evaluated accuracy of deformable image registration (DIR) with twelve parameter settings for thoracic images. We used peak-inhale and peak-exhale images for ten patients provided by DIR-lab. We used a prototype version of iCView software (ITEM Corporation) with DIR to perform intensity, structure, and hybrid-based DIR with the twelve parameter settings. DIR accuracy was evaluated by a target registration error (TRE) using 300 bronchial bifurcations and the Dice similarity coefficient (DSC) of the lungs. For twelve parameter settings, TRE ranged from 2.83 mm to 5.27 mm, whereas DSC ranged from 0.96 to 0.98. These results demonstrated that DIR accuracy differed among parameter settings and show that appropriate parameter settings are required for clinical practice.

The application of multiple metrics in deformable image registration for target volume delineation of breast tumor bed

BACKGROUND AND PURPOSE

For postoperative breast cancer patients, deformable image registration (DIR) is challenged due to the large deformations and non-correspondence caused by tumor resection and clip insertion. To deal with it, three metrics (fiducial-, region-, and intensity-based) were jointly used in DIR algorithm for improved accuracy.

MATERIALS AND METHODS

Three types of metrics were combined to form a single-objective function in DIR algorithm. Fiducial-based metric was used to minimize the distance between the corresponding point sets of two images. Region-based metric was used to improve the overlap between the corresponding areas of two images. Intensity-based metric was used to maximize the correlation between the corresponding voxel intensities of two images. The two CT images, one before surgery and the other after surgery, were acquired from the same patient in the same radiotherapy treatment position. Twenty patients who underwent breast-conserving surgery and postoperative radiotherapy were enrolled in this study.

RESULTS

For target registration error, the difference between the proposed and the conventional registration methods was statistically significant for soft tissue (2.06 vs. 7.82, p = 0.00024 < 0.05) and body boundary (3.70 vs. 6.93, p = 0.021 < 0.05). For visual assessment, the proposed method achieved better matching result for soft tissue and body boundary.

CONCLUSIONS

Comparing to the conventional method, the registration accuracy of the proposed method was significantly improved. This method provided a feasible way for target volume delineation of tumor bed in postoperative radiotherapy of breast cancer patients.

Feasibility of 4D VMAT-CT

Objective.Feasibility of three-dimensional (3D) tracking of volumetric modulated arc therapy (VMAT) based on VMAT-computed tomography (VMAT-CT) has been shown previously by our group. However, 3D VMAT-CT is not suitable for treatments that involve significant target movement due to patient breathing. The goal of this study was to reconstruct four-dimensional (4D) VMAT-CT and evaluate the feasibility of tracking based on 4D VMAT-CT.Approach.Synchronized portal images of phantoms and linac log were both sorted into four phases, and VMAT-CT+ was generated in each phase by fusing reconstructed VMAT-CT and planning CT using rigid or deformable registration. Dose was calculated in each phase and was registered to the mean position planning CT for 4D dose reconstruction. Trackings based on 4D VMAT-CT+ and 4D cone beam CT (CBCT) were compared. Potential uncertainties were also evaluated.Main results.Tracking based on 4D VMAT-CT+ was accurate, could detect phantom deformation and/or change of breathing pattern, and was superior to that based on 4D CBCT. The impact of uncertainties on tracking was minimal.Significance.Our study shows it is feasible to accurately track position and dose based on 4D VMAT-CT for patients whose VMAT treatments are subject to respiratory motion. It will significantly increase the confidence of VMAT and is a clinically viable solution to daily patient positioning,in vivodosimetry and treatment monitoring.

COMP Report: CPQR technical quality control guidelines for use of positron emission tomography/computed tomography in radiation treatment planning

Positron emission tomography with x-ray computed tomography (PET/CT) is increasingly being utilized for radiation treatment planning (RTP). Accurate delivery of RT therefore depends on quality PET/CT data. This study covers quality control (QC) procedures required for PET/CT for diagnostic imaging and incremental QC required for RTP. Based on a review of the literature, it compiles a list of recommended tests, performance frequencies, and tolerances, as well as references to documents detailing how to perform each test. The report was commissioned by the Canadian Organization of Medical Physicists as part of the Canadian Partnership for Quality Radiotherapy initiative.

A Prior Knowledge-Guided, Deep Learning-Based Semiautomatic Segmentation for Complex Anatomy on Magnetic Resonance Imaging

PURPOSE

Despite recent substantial improvement in autosegmentation using deep learning (DL) methods, labor-intensive and time-consuming slice-by-slice manual editing is often needed, particularly for complex anatomy (eg, abdominal organs). This work aimed to develop a fast, prior knowledge-guided DL semiautomatic segmentation (DL-SAS) method for complex structures on abdominal magnetic resonance imaging (MRI) scans.

METHODS AND MATERIALS

A novel application using contours on an adjacent slice as a prior knowledge informant in a 2-dimensional UNet DL model to guide autosegmentation for a subsequent slice was implemented for DL-SAS. A generalized, instead of organ-specific, DL-SAS model was trained and tested for abdominal organs on T2-weighted MRI scans collected from 75 patients (65 for training and 10 for testing). The DL-SAS model performance was compared with 3 common autocontouring methods (linear interpolation, rigid propagation, and a full 3-dimensional DL autosegmentation model trained with the same training data set) based on various quantitative metrics including the Dice similarity coefficient (DSC) and ratio of acceptable slices (ROA) using paired t tests.

RESULTS

For the 10 testing cases, the DL-SAS model performed best with the slice interval (SI) of 1, resulting in an average DSC of 0.93 ± 0.02, 0.92 ± 0.02, 0.91 ± 0.02, 0.88 ± 0.03, and 0.87 ± 0.02 for the large bowel, stomach, small bowel, duodenum, and pancreas, respectively. The performance decreased with increased SIs from the guidance slice. The DL-SAS method performed significantly better (P < .05) than the other 3 methods. The ROA values were in the range of 48% to 66% for all the organs with an SI of 1 for DL-SAS, higher than those for linear interpolation (31%-57% for an SI of 1) and DL auto-segmentation (16%-51%).

CONCLUSIONS

The developed DL-SAS model segmented complex abdominal structures on MRI with high accuracy and efficiency and may be implemented as an interactive manual contouring tool or a contour editing tool in conjunction with a full autosegmentation process, facilitating fast and accurate segmentation for MRI-guided online adaptive radiation therapy.

Estimating medical image registration error and confidence: A taxonomy and scoping review

Given that image registration is a fundamental and ubiquitous task in both clinical and research domains of the medical field, errors in registration can have serious consequences. Since such errors can mislead clinicians during image-guided therapies or bias the results of a downstream analysis, methods to estimate registration error are becoming more popular. To give structure to this new heterogenous field we developed a taxonomy and performed a scoping review of methods that quantitatively and automatically provide a dense estimation of registration error. The taxonomy breaks down error estimation methods into Approach (Image- or Transformation-based), Framework (Machine Learning or Direct) and Measurement (error or confidence) components. Following the PRISMA guidelines for scoping reviews, the 570 records found were reduced to twenty studies that met inclusion criteria, which were then reviewed according to the proposed taxonomy. Trends in the field, advantages and disadvantages of the methods, and potential sources of bias are also discussed. We provide suggestions for best practices and identify areas of future research.

Evaluation of dose accumulation methods and workflows utilising cone beam computed tomography images

INTRODUCTION

Various adaptive radiation therapy (ART) methods have emerged, with little consensus amongst the literature as to which is most appropriate. This study aimed to compare dose mapping (DM) versus Monte Carlo recalculation (MCR), using cone beam computed tomography (CBCT) images when utilised in automated ART dose accumulation workflows in the MIM Maestro software package.

METHODS

The treatment plans for 38 cancer patients (19 prostate and 19 head and neck cases) were used to perform DM or MCR retrospectively upon CBCTs acquired during treatment, which were then deformably registered to the planning CT (DR-pCT) to facilitate dose accumulation. Dose-volume and region-of-interest data were extracted for the planning target volumes and organs at risk. Intraclass correlation (ICC) values and Bland-Altman plots were utilised to compare DM versus MCR doses on the CBCT images as well as CBCT versus DR-pCT doses.

RESULTS

When comparing DM and MCR on CBCTs, the differences across dose level mean dose differences were mostly within a ±5% level of agreement based on the Bland-Altman plots, with over 67% of ICC values over 0.9 and indicative of good correlation. When these distributions were deformed back to the planning CT, the agreement was reduced considerably, with larger differences (exceeding ±5%) resulting from workflow-related issues.

CONCLUSION

The results emphasise the need to consider and make adaptations to minimise the effect of workflows on algorithm performance. Manual user intervention, refined departmental protocols and further developments to the MIM Maestro software will enhance the use of this tool.

Treatment-integrated imaging, radiomics, and personalised radiotherapy: the future is at hand

Since the introduction of computed tomography for planning purposes in the 1970s, we have been observing a continuous development of different imaging methods in radiotherapy. The current achievements of imaging technologies in radiotherapy enable more than just improvement of accuracy on the planning stage. Through integrating imaging with treatment machines, they allow advanced control methods of dose delivery during the treatment. This article reviews how the integration of existing and novel forms of imaging changes radiotherapy and how these advances can allow a more individualised approach to cancer therapy. We believe that the significant challenge for the next decade is the continued integration of a range of different imaging devices into linear accelerators. These imaging modalities should show intra-fraction changes in body morphology and inter-fraction metabolic changes. As the use of these more advanced, integrated machines grows, radiotherapy delivery will become more accurate, thus resulting in better clinical outcomes: higher cure rates with fewer side effects.

The impact of an educational tool in cervix image registration across three imaging modalities

OBJECTIVES

Accurate image registration is vital in cervical cancer where changes in both planning target volume (PTV) and organs at risk (OARs) can make decisions regarding image registration complicated. This work aims to determine the impact of a dedicated educational tool compared with experience gained in MR-guided radiotherapy (MRgRT).

METHODS

10 therapeutic radiographers acted as observers and were split into two groups based on previous experience with MRgRT and Monaco treatment planning system. Three CBCT-CT, three MR-CT and two MR-MR registrations were completed per patient by each observer. Observers recorded translations, time to complete image registration and confidence. Data were collected in two phases; prior to and following the introduction of a cervix registration guide.

RESULTS

No statistically significant differences were noted between imaging modalities. Each group was assessed independently pre- and post-education, no statistically significant differences were noted in either CBCT-CT or MR-CT imaging. Group 1 MR-MR imaging showed a statistically significant reduction in interobserver variability (p=0.04), in Group 2, the result was not statistically significant (p=0.06). Statistically significant increases in confidence were seen in all three modalities (p≤0.05).

CONCLUSIONS

At The Christie NHS Foundation Trust, radiographers consistently registered images across three different imaging modalities regardless of their previous experience. The implementation of an image registration guide had limited impact on inter- and intraobserver variability. Radiographers' confidence showed statistically significant improvements following the use of the registration manual.

ADVANCES IN KNOWLEDGE

This work helps evaluate training methods for novel roles that are developing in MRgRT.

CDFRegNet: A cross-domain fusion registration network for CT-to-CBCT image registration

BACKGROUND AND OBJECTIVE

Computer tomography (CT) to cone-beam computed tomography (CBCT) image registration plays an important role in radiotherapy treatment placement, dose verification, and anatomic changes monitoring during radiotherapy. However, fast and accurate CT-to-CBCT image registration is still very challenging due to the intensity differences, the poor image quality of CBCT images, and inconsistent structure information.

METHODS

To address these problems, a novel unsupervised network named cross-domain fusion registration network (CDFRegNet) is proposed. First, a novel edge-guided attention module (EGAM) is designed, aiming at capturing edge information based on the gradient prior images and guiding the network to model the spatial correspondence between two image domains. Moreover, a novel cross-domain attention module (CDAM) is proposed to improve the network's ability to guide the network to effectively map and fuse the domain-specific features.

RESULTS

Extensive experiments on a real clinical dataset were carried out, and the experimental results verify that the proposed CDFRegNet can register CT to CBCT images effectively and obtain the best performance, while compared with other representative methods, with a mean DSC of 80.01±7.16%, a mean TRE of 2.27±0.62 mm, and a mean MHD of 1.50±0.32 mm. The ablation experiments also proved that our EGAM and CDAM can further improve the accuracy of the registration network and they can generalize well to other registration networks.

CONCLUSION

This paper proposed a novel CT-to-CBCT registration method based on EGAM and CDAM, which has the potential to improve the accuracy of multi-domain image registration.