Adam Wang

Non-local total-variation (NLTV) minimization combined with reweighted L1-norm for compressed sensing CT reconstruction

The compressed sensing (CS) technique has been employed to reconstruct CT/CBCT images from fewer projections as it is designed to recover a sparse signal from highly under-sampled measurements. Since the CT image itself cannot be sparse, a variety of transforms were developed to make the image sufficiently sparse. The total-variation (TV) transform with local image gradient in L1-norm was adopted in most cases. This approach, however, which utilizes very local information and penalizes the weight at a constant rate regardless of different degrees of spatial gradient, may not produce qualified reconstructed images from noise-contaminated CT projection data. This work presents a new non-local operator of total-variation (NLTV) to overcome the deficits stated above by utilizing a more global search and non-uniform weight penalization in reconstruction. To further improve the reconstructed results, a reweighted L1-norm that approximates the ideal sparse signal recovery of the L0-norm is incorporated into the NLTV reconstruction with additional iterates. This study tested the proposed reconstruction method (reweighted NLTV) from under-sampled projections of 4 objects and 5 experiments (1 digital phantom with low and high noise scenarios, 1 pelvic CT, and 2 CBCT images). We assessed its performance against the conventional TV, NLTV and reweighted TV transforms in the tissue contrast, reconstruction accuracy, and imaging resolution by comparing contrast-noise-ratio (CNR), normalized root-mean square error (nRMSE), and profiles of the reconstructed images. Relative to the conventional NLTV, combining the reweighted L1-norm with NLTV further enhanced the CNRs by 2-4 times and improved reconstruction accuracy. Overall, except for the digital phantom with low noise simulation, our proposed algorithm produced the reconstructed image with the lowest nRMSEs and the highest CNRs for each experiment.

A novel multilayer MV imager computational model for component optimization

Purpose: A novel Megavoltage (MV) multilayer imager (MLI) design featuring higher detective quantum efficiency and lower noise than current conventional MV imagers in clinical use has been recently reported. Optimization of the MLI design for multiple applications including tumor tracking, MV‐CBCT and portal dosimetry requires a computational model that will provide insight into the physics processes that affect the overall and individual components performance. The purpose of the current work was to develop and validate a comprehensive computational model that can be used for MLI optimization. Methods: The MLI model was built using the Geant4 Application for Tomographic Emission (GATE) application. The model includes x‐ray and charged‐particle interactions as well as the optical transfer within the phosphor. A first prototype MLI device featuring a stack of four detection layers was used for model validation. Each layer of the prototype contains a copper buildup plate, a phosphor screen and photodiode array. The model was validated against measured data of Modulation Transfer Function (MTF), Noise‐Power Spectrum (NPS), and Detective Quantum Efficiency (DQE). MTF was computed using a slanted slit with 2.3° angle and 0.1 mm width. NPS was obtained using the autocorrelation function technique. DQE was calculated from MTF and NPS data. The comparison metrics between simulated and measured data were the Pearsons correlation coefficient (r) and the normalized root‐mean‐square error (NRMSE). Results: Good agreement between measured and simulated MTF and NPS values was observed. Pearsons correlation coefficient for the combined signal from all layers of the MLI was equal to 0.9991 for MTF and 0.9992 for NPS; NRMSE was 0.0121 for MTF and 0.0194 for NPS. Similarly, the DQE correlation coefficient for the combined signal was 0.9888 and the NRMSE was 0.0686. Conclusions: A comprehensive model of the novel MLI design was developed using the GATE toolkit and validated against measured MTF, NPS, and DQE data acquired with a prototype device featuring four layers. This model will be used for further optimization of the imager components and configuration for clinical radiotherapy applications.

A piecewise-focused high DQE detector for MV imaging.

PURPOSEnElectronic portal imagers (EPIDs) with high detective quantum efficiencies (DQEs) are sought to facilitate the use of the megavoltage (MV) radiotherapy treatment beam for image guidance. Potential advantages include high quality (treatment) beams eye view imaging, and improved cone-beam computed tomography (CBCT) generating images with more accurate electron density maps with immunity to metal artifacts. One approach to increasing detector sensitivity is to couple a thick pixelated scintillator array to an active matrix flat panel imager (AMFPI) incorporating amorphous silicon thin film electronics. Cadmium tungstate (CWO) has many desirable scintillation properties including good light output, a high index of refraction, high optical transparency, and reasonable cost. However, due to the 0 1 0 cleave plane inherent in its crystalline structure, the difficulty of cutting and polishing CWO has, in part, limited its study relative to other scintillators such as cesium iodide and bismuth germanate (BGO). The goal of this work was to build and test a focused large-area pixelated strip CWO detector.nnnMETHODSnA 361 × 52 mm scintillator assembly that contained a total of 28 072 pixels was constructed. The assembly comprised seven subarrays, each 15 mm thick. Six of the subarrays were fabricated from CWO with a pixel pitch of 0.784 mm, while one array was constructed from BGO for comparison. Focusing was achieved by coupling the arrays to the Varian AS1000 AMFPI through a piecewise linear arc-shaped fiber optic plate. Simulation and experimental studies of modulation transfer function (MTF) and DQE were undertaken using a 6 MV beam, and comparisons were made between the performance of the pixelated strip assembly and the most common EPID configuration comprising a 1 mm-thick copper build-up plate attached to a 133 mg/cm(2) gadolinium oxysulfide scintillator screen (Cu-GOS). Projection radiographs and CBCT images of phantoms were acquired. The work also introduces the use of a lightweight edge phantom to generate MTF measurements at MV energies and shows its functional equivalence to the more cumbersome slit-based method.nnnRESULTSnMeasured and simulated DQE(0)s of the pixelated CWO detector were 22% and 26%, respectively. The average measured and simulated ratios of CWO DQE(f) to Cu-GOS DQE(f) across the frequency range of 0.0-0.62 mm(-1) were 23 and 29, respectively. 2D and 3D imaging studies confirmed the large dose efficiency improvement and that focus was maintained across the field of view. In the CWO CBCT images, the measured spatial resolution was 7 lp/cm. The contrast-to-noise ratio was dramatically improved reflecting a 22 × sensitivity increase relative to Cu-GOS. The CWO scintillator material showed significantly higher stability and light yield than the BGO material.nnnCONCLUSIONSnAn efficient piecewise-focused pixelated strip scintillator for MV imaging is described that offers more than a 20-fold dose efficiency improvement over Cu-GOS.

Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter – Part II: System modeling, scatter correction, and optimization

PURPOSEnTo correct for scatter in kV cone-beam CT (CBCT) projection data on a clinical system using a new tool, Acuros® CTS, that estimates scatter images rapidly and accurately by deterministically solving the linear Boltzmann transport equation.nnnMETHODSnPhantom and patient CBCT scans were acquired on TrueBeam® radiotherapy machines. A first-pass reconstruction was used to create water and bone density maps of the imaged object, which was updated to include a more accurate representation of the patient couch. The imaging system model accounted for the TrueBeam x-ray source (polychromatic spectrum, beam filtration, bowtie filter, and collimation hardware) and x-ray detection system (antiscatter grid, flat-panel imager). Acuros CTS then used the system and object models to estimate the scatter component of each projection image, which was subtracted from the measured projections. The corrected projections were then reconstructed to produce the final result. We examined the tradeoff between run time and accuracy using a Pareto optimization of key parameters, including the voxel size of the down-sampled object model, the number of pixels in the down-sampled detector, and the number of scatter images (angular down-sampling). All computations and reconstructions were performed on a research workstation containing two graphics processing units (GPUs). In addition, we established a method for selecting a subset of projections for which scatter images were calculated. The projections were selected to minimize interpolation errors in the remaining projections. Image quality improvement was assessed by measuring the accuracy of the reconstructed phantom and patient images.nnnRESULTSnThe Pareto optimization yielded a set of parameters with an average run time of 26 seconds for scatter correction while maintaining high accuracy of scatter estimation. This was achieved in part by means of optimizing the projection angles that were processed, thus favoring the use of more angles in the lateral (i.e., horizontal) direction and fewer angles in the AP direction. In a 40 cm solid water phantom reconstruction, nonuniformities were decreased from 217 HU without scatter correction to 51 HU with conventional (kernel-based) scatter correction to 17 HU with Acuros CTS-based scatter correction. In clinical pelvis scans, nonuniformities in the bladder were reduced from 85 HU with conventional scatter correction to 14 HU with Acuros CTS.nnnCONCLUSIONSnAcuros CTS is a promising new tool for fast and accurate scatter correction for CBCT imaging. By carefully modeling the imaging chain and optimizing several parameters, we achieved high correction accuracies with computation times compatible with the clinical workflow. The improvement in image quality enables better soft-tissue visualization and potentially enables applications such as adaptive radiotherapy.

Asymmetric scatter kernels for software-based scatter correction of gridless mammography

Scattered radiation remains one of the primary challenges for digital mammography, resulting in decreased image contrast and visualization of key features. While anti-scatter grids are commonly used to reduce scattered radiation in digital mammography, they are an incomplete solution that can add radiation dose, cost, and complexity. Instead, a software-based scatter correction method utilizing asymmetric scatter kernels is developed and evaluated in this work, which improves upon conventional symmetric kernels by adapting to local variations in object thickness and attenuation that result from the heterogeneous nature of breast tissue. This fast adaptive scatter kernel superposition (fASKS) method was applied to mammography by generating scatter kernels specific to the object size, x-ray energy, and system geometry of the projection data. The method was first validated with Monte Carlo simulation of a statistically-defined digital breast phantom, which was followed by initial validation on phantom studies conducted on a clinical mammography system. Results from the Monte Carlo simulation demonstrate excellent agreement between the estimated and true scatter signal, resulting in accurate scatter correction and recovery of 87% of the image contrast originally lost to scatter. Additionally, the asymmetric kernel provided more accurate scatter correction than the conventional symmetric kernel, especially at the edge of the breast. Results from the phantom studies on a clinical system further validate the ability of the asymmetric kernel correction method to accurately subtract the scatter signal and improve image quality. In conclusion, software-based scatter correction for mammography is a promising alternative to hardware-based approaches such as anti-scatter grids.

Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter – Part I: Core algorithms and validation

PURPOSEnTo describe Acuros® CTS, a new software tool for rapidly and accurately estimating scatter in x-ray projection images by deterministically solving the linear Boltzmann transport equationxa0(LBTE).nnnMETHODSnThe LBTE describes the behavior of particles as they interact with an object across spatial, energy, and directional (propagation) domains. Acuros CTS deterministically solves the LBTE by modeling photon transport associated with an x-ray projection in three main steps: (a) Ray tracing photons from the x-ray source into the object where they experience their first scattering event and form scattering sources. (b) Propagating photons from their first scattering sources across the object in all directions to form second scattering sources, then repeating this process until all high-order scattering sources are computed using the source iteration method. (c) Ray-tracing photons from scattering sources within the object to the detector, accounting for the detectors energy and anti-scatter grid responses. To make this process computationally tractable, a combination of analytical and discrete methods is applied. The three domains are discretized using the Linear Discontinuous Finite Elements, Multigroup, and Discrete Ordinates methods, respectively, which confer the ability to maintain the accuracy of a continuous solution. Furthermore, through the implementation in CUDA, we sought to exploit the parallel computing capabilities of graphics processing units (GPUs) to achieve the speeds required for clinical utilization. Acuros CTS was validated against Geant4 Monte Carlo simulations using two digital phantoms: (a) a water phantom containing lung, air, and bone inserts (WLAB phantom) and (b) a pelvis phantom derived from a clinical CT dataset. For these studies, we modeled the TrueBeam® (Varian Medical Systems, Palo Alto, CA) kV imaging system with a source energy of 125xa0kVp. The imager comprised a 600xa0μm-thick Cesium Iodide (CsI) scintillator and a 10:1 one-dimensional anti-scatter grid. For the WLAB studies, the full-fan geometry without a bowtie filter was used (with and without the anti-scatter grid). For the pelvis phantom studies, a half-fan geometry with bowtie was used (with the anti-scatter grid). Scattered and primary photon fluences and energies deposited in the detector were recorded.nnnRESULTSnThe Acuros CTS and Monte Carlo results demonstrated excellent agreement. For the WLAB studies, the average percent difference between the Monte Carlo- and Acuros-generated scattered photon fluences at the face of the detector was -0.7%. After including the detector response, the average percent differences between the Monte Carlo- and Acuros-generated scatter fractions (SF) were -0.1% without the grid and 0.6% with the grid. For the digital pelvis simulation, the Monte Carlo- and Acuros-generated SFs agreed to within 0.1% on average, despite the scatter-to-primary ratios (SPRs) being as high as 5.5. The Acuros CTS computation time for each scatter image was ~1xa0s using a single GPU.nnnCONCLUSIONSnAcuros CTS enables a fast and accurate calculation of scatter images by deterministically solving the LBTE thus offering a computationally attractive alternative to Monte Carlo methods. Part II describes the application of Acuros CTS to scatter correction of CBCT scans on the TrueBeam system.

WE-AB-207A-08: BEST IN PHYSICS (IMAGING): Advanced Scatter Correction and Iterative Reconstruction for Improved Cone-Beam CT Imaging On the TrueBeam Radiotherapy Machine

PURPOSEnTo improve CBCT image quality for image-guided radiotherapy by applying advanced reconstruction algorithms to overcome scatter, noise, and artifact limitations METHODS: CBCT is used extensively for patient setup in radiotherapy. However, image quality generally falls short of diagnostic CT, limiting soft-tissue based positioning and potential applications such as adaptive radiotherapy. The conventional TrueBeam CBCT reconstructor uses a basic scatter correction and FDK reconstruction, resulting in residual scatter artifacts, suboptimal image noise characteristics, and other artifacts like cone-beam artifacts. We have developed an advanced scatter correction that uses a finite-element solver (AcurosCTS) to model the behavior of photons as they pass (and scatter) through the object. Furthermore, iterative reconstruction is applied to the scatter-corrected projections, enforcing data consistency with statistical weighting and applying an edge-preserving image regularizer to reduce image noise. The combined algorithms have been implemented on a GPU. CBCT projections from clinically operating TrueBeam systems have been used to compare image quality between the conventional and improved reconstruction methods. Planning CT images of the same patients have also been compared.nnnRESULTSnThe advanced scatter correction removes shading and inhomogeneity artifacts, reducing the scatter artifact from 99.5 HU to 13.7 HU in a typical pelvis case. Iterative reconstruction provides further benefit by reducing image noise and eliminating streak artifacts, thereby improving soft-tissue visualization. In a clinical head and pelvis CBCT, the noise was reduced by 43% and 48%, respectively, with no change in spatial resolution (assessed visually). Additional benefits include reduction of cone-beam artifacts and reduction of metal artifacts due to intrinsic downweighting of corrupted rays.nnnCONCLUSIONnThe combination of an advanced scatter correction with iterative reconstruction substantially improves CBCT image quality. It is anticipated that clinically acceptable reconstruction times will result from a multi-GPU implementation (the algorithms are under active development and not yet commercially available). All authors are employees of and (may) own stock of Varian Medical Systems.

TH-EF-BRA-02: Patient-Specific Dose Maps for CT Scans Using a Fast, Deterministic Boltzmann Transport Equation Solver

Purpose: To develop a rapid and accurate software tool for computing patient-specific radiation dose maps of dose delivered from kV computed tomography (CT) scans. Methods: Monte Carlo methods currently provide the gold-standard for calculating patient-specific dose maps, but require immense computational resources to achieve sufficiently high statistical accuracy. To overcome this limitation, a deterministic method was implemented to solve the same underlying Boltzmann transport equation (BTE) that governs particle interactions and transport. Phase-space was discretized according to spatial location, energy, and angle, and a deterministic finite element algorithm was applied to compute the object’s photon fluence distribution, which does not exhibit stochastic noise. A computationally efficient GPU implementation for a standard workstation was developed, and comparison was made between the performance of the deterministic BTE solver and a standard Monte Carlo algorithm for a cone-beam projection of a virtual anthropomorphic chest phantom. Results: The BTE solution and Monte Carlo results were in strong agreement with a relative root-mean square error (RMSE) of 3.47%. Some larger differences existed at high-contrast boundaries (e.g., air/water) and within the bone, and are under further investigation. Notably, the computation time of the BTE solver was 8 seconds, while to obtain the same level of statistical uncertainty with conventional Monte Carlo required 1200 CPU-hours. Additionally, unlike Monte Carlo, the BTE computation time is only weakly dependent on the number of sources, making it extremely well-suited for CT dose calculations. Therefore, the BTE-based method is expected to offer a >30,000x speed increase compared to Monte Carlo for entire CT scans, even after application of variance reduction techniques and GPU implementation. Conclusion: The novel deterministic BTE solver offers a significantly faster alternative to Monte Carlo-based methods for computing dose delivered by CT scans, which can enable estimation of patient-specific organ doses for each CT examination performed. Adam Wang, Alex Maslowski, Todd Wareing, and Josh Star-Lack are employees of Varian Medical Systems.

Multi-layer imager design for mega-voltage spectral imaging

The architecture of multi-layer imagers (MLIs) can be exploited to provide megavoltage spectral imaging (MVSPI) for specific imaging tasks. In the current work, we investigated bone suppression and gold fiducial contrast enhancement as two clinical tasks which could be improved with spectral imaging. A method based on analytical calculations that enables rapid investigation of MLI component materials and thicknesses was developed and validated against Monte Carlo computations. The figure of merit for task-specific imaging performance was the contrast-to-noise ratio (CNR) of the gold fiducial when the CNR of bone was equal to zero after a weighted subtraction of the signals obtained from each MLI layer. Results demonstrated a sharp increase in the CNR of gold when the build-up component or scintillation materials and thicknesses were modified. The potential for low-cost, prompt implementation of specific modifications (e.g. composition of the build-up component) could accelerate clinical translation of MVSPI.

Physics considerations in MV-CBCT multi-layer imager design

Megavoltage (MV) cone-beam computed tomography (CBCT) using an electronic portal imaging (EPID) offers advantageous features, including 3D mapping, treatment beam registration, high-z artifact suppression, and direct radiation dose calculation. Adoption has been slowed by image quality limitations and concerns about imaging dose. Developments in imager design, including pixelated scintillators, structured phosphors, inexpensive scintillation materials, and multi-layer imager (MLI) architecture have been explored to improve EPID image quality and reduce imaging dose. The present study employs a hybrid Monte Carlo and linear systems model to determine the effect of detector design elements, such as multi-layer architecture and scintillation materials. We follow metrics of image quality including modulation transfer function (MTF) and noise power spectrum (NPS) from projection images to 3D reconstructions to in-plane slices and apply a task based figure-of-merit, the ideal observer signal-to-noise ratio (d) to determine the effect of detector design on object detectability. Generally, detectability was limited by detector noise performance. Deploying an MLI imager with a single scintillation material for all layers yields improvement in noise performance and d linear with the number of layers. In general, improving x-ray absorption using thicker scintillators results in improved DQE(0). However, if light yield is low, performance will be affected by electronic noise at relatively high doses, resulting in rapid image quality degradation. Maximizing image quality in a heterogenous MLI detector (i.e. multiple different scintillation materials) is most affected by limiting total noise. However, while a second-order effect, maximizing total spatial resolution of the MLI detector is a balance between the intensity contribution of each layer against its individual MTF. So, while a thinner scintillator may yield a maximal individual-layer MTF, its quantum efficiency will be relatively low in comparison to a thicker scintillator and thus, intensity contribution may be insufficient to noticeably improve the total detector MTF.