Josh Star-Lack

Improved scatter correction using adaptive scatter kernel superposition

Accurate scatter correction is required to produce high-quality reconstructions of x-ray cone-beam computed tomography (CBCT) scans. This paper describes new scatter kernel superposition (SKS) algorithms for deconvolving scatter from projection data. The algorithms are designed to improve upon the conventional approach whose accuracy is limited by the use of symmetric kernels that characterize the scatter properties of uniform slabs. To model scatter transport in more realistic objects, nonstationary kernels, whose shapes adapt to local thickness variations in the projection data, are proposed. Two methods are introduced: (1) adaptive scatter kernel superposition (ASKS) requiring spatial domain convolutions and (2) fast adaptive scatter kernel superposition (fASKS) where, through a linearity approximation, convolution is efficiently performed in Fourier space. The conventional SKS algorithm, ASKS, and fASKS, were tested with Monte Carlo simulations and with phantom data acquired on a table-top CBCT system matching the Varian On-Board Imager (OBI). All three models accounted for scatter point-spread broadening due to object thickening, object edge effects, detector scatter properties and an anti-scatter grid. Hounsfield unit (HU) errors in reconstructions of a large pelvis phantom with a measured maximum scatter-to-primary ratio over 200% were reduced from -90 ± 58 HU (mean ± standard deviation) with no scatter correction to 53 ± 82 HU with SKS, to 19 ± 25 HU with fASKS and to 13 ± 21 HU with ASKS. HU accuracies and measured contrast were similarly improved in reconstructions of a body-sized elliptical Catphan phantom. The results show that the adaptive SKS methods offer significant advantages over the conventional scatter deconvolution technique.

Scatter correction method for x-ray CT using primary modulation: Phantom studies

PURPOSE Scatter correction is a major challenge in x-ray imaging using large area detectors. Recently, the authors proposed a promising scatter correction method for x-ray computed tomography (CT) using primary modulation. Proof of concept was previously illustrated by Monte Carlo simulations and physical experiments on a small phantom with a simple geometry. In this work, the authors provide a quantitative evaluation of the primary modulation technique and demonstrate its performance in applications where scatter correction is more challenging. METHODS The authors first analyze the potential errors of the estimated scatter in the primary modulation method. On two tabletop CT systems, the method is investigated using three phantoms: A Catphan 600 phantom, an anthropomorphic chest phantom, and the Catphan 600 phantom with two annuli. Two different primary modulators are also designed to show the impact of the modulator parameters on the scatter correction efficiency. The first is an aluminum modulator with a weak modulation and a low modulation frequency, and the second is a copper modulator with a strong modulation and a high modulation frequency. RESULTS On the Catphan 600 phantom in the first study, the method reduces the error of the CT number in the selected regions of interest (ROIs) from 371.4 to 21.9 Hounsfield units (HU); the contrast to noise ratio also increases from 10.9 to 19.2. On the anthropomorphic chest phantom in the second study, which represents a more difficult case due to the high scatter signals and object heterogeneity, the method reduces the error of the CT number from 327 to 19 HU in the selected ROIs and from 31.4% to 5.7% on the overall average. The third study is to investigate the impact of object size on the efficiency of our method. The scatter-to-primary ratio estimation error on the Catphan 600 phantom without any annulus (20 cm in diameter) is at the level of 0.04, it rises to 0.07 and 0.1 on the phantom with an elliptical annulus (30 cm in the minor axis and 38 cm in the major axis) and with a circular annulus (38 cm in diameter). CONCLUSIONS On the three phantom studies, good scatter correction performance of the proposed method has been demonstrated using both image comparisons and quantitative analysis. The theory and experiments demonstrate that a strong primary modulation that possesses a low transmission factor and a high modulation frequency is preferred for high scatter correction accuracy.

Scanning-beam digital x-ray (SBDX) technology for interventional and diagnostic cardiac angiography

The scanning-beam digital x-ray (SBDX) system is designed for x-ray dose reduction in cardiac angiographic applications. Scatter reduction, efficient detection of primary x-rays, and an inverse beam geometry are the main components of the entrance dose reduction strategy. This paper reports the construction of an SBDX prototype, image reconstruction techniques, and measurements of spatial resolution and x-ray output. The x-ray source has a focal spot that is electronically scanned across a large-area transmission target. A multihole collimator beyond the target defines a series of x-ray beams directed at a distant small-area detector array. The prototype has a 23 cm X 23 cm target, 100 X 100 focal spot positions, and a 5 cm X 5 cm CdTe detector positioned 150 cm from the target. With this nonmechanical method of beam scanning, patient images with low detected scatter are generated at up to 30 frame/s. SBDX data acquisition is tomosynthetic. The prototype simultaneously reconstructs 16 planes spaced throughout the cardiac volume using shift-and-add backprojection. Image frames analogous to conventional projection images are generated with a multiplane compositing algorithm. Single-plane versus multiplane reconstruction of contrast-filled coronary arteries is demonstrated with images of the porcine heart. Phantom and porcine imaging studies show multiplane reconstruction is practicable under clinically realistic levels of patient attenuation and cardiac motion. The modulation transfer function for an in-plane slit at mechanical isocenter measured 0.41-0.56 at 1 cycle/mm, depending on the detector element to image pixel interpolation technique. Modeling indicates that desired gains in spatial resolution are achievable by halving the detector element width. The x-ray exposure rate 15 cm below isocenter, without table or patient in the beam, measured 11.5 R/min at 120 kVp, 24.3 kWp and 3.42 R/min at 70 kVp, 14.2 kWp.

Efficient scatter correction using asymmetric kernels

X-ray cone-beam (CB) projection data often contain high amounts of scattered radiation, which must be properly modeled in order to produce accurate computed tomography (CT) reconstructions. A well known correction technique is the scatter kernel superposition (SKS) method that involves deconvolving projection data with kernels derived from pencil beam-generated scatter point-spread functions. The method has the advantages of being practical and computationally efficient but can suffer from inaccuracies. We show that the accuracy of the SKS algorithm can be significantly improved by replacing the symmetric kernels that traditionally have been used with nonstationary asymmetric kernels. We also show these kernels can be well approximated by combinations of stationary kernels thus allowing for efficient implementation of convolution via FFT. To test the new algorithm, Monte Carlo simulations and phantom experiments were performed using a table-top system with geometry and components matching those of the Varian On-Board Imager (OBI). The results show that asymmetric kernels produced substantially improved scatter estimates. For large objects with scatter-to-primary ratios up to 2.0, scatter profiles were estimated to within 10% of measured values. With all corrections applied, including beam hardening and lag, the resulting accuracies of the CBCT reconstructions were within ±25 Hounsfield Units (±2.5%).

Rapid Monte Carlo simulation of detector DQE(f)

PURPOSE Performance optimization of indirect x-ray detectors requires proper characterization of both ionizing (gamma) and optical photon transport in a heterogeneous medium. As the tool of choice for modeling detector physics, Monte Carlo methods have failed to gain traction as a design utility, due mostly to excessive simulation times and a lack of convenient simulation packages. The most important figure-of-merit in assessing detector performance is the detective quantum efficiency (DQE), for which most of the computational burden has traditionally been associated with the determination of the noise power spectrum (NPS) from an ensemble of flood images, each conventionally having 10(7) - 10(9) detected gamma photons. In this work, the authors show that the idealized conditions inherent in a numerical simulation allow for a dramatic reduction in the number of gamma and optical photons required to accurately predict the NPS. METHODS The authors derived an expression for the mean squared error (MSE) of a simulated NPS when computed using the International Electrotechnical Commission-recommended technique based on taking the 2D Fourier transform of flood images. It is shown that the MSE is inversely proportional to the number of flood images, and is independent of the input fluence provided that the input fluence is above a minimal value that avoids biasing the estimate. The authors then propose to further lower the input fluence so that each event creates a point-spread function rather than a flood field. The authors use this finding as the foundation for a novel algorithm in which the characteristic MTF(f), NPS(f), and DQE(f) curves are simultaneously generated from the results of a single run. The authors also investigate lowering the number of optical photons used in a scintillator simulation to further increase efficiency. Simulation results are compared with measurements performed on a Varian AS1000 portal imager, and with a previously published simulation performed using clinical fluence levels. RESULTS On the order of only 10-100 gamma photons per flood image were required to be detected to avoid biasing the NPS estimate. This allowed for a factor of 10(7) reduction in fluence compared to clinical levels with no loss of accuracy. An optimal signal-to-noise ratio (SNR) was achieved by increasing the number of flood images from a typical value of 100 up to 500, thereby illustrating the importance of flood image quantity over the number of gammas per flood. For the point-spread ensemble technique, an additional 2× reduction in the number of incident gammas was realized. As a result, when modeling gamma transport in a thick pixelated array, the simulation time was reduced from 2.5 × 10(6) CPU min if using clinical fluence levels to 3.1 CPU min if using optimized fluence levels while also producing a higher SNR. The AS1000 DQE(f) simulation entailing both optical and radiative transport matched experimental results to within 11%, and required 14.5 min to complete on a single CPU. CONCLUSIONS The authors demonstrate the feasibility of accurately modeling x-ray detector DQE(f) with completion times on the order of several minutes using a single CPU. Convenience of simulation can be achieved using GEANT4 which offers both gamma and optical photon transport capabilities.

Inverse-geometry volumetric CT system with multiple detector arrays for wide field-of-view imaging

Current volumetric computed tomography (CT) methods require seconds to acquire a thick volume (>8 cm) with high resolution. Inverse-geometry CT (IGCT) is a new system geometry under investigation that is anticipated to be able to image a thick volume in a single gantry rotation with isotropic resolution and no cone-beam artifacts. IGCT employs a large array of source spots opposite a smaller detector array. The in-plane field of view (FOV) is primarily determined by the size of the source array, in much the same way that the FOV is determined by the size of the detector array in a conventional CT system. Thus, the size of the source array can be a limitation on the achievable FOV. We propose adding additional detector arrays, spaced apart laterally, to increase the in-plane FOV while still using a modestly sized source array. We determine optimal detector placement to maximize the FOV while obtaining relatively uniform sampling. We also demonstrate low wasted radiation of the proposed system through design and simulation of a pre-patient collimator. Reconstructions from simulated projection data show no artifacts when combining the data from the detector arrays. Finally, to demonstrate feasibility of the concept, an anthropomorphic thorax phantom containing a porcine heart was scanned on a prototype table-top system. The reconstructed axial images demonstrate a 45 cm in-plane FOV using a 23 cm source array.

Comparison of entrance exposure and signal-to-noise ratio between an SBDX prototype and a wide-beam cardiac angiographic system

The scanning-beam digital x-ray (SBDX) system uses an inverse geometry, narrow x-ray beam, and a 2-mm thick CdTe detector to improve the dose efficiency of the coronary angiographic procedure. Entrance exposure and large-area iodine signal-to-noise ratio (SNR) were measured with the SBDX prototype and compared to that of a clinical cardiac interventional system with image intensifier (II) and charge coupled device (CCD) camera (Philips H5000, MRC-200 x-ray tube, 72 kWp max). Phantoms were 18.6-35.0 cm acrylic with an iohexol-equivalent disk placed at midthickness (35 mg/cm2 iodine radiographic density). Imaging was performed at 15 frame/s, with the disk at mechanical isocenter and an 11-cm object-plane field width. The II/CCD system was operated in cine mode with automatic exposure control. With the SBDX prototype at maximum x-ray output (120 kVp, 24.3 kWp), the SBDX SNR was 107%-69% of the II/CCD SNR, depending on phantom thickness, and the SBDX entrance exposure rate was 10.7-9.3 R/min (9.4-8.2 cGy/min air kerma). For phantoms where an equal-kVp imaging comparison was possible (> or = 23.3 cm), the SBDX SNR ranged from 47% to 69% of the II/CCD SNR while delivering 6% to 9% of the II/CCD entrance exposure rate. From these measurements it was determined that the relative SBDX entrance exposure at equal SNR would be 31%-16%. Results were consistent with a model for relative entrance exposure at equal SNR, which predicted a 3-7 times reduction in entrance exposure due to SBDXs comparatively low scatter fraction (5.5%-8.1% measured, including off-focus radiation), high detector detective quantum efficiency (66%-73%, measured from 70 to 120 kVp), and large entrance field area (1.7x - 2.3x, for the same object-plane field width). With improvements to the system geometry, detector, and x-ray source, SBDX technology is projected to achieve conventional cine-quality SNR over a full range of patient thicknesses, with 5-10 times lower skin dose.

A practical reconstruction algorithm for CT noise variance maps using FBP reconstruction

Accurate prediction of reconstructed noise in computed tomography (CT) images is important for purposes of system design, optimization and evaluation. A large body of work describes noise prediction methods for CT, the vast majority of which assume stationarity of both noise and signal processes. Consequently, these methods are usually applied to and evaluated using simple phantoms, and only a portion of the image is scrutinized. In this work, we derive a practical method for reconstructing CT noise variance maps for arbitrary objects. Photon Poisson noise and system electronic noise are considered. The final formula has the same structure as that of the filtered backprojection (FBP) formula, but with different weighting factors and convolution kernels. The algorithm is verified using computer simulations of the Shepp-Logan phantom, and a good match is found between the predicted noise map from one single noisy scan and the measured noise using 128 noisy scans. As compared to other proposed noise models, our complementary work provides a method of noise prediction by simple adaptation of FBP reconstruction algorithms. The result is a tool that can be useful for system optimization and evaluation tasks as well as the design of reconstruction filters.

A novel EPID design for enhanced contrast and detective quantum efficiency

Beams-eye-view imaging applications such as real-time soft-tissue motion estimation are hindered by the inherently low image contrast of electronic portal imaging devices (EPID) currently available for clinical use. We introduce and characterize a novel EPID design that provides substantially increased detective quantum efficiency (DQE), contrast-to-noise ratio (CNR) and sensitivity without degradation in spatial resolution. The prototype design features a stack of four conventional EPID layers combined with low noise integrated readout electronics. Each layer consists of a copper plate, a scintillator ([Formula: see text]) and a photodiode/TFT-switch (aSi:H). We characterize the prototypes signal response to a 6 MV photon beam in terms of modulation transfer function (MTF), DQE and CNR. The presampled MTF is estimated using a slanted slit technique, the DQE is calculated from measured normalized noise power spectra (nNPS) and the MTF and CNR is estimated using a Las Vegas contrast phantom. The prototype has been designed and built to be interchangeable with the current clinical EPID on the Varian TrueBeam platform (AS-1200) in terms of size and data output specifications. Performance evaluation is conducted in absolute values as well as in relative terms using the Varian AS-1200 EPID as a reference detector. A fivefold increase of DQE(0) to about 6.7% was observed by using the four-layered design versus the AS-1200 reference detector. No substantial differences are observed between each layers individual MTF and the one for all four layers operating combined indicating that defocusing due to beam divergence is negligible. Also, using four layers instead of one increases the signal to noise ratio by a factor of 1.7.

Metal artifact correction for x-ray computed tomography using kV and selective MV imaging.

PURPOSE The overall goal of this work is to improve the computed tomography (CT) image quality for patients with metal implants or fillings by completing the missing kilovoltage (kV) projection data with selectively acquired megavoltage (MV) data that do not suffer from photon starvation. When both of these imaging systems, which are available on current radiotherapy devices, are used, metal streak artifacts are avoided, and the soft-tissue contrast is restored, even for regions in which the kV data cannot contribute any information. METHODS Three image-reconstruction methods, including two filtered back-projection (FBP)-based analytic methods and one iterative method, for combining kV and MV projection data from the two on-board imaging systems of a radiotherapy device are presented in this work. The analytic reconstruction methods modify the MV data based on the information in the projection or image domains and then patch the data onto the kV projections for a FBP reconstruction. In the iterative reconstruction, the authors used dual-energy (DE) penalized weighted least-squares (PWLS) methods to simultaneously combine the kV/MV data and perform the reconstruction. RESULTS The authors compared kV/MV reconstructions to kV-only reconstructions using a dental phantom with fillings and a hip-implant numerical phantom. Simulation results indicated that dual-energy sinogram patch FBP and the modified dual-energy PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in the kV projections. The root-mean-square errors of soft-tissue patterns obtained using combined kV/MV data are 10-15 Hounsfield units smaller than those of the kV-only images, and the structural similarity index measure also indicates a 5%-10% improvement in the image quality. The added dose from the MV scan is much less than the dose from the kV scan if a high efficiency MV detector is assumed. CONCLUSIONS The authors have shown that it is possible to improve the image quality of kV CTs for patients with metal implants or fillings by completing the missing kV projection data with selectively acquired MV data that do not suffer from photon starvation. Numerical simulations demonstrated that dual-energy sinogram patch FBP and a modified kV/MV PWLS method can successfully suppress metal streak artifacts and restore information lost due to photon starvation in kV projections. Combined kV/MV images may permit the improved delineation of structures of interest in CT images for patients with metal implants or fillings.