Mammographic density. Measurement of mammographic density
© BioMed Central Ltd 2008
Published: 19 June 2008
Skip to main content
© BioMed Central Ltd 2008
Published: 19 June 2008
Mammographic density has been strongly associated with increased risk of breast cancer. Furthermore, density is inversely correlated with the accuracy of mammography and, therefore, a measurement of density conveys information about the difficulty of detecting cancer in a mammogram. Initial methods for assessing mammographic density were entirely subjective and qualitative; however, in the past few years methods have been developed to provide more objective and quantitative density measurements. Research is now underway to create and validate techniques for volumetric measurement of density. It is also possible to measure breast density with other imaging modalities, such as ultrasound and MRI, which do not require the use of ionizing radiation and may, therefore, be more suitable for use in young women or where it is desirable to perform measurements more frequently. In this article, the techniques for measurement of density are reviewed and some consideration is given to their strengths and limitations.
In 1976, John Wolfe, a radiologist who specialized in mammography, first proposed that there was a strong association between the "parenchymal patterns" seen in the mammogram and the risk that a women would later develop breast cancer [2, 3]. He defined four patterns (later known as Wolfe grades) to characterize the breast. The N pattern, which represented a fatty radiolucent breast, connoted the lowest breast cancer risk. The P1 and P2 patterns indicated progressively greater levels of prominence of fibrous tissue surrounding the ducts and correspondingly higher risk, while the DY pattern indicated the highest risk with a breast that contained dense sheets of fibroglandular tissue. The association of the Wolfe patterns with risk of breast cancer has been reviewed by Saftlas and Szklo  and by Goodwin and Boyd , who concluded that there is a two- to three-fold increase in risk between the N and DY patterns. Because it appears that it is the increasing prevalence of fibroglandular tissue in the breast that gives rise to the increased risk, most subsequent work in this field has attempted to measure mammographic density explicitly.
Each of the images in Figure 1 was selected as representative of one of the categories of a six-category classification (SCC) scheme, which is quantitative for the proportion of the breast appearing as mammographically dense tissue. The six categories range from an absence of density to extensive density (the exact categories for the classification of Figure 1 are summarized in the legend).
Currently, a widely used density classification scheme is the one associated with the Breast Imaging Reporting and Data System (BIRADS)  for reporting findings on mammography. This density system has four categories: BIRADS-1 indicates a predominantly fatty breast; BIRADS-2 scattered fibroglandular densities; BIRADS-3 a breast that is heterogeneously dense; and BIRADS-4, the highest level, an extremely dense breast that could obscure a lesion. This qualitative system was not developed to quantify risk, but to allow an interpreting radiologist to indicate the level of concern that a cancer in the breast might be missed on mammography due to masking by dense tissue. It is well known that the sensitivity of mammography is decreased in the dense breast [7, 8] and a high BIRADS score tells a referring physician who is concerned about breast cancer that other tests less affected by density, such as ultrasound or magnetic resonance imaging (MRI), might be warranted. More recently, in an attempt to make the BIRADS density system more quantitative, it has been recommended that mammograms be classified into four density categories with upper bounds of 24%, 49%, 74% and 100%.
Planimetry refers to the direct measurement of the area of dense tissue seen on the mammogram. Typically, it is performed by tracing around the regions of dense tissue on the mammogram using an instrument called a planimeter. This integrates the total enclosed area. A similar measurement of the total projected area of the breast on the mammogram is also made and the first measurement is divided by the second to obtain the fractional area of the breast that is considered to be dense. This measurement is straightforward to perform, but becomes increasingly labor intensive if it is attempted to separately measure the individual 'islands' of dense tissue that are frequently present in the image. This method was used in the work by Wolfe and colleagues [9, 10].
For many of the quantitative density measurement techniques (but not planimetry) the image must first be digitized. Generally, this is accomplished using a device that scans the film point by point (raster scanning) or line by line with an intense, highly collimated light source such as a laser. In the digitized image, the brightness of each picture element (pixel) is represented by a numerical value ranging from 0 to 2n - 1, where n is the number of bits of digitization. This value can be either linearly or logarithmically related to the brightness. It is important that the digitizer be capable of registering signals over the full range of film opacity from clear to fully black without saturating. Generally, this requires a digitizer with at least 12 bits of precision (4,096 grey levels) combined with an optical design that supports this range. For example, the digitizer must be free from sources of extraneous glare light that would interfere with the measurement.
The digitizer must have adequate spatial resolution to allow local density changes to be tracked. For most measurements the resolution requirement is moderate and even fairly coarse digitization (for example, 0.25 mm pixels) is adequate. For more advanced measures such as texture analysis [11–13] or if the digitized image is also to be used for diagnostic purposes, it may be necessary to digitize to pixels as small as 0.050 mm (50 μm).
As a less time-consuming alternative to planimetry for providing a quantitative estimate of mammographically dense tissue, a simple observer-assisted technique called interactive thresholding was developed by our group . This technique can easily be applied to a digital representation of the mammogram.
Several investigators have developed methods for analyzing mammographic patterns according to texture and found that these texture measures were associated, to varying degrees, with risk [11, 26–31]. For example, Caldwell and colleagues  tested the correlation of the fractal dimension of the digitized mammogram with the Wolfe parenchymal patterns. Magnin and colleagues  in France and Giger's group at The University of Chicago [29, 30] have evaluated the ability of a number of computer-calculated image texture measures to predict risk. While these ideas are intriguing and may lead to more powerful analytical tools in the future, none has yet been demonstrated to provide as strong an association with breast cancer risk as have more direct measures of mammographic density.
While a strong association has been demonstrated between percent mammographic density by area and breast cancer risk, it is more logical that risk is related more directly to the number of target cells, which in turn will be proportional to the 'volume' or fractional volume of dense tissue in the breast.
The most straightforward method of measuring volumetric radiological density is from computed tomography (CT). The CT scan is actually a three-dimensional reconstruction of the X-ray attenuation coefficient of tissues presented as a series of planar images. The values of each image pixel characterize the tissue in terms of its effective atomic number and electron density in a more or less continuous manner. Alternatively, if desired, a simple binary threshold can be reliably set to differentiate between fat-like and water-like tissues and the volume of each type of tissue, the total breast volume and the fraction by volume of each tissue type can be calculated.
Such data for the breast can be obtained from thoracic CT imaging performed without contrast media. There are also dedicated breast CT systems now under development in which only the breast is irradiated [32, 33]. These can provide X-ray attenuation data corresponding to volume elements within the breast. One possible limitation of the latter is that because the breast is imaged when it is pendant into the imaging system with the woman lying prone on a table, some tissue near the chest wall may be excluded from the image and, therefore, from the calculation.
where l and h represent the low and high energies. From these two equations in two unknowns, the effective values of t bone and t soft can be estimated. The same approach has been used for breast density measurement by Shepherd and colleagues , who have built a dedicated system for this purpose. Instead of bone and soft tissue, transmission through the breast is analyzed in terms of effective thicknesses of fibroglandular tissue and fat. Such a system should offer very precise results and, because of the narrow X-ray beams used, should be free of the effects of scattered radiation. One negative aspect is that even though the required radiation dose is very low, it does require a separate procedure be undertaken by the woman while most other methods simply make use of a mammogram that was obtained for other purposes.
Until three-dimensional X-ray breast imaging techniques become widely used, it is most practical to obtain volumetric density information from images produced by two-dimensional mammography systems. Several authors have suggested methods for doing this [37–39]. All methods are based on the known exponential attenuation properties of X-rays. If I 0 X-rays of energy E are incident upon a breast of thickness T, with effective X-ray attenuation coefficient μ(E), then the number transmitted that can be measured by an imaging system is:I tr (E) = I 0(E)e -μT
This relationship is based on two simplifying assumptions: first, that the X-rays are monoenergetic; and second, that no X-rays scattered in the breast reach the imaging system.
Proceeding further, we can consider the breast to be composed of only two materials, fibroglandular tissue and fat, of thicknesses t fib and t fat such that for any path of the X-ray beam though the breast of length T:T = t fib + t fat
where is the fractional density along the measured path.
If I 0/I tr is measured and T is known, then m can be calculated using the known attenuation coefficients of fibroglandular tissue and fat.
In practice, X-ray beams available for clinical mammography are polyenergetic, so the problem becomes more complicated. Researchers have taken different approaches to solve this problem. For example, Highnam and colleagues [39, 40] have chosen to create a physics model of the complete image forming system, including the X-ray source, X-ray scattering and scatter removal and the image receptor, and have calculated a quantity referred to as hint, the thickness of "interesting" (that is, fibroglandular) tissue. Van Engeland and colleagues  developed a physical model to describe image acquisition of full-field digital mammograms and demonstrated good correlation of their volumetric density measurements with three-dimensional data from breast MRI.
In digital mammography the screen-film image receptor is replaced by a detector that produces an electronic signal that precisely and predictably (with generally linear or logarithmic response) tracks the fluence of X-rays transmitted by the breast over a very wide range. This signal is digitized and the image is stored as a matrix in computer memory. This greatly facilitates quantitative density measurement, both because of the improved quality of the signal and because it is no longer necessary to scan the mammogram to digitize it. Furthermore, the modern X-ray systems used for digital mammography give highly reproducible X-ray outputs, largely eliminating the need to monitor drifting of signals from image to image.
There are a few important considerations associated with density assessment from digital mammograms. Many digital mammography systems produce images in two forms, commonly referred to as the 'for processing' or 'raw' image and the 'for presentation' or 'processed' images. The raw image data are based on the detector signal, which is normally proportional to X-ray transmission through the breast and, therefore, should relate closely to breast composition. This image would be subjected to only slight corrections, for example, to compensate for detector flaws.
In order to make the information more suitable for display on a computer screen or for laser printing on film, these images subsequently undergo extensive processing. The image processing operations may be linear or non-linear and may be applied globally (that is, over the entire image in a consistent manner) or locally. These algorithms are largely proprietary to the manufacturers of the digital mammography systems so that the exact details on what they do to the image data are not known. Such processing is likely to distort the relationship between the image signal and X-ray transmission and, thereby, interfere with the ability to derive density information from the images. Radiologists report that when viewing these processed mammograms, breasts appear to be less dense than when imaged with film mammography. For the purpose of cancer detection this is generally considered to be advantageous. In particular, attempting to measure density using thresholding algorithms or physics-based modeling algorithms is likely to be problematic, especially if such measurements are to be compared to those obtained from film mammograms.
It is strongly recommended that density analysis from digital mammograms be performed using the raw image data. Certainly volumetric analysis should be more accurate when done in this way. If it is desired to perform two-dimensional thresholding on digital mammograms, the best approach may be to transform the raw image using a clearly defined global processing algorithm that emulates the characteristics of mammography film before utilizing the thresholding algorithm to measure density. Work to evaluate the performance of density measurements made in this way is currently in progress in my laboratory.
There has been relatively little work done in comparing the measurements provided by different breast density measurement techniques and even less on comparing their performance in predicting breast cancer risk. In a study utilizing mammograms from 65 women, Martin and colleagues  compared two-dimensional density measurements derived from several qualitative, quantitative and semi-automated methods. These included a ten-category subjective scale based on percent density, the qualitative BIRADS scale, a newly introduced quantitative BIRADS scale (four quartiles), and a semi-automated version of a system similar to that described by Byng and colleagues . Consistent with the observations of Warner and colleagues , they found large differences between assessments based on qualitative and quantitative methods. Qualitative assessments were also less reproducible. The authors also observed that qualitative assessments tended to overestimate the degree of density.
There have also been a limited number of comparisons between volumetric and area-based methods. While there is reason to presume that the latter should better correlate with the biological factors responsible for breast cancer risk, volumetric methods depend critically on knowledge of breast thickness, which is difficult to determine accurately in the clinical environment. This may be responsible for the recent findings that the volumetric technique developed by Highnam was less reliable than threshold-based two-dimensional thresholding  and did not provide a stronger predictor of breast cancer risk .
Although most of the work on breast density measurement has been done with mammography, other medical breast imaging modalities also provide information about tissue composition. These have the advantages of providing three-dimensional images and do not involve exposure of the breast to ionizing radiation. One of these is ultrasound. Although the images primarily are sensitive to acoustic reflections at tissue boundaries, the signals are also dependent on the speed of sound and its attenuation, and all three of these factors are, in turn, dependent on tissue composition. There is indication that measurements with ultrasound could provide equivalent density information to that from mammography [46–48]. One of the current limitations of ultrasound, however, is that imaging is highly operator dependent, and this will likely lead to variability in density measurement. Nevertheless, it should be possible to produce an automated volume ultrasound system that would be reproducible and produce reliable quantitative results.
MRI images can be produced that provide signals related to the fat and water composition of the breast. Since the water composition is highly correlated with the prevalence of fibroglandular tissue, these images should be useful for density assessment. Several groups are developing approaches to quantifying density using MRI [49, 50].
Several methods are available for the measurement of breast density. Generally, quantitative approaches that use data extracted from the digitized mammogram allow more precise and reliable measurement than possible with subjective and qualitative techniques. Methods for volumetric assessment of density are currently being developed and evaluated as well as techniques that do not require the use of ionizing radiation.
Breast Imaging Reporting and Data System
magnetic resonance imaging
I would like to express my appreciation to Professor Norman Boyd, with whom I have enjoyed a collaboration over the past 18 years in the investigation of mammographic density and its connection with risk and the accuracy of breast cancer detection. Some of the work described here has been supported by funding from a Terry Fox Program Project Award through the National Cancer Institute of Canada.