MR-guided stereotactic radiation therapy for head and neck cancers.

Purpose: MR-guided radiotherapy (MRgRT) has the advantage of utilizing high soft tissue contrast imaging to track daily changes in target and critical organs throughout the entire radiation treatment course. Head and neck (HN) stereotactic body radiation therapy (SBRT) has been increasingly used to treat localized lesions within a shorter timeframe. The purpose of this study is to examine the dosimetric difference between the step-and-shot intensity modulated radiation therapy (IMRT) plans on Elekta Unity and our clinical volumetric modulated arc therapy (VMAT) plans on Varian TrueBeam for HN SBRT. Method: Fourteen patients treated on TrueBeam sTx with VMAT treatment plans were re-planned in the Monaco treatment planning system for Elekta Unity MR-Linac (MRL). The plan qualities, including target coverage, conformity, homogeneity, nearby critical organ doses, gradient index and low dose bath volume, were compared between VMAT and Monaco IMRT plans. Additionally, we evaluated the Unity adaptive plans of adapt-to-position (ATP) and adapt-to-shape (ATS) workflows using simulated setup errors for five patients and assessed the outcomes of our treated patients. Results: Monaco IMRT plans achieved comparable results to VMAT plans in terms of target coverage, uniformity and homogeneity, with slightly higher target maximum and mean doses. The critical organ doses in Monaco IMRT plans all met clinical goals; however, the mean doses and low dose bath volumes were higher than in VMAT plans. The adaptive plans demonstrated that the ATP workflow may result in degraded target coverage and OAR doses for HN SBRT, while the ATS workflow can maintain the plan quality. Conclusion: The use of Monaco treatment planning and online adaptation can achieve dosimetric results comparable to VMAT plans, with the additional benefits of real-time tracking of target volume and nearby critical structures. This offers the potential to treat aggressive and variable tumors in HN SBRT and improve local control and treatment toxicity.

Quantifying the Effect of 4-Dimensional Computed Tomography-Based Deformable Dose Accumulation on Representing Radiation Damage for Patients with Locally Advanced Non-Small Cell Lung Cancer Treated with Standard-Fractionated Intensity-Modulated Radiation Therapy.

PURPOSE: The aim of this study was to investigate the dosimetric and clinical effects of 4-dimensional computed tomography (4DCT)-based longitudinal dose accumulation in patients with locally advanced non-small cell lung cancer treated with standard-fractionated intensity-modulated radiation therapy (IMRT). METHODS AND MATERIALS: Sixty-seven patients were retrospectively selected from a randomized clinical trial. Their original IMRT plan, planning and verification 4DCTs, and ∼4-month posttreatment follow-up CTs were imported into a commercial treatment planning system. Two deformable image registration algorithms were implemented for dose accumulation, and their accuracies were assessed. The planned and accumulated doses computed using average-intensity images or phase images were compared. At the organ level, mean lung dose and normal-tissue complication probability (NTCP) for grade ≥2 radiation pneumonitis were compared. At the region level, mean dose in lung subsections and the volumetric overlap between isodose intervals were compared. At the voxel level, the accuracy in estimating the delivered dose was compared by evaluating the fit of a dose versus radiographic image density change (IDC) model. The dose-IDC model fit was also compared for subcohorts based on the magnitude of NTCP difference (|ΔNTCP|) between planned and accumulated doses. RESULTS: Deformable image registration accuracy was quantified, and the uncertainty was considered for the voxel-level analysis. Compared with planned doses, accumulated doses on average resulted in <1-Gy lung dose increase and <2% NTCP increase (up to 8.2 Gy and 18.8% for a patient, respectively). Volumetric overlap of isodose intervals between the planned and accumulated dose distributions ranged from 0.01 to 0.93. Voxel-level dose-IDC models demonstrated a fit improvement from planned dose to accumulated dose (pseudo-R2 increased 0.0023) and a further improvement for patients with ≥2% |ΔNTCP| versus for patients with <2% |ΔNTCP|. CONCLUSIONS: With a relatively large cohort, robust image registrations, multilevel metric comparisons, and radiographic image-based evidence, we demonstrated that dose accumulation more accurately represents the delivered dose and can be especially beneficial for patients with greater longitudinal response.

Validation of an automated contouring and treatment planning tool for pediatric craniospinal radiation therapy.

Purpose: Treatment planning for craniospinal irradiation (CSI) is complex and time-consuming, especially for resource-constrained centers. To alleviate demanding workflows, we successfully automated the pediatric CSI planning pipeline in previous work. In this work, we validated our CSI autosegmentation and autoplanning tool on a large dataset from St. Jude Children's Research Hospital. Methods: Sixty-three CSI patient CT scans were involved in the study. Pre-planning scripts were used to automatically verify anatomical compatibility with the autoplanning tool. The autoplanning pipeline generated 15 contours and a composite CSI treatment plan for each of the compatible test patients (n=51). Plan quality was evaluated quantitatively with target coverage and dose to normal tissue metrics and qualitatively with physician review, using a 5-point Likert scale. Three pediatric radiation oncologists from 3 institutions reviewed and scored 15 contours and a corresponding composite CSI plan for the final 51 test patients. One patient was scored by 3 physicians, resulting in 53 plans scored total. Results: The algorithm automatically detected 12 incompatible patients due to insufficient junction spacing or head tilt and removed them from the study. Of the 795 autosegmented contours reviewed, 97% were scored as clinically acceptable, with 92% requiring no edits. Of the 53 plans scored, all 51 brain dose distributions were scored as clinically acceptable. For the spine dose distributions, 92%, 100%, and 68% of single, extended, and multiple-field cases, respectively, were scored as clinically acceptable. In all cases (major or minor edits), the physicians noted that they would rather edit the autoplan than create a new plan. Conclusions: We successfully validated an autoplanning pipeline on 51 patients from another institution, indicating that our algorithm is robust in its adjustment to differing patient populations. We automatically generated 15 contours and a comprehensive CSI treatment plan for each patient without physician intervention, indicating the potential for increased treatment planning efficiency and global access to high-quality radiation therapy.

Identifying the optimal deep learning architecture and parameters for automatic beam aperture definition in 3D radiotherapy.

PURPOSE: Two-dimensional radiotherapy is often used to treat cervical cancer in low- and middle-income countries, but treatment planning can be challenging and time-consuming. Neural networks offer the potential to greatly decrease planning time through automation, but the impact of the wide range of hyperparameters to be set during training on model accuracy has not been exhaustively investigated. In the current study, we evaluated the effect of several convolutional neural network architectures and hyperparameters on 2D radiotherapy treatment field delineation. METHODS: Six commonly used deep learning architectures were trained to delineate four-field box apertures on digitally reconstructed radiographs for cervical cancer radiotherapy. A comprehensive search of optimal hyperparameters for all models was conducted by varying the initial learning rate, image normalization methods, and (when appropriate) convolutional kernel size, the number of learnable parameters via network depth and the number of feature maps per convolution, and nonlinear activation functions. This yielded over 1700 unique models, which were all trained until performance converged and then tested on a separate dataset. RESULTS: Of all hyperparameters, the choice of initial learning rate was most consistently significant for improved performance on the test set, with all top-performing models using learning rates of 0.0001. The optimal image normalization was not consistent across architectures. High overlap (mean Dice similarity coefficient = 0.98) and surface distance agreement (mean surface distance < 2 mm) were achieved between the treatment field apertures for all architectures using the identified best hyperparameters. Overlap Dice similarity coefficient (DSC) and distance metrics (mean surface distance and Hausdorff distance) indicated that DeepLabv3+ and D-LinkNet architectures were least sensitive to initial hyperparameter selection. CONCLUSION: DeepLabv3+ and D-LinkNet are most robust to initial hyperparameter selection. Learning rate, nonlinear activation function, and kernel size are also important hyperparameters for improving performance.

Feasibility and Clinical Acceptability of Automation-Assisted 3D Conformal Radiotherapy Planning for Patients With Cervical Cancer in a Resource-Constrained Setting.

PURPOSE: The Ocean Road Cancer Institute (ORCI) in Tanzania began offering 3D conformal radiation therapy (3DCRT) in 2018. Steep learning curves, high patient volume, and a limited workforce resulted in long radiation therapy (RT) planning workflows. We aimed to establish the feasibility of implementing an automation-assisted cervical cancer 3DCRT planning system. MATERIALS AND METHODS: We performed chart abstractions on 30 patients with cervical cancer treated with 3DCRT at ORCI. The Radiation Planning Assistant (RPA) generated a new automated set of contours and plans on the basis of anonymized computed tomography images. Each were assessed for edit time requirements, dose-volume safety metrics, and clinical acceptability by two ORCI physician investigators. Dice similarity coefficient (DSC) agreement analysis was conducted between original and new contour sets. RESULTS: The average time to manually develop treatment plans was 7 days. Applying RPA, automated same-day contours and plans were developed for 29 of 30 patients (97%). Of the 29 evaluable contours, all were approved with <2 minutes of edit time. Agreement between clinical and RPA contours was highest for the rectum (median DSC, 0.72) and bladder (DSC, 0.90). Agreement was lower with the primary tumor clinical target volume (CTVp; DSC, 0.69) and elective nodal clinical target volume (CTVn; DSC, 0.63). All RPA plans were approved with <4 minutes of edit time. RPA target coverage was excellent, covering the CTVp with median V45 Gy 100% and CTVn with median V45 Gy 99.9%. CONCLUSION: Automation-assisted 3DCRT contouring yielded high levels of agreement for normal structures. The RPA met all planning safety metrics and sustained high levels of clinical acceptability with minimal edit times. This tool offers the potential to significantly decrease RT planning timelines while maintaining high-quality RT delivery in resource-constrained settings.

A real-time contouring feedback tool for consensus-based contour training.

Purpose: Variability in contouring structures of interest for radiotherapy continues to be challenging. Although training can reduce such variability, having radiation oncologists provide feedback can be impractical. We developed a contour training tool to provide real-time feedback to trainees, thereby reducing variability in contouring. Methods: We developed a novel metric termed localized signed square distance (LSSD) to provide feedback to the trainee on how their contour compares with a reference contour, which is generated real-time by combining trainee contour and multiple expert radiation oncologist contours. Nine trainees performed contour training by using six randomly assigned training cases that included one test case of the heart and left ventricle (LV). The test case was repeated 30 days later to assess retention. The distribution of LSSD maps of the initial contour for the training cases was combined and compared with the distribution of LSSD maps of the final contours for all training cases. The difference in standard deviations from the initial to final LSSD maps, ΔLSSD, was computed both on a per-case basis and for the entire group. Results: For every training case, statistically significant ΔLSSD were observed for both the heart and LV. When all initial and final LSSD maps were aggregated for the training cases, before training, the mean LSSD ([range], standard deviation) was -0.8 mm ([-37.9, 34.9], 4.2) and 0.3 mm ([-25.1, 32.7], 4.8) for heart and LV, respectively. These were reduced to -0.1 mm ([-16.2, 7.3], 0.8) and 0.1 mm ([-6.6, 8.3], 0.7) for the final LSSD maps during the contour training sessions. For the retention case, the initial and final LSSD maps of the retention case were aggregated and were -1.5 mm ([-22.9, 19.9], 3.4) and -0.2 mm ([-4.5, 1.5], 0.7) for the heart and 1.8 mm ([-16.7, 34.5], 5.1) and 0.2 mm ([-3.9, 1.6],0.7) for the LV. Conclusions: A tool that uses real-time contouring feedback was developed and successfully used for contour training of nine trainees. In all cases, the utility was able to guide the trainee and ultimately reduce the variability of the trainee's contouring.

Automated Contouring and Planning in Radiation Therapy: What Is 'Clinically Acceptable'?

Developers and users of artificial-intelligence-based tools for automatic contouring and treatment planning in radiotherapy are expected to assess clinical acceptability of these tools. However, what is 'clinical acceptability'? Quantitative and qualitative approaches have been used to assess this ill-defined concept, all of which have advantages and disadvantages or limitations. The approach chosen may depend on the goal of the study as well as on available resources. In this paper, we discuss various aspects of 'clinical acceptability' and how they can move us toward a standard for defining clinical acceptability of new autocontouring and planning tools.

Dose Escalation for Pancreas SBRT: Potential and Limitations of using Daily Online Adaptive Radiation Therapy and an Iterative Isotoxicity Automated Planning Approach.

Purpose: To determine the dosimetric limitations of daily online adaptive pancreas stereotactic body radiation treatment by using an automated dose escalation approach. Methods and Materials: We collected 108 planning and daily computed tomography (CT) scans from 18 patients (18 patients × 6 CT scans) who received 5-fraction pancreas stereotactic body radiation treatment at MD Anderson Cancer Center. Dose metrics from the original non-dose-escalated clinical plan (non-DE), the dose-escalated plan created on the original planning CT (DE-ORI), and the dose-escalated plan created on daily adaptive radiation therapy CT (DE-ART) were analyzed. We developed a dose-escalation planning algorithm within the radiation treatment planning system to automate the dose-escalation planning process for efficiency and consistency. In this algorithm, the prescription dose of the dose-escalation plan was escalated before violating any organ-at-risk (OAR) dose constraint. Dose metrics for 3 targets (gross target volume [GTV], tumor vessel interface [TVI], and dose-escalated planning target volume [DE-PTV]) and 9 OARs (duodenum, large bowel, small bowel, stomach, spinal cord, kidneys, liver, and skin) for the 3 plans were compared. Furthermore, we evaluated the effectiveness of the online adaptive dose-escalation planning process by quantifying the effect of the interfractional dose distribution variations among the DE-ART plans. Results: The median D95% dose to the GTV/TVI/DE-PTV was 33.1/36.2/32.4 Gy, 48.5/50.9/40.4 Gy, and 53.7/58.2/44.8 Gy for non-DE, DE-ORI, and DE-ART, respectively. Most OAR dose constraints were not violated for the non-DE and DE-ART plans, while OAR constraints were violated for the majority of the DE-ORI patients due to interfractional motion and lack of adaptation. The maximum difference per fraction in D95%, due to interfractional motion, was 2.5 ± 2.7 Gy, 3.0 ± 2.9 Gy, and 2.0 ± 1.8 Gy for the TVI, GTV, and DE-PTV, respectively. Conclusions: Most patients require daily adaptation of the radiation planning process to maximally escalate delivered dose to the pancreatic tumor without exceeding OAR constraints. Using our automated approach, patients can receive higher target dose than standard of care without violating OAR constraints.

Automating the treatment planning process for 3D-conformal pediatric craniospinal irradiation therapy.

PURPOSE: Pediatric patients with medulloblastoma in low- and middle-income countries (LMICs) are most treated with 3D-conformal photon craniospinal irradiation (CSI), a time-consuming, complex treatment to plan, especially in resource-constrained settings. Therefore, we developed and tested a 3D-conformal CSI autoplanning tool for varying patient lengths. METHODS AND MATERIALS: Autocontours were generated with a deep learning model trained:tested (80:20 ratio) on 143 pediatric medulloblastoma CT scans (patient ages: 2-19 years, median = 7 years). Using the verified autocontours, the autoplanning tool generated two lateral brain fields matched to a single spine field, an extended single spine field, or two matched spine fields. Additional spine subfields were added to optimize the corresponding dose distribution. Feathering was implemented (yielding nine to 12 fields) to give a composite plan. Each planning approach was tested on six patients (ages 3-10 years). A pediatric radiation oncologist assessed clinical acceptability of each autoplan. RESULTS: The autocontoured structures' average Dice similarity coefficient ranged from .65 to .98. The average V95 for the brain/spinal canal for single, extended, and multi-field spine configurations was 99.9% ± 0.06%/99.9% ± 0.10%, 99.9% ± 0.07%/99.4% ± 0.30%, and 99.9% ± 0.06%/99.4% ± 0.40%, respectively. The average maximum dose across all field configurations to the brainstem, eyes (L/R), lenses (L/R), and spinal cord were 23.7 ± 0.08, 24.1 ± 0.28, 13.3 ± 5.27, and 25.5 ± 0.34 Gy, respectively (prescription = 23.4 Gy/13 fractions). Of the 18 plans tested, all were scored as clinically acceptable as-is or clinically acceptable with minor, time-efficient edits preferred or required. No plans were scored as clinically unacceptable. CONCLUSION: The autoplanning tool successfully generated pediatric CSI plans for varying patient lengths in 3.50 ± 0.4 minutes on average, indicating potential for an efficient planning aid in a resource-constrained settings.

Customizable landmark-based field aperture design for automated whole-brain radiotherapy treatment planning.

PURPOSE: To develop and evaluate an automated whole-brain radiotherapy (WBRT) treatment planning pipeline with a deep learning-based auto-contouring and customizable landmark-based field aperture design. METHODS: The pipeline consisted of the following steps: (1) Auto-contour normal structures on computed tomography scans and digitally reconstructed radiographs using deep learning techniques, (2) locate the landmark structures using the beam's-eye-view, (3) generate field apertures based on eight different landmark rules addressing different clinical purposes and physician preferences. Two parallel approaches for generating field apertures were developed for quality control. The performance of the generated field shapes and dose distributions were compared with the original clinical plans. The clinical acceptability of the plans was assessed by five radiation oncologists from four hospitals. RESULTS: The performance of the generated field apertures was evaluated by the Hausdorff distance (HD) and mean surface distance (MSD) from 182 patients' field apertures used in the clinic. The average HD and MSD for the generated field apertures were 16 ± 7 and 7 ± 3 mm for the first approach, respectively, and 17 ± 7 and 7 ± 3 mm, respectively, for the second approach. The differences regarding HD and MSD between the first and the second approaches were 1 ± 2 and 1 ± 3 mm, respectively. A clinical review of the field aperture design, conducted using 30 patients, achieved a 100% acceptance rate for both the first and second approaches, and the plan review achieved a 100% acceptance rate for the first approach and a 93% acceptance rate for the second approach. The average acceptance rate for meeting lens dosimetric recommendations was 80% (left lens) and 77% (right lens) for the first approach, and 70% (both left and right lenses) for the second approach, compared with 50% (left lens) and 53% (right lens) for the clinical plans. CONCLUSION: This study provided an automated pipeline with two field aperture generation approaches to automatically generate WBRT treatment plans. Both quantitative and qualitative evaluations demonstrated that our novel pipeline was comparable with the original clinical plans.

Multi-organ segmentation of abdominal structures from non-contrast and contrast enhanced CT images.

Manually delineating upper abdominal organs at risk (OARs) is a time-consuming task. To develop a deep-learning-based tool for accurate and robust auto-segmentation of these OARs, forty pancreatic cancer patients with contrast-enhanced breath-hold computed tomographic (CT) images were selected. We trained a three-dimensional (3D) U-Net ensemble that automatically segments all organ contours concurrently with the self-configuring nnU-Net framework. Our tool's performance was assessed on a held-out test set of 30 patients quantitatively. Five radiation oncologists from three different institutions assessed the performance of the tool using a 5-point Likert scale on an additional 75 randomly selected test patients. The mean (± std. dev.) Dice similarity coefficient values between the automatic segmentation and the ground truth on contrast-enhanced CT images were 0.80 ± 0.08, 0.89 ± 0.05, 0.90 ± 0.06, 0.92 ± 0.03, 0.96 ± 0.01, 0.97 ± 0.01, 0.96 ± 0.01, and 0.96 ± 0.01 for the duodenum, small bowel, large bowel, stomach, liver, spleen, right kidney, and left kidney, respectively. 89.3% (contrast-enhanced) and 85.3% (non-contrast-enhanced) of duodenum contours were scored as a 3 or above, which required only minor edits. More than 90% of the other organs' contours were scored as a 3 or above. Our tool achieved a high level of clinical acceptability with a small training dataset and provides accurate contours for treatment planning.

An Automated Treatment Planning Framework for Spinal Radiation Therapy and Vertebral-Level Second Check.

PURPOSE: Complicating factors such as time pressures, anatomic variants in the spine, and similarities in adjacent vertebrae are associated with incorrect level treatments of the spine. The purpose of this work was to mitigate such challenges by fully automating the treatment planning process for diagnostic and simulation computed tomography (CT) scans. METHODS AND MATERIALS: Vertebral bodies are labeled on CT scans of any length using 2 intendent deep-learning models-mirroring 2 different experts labeling the spine. Then, a U-Net++ architecture was trained, validated, and tested to contour each vertebra (n = 220 CT scans). Features from the CT and auto-contours were input into a random forest classifier to predict whether vertebrae were correctly labeled. This classifier was trained using auto-contours from cone beam computed tomography, positron emission tomography/CT, simulation CT, and diagnostic CT images (n = 56 CT scans, 751 contours). Auto-plans were generated via scripting. Each model was combined into a framework to make a fully automated clinical tool. A retrospective planning study was conducted in which 3 radiation oncologists scored auto-plan quality on an unseen patient cohort (n = 60) on a 5-point scale. CT scans varied in scan length, presence of surgical implants, imaging protocol, and metastatic burden. RESULTS: The results showed that the uniquely designed convolutional neural networks accurately labeled and segmented vertebral bodies C1-L5 regardless of imaging protocol or metastatic burden. Mean dice-similarity coefficient was 85.0% (cervical), 90.3% (thoracic), and 93.7% (lumbar). The random forest classifier predicted mislabeling across various CT scan types with an area under the curve of 0.82. All contouring and labeling errors within treatment regions (11 of 11), including errors from patient plans with atypical anatomy (eg, T13, L6) were detected. Radiation oncologists scored 98% of simulation CT-based plans and 92% of diagnostic CT-based plans as clinically acceptable or needing minor edits for patients with typical anatomy. On average, end-to-end treatment planning time of the clinical tool was less than 8 minutes. CONCLUSIONS: This novel method to automatically verify, contour, and plan palliative spine treatments is efficient and effective across various CT scan types. Furthermore, it is the first to create a clinical tool that can automatically verify vertebral level in CT images.

Automation of radiation treatment planning for rectal cancer.

PURPOSE: To develop an automated workflow for rectal cancer three-dimensional conformal radiotherapy (3DCRT) treatment planning that combines deep learning (DL) aperture predictions and forward-planning algorithms. METHODS: We designed an algorithm to automate the clinical workflow for 3DCRT planning with field aperture creations and field-in-field (FIF) planning. DL models (DeepLabV3+ architecture) were trained, validated, and tested on 555 patients to automatically generate aperture shapes for primary (posterior-anterior [PA] and opposed laterals) and boost fields. Network inputs were digitally reconstructed radiographs, gross tumor volume (GTV), and nodal GTV. A physician scored each aperture for 20 patients on a 5-point scale (>3 is acceptable). A planning algorithm was then developed to create a homogeneous dose using a combination of wedges and subfields. The algorithm iteratively identifies a hotspot volume, creates a subfield, calculates dose, and optimizes beam weight all without user intervention. The algorithm was tested on 20 patients using clinical apertures with varying wedge angles and definitions of hotspots, and the resulting plans were scored by a physician. The end-to-end workflow was tested and scored by a physician on another 39 patients. RESULTS: The predicted apertures had Dice scores of 0.95, 0.94, and 0.90 for PA, laterals, and boost fields, respectively. Overall, 100%, 95%, and 87.5% of the PA, laterals, and boost apertures were scored as clinically acceptable, respectively. At least one auto-plan was clinically acceptable for all patients. Wedged and non-wedged plans were clinically acceptable for 85% and 50% of patients, respectively. The hotspot dose percentage was reduced from 121% (σ = 14%) to 109% (σ = 5%) of prescription dose for all plans. The integrated end-to-end workflow of automatically generated apertures and optimized FIF planning gave clinically acceptable plans for 38/39 (97%) of patients. CONCLUSION: We have successfully automated the clinical workflow for generating radiotherapy plans for rectal cancer for our institution.

Clinical acceptability of fully automated external beam radiotherapy for cervical cancer with three different beam delivery techniques.

PURPOSE: To fully automate CT-based cervical cancer radiotherapy by automating contouring and planning for three different treatment techniques. METHODS: We automated three different radiotherapy planning techniques for locally advanced cervical cancer: 2D 4-field-box (4-field-box), 3D conformal radiotherapy (3D-CRT), and volumetric modulated arc therapy (VMAT). These auto-planning algorithms were combined with a previously developed auto-contouring system. To improve the quality of the 4-field-box and 3D-CRT plans, we used an in-house, field-in-field (FIF) automation program. Thirty-five plans were generated for each technique on CT scans from multiple institutions and evaluated by five experienced radiation oncologists from three different countries. Every plan was reviewed by two of the five radiation oncologists and scored using a 5-point Likert scale. RESULTS: Overall, 87%, 99%, and 94% of the automatically generated plans were found to be clinically acceptable without modification for the 4-field-box, 3D-CRT, and VMAT plans, respectively. Some customizations of the FIF configuration were necessary on the basis of radiation oncologist preference. Additionally, in some cases, it was necessary to renormalize the plan after it was generated to satisfy radiation oncologist preference. CONCLUSION: Approximately, 90% of the automatically generated plans were clinically acceptable for all three planning techniques. This fully automated planning system has been implemented into the radiation planning assistant for further testing in resource-constrained radiotherapy departments in low- and middle-income countries.

Automatic contouring QA method using a deep learning-based autocontouring system.

PURPOSE: To determine the most accurate similarity metric when using an independent system to verify automatically generated contours. METHODS: A reference autocontouring system (primary system to create clinical contours) and a verification autocontouring system (secondary system to test the primary contours) were used to generate a pair of 6 female pelvic structures (UteroCervix [uterus + cervix], CTVn [nodal clinical target volume (CTV)], PAN [para-aortic lymph nodes], bladder, rectum, and kidneys) on 49 CT scans from our institution and 38 from other institutions. Additionally, clinically acceptable and unacceptable contours were manually generated using the 49 internal CT scans. Eleven similarity metrics (volumetric Dice similarity coefficient (DSC), Hausdorff distance, 95% Hausdorff distance, mean surface distance, and surface DSC with tolerances from 1 to 10 mm) were calculated between the reference and the verification autocontours, and between the manually generated and the verification autocontours. A support vector machine (SVM) was used to determine the threshold that separates clinically acceptable and unacceptable contours for each structure. The 11 metrics were investigated individually and in certain combinations. Linear, radial basis function, sigmoid, and polynomial kernels were tested using the combinations of metrics as inputs for the SVM. RESULTS: The highest contouring error detection accuracies were 0.91 for the UteroCervix, 0.90 for the CTVn, 0.89 for the PAN, 0.92 for the bladder, 0.95 for the rectum, and 0.97 for the kidneys and were achieved using surface DSCs with a thickness of 1, 2, or 3 mm. The linear kernel was the most accurate and consistent when a combination of metrics was used as an input for the SVM. However, the best model accuracy from the combinations of metrics was not better than the best model accuracy from a surface DSC as an input. CONCLUSIONS: We distinguished clinically acceptable contours from clinically unacceptable contours with an accuracy higher than 0.9 for the targets and critical structures in patients with cervical cancer; the most accurate similarity metric was surface DSC with a thickness of 1, 2, or 3 mm.

Clinical implementation of automated treatment planning for whole-brain radiotherapy.

The purpose of the study was to develop and clinically deploy an automated, deep learning-based approach to treatment planning for whole-brain radiotherapy (WBRT). We collected CT images and radiotherapy treatment plans to automate a beam aperture definition from 520 patients who received WBRT. These patients were split into training (n = 312), cross-validation (n = 104), and test (n = 104) sets which were used to train and evaluate a deep learning model. The DeepLabV3+ architecture was trained to automatically define the beam apertures on lateral-opposed fields using digitally reconstructed radiographs (DRRs). For the beam aperture evaluation, 1st quantitative analysis was completed using a test set before clinical deployment and 2nd quantitative analysis was conducted 90 days after clinical deployment. The mean surface distance and the Hausdorff distances were compared in the anterior-inferior edge between the clinically used and the predicted fields. Clinically used plans and deep-learning generated plans were evaluated by various dose-volume histogram metrics of brain, cribriform plate, and lens. The 1st quantitative analysis showed that the average mean surface distance and Hausdorff distance were 7.1 mm (±3.8 mm) and 11.2 mm (±5.2 mm), respectively, in the anterior-inferior edge of the field. The retrospective dosimetric comparison showed that brain dose coverage (D99%, D95%, D1%) of the automatically generated plans was 29.7, 30.3, and 32.5 Gy, respectively, and the average dose of both lenses was up to 19.0% lower when compared to the clinically used plans. Following the clinical deployment, the 2nd quantitative analysis showed that the average mean surface distance and Hausdorff distance between the predicted and clinically used fields were 2.6 mm (±3.2 mm) and 4.5 mm (±5.6 mm), respectively. In conclusion, the automated patient-specific treatment planning solution for WBRT was implemented in our clinic. The predicted fields appeared consistent with clinically used fields and the predicted plans were dosimetrically comparable.

Technical Note: Dose prediction for head and neck radiotherapy using a three-dimensional dense dilated U-net architecture.

PURPOSE: Radiation therapy treatment planning is a time-consuming and iterative manual process. Consequently, plan quality varies greatly between and within institutions. Artificial intelligence shows great promise in improving plan quality and reducing planning times. This technical note describes our participation in the American Association of Physicists in Medicine Open Knowledge-Based Planning Challenge (OpenKBP), a competition to accurately predict radiation therapy dose distributions. METHODS: A three-dimensional (3D) densely connected U-Net with dilated convolutions was developed to predict 3D dose distributions given contoured CT images of head and neck patients as input. While traditional augmentation techniques such as rotations and translations were explored, it was found that training on random patches alone resulted in the greatest model performance. A custom-weighted mean squared error loss function was employed. Finally, an ensemble of best-performing networks was used to generate the final challenge predictions. RESULTS: Our team (SuperPod) placed second in the dose stream of the OpenKBP challenge. The average mean absolute difference between the predicted and clinical dose distributions of the testing dataset was 2.56 Gy. On average, the predicted normalized target DVH metrics were within 3% of the clinical plans, and the predicted organ at risk DVH metrics were within 2 Gy of the clinical plans. CONCLUSIONS: The developed 3D dense dilated U-Net architecture can accurately predict 3D radiotherapy dose distributions and can be used as part of a fully automated radiation therapy planning pipeline.

Development and dosimetric assessment of an automatic dental artifact classification tool to guide artifact management techniques in a fully automated treatment planning workflow.

PURPOSE: We conducted our study to develop a tool capable of automatically detecting dental artifacts in a CT scan on a slice-by-slice basis and to assess the dosimetric impact of implementing the tool into the Radiation Planning Assistant (RPA), a web-based platform designed to fully automate the radiation therapy treatment planning process. METHODS: We developed an automatic dental artifact identification tool and assessed the dosimetric impact of its use in the RPA. Three users manually annotated 83,676 head-and-neck (HN) CT slices (549 patients). Majority-voting was applied to the individual annotations to determine the presence or absence of dental artifacts. The patients were divided into train, cross-validation, and test data sets (ratio: 3:1:1, respectively). A random subset of images without dental artifacts was used to balance classes (1:1) in the training data set. The Inception-V3 deep learning model was trained with the binary cross-entropy loss function. With use of this model, we automatically identified artifacts on 15 RPA HN plans on a slice-by-slice basis and investigated three dental artifact management methods applied before and after volumetric modulated arc therapy (VMAT) plan optimization. The resulting dose distributions and target coverage were quantified. RESULTS: Per-slice accuracy, sensitivity, and specificity were 99 %, 91 %, and 99 %, respectively. The model identified all patients with artifacts. Small dosimetric differences in total plan dose were observed between the various density-override methods (±1 Gy). For the pre- and post-optimized plans, 90 % and 99 %, respectively, of dose comparisons resulted in normal structure dose differences of ±1 Gy. Differences in the volume of structures receiving 95 % of the prescribed dose (V95[%]) were ≤0.25 % for 100 % of plans. CONCLUSION: The dosimetric impact of applying dental artifact management before and after artifact plan optimization was minor. Our results suggest that not accounting for dental artifacts in the current RPA workflow (where only post-optimization dental artifact management is possible) may result in minor dosimetric differences. If RPA users choose to override CT densities as a solution to managing dental artifacts, our results suggest segmenting the volume of the artifact and overriding its density to water is a safe option.

Impact of slice thickness, pixel size, and CT dose on the performance of automatic contouring algorithms.

PURPOSE: To investigate the impact of computed tomography (CT) image acquisition and reconstruction parameters, including slice thickness, pixel size, and dose, on automatic contouring algorithms. METHODS: Eleven scans from patients with head-and-neck cancer were reconstructed with varying slice thicknesses and pixel sizes. CT dose was varied by adding noise using low-dose simulation software. The impact of these imaging parameters on two in-house auto-contouring algorithms, one convolutional neural network (CNN)-based and one multiatlas-based system (MACS) was investigated for 183 reconstructed scans. For each algorithm, auto-contours for organs-at-risk were compared with auto-contours from scans with 3 mm slice thickness, 0.977 mm pixel size, and 100% CT dose using Dice similarity coefficient (DSC), Hausdorff distance (HD), and mean surface distance (MSD). RESULTS: Increasing the slice thickness from baseline value of 3 mm gave a progressive reduction in DSC and an increase in HD and MSD on average for all structures. Reducing the CT dose only had a relatively minimal effect on DSC and HD. The rate of change with respect to dose for both auto-contouring methods is approximately 0. Changes in pixel size had a small effect on DSC and HD for CNN-based auto-contouring with differences in DSC being within 0.07. Small structures had larger deviations from the baseline values than large structures for DSC. The relative differences in HD and MSD between the large and small structures were small. CONCLUSIONS: Auto-contours can deviate substantially with changes in CT acquisition and reconstruction parameters, especially slice thickness and pixel size. The CNN was less sensitive to changes in pixel size, and dose levels than the MACS. The results contraindicated more restrictive values for the parameters should be used than a typical imaging protocol for head-and-neck.

Training deep-learning segmentation models from severely limited data.

PURPOSE: To enable generation of high-quality deep learning segmentation models from severely limited contoured cases (e.g., ~10 cases). METHODS: Thirty head and neck computed tomography (CT) scans with well-defined contours were deformably registered to 200 CT scans of the same anatomic site without contours. Acquired deformation vector fields were used to train a principal component analysis (PCA) model for each of the 30 contoured CT scans by capturing the mean deformation and most prominent variations. Each PCA model can produce an infinite number of synthetic CT scans and corresponding contours by applying random deformations. We used 300, 600, 1000, and 2000 synthetic CT scans and contours generated from one PCA model to train V-Net, a 3D convolutional neural network architecture, to segment parotid and submandibular glands. We repeated the training using same numbers of training cases generated from 7, 10, 20, and 30 PCA models, with the data distributed evenly between each PCA model. Performance of the segmentation models was evaluated with Dice similarity coefficients between auto-generated contours and physician-drawn contours on 162 test CT scans for parotid glands and another 21 test CT scans for submandibular glands. RESULTS: Dice values varied with the number of synthetic CT scans and the number of PCA models used to train the network. By using 2000 synthetic CT scans generated from 10 PCA models, we achieved Dice values of 82.8% ± 6.8% for right parotid, 82.0% ± 6.9% for left parotid, and 74.2% ± 6.8% for submandibular glands. These results are comparable with those obtained from state-of-the-art auto-contouring approaches, including a deep learning network trained from more than 1000 contoured patients and a multi-atlas algorithm from 12 well-contoured atlases. Improvement was marginal when >10 PCA models or >2000 synthetic CT scans were used. CONCLUSIONS: We demonstrated an effective data augmentation approach to train high-quality deep learning segmentation models from a limited number of well-contoured patient cases.