Development of automated region of interest selection algorithms for surface-guided radiation therapy of breast cancer.

To investigate automation of the preparation of the region of interest (ROI) for surface-guided radiotherapy (SGRT) of the whole breast with two algorithms based on contour anatomies: using the body contour, and using the breast contour. The patient dataset used for modeling consisted of 39 breast cancer patients previously treated with SGRT. The patient's anatomical structures (body and ipsilateral breast) were retrieved from the planning system, and the clinical ROI (cROI) drawn by the planners was retrieved from the SGRT system for comparison. For the body-contour-based algorithm, a convolutional neural network (MobileNet-v2) was utilized to train a synthetic human model dataset to predict body joint locations. With the body joint location knowledge, an automated ROI (aROIbody ) can be created based on: (1) the superior-inferior (S-I) borders defined by the joint locations, (2) the left-right (L-R) borders defined with 3/4 of chest width, and (3) a curation of the ROI to avoid the ipsilateral armpit. For the breast-contour-based algorithm, an aROIbreast was created by first defining the ROI in the S-I direction with the ipsilateral breast boundaries. Other steps are the same as with the body-contour-based algorithm. Among the 39 patients, 24 patients were used to fine-tune the algorithm parameters, and the remaining 15 patients were used to evaluate the quality of the aROIs against the cROIs. A blinded evaluation was performed by three SGRT expert physicists to rate the acceptability and the quality (1-10 scale) of the aROIs and cROIs, and the dice similarity coefficient (DSC) was also calculated to compare the similarity between the aROIs and cROIs. The results showed that the average acceptability was 14/15 (range: 13/15-15/15) for cROIs, 13.3/15 (range: 13/15-14/15) for aROIbody , and 14.6/15 (range: 14/15-15/15) for aROIbreast . The average quality was 7.4 ± 0.8 for cROIs, 8.1 ± 1.2 for aROIbody , and 8.2 ± 0.9 for aROIbreast . The DSC with cROIs was 0.81 ± 0.06 for aROIbody , and 0.83 ± 0.04 for aROIbreast . The ROI creation time was ∼120 s for clinical, 1.3 s for aROIbody , and 1.2 s for aROIbreast . The proposed automated algorithms can improve the ROI compliance with the SGRT protocol, with a shortened preparation time. It is ready to be integrated into the clinical workflow for automated ROI preparation.

Reproducibility of chestwall and heart position using surface-guided versus RPM-guided DIBH radiotherapy for left breast cancer.

This study compared the reproducibility of chestwall and heart position using surface-guided versus RPM (real-time position management)-guided deep inspiration breath hold (DIBH) radiotherapy for left sided breast cancer. Forty DIBH patients under either surface-guided radiotherapy (SGRT) or RPM guidance were studied. For patients treated with tangential fields, reproducibility was measured as the displacements in central lung distance (CLD) and heart shadow to field edge distance (HFD) between pretreatment MV (megavoltage) images and planning DRRs (digitally reconstructed radiographs). For patients treated with volumetric modulated arc therapy (VMAT), sternum to isocenter (ISO) distance (StID), spine to rib edge distance (SpRD), and heart shadow to central axis (CAX) distance (HCD) between pretreatment kV images and planning DRRs were measured. These displacements were compared between SGRT and RPM-guided DIBH. In tangential patients, the mean absolute displacements of SGRT versus RPM guidance were 0.19 versus 0.23 cm in CLD, and 0.33 versus 0.62 cm in HFD. With respect to planning DRR, heart appeared closer to the field edge by 0.04 cm with surface imaging versus 0.62 cm with RPM. In VMAT patients, the displacements of surface imaging versus RPM guidance were 0.21 versus 0.15 cm in StID, 0.24 versus 0.19 cm in SpRD, and 0.72 versus 0.41 cm in HCD. Heart appeared 0.41 cm further away from CAX with surface imaging, whereas 0.10 cm closer to field CAX with RPM. None of the differences between surface imaging and RPM guidance was statistically significant. In conclusion, the displacements of chestwall were small and were comparable with SGRT- or RPM-guided DIBH. The position deviations of heart were larger than those of chestwall with SGRT or RPM. Although none of the differences between SGRT and RPM guidance were statistically significant, there was a trend that the position deviations of heart were smaller and more favorable with SGRT than with RPM guidance in tangential patients.

Comparison of setup accuracy of optical surface image versus orthogonal x-ray images for VMAT of the left breast using deep-inspiration breath-hold.

To compare the setup accuracy of optical surface image (OSI) versus orthogonal x-ray images (2DkV) using cone beam computed tomography (CBCT) as ground truth for radiotherapy of left breast cancer in deep-inspiration breath-hold (DIBH). Ten left breast DIBH patients treated with volumetric modulated arc therapy (VMAT) were studied retrospectively. OSI, 2DkV, and CBCT were acquired weekly at treatment setup. OSI, 2DkV, and CBCT were registered to planning CT or planning DRR based on a breast surface region of interest (ROI), bony anatomy (chestwall and sternum), and both bony anatomy and breast surface, respectively. These registrations provided couch shifts for each imaging system. The setup errors, or the difference in couch shifts between OSI and CBCT were compared to those between 2DkV and CBCT. A second OSI was acquired during last beam delivery to evaluate intrafraction motion. The median absolute setup errors were (0.21, 0.27, 0.23 cm, 0.6°, 1.3°, 1.0°) for OSI, and (0.26, 0.24, 0.18 cm, 0.9°, 1.0°, 0.6°) for 2DkV in vertical, longitudinal and lateral translations, and in rotation, roll and pitch, respectively. None of the setup errors was significantly different between OSI and 2DkV. For both systems, the systematic and random setup errors were ≤0.6 cm and ≤1.5° in all directions. Nevertheless, larger setup errors were observed in some sessions in both systems. There was no correlation between OSI and CBCT whereas there was modest correlation between 2DkV and CBCT. The intrafraction motion in DIBH detected by OSI was small with median absolute translations <0.2 cm, and rotations ≤0.4°. Though OSI showed comparable and small setup errors as 2DkV, it showed no correlation with CBCT. We concluded that to achieve accurate setup for both bony anatomy and breast surface, daily 2DkV can't be omitted following OSI for left breast patients treated with DIBH VMAT.

Automated tracking of morphologic changes in weekly magnetic resonance imaging during head and neck radiotherapy.

BACKGROUND AND PURPOSE: Anatomic changes during head and neck radiotherapy can impact dose delivery, necessitate adaptive replanning, and indicate patient-specific response to treatment. We have developed an automated system to track these changes through longitudinal MRI scans to aid identification and clinical intervention. The purpose of this article is to describe this tracking system and present results from an initial cohort of patients. MATERIALS AND METHODS: The Automated Watchdog in Adaptive Radiotherapy Environment (AWARE) was developed to process longitudinal MRI data for radiotherapy patients. AWARE automatically identifies and collects weekly scans, propagates radiotherapy planning structures, computes structure changes over time, and reports important trends to the clinical team. AWARE also incorporates manual structure review and revision from clinical experts and dynamically updates tracking statistics when necessary. AWARE was applied to patients receiving weekly T2-weighted MRI scans during head and neck radiotherapy. Changes in nodal gross tumor volume (GTV) and parotid gland delineations were tracked over time to assess changes during treatment and identify early indicators of treatment response. RESULTS: N = 91 patients were tracked and analyzed in this study. Nodal GTVs and parotids both shrunk considerably throughout treatment (-9.7 ± 7.7% and -3.7 ± 3.3% per week, respectively). Ipsilateral parotids shrunk significantly faster than contralateral (-4.3 ± 3.1% vs. -2.9 ± 3.3% per week, p = 0.005) and increased in distance from GTVs over time (+2.7 ± 7.2% per week, p < 1 × 10-5 ). Automatic structure propagations agreed well with manual revisions (Dice = 0.88 ± 0.09 for parotids and 0.80 ± 0.15 for GTVs), but for GTVs the agreement degraded 4-5 weeks after the start of treatment. Changes in GTV volume observed by AWARE as early as one week into treatment were predictive of large changes later in the course (AUC = 0.79). CONCLUSION: AWARE automatically identified longitudinal changes in GTV and parotid volumes during radiotherapy. Results suggest that this system may be useful for identifying rapidly responding patients as early as one week into treatment.

Intrafractional accuracy and efficiency of a surface imaging system for deep inspiration breath hold during ablative gastrointestinal cancer treatment.

PURPOSE: Beam gating with deep inspiration breath hold (DIBH) usually depends on some external surrogate to infer internal target movement, and the exact internal movement is unknown. In this study, we tracked internal targets and characterized residual motion during DIBH treatment, guided by a surface imaging system, for gastrointestinal cancer. We also report statistics on treatment time. METHODS AND MATERIALS: We included 14 gastrointestinal cancer patients treated with surface imaging-guided DIBH volumetrically modulated arc therapy, each with at least one radiopaque marker implanted near or within the target. They were treated in 25, 15, or 10 fractions. Thirteen patients received treatment for pancreatic cancer, and one underwent separate treatments for two liver metastases. The surface imaging system monitored a three-dimensional surface with ± 3 mm translation and ± 3° rotation threshold. During delivery, a kilovolt image was automatically taken every 20° or 40° gantry rotation, and the internal marker was identified from the image. The displacement and residual motion of the markers were calculated. To analyze the treatment efficiency, the treatment time of each fraction was obtained from the imaging and treatment timestamps in the record and verify system. RESULTS: Although the external surface was monitored and limited to ± 3 mm and ± 3°, significant residual internal target movement was observed in some patients. The range of residual motion was 3-21 mm. The average displacement for this cohort was 0-3 mm. In 19% of the analyzed images, the magnitude of the instantaneous displacement was > 5 mm. The mean treatment time was 17 min with a standard deviation of 4 min. CONCLUSIONS: Precaution is needed when applying surface image guidance for gastrointestinal cancer treatment. Using it as a solo DIBH technique is discouraged when the correlation between internal anatomy and patient surface is limited. Real-time radiographic verification is critical for safe treatments.

Intrafractional relationship changes between an external breathing signal and fiducial marker positions in pancreatic cancer patients.

BACKGROUND AND PURPOSE: The purpose of this study of pancreatic cancer patients treated with respiratory-guided stereotactic body radiotherapy (SBRT) on a standard linac was to investigate (a) the intrafractional relationship change (IRC) between a breathing signal and the tumor position, (b) the impact of IRC on the delivered dose, and (c) potential IRC predictors. MATERIALS AND METHODS: We retrospectively investigated 10 pancreatic cancer patients with 2-4 implanted fiducial markers in the tumor treated with SBRT. Fluoroscopic images were acquired before and after treatment delivery simultaneously with the abdominal breathing motion. We quantified the IRC as the change in fiducial location for a given breathing amplitude in the left-right (LR), anterior-posterior (AP), and superior-inferior (SI) directions from before to after treatment delivery. The treatment plans were re-calculated after changing the isocenter coordinates according to the IRCs. Four treatment- or patient-related factors were investigated as potential predictors for IRC using linear models. RESULTS: The average (±1 SD) absolute IRCs in the LR, AP, and SI directions were 1.2 ± 1.2 mm, 0.7 ± 0.7 mm, and 1.1 ± 0.8 mm, respectively. The average 3D IRC was 2.0 ± 1.3 mm (range: 0.4-5.3 mm) for a median treatment delivery time of 8.5 min (range: 5.7-19.9 min; n = 31 fractions). The dose coverage of the internal target volume (ITV) decreased by more than 3% points in three of 31 fractions. In those cases, the 3D IRC had been larger than 4.3 mm. The 3D IRC was found to correlate with changes in the minimum breathing amplitude during treatment delivery. CONCLUSION: On average, 2 mm of treatment delivery accuracy was lost due to IRC. Periodical intrafractional imaging is needed to safely deliver respiratory-guided SBRT.

Assessment of the use of different imaging and delivery techniques for cranial treatments on the Halcyon linac.

PURPOSE: In this work, we investigated the effect on the workflow and setup accuracy of using surface guided radiation therapy (SGRT) for patient setup, megavoltage cone beam CT (MVCBCT) or kilovoltage cone beam CT (kVCBCT) for imaging and fixed IMRT or volumetric-modulated arc therapy (VMAT) for treatment delivery with the Halcyon linac. METHODS: We performed a retrospective investigation of 272 treatment fractions, using three different workflows. The first and second workflows used MVCBCT and fixed IMRT for imaging and treatment delivery, and the second one also used SGRT for patient setup. The third workflow used SGRT for setup, kVCBCT for imaging and VMAT for delivery. Workflows were evaluated by comparing the number of fractions requiring repeated imaging acquisitions and the time required for setup, imaging and treatment delivery. Setup position accuracy was assessed by comparing the daily kV- or MV- CBCT with the planning CT and measuring the residual rotational errors for pitch, yaw and roll angles. RESULTS: Without the use of SGRT, the imaging fields were delivered more than once on 11.1% of the fractions, while re-imaging was necessary in 5.5% of the fractions using SGRT. The total treatment time, including setup, imaging, and delivery, for the three workflows was 531 ± 157 s, 503 ± 130 s and 457 ± 91 s, respectively. A statistically significant difference was observed when comparing the third workflow with the first two. The total residual rotational errors were 1.96 ± 1.29°, 1.28 ± 0.67° and 1.22 ± 0.76° and statistically significant differences were observed when comparing workflows with and without SGRT. CONCLUSIONS: The use of SGRT allowed for a reduction of re-imaging during patient setup and improved patient position accuracy by reducing residual rotational errors. A reduction in treatment time using kVCBCT with SGRT was observed. The most efficient workflow was the one including kVCBCT and SGRT for setup and VMAT for delivery.

A survey of surface imaging use in radiation oncology in the United States.

Surface imaging (SI) has been rapidly integrated into radiotherapy clinics across the country without specific guidelines and recommendations on its commissioning and use aside from vendor-provided information. A survey was created under the auspices of AAPM TG-302 to assess the current status of SI to identify if there is need for formal guidance. The survey was designed to determine the institutional setting of responders, availability and length of its use, commissioning procedures, and clinical applications. This survey was created in REDCap, and approved as IRB exempt to collect anonymized data. Questions were reviewed by multiple physicists to ensure concept validity and piloted by a small group of independent physicists to ensure process validity. All full members of AAPM self-identified as "therapy" or "other" were sent the survey link by email. The survey was active from February to March 2018. Of 3677 members successfully contacted, 439 completed responses; the summary of these responses provides insight on current surface imaging clinical practices, though they should not be assumed to be representative of radiation oncology as a whole. Results showed that 53.3% of respondents have SI in their clinics, mostly in treatment rooms, rarely in simulation rooms. Half of those without SI plan on purchasing it within 3 years. Over 10% have SI but do not use it clinically, 36.8% classify themselves as "expert" users, and 85.5% agreed/strongly agreed that SI guidelines are needed. Initial positioning with SI is most common for breast/chestwall and SRS/SBRT treatments, least common for pediatrics. Use of SI for intra-fraction monitoring follows a similar distribution. Gating with SI is most prevalent for breast/chestwall (66.0%) but also used in SBRT (33.0%), and non-SBRT lung/abdomen (<30%) treatments. SI is a rapidly growing technology in the field with widespread use for several anatomic sites. Guidelines and recommendations on commissioning and clinical use are warranted.

Dosimetric study of the plan quality and dose to organs at risk on tangential breast treatments using the Halcyon linac.

PURPOSE: To investigate the plan quality and doses to the heart, contralateral breast (CB), ipsilateral lung (IL), and contralateral lung (CL) in tangential breast treatments using the Halcyon linac with megavoltage setup fields. METHODS: Radiotherapy treatment plans with tangential beams from 25 breast cancer patients previously treated on a C-arm linac were replanned for Halcyon. Thirteen corresponded to right-sided breasts and 12 to left-sided breasts, all with a dose prescription of 50 Gy in 25 fractions. Plans were created with the following setup imaging techniques: low-dose (LD) MVCBCT, high-quality (HQ) MVCBCT, LD-MV and HQ-MV pairs and the imaging dose was included in the plans. Plan quality metric values for the lumpectomy cavity, whole-breast and doses to the organs at risk (OARs) were measured and compared with those from the original plans. RESULTS: No significant differences in plan quality were observed between the original and Halcyon plans. An increase in the mean dose (Mean) for all the organs was observed for the Halcyon plans. For right-sided plans, the accumulated Mean over the 25 fractions in the C-arm plans was 0.4 ± 0.3, 0.2 ± 0.2, 5.4 ± 1.3, and 0.1 ± 0.1 Gy for the heart, CB, IL, and CL, respectively, while values in the MVCBCT-LD Halcyon plans were 1.2 ± 0.2, 0.6 ± 0.1, 6.5 ± 1.4, and 0.4 ± 0.1 Gy, respectively. For left-sided treatments, Mean in the original plans was 0.9 ± 0.2, 0.1 ± 0.0, 4.2 ± 1.2, and 0.0 ± 0.0 Gy, while for the MVCBCT-LD Halcyon plans values were 1.9 ± 0.2, 0.6 ± 0.2, 5.1 ± 1.2, and 0.5 ± 0.2 Gy, respectively. CONCLUSIONS: Plan quality for breast treatments using Halcyon is similar to the quality for a 6 MV, C-arm plan. For treatments using megavoltage setup fields, the dose contribution to OARs from the imaging fields can be equal or higher than the dose from treatment fields.

Automatic patient positioning and gating window settings in respiratory-gated stereotactic body radiation therapy for pancreatic cancer using fluoroscopic imaging.

Before treatment delivery of respiratory-gated radiation therapy (RT) in patients with implanted fiducials, both the patient position and the gating window thresholds must be set. In linac-based RT, this is currently done manually and setup accuracy will therefore be dependent on the skill of the user. In this study, we present an automatic method for finding the patient position and the gating window thresholds. Our method uses sequentially acquired anterior-posterior (AP) and lateral fluoroscopic imaging with simultaneous breathing amplitude monitoring and intends to reach 100% gating accuracy while keeping the duty cycle as high as possible. We retrospectively compared clinically used setups to the automatic setups by our method in five pancreatic cancer patients treated with hypofractionated RT. In 15 investigated fractions, the average (±standard deviation) differences between the clinical and automatic setups were -0.4 ± 0.8 mm, -1.0 ± 1.1 mm, and 1.8 ± 1.3 mm in the left-right (LR), the AP, and the superior-inferior (SI) direction, respectively. For the clinical setups, typical interfractional setup variations were 1-2 mm in the LR and AP directions, and 2-3 mm in the SI direction. Using the automatic method, the duty cycle could be improved in six fractions, in four fractions the duty cycle had to be lowered to improve gating accuracy, and in five fractions both duty cycle and gating accuracy could be improved. Our automatic method has the potential to increase accuracy and decrease user dependence of setup for patients with implanted fiducials treated with respiratory-gated RT. After fluoroscopic image acquisition, the calculated patient shifts and gating window thresholds are calculated in 1-2 s. The method gives the user the possibility to evaluate the effect of different patient positions and gating window thresholds on gating accuracy and duty cycle. If deemed necessary, it can be used at any time during treatment delivery.

Clinical efficacy and safety of surface imaging guided radiosurgery (SIG-RS) in the treatment of benign skull base tumors.

Frameless, surface imaging guided radiosurgery (SIG-RS) is a novel platform for stereotactic radiosurgery (SRS) wherein patient positioning is monitored in real-time through infra-red camera tracking of facial topography. Here we describe our initial clinical experience with SIG-RS for the treatment of benign neoplasms of the skull base. We identified 48 patients with benign skull base tumors consecutively treated with SIG-RS at a single institution between 2009 and 2011. Patients were diagnosed with meningioma (n = 22), vestibular schwannoma (n = 20), or nonfunctional pituitary adenoma (n = 6). Local control and treatment-related toxicity were retrospectively assessed. Median follow-up was 65 months (range 61-72 months). Prescription doses were 12-13 Gy in a single fraction (n = 18), 8 Gy × 3 fractions (n = 6), and 5 Gy × 5 fractions (n = 24). Actuarial tumor control rate at 5 years was 98%. No grade ≥3 treatment-related toxicity was observed. Grade ≤2 toxicity was associated with symptomatic lesions (p = 0.049) and single fraction treatment (p = 0.005). SIG-RS for benign skull base tumors produces clinical outcomes comparable to conventional frame-based SRS techniques while enhancing patient comfort.

Dosimetric assessment of brass mesh bolus for postmastectomy photon radiotherapy.

Brass mesh bolus has been shown to be an acceptable substitute for tissue-equivalent bolus to increase superficial dose for chest wall tangent photon radiotherapy. This work investigated the increase in surface dose, the change in the dose at depth, and the safety implications of higher energy photon beams when using brass mesh bolus for postmastectomy chest wall radiotherapy. A photon tangent plan was delivered to a thorax phantom, and the superficial dose ranged from 40%-72% of prescription dose with no bolus. The surface dose increased to 75%-110% of prescription dose with brass mesh bolus and 85%-109% of prescription dose with tissue-equivalent bolus. It was also found that the dose at depth when using brass mesh bolus is comparable to that measured with no bolus for en face and oblique incidence. Monte Carlo calculations were used to assess the photoneutron production from brass mesh bolus used with 15 MV and 24 MV photon beams. The effective dose from photoneutrons was approximated and found to be relatively small, yet not negligible. Activation products generated by these photoneutrons, the surface dose rate due to the activation products, and the half-life of the activation products were also considered in this work. The authors conclude that brass mesh bolus is a reasonable alternative to tissue-equivalent bolus, and it may be used with high-energy beam; but one should be aware of the potential increased effective dose to staff and patients due to the activation products produced by photoneutrons.

SPARSE: Seed Point Auto-Generation for Random Walks Segmentation Enhancement in medical inhomogeneous targets delineation of morphological MR and CT images.

In medical image processing, robust segmentation of inhomogeneous targets is a challenging problem. Because of the complexity and diversity in medical images, the commonly used semiautomatic segmentation algorithms usually fail in the segmentation of inhomogeneous objects. In this study, we propose a novel algorithm imbedded with a seed point autogeneration for random walks segmentation enhancement, namely SPARSE, for better segmentation of inhomogeneous objects. With a few user-labeled points, SPARSE is able to generate extended seed points by estimating the probability of each voxel with respect to the labels. The random walks algorithm is then applied upon the extended seed points to achieve improved segmentation result. SPARSE is implemented under the compute unified device architecture (CUDA) programming environment on graphic processing unit (GPU) hardware platform. Quantitative evaluations are performed using clinical homogeneous and inhomogeneous cases. It is found that the SPARSE can greatly decrease the sensitiveness to initial seed points in terms of location and quantity, as well as the freedom of selecting parameters in edge weighting function. The evaluation results of SPARSE also demonstrate substantial improvements in accuracy and robustness to inhomogeneous target segmentation over the original random walks algorithm.

A method for generating large datasets of organ geometries for radiotherapy treatment planning studies.

BACKGROUND: With the rapidly increasing application of adaptive radiotherapy, large datasets of organ geometries based on the patient's anatomy are desired to support clinical application or research work, such as image segmentation, re-planning, and organ deformation analysis. Sometimes only limited datasets are available in clinical practice. In this study, we propose a new method to generate large datasets of organ geometries to be utilized in adaptive radiotherapy. METHODS: Given a training dataset of organ shapes derived from daily cone-beam CT, we align them into a common coordinate frame and select one of the training surfaces as reference surface. A statistical shape model of organs was constructed, based on the establishment of point correspondence between surfaces and non-uniform rational B-spline (NURBS) representation. A principal component analysis is performed on the sampled surface points to capture the major variation modes of each organ. RESULTS: A set of principal components and their respective coefficients, which represent organ surface deformation, were obtained, and a statistical analysis of the coefficients was performed. New sets of statistically equivalent coefficients can be constructed and assigned to the principal components, resulting in a larger geometry dataset for the patient's organs. CONCLUSIONS: These generated organ geometries are realistic and statistically representative.

Creating patient-specific digital phantoms with a longitudinal atlas for evaluating deformable CT-CBCT registration in adaptive lung radiotherapy.

BACKGROUND: Quality assurance of deformable image registration (DIR) is challenging because the ground truth is often unavailable. In addition, current approaches that rely on artificial transformations do not adequately resemble clinical scenarios encountered in adaptive radiotherapy. PURPOSE: We developed an atlas-based method to create a variety of patient-specific serial digital phantoms with CBCT-like image quality to assess the DIR performance for longitudinal CBCT imaging data in adaptive lung radiotherapy. METHODS: A library of deformations was created by extracting the longitudinal changes observed between a planning CT and weekly CBCT from an atlas of lung radiotherapy patients. The planning CT of an inquiry patient was first deformed by mapping the deformation pattern from a matched atlas patient, and subsequently appended with CBCT artifacts to imitate a weekly CBCT. Finally, a group of digital phantoms around an inquiry patient was produced to simulate a series of possible evolutions of tumor and adjacent normal structures. We validated the generated deformation vector fields (DVFs) to ensure numerically and physiologically realistic transformations. The proposed framework was applied to evaluate the performance of the DIR algorithm implemented in the commercial Eclipse treatment planning system in a retrospective study of eight inquiry patients. RESULTS: The generated DVFs were inverse consistent within less than 3 mm and did not exhibit unrealistic folding. The deformation patterns adequately mimicked the observed longitudinal anatomical changes of the matched atlas patients. Worse Eclipse DVF accuracy was observed in regions of low image contrast or artifacts. The structure volumes exhibiting a DVF error magnitude of equal or more than 2 mm ranged from 24.5% (spinal cord) to 69.2% (heart) and the maximum DVF error exceeded 5 mm for all structures except the spinal cord. Contour-based evaluations showed a high degree of alignment with dice similarity coefficients above 0.8 in all cases, which underestimated the overall DVF accuracy within the structures. CONCLUSIONS: It is feasible to create and augment digital phantoms based on a particular patient of interest using multiple series of deformation patterns from matched patients in an atlas. This can provide a semi-automated procedure to complement the quality assurance of CT-CBCT DIR and facilitate the clinical implementation of image-guided and adaptive radiotherapy that involve longitudinal CBCT imaging studies.

Toward quantitative intrafractional monitoring in paraspinal SBRT using a proprietary software application: clinical implementation and patient results.

Objective.We report on paraspinal motion and the clinical implementation of our proprietary software that leverages Varian's intrafraction motion review (IMR) capability for quantitative tracking of the spine during paraspinal SBRT. The work is based on our prior development and analysis on phantoms.Approach.To address complexities in patient anatomy, digitally reconstructed radiographs (DRR's) that highlight only the spine or hardware were constructed as tracking reference. Moreover, a high-pass filter and first-pass coarse search were implemented to enhance registration accuracy and stability. For evaluation, 84 paraspinal SBRT patients with sites spanning across the entire vertebral column were enrolled with prescriptions ranging from 24 to 40 Gy in one to five fractions. Treatments were planned and delivered with 9 IMRT beams roughly equally distributed posteriorly. IMR was triggered every 200 or 500 MU for each beam. During treatment, the software grabbed the IMR image, registered it with the corresponding DRR, and displayed the motion result in near real-time on auto-pilot mode. Four independent experts completed offline manual registrations as ground truth for tracking accuracy evaluation.Main results.Our software detected ≥1.5 mm and ≥2 mm motions among 17.1% and 6.6% of 1371 patient images, respectively, in either lateral or longitudinal direction. In the validation set of 637 patient images, 91.9% of the tracking errors compared to manual registration fell within ±0.5 mm in either direction. Given a motion threshold of 2 mm, the software accomplished a 98.7% specificity and a 93.9% sensitivity in deciding whether to interrupt treatment for patient re-setup.Significance.Significant intrafractional motion exists in certain paraspinal SBRT patients, supporting the need for quantitative motion monitoring during treatment. Our improved software achieves high motion tracking accuracy clinically and provides reliable guidance for treatment intervention. It offers a practical solution to ensure accurate delivery of paraspinal SBRT on a conventional Linac platform.

Automated planning of stereotactic spine re-irradiation using cumulative dose limits.

Background and Purpose: The lack of dedicated tools in commercial planning systems currently restricts efficient review and planning for re-irradiation. The aim of this study was to develop an automated re-irradiation planning framework based on cumulative doses. Materials and Methods: We performed a retrospective study of 14 patients who received spine SBRT re-irradiation near a previously irradiated treatment site. A fully-automated workflow, DART (Dose Accumulation-based Re-irradiation Tool), was implemented within Eclipse by leveraging a combination of a dose accumulation script and a proprietary automated optimization algorithm. First, we converted the prior treatment dose into equivalent dose in 2 Gy fractions (EQD2) and mapped it to the current anatomy, utilizing deformable image registration. Subsequently, the intersection of EQD2 isodose lines with relevant organs at risk defines a series of optimization structures. During plan optimization, the residual allowable dose at a specified tissue tolerance was treated as a hard constraint. Results: All DART plans met institutional physical and cumulative constraints and passed plan checks by qualified medical physicists. DART demonstrated significant improvements in target coverage over clinical plans, with an average increase in PTV D99% and V100% of 2.3 Gy [range -0.3-7.7 Gy] and 3.4 % [range -0.4 %-7.6 %] (p < 0.01, paired t-test), respectively. Moreover, high-dose spillage (>105 %) outside the PTV was reduced by up to 7 cm3. The homogeneity index for DART plans was improved by 19 % (p < 0.001). Conclusions: DART provides a powerful framework to achieve more tailored re-irradiation plans by accounting for dose distributions from the previous treatments. The superior plan quality could improve the therapeutic ratio for re-irradiation patients.

An automated technique for global noise level measurement in CT image with a conjunction of image gradient.

Automated assessment of noise level in clinical computed tomography (CT) images is a crucial technique for evaluating and ensuring the quality of these images. There are various factors that can impact CT image noise, such as statistical noise, electronic noise, structure noise, texture noise, artifact noise, etc. In this study, a method was developed to measure the global noise index (GNI) in clinical CT scans due to the fluctuation of x-ray quanta. Initially, a noise map is generated by sliding a 10 × 10 pixel for calculating Hounsfield unit (HU) standard deviation and the noise map is further combined with the gradient magnitude map. By employing Boolean operation, pixels with high gradients are excluded from the noise histogram generated with the noise map. By comparing the shape of the noise histogram from this method with Christianson's tissue-type global noise measurement algorithm, it was observed that the noise histogram computed in anthropomorphic phantoms had a similar shape with a close GNI value. In patient CT images, excluding the HU deviation due the structure change demonstrated to have consistent GNI values across the entire CT scan range with high heterogeneous tissue compared to the GNI values using Christianson's tissue-type method. The proposed GNI was evaluated in phantom scans and was found to be capable of comparing scan protocols between different scanners. The variation of GNI when using different reconstruction kernels in clinical CT images demonstrated a similar relationship between noise level and kernel sharpness as observed in uniform phantom: sharper kernel resulted in noisier images. This indicated that GNI was a suitable index for estimating the noise level in clinical CT images with either a smooth or grainy appearance. The study's results suggested that the algorithm can be effectively utilized to screen the noise level for a better CT image quality control.

Deep learning-based target decomposition for markerless lung tumor tracking in radiotherapy.

BACKGROUND: In radiotherapy, real-time tumor tracking can verify tumor position during beam delivery, guide the radiation beam to target the tumor, and reduce the chance of a geometric miss. Markerless kV x-ray image-based tumor tracking is challenging due to the low tumor visibility caused by tumor-obscuring structures. Developing a new method to enhance tumor visibility for real-time tumor tracking is essential. PURPOSE: To introduce a novel method for markerless kV image-based tracking of lung tumors via deep learning-based target decomposition. METHODS: We utilized a conditional Generative Adversarial Network (cGAN), known as Pix2Pix, to build a patient-specific model and generate the synthetic decomposed target image (sDTI) to enhance tumor visibility on the real-time kV projection images acquired by the onboard kV imager equipped on modern linear accelerators. We used 4DCT simulation images to generate the digitally reconstructed radiograph (DRR) and DTI image pairs for model training. We augmented the training dataset by randomly shifting the 4DCT in the superior-inferior, anterior-posterior, and left-right directions during the DRR and DTI generation process. We performed real-time 2D tumor tracking via template matching between the DTI generated from the CT simulation and the sDTI generated from the real-time kV projection images. We validated the proposed method using nine patients' datasets with implanted beacons near the tumor. RESULTS: The sDTI can effectively improve the image contrast around the lung tumors on the kV projection images for the nine patients. With the beacon motion as ground truth, the tracking errors were on average 0.8 ± 0.7 mm in the superior-inferior (SI) direction and 0.9 ± 0.8 mm in the in-plane left-right (IPLR) direction. The percentage of successful tracking, defined as a tracking error less than 2 mm in the SI direction, is 92.2% on the 4312 tested images. The patient-specific model took approximately 12 h to train. During testing, it took approximately 35 ms to generate one sDTI, and 13 ms to perform the tumor tracking using template matching. CONCLUSIONS: Our method offers the potential solution for nearly real-time markerless lung tumor tracking. It achieved a high level of accuracy and an impressive tracking rate. Further development of 3D lung tumor tracking is warranted.

Clinical Experience With an Offline Adaptive Radiation Therapy Head and Neck Program: Dosimetric Benefits and Opportunities for Patient Selection.

PURPOSE: The objective of this study was to develop a linear accelerator (LINAC)-based adaptive radiation therapy (ART) workflow for the head and neck that is informed by automated image tracking to identify major anatomic changes warranting adaptation. In this study, we report our initial clinical experience with the program and an investigation into potential trigger signals for ART. METHODS AND MATERIALS: Offline ART was systematically performed on patients receiving radiation therapy for head and neck cancer on C-arm LINACs. Adaptations were performed at a single time point during treatment with resimulation approximately 3 weeks into treatment. Throughout treatment, all patients were tracked using an automated image tracking system called the Automated Watchdog for Adaptive Radiotherapy Environment (AWARE). AWARE measures volumetric changes in gross tumor volumes (GTVs) and selected normal tissues via cone beam computed tomography scans and deformable registration. The benefit of ART was determined by comparing adaptive plan dosimetry and normal tissue complication probabilities against the initial plans recalculated on resimulation computed tomography scans. Dosimetric differences were then correlated with AWARE-measured volume changes to identify patient-specific triggers for ART. Candidate trigger variables were evaluated using receiver operator characteristic analysis. RESULTS: In total, 46 patients received ART in this study. Among these patients, we observed a significant decrease in dose to the submandibular glands (mean ± standard deviation: -219.2 ± 291.2 cGy, P < 10-5), parotids (-68.2 ± 197.7 cGy, P = .001), and oral cavity (-238.7 ± 206.7 cGy, P < 10-5) with the adaptive plan. Normal tissue complication probabilities for xerostomia computed from mean parotid doses also decreased significantly with the adaptive plans (P = .008). We also observed systematic intratreatment volume reductions (ΔV) for GTVs and normal tissues. Candidate triggers were identified that predicted significant improvement with ART, including parotid ΔV = 7%, neck ΔV = 2%, and nodal GTV ΔV = 29%. CONCLUSIONS: Systematic offline head and neck ART was successfully deployed on conventional LINACs and reduced doses to critical salivary structures and the oral cavity. Automated cone beam computed tomography tracking provided information regarding anatomic changes that may aid patient-specific triggering for ART.