Quantification and normalization of x-ray mammograms.

The analysis of (x-ray) mammograms remains qualitative, relying on the judgement of clinicians. We present a novel method to compute a quantitative, normalized measure of tissue radiodensity traversed by the primary beam incident on each pixel of a mammogram, a measure we term the standard attenuation rate (SAR). SAR enables: the estimation of breast density which is linked to cancer risk; direct comparison between images; the full potential of computer aided diagnosis to be utilized; and a basis for digital breast tomosynthesis reconstruction. It does this by removing the effects of the imaging conditions under which the mammogram is acquired. First, the x-ray spectrum incident upon the breast is calculated, and from this, the energy exiting the breast is calculated. The contribution of scattered radiation is calculated and subtracted. The SAR measure is the scaling factor that must be applied to the reference material in order to match the primary attenuation of the breast. Specifically, this is the scaled reference material attenuation which when traversed by an identical beam to that traversing the breast, and when subsequently detected, results in the primary component of the pixel intensity observed in the breast image. We present results using two tissue equivalent phantoms, as well as a sensitivity analysis to detector response changes over time and possible errors in compressed thickness measurement.

A model of primary and scattered photon fluence for mammographic x-ray image quantification.

We present an efficient method to calculate the primary and scattered x-ray photon fluence component of a mammographic image. This can be used for a range of clinically important purposes, including estimation of breast density, personalized image display, and quantitative mammogram analysis. The method is based on models of: the x-ray tube; the digital detector; and a novel ray tracer which models the diverging beam emanating from the focal spot. The tube model includes consideration of the anode heel effect, and empirical corrections for wear and manufacturing tolerances. The detector model is empirical, being based on a family of transfer functions that cover the range of beam qualities and compressed breast thicknesses which are encountered clinically. The scatter estimation utilizes optimal information sampling and interpolation (to yield a clinical usable computation time) of scatter calculated using fundamental physics relations. A scatter kernel arising around each primary ray is calculated, and these are summed by superposition to form the scatter image. Beam quality, spatial position in the field (in particular that arising at the air-boundary due to the depletion of scatter contribution from the surroundings), and the possible presence of a grid, are considered, as is tissue composition using an iterative refinement procedure. We present numerous validation results that use a purpose designed tissue equivalent step wedge phantom. The average differences between actual acquisitions and modelled pixel intensities observed across the adipose to fibroglandular attenuation range vary between 5% and 7%, depending on beam quality and, for a single beam quality are 2.09% and 3.36% respectively with and without a grid.