Effect of scatter suppression with 2D antiscatter grids in photon counting compact CBCT.

Energy sensitive and photon counting detectors can provide improved tissue visualization and material quantification capabilities in Cone Beam Computed Tomography (CBCT) systems. However, their implementation in CBCT systems is more challenging, which is in part due to high fluence of scattered X-rays in wide cone angle CBCT geometry. Specifically, high scatter contamination in lower energy spectrum challenges reconstruction of high fidelity CBCT images by using lower energy X-rays. To address this problem, we investigated a robust scatter rejection with 2D antiscatter grids in a benchtop photon counting and compact CBCT system. The benchtop system employs a 35 cm wide CdTe photon counting detector with two energy thresholds. To reject scatter, a dedicated 2D antiscatter grid (2D grid) prototype made from tungsten was developed and mounted directly on the detector. To correct residual scatter not stopped by the 2D grid, a measurement-based scatter correction method, referred to as Grid-based Scatter Sampling (GSS), was utilized. Without 2D grid, scatter to primary ratio (SPR) reached 2.3 in the 15-40 keV energy bin. SPR was factor of 3 higher in the lowest energy bin when compared to the highest energy bin (90-120 keV). With the 2D grid, SPR was reduced below 0.14, and SPR values were more homogenous across the energy spectrum. CT number nonuniformity was factor of 3 lower in both low and high energy bin CBCT reconstructions. Improvement in contrast to noise ratio and contrast was more pronounced in the low energy bin CBCT images. This work indicates that 2D grids can significantly reduce spectral contamination caused by scatter in photon counting compact CBCT, and potentially enable higher fidelity CBCT image reconstructions.

2D antiscatter grid and scatter sampling based CBCT method for online dose calculations during CBCT guided radiation therapy of pelvis.

BACKGROUND: Online dose calculations before the delivery of radiation treatments have applications in dose delivery verification, online adaptation of treatment plans, and simulation-free treatment planning. While dose calculations by directly utilizing CBCT images are desired, dosimetric accuracy can be compromised due to relatively lower HU accuracy in CBCT images. PURPOSE: In this work, we propose a novel CBCT imaging pipeline to enhance the accuracy of CBCT-based dose calculations in the pelvis region. Our approach aims to improve the HU accuracy in CBCT images, thereby improving the overall accuracy of CBCT-based dose calculations prior to radiation treatment delivery. METHODS: An in-house developed quantitative CBCT pipeline was implemented to address the CBCT raw data contamination problem. The pipeline combines algorithmic data correction strategies and 2D antiscatter grid-based scatter rejection to achieve high CT number accuracy. To evaluate the effect of the quantitative CBCT pipeline on CBCT-based dose calculations, phantoms mimicking pelvis anatomy were scanned using a linac-mounted CBCT system, and a gold standard multidetector CT used for treatment planning (pCT). A total of 20 intensity-modulated treatment plans were generated for five targets, using 6 and 10 MV flattening filter-free beams, and utilizing small and large pelvis phantom images. For each treatment plan, four different dose calculations were performed using pCT images and three CBCT imaging configurations: quantitative CBCT, clinical CBCT protocol, and a high-performance 1D antiscatter grid (1D ASG). Subsequently, dosimetric accuracy was evaluated for both targets and organs at risk as a function of patient size, target location, beam energy, and CBCT imaging configuration. RESULTS: When compared to the gold-standard pCT, dosimetric errors in quantitative CBCT-based dose calculations were not significant across all phantom sizes, beam energies, and treatment sites. The largest error observed was 0.6% among all dose volume histogram metrics and evaluated dose calculations. In contrast, dosimetric errors reached up to 7% and 97% in clinical CBCT and high-performance ASG CBCT-based treatment plans, respectively. The largest dosimetric errors were observed in bony targets in the large phantom treated with 6 MV beams. The trends of dosimetric errors in organs at risk were similar to those observed in the targets. CONCLUSIONS: The proposed quantitative CBCT pipeline has the potential to provide comparable dose calculation accuracy to the gold-standard planning CT in photon radiation therapy for the abdomen and pelvis. These robust dose calculations could eliminate the need for density overrides in CBCT images and enable direct utilization of CBCT images for dose delivery monitoring or online treatment plan adaptations before the delivery of radiation treatments.

Feasibility of dual-energy CBCT material decomposition in the human torso with 2D anti-scatter grids and grid-based scatter sampling.

BACKGROUND: Dual-energy (DE) imaging techniques in cone-beam computed tomography (CBCT) have potential clinical applications, including material quantification and improved tissue visualization. However, the performance of DE CBCT is limited by the effects of scattered radiation, which restricts its use to small object imaging. PURPOSE: This study investigates the feasibility of DE CBCT material decomposition by reducing scatter with a 2D anti-scatter grid and a measurement-based scatter correction method. Specifically, the investigation focuses on iodine quantification accuracy and virtual monoenergetic (VME) imaging in phantoms that mimic head, thorax, abdomen, and pelvis anatomies. METHODS: A 2D anti-scatter grid prototype was utilized with a residual scatter correction method in a linac-mounted CBCT system to investigate the effects of robust scatter suppression in DE CBCT. Scans were acquired at 90 and 140 kVp using phantoms that mimic head, thorax, and abdomen/pelvis anatomies. Iodine vials with varying concentrations were placed in each phantom, and CBCT images were decomposed into iodine and water basis material images. The effect of a 2D anti-scatter grid with and without residual scatter correction on iodine concentration quantification and contrast visualization in VME images was evaluated. To benchmark iodine concentration quantification accuracy, a similar set of experiments and DE processing were also performed with a conventional multidetector CT scanner. RESULTS: In CBCT images, a 2D grid with or without scatter correction can differentiate iodine and water after DE processing in human torso-sized phantom images. However, iodine quantification errors were up to 10 mg/mL in pelvis phantoms when only the 2D grid was used. Adding scatter correction to 2D-grid CBCT reduced iodine quantification errors below 1.5 mg/mL in pelvis phantoms, comparable to iodine quantification errors in multidetector CT. While a noticeable contrast-to-noise ratio improvement was not observed in VME CBCT images, contrast visualization was substantially better in 40 keV VME images in visual comparisons with 90 and 140 kVp CBCT images across all phantom sizes investigated. CONCLUSIONS: This study indicates that accurate DE decomposition is potentially feasible in DE CBCT of the human torso if robust scatter suppression is achieved with 2D anti-scatter grids and residual scatter correction. This approach can potentially enable better contrast visualization and tissue and contrast agent quantification in various CBCT applications.

A 2D antiscatter grid and scatter sampling based CBCT method for online dose calculations during CBCT guided radiation therapy of pelvis.

Online dose calculations before radiation treatment have applications in dose delivery verification, plan adaptation, and treatment planning. We propose a novel CBCT imaging pipeline to enhance accuracy. Our approach aims to improve HU accuracy in CBCT images for more precise dose calculations. A quantitative CBCT pipeline was implemented, combining data correction strategies and scatter rejection, achieving high CT number accuracy. We evaluated the pipeline's effect using pelvis anatomy phantoms and found that dosimetric errors in quantitative CBCT-based dose calculations were minimal. In contrast, clinical CBCT and high-performance ASG CBCT-based plans showed significant errors. The proposed quantitative CBCT pipeline offers comparable dose calculation accuracy to the gold-standard planning CT, eliminating the need for density overrides and enabling precise dose delivery monitoring or online plan adaptations in radiation therapy.

A simulation study to evaluate the effect of 2D antiscatter grid primary transmission on flat panel detector based CBCT image quality.

Purpose. Two-dimensional antiscatter grids' (2D-ASGs) septal shadows and their impact on primary transmission play a critical role in cone-beam computed tomography (CBCT) image noise and artifact characteristics. Therefore, a numerical simulation platform was developed to evaluate the effect of 2D-ASG's primary transmission on image quality, as a function of grid geometry and CBCT system properties.Methods. To study the effect of 2D-ASG's septal shadows on primary transmission and CBCT image quality, two new methods were introduced; one to simulate projection signal gradients in septal shadows, and the other to simulate septal shadow variations due to gantry flex. Signal gradients in septal shadows were simulated by generating a system point spread function that was directly extracted from projection images of 2D-ASG prototypes in experiments. Variations in septal shadows due to gantry flex were simulated by generating oversampled shadow profiles extracted from experiments. Subsequently, the effect of 2D-ASG's septal shadows on primary transmission and image quality was evaluated.Results.For an apparent septal thickness of 0.15 mm, the primary transmission of 2D-ASG varied between 72%-90% for grid pitches 1-3 mm. In low-contrast phantoms, the effect of 2D-ASG's radiopaque footprint on information loss was subtle. At high spatial frequencies, information loss manifested itself as undersampling artifacts, however, its impact on image quality is subtle when compared to quantum noise. Effects of additive electronic noise and gantry flex induced ring artifacts on image quality varied as a function of grid pitch and septal thickness. Such artifacts were substantially less in lower resolution images.Conclusion. The proposed simulation platform allowed successful evaluation of CBCT image quality variations as a function of 2D-ASG primary transmission properties and CBCT system characteristics. This platform can be potentially used for optimizing 2D-ASG design properties based on the imaging task and properties of the CBCT system.

A quantitative CBCT pipeline based on 2D antiscatter grid and grid-based scatter sampling for image-guided radiation therapy.

BACKGROUND: Quantitative accuracy is critical for expanding the role of cone beam CT (CBCT) imaging from target localization to quantitative treatment monitoring and plan adaptations in radiation therapy. Despite advances in CBCT image quality improvement methods, quantitative accuracy gap between CBCT and multi-detector CT (MDCT) remains. PURPOSE: In this work, a physics-driven approach was investigated that combined robust scatter rejection, raw data correction and iterative image reconstruction to further improve CBCT image quality and quantitative accuracy, referred to as quantitative CBCT (qCBCT). METHODS: QCBCT approach includes tungsten 2D antiscatter grid hardware, residual scatter correction with grid-based scatter sampling, image lag, and beam hardening correction for offset detector geometry linac-mounted CBCT. Images were reconstructed with iterative image reconstruction to reduce image noise. qCBCT was evaluated using a variety of phantoms to investigate the effect of object size and its composition on image quality, and image quality was benchmarked against clinical CBCT and gold standard MDCT images used for treatment planning. RESULTS: QCBCT provided statistically significant improvement in CT number accuracy and reduced image artifacts when compared to clinical CBCT images. When compared to gold standard MDCT, mean HU errors in qCBCT and clinical CBCT were 17 ± 9 and 38 ± 29 HU, respectively. Magnitude of phantom size dependent HU variations were comparable between MDCT and qCBCT images. With iterative reconstruction, contrast-to-noise ratio improved by 25% when compared to clinical CBCT protocols. CONCLUSIONS: Combination of novel scatter suppression techniques and other data correction methods in qCBCT provided CT number accuracy comparable to gold standard MDCT used for treatment planning. This approach may potentially improve CBCT's promise in fulfilling the tasks that demand high quantitative accuracy, such as online dose calculations and treatment response assessment, in image guided radiation therapy.

2D antiscatter grid and scatter sampling based CBCT pipeline for image guided radiation therapy.

Poor tissue visualization and quantitative accuracy in CBCT is a major barrier in expanding the role of CBCT imaging from target localization to quantitative treatment monitoring and plan adaptations in radiation therapy sessions. To further improve image quality in CBCT, 2D antiscatter grid based scatter rejection was combined with a raw data processing pipeline and iterative image reconstruction. The culmination of these steps was referred as quantitative CBCT, qCBCT. qCBCT data processing steps include 2D antiscatter grid implementation, measurement based residual scatter, image lag, and beam hardening correction for offset detector geometry CBCT with a bow tie filter. Images were reconstructed with iterative image reconstruction to reduce image noise. To evaluate image quality, qCBCT acquisitions were performed using a variety of phantoms to investigate the effect of object size and its composition on image quality. qCBCT image quality was benchmarked against clinical CBCT and MDCT images. Addition of image lag and beam hardening correction to scatter suppression reduced HU degradation in qCBCT by 10 HU and 40 HU, respectively. When compared to gold standard MDCT, mean HU errors in qCBCT and clinical CBCT were 10 HU and 27 HU, respectively. HU inaccuracy due to change in phantom size was 22 HU and 85 HU in qCBCT and clinical CBCT images, respectively. With iterative reconstruction, contrast to noise ratio improved by a factor of 1.25 when compared to clinical CBCT protocols. Robust artifact and noise suppression in qCBCT images can reduce the image quality gap between CBCT and MDCT, improving the promise of qCBCT in fulfilling the tasks that demand high quantitative accuracy, such as CBCT based dose calculations and treatment response assessment in image guided radiation therapy.

Concurrent kilovoltage CBCT imaging and megavoltage beam delivery: suppression of cross-scatter with 2D antiscatter grids and grid-based scatter sampling.

Objective. The concept of using kilovoltage (kV) and megavoltage (MV) beams concurrently has potential applications in cone beam computed tomography (CBCT) guided radiation therapy, such as single breath hold scans, metal artifact reduction, and simultaneous imaging during MV treatment delivery. However, MV cross-scatter generated during MV beam delivery degrades CBCT image quality. To address this, a 2D antiscatter grid and a cross-scatter correction method were investigated in the context of high dose MV treatment delivery.Approach. A 3D printed, tungsten 2D antiscatter grid prototype was utilized in kV CBCT scans to reduce MV cross-scatter fluence during concurrent MV beam delivery. Remaining cross-scatter in projections was corrected by using the 2D grid itself as a cross-scatter intensity sampling device, referred to as grid-based scatter sampling (GSS). To test this approach, kV CBCT acquisitions were performed while delivering 6 and 10 MV beams, mimicking high dose rate treatment delivery scenarios. kV and MV beam deliveries were not synchronized to eliminate MV beam delivery interruption. MV cross-scatter suppression performance of the proposed approach was evaluated in projections and CBCT images of phantoms.Main results. 2D grid reduced the intensity of MV cross-scatter in kV projections by a factor of 3 on the average, when compared to conventional antiscatter grid. Remaining cross scatter as measured by the GSS method was within 7% of measured reference intensity values, and subsequently corrected. CBCT image quality was improved substantially during concurrent kV-MV beam delivery. Median Hounsfield Unit (HU) inaccuracy was up to 182 HU without our methods, and it was reduced to a median 6.5 HU with our 2D grid and scatter correction approach. Our methods provided a factor of 2-6 improvement in contrast-to-noise ratio.Significance. This investigation demonstrates the utility of 2D antiscatter grids and grid-based scatter sampling in suppressing MV cross-scatter. Our approach successfully minimized the effects of MV cross-scatter in concurrent kV CBCT imaging and high dose MV treatment delivery scenarios. Hence, robust MV cross-scatter suppression is potentially feasible without MV beam delivery interruption or compromising kV image acquisition rate.

Megavoltage cross-scatter rejection and correction using 2D antiscatter grids in kilovoltage CBCT imaging.

Simultaneous use of kilovoltage (kV) and megavoltage (MV) beams has numerous potential applications in cone beam computed tomography (CBCT)-guided radiotherapy, such as fast MV+kV CBCT for single breath-hold scan, tumor localization with kV CBCT imaging during MV therapy delivery, and metal artifact suppression. However, the introduction of MV beams results in a large MV-cross scatter fluence incident on the kV Flat Panel Detector (FPD), and thus, deteriorating the low contrast visualization and Hounsfield Unit (HU) accuracy. In this work, we introduced a novel and robust method for reducing the effects of MV cross scatter. First, we implemented a 2D antiscatter grid atop the detector which rejects a large section of MV cross scatter. This hardware-based approach, while effective, allows a fraction of MV cross scatter to be transmitted to the FPD, resulting in artifacts and degraded HU accuracy in CBCT images. We thus introduced a data correction step, which aimed to estimate and correct the remaining MV cross scatter. This approach, referred to as Grid-Based Scatter Sampling, utilized 2D antiscatter grid itself to measure and correct remaining MV cross scatter in projections. We investigated the performance of the proposed approach in experiments by simultaneously acquiring kV CBCT and delivering MV beams with a clinical linac. The results show that the proposed method can substantially reduce HU inaccuracy and increase contrast-to-noise ratio (CNR). Our method does not require synchronization of kV and MV beam pulses, reduction of kV frame acquisition rate, or MV dose rate, and therefore, it is more practical to implement in radiation therapy clinical setting.

Evaluation of scatter rejection and correction performance of 2D antiscatter grids in cone beam computed tomography.

PURPOSE: We have been investigating two-dimensional (2D) antiscatter grids (2D ASGs) to reduce scatter fluence and improve image quality in cone beam computed tomography (CBCT). In this work, two different aspects of 2D ASGs, their scatter rejection and correction capability, were investigated in CBCT experiments. To correct residual scatter transmitted through the 2D ASG, it was used as a scatter measurement device with a novel method: grid-based scatter sampling. METHODS: Three focused 2D ASG prototypes with grid ratios of 8, 12, and 16 were developed for linac-mounted offset detector CBCT geometry. In the first phase, 2D ASGs were used as a scatter rejection device, and the effect of grid ratio on CT number accuracy and contrast-to-noise ratio (CNR) evaluated in CBCT images. In the second phase, a grid-based scatter sampling method which exploits the signal modulation characteristics of the 2D ASG's septal shadows to measure and correct residual scatter transmitted through the grid was implemented. To evaluate CT number accuracy, the percent change in CT numbers was measured by changing the phantom from head to pelvis size and configuration. RESULTS: When 2D ASG was used as a scatter rejection device, CT number accuracy increased and the CT number variation due to change in phantom dimensions was reduced from 23% to 2-6%. A grid ratio of 16 yielded the lowest CT number variation. All three 2D ASGs yielded improvement in CNR, up to a factor of two in pelvis-sized phantoms. When 2D ASG prototypes were used for both scatter rejection and correction, CT number variations were reduced further, to 1.3-2.6%. In comparisons with a clinical CBCT system and a high-performance radiographic ASG, 2D ASG provided higher CT number accuracy under the same imaging conditions. CONCLUSIONS: When 2D ASG is used solely as a scatter rejection device, substantial improvement in CT number accuracy can be achieved by increasing the grid ratio. Two-dimensional ASGs also provided significant CNR improvement even at lower grid ratios. Two-dimensional ASGs used in conjunction with the grid-based scatter sampling method provided further improvement in CT number accuracy, irrespective of the grid ratio, while preserving 2D ASGs' capacity to improve CNR. The combined effect of scatter rejection and residual scatter correction by 2D ASG may accelerate implementation of new techniques in CBCT that require high quantitative accuracy, such as radiotherapy dose calculation and dual energy CBCT.

A unified scatter rejection and correction method for cone beam computed tomography.

PURPOSE: Scattered radiation is a major cause of image quality degradation in flat panel detector-based cone beam CT (CBCT). While recently introduced 2D antiscatter grids reject the majority of scatter fluence, the small percentage of scatter fluence still transmitted to the detector remains a major challenge for implementation of quantitative imaging techniques such as dual energy imaging in CBCT. Additionally, this residual scatter is also a major source of grid-induced artifacts, which impedes implementation of 2D grids in CBCT. We therefore present a new method to achieve both robust scatter rejection and residual scatter correction using a 2D antiscatter grid; in doing so, we expand the role of 2D grids from mere scatter rejection devices to scatter measurement devices. METHOD: In our method, the radiopaque septa of the 2D grid emulate a micro array of beam-stops placed on the detector which introduce spatially periodic septal shadows. By selecting sufficiently thin grid septa, the primary intensity can be reduced while preserving the uniformity of scatter intensity. This enables us to correlate the modulated pixel signal intensity in septal shadows with local scatter intensity. Our method then exploits this correlation to measure and remove residual scatter intensity from projections. No assumptions are made about the object being imaged. We refer to this as Grid-based Scatter Sampling (GSS). In this work, we evaluate the principle of signal modulation with grid septa, the accuracy of scatter estimates, and the effect of the GSS method on image quality using simulations and measurements. We also implement the GSS method experimentally using a 2D grid prototype. RESULTS: Our results demonstrate that the GSS method increased CT number accuracy and reduced image artifacts associated with scatter. With 2D grid and residual scatter correction, HU nonuniformity was reduced from 65 HU to 30 HU in pelvis sized phantoms, and HU variations due to change in phantom size were reduced from 59 HU to 20 HU, when compared to use of only a 2D grid. With residual scatter correction via GSS method, grid-induced ring artifacts were suppressed, leading to a 41% reduction in noise. The shape of the modulation transfer function (MTF) was preserved before and after suppression of ring artifacts. CONCLUSIONS: Our grid-based scatter sampling method enables utilization of a 2D grid as a scatter measurement and correction device. This method significantly improves quantitative accuracy in CBCT, further reducing the image quality gap between CBCT and multi-detector CT. By correcting residual scatter with the proposed method, grid-induced line artifacts in projections and associated ring artifacts in CBCT images were also suppressed with no compromise of spatial resolution.

Simultaneous scatter rejection and correction method using 2D antiscatter grids for CBCT.

While two-dimensional antiscatter grids (2D grid) reduce scatter intensity substantially in Cone Beam Computed Tomography (CBCT), a small fraction of scattered radiation is transmitted through the 2D grid to the detector. Residual scatter limits the accuracy of CT numbers and interferes with the correction of grid's septal shadows, or footprint, in projections. If grid's septal shadows are not adequately suppressed in projections, it will lead to ring artifacts in CBCT images. In this work, we present a new method to correct residual scatter transmitted through the grid by employing the 2D grid itself as a residual scatter measurement device. Our method, referred as grid-based scatter sampling (GSS), exploits the spatial modulation of primary x-ray fluence by 2D grid's septal shadows. The shape of the signal modulation pattern varies as a function of residual scatter intensity registered by detector pixels. Such a variation in signal pattern was employed to measure residual scatter intensity in each projection, and subsequently, residual scatter was subtracted from each projection. To validate the GSS method, CBCT imaging experiments were conducted using a 2D antiscatter grid prototype in a linac mounted CBCT system. The effect of GSS method on the ring artifact reduction was quantified by measuring noise in CBCT images. In addition, the nonuniformity of Hounsfield Units (HU) and HU accuracy was measured in both head and pelvis-sized phantoms. In qualitative evaluations, GSS method successfully reduced ring artifacts caused by 2D grid's footprint. Image noise reduced by 23% due to reduction of ring artifacts. HU nonuniformity in water-equivalent sections was reduced from 20 HU to 10 HU, and streak artifacts between high density inserts were reduced. The phantom size dependent variations in HU was also reduced after application of GSS method. Without GSS method, HU of density inserts reduced by 9% on the average when phantom size was increased from head to pelvis. With GSS method, HU values reduced only by 5% under the same conditions. In summary, GSS method complements the 2D grid's scatter suppression performance, by correcting the scatter transmitted through the grid. This approach does not require additional scatter-measurement hardware, such as beam-stop arrays, since the grid itself is employed as the scatter measurement device. By suppressing residual scatter in projections, our proposed method successfully reduced artifacts caused by 2D grid's footprint, and further improved CT number accuracy.

Simulation of x-ray-induced acoustic imaging for absolute dosimetry: Accuracy of image reconstruction methods.

PURPOSE: Three-dimensional in-vivo dose verification is one of the standing challenges in radiation therapy. X-ray-induced acoustic tomography has recently been proposed as an imaging method for use in in-vivo dosimetry. The aim of this study was to investigate the accuracy of reconstructing three-dimensional (3D) absolute dose using x-ray-induced acoustic tomography. We performed this investigation using two different tomographic dose reconstruction techniques. METHODS: Two examples of 3D dose reconstruction techniques for x-ray acoustic imaging are investigated. Dose distributions are calculated for varying field sizes using a clinical treatment planning system. The induced acoustic pressure waves which are generated by the increase in temperature due to the absorption of pulsed MV x-rays are simulated using an advanced numerical modeling package for acoustic wave propagation in the time domain. Two imaging techniques, back projection and iterative time reversal, are used to reconstruct the 3D dose distribution in a water phantom with open fields. Image analysis is performed and reconstructed depth dose curves from x-ray acoustic imaging are compared to the depth dose curves calculated from the treatment planning system. Calculated field sizes from the reconstructed dose profiles by back projection and time reversal are compared to the planned field size to determine their accuracy. The iterative time reversal imaging technique is also used to reconstruct dose in an example clinical dose distribution. Image analysis of this clinical test case is performed using the gamma passing rate. In addition, gamma passing rates are used to validate the stopping criteria in the iterative time reversal method. RESULTS: Water phantom simulations showed that back projection does not adequately reconstruct the shape and intensity of the depth dose. When compared to the depth of maximum dose calculated by a treatment planning system, the maximum dose depth by back projection is shifted deeper by 55 and 75 mm for 4 × 4 cm and 10 × 10 cm field sizes, respectively. The reconstructed depth dose by iterative time reversal accurately agrees with the planned depth dose for a 4 × 4 cm field size and is shifted deeper by 12 mm for the 10 × 10 cm field size. When reconstructing field sizes, the back projection method leads to 18% and 35% larger sizes for the 4 × 4 cm and 10 × 10 cm fields, respectively, whereas the iterative time reversal method reconstructs both field sizes with < 2% error. For the clinical dose distribution, we were able to reconstruct the dose delivered by a 1 degree sub-arc with a good accuracy. The reconstructed and planned doses were compared using gamma analysis, with> 96% gamma passing rate at 3%/2 mm. CONCLUSIONS: Our results show that the 3D x-ray acoustic reconstructed dose by iterative time reversal is considerably more accurate than the dose reconstructed by back projection. Iterative time reversal imaging has a potential for use in 3D absolute dosimetry.

Effect of grid geometry on the transmission properties of 2D grids for flat detectors in CBCT.

To suppress scatter in cone beam computed tomography (CBCT), two-dimensional antiscatter grids (2D grid) have been recently proposed. In this work, we developed several grid prototypes with higher grid ratios and smaller grid pitches than previous designs, and quantified their primary and scatter transmission properties in the context of CBCT for radiation therapy. Three focused 2D grid prototypes were developed with grid ratios at 12 and 16, and grid pitches at 2 and 3 mm. Their scatter transmission properties were measured between 80-140 kVp, and benchmarked against a high performance radiographic grid (1D grid) using a Varian TrueBeam CBCT system. The effect of source-grid misalignment on the primary transmission and the improvement in contrast-to-noise ratio (CNR) were also evaluated. Changing the grid pitch from two to three mm increased the average primary transmission from 84% to 89%. Maximum scatter-to-primary ratio (SPR) with grid ratio of 12 was 0.3, and increasing the grid ratio to 16 reduced SPR by 30%. A 10 mm misalignment in 2D grid position led to a 6%-8% reduction in average primary transmission, and reduction was more pronounced for the higher grid ratio. 2D grids provided up to factor of seven lower SPR and 21% better primary transmission than the 1D grid, and their scatter transmission exhibited lower energy dependence. While 2D grids provided up to factor of 2.3 higher CNR improvement, a significant variation in CNR improvement was not observed among different grid pitch and ratios. In summary, grid ratio of 16 and grid pitch of 2 mm can keep SPRs below 0.2 even in high scatter conditions, while keeping primary transmission fractions above 80%, key benefits of the investigated 2D grids in improving image quality of CBCT. However, such grids require precise alignment in source-grid geometry during CBCT acquisitions. This study also implies that 2D grids can provide substantially better scatter suppression and primary transmission than high-performance 1D grids currently available.

A novel total variation based ring artifact suppression method for CBCT imaging with two-dimensional antiscatter grids.

PURPOSE: Two-dimensional antiscatter grids (2DASG) for cone beam computed tomography (CBCT) is a new area of research to reduce scatter intensity, and consequently improve CBCT image quality. One of the challenges in implementation of 2DASGs is their septa shadows that are impinged on the projections. If these artifacts are not corrected, they may lead to ring artifacts in CBCT images. In this work, we present a novel method to suppress ring artifacts in FPD-based CBCT images. METHODS: Briefly, our method first detects the locations of 2DASG's septa shadows in projections and then, reduces projection pixel values in septa shadows iteratively until a residual-based convergence criterion is met. To suppress the 2DASG's septa shadows, we developed a total variation minimization (TVM) formulation, referred to as adaptive-diffusive total variation minimization (adTVM), where the diffusivity regularization parameter was adaptively adjusted during each iteration based on the magnitude of the local pixel gradients. To test our method, we have acquired CBCT scans of phantoms using three 2DASG prototypes with different grid geometries. Projections were acquired with a linac mounted CBCT system, operated in offset detector geometry. These projections were then corrected in the following steps: first, projections were corrected using a gantry angle-specific gain correction map; next, projections were corrected by applying our adTVM method. CBCT images were reconstructed using FDK filtered backprojection algorithm. To evaluate adTVM's performance, pixel value statistics and contrast-to-noise ratio (CNR) were compared in CBCT images corrected with and without our adTVM method. RESULTS: Without our adTVM method, all three 2DASG prototypes introduced ring artifacts with varying intensities in CBCT images. With our method, significant reduction in ring artifacts was observed in all test cases. Standard deviation of CT numbers was reduced by 7-74% in uniform density phantom CBCT images, CNR was increased by 8-67%, and CT number accuracy of contrast objects embedded in the phantom was preserved. CONCLUSION: We propose a new method to suppress ring artifacts caused by the 2DASG's septa shadows in CBCT images. Our initial investigations indicated that adTVM method could substantially reduce such ring artifacts while preserving CT number accuracy and maintaining good spatial resolution. Therefore, our method may potentially play an important role in enabling the implementation of 2DASGs in flat panel detector based CBCT systems.

Two-dimensional antiscatter grid: A novel scatter rejection device for Cone-beam computed tomography.

PURPOSE: Scattered radiation remains to be a major cause of image quality degradation in Flat Panel Detector (FPD)-based Cone-beam computed tomography (CBCT). We have been investigating a novel two-dimensional antiscatter grid (2D-ASG) concept to reduce scatter intensity, and hence improve CBCT image quality. We present the first CBCT imaging experiments performed with the 2D-ASG prototype, and demonstrate its efficacy in improving CBCT image quality. METHODS: A 2D-ASG prototype with septa focused to x-ray source was additively manufactured from tungsten and mounted on a Varian TrueBeam CBCT system. CBCT projections of phantoms were acquired with an offset detector geometry using TrueBeam's "developer" mode. To minimize the effect of gantry flex, projections were gain corrected on angle-specific bases. CBCT images were reconstructed using a filtered backprojection algorithm and image quality improvement was quantified by measuring contrast-to-noise ratio (CNR) and CT number accuracy in images acquired with no antiscatter grid (NO-ASG), conventional one dimensional antiscatter grid (1D-ASG), and the 2D-ASG prototype. RESULTS: A significant improvement in contrast resolution was achieved using our 2D-ASG prototype compared to results of 1D-ASG and NO-ASG acquisitions. Compared to NO-ASG and 1D-ASG experiments, the CNR of material inserts improved by as much as 86% and 54% respectively. Using 2D-ASG, CT number underestimation in water equivalent material section of the phantom was reduced by up to 325 HU when compared to NO-ASG and up to 179 HU when compared to 1D-ASG. CONCLUSION: We successfully performed the first CBCT imaging experiments with a 2D-ASG prototype. 2D-ASG provided significantly higher CT number accuracy, higher CNR, and diminished scatter-induced image artifacts in qualitative evaluations. We strongly believe that utilization of a 2D-ASG may potentially lead to better soft tissue visualization in CBCT and may enable novel clinical applications that require high CT number accuracy.

Transmission characteristics of a two dimensional antiscatter grid prototype for CBCT.

AIM: High fraction of scattered radiation in cone-beam CT (CBCT) imaging degrades CT number accuracy and visualization of low contrast objects. To suppress scatter in CBCT projections, we developed a focused, two-dimensional antiscatter grid (2DASG) prototype. In this work, we report on the primary and scatter transmission characteristics of the 2DASG prototype aimed for linac mounted, offset detector geometry CBCT systems in radiation therapy, and compared its performance to a conventional one-dimensional ASG (1DASG). METHODS: The 2DASG is an array of through-holes separated by 0.1 mm septa that was fabricated from tungsten using additive manufacturing techniques. Through-holes' focusing geometry was designed for offset detector CBCT in Varian TrueBeam system. Two types of ASGs were evaluated: (a) a conventional 1DASG with a grid ratio of 10, (b) the 2DASG prototype with a grid ratio of 8.2. To assess the scatter suppression performance of both ASGs, Scatter-to-primary ratio (SPR) and scatter transmission fraction (Ts) were measured using the beam stop method. Scatter and primary intensities were modulated by varying the phantom thickness between 10 and 40 cm. Additionally, the effect of air gap and bow tie (BT) filter on SPR and Ts were evaluated. Average primary transmission fraction (TP ) and pixel specific primary transmission were also measured for both ASGs. To assess the effect of transmission characteristics on projection image signal-to-noise ratio (SNR), SNR improvement factor was calculated. Improvement in contrast to noise ratio (CNR) was demonstrated using a low contrast object. RESULTS: In comparison to 1DASG, 2DASG reduced SPRs by a factor of 3 to 6 across the range of phantom setups investigated. Ts values for 1D and 2DASGs were in the range of 21 to 29%, and 5 to 14% respectively. 2DASG continued to provide lower SPR and Ts at increased air gap and with BT filter. Tp of 1D and 2DASGs were 70.6% and 84.7% respectively. Due to the septal shadow of the 2DASG, its pixel specific primary transmission values varied between 32.5% and 99.1%. With respect to 1DASG, 2DASG provided up to factor of 1.7 more improvement in SNR across the SPR range investigated. Moreover, 2DASG provided improved visualization of low contrast objects with respect to 1DASG and NOASG setups. CONCLUSIONS: When compared to a conventional 1DASG, 2DASG prototype provided noticeably lower SPR and Ts values, indicating its superior scatter suppression performance. 2DASG also provided 19% higher average primary transmission that was attributed to the absence of interseptal spacers and optimized grid geometry. Our results indicate that the combined effect of lower scatter and higher primary transmission provided by 2DASG may potentially translate into more accurate CT numbers and improved contrast resolution in CBCT images.

Should we customize PTV expansions for BMI? Daily cone beam computerized tomography to assess organ motion in postoperative endometrial and cervical cancer patients.

AIM: A single-institution review assessing patient characteristics contributing to daily organ motion in postoperative endometrial and cervical cancer patients treated with intensity-modulated radiotherapy (IMRT). BACKGROUND: The Radiation Therapy Oncology Group has established consensus guidelines for postoperative pelvic IMRT, recommending a 7 mm margin on all three axes of the target volume. MATERIALS AND METHODS: Daily shifts on 457 radiation setups for 18 patients were recorded in the x axis (lateral), y axis (superior-inferior) and z axis (anterior-posterior); daily positions of the planning tumor volume were referenced with the initial planning scan to quantify variations. RESULTS: Of the 457 sessions, 85 (18.6%) had plan shifts of at least 7 mm in one of the three dimensions. For obese patients (body mass index [BMI] ≥ 30), 75/306 (24.5%) sessions had plan shifts ≥7 mm. Odds of having a shift ≥7 mm in any direction was greater for obese patients under both univariate (OR 4.227, 95% CI 1.235-14.466, p = 0.021) and multivariate (OR 5.000, 95% CI 1.341-18.646, p = 0.016) analyses (MVA). Under MVA, having a BMI ≥ 30 was associated with increased odds of shifts in the anterior-posterior (1.173 mm, 95% CI 0.281-2.065, p = 0.001) and lateral (2.074 mm, 95% CI 1.284-2.864, p < 0.000) directions but not in the superior-inferior axis (0.298 mm, 95% CI -0.880 to 1.475, p = 0.619) exceeding 7 mm. CONCLUSIONS: Based on these findings, the standard planned tumor volume expansion of 7 mm is less likely to account for daily treatment changes in obese patients.

Radiation dose uncertainty and correction for a mouse orthotopic and xenograft irradiation model.

PURPOSE: In animal irradiation models, reported dose can vary significantly from the actual doses delivered. We describe an effective method for in vivo dose verification. MATERIALS AND METHODS: Mice bearing commercially-available cell line or patient-derived tumor cell orthotopic or flank xenografts were irradiated using a 160 kVp, 25 mA X-ray source. Entrance dose was evaluated using optically-stimulated luminescence dosimeters (OSLD) and exit dose was assessed using radiochromic film dosimetry. RESULTS: Tumor position within the irradiation field was validated using external fiducial markers. The average entrance dose in orthotopic tumors from 10 OSLDs placed on two different animal irradiation days was 514 ± 37 cGy (range: 437-545). Exit dose measurements taken from seven radiochromic films on two separate days were 341 ± 21 cGy (a 34% attenuation). Flank tumor irradiation doses measured by OSLD were 368 ± 9 cGy compared to exit doses of 330 cGy measured by radiochromic film. CONCLUSION: Variations related to the irradiation model can lead to significant under or overdosing in vivo which can affect tumor control and/or biologic endpoints that are dose-dependent. We recommend that dose measurements be determined empirically based on the mouse model and irradiator used and dose compensation adjustments performed to ensure correct and appropriate doses.

Hedgehog signaling drives radioresistance and stroma-driven tumor repopulation in head and neck squamous cancers.

Local control and overall survival in patients with advanced head and neck squamous cell cancer (HNSCC) remains dismal. Signaling through the Hedgehog (Hh) pathway is associated with epithelial-to-mesenchymal transition, and activation of the Hh effector transcription factor Gli1 is a poor prognostic factor in this disease setting. Here, we report that increased GLI1 expression in the leading edge of HNSCC tumors is further increased by irradiation, where it contributes to therapeutic inhibition. Hh pathway blockade with cyclopamine suppressed GLI1 activation and enhanced tumor sensitivity to radiotherapy. Furthermore, radiotherapy-induced GLI1 expression was mediated in part by the mTOR/S6K1 pathway. Stroma exposed to radiotherapy promoted rapid tumor repopulation, and this effect was suppressed by Hh inhibition. Our results demonstrate that Gli1 that is upregulated at the tumor-stroma intersection in HNSCC is elevated by radiotherapy, where it contributes to stromal-mediated resistance, and that Hh inhibitors offer a rational strategy to reverse this process to sensitize HNSCC to radiotherapy.