Power spectrum analysis of the x-ray scatter signal in mammography and breast tomosynthesis projections.

PURPOSE To analyze the frequency domain characteristics of the signal in mammography images and breast tomosynthesis projections with patient tissue texture due to detected scattered x-rays. METHODS Acquisitions of x-ray projection images of 19 different patient breasts were simulated using previously acquired volumetric patient images. Acquisition of these images was performed with a dedicated breast CT prototype system, and the images were classified into voxels representing skin, adipose, and glandular tissue with a previously validated automated algorithm. The classified three dimensional images then underwent simulated mechanical compression representing that which is performed during acquisition of mammography and breast tomosynthesis images. The acquisition of projection images of each patient breast was simulated using Monte Carlo methods with each simulation resulting in two images: one of the primary (non-scattered) signal and one of the scatter signal. To analyze the scatter signal for both mammography and breast tomosynthesis, two projections images of each patient breast were simulated, one with the x-ray source positioned at 0° (mammography and central tomosynthesis projection) and at 30° (wide tomosynthesis projection). The noise power spectra (NPS) for both the scatter signal alone and the total signal (primary + scatter) for all images were obtained and the combined results of all patients analyzed. The total NPS was fit to the expected power-law relationship NPS(f) = k/f β and the results were compared with those previously published on the power spectrum characteristics of mammographic texture. The scatter signal alone was analyzed qualitatively and a power-law fit was also performed. RESULTS The mammography and tomosynthesis projections of three patient breasts were too small to analyze, so a total of 16 patient breasts were analyzed. The values of β for the total signal of the 0° projections agreed well with previously published results. As expected, the scatter power spectrum reflected a fast drop-off with increasing spatial frequency, with a reduction of four orders of magnitude by 0.1 lp/mm. The β values for the scatter signal were 6.14 and 6.39 for the 0° and 30° projections, respectively. CONCLUSIONS Although the low-frequency characteristics of scatter in mammography and breast tomosynthesis were known, a quantitative analysis of the frequency domain characteristics of this signal was needed in order to optimize previously proposed software-based x-ray scatter reduction algorithms for these imaging modalities.

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

[2]  Ingrid Reiser,et al.  Comparison of power spectra for tomosynthesis projections and reconstructed images. , 2009, Medical physics.

[3]  Yong Yu,et al.  A novel cone beam breast CT scanner: system evaluation , 2007, SPIE Medical Imaging.

[4]  Arthur E. Burgess,et al.  Mammographic structure: data preparation and spatial statistics analysis , 1999, Medical Imaging.

[5]  John J Heine,et al.  Spectral analysis of full field digital mammography data. , 2002, Medical physics.

[6]  A. Burgess,et al.  Human observer detection experiments with mammograms and power-law noise. , 2001, Medical physics.

[7]  Ioannis Sechopoulos,et al.  Dosimetric characterization of a dedicated breast computed tomography clinical prototype. , 2010, Medical physics.

[8]  D. Kopans.  Digital breast tomosynthesis. , 2008, Health devices.

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

[10]  J. Boone,et al.  Evaluation of x-ray scatter properties in a dedicated cone-beam breast CT scanner. , 2005, Medical physics.

[11]  Polad M Shikhaliev,et al.  Beam hardening artefacts in computed tomography with photon counting, charge integrating and energy weighting detectors: a simulation study , 2005, Physics in medicine and biology.

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

[13]  John M Boone,et al.  Characterizing anatomical variability in breast CT images. , 2008, Medical physics.

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

[15]  C Abbey,et al.  Statistical texture synthesis of mammographic images with super-blob lumpy backgrounds. , 1999, Optics express.

[16]  Ruola Ning,et al.  Composite modulation transfer function evaluation of a cone beam computed tomography breast imaging system , 2009 .

[17]  J. Boone,et al.  Evaluation of x-ray scatter properties in a dedicated cone-beam breast CT scanner. , 2005, Medical physics.

[18]  I. Sechopoulos,et al.  A software-based x-ray scatter correction method for breast tomosynthesis. , 2011, Medical physics.

[19]  Roberto Manduchi,et al.  Bilateral filtering for gray and color images , 1998, Sixth International Conference on Computer Vision (IEEE Cat. No.98CH36271).

[20]  Stephen J. Glick,et al.  Characterization of scatter radiation in cone beam CT mammography , 2005, SPIE Medical Imaging.

[21]  C. D'Orsi,et al.  Computation of the glandular radiation dose in digital tomosynthesis of the breast. , 2006, Medical physics.

[22]  P P Fatouros,et al.  Quality assurance in mammography: artifact analysis. , 1999, Radiographics : a review publication of the Radiological Society of North America, Inc.

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

[24]  J A Seibert,et al.  Characterization of the point spread function and modulation transfer function of scattered radiation using a digital imaging system. , 1986, Medical physics.

[25]  John M. Boone,et al.  Evaluation of scatter effects on image quality for breast tomosynthesis , 2009 .

[26]  I. Sechopoulos,et al.  Cupping artifact correction and automated classification for high-resolution dedicated breast CT images. , 2012, Medical physics.

[27]  Anita Nosratieh,et al.  Anatomical complexity in breast parenchyma and its implications for optimal breast imaging strategies. , 2012, Medical physics.

[28]  John M. Boone,et al.  Simulation of Mechanical Compression of Breast Tissue , 2007, IEEE Transactions on Biomedical Engineering.

[29]  Ruola Ning,et al.  Cone-beam CT for breast imaging: Radiation dose, breast coverage, and image quality. , 2010, AJR. American journal of roentgenology.

[30]  I. Sechopoulos,et al.  Characterization of the homogeneous tissue mixture approximation in breast imaging dosimetry. , 2012, Medical physics.

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

[32]  Ruola Ning,et al.  NPS characterization and evaluation of a cone beam CT breast imaging system. , 2009, Journal of X-ray science and technology.

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

[34]  Kristina Bliznakova,et al.  A novel simulation algorithm for soft tissue compression , 2007, Medical & Biological Engineering & Computing.

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