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

Ioannis Sechopoulos, Kristina Bliznakova, Baowei Fei

Research output: Contribution to journalArticle

4 Citations (Scopus)

Abstract

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.

Original languageEnglish (US)
Article number101905
JournalMedical physics
Volume40
Issue number10
DOIs
StatePublished - Jan 1 2013
Externally publishedYes

Fingerprint

Mammography
Spectrum Analysis
Breast
X-Rays
Monte Carlo Method
Three-Dimensional Imaging
Adipose Tissue
Software
Skin

Keywords

  • Breast
  • Mammography
  • Monte Carlo
  • Noise power spectrum
  • Tomosynthesis
  • X-ray scatter

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

Power spectrum analysis of the x-ray scatter signal in mammography and breast tomosynthesis projections. / Sechopoulos, Ioannis; Bliznakova, Kristina; Fei, Baowei.

In: Medical physics, Vol. 40, No. 10, 101905, 01.01.2013.

Research output: Contribution to journalArticle

@article{e4e7f0c2e3eb4d83823c28714cb0a1ca,
title = "Power spectrum analysis of the x-ray scatter signal in mammography and breast tomosynthesis projections",
abstract = "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.",
keywords = "Breast, Mammography, Monte Carlo, Noise power spectrum, Tomosynthesis, X-ray scatter",
author = "Ioannis Sechopoulos and Kristina Bliznakova and Baowei Fei",
year = "2013",
month = "1",
day = "1",
doi = "10.1118/1.4820442",
language = "English (US)",
volume = "40",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "10",

}

TY - JOUR

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

AU - Sechopoulos, Ioannis

AU - Bliznakova, Kristina

AU - Fei, Baowei

PY - 2013/1/1

Y1 - 2013/1/1

N2 - 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.

AB - 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.

KW - Breast

KW - Mammography

KW - Monte Carlo

KW - Noise power spectrum

KW - Tomosynthesis

KW - X-ray scatter

UR - http://www.scopus.com/inward/record.url?scp=84885719144&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84885719144&partnerID=8YFLogxK

U2 - 10.1118/1.4820442

DO - 10.1118/1.4820442

M3 - Article

C2 - 24089907

AN - SCOPUS:84885719144

VL - 40

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 10

M1 - 101905

ER -