Purpose: X-ray scattering leads to CT images with a reduced contrast, inaccurate CT values as well as streak and cupping artifacts. Therefore, scatter correction is crucial to maintain the diagnostic value of CT and CBCT examinations. However, existing approaches are not able to combine both high accuracy and high computational performance. Therefore, we propose the deep scatter estimation (DSE): a deep convolutional neural network to derive highly accurate scatter estimates in real time. Methods: Gold standard scatter estimation approaches rely on dedicated Monte Carlo (MC) photon transport codes. However, being computationally expensive, MC methods cannot be used routinely. To enable real-time scatter correction with similar accuracy, DSE uses a deep convolutional neural network that is trained to predict MC scatter estimates based on the acquired projection data. Here, the potential of DSE is demonstrated using simulations of CBCT head, thorax, and abdomen scans as well as measurements at an experimental table-top CBCT. Two conventional computationally efficient scatter estimation approaches were implemented as reference: a kernel-based scatter estimation (KSE) and the hybrid scatter estimation (HSE). Results: The simulation study demonstrates that DSE generalizes well to varying tube voltages, varying noise levels as well as varying anatomical regions as long as they are appropriately represented within the training data. In any case the deviation of the scatter estimates from the ground truth MC scatter distribution is less than 1.8% while it is between 6.2% and 293.3% for HSE and between 11.2% and 20.5% for KSE. To evaluate the performance for real data, measurements of an anthropomorphic head phantom were performed. Errors were quantified by a comparison to a slit scan reconstruction. Here, the deviation is 278 HU (no correction), 123 HU (KSE), 65 HU (HSE), and 6 HU (DSE), respectively. Conclusions: The DSE clearly outperforms conventional scatter estimation approaches in terms of accuracy. DSE is nearly as accurate as Monte Carlo simulations but is superior in terms of speed (%10 ms/projection) by orders of magnitude.

X-ray scatter is a major cause of image quality degradation in dimensional CT. Especially, in case of highly attenuating components scatter-to-primary ratios may easily be higher than 1. The corresponding artifacts which appear as cupping or dark streaks in the CT reconstruction may impair a metrological assessment. Therefore, an appropriate scatter correction is crucial. Thereby, the gold standard is to predict the scatter distribution using a Monte Carlo (MC) code and subtract the corresponding scatter estimate from the measured raw data. MC, however, is too slow to be used routinely. To correct for scatter in real-time, we developed the deep scatter estimation (DSE). It uses a deep convolutional neural network which is trained to reproduce the output of MC simulations using only the acquired projection data as input. Once trained, DSE can be applied in real-time. The present study demonstrates the potential of the proposed approach using simulations and measurements. In both cases the DSE yields highly accurate scatter estimates that differ by < 3% from our MC scatter predictions. Further, DSE clearly outperforms kernel-based scatter estimation techniques and hybrid approaches, as they are in use today.

Purpose: Scattered radiation is one of the major problems facing image quality in flat detector conebeam computed tomography (CBCT). Previously, a new scatter estimation and correction method using primary beam modulation has been proposed. The original image processing technique used a frequency-domain-based analysis, which proved to be sensitive to the accuracy of the modulator pattern both spatially and in amplitude as well as to the frequency of the modulation pattern. In addition, it cannot account for penumbra effects that occur, for example, due to the finite focal spot size and the scatter estimate can be degraded by high-frequency components of the primary image. Methods: In this paper, the authors present a new way to estimate the scatter using primary modulation. It is less sensitive to modulator nonidealities and most importantly can handle arbitrary modulator shapes and changes in modulator attenuation. The main idea is that the scatter estimation can be expressed as an optimization problem, which yields a separation of the scatter and the primary image. The method is evaluated using simulated and experimental CBCT data. The scattering properties of the modulator itself are analyzed using a Monte Carlo simulation.