Embodiments disclose a method performed by at least one processor for processing a plurality of x-ray projection images of a subject, the method comprising a plurality of operations including reconstructing the projection images to yield a volume reconstruction； segmenting the reconstructed volume by assigning a material type to each voxel； estimating a first set of scatter images corresponding to a subset of the projection images by calculating probabilistic predictions of interactions of x-rays with the subject and applying a low pass spatial filter to the scatter images； estimating a second set of scatter images corresponding to projection images not included in the subset, based on the first set of scatter images； and subtracting, for each projection image, the corresponding scatter image to yield a corrected projection image. Applying the plurality of operations of the method in an iterative process on graphics-processor-quality computational hardware enables real-time scatter correction and reconstruction of CT volume images.
