The effect of common dental fixtures on treatment planning and delivery for head and neck intensity modulated proton therapy.

PURPOSE: Proton treatment plan perturbation by common dental fixtures such as amalgams (Am) and porcelain-fused-to-metal (PFM) crowns has, to date, been uncharacterized. Previous studies have been conducted to determine the physical effect of these materials within the beam path for single spots, but their effects on complex treatment plans and clinical anatomy have not yet been quantified. The present manuscript aims to study the effect of Am and PFM fixtures on proton treatment planning in a clinical setting. METHODS: An anthropomorphic phantom with removable tongue, maxilla, and mandible modules was simulated on a clinical computed tomography (CT) scanner. Spare maxilla modules were modified to include either a 1.5 mm depth central groove occlusal amalgam (Am) or a porcelain-fused-to-metal (PFM) crown, implanted on the first right molar. Modified tongue modules were 3D printed to accommodate several axial or sagittal oriented pieces of EBT-3 film. Clinically representative spot-scanning proton plans were generated in Eclipse v.15.6 using the proton convolution superposition (PCS) algorithm v.15.6.06 using a multi-field optimization (MFO) technique with the goal of delivering a uniform 54 Gy dose to a clinical target volume (CTV) typical of a base-of-tongue (BoT) treatment. A typical geometric beam arrangement of two anterior oblique (AO) beams and a posterior beam was employed. Plans optimized without any material overrides were delivered to the phantom A) without implants; B) with Am fixture; or C) with PFM crown. Plans were also reoptimized and delivered with inclusion of material overrides to equate relative stopping power of the fixture with that of a previously measured result. RESULTS: Plans exhibit slightly greater dose weight towards AO beams. The optimizer accounted for inclusion of fixture overrides by increasing beam weights to the beam closest to the implant. Film measurements exhibited cold spots directly within the beam path through the fixture in plans with and without overridden materials. Cold spots were somewhat mitigated in plans including overridden materials in the structure set but were not entirely eliminated. Cold spots associated with Am and PFM fixtures were quantified at 17% and 14% for plans without overrides, respectively, and 11% and 9% with using Monte Carlo simulation. Compared with film measurements and Monte Carlo simulation, the treatment planning system underestimates the dose shadowing effect in plans including material overrides. CONCLUSIONS: Dental fixtures create a dose shadowing effect directly in line with the beam path through the material. This cold spot is partially mitigated by overriding the material to measured relative stopping powers. Due to uncertainties in modeling perturbation through the fixture, the magnitude of the cold spot is underestimated using the institutional TPS when compared to measurement and MC simulation.

Comparison of MR-soft tissue based versus biliary stent based alignment for image guidance in pancreatic SBRT.

PURPOSE: The role of biliary stents in image-guided localization for pancreatic cancer has been inconclusive. To date, stent accuracy has been largely evaluated against implanted fiducials on cone beam computed tomography. We aim to use magnetic resonance (MR) soft tissue as a direct reference to examine the geometric and dosimetric impacts of stent-based localization on the newly available MR linear accelerator. METHODS: Thirty pancreatic cancer patients (132 fractions) treated on our MR linear accelerator were identified to have a biliary stent. In our standard adaptive workflow, patients were set up to the target using soft tissue for image registration and structures were re-contoured on daily MR images. The original plan was then projected on treatment anatomy and dose predicted, followed by plan re-optimization and treatment delivery. These online predicted plans were soft tissue-based and served as reference plans. Retrospective image registration to the stent was performed offline to simulate stent-based localization and the magnitude of shifts was taken as the geometric accuracy of stent localization. New predicted plans were generated based on stent-alignment for dosimetric comparison. RESULTS: Shifts were within 3 mm for 90% of the cases (mean = 1.5 mm); however, larger shifts up to 7.2 mm were observed. Average PTV coverage dropped by 1.1% with a maximum drop of 26.8%. The mean increase in V35Gy was 0.15, 0.05, 0.02, and 0.02 cc for duodenum, stomach, small bowel and large bowel, respectively. Stent alignment was significantly worse for all metrics except for small bowel (p = 0.07). CONCLUSIONS: Overall discrepancy between stent- and soft tissue-alignment was modest; however, large discrepancies were observed for select cases. While PTV coverage loss may be compensated for by using a larger margin, the increase in dose to gastrointestinal organs at risk may limit the role of biliary stents in image-guided localization.

Characterization and commissioning of a Leksell Gamma Knife ICON system for framed and frameless stereotactic radiosurgery.

PURPOSE: The Leksell Gamma Knife Icon unit (IU) was introduced recently as an upgrade to the Perfexion unit (PU) at our Gamma Knife practice. In the current study, we sought mainly to characterize dosimetry and targeting accuracy of the IU treatment deliveries using both invasive frame and frameless treatment workflows. METHODS: Relative output factors were measured by delivering single-shot 4, 8 and 16 mm radiation profiles in the manufacturer's acrylonitrile butadiene styrene spherical phantom in coronal and sagittal planes using EBT3 film. Resultant dosimetry was compared with the manufacturer's dose calculation and derived output factors were compared with the manufacturer's published value. Geometric consistency of stereotactic coordinates based on cone-beam computed tomography (CBCT) versus the traditional conventional CT-based method was characterized using a rigid phantom containing nine fiducial indicators over four separate trials. End-to-end (E2E) testing using EBT3 film was designed to evaluate both dosimetric and geometric accuracy for hypothetical framed and frameless workflows. RESULTS: Relative output factors as measured by the manufacturer were independently confirmed using EBT3 film measurements to within 2%. The mean 3D radial discrepancy in stereotactic space between CBCT and CT-based definition over the sampled locations in our rigid geometry phantom was demonstrated to be between 0.40 mm and 0.56 mm over the set of trials, larger than prior reported values. E2E performed in 2D demonstrates sub-mm (and typically < 0.5 mm) accuracy for framed and frameless workflows; geometric accuracy of framed treatments using CBCT-defined stereotactic coordinates was shown to be slightly improved in comparison with those defined using conventional CT. Furthermore, in phantom, frameless workflows exhibited better accuracy than framed workflows for fractionated treatments, despite large magnitudes of introduced interfraction setup error. Accuracy of dosimetric delivery was confirmed in terms of qualitative comparisons of dose profiles and in terms of 2D gamma pass rates based on 1%/1 mm criteria. CONCLUSION: The IU was commissioned for clinical use of frameless and framed treatment protocols. The present study outlines an extensive E2E methodology for confirmation of dosimetric and geometric treatment accuracy.

MV-based relative electron density estimation (iMREDe) for MR-LINAC dose calculation.

BACKGROUND: MR-LINAC systems have been increasingly utilized for real-time imaging in adaptive treatments worldwide. Challenges in MR representation of air cavities and subsequent estimation of electron density maps impede planning efficiency and may lead to dose calculation uncertainties. PURPOSE: To demonstrate the generation of accurate electron density maps using the primary MV beam with a flat-panel imager. METHODS: The ViewRay MRIdian MR-LINAC system was modeled digitally for Monte Carlo simulations. Iron shimming, the magnetic field, and the proposed flat panel detector were included in the model. The effect of the magnetic field on the detector response was investigated. Acquisition of projections over 360 degrees was simulated for digital phantoms of the Catphan 505 phantom and a patient treated for Head and Neck cancer. Shim patterns on the projections were removed and detector noise linearity was assessed. Electron density maps were generated for the digital patient phantom using the flat-panel detector and compared with actual treatment planning CT generated electron density maps of the same patient. RESULTS: The effect of the magnetic field on the detector point-spread function (PSF) was found to be substantial for field strengths above 50 mT. Shims correction in the projection images using air normalization and in-painting effectively removed reconstruction artifacts without affecting noise linearity. The relative difference between reconstructed electron density maps from the proposed method and electron density maps generated from the treatment planning CT was 11% on average along all slices included in the iMREDe reconstruction. CONCLUSIONS: The proposed iMREDe technique demonstrated the feasibility of generating accurate electron densities for the ViewRay MRIdian MR-LINAC system with a flat-panel imager and the primary MV beam. This work is a step towards reducing the time and effort required for adaptive radiotherapy in the current ViewRay MR-LINAC systems.

A kV-MV approach to CBCT metal artifact reduction using multi-layer MV-CBCT.

Objective. To demonstrate that complete cone beam CT (CBCT) scans from both MV-energy and kV-energy LINAC sources can reduce metal artifacts in radiotherapy guidance, while maintaining standard-of-care x-ray doses levels.Approach. MV-CBCT and kV-CBCT scans are acquired at half normal dose. The impact of lowered dose on MV-CBCT data quality is mitigated by the use of a 4-layer MV-imager prototype and reduced LINAC energy settings (2.5 MV) to improve photon capture. Additionally, the MV-CBCT is used to determine the 3D position and pose of metal implants, which in turn is used to guide model-based poly-energetic correction and interleaving of the kV-CBCT and MV-CBCT data. Certain edge-preserving regularization steps incorporated into the model-based correction algorithm further reduce MV data noise.Main results. The method was tested in digital phantoms and a real pelvis phantom with large 2.5″ spherical inserts, emulating hip replacements of different materials. The proposed method demonstrated an appealing compromise between the high contrast of kV-CBCT and low artifact content of MV-CBCT. Contrast-to-noise improved 3-fold compared to MV-CBCT with a clinical 1-layer architecture at matched dose (37 mGy) and edge blur levels. Visual delineation of the bladder and prostate improved noteably over kV- or MV-CBCT alone.Significance. The proposed method demonstrates that a full MV-CBCT scan can be combined with kV-CBCT to reduce metal artifacts without resorting to complicated beam collimation strategies to limit the MV-CBCT dose contribution. Additionally, significant improvements in CNR can be achieved as compared to metal artifact reduction through current clinical MV-CBCT practices.

Lung sparing in MR-guided non-adaptive SBRT treatment of peripheral lung tumors.

Objective.We aim to: (1) quantify the benefits of lung sparing using non-adaptive magnetic resonance guided stereotactic body radiotherapy (MRgSBRT) with advanced motion management for peripheral lung cancers compared to conventional x-ray guided SBRT (ConvSBRT); (2) establish a practical decision-making guidance metric to assist a clinician in selecting the appropriate treatment modality.Approach.Eleven patients with peripheral lung cancer who underwent breath-hold, gated MRgSBRT on an MR-guided linear accelerator (MR linac) were studied. Four-dimensional computed tomography (4DCT)-based retrospective planning using an internal target volume (ITV) was performed to simulate ConvSBRT, which were evaluated against the original MRgSBRT plans. Metrics analyzed included planning target volume (PTV) coverage, various lung metrics and the generalized equivalent unform dose (gEUD). A dosimetric predictor for achievable lung metrics was derived to assist future patient triage across modalities.Main results.PTV coverage was high (median V100% > 98%) and comparable for both modalities. MRgSBRT had significantly lower lung doses as measured by V20 (median 3.2% versus 4.2%), mean lung dose (median 3.3 Gy versus 3.8 Gy) and gEUD. Breath-hold, gated MRgSBRT resulted in an average reduction of 47% in PTV volume and an average increase of 19% in lung volume. Strong correlation existed between lung metrics and the ratio of PTV to lung volumes (RPTV/Lungs) for both modalities, indicating that RPTV/Lungsmay serve as a good predictor for achievable lung metrics without the need for pre-planning. A threshold value of RPTV/Lungs< 0.035 is suggested to achieve V20 < 10% using ConvSBRT. MRgSBRT should otherwise be considered if the threshold cannot be met.Significance.The benefits of lung sparing using MRgSBRT were quantified for peripheral lung tumors; RPTV/Lungswas found to be an effective predictor for achievable lung metrics across modalities. RPTV/Lungscan assist a clinician in selecting the appropriate modality without the need for labor-intensive pre-planning, which has significant practical benefit for a busy clinic.

Monte Carlo model of a prototype flat-panel detector for multi-energy applications in radiotherapy.

BACKGROUND: The incorporation of multi-energy capabilities into radiotherapy flat-panel detectors offers advantages including enhanced soft tissue visualization by reduction of signal from overlapping anatomy such as bone in 2D image projections; creation of virtual monoenergetic images for 3D contrast enhancement, metal artefact reduction and direct acquisition of relative electron density. A novel dual-layer on-board imager offering dual energy processing capabilities is being designed. As opposed to other dual-energy implementation techniques which require separate acquisition with two different x-ray spectra, the dual-layer detector design enables simultaneous acquisition of high and low energy images with a single exposure. A computational framework is required to optimize the design parameters and evaluate detector performance for specific clinical applications. PURPOSE: In this study, we report on the development of a Monte Carlo (MC) model of the imager including model validation. METHODS: The stack-up of the dual-layer imager (DLI) was implemented in GEANT4 Application for Tomographic Emission (GATE). The DLI model has an active area of 43×43 cm2 , with top and bottom Cesium Iodide (CsI) scintillators of 600 and 800 µm thickness, respectively. Measurement of spatial resolution and imaging of dedicated multi-material dual-energy (DE) phantoms were used to validate the model. The modulation transfer function (MTF) of the detector was calculated for a 120 kVp x-ray spectrum using a 0.5 mm thick tantalum edge rotated by 2.5o . For imaging validation, the DE phantom was imaged using a 140 kVp x-ray spectrum. For both validation simulations, corresponding measurements were done using an initial prototype of the imager. Agreement between simulations and measurement was assessed using normalized root mean square error (NRMSE) and 1D profile difference for the MTF and phantom images respectively. Further comparison between measurement and simulation was made using virtual monoenergetic images (VMIs) generated from basis material images derived using precomputed look-up tables. RESULTS: The MTF of the bottom layer of the dual-layer model shows values decreasing more quickly with spatial frequency, compared to the top layer, due to the thicker bottom scintillator thickness and scatter from the top layer. A comparison with measurement shows NRMSE of 0.013 and 0.015 as well as identical MTF50 of 0.8 mm1 and 1.0 mm1 for the top and bottom layer respectively. For the DE imaging of the DE-phantom, although a maximum deviation of 3.3% is observed for the 10 mm aluminum and Teflon inserts at the top layer, the agreement for all other inserts is less than 2.2% of the measured value at both layers. Material decomposition of simulated scatter-free DE images gives an average accuracy in PMMA and aluminum composition of 4.9% and 10.3% for 11-30 mm PMMA and 1-10 mm aluminum objects respectively. A comparison of decomposed values using scatter containing measured and simulated DE images shows good agreement within statistical uncertainty. CONCLUSION: Validation using both MTF and phantom imaging shows good agreement between simulation and measurements. With the present configuration of the digital prototype, the model can generate material decomposed images and virtual monoenergetic images.

Physical characterization of therapeutic proton delivery through common dental materials.

PURPOSE: Dental fixtures are commonplace in an aging, radiation treatment population. The current, local standard of practice in particle therapy is to employ treatment geometries to avoid delivery through implanted dental fixtures. The present study aims to observe the physical effect of delivering therapeutic proton beams through common dental fixture materials as prelude to an eventual goal of assessing the feasibility of using treatment geometries not specified for avoidance of oral implants. A sampling of common dental materials was selected based on prosthodontic consult and was evaluated in terms of relative stopping power and three-dimensional (3D) dose perturbation. METHODS: Amalgams, porcelain-fused-to-metal (PFM) crowns consisting of zirconia and non-noble base metals, and lithium disilicate implants were chosen for analysis. Theoretical stopping power (S) and mass stopping power (S/ρ) were calculated using the Stopping and Range of Ions in Matter (SRIM) application, basing stoichiometric compositions of each fixture on published materials data. S and S/ρ were calculated for a range of historically available compositions of amalgams from 1900 until the current era. The perturbance of S and S/ρ as a function of clinically relevant ranges of amalgam compositions for the modern era was analyzed. Water equivalent thickness (WET) and relative stopping power (Srel ) of each material was measured for a clinical spot-scanning proton beam with monoenergies of 159.9 and 228.8 MeV with a multi-layer ionization chamber (MLIC). Subsequently, 3D dose perturbation was assessed by delivering proton beams through a custom phantom designed to simulate both en-face and on-edge treatment geometries through the selected materials. A treatment plan mimicking the experimental delivery was constructed in the institutional treatment planning system and calculated using TOPAS-based Monte Carlo simulation (MCS). Experimental results were used to validate the MCS. Finally, treatment planning system (TPS) outputs were compared to MCS to determine the accuracy of the dose calculation model. RESULTS: Historical compositions of amalgams ranged in S from 44.8 to 42.9 MeV/cm, with the greatest deviation being observed for the 1900-1959 era. Deviation as a function of amalgam composition from the modern era was most sensitive to proportion of Hg, accounting for deviations up to -4.2% at the greatest clinically relevant concentration. S/ρ was not found to vary greatly between each porcelain and metal alloy material for PFM type crowns. Relative stopping powers ranged between 1.3 and 5.4 for all studied materials, suggesting substantial changes in proton range with respect to water. Film measurements of pristine spots confirm dose perturbance and shortening of proton range, with an upstream shift of each Bragg peak being observed directly behind the installed fixture. At high energies, cold spots were found in all cases directly behind each material feature with a medial fill-in of dose occurring distally. Qualitative agreement of spot perturbance was confirmed between film measurements and MCS. Finally, when comparing integrated depth doses (IDD) by summing over all axial directions, good agreement is observed between TPS and MCS. CONCLUSIONS: All dental materials studied substantially perturbed the dosimetry of pristine proton spots both in terms of WET/Srel as well as the spatial distribution of dose. Proton range was quantifiably shortened, and each dental material affected a cold spot directly behind the object with medial dose back-filling was observed distally. MCS and Eclipse dose calculations exhibited good agreement with measurements, suggesting that treatment planning without employing avoidance strategies may be possible with further investigation.

Synthetic CT generation for MRI-guided adaptive radiotherapy in prostate cancer.

Current MRI-guided adaptive radiotherapy (MRgART) workflows require fraction-specific electron and/or mass density maps, which are created by deformable image registration (DIR) between the simulation CT images and daily MR images. Manual density overrides may also be needed where DIR-produced results are inaccurate. This approach slows the adaptive radiotherapy workflow and introduces additional dosimetric uncertainties, especially in the presence of the magnetic field. This study investigated a method based on a conditional generative adversarial network (cGAN) with a multi-planar method to generate synthetic CT images from low-field MR images to improve efficiency in MRgART workflows for prostate cancer. Fifty-seven male patients, who received MRI-guided radiation therapy to the pelvis using the ViewRay MRIdian Linac, were selected. Forty-five cases were randomly assigned to the training cohort with the remaining twelve cases assigned to the validation/testing cohort. All patient datasets had a semi-paired DIR-deformed CT-sim image and 0.35T MR image acquired using a true fast imaging with steady-state precession (TrueFISP) sequence. Synthetic CT images were compared with deformed CT images to evaluate image quality and dosimetric accuracy. To evaluate the dosimetric accuracy of this method, clinical plans were recalculated on synthetic CT images in the MRIdian treatment planning system. Dose volume histograms for planning target volumes (PTVs) and organs-at-risk (OARs) and dose distributions using gamma analyses were evaluated. The mean-absolute-errors (MAEs) in CT numbers were 30.1 ± 4.2 HU, 19.6 ± 2.3 HU and 158.5 ± 26.0 HU for the whole pelvis, soft tissue, and bone, respectively. The peak signal-to-noise ratio was 35.2 ± 1.7 and the structural index similarity measure was 0.9758 ± 0.0035. The dosimetric difference was on average less than 1% for all PTV and OAR metrics. Plans showed good agreement with gamma pass rates of 99% and 99.9% for 1%/1 mm and 2%/2 mm, respectively. Our study demonstrates the potential of using synthetic CT images created with a multi-planar cGAN method from 0.35T MRI TrueFISP images for the MRgART treatment of prostate radiotherapy. Future work will validate the method in a large cohort of patients and investigate the limitations of the method in the adaptive workflow.

Analysis of the Rate of Re-planning in Spot-Scanning Proton Therapy.

Purpose: Finite proton range affords improved dose conformality of radiation therapy when patient regions-of-interest geometries are well characterized. Substantial changes in patient anatomy necessitate re-planning (RP) to maintain effective, safe treatment. Regularly planned verification scanning (VS) is performed to ensure consistent treatment quality. Substantial resources, however, are required to conduct an effective proton plan verification program, which includes but is not limited to, additional computed tomography (CT) scanner time and dedicated personnel: radiation therapists, medical physicists, physicians, and medical dosimetrists. Materials and Methods: Verification scans (VSs) and re-plans (RPs) of 711 patients treated with proton therapy between June 2015 and June 2018 were studied. All treatment RP was performed with the intent to maintain original plan integrity and coverage. The treatments were classified by anatomic site: brain, craniospinal, bone, spine, head and neck (H&N), lung or chest, breast, prostate, rectum, anus, pelvis, esophagus, liver, abdomen, and extremity. Within each group, the dates of initial simulation scan, number of VSs, number of fractions completed at the time of VS, and the frequency of RP were collected. Data were analyzed in terms of rate of RP and individual likelihood of RP. Results: A total of 2196 VSs and 201 RPs were performed across all treatment sites. H&N and lung or chest disease sites represented the largest proportion of plan modifications in terms of rate of re-plan (RoR: 54% and 58%, respectively) and individual likelihood of RP on a per patient basis (likelihood of RP [RP%]: 46% and 39%, respectively). These sites required RP beyond 4 weeks of treatment, suggesting continued benefit for frequent, periodic VS. Disease sites in the lower pelvis demonstrated a low yield for RP per VS (0.01-0.02), suggesting that decreasing VS frequency, particularly late in treatment, may be reasonable. Conclusions: A large degree of variation in RoR and individual RP% was observed between anatomic treatment sites. The present retrospective analysis provides data to help develop anatomic site-based VS protocols.

The effect of angular dose distribution on the detection of microcalcifications in digital breast tomosynthesis.

PURPOSE: Substantial effort has been devoted to the clinical development of digital breast tomosynthesis (DBT). DBT is a three-dimensional (3D) x-ray imaging modality that reconstructs a number of thin image slices parallel to a stationary detector plane. Preliminary clinical studies have shown that the removal of overlapping breast tissue reduces image clutter and increases detectability of large, low contrast lesions. However, some studies, as well as anecdotal evidence, suggested decreased conspicuity of small, high contrast objects such as microcalcifications. Several investigators have proposed alternative imaging methods for improving microcalcification detection by delivering half of the total dose to the central view in addition to a separate DBT scan. Preliminary observer studies found possible improvement by either viewing the central projection alone or combining all views with a reconstruction algorithm. METHODS: In this paper, we developed a generalized imaging theory based on a cascaded linear-system model for DBT to calculate the effect of variable angular dose distribution on the 3D modulation transfer function (MTF) and noise power spectrum (NPS). Using the ideal observer signal-to-noise ratio (SNR), d', as a figure-of-merit (FOM) for a signal embedded in a uniform background, we compared the detectability of objects with different sizes under different imaging conditions (e.g., angular dose distribution and reconstruction filters). Experimental investigation was conducted for three different angular dose schemes (ADS) using a Siemens Novation(TOMO) prototype unit. RESULTS: Our results show excellent agreement between modeled and experimental measurements of 3D NPS with different angular dose distribution. The ideal observer detectability index for the detection of Gaussian objects with different angular dose distributions depends strongly on the applied reconstruction filter as well as the imaging task. For detection tasks of small calcifications with reconstruction filters used typically in a clinical setting, variable angular dose distribution with more dose delivered to the central views may lead to higher d' than a uniform angular dose distribution. CONCLUSIONS: The conspicuity of the detection of small calcifications may be improved, under certain imaging conditions, by delivering higher dose toward the central views of a tomosynthesis scan, while also reducing the dose at peripheral angles to keep total administered radiation dose equivalent. The degree of improvement depends on the choice of reconstruction filters as well as the imaging task. The improvement is more substantial for high-frequency imaging tasks and when an aggressive slice-thickness (ST) filter is applied to reduced the high-frequency noise at peripheral angles.

Abbreviated on-treatment CBCT using roughness penalized mono-energization of kV-MV data and a multi-layer MV imager.

Simultaneous acquisition of cone beam CT (CBCT) projections using both the kV and MV imagers of an image guided radiotherapy system reduces set-up scan times-a benefit to lung cancer radiation oncology patients-but increases noise in the 3D reconstruction. In this article, we present a kV-MV scan time reduction technique that uses two noise-reducing measures to achieve superior performance. The first is a high-DQE multi-layer MV imager prototype. The second is a beam hardening correction algorithm which combines poly-energetic modeling with edge-preserving, regularized smoothing of the projections. Performance was tested in real acquisitions of the Catphan 604 and a thorax phantom. Percent noise was quantified from voxel values in a soft tissue volume of interest (VOI) while edge blur was quantified from a VOI straddling a boundary between air and soft material. Comparisons in noise/resolution performance trade-off were made between our proposed approach, a dose-equivalent kV-only scan, and a kV-MV reconstruction technique previously published by Yinet al(2005Med. Phys.329). The proposed technique demonstrated lower noise as a function of spatial resolution than the baseline kV-MV method, notably a 50% noise reduction at typical edge blur levels. Our proposed method also exhibited fainter non-uniformity artifacts and in some cases superior contrast. Overall, we find that the combination of a multi-layer MV imager, acquiring at a LINAC source energy of 2.5 MV, and a denoised beam hardening correction algorithm enables noise, resolution, and dose performance comparable to standard kV-imager only set-up CBCT, but with nearly half the gantry rotation time.

Low-dose megavoltage cone-beam computed tomography using a novel multi-layer imager (MLI).

PURPOSE: The feasibility of low-dose megavoltage cone-beam acquisition (MVCBCT) using a novel, high detective quantum efficiency (DQE) multi-layer imager (MLI) was investigated. The aim of this work was to reconstruct MVCBCT images using the MLI at different total dose levels, and assess Hounsfield Unit (HU) accuracy, noise and contrast-to-noise ratio (CNR) for low-dose megavoltage cone-beam acquisition. METHODS: The MLI has four stacked layers; each layer contains a combination of copper filter/converter, gadolinium oxysulfide (GOS) scintillator and a-Si detector array. In total, 720 projections of a CATPHAN® phantom were acquired over 360° at 2.5, 6, and 6 MV flattening filter free (FFF) beam energies on a Varian TrueBeam LINAC. The dose per projection was 0.01, 0.0167, and 0.05 MU for 2.5, 6, and 6 MV FFF, respectively. MVCBCT images were reconstructed with varying numbers of projections to provide a range of doses for evaluation. Hounsfield Unit uniformity, accuracy, noise and CNR were estimated. Improvements were quantified relative to the standard AS1200 single-layer imager. RESULTS: Average HU uniformity for the MLI reconstructions was within a range of 95%-99% for all of the energies studied. Relative electron density estimation from HU values was within 0.4% ± 1.8% from nominal values. The CNR for MVCBCT based on MLI projections was 2-4× greater than from AS1200 projections. The 2.5 MV beam acquisition with the MLI exhibited the lowest noise and the best balance between CNR and dose for low-dose reconstructions. CONCLUSIONS: Megavoltage cone-beam acquisition imaging with a novel MLI prototype mounted on a clinical linear accelerator demonstrated substantial improvement over the standard AS1200 EPID. Further optimization of MVCBCT reconstruction, particularly for 2.5 MV acquisitions, will improve image metrics. Overall, the MLI improves CNR at substantially lower doses than currently required by conventional detectors. This new high DQE detector could provide high-quality MVCBCT at clinically acceptable doses.

Low-Dose Hypofractionated Total Skin Electron Beam Therapy for Adult Cutaneous T-Cell Lymphoma.

PURPOSE: Historically, the standard of care for total skin electron beam therapy (TSEBT) delivered 30 to 36 Gy over 5 to 10 weeks. Given the high risk of relapse, a majority of patients require additional treatments. Therefore, attempts to use a shortened course of TSEBT have been investigated. METHODS AND MATERIALS: We conducted a single-institution retrospective review to evaluate disease response, control, and toxicity using a low-dose, hypofractionated course of TSEBT (HTSEBT) in patients with mycosis fungoides. RESULTS: Forty patients received 57 courses of HTSEBT. Median dose (Gy)/fractionation was 12/3, spanning a median time of 2.4 weeks. Overall response rate of patients assessed (n = 54) was 100%. Thirty-one courses (57.4%) resulted in a complete response and 23 courses (42.6%) resulted in a partial response. Cumulative incidence of progressive skin disease at 3 months was 37.2%, at 6 months, 56.9%, and at 1 year, 81.5%. Of the 40 patients treated with a first course of HTSEBT, 31 received subsequent courses of radiotherapy. Cumulative incidence of subsequent treatment was 28.0% at 3 months, 46.8% at 6 months, and 70.0% at 1 year. Patients who underwent repeat courses of HTSEBT continued to have similar treatment responses to repeat courses without increased toxicities. Toxicities from all courses were acceptable with the exception of 1 patient, who experienced grade 4 skin toxicity (moist desquamation requiring hospitalization). CONCLUSIONS: Low-dose HTSEBT provides good palliation in patients with cutaneous T-cell lymphoma with a satisfactory response and toxicity profile. HTSEBT allows therapy to be completed in far fewer treatments. Low-dose HTSEBT is an appropriate treatment option for patients unable to come for daily treatment. HTSEBT provides a way to decrease exposure to other patients and staff during public health emergencies such as the coronavirus disease 2019 (COVID-19) pandemic.

Experimental validation of a three-dimensional linear system model for breast tomosynthesis.

A three-dimensional (3D) linear model for digital breast tomosynthesis (DBT) was developed to investigate the effects of different imaging system parameters on the reconstructed image quality. In the present work, experimental validation of the model was performed on a prototype DBT system equipped with an amorphous selenium (a-Se) digital mammography detector and filtered back-projection (FBP) reconstruction methods. The detector can be operated in either full resolution with 85 microm pixel size or 2 x 1 pixel binning mode to reduce acquisition time. Twenty-five projection images were acquired with a nominal angular range of +/- 20 degrees. The images were reconstructed using a slice thickness of 1 mm with 0.085 x 0.085 mm in-plane pixel dimension. The imaging performance was characterized by spatial frequency-dependent parameters including a 3D noise power spectrum (NPS) and in-plane modulation transfer function (MTF). Scatter-free uniform x-ray images were acquired at four different exposure levels for noise analysis. An aluminum (Al) edge phantom with 0.2 mm thickness was imaged to measure the in-plane presampling MTF. The measured in-plane MTF and 3D NPS were both in good agreement with the model. The dependence of DBT image quality on reconstruction filters was investigated. It was found that the slice thickness (ST) filter, a Hanning window to limit the high-frequency components in the slice thickness direction, reduces noise aliasing and improves 3D DQE. An ACR phantom was imaged to investigate the effects of angular range and detector operational modes on reconstructed image quality. It was found that increasing the angular range improves the MTF at low frequencies, resulting in better detection of large-area, low-contrast mass lesions in the phantom. There is a trade-off between noise and resolution for pixel binning and full resolution modes, and the choice of detector mode will depend on radiation dose and the targeted lesion.

A novel method for fast image simulation of flat panel detectors.

We have developed a novel method for fast image simulation of flat panel detectors, based on the photon energy deposition efficiency and the optical spread function (OSF). The proposed method, FastEPID, determines the photon detection using photon energy deposition and replaces particle transport within the detector with precalculated OSFs. The FastEPID results are validated against experimental measurement and conventional Monte Carlo simulation in terms of modulation transfer function (MTF), signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), contrast, and relative difference of pixel value, obtained with a slanted slit image, Las Vegas phantom, and anthropomorphic pelvis phantom. Excellent agreement is observed between simulation and measurement in all cases. Without degrading image quality, the FastEPID method is capable of reducing simulation time up to a factor of 150. Multiple applications, such as imager design optimization for planar and volumetric imaging, are expected to benefit from the implementation of the FastEPID method.

Characterizing a novel scintillating glass for application to megavoltage cone-beam computed tomography.

PURPOSE: The purpose of this study was to evaluate the performance of a prototype electric portal imaging device (EPID) with a high detective quantum efficiency (DQE) scintillator, LKH-5. Specifically, image quality in context of both planar and megavoltage (MV) cone-beam computed tomography (CBCT) is analyzed. METHODS: Planar image quality in terms of modulation transfer function (MTF), noise power spectrum (NPS), and DQE are measured and compared to an existing EPID (AS-1200) using the 6 MV beamline for a Varian TrueBeam linac. Imager performance is contextualized for three-dimensional (3D), MV-CBCT performance by measuring imager lag and analyzing the expected degradation of the DQE as a function of dose. Finally, comparisons between reconstructed images of the Catphan phantom in terms of qualitative quality and signal-difference-to-noise ratio (SDNR) are made for 6 MV images using both conventional and LKH-5 EPIDs as well as for the kilovoltage (kV) on-board imager (OBI). RESULTS: Analysis of the NPS reveals linearity at all measured doses using the prototype LKH-5 detector. While the first zero of the MTF is much lower for the LKH-5 detector than the conventional EPID (0.6 cycles/mm vs 1.6 cycles/mm), the normalized NPS (NNPS) multiplied by total quanta (qNNPS) of the LKH-5 detector is roughly a factor of seven to eight times lower, yielding a DQE(0) of approximately 8%. First, second, and third frame lag were measured at approximately 23%, 5%, and 1%, respectively, although no noticeable image artifacts were apparent in reconstructed volumes. Analysis of low-dose performance reveals that DQE(0) remains at 80% of its maximum value at a dose as low as 7.5 × 10-6  MU. For a 400 projection technique, this represents a total scan dose of 0.0030 MU, suggesting that if imaging doses are increased to a value typical of kV-CBCT scans (~2.7 cGy), the LKH-5 detector will retain quantum noise limited performance. Finally, comparing Catphan scans, the prototype detector exhibits much lower image noise than the conventional EPID, resulting in improved small object representation. Furthermore, SDNR of H2 O and polystyrene cylinders improved from -1.95 and 2.94 to -15 and 18.7, respectively. CONCLUSIONS: Imaging performance of the prototype LKH-5 detector was measured and analyzed for both planar and 3D contexts. Improving noise transfer of the detector results in concurrent improvement of DQE(0). For 3D imaging, temporal characteristics were adequate for artifact-free performance and at relevant doses, the detector retained quantum noise limited performance. Although quantitative MTF measurements suggest poorer resolution, small object representation of the prototype imager is qualitatively improved over the conventional detector due to the measured reduction in noise.

Image artifacts in digital breast tomosynthesis: investigation of the effects of system geometry and reconstruction parameters using a linear system approach.

Digital breast tomosynthesis (DBT) is a three-dimensional (3D) x-ray imaging modality that reconstructs image slices parallel to the detector plane. Image acquisition is performed using a limited angular range (less than 50 degrees) and a limited number of projection views (less than 50 views). Due to incomplete data sampling, image artifacts are unavoidable in DBT. In this preliminary study, the image artifacts in DBT were investigated systematically using a linear system approximation. A cascaded linear system model of DBT was developed to calculate the 3D presampling modulation transfer function (MTF) with different image acquisition geometries and reconstruction filters using a filtered backprojection (FBP) algorithm. A thin, slanted tungsten (W) wire was used to measure the presampling MTF of the DBT system in the cross-sectional plane defined by the thickness (z-) and tube travel (x-) directions. The measurement was in excellent agreement with the calculation using the model. A small steel bead was used to calculate the artifact spread function (ASF) of the DBT system. The ASF was correlated with the convolution of the two-dimensional (2D) point spread function (PSF) of the system and the object function of the bead. The results showed that the cascaded linear system model can be used to predict the magnitude of image artifacts of small, high-contrast objects with different image acquisition geometry and reconstruction filters.

A Monte Carlo study of the impact of phosphor optical properties on EPID imaging performance.

We have developed a Monte Carlo computational model of a clinically employed electronic portal imaging device (EPID), and demonstrated the impact of phosphor optical properties on the imaging performance. The EPID model was built with Geant4 application for tomographic emission. Both radiative and optical transport were included in the model. Modulation transfer function (MTF), normalized noise-power spectrum times the incident x-ray fluence (qNNPS), and detective quantum efficiency (DQE) were calculated for simulated and measured data, and their agreement was quantified by the normalized root-mean-square error (NRMSE). MTF was computed using a 100 µm wide slit tilted by 1.5° and qNNPS was estimated using the Fujita-Lubberts-Swank method. DQE was calculated from MTF and qNNPS data. The NRMSE value was 0.0467 for MTF, 0.0217 for qNNPS, and 0.0885 for DQE, showing good agreement between measurement and simulation. Five major optical properties, phosphor grain size, phosphor thickness, phosphor refractive index, binder refractive index, and packing ratio were tested for their influence on the qNNPS, MTF, and DQE(0) of the model. Generally, the effect on the qNNPS is greater than MTF, and no impact on DQE(0), except from phosphor thickness, was observed. Multiple applications, such as imager design optimization and investigations of the dosimetric performance, are expected to benefit from the validated model.

Feasibility of closed-MLC tracking using high sensitivity and multi-layer electronic portal imagers.

In radiation therapy, improvements in treatment conformality are often limited by movement of target tissue. To better treat the target, tumor tracking strategies involving beam's-eye-view (BEV) have been explored. However, localization surrogates like implanted fiducial markers may sometimes leave the field-of-view (FOV), as defined by the linear accelerator (LINAC) multi-leaf collimator (MLC). Radiation leakage through the MLC has been measured previously at approximately 1%-2%. High sensitivity prototype detectors imagers may improve the ability to visualize objects outside of the MLC FOV during treatment. The present study presents a proof-of-concept for tracking fiducial markers outside the MLC FOV by employing high sensitivity detectors using a high-efficiency, prototype scintillating glass called LKH-5 and also investigates the impact of multi-layer imager (MLI) architecture. It was found that by improving the detector efficiency, using either of these methods results in a reduction of dose required for fiducial marker visibility. Further, image correction by a rectangular median filter will improve fiducial marker representation in the MLC blocked images. Quantified by measuring the peak-to-sidelobe ratio (PSR) of the normalized cross correlation (NCC) between a template of the fiducial marker with the blocked MLC acquisition, visibility has been found at a threshold of roughly 5 for all configurations with a 3 × 3 cm2 ROI. For typical gadolinium oxysulfide (GOS) detectors in single and simulated 4-layer configurations, the minimum dose required for visualization was 20 and 10 MU, respectively. For LKH-5 detectors in single and simulated 4-layer configurations, this minimum dose was reduced to 4 and 2 MU, respectively. With a 6 MV flattening filter free (FFF) beam dose rate of 1400 MU min-1, the maximum detector frame rate while maintaining fiducial visibility is approximately 12 fps for a 4-layer LKH-5 configuration.