Estimation of scattered radiation in digital breast tomosynthesis

Digital breast tomosynthesis (DBT) is a promising technique to overcome the tissue superposition limitations found in planar 2D x-ray mammography. However, as most DBT systems do not employ an anti-scatter grid, the levels of scattered radiation recorded within the image receptor are significantly higher than that observed in planar 2D x-ray mammography. Knowledge of this field is necessary as part of any correction scheme and for computer modelling and optimisation of this examination. Monte Carlo (MC) simulations are often used for this purpose, however they are computationally expensive and a more rapid method of calculation is desirable. This issue is addressed in this work by the development of a fast kernel-based methodology for scatter field estimation using a detailed realistic DBT geometry. Thickness-dependent scatter kernels, which were validated against the literature with a maximum discrepancy of 4% for an idealised geometry, have been calculated and a new physical parameter (air gap distance) was used to estimate more accurately the distribution of scattered radiation for a series of anthropomorphic breast phantom models. The proposed methodology considers, for the first time, the effects of scattered radiation from the compression paddle and breast support plate, which can represent more than 30% of the total scattered radiation recorded within the image receptor. The results show that the scatter field estimator can calculate scattered radiation images in an average of 80 min for projection angles up to 25° with equal to or less than a 10% error across most of the breast area when compared with direct MC simulations.

[1]  John M Boone,et al.  Evaluation of scatter effects on image quality for breast tomosynthesis , 2007, SPIE Medical Imaging.

[2]  J A Seibert,et al.  X-ray scatter removal by deconvolution. , 1988, Medical physics.

[3]  G Panayiotakis,et al.  A multiple projection method for digital tomosynthesis. , 1992, Medical physics.

[4]  J A Seibert,et al.  An edge spread technique for measurement of the scatter-to-primary ratio in mammography. , 2000, Medical physics.

[5]  D M Cunha,et al.  Evaluation of scatter-to-primary ratio, grid performance and normalized average glandular dose in mammography by Monte Carlo simulation including interference and energy broadening effects , 2010, Physics in medicine and biology.

[6]  O. Díaz Montesdeoca,et al.  Scattered radiation in projection X-raymammography and digital breast tomosynthesis , 2013 .

[7]  S. Incerti,et al.  Geant4 developments and applications , 2006, IEEE Transactions on Nuclear Science.

[8]  Ruola Ning,et al.  Why should breast tumour detection go three dimensional? , 2003, Physics in medicine and biology.

[9]  D. Kopans,et al.  Digital tomosynthesis in breast imaging. , 1997, Radiology.

[10]  D R Dance,et al.  The computation of scatter in mammography by Monte Carlo methods. , 1984, Physics in medicine and biology.

[11]  James A. Scott Photon, Electron, Proton and Neutron Interaction Data for Body Tissues ICRU Report 46. International Commission on Radiation Units and Measurements, Bethesda, 1992, $40.00 , 1993 .

[12]  Olivier Alonzo-Proulx,et al.  Effect of Tissue Thickness Variation in Volumetric Breast Density Estimation , 2008, Digital Mammography / IWDM.

[13]  R. Siddon Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[14]  Andrew D. A. Maidment,et al.  Development and characterization of an anthropomorphic breast software phantom based upon region-growing algorithm. , 2011, Medical physics.

[15]  T. R. Fewell,et al.  Molybdenum, rhodium, and tungsten anode spectral models using interpolating polynomials with application to mammography. , 1997, Medical physics.

[16]  Gyula Faigel,et al.  X-Ray Holography , 1999 .

[17]  C. Mathers,et al.  Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008 , 2010, International journal of cancer.

[18]  J M Boone,et al.  Scatter/primary in mammography: Monte Carlo validation. , 2000, Medical physics.

[19]  J. Baker,et al.  Breast tomosynthesis: state-of-the-art and review of the literature. , 2011, Academic radiology.

[20]  J. Boone,et al.  Scatter/primary in mammography: comprehensive results. , 2000, Medical physics.

[21]  Federica Zanca,et al.  Two-view and single-view tomosynthesis versus full-field digital mammography: high-resolution X-ray imaging observer study. , 2012, Radiology.

[22]  R. Kruger,et al.  Scatter estimation for a digital radiographic system using convolution filtering. , 1987, Medical physics.

[23]  I. Sechopoulos A review of breast tomosynthesis. Part I. The image acquisition process. , 2013, Medical physics.

[24]  S. Molloi,et al.  Scatter correction in digital mammography based on image deconvolution , 2010, Physics in medicine and biology.

[25]  James T Dobbins,et al.  Digital x-ray tomosynthesis: current state of the art and clinical potential. , 2003, Physics in medicine and biology.

[26]  John M Boone,et al.  Grid and slot scan scatter reduction in mammography: comparison by using Monte Carlo techniques. , 2002, Radiology.

[27]  Andriy I. Bandos,et al.  Comparison of digital mammography alone and digital mammography plus tomosynthesis in a population-based screening program. , 2013, Radiology.

[28]  C. D'Orsi,et al.  Scatter radiation in digital tomosynthesis of the breast. , 2007, Medical physics.

[29]  Kenneth C. Young,et al.  A fast scatter field estimator for digital breast tomosynthesis , 2012, Medical Imaging.

[30]  E. W. Shrigley Medical Physics , 1944, British medical journal.

[31]  H. O. Wyckoff,et al.  The International Commission on Radiation Units and Measurements , 2001, Journal of the ICRU.

[32]  J. S. Laughlin,et al.  Absorbed radiation dose in mammography. , 1979, Radiology.

[33]  I. Sechopoulos A review of breast tomosynthesis. Part II. Image reconstruction, processing and analysis, and advanced applications. , 2013, Medical physics.