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.