Quick search Find article
Quick search
Find article

Rayleigh scatter in kilovoltage x-ray imaging: is the independent atom approximation good enough?

G Poludniowski, P M Evans and S Webb

Show affiliations


Monte Carlo simulation is the gold standard method for modelling scattering processes in medical x-ray imaging. General-purpose Monte Carlo codes, however, typically use the independent atom approximation (IAA). This is known to be inaccurate for Rayleigh scattering, for many materials, in the forward direction. This work addresses whether the IAA is sufficient for the typical modelling tasks in medical kilovoltage x-ray imaging. As a means of comparison, we incorporate a more realistic 'interference function' model into a custom-written Monte Carlo code. First, we conduct simulations of scatter from isolated voxels of soft tissue, adipose, cortical bone and spongiosa. Then, we simulate scatter profiles from a cylinder of water and from phantoms of a patient's head, thorax and pelvis, constructed from diagnostic-quality CT data sets. Lastly, we reconstruct CT numbers from simulated sets of projection images and investigate the quantitative effects of the approximation. We show that the IAA can produce errors of several per cent of the total scatter, across a projection image, for typical x-ray beams and patients. The errors in reconstructed CT number, however, for the phantoms simulated, were small (typically < 10 HU). The IAA can therefore be considered sufficient for the modelling of scatter correction in CT imaging. Where accurate quantitative estimates of scatter in individual projection images are required, however, the appropriate interference functions should be included.


PACS

87.59.bd Computed radiography

87.55.K- Monte Carlo methods

87.55.-x Treatment strategy

Subjects

Medical physics

Dates

Issue 22 (21 November 2009)

Received 5 August 2009, in final form 10 October 2009

Published 4 November 2009



  1. Rayleigh scatter in kilovoltage x-ray imaging: is the independent atom approximation good enough?

    G Poludniowski et al 2009 Phys. Med. Biol. 54 6931

  2. Relation between polymer and Fock excitations

    Abhay Ashtekar and Jerzy Lewandowski 2001 Class. Quantum Grav. 18 L117

  3. Fundamentals of zinc oxide as a semiconductor

    Anderson Janotti and Chris G Van de Walle 2009 Rep. Prog. Phys. 72 126501

  4. Magneto-dilatonic Bianchi-I cosmology: isotropization and singularity problems

    K A Bronnikov et al 2004 Class. Quantum Grav. 21 3389

  5. Error analysis of exponential integrators for oscillatory second-order differential equations

    Volker Grimm and Marlis Hochbruck 2006 J. Phys. A: Math. Gen. 39 5495

  6. Completely localized gravity with higher curvature terms

    Ishwaree P Neupane 2002 Class. Quantum Grav. 19 5507

  7. The influence of magnetic islands on drift mode stability in magnetized plasma

    H R Wilson and J W Connor 2009 Plasma Phys. Control. Fusion 51 115007

  8. Luminous Infrared Galaxies with the Submillimeter Array. I. Survey Overview and the Central Gas to Dust Ratio

    Christine D. Wilson et al. 2008 ApJS 178 189

  9. A canonical reduction of order for the Kepler problem

    S E Godfrey and G E Prince 1991 J. Phys. A: Math. Gen. 24 5465

  10. Combining Strong and Weak Gravitational Lensing in Abell 1689

    Marceau Limousin et al. 2007 ApJ 668 643

Related review articles

What's this?
View review articles related to this research to gain an insight into the key trends in this subject area. Related review articles are selected based on PACS/MSC codes, and are no more than three years old.

  1. Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction?

View by subject




Export








Please login to access our web services, or create an account if you don't yet have one.

You must have cookies enabled in your web browser to be able to login.

Username
Password

Forgotten your password? Get a new one here.