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

Adaptive Radiation Therapy: A Review of CT-based Techniques

Adaptive radiation therapy is a feedback process by which imaging information acquired over the course of treatment, such as changes in patient anatomy, can be used to reoptimize the treatment plan, with the end goal of improving target coverage and reducing treatment toxicity. This review describes different types of adaptive radiation therapy and their clinical implementation with a focus on CT-guided online adaptive radiation therapy. Depending on local anatomic changes and clinical context, different anatomic sites and/or disease stages and presentations benefit from different adaptation strategies. Online adaptive radiation therapy, where images acquired in-room before each fraction are used to adjust the treatment plan while the patient remains on the treatment table, has emerged to address unpredictable anatomic changes between treatment fractions. Online treatment adaptation places unique pressures on the radiation therapy workflow, requiring high-quality daily imaging and rapid recontouring, replanning, plan review, and quality assurance. Generating a new plan with every fraction is resource intensive and time sensitive, emphasizing the need for workflow efficiency and clinical resource allocation. Cone-beam CT is widely used for image-guided radiation therapy, so implementing cone-beam CT-guided online adaptive radiation therapy can be easily integrated into the radiation therapy workflow and potentially allow for rapid imaging and replanning. The major challenge of this approach is the reduced image quality due to poor resolution, scatter, and artifacts. Keywords: Adaptive Radiation Therapy, Cone-Beam CT, Organs at Risk, Oncology © RSNA, 2023.

Improving accelerated 3D imaging in MRI-guided radiotherapy for prostate cancer using a deep learning method

PURPOSE

This study was to improve image quality for high-speed MR imaging using a deep learning method for online adaptive radiotherapy in prostate cancer. We then evaluated its benefits on image registration.

METHODS

Sixty pairs of 1.5 T MR images acquired with an MR-linac were enrolled. The data included low-speed, high-quality (LSHQ), and high-speed low-quality (HSLQ) MR images. We proposed a CycleGAN, which is based on the data augmentation technique, to learn the mapping between the HSLQ and LSHQ images and then generate synthetic LSHQ (synLSHQ) images from the HSLQ images. Five-fold cross-validation was employed to test the CycleGAN model. The normalized mean absolute error (nMAE), peak signal-to-noise ratio (PSNR), structural similarity index measurement (SSIM), and edge keeping index (EKI) were calculated to determine image quality. The Jacobian determinant value (JDV), Dice similarity coefficient (DSC), and mean distance to agreement (MDA) were used to analyze deformable registration.

RESULTS

Compared with the LSHQ, the proposed synLSHQ achieved comparable image quality and reduced imaging time by ~ 66%. Compared with the HSLQ, the synLSHQ had better image quality with improvement of 57%, 3.4%, 26.9%, and 3.6% for nMAE, SSIM, PSNR, and EKI, respectively. Furthermore, the synLSHQ enhanced registration accuracy with a superior mean JDV (6%) and preferable DSC and MDA values compared with HSLQ.

CONCLUSION

The proposed method can generate high-quality images from high-speed scanning sequences. As a result, it shows potential to shorten the scan time while ensuring the accuracy of radiotherapy.

Dose outside of the prostate is associated with improved outcomes for high-risk prostate cancer patients treated with brachytherapy boost

Background

One in three high-risk prostate cancer patients treated with radiotherapy recur. Detection of lymph node metastasis and microscopic disease spread using conventional imaging is poor, and many patients are under-treated due to suboptimal seminal vesicle or lymph node irradiation. We use Image Based Data Mining (IBDM) to investigate association between dose distributions, and prognostic variables and biochemical recurrence (BCR) in prostate cancer patients treated with radiotherapy. We further test whether including dose information in risk-stratification models improves performance.

Method

Planning CTs, dose distributions and clinical information were collected for 612 high-risk prostate cancer patients treated with conformal hypo-fractionated radiotherapy, intensity modulated radiotherapy (IMRT), or IMRT plus a single fraction high dose rate (HDR) brachytherapy boost. Dose distributions (including HDR boost) of all studied patients were mapped to a reference anatomy using the prostate delineations. Regions where dose distributions significantly differed between patients that did and did-not experience BCR were assessed voxel-wise using 1) a binary endpoint of BCR at four-years (dose only) and 2) Cox-IBDM (dose and prognostic variables). Regions where dose was associated with outcome were identified. Cox proportional-hazard models with and without region dose information were produced and the Akaike Information Criterion (AIC) was used to assess model performance.

Results

No significant regions were observed for patients treated with hypo-fractionated radiotherapy or IMRT. Regions outside the target where higher dose was associated with lower BCR were observed for patients treated with brachytherapy boost. Cox-IBDM revealed that dose response was influenced by age and T-stage. A region at the seminal vesicle tips was identified in binary- and Cox-IBDM. Including the mean dose in this region in a risk-stratification model (hazard ratio=0.84, p=0.005) significantly reduced AIC values (p=0.019), indicating superior performance, compared with prognostic variables only. The region dose was lower in the brachytherapy boost patients compared with the external beam cohorts supporting the occurrence of marginal misses.

Conclusion

Association was identified between BCR and dose outside of the target region in high-risk prostate cancer patients treated with IMRT plus brachytherapy boost. We show, for the first-time, that the importance of irradiating this region is linked to prognostic variables.

Comparison of different dose accumulation strategies to estimate organ doses after stereotactic magnetic resonance-guided adaptive radiotherapy

INTRODUCTION

Re-irradiation is frequently performed in the era of precision oncology, but previous doses to organs-at-risk (OAR) must be assessed to avoid cumulative overdoses. Stereotactic magnetic resonance-guided online adaptive radiotherapy (SMART) enables highly precise ablation of tumors close to OAR. However, OAR doses may change considerably during adaptive treatment, which complicates potential re-irradiation. We aimed to compare the baseline plan with different dose accumulation techniques to inform re-irradiation.

PATIENTS & METHODS

We analyzed 18 patients who received SMART to lung or liver tumors inside prospective databases. Cumulative doses were calculated inside the planning target volumes (PTV) and OAR for the adapted plans and theoretical non-adapted plans via (1) cumulative dose volume histograms (DVH sum plan) and (2) deformable image registration (DIR)-based dose accumulation to planning images (DIR sum plan). We compared cumulative dose parameters between the baseline plan, DVH sum plan and DIR sum plan using equivalent doses in 2 Gy fractions (EQD2).

RESULTS

Individual patients presented relevant increases of near-maximum doses inside the proximal bronchial tree, spinal cord, heart and gastrointestinal OAR when comparing adaptive treatment to the baseline plans. The spinal cord near-maximum doses were significantly increased in the liver patients (D2% median: baseline 6.1 Gy, DIR sum 8.1 Gy, DVH sum 8.4 Gy, p = 0.04; D0.1 cm³ median: baseline 6.1 Gy, DIR sum 8.1 Gy, DVH sum 8.5 Gy, p = 0.04). Three OAR overdoses occurred during adaptive treatment (DIR sum: 1, DVH sum: 2), and four more intense OAR overdoses would have occurred during non-adaptive treatment (DIR sum: 4, DVH sum: 3). Adaptive treatment maintained similar PTV coverages to the baseline plans, while non-adaptive treatment yielded significantly worse PTV coverages in the lung (D95% median: baseline 86.4 Gy, DIR sum 82.4 Gy, DVH sum 82.2 Gy, p = 0.006) and liver patients (D95% median: baseline 87.4 Gy, DIR sum 82.1 Gy, DVH sum 81.1 Gy, p = 0.04).

CONCLUSION

OAR doses can increase during SMART, so that re-irradiation should be planned based on dose accumulations of the adapted plans instead of the baseline plan. Cumulative dose volume histograms represent a simple and conservative dose accumulation strategy.

Optimizing 3DCT image registration for interfractional changes in carbon-ion prostate radiotherapy

To perform setup procedures including both positional and dosimetric information, we developed a CT-CT rigid image registration algorithm utilizing water equivalent pathlength (WEPL)-based image registration and compared the resulting dose distribution with those of two other algorithms, intensity-based image registration and target-based image registration, in prostate cancer radiotherapy using the carbon-ion pencil beam scanning technique. We used the data of the carbon ion therapy planning CT and the four-weekly treatment CTs of 19 prostate cancer cases. Three CT-CT registration algorithms were used to register the treatment CTs to the planning CT. Intensity-based image registration uses CT voxel intensity information. Target-based image registration uses target position on the treatment CTs to register it to that on the planning CT. WEPL-based image registration registers the treatment CTs to the planning CT using WEPL values. Initial dose distributions were calculated using the planning CT with the lateral beam angles. The treatment plan parameters were optimized to administer the prescribed dose to the PTV on the planning CT. Weekly dose distributions using the three different algorithms were calculated by applying the treatment plan parameters to the weekly CT data. Dosimetry, including the dose received by 95% of the clinical target volume (CTV-D95), rectal volumes receiving > 20 Gy (RBE) (V20), > 30 Gy (RBE) (V30), and > 40 Gy (RBE) (V40), were calculated. Statistical significance was assessed using the Wilcoxon signed-rank test. Interfractional CTV displacement over all patients was 6.0 ± 2.7 mm (19.3 mm maximum standard amount). WEPL differences between the planning CT and the treatment CT were 1.2 ± 0.6 mm-H₂O (< 3.9 mm-H₂O), 1.7 ± 0.9 mm-H₂O (< 5.7 mm-H₂O) and 1.5 ± 0.7 mm-H₂O (< 3.6 mm-H₂O maxima) with the intensity-based image registration, target-based image registration, and WEPL-based image registration, respectively. For CTV coverage, the D95 values on the planning CT were > 95% of the prescribed dose in all cases. The mean CTV-D95 values were 95.8 ± 11.5% and 98.8 ± 1.7% with the intensity-based image registration and target-based image registration, respectively. The WEPL-based image registration was CTV-D95 to 99.0 ± 0.4% and rectal Dmax to 51.9 ± 1.9 Gy (RBE) compared to 49.4 ± 9.1 Gy (RBE) with intensity-based image registration and 52.2 ± 1.8 Gy (RBE) with target-based image registration. The WEPL-based image registration algorithm improved the target coverage from the other algorithms and reduced rectal dose from the target-based image registration, even though the magnitude of the interfractional variation was increased.

Intra-Patient Lung CT Registration through Large Deformation Decomposition and Attention-Guided Refinement

Deformable lung CT image registration is an essential task for computer-assisted interventions and other clinical applications, especially when organ motion is involved. While deep-learning-based image registration methods have recently achieved promising results by inferring deformation fields in an end-to-end manner, large and irregular deformations caused by organ motion still pose a significant challenge. In this paper, we present a method for registering lung CT images that is tailored to the specific patient being imaged. To address the challenge of large deformations between the source and target images, we break the deformation down into multiple continuous intermediate fields. These fields are then combined to create a spatio-temporal motion field. We further refine this field using a self-attention layer that aggregates information along motion trajectories. By leveraging temporal information from a respiratory cycle, our proposed methods can generate intermediate images that facilitate image-guided tumor tracking. We evaluated our approach extensively on a public dataset, and our numerical and visual results demonstrate the effectiveness of the proposed method.

Applicability and usage of dose mapping/accumulation in radiotherapy

Dose mapping/accumulation (DMA) is a topic in radiotherapy (RT) for years, but has not yet found its widespread way into clinical RT routine. During the ESTRO Physics workshop 2021 on "commissioning and quality assurance of deformable image registration (DIR) for current and future RT applications", we built a working group on DMA from which we present the results of our discussions in this article. Our aim in this manuscript is to shed light on the current situation of DMA in RT and to highlight the issues that hinder consciously integrating it into clinical RT routine. As a first outcome of our discussions, we present a scheme where representative RT use cases are positioned, considering expected anatomical variations and the impact of dose mapping uncertainties on patient safety, which we have named the DMA landscape (DMAL). This tool is useful for future reference when DMA applications get closer to clinical day-to-day use. Secondly, we discussed current challenges, lightly touching on first-order effects (related to the impact of DIR uncertainties in dose mapping), and focusing in detail on second-order effects often dismissed in the current literature (as resampling and interpolation, quality assurance considerations, and radiobiological issues). Finally, we developed recommendations, and guidelines for vendors and users. Our main point include: Strive for context-driven DIR (by considering their impact on clinical decisions/judgements) rather than perfect DIR; be conscious of the limitations of the implemented DIR algorithm; and consider when dose mapping (with properly quantified uncertainties) is a better alternative than no mapping.

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.

CBCT-guided adaptive radiotherapy using self-supervised sequential domain adaptation with uncertainty estimation

Adaptive radiotherapy (ART) is an advanced technology in modern cancer treatment that incorporates progressive changes in patient anatomy into active plan/dose adaption during the fractionated treatment. However, the clinical application relies on the accurate segmentation of cancer tumors on low-quality on-board images, which has posed challenges for both manual delineation and deep learning-based models. In this paper, we propose a novel sequence transduction deep neural network with an attention mechanism to learn the shrinkage of the cancer tumor based on patients' weekly cone-beam computed tomography (CBCT). We design a self-supervised domain adaption (SDA) method to learn and adapt the rich textural and spatial features from pre-treatment high-quality computed tomography (CT) to CBCT modality in order to address the poor image quality and lack of labels. We also provide uncertainty estimation for sequential segmentation, which aids not only in the risk management of treatment planning but also in the calibration and reliability of the model. Our experimental results based on a clinical non-small cell lung cancer (NSCLC) dataset with sixteen patients and ninety-six longitudinal CBCTs show that our model correctly learns weekly deformation of the tumor over time with an average dice score of 0.92 on the immediate next step, and is able to predict multiple steps (up to 5 weeks) for future patient treatments with an average dice score reduction of 0.05. By incorporating the tumor shrinkage predictions into a weekly re-planning strategy, our proposed method demonstrates a significant decrease in the risk of radiation-induced pneumonitis up to 35% while maintaining the high tumor control probability.

Cross-modality image translation: CT image synthesis of MR brain images using multi generative network with perceptual supervision

BACKGROUND

Computed tomography (CT) and magnetic resonance imaging (MRI) are the mainstream imaging technologies for clinical practice. CT imaging can reveal high-quality anatomical and physiopathological structures, especially bone tissue, for clinical diagnosis. MRI provides high resolution in soft tissue and is sensitive to lesions. CT combined with MRI diagnosis has become a regular image-guided radiation treatment plan.

METHODS

In this paper, to reduce the dose of radiation exposure in CT examinations and ameliorate the limitations of traditional virtual imaging technologies, we propose a Generative MRI-to-CT transformation method with structural perceptual supervision. Even though structural reconstruction is structurally misaligned in the MRI-CT dataset registration, our proposed method can better align structural information of synthetic CT (sCT) images to input MRI images while simulating the modality of CT in the MRI-to-CT cross-modality transformation.

RESULTS

We retrieved a total of 3416 brain MRI-CT paired images as the train/test dataset, including 1366 train images of 10 patients and 2050 test images of 15 patients. Several methods (the baseline methods and the proposed method) were evaluated by the HU difference map, HU distribution, and various similarity metrics, including the mean absolute error (MAE), structural similarity index (SSIM), peak signal-to-noise ratio (PSNR), and normalized cross-correlation (NCC). In our quantitative experimental results, the proposed method achieves the lowest MAE mean of 0.147, highest PSNR mean of 19.27, and NCC mean of 0.431 in the overall CT test dataset.

CONCLUSIONS

In conclusion, both qualitative and quantitative results of synthetic CT validate that the proposed method can preserve higher similarity of structural information of the bone tissue of target CT than the baseline methods. Furthermore, the proposed method provides better HU intensity reconstruction for simulating the distribution of the CT modality. The experimental estimation indicates that the proposed method is worth further investigation.

A ready-to-use machine learning tool for symmetric multi-modality registration of brain MRI

Volumetric registration of brain MRI is routinely used in human neuroimaging, e.g., to align different MRI modalities, to measure change in longitudinal analysis, to map an individual to a template, or in registration-based segmentation. Classical registration techniques based on numerical optimization have been very successful in this domain, and are implemented in widespread software suites like ANTs, Elastix, NiftyReg, or DARTEL. Over the last 7-8 years, learning-based techniques have emerged, which have a number of advantages like high computational efficiency, potential for higher accuracy, easy integration of supervision, and the ability to be part of a meta-architectures. However, their adoption in neuroimaging pipelines has so far been almost inexistent. Reasons include: lack of robustness to changes in MRI modality and resolution; lack of robust affine registration modules; lack of (guaranteed) symmetry; and, at a more practical level, the requirement of deep learning expertise that may be lacking at neuroimaging research sites. Here, we present EasyReg, an open-source, learning-based registration tool that can be easily used from the command line without any deep learning expertise or specific hardware. EasyReg combines the features of classical registration tools, the capabilities of modern deep learning methods, and the robustness to changes in MRI modality and resolution provided by our recent work in domain randomization. As a result, EasyReg is: fast; symmetric; diffeomorphic (and thus invertible); agnostic to MRI modality and resolution; compatible with affine and nonlinear registration; and does not require any preprocessing or parameter tuning. We present results on challenging registration tasks, showing that EasyReg is as accurate as classical methods when registering 1 mm isotropic scans within MRI modality, but much more accurate across modalities and resolutions. EasyReg is publicly available as part of FreeSurfer; see https://surfer.nmr.mgh.harvard.edu/fswiki/EasyReg .

Management of reirradiations: A clinical and technical overview based on a French survey

INTRODUCTION

The reirradiation number increased due to systemic therapies and patient survival. Few guidelines regarding acceptable cumulative doses to organs at risk (OARs) and appropriate dose accumulation tools need, made reirradiation challenging. The survey objective was to present the French current technical and clinical practices in reirradiations.

METHODS

A group of physician and physicists developed a survey gathering major issues of the topic. The questionnaire consisted in 4 parts: data collection, demographic, clinical and technical aspects. It was delivered through the SFRO and the SFPM. Data collection lasted 2 months and were gathered to compute statistical analysis.

RESULTS

48 institutions answered the survey. Difficulties about patient data collection were related to patient safety, administrative and technical limitations. Half of the institutions discussed reirradiation cases during a multidisciplinary meeting. It mainly aimed at discussing the indication and the new treatment total dose (92%). 79% of the respondents used various references but only 6% of them were specific to reirradiations. Patients with pain and clinical deficit were ranked as best inclusion criteria. 54.2% of the institutions considered OARs recovery, especially for spinal cord and brainstem. A commercial software was used for dose accumulation for 52% of respondents. Almost all institutions performed equivalent dose conversion (94%). A quarter of the institutions estimated not to have the appropriate equipment for reirradiation.

CONCLUSION

This survey showed the various approaches and tools used in reirradiation management. It highlighted issues in collecting data, and the guidelines necessity for safe practices, to increase clinicians confidence in retreating patients.

STereotactic Arrhythmia Radioablation (STAR): the Standardized Treatment and Outcome Platform for Stereotactic Therapy Of Re-entrant tachycardia by a Multidisciplinary consortium (STOPSTORM.eu) and review of current patterns of STAR practice in Europe

The EU Horizon 2020 Framework-funded Standardized Treatment and Outcome Platform for Stereotactic Therapy Of Re-entrant tachycardia by a Multidisciplinary (STOPSTORM) consortium has been established as a large research network for investigating STereotactic Arrhythmia Radioablation (STAR) for ventricular tachycardia (VT). The aim is to provide a pooled treatment database to evaluate patterns of practice and outcomes of STAR and finally to harmonize STAR within Europe. The consortium comprises 31 clinical and research institutions. The project is divided into nine work packages (WPs): (i) observational cohort; (ii) standardization and harmonization of target delineation; (iii) harmonized prospective cohort; (iv) quality assurance (QA); (v) analysis and evaluation; (vi, ix) ethics and regulations; and (vii, viii) project coordination and dissemination. To provide a review of current clinical STAR practice in Europe, a comprehensive questionnaire was performed at project start. The STOPSTORM Institutions' experience in VT catheter ablation (83% ≥ 20 ann.) and stereotactic body radiotherapy (59% > 200 ann.) was adequate, and 84 STAR treatments were performed until project launch, while 8/22 centres already recruited VT patients in national clinical trials. The majority currently base their target definition on mapping during VT (96%) and/or pace mapping (75%), reduced voltage areas (63%), or late ventricular potentials (75%) during sinus rhythm. The majority currently apply a single-fraction dose of 25 Gy while planning techniques and dose prescription methods vary greatly. The current clinical STAR practice in the STOPSTORM consortium highlights potential areas of optimization and harmonization for substrate mapping, target delineation, motion management, dosimetry, and QA, which will be addressed in the various WPs.

Analytic calculations and synthetic shapes for validation of quantitative contour comparison software

A high level of variability in reported values was observed in a recent survey of contour similarity measures (CSMs) calculation tools. Such variations in the output measurements prevent meaningful comparison between studies. The purpose of this study was to develop a dataset with analytically calculated gold standard values to facilitate standardization and ensure accuracy of CSM implementations. The dataset was generated in the Digital Imaging and Communications in Medicine (DICOM) format. Both the dataset and the software used for its generation are made publicly available to encourage robust testing of CSM implementations for accuracy, improving consistency between different implementations.

On the implant stability in adaptive multi-catheter breast brachytherapy: Establishment of a decision-tree for treatment re-planning

BACKGROUND AND PURPOSE

To assess implant stability and identify causes of implant variations during high-dose-rate multi-catheter breast brachytherapy.

MATERIALS AND METHODS

Planning-CTs were compared to control-CTs acquired halfway through the treatment for 100 patients. For assessing geometric stability, Fréchet-distance and button-to-button distance changes of all catheters as well as variations of Euclidean distances and convex hulls of all dwell positions were determined. The CTs were inspected to identify the causes of geometric changes. Dosimetric effects were evaluated by target volume transfers and re-contouring of organs at risk. The dose non-uniformity ratio (DNR), 100% and 150% isodose volumes (V₁₀₀ and V₁₅₀), coverage index (CI), and organ doses were calculated. Correlations between the examined geometric and dosimetric parameters were assessed.

RESULTS

Fréchet-distance and dwell position deviations >2.5 mm as well as button-to-button distance changes >5 mm were detected for 5%, 2%, and 6.3% of catheters, but for 32, 17, and 37 patients, respectively. Variations occurred enhanced in the lateral breast and close to the ribs, e.g. due to different arm positions. Only small dosimetric effects with median DNR, V₁₀₀, and CI variations of -0.01 ± 0.02, (-0.5 ± 1.3)ccm, and (-1.4 ± 1.8)% were observed in general. Skin dose exceeded recommended levels for 12 of 100 patients. Various correlations between geometric and dosimetric implant stability were found, based on which decision-tree regarding treatment re-planning was established.

CONCLUSION

Multi-catheter breast brachytherapy shows a high implant stability in general, but considering skin dose changes is important. To increase implant stability for individual patients, we plan to investigate patient immobilization aids during treatments.

Development of a multi-layer quality assurance program to evaluate the uncertainty of deformable dose accumulation in adaptive radiotherapy

PURPOSE

Deformable dose accumulation (DDA) has uncertainties which impede the implementation of DDA-based adaptive radiotherapy (ART) in clinic. The purpose of this study is to develop a multi-layer quality assurance (MLQA) program to evaluate uncertainties in DDA.

METHODS

A computer program is developed to generate a pseudo-inverse displacement vector field (DVF) for each deformable image registration (DIR) performed in Accuray's PreciseART. The pseudo-inverse DVF is first used to calculate a pseudo-inverse consistency error (PICE) and then implemented in an energy and mass congruent mapping (EMCM) method to reconstruct a deformed dose. The PICE is taken as a metric to estimate DIR uncertainties. A pseudo-inverse dose agreement rate (PIDAR) is used to evaluate the consequence of the DIR uncertainties in DDA and the principle of energy conservation is used to validate the integrity of dose mappings. The developed MLQA program was tested using the data collected from five representative cancer patients treated with tomotherapy.

RESULTS

DIRs were performed in PreciseART to generate primary DVFs for the five patients. The fidelity index and PICE of these DVFs on average are equal to 0.028 mm and 0.169 mm, respectively. With the criteria of 3 mm/3% and 5 mm/5%, the PIDARs of the PreciseART-reconstructed doses are 73.9 ± 4.4% and 87.2 ± 3.3%, respectively. The PreciseART and EMCM-based dose reconstructions have their deposited energy changed by 5.6 ± 3.9% and 2.6 ± 1.5% in five GTVs, and by 9.2 ± 7.8% and 4.7 ± 3.6% in 30 OARs, respectively.

CONCLUSIONS

A pseudo-inverse map-based EMCM program has been developed to evaluate DIR and dose mapping uncertainties. This program could also be used as a sanity check tool for DDA-based ART.

Evaluating the accumulated dose distribution of organs at risk in combined radiotherapy for cervical carcinoma based on deformable image registration

OBJECTIVE

To evaluate the feasibility and value of deformable image registration (DIR) in calculating the cumulative doses of organs at risk (OARs) in the combined radiotherapy of cervical cancer.

PATIENTS AND METHODS

Thirty cervical cancer patients treated with external beam radiotherapy (EBRT) combined with intracavitary brachytherapy (ICBT) were reviewed. The simulation CT images of EBRT and ICBT were imported into Varian Velocity 4.1 for the DIR-based dose accumulation. Cumulative dose-volume parameters of D2cc for rectum and bladder were compared between the direct addition (DA) and DIR methods. The quantitative parameters were measured to evaluate the accuracy of DIR.

RESULTS

The three-dimensional cumulative dose distribution of the tumor and OARs were graphically well illustrated by composite isodose lines. In combined EBRT and ICBT, the mean cumulative bladder D2cc calculated by DIR and DA was 86.13 Gy and 86.27 Gy, respectively. The mean cumulative rectal D2cc calculated by DIR and DA was 72.97 Gy and 73.90 Gy, respectively. No significant differences were noted between these two methods (p > 0.05). As to the parameters used to evaluate the DIR accuracy, the mean DSC, Jacobian, MDA (mm) and Hausdorff distance (mm) were 0.79, 1.0, 3.84, and 22.01 respectively for the bladder and 0.53, 1.2, 7.31, and 29.58 respectively for the rectum. In this study, the DSC seemed to be slightly lower compared with previous studies.

CONCLUSION

Dose accumulation based on DIR might be an alternative method to illustrate and evaluate the cumulative doses of the OARs in combined radiotherapy for cervical cancer. However, DIR should be used with caution before overcoming the relevant limitations.

HaN-Seg: The head and neck organ-at-risk CT and MR segmentation dataset

PURPOSE

For the cancer in the head and neck (HaN), radiotherapy (RT) represents an important treatment modality. Segmentation of organs-at-risk (OARs) is the starting point of RT planning, however, existing approaches are focused on either computed tomography (CT) or magnetic resonance (MR) images, while multimodal segmentation has not been thoroughly explored yet. We present a dataset of CT and MR images of the same patients with curated reference HaN OAR segmentations for an objective evaluation of segmentation methods.

ACQUISITION AND VALIDATION METHODS

The cohort consists of HaN images of 56 patients that underwent both CT and T1-weighted MR imaging for image-guided RT. For each patient, reference segmentations of up to 30 OARs were obtained by experts performing manual pixel-wise image annotation. By maintaining the distribution of patient age and gender, and annotation type, the patients were randomly split into training Set 1 (42 cases or 75%) and test Set 2 (14 cases or 25%). Baseline auto-segmentation results are also provided by training the publicly available deep nnU-Net architecture on Set 1, and evaluating its performance on Set 2.

DATA FORMAT AND USAGE NOTES

The data are publicly available through an open-access repository under the name HaN-Seg: The Head and Neck Organ-at-Risk CT & MR Segmentation Dataset. Images and reference segmentations are stored in the NRRD file format, where the OAR filenames correspond to the nomenclature recommended by the American Association of Physicists in Medicine, and OAR and demographics information is stored in separate comma-separated value  files.

POTENTIAL APPLICATIONS

The HaN-Seg: The Head and Neck Organ-at-Risk CT & MR Segmentation Challenge is launched in parallel with the dataset release to promote the development of automated techniques for OAR segmentation in the HaN. Other potential applications include out-of-challenge algorithm development and benchmarking, as well as external validation of the developed algorithms.

AAPM medical physics practice guideline 13.a: HDR brachytherapy, part A

The American Association of Physicists in Medicine (AAPM) is a nonprofit professional society whose primary purposes are to advance the science, education, and professional practice of medical physics. The AAPM has more than 8000 members and is the principal organization of medical physicists in the United States. The AAPM will periodically define new practice guidelines for medical physics practice to help advance the science of medical physics and to improve the quality of service to patients throughout the United States. Existing medical physics practice guidelines (MPPGs) will be reviewed for the purpose of revision or renewal, as appropriate, on their fifth anniversary or sooner. Each medical physics practice guideline represents a policy statement by the AAPM, has undergone a thorough consensus process in which it has been subjected to extensive review, and requires the approval of the Professional Council. The medical physics practice guidelines recognize that the safe and effective use of diagnostic and therapeutic radiology requires specific training, skills, and techniques, as described in each document. Reproduction or modification of the published practice guidelines and technical standards by those entities not providing these services is not authorized. The following terms are used in the AAPM practice guidelines: (1) Must and must not: Used to indicate that adherence to the recommendation is considered necessary to conform to this practice guideline. (2) Should and should not: Used to indicate a prudent practice to which exceptions may occasionally be made in appropriate circumstances. Approved by AAPM's Executive Committee April 28, 2022.

Quality and Safety Considerations in Image Guided Radiation Therapy: An ASTRO Safety White Paper Update

PURPOSE

This updated report on image guided radiation therapy (IGRT) is part of a series of consensus-based white papers previously published by the American Society for Radiation Oncology addressing patient safety. Since the first white papers were published, IGRT technology and procedures have progressed significantly such that these procedures are now more commonly used. The use of IGRT has now extended beyond high-precision treatments, such as stereotactic radiosurgery and stereotactic body radiation therapy, and into routine clinical practice for many treatment techniques and anatomic sites. Therefore, quality and patient safety considerations for these techniques remain an important area of focus.

METHODS AND MATERIALS

The American Society for Radiation Oncology convened an interdisciplinary task force to assess the original IGRT white paper and update content where appropriate. Recommendations were created using a consensus-building methodology, and task force members indicated their level of agreement based on a 5-point Likert scale from "strongly agree" to "strongly disagree." A prespecified threshold of ≥75% of raters who selected "strongly agree" or "agree" indicated consensus.

SUMMARY

This IGRT white paper builds on the previous version and uses other guidance documents to primarily focus on processes related to quality and safety. IGRT requires an interdisciplinary team-based approach, staffed by appropriately trained specialists, as well as significant personnel resources, specialized technology, and implementation time. A thorough feasibility analysis of resources is required to achieve the clinical and technical goals and should be discussed with all personnel before undertaking new imaging techniques. A comprehensive quality-assurance program must be developed, using established guidance, to ensure IGRT is performed in a safe and effective manner. As IGRT technologies continue to improve or emerge, existing practice guidelines should be reviewed or updated regularly according to the latest American Association of Physicists in Medicine Task Group reports or guidelines. Patient safety in the application of IGRT is everyone's responsibility, and professional organizations, regulators, vendors, and end-users must demonstrate a clear commitment to working together to ensure the highest levels of safety.

An unsupervised image registration method employing chest computed tomography images and deep neural networks

BACKGROUND

Deformable image registration is crucial for multiple radiation therapy applications. Fast registration of computed tomography (CT) lung images is challenging because of the large and nonlinear deformation between inspiration and expiration. With advancements in deep learning techniques, learning-based registration methods are considered efficient alternatives to traditional methods in terms of accuracy and computational cost.

METHOD

In this study, an unsupervised lung registration network (LRN) with cycle-consistent training is proposed to align two acquired CT-derived lung datasets during breath-holds at inspiratory and expiratory levels without utilizing any ground-truth registration results. Generally, the LRN model uses three loss functions: image similarity, regularization, and Jacobian determinant. Here, LRN was trained on the CT datasets of 705 subjects and tested using 10 pairs of public CT DIR-Lab datasets. Furthermore, to evaluate the effectiveness of the registration technique, target registration errors (TREs) of the LRN model were compared with those of the conventional algorithm (sum of squared tissue volume difference; SSTVD) and a state-of-the-art unsupervised registration method (VoxelMorph).

RESULTS

The results showed that the LRN with an average TRE of 1.78 ± 1.56 mm outperformed VoxelMorph with an average TRE of 2.43 ± 2.43 mm, which is comparable to that of SSTVD with an average TRE of 1.66 ± 1.49 mm. In addition, estimating the displacement vector field without any folding voxel consumed less than 2 s, demonstrating the superiority of the learning-based method with respect to fiducial marker tracking and the overall soft tissue alignment with a nearly real-time speed.

CONCLUSIONS

Therefore, this proposed method shows significant potential for use in time-sensitive pulmonary studies, such as lung motion tracking and image-guided surgery.

Effects of respiratory and cardiac motion on estimating radiation dose to the left ventricle during radiotherapy for lung cancer

PURPOSE

Establish a workflow to evaluate radiotherapy (RT) dose variation induced by respiratory and cardiac motion on the left ventricle (LV) and left ventricular myocardium (LVM).

METHODS

Eight lung cancer patients underwent 4D-CT, expiratory T1-volumetric-interpolated-breath-hold-examination (VIBE), and cine MRI scans in expiration. Treatment plans were designed on the average intensity projection (AIP) datasets from 4D-CTs. RT dose from AIP was transferred onto 4D-CT respiratory phases. About 50% 4D-CT dose was mapped onto T1-VIBE (following registration) and from there onto average cine MRI datasets. Dose from average cine MRI was transferred onto all cardiac phases. Cumulative cardiac dose was estimated by transferring dose from each cardiac phase onto a reference cine phase following deformable image registration. The LV was contoured on each 4D-CT breathing phase and was called clinical LV (cLV); this structure is blurred by cardiac motion. Additionally, LV, LVM, and an American Heart Association (AHA) model were contoured on all cardiac phases. Relative maximum/mean doses for contoured regions were calculated with respect to each patient's maximum/mean AIP dose.

RESULTS

During respiration, relative maximum and mean doses on the cLV ranged from -4.5% to 5.6% and -14.2% to 16.5%, respectively, with significant differences in relative mean doses between inspiration and expiration (P < 0.0145). During cardiac motion at expiration, relative maximum and mean doses on the LV ranged from 1.6% to 59.3%, 0.5% to 27.4%, respectively. Relative mean doses were significantly different between diastole and systole (P = 0.0157). No significant differences were noted between systolic, diastolic, or cumulative cardiac doses compared to the expiratory 4D-CT (P > 0.14). Significant differences were observed in AHA segmental doses depending on tumour proximity compared to global LV doses on expiratory 4D-CT (P < 0.0117).

CONCLUSION

In this study, the LV dose was highest during expiration and diastole. Segmental evaluation suggested that future cardiotoxicity evaluations may benefit from regional assessments of dose that account for cardiopulmonary motion.

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.