A calibration approach to glandular tissue composition estimation in digital mammography

Studies of parenchymal texture added to mammographic breast density and risk of breast cancer: a systematic review of the methods used in the literature

This systematic review aimed to assess the methods used to classify mammographic breast parenchymal features in relation to the prediction of future breast cancer. The databases including Medline (Ovid) 1946-, Embase.com 1947-, CINAHL Plus 1937-, Scopus 1823-, Cochrane Library (including CENTRAL), and Clinicaltrials.gov were searched through October 2021 to extract published articles in English describing the relationship of parenchymal texture features with the risk of breast cancer. Twenty-eight articles published since 2016 were included in the final review. The identification of parenchymal texture features varied from using a predefined list to machine-driven identification. A reduction in the number of features chosen for subsequent analysis in relation to cancer incidence then varied across statistical approaches and machine learning methods. The variation in approach and number of features identified for inclusion in analysis precluded generating a quantitative summary or meta-analysis of the value of these features to improve predicting risk of future breast cancers. This updated overview of the state of the art revealed research gaps; based on these, we provide recommendations for future studies using parenchymal features for mammogram images to make use of accumulating image data, and external validation of prediction models that extend to 5 and 10 years to guide clinical risk management. Following these recommendations could enhance the applicability of models, helping improve risk classification and risk prediction for women to tailor screening and prevention strategies to the level of risk.

Quantification of water and lipid density with dual-energy mammography: validation in postmortem breasts

OBJECTIVES

Breast cancer is the most common cancer in women and the second leading cause of cancer death. It is well known that breast density is an important risk factor for breast cancer and also can be used to personalize screening and for assessment of treatment response. Breast density has previously been correlated to volumetric water density. The purpose of this study is to validate the accuracy and precision of dual-energy mammography in measuring water density in postmortem breasts.

METHODS

Twenty pairs of postmortem breasts were imaged using dual-energy mammography with energy-sensitive photon-counting detectors. Chemical analysis was used as the reference standard to assess the accuracy of dual-energy mammography in measuring volumetric water and lipid density. Images from different views and contralateral breasts were used to assess estimate of precision for water and lipid volumetric density measurements.

RESULTS

The measured volumetric water and lipid density from dual-energy mammography and chemical analysis were in good agreement, where the standard errors of estimates (SEE) of both were calculated to be 2.1%. Volumetric water and lipid density measurements from different views were also in good agreement, with a SEE of 1.3% and 1.1%, respectively.

CONCLUSIONS

The results indicate that dual-energy mammography can be used to accurately measure volumetric water and lipid density in breast tissue. Accurate quantification of volumetric water density is expected to enhance its utility as a risk factor for breast cancer and for assessment of response to therapy.

KEY POINTS

• Dual-energy mammography can be used to accurately measure water and lipid volumetric density in breast tissue. • Improved quantification of volumetric water density is expected to enhance its utility for assessment of response to therapy and as a risk factor for breast cancer.

Calibrated Breast Density Measurements

RATIONALE AND OBJECTIVES

Mammographic density is an important risk factor for breast cancer, but translation to the clinic requires assurance that prior work based on mammography is applicable to current technologies. The purpose of this work is to evaluate whether a calibration methodology developed previously produces breast density metrics predictive of breast cancer risk when applied to a case-control study.

MATERIALS AND METHODS

A matched case control study (n = 319 pairs) was used to evaluate two calibrated measures of breast density. Two-dimensional mammograms were acquired from six Hologic mammography units: three conventional Selenia two-dimensional full-field digital mammography systems and three Dimensions digital breast tomosynthesis systems. We evaluated the capability of two calibrated breast density measures to quantify breast cancer risk: the mean (PGm) and standard deviation (PGsd) of the calibrated pixels. Matching variables included age, hormone replacement therapy usage/duration, screening history, and mammography unit. Calibrated measures were compared to the percentage of breast density (PD) determined with the operator-assisted Cumulus method. Conditional logistic regression was used to generate odds ratios (ORs) from continuous and quartile (Q) models with 95% confidence intervals. The area under the receiver operating characteristic curve (Az) was also used as a comparison metric. Both univariate models and models adjusted for body mass index and ethnicity were evaluated.

RESULTS

In adjusted models, both PGsd and PD were statistically significantly associated with breast cancer with similar Az of 0.61-0.62. The corresponding ORs and confidence intervals were also similar. For PGsd, the OR was 1.34 (1.09, 1.66) for the continuous measure and 1.83 (1.11, 3.02), 2.19 (1.28, 3.73), and 2.20 (1.26, 3.85) for Q2-Q4. For PD, the OR was 1.43 (1.16, 1.76) for the continuous measure and 0.84 (0.52, 1.38), 1.96 (1.19, 3.23), and 2.27 (1.29, 4.00) for Q2-Q4. The results for PGm were slightly attenuated and not statistically significant. The OR was 1.22 (0.99, 1.51) with Az = 0.60 for the continuous measure and 1.24 (0.78, 1.97), 0.98 (0.60, 1.61), and 1.26, (0.77, 2.07) for Q2-Q4 with Az = 0.60.

CONCLUSION

The calibrated PGsd measure provided significant associations with breast cancer comparable to those given by PD. The calibrated PGm performed slightly worse. These findings indicate that the calibration approach developed previously replicates under more general conditions.

Radiomic Phenotypes of Mammographic Parenchymal Complexity: Toward Augmenting Breast Density in Breast Cancer Risk Assessment

Purpose To identify phenotypes of mammographic parenchymal complexity by using radiomic features and to evaluate their associations with breast density and other breast cancer risk factors. Materials and Methods Computerized image analysis was used to quantify breast density and extract parenchymal texture features in a cross-sectional sample of women screened with digital mammography from September 1, 2012, to February 28, 2013 (n = 2029; age range, 35-75 years; mean age, 55.9 years). Unsupervised clustering was applied to identify and reproduce phenotypes of parenchymal complexity in separate training (n = 1339) and test sets (n = 690). Differences across phenotypes by age, body mass index, breast density, and estimated breast cancer risk were assessed by using Fisher exact, χ², and Kruskal-Wallis tests. Conditional logistic regression was used to evaluate preliminary associations between the detected phenotypes and breast cancer in an independent case-control sample (76 women diagnosed with breast cancer and 158 control participants) matched on age. Results Unsupervised clustering in the screening sample identified four phenotypes with increasing parenchymal complexity that were reproducible between training and test sets (P = .001). Breast density was not strongly correlated with phenotype category (R² = 0.24 for linear trend). The low- to intermediate-complexity phenotype (prevalence, 390 of 2029 [19%]) had the lowest proportion of dense breasts (eight of 390 [2.1%]), whereas similar proportions were observed across other phenotypes (from 140 of 291 [48.1%] in the high-complexity phenotype to 275 of 511 [53.8%] in the low-complexity phenotype). In the independent case-control sample, phenotypes showed a significant association with breast cancer (P = .001), resulting in higher discriminatory capacity when added to a model with breast density and body mass index (area under the curve, 0.84 vs 0.80; P = .03 for comparison). Conclusion Radiomic phenotypes capture mammographic parenchymal complexity beyond conventional breast density measures and established breast cancer risk factors. © RSNA, 2018 Online supplemental material is available for this article. See also the editorial by Pinker in this issue.

Breast density estimation from high spectral and spatial resolution MRI

A three-dimensional breast density estimation method is presented for high spectral and spatial resolution (HiSS) MR imaging. Twenty-two patients were recruited (under an Institutional Review Board--approved Health Insurance Portability and Accountability Act-compliant protocol) for high-risk breast cancer screening. Each patient received standard-of-care clinical digital x-ray mammograms and MR scans, as well as HiSS scans. The algorithm for breast density estimation includes breast mask generating, breast skin removal, and breast percentage density calculation. The inter- and intra-user variabilities of the HiSS-based density estimation were determined using correlation analysis and limits of agreement. Correlation analysis was also performed between the HiSS-based density estimation and radiologists' breast imaging-reporting and data system (BI-RADS) density ratings. A correlation coefficient of 0.91 ([Formula: see text]) was obtained between left and right breast density estimations. An interclass correlation coefficient of 0.99 ([Formula: see text]) indicated high reliability for the inter-user variability of the HiSS-based breast density estimations. A moderate correlation coefficient of 0.55 ([Formula: see text]) was observed between HiSS-based breast density estimations and radiologists' BI-RADS. In summary, an objective density estimation method using HiSS spectral data from breast MRI was developed. The high reproducibility with low inter- and low intra-user variabilities shown in this preliminary study suggest that such a HiSS-based density metric may be potentially beneficial in programs requiring breast density such as in breast cancer risk assessment and monitoring effects of therapy.

Improvement in risk prediction, early detection and prevention of breast cancer in the NHS Breast Screening Programme and family history clinics: a dual cohort study

Background

In the UK, women are invited for 3-yearly mammography screening, through the NHS Breast Screening Programme (NHSBSP), from the ages of 47–50 years to the ages of 69–73 years. Women with family histories of breast cancer can, from the age of 40 years, obtain enhanced surveillance and, in exceptionally high-risk cases, magnetic resonance imaging. However, no NHSBSP risk assessment is undertaken. Risk prediction models are able to categorise women by risk using known risk factors, although accurate individual risk prediction remains elusive. The identification of mammographic breast density (MD) and common genetic risk variants [single nucleotide polymorphisms (SNPs)] has presaged the improved precision of risk models.

Objectives

To (1) identify the best performing model to assess breast cancer risk in family history clinic (FHC) and population settings; (2) use information from MD/SNPs to improve risk prediction; (3) assess the acceptability and feasibility of offering risk assessment in the NHSBSP; and (4) identify the incremental costs and benefits of risk stratified screening in a preliminary cost-effectiveness analysis.

Design

Two cohort studies assessing breast cancer incidence.

Setting

High-risk FHC and the NHSBSP Greater Manchester, UK.

Participants

A total of 10,000 women aged 20–79 years [Family History Risk Study (FH-Risk); UK Clinical Research Network identification number (UKCRN-ID) 8611] and 53,000 women from the NHSBSP [aged 46–73 years; Predicting the Risk of Cancer At Screening (PROCAS) study; UKCRN-ID 8080].

Interventions

Questionnaires collected standard risk information, and mammograms were assessed for breast density by a number of techniques. All FH-Risk and 10,000 PROCAS participants participated in deoxyribonucleic acid (DNA) studies. The risk prediction models Manual method, Tyrer–Cuzick (TC), BOADICEA (Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm) and Gail were used to assess risk, with modelling based on MD and SNPs. A preliminary model-based cost-effectiveness analysis of risk stratified screening was conducted.

Main outcome measures

Breast cancer incidence.

Data sources

The NHSBSP; cancer registration.

Results

A total of 446 women developed incident breast cancers in FH-Risk in 97,958 years of follow-up. All risk models accurately stratified women into risk categories. TC had better risk precision than Gail, and BOADICEA accurately predicted risk in the 6268 single probands. The Manual model was also accurate in the whole cohort. In PROCAS, TC had better risk precision than Gail [area under the curve (AUC) 0.58 vs. 0.54], identifying 547 prospective breast cancers. The addition of SNPs in the FH-Risk case–control study improved risk precision but was not useful inBRCA1(breast cancer 1 gene) families. Risk modelling of SNPs in PROCAS showed an incremental improvement from using SNP18 used in PROCAS to SNP67. MD measured by visual assessment score provided better risk stratification than automatic measures, despite wide intra- and inter-reader variability. Using a MD-adjusted TC model in PROCAS improved risk stratification (AUC = 0.6) and identified significantly higher rates (4.7 per 10,000 vs. 1.3 per 10,000;p < 0.001) of high-stage cancers in women with above-average breast cancer risks. It is not possible to provide estimates of the incremental costs and benefits of risk stratified screening because of lack of data inputs for key parameters in the model-based cost-effectiveness analysis.

Conclusions

Risk precision can be improved by using DNA and MD, and can potentially be used to stratify NHSBSP screening. It may also identify those at greater risk of high-stage cancers for enhanced screening. The cost-effectiveness of risk stratified screening is currently associated with extensive uncertainty. Additional research is needed to identify data needed for key inputs into model-based cost-effectiveness analyses to identify the impact on health-care resource use and patient benefits.

Future work

A pilot of real-time NHSBSP risk prediction to identify women for chemoprevention and enhanced screening is required.

Funding

The National Institute for Health Research Programme Grants for Applied Research programme. The DNA saliva collection for SNP analysis for PROCAS was funded by the Genesis Breast Cancer Prevention Appeal.

Volumetric breast density measurement: sensitivity analysis of a relative physics approach

OBJECTIVE

To investigate the sensitivity and robustness of a volumetric breast density (VBD) measurement system to errors in the imaging physics parameters including compressed breast thickness (CBT), tube voltage (kVp), filter thickness, tube current-exposure time product (mAs), detector gain, detector offset and image noise.

METHODS

3317 raw digital mammograms were processed with Volpara(®) (Matakina Technology Ltd, Wellington, New Zealand) to obtain fibroglandular tissue volume (FGV), breast volume (BV) and VBD. Errors in parameters including CBT, kVp, filter thickness and mAs were simulated by varying them in the Digital Imaging and Communications in Medicine (DICOM) tags of the images up to ±10% of the original values. Errors in detector gain and offset were simulated by varying them in the Volpara configuration file up to ±10% from their default values. For image noise, Gaussian noise was generated and introduced into the original images.

RESULTS

Errors in filter thickness, mAs, detector gain and offset had limited effects on FGV, BV and VBD. Significant effects in VBD were observed when CBT, kVp, detector offset and image noise were varied (p < 0.0001). Maximum shifts in the mean (1.2%) and median (1.1%) VBD of the study population occurred when CBT was varied.

CONCLUSION

Volpara was robust to expected clinical variations, with errors in most investigated parameters giving limited changes in results, although extreme variations in CBT and kVp could lead to greater errors.

ADVANCES IN KNOWLEDGE

Despite Volpara's robustness, rigorous quality control is essential to keep the parameter errors within reasonable bounds. Volpara appears robust within those bounds, albeit for more advanced applications such as tracking density change over time, it remains to be seen how accurate the measures need to be.

Consistency of breast density categories in serial screening mammograms: A comparison between automated and human assessment

Reliable breast density measurement is needed to personalize screening by using density as a risk factor and offering supplemental screening to women with dense breasts. We investigated the categorization of pairs of subsequent screening mammograms into density classes by human readers and by an automated system. With software (VDG) and by four readers, including three specialized breast radiologists, 1000 mammograms belonging to 500 pairs of subsequent screening exams were categorized into either two or four density classes. We calculated percent agreement and the percentage of women that changed from dense to non-dense and vice versa. Inter-exam agreement (IEA) was calculated with kappa statistics. Results were computed for each reader individually and for the case that each mammogram was classified by one of the four readers by random assignment (group reading). Higher percent agreement was found with VDG (90.4%, CI 87.9-92.9%) than with readers (86.2-89.2%), while less plausible changes from non-dense to dense occur less often with VDG (2.8%, CI 1.4-4.2%) than with group reading (4.2%, CI 2.4-6.0%). We found an IEA of 0.68-0.77 for the readers using two classes and an IEA of 0.76-0.82 using four classes. IEA is significantly higher with VDG compared to group reading. The categorization of serial mammograms in density classes is more consistent with automated software than with a mixed group of human readers. When using breast density to personalize screening protocols, assessment with software may be preferred over assessment by radiologists.

Postmortem validation of breast density using dual-energy mammography

PURPOSE

Mammographic density has been shown to be an indicator of breast cancer risk and also reduces the sensitivity of screening mammography. Currently, there is no accepted standard for measuring breast density. Dual energy mammography has been proposed as a technique for accurate measurement of breast density. The purpose of this study is to validate its accuracy in postmortem breasts and compare it with other existing techniques.

METHODS

Forty postmortem breasts were imaged using a dual energy mammography system. Glandular and adipose equivalent phantoms of uniform thickness were used to calibrate a dual energy basis decomposition algorithm. Dual energy decomposition was applied after scatter correction to calculate breast density. Breast density was also estimated using radiologist reader assessment, standard histogram thresholding and a fuzzy C-mean algorithm. Chemical analysis was used as the reference standard to assess the accuracy of different techniques to measure breast composition.

RESULTS

Breast density measurements using radiologist reader assessment, standard histogram thresholding, fuzzy C-mean algorithm, and dual energy were in good agreement with the measured fibroglandular volume fraction using chemical analysis. The standard error estimates using radiologist reader assessment, standard histogram thresholding, fuzzy C-mean, and dual energy were 9.9%, 8.6%, 7.2%, and 4.7%, respectively.

CONCLUSIONS

The results indicate that dual energy mammography can be used to accurately measure breast density. The variability in breast density estimation using dual energy mammography was lower than reader assessment rankings, standard histogram thresholding, and fuzzy C-mean algorithm. Improved quantification of breast density is expected to further enhance its utility as a risk factor for breast cancer.

A Review on Automatic Mammographic Density and Parenchymal Segmentation

Breast cancer is the most frequently diagnosed cancer in women. However, the exact cause(s) of breast cancer still remains unknown. Early detection, precise identification of women at risk, and application of appropriate disease prevention measures are by far the most effective way to tackle breast cancer. There are more than 70 common genetic susceptibility factors included in the current non-image-based risk prediction models (e.g., the Gail and the Tyrer-Cuzick models). Image-based risk factors, such as mammographic densities and parenchymal patterns, have been established as biomarkers but have not been fully incorporated in the risk prediction models used for risk stratification in screening and/or measuring responsiveness to preventive approaches. Within computer aided mammography, automatic mammographic tissue segmentation methods have been developed for estimation of breast tissue composition to facilitate mammographic risk assessment. This paper presents a comprehensive review of automatic mammographic tissue segmentation methodologies developed over the past two decades and the evidence for risk assessment/density classification using segmentation. The aim of this review is to analyse how engineering advances have progressed and the impact automatic mammographic tissue segmentation has in a clinical environment, as well as to understand the current research gaps with respect to the incorporation of image-based risk factors in non-image-based risk prediction models.

A population-based tissue probability map-driven level set method for fully automated mammographic density estimations

PURPOSE

A major challenge when distinguishing glandular tissues on mammograms, especially for area-based estimations, lies in determining a boundary on a hazy transition zone from adipose to glandular tissues. This stems from the nature of mammography, which is a projection of superimposed tissues consisting of different structures. In this paper, the authors present a novel segmentation scheme which incorporates the learned prior knowledge of experts into a level set framework for fully automated mammographic density estimations.

METHODS

The authors modeled the learned knowledge as a population-based tissue probability map (PTPM) that was designed to capture the classification of experts' visual systems. The PTPM was constructed using an image database of a selected population consisting of 297 cases. Three mammogram experts extracted regions for dense and fatty tissues on digital mammograms, which was an independent subset used to create a tissue probability map for each ROI based on its local statistics. This tissue class probability was taken as a prior in the Bayesian formulation and was incorporated into a level set framework as an additional term to control the evolution and followed the energy surface designed to reflect experts' knowledge as well as the regional statistics inside and outside of the evolving contour.

RESULTS

A subset of 100 digital mammograms, which was not used in constructing the PTPM, was used to validate the performance. The energy was minimized when the initial contour reached the boundary of the dense and fatty tissues, as defined by experts. The correlation coefficient between mammographic density measurements made by experts and measurements by the proposed method was 0.93, while that with the conventional level set was 0.47.

CONCLUSIONS

The proposed method showed a marked improvement over the conventional level set method in terms of accuracy and reliability. This result suggests that the proposed method successfully incorporated the learned knowledge of the experts' visual systems and has potential to be used as an automated and quantitative tool for estimations of mammographic breast density levels.

Calibrated breast density methods for full field digital mammography: a system for serial quality control and inter-system generalization

PURPOSE

The authors are developing a system for calibrated breast density measurements using full field digital mammography (FFDM). Breast tissue equivalent (BTE) phantom images are used to establish baseline (BL) calibration curves at time zero. For a given FFDM unit, the full BL dataset is comprised of approximately 160 phantom images, acquired prior to calibrating prospective patient mammograms. BL curves are monitored serially to ensure they produce accurate calibration and require updating when calibration accuracy degrades beyond an acceptable tolerance, rather than acquiring full BL datasets repeatedly. BL updating is a special case of generalizing calibration datasets across FFDM units, referred to as cross-calibration. Serial monitoring, BL updating, and cross-calibration techniques were developed and evaluated.

METHODS

BL curves were established for three Hologic Selenia FFDM units at time zero. In addition, one set of serial phantom images, comprised of equal proportions of adipose and fibroglandular BTE materials (50/50 compositions) of a fixed height, was acquired biweekly and monitored with the cumulative sum (Cusum) technique. These 50/50 composition images were used to update the BL curves when the calibration accuracy degraded beyond a preset tolerance of ±4 standardized units. A second set of serial images, comprised of a wide-range of BTE compositions, was acquired biweekly to evaluate serial monitoring, BL updating, and cross-calibration techniques.

RESULTS

Calibration accuracy can degrade serially and is a function of acquisition technique and phantom height. The authors demonstrated that all heights could be monitored simultaneously while acquiring images of a 50/50 phantom with a fixed height for each acquisition technique biweekly, translating into approximately 16 image acquisitions biweekly per FFDM unit. The same serial images are sufficient for serial monitoring, BL updating, and cross-calibration. Serial calibration accuracy was maintained within ±4 standardized unit variation from the ideal when applying BL updating. BL updating is a special case of cross-calibration; the BL dataset of unit 1 can be converted to the BL dataset for another similar unit (i.e., unit 2) at any given time point using the 16 serial monitoring 50/50 phantom images of unit 2 (or vice versa) acquired near this time point while maintaining the ±4 standardized unit tolerance.

CONCLUSIONS

A methodology for monitoring and maintaining serial calibration accuracy for breast density measurements was evaluated. Calibration datasets for a given unit can be translated forward in time with minimal phantom imaging effort. Similarly, cross-calibration is a method for generalizing calibration datasets across similar units without additional phantom imaging. This methodology will require further evaluation with mammograms for complete validation.

The simulation of 3D mass models in 2D digital mammography and breast tomosynthesis

PURPOSE

This work proposes a new method of building 3D breast mass models with different morphological shapes and describes the validation of the realism of their appearance after simulation into 2D digital mammograms and breast tomosynthesis images.

METHODS

Twenty-five contrast enhanced MRI breast lesions were collected and each mass was manually segmented in the three orthogonal views: sagittal, coronal, and transversal. The segmented models were combined, resampled to have isotropic voxel sizes, triangularly meshed, and scaled to different sizes. These masses were referred to as nonspiculated masses and were then used as nuclei onto which spicules were grown with an iterative branching algorithm forming a total of 30 spiculated masses. These 55 mass models were projected into 2D projection images to obtain mammograms after image processing and into tomographic sequences of projection images, which were then reconstructed to form 3D tomosynthesis datasets. The realism of the appearance of these mass models was assessed by five radiologists via receiver operating characteristic (ROC) analysis when compared to 54 real masses. All lesions were also given a breast imaging reporting and data system (BIRADS) score. The data sets of 2D mammography and tomosynthesis were read separately. The Kendall's coefficient of concordance was used for the interrater observer agreement assessment for the BIRADS scores per modality. Further paired analysis, using the Wilcoxon signed rank test, of the BIRADS assessment between 2D and tomosynthesis was separately performed for the real masses and for the simulated masses.

RESULTS

The area under the ROC curves, averaged over all observers, was 0.54 (95% confidence interval [0.50, 0.66]) for the 2D study, and 0.67 (95% confidence interval [0.55, 0.79]) for the tomosynthesis study. According to the BIRADS scores, the nonspiculated and the spiculated masses varied in their degrees of malignancy from normal (BIRADS 1) to highly suggestive for malignancy (BIRADS 5) indicating the required variety of shapes and margins of these models. The assessment of the BIRADS scores for all observers indicated good agreement based on Kendall's coefficient for both the 2D and the tomosynthesis evaluations. The paired analysis of the BIRADS scores between 2D and tomosynthesis for each observer revealed consistent behavior for the real and simulated masses.

CONCLUSIONS

A database of 3D mass models, with variety of shapes and margins, was validated for the realism of their appearance for 2D digital mammography and for breast tomosynthesis. This database is suitable for use in future observer performance studies whether in virtual clinical trials or in patient images with simulated lesions.

Comparison of volumetric breast density estimations from mammography and thorax CT

Breast density has become an important issue in current breast cancer screening, both as a recognized risk factor for breast cancer and by decreasing screening efficiency by the masking effect. Different qualitative and quantitative methods have been proposed to evaluate area-based breast density and volumetric breast density (VBD). We propose a validation method comparing the computation of VBD obtained from digital mammographic images (VBDMX) with the computation of VBD from thorax CT images (VBDCT). We computed VBDMX by applying a conversion function to the pixel values in the mammographic images, based on models determined from images of breast equivalent material. VBDCT is computed from the average Hounsfield Unit (HU) over the manually delineated breast volume in the CT images. This average HU is then compared to the HU of adipose and fibroglandular tissues from patient images. The VBDMX method was applied to 663 mammographic patient images taken on two Siemens Inspiration (hospL) and one GE Senographe Essential (hospJ). For the comparison study, we collected images from patients who had a thorax CT and a mammography screening exam within the same year. In total, thorax CT images corresponding to 40 breasts (hospL) and 47 breasts (hospJ) were retrieved. Averaged over the 663 mammographic images the median VBDMX was 14.7% . The density distribution and the inverse correlation between VBDMX and breast thickness were found as expected. The average difference between VBDMX and VBDCT is smaller for hospJ (4%) than for hospL (10%). This study shows the possibility to compare VBDMX with the VBD from thorax CT exams, without additional examinations. In spite of the limitations caused by poorly defined breast limits, the calibration of mammographic images to local VBD provides opportunities for further quantitative evaluations.

Research in digital mammography and tomosynthesis at the University of Toronto

There have been major advances in the field of breast cancer imaging since the early 1970s, both in technological improvements and in the use of the methods of medical physics and image analysis to optimize image quality. The introduction of digital mammography in 2000 provided a marked improvement in imaging of dense breasts. In addition, it became possible to produce tomographic and functional images on modified digital mammography systems. Digital imaging also greatly facilitated the extraction of quantitative information from images. My laboratory has been fortunate in being able to participate in some of these exciting developments. I will highlight some of the areas of our research interest which include modeling of the image formation process, development of high-resolution X-ray detectors for digital mammography and investigating new methods for analyzing image quality. I will also describe our more recent work on developing new applications of digital mammography including tomosynthesis, contrast-enhanced mammography, and measurement of breast density. Finally, I will point to a new area for our research--the application of the techniques of medical imaging to making pathology more quantitative to contribute to use of biomarkers for better characterizing breast cancer and directing therapeutic decisions.

Automated Volumetric Breast Density derived by Shape and Appearance Modeling

The image shape and texture (appearance) estimation designed for facial recognition is a novel and promising approach for application in breast imaging. The purpose of this study was to apply a shape and appearance model to automatically estimate percent breast fibroglandular volume (%FGV) using digital mammograms. We built a shape and appearance model using 2000 full-field digital mammograms from the San Francisco Mammography Registry with known %FGV measured by single energy absorptiometry method. An affine transformation was used to remove rotation, translation and scale. Principal Component Analysis (PCA) was applied to extract significant and uncorrelated components of %FGV. To build an appearance model, we transformed the breast images into the mean texture image by piecewise linear image transformation. Using PCA the image pixels grey-scale values were converted into a reduced set of the shape and texture features. The stepwise regression with forward selection and backward elimination was used to estimate the outcome %FGV with shape and appearance features and other system parameters. The shape and appearance scores were found to correlate moderately to breast %FGV, dense tissue volume and actual breast volume, body mass index (BMI) and age. The highest Pearson correlation coefficient was equal 0.77 for the first shape PCA component and actual breast volume. The stepwise regression method with ten-fold cross-validation to predict %FGV from shape and appearance variables and other system outcome parameters generated a model with a correlation of r² = 0.8. In conclusion, a shape and appearance model demonstrated excellent feasibility to extract variables useful for automatic %FGV estimation. Further exploring and testing of this approach is warranted.

Volumetric breast density estimation from full-field digital mammograms: a validation study

OBJECTIVES

To objectively evaluate automatic volumetric breast density assessment in Full-Field Digital Mammograms (FFDM) using measurements obtained from breast Magnetic Resonance Imaging (MRI).

MATERIAL AND METHODS

A commercially available method for volumetric breast density estimation on FFDM is evaluated by comparing volume estimates obtained from 186 FFDM exams including mediolateral oblique (MLO) and cranial-caudal (CC) views to objective reference standard measurements obtained from MRI.

RESULTS

Volumetric measurements obtained from FFDM show high correlation with MRI data. Pearson's correlation coefficients of 0.93, 0.97 and 0.85 were obtained for volumetric breast density, breast volume and fibroglandular tissue volume, respectively.

CONCLUSIONS

Accurate volumetric breast density assessment is feasible in Full-Field Digital Mammograms and has potential to be used in objective breast cancer risk models and personalized screening.

Tissue Cancellation in Dual Energy Mammography Using a Calibration Phantom Customized for Direct Mapping

An easily implementable tissue cancellation method for dual energy mammography is proposed to reduce anatomical noise and enhance lesion visibility. For dual energy calibration, the images of an imaging object are directly mapped onto the images of a customized calibration phantom. Each pixel pair of the low and high energy images of the imaging object was compared to pixel pairs of the low and high energy images of the calibration phantom. The correspondence was measured by absolute difference between the pixel values of imaged object and those of the calibration phantom. Then the closest pixel pair of the calibration phantom images is marked and selected. After the calibration using direct mapping, the regions with lesion yielded different thickness from the background tissues. Taking advantage of the different thickness, the visibility of cancerous lesions was enhanced with increased contrast-to-noise ratio, depending on the size of lesion and breast thickness. However, some tissues near the edge of imaged object still remained after tissue cancellation. These remaining residuals seem to occur due to the heel effect, scattering, nonparallel X-ray beam geometry and Poisson distribution of photons. To improve its performance further, scattering and the heel effect should be compensated.

A comparison of calibration data from full field digital mammography units for breast density measurements

Background

Breast density is a significant breast cancer risk factor measured from mammograms. The most appropriate method for measuring breast density for risk applications is still under investigation. Calibration standardizes mammograms to account for acquisition technique differences prior to making breast density measurements. We evaluated whether a calibration methodology developed for an indirect x-ray conversion full field digital mammography (FFDM) technology applies to direct x-ray conversion FFDM systems.

Methods

Breast tissue equivalent (BTE) phantom images were used to establish calibration datasets for three similar direct x-ray conversion FFDM systems. The calibration dataset for each unit is a function of the target/filter combination, x-ray tube voltage, current × time (mAs), phantom height, and two detector fields of view (FOVs). Methods were investigated to reduce the amount of calibration data by restricting the height, mAs, and FOV sampling. Calibration accuracy was evaluated with mixture phantoms. We also compared both intra- and inter-system calibration characteristics and accuracy.

Results

Calibration methods developed previously apply to direct x-ray conversion systems with modification. Calibration accuracy was largely within the acceptable range of ± 4 standardized units from the ideal value over the entire acquisition parameter space for the direct conversion units. Acceptable calibration accuracy was maintained with a cubic-spline height interpolation, representing a modification to previous work. Calibration data is unit specific, can be acquired with the large FOV, and requires a minimum of one reference mAs sample. The mAs sampling, calibration accuracy, and the necessity for machine specific calibration data are common characteristics and in agreement with our previous work.

Conclusion

The generality of our calibration approach was established under ideal conditions. Evaluation with patient data using breast cancer status as the endpoint is required to demonstrate that the approach produces a breast density measure associated with breast cancer.

Personalized estimates of radiation dose from dedicated breast CT in a diagnostic population and comparison with diagnostic mammography

This study retrospectively analyzed the mean glandular dose (MGD) to 133 breasts from 132 subjects, all women, who participated in a clinical trial evaluating dedicated breast CT in a diagnostic population. The clinical trial was conducted in adherence to a protocol approved by institutional review boards and the study participants provided written informed consent. Individual estimates of MGD to each breast from dedicated breast CT was obtained by combining x-ray beam characteristics with estimates of breast dimensions and fibroglandular fraction from volumetric breast CT images, and using normalized glandular dose coefficients. For each study participant and for the breast corresponding to that imaged with breast CT, an estimate of the MGD from diagnostic mammography (including supplemental views) was obtained from the DICOM image headers for comparison. This estimate uses normalized glandular dose coefficients corresponding to a breast with 50% fibroglandular weight fraction. The median fibroglandular weight fraction for the study cohort determined from volumetric breast CT images was 15%. Hence, the MGD from diagnostic mammography was corrected to be representative of the study cohort. Individualized estimates of MGD from breast CT ranged from 5.7 to 27.8 mGy. Corresponding to the breasts imaged with breast CT, the MGD from diagnostic mammography ranged from 2.6 to 31.6 mGy. The mean (± inter-breast SD) and the median MGD (mGy) from dedicated breast CT exam were 13.9 ± 4.6 and 12.6, respectively. For the corresponding breasts, the mean (± inter-breast SD) and the median MGD (mGy) from diagnostic mammography were 12.4 ± 6.3 and 11.1, respectively. Statistical analysis indicated that at the 0.05 level, the distributions of MGD from dedicated breast CT and diagnostic mammography were significantly different (Wilcoxon signed ranks test, p = 0.007). While the interquartile range and the range (maximum-minimum) of MGD from dedicated breast CT was lower than diagnostic mammography, the median MGD from dedicated breast CT was approximately 13.5% higher than that from diagnostic mammography. The MGD for breast CT is based on a 1.45 mm skin layer and that for diagnostic mammography is based on a 4 mm skin layer; thus, favoring a lower estimate for MGD from diagnostic mammography. The median MGD from dedicated breast CT corresponds to the median MGD from four to five diagnostic mammography views. In comparison, for the same 133 breasts, the mean and the median number of views per breast during diagnostic mammography were 4.53 and 4, respectively. Paired analysis showed that there was approximately equal likelihood of receiving lower MGD from either breast CT or diagnostic mammography. Future work will investigate methods to reduce and optimize radiation dose from dedicated breast CT.

Wavelet-based noise-model driven denoising algorithm for differential phase contrast mammography

Traditional mammography can be positively complemented by phase contrast and scattering x-ray imaging, because they can detect subtle differences in the electron density of a material and measure the local small-angle scattering power generated by the microscopic density fluctuations in the specimen, respectively. The grating-based x-ray interferometry technique can produce absorption, differential phase contrast (DPC) and scattering signals of the sample, in parallel, and works well with conventional X-ray sources; thus, it constitutes a promising method for more reliable breast cancer screening and diagnosis. Recently, our team proved that this novel technology can provide images superior to conventional mammography. This new technology was used to image whole native breast samples directly after mastectomy. The images acquired show high potential, but the noise level associated to the DPC and scattering signals is significant, so it is necessary to remove it in order to improve image quality and visualization. The noise models of the three signals have been investigated and the noise variance can be computed. In this work, a wavelet-based denoising algorithm using these noise models is proposed. It was evaluated with both simulated and experimental mammography data. The outcomes demonstrated that our method offers a good denoising quality, while simultaneously preserving the edges and important structural features. Therefore, it can help improve diagnosis and implement further post-processing techniques such as fusion of the three signals acquired.

X-ray spectral measurements for tungsten-anode from 20 to 49 kVp on a digital breast tomosynthesis system

PURPOSE

This paper presents new spectral measurements of a tungsten-target digital breast tomosynthesis (DBT) system, including spectra of 43-49 kVp.

METHODS

Raw x-ray spectra of 20-49 kVp were directly measured from the tube port of a Selenia Dimensions DBT system using a CdTe based spectrometer. Two configurations of collimation were employed: one with two tungsten pinholes of 25 μm and 200 μm diameters, and the other with a single pinhole of 25 μm diameter, for acquiring spectra from the focal spot and from the focal spot as well as its vicinity. Stripping correction was applied to the measured spectra to compensate distortions due to escape events. The measured spectra were compared with the existing mammographic spectra of the TASMIP model in terms of photon fluence per exposure, spectral components, and half-value layer (HVL). HVLs were calculated from the spectra with a numerical filtration of 0.7 mm aluminum and were compared against actual measurements on the DBT system using W/Al (target-filter) combination, without paddle in the beam.

RESULTS

The spectra from the double-pinhole configuration, in which the acceptance aperture pointed right at the focal spot, were harder than the single-pinhole spectra which include both primary and off-focus radiation. HVL calculated from the single-pinhole setup agreed with the measured HVL within 0.04 mm aluminum, while the HVL values from the double-pinhole setup were larger than the single-pinhole HVL by at most 0.1 mm aluminum. The spectra from single-pinhole setup agreed well with the TASMIP mammographic spectra, and are more relevant for clinical purpose.

CONCLUSIONS

The spectra data would be useful for future research on DBT system with tungsten targets.

A novel automated mammographic density measure and breast cancer risk

BACKGROUND

Mammographic breast density is a strong breast cancer risk factor but is not used in the clinical setting, partly because of a lack of standardization and automation. We developed an automated and objective measurement of the grayscale value variation within a mammogram, evaluated its association with breast cancer, and compared its performance with that of percent density (PD).

METHODS

Three clinic-based studies were included: a case-cohort study of 217 breast cancer case subjects and 2094 non-case subjects and two case-control studies comprising 928 case subjects and 1039 control subjects and 246 case subjects and 516 control subjects, respectively. Percent density was estimated from digitized mammograms using the computer-assisted Cumulus thresholding program, and variation was estimated from an automated algorithm. We estimated hazards ratios (HRs), odds ratios (ORs), the area under the receiver operating characteristic curve (AUC), and 95% confidence intervals (CIs) using Cox proportional hazards models for the cohort and logistic regression for case-control studies, with adjustment for age and body mass index. We performed a meta-analysis using random study effects to obtain pooled estimates of the associations between the two mammographic measures and breast cancer. All statistical tests were two-sided.

RESULTS

The variation measure was statistically significantly associated with the risk of breast cancer in all three studies (highest vs lowest quartile: HR = 2.0 [95% CI = 1.3 to 3.1]; OR = 2.7 [95% CI = 2.1 to 3.6]; OR = 2.4 [95% CI = 1.4 to 3.9]; [corrected] all P (trend) < .001). [corrected]. The risk estimates and AUCs for the variation measure were similar to [corrected] those for percent density (AUCs for variation = 0.60-0.62 and [corrected] AUCs for percent density = 0.61-0.65). [corrected]. A meta-analysis of the three studies demonstrated similar associations [corrected] between variation and breast cancer (highest vs lowest quartile: RR = 1.8, 95% CI = 1.4 to 2.3) and [corrected] percent density and breast cancer (highest vs lowest quartile: RR = 2.3, 95% CI = 1.9 to 2.9).

CONCLUSION

The association between the automated variation measure and the risk of breast cancer is at least as strong as that for percent density. Efforts to further evaluate and translate the variation measure to the clinical setting are warranted.

[Computerized estimation of a percent glandular tissue composition in computed radiography mammography]

Measurement of a percent glandular tissue composition (%GTC) is important in terms of the estimation of individual patient exposure dose and the prediction of malignancy, and thus a number of reports for estimating %GTC by use of a mammogram have been published. In this study, we propose a method for estimating individual %GTC by use of computed radiography (CR) mammograms. By employing breast-equivalent phantoms that are able to create breast phantom images with various combinations of fat and glandular tissue, as well as the thickness of whole breast, we determined a reference table for converting an each pixel value on CR mammography to the glandular tissue ratio. Therefore, the %GTC for individual breast was estimated by averaging glandular tissue ratio for a whole region. The clinical image data set that consisted of 49 CR mammograms were used for estimating %GTC. A paired comparison method for determining subjective ranking of the degree of breast density was employed in order to demonstrate the validity of our method. The results indicate that the average estimated %GTC was 35.0% (ranged from 12.0% to 67.0%) and they had a increased correlation with the ranking of those obtained by observer test. Therefore, it was suggested that our proposed method would be utilized for estimating the %GTC in objective manner.

Full field digital mammography and breast density: comparison of calibrated and noncalibrated measurements

RATIONALE AND OBJECTIVES

Mammographic breast density is an important and widely accepted risk factor for breast cancer. A statement about breast density in the mammographic report is becoming a requirement in many States. However, there is significant inter-observer variation between radiologists in their interpretation of breast density. A properly designed automated system could provide benefits in maintaining consistency and reproducibility. We have developed a new automated and calibrated measure of breast density using full field digital mammography (FFDM). This new measure assesses spatial variation within a mammogram and produced significant associations with breast cancer in a small study. The costs of this automation are delays from advanced image and data analyses before the study can be processed. We evaluated this new calibrated variation measure using a larger dataset than previously. We also explored the possibility of developing an automated measure from unprocessed (raw data) mammograms as an approximation for this calibrated breast density measure.

MATERIALS AND METHODS

A case-control study comprised of 160 cases and 160 controls matched by age, screening history, and hormone replacement therapy was used to compare the calibrated variation measure of breast density with three variants of a noncalibrated measure of spatial variation. The operator-assisted percentage of breast density measure (PD) was used as a standard reference for comparison. Odds ratio (OR) quartile analysis was used to compare these measures. Linear regression analysis was applied to assess the calibration's impact on the raw pixel distribution.

RESULTS

All breast density measures showed significant breast cancer associations. The calibrated spatial variation measure produced the strongest associations (OR: 1.0 [ref.], 4.6, 4.3, 7.4). The associations for PD were diminished in comparison (OR: 1.0 [ref.], 2.7, 2.9, 5.2). Two additional non-calibrated measures restricted in region size also showed significant associations (OR: 1.0 [ref.], 2.9, 4.4, 5.4), and (OR: 1.0 [ref.], 3.5, 3.1, 4.9). Regression analyses indicated the raw image mean is influenced by the calibration more so than its standard deviation.

CONCLUSION

Breast density measures can be automated. The associated calibration produced risk information not retrievable from the raw data representation. Although the calibrated measure produced the stronger association, the non-calibrated measures may offer an alternative to PD and other operator based methods after further evaluation, because they can be implemented automatically with a simple processing algorithm.

Volume of mammographic density and risk of breast cancer

BACKGROUND

Assessing the volume of mammographic density might more accurately reflect the amount of breast volume at risk of malignant transformation and provide a stronger indication of risk of breast cancer than methods based on qualitative scores or dense breast area.

METHODS

We prospectively collected mammograms for women undergoing screening mammography. We determined the diagnosis of subsequent invasive or ductal carcinoma in situ for 275 cases, selected 825 controls matched for age, ethnicity, and mammography system, and assessed three measures of breast density: percent dense area, fibroglandular volume, and percent fibroglandular volume.

RESULTS

After adjustment for familial breast cancer history, body mass index, history of breast biopsy, and age at first live birth, the ORs for breast cancer risk in the highest versus lowest measurement quintiles were 2.5 (95% CI: 1.5-4.3) for percent dense area, 2.9 (95% CI: 1.7-4.9) for fibroglandular volume, and 4.1 (95% CI: 2.3-7.2) for percent fibroglandular volume. Net reclassification indexes for density measures plus risk factors versus risk factors alone were 9.6% (P = 0.07) for percent dense area, 21.1% (P = 0.0001) for fibroglandular volume, and 14.8% (P = 0.004) for percent fibroglandular volume. Fibroglandular volume improved the categorical risk classification of 1 in 5 women for both women with and without breast cancer.

CONCLUSION

Volumetric measures of breast density are more accurate predictors of breast cancer risk than risk factors alone and than percent dense area.

IMPACT

Risk models including dense fibroglandular volume may more accurately predict breast cancer risk than current risk models.

Calibrated measures for breast density estimation

RATIONALE AND OBJECTIVES

Breast density is a significant breast cancer risk factor measured from mammograms. Evidence suggests that the spatial variation in mammograms may also be associated with risk. We investigated the variation in calibrated mammograms as a breast cancer risk factor and explored its relationship with other measures of breast density using full field digital mammography (FFDM).

MATERIALS AND METHODS

A matched case-control analysis was used to assess a spatial variation breast density measure in calibrated FFDM images, normalized for the image acquisition technique variation. Three measures of breast density were compared between cases and controls: (a) the calibrated average measure, (b) the calibrated variation measure, and (c) the standard percentage of breast density (PD) measure derived from operator-assisted labeling. Linear correlation and statistical relationships between these three breast density measures were also investigated.

RESULTS

Risk estimates associated with the lowest to highest quartiles for the calibrated variation measure were greater in magnitude (odds ratios: 1.0 [ref.], 3.5, 6.3, and 11.3) than the corresponding risk estimates for quartiles of the standard PD measure (odds ratios: 1.0 [ref.], 2.3, 5.6, and 6.5) and the calibrated average measure (odds ratios: 1.0 [ref.], 2.4, 2.3, and 4.4). The three breast density measures were highly correlated, showed an inverse relationship with breast area, and related by a mixed distribution relationship.

CONCLUSION

The three measures of breast density capture different attributes of the same data field. These preliminary findings indicate the variation measure is a viable automated method for assessing breast density. Insights gained by this work may be used to develop a standard for measuring breast density.

A quantitative description of the percentage of breast density measurement using full-field digital mammography

RATIONALE AND OBJECTIVES

Breast density is a significant breast cancer risk factor that is measured from mammograms. However, uncertainty remains in both understanding its underlying physical properties as it relates to the breast and determining the optimal method for its measurement. A quantitative description of the information captured by the standard operator-assisted percentage of breast density (PD) measure was developed using full-field digital mammography (FFDM) images that were calibrated to adjust for interimage acquisition technique differences.

MATERIALS AND METHODS

The information captured by the standard PD measure was quantified by developing a similar measure of breast density (PD(c)) from calibrated mammograms automatically by applying a static threshold to each image. The specific threshold was estimated by first sampling the probability distributions for breast tissue in calibrated mammograms. A percent glandular (PG) measure of breast density was also derived from calibrated mammograms. The PD, PD(c), and PG breast density measures were compared using both linear correlation (R) and quartile odds ratio measures derived from a matched case-control study.

RESULTS

The standard PD measure is an estimate of the number of pixel values above a fixed idealized x-ray attenuation fraction. There was significant correlation (P < .0001) between the PD(c)-PD (r = 0.78), PD(c)-PG (r = 0.87), and PD-PG (r = 0.71) measures of breast density. Risk estimates associated with the lowest to highest quartiles for the PD(c) measure (odds ratio [OR]: 1.0 ref., 3.4, 3.6, and 5.6), and the standard PD measure (OR 1.0 ref., 2.9, 4.8, and 5.1) were similar and greater than that of the calibrated PG measure (OR 1.0 ref., 2.0, 2.4, and 2.4).

CONCLUSIONS

The information captured by the standard PD measure was quantified as it relates to calibrated mammograms and used to develop an automated method for measuring breast density. These findings represent an initial step for developing an automated measure built on an established calibration platform. A fully developed automated measure may be useful for both research- and clinical-based risk applications.

Statistical analysis of mammographic breast composition measurements: towards a quantitative measure of relative breast cancer risk

OBJECTIVE

A number of studies have identified the relationship between the visual appearance of high breast density at mammography and an increased risk of breast cancer. With the advent of digital mammography and the promise of routine measurements of parameters associated with breast composition, the possibility arises of using breast composition in a quantitative manner to predict relative breast cancer risk. Previous measurements have shown that the average proportion of glandular and adipose tissue within the breast varies with both age and breast size. In order to be able to identify individual women with an unusually high volume of glandular tissue, it will therefore be necessary to make comparisons with a disease-free population matched for age and breast size.

METHODS

A large number of breast glandular thickness measurements were analysed to investigate the statistics of breast composition across a disease-free population as a test of a suitable methodology for relative risk estimation. The large data set is also used to revisit the trends in breast composition used in the current UK method of breast radiation dosimetry.

RESULTS

It is demonstrated that a non-linear transformation can be used to produce normal statistical distributions, suitable for producing a standardised "Z-score" for breast composition.

CONCLUSION

A standard "Z-score" approach to identify women with unusually glandular breasts is recommended and so provide a basis for cancer risk estimations.

Effective radiation attenuation calibration for breast density: compression thickness influences and correction

Background

Calibrating mammograms to produce a standardized breast density measurement for breast cancer risk analysis requires an accurate spatial measure of the compressed breast thickness. Thickness inaccuracies due to the nominal system readout value and compression paddle orientation induce unacceptable errors in the calibration.

Method

A thickness correction was developed and evaluated using a fully specified two-component surrogate breast model. A previously developed calibration approach based on effective radiation attenuation coefficient measurements was used in the analysis. Water and oil were used to construct phantoms to replicate the deformable properties of the breast. Phantoms consisting of measured proportions of water and oil were used to estimate calibration errors without correction, evaluate the thickness correction, and investigate the reproducibility of the various calibration representations under compression thickness variations.

Results

The average thickness uncertainty due to compression paddle warp was characterized to within 0.5 mm. The relative calibration error was reduced to 7% from 48-68% with the correction. The normalized effective radiation attenuation coefficient (planar) representation was reproducible under intra-sample compression thickness variations compared with calibrated volume measures.

Conclusion

Incorporating this thickness correction into the rigid breast tissue equivalent calibration method should improve the calibration accuracy of mammograms for risk assessments using the reproducible planar calibration measure.

Validation of a method for measuring the volumetric breast density from digital mammograms

The purpose of this study was to evaluate the performance of an algorithm used to measure the volumetric breast density (VBD) from digital mammograms. The algorithm is based on the calibration of the detector signal versus the thickness and composition of breast-equivalent phantoms. The baseline error in the density from the algorithm was found to be 1.25 +/- 2.3% VBD units (PVBD) when tested against a set of calibration phantoms, of thicknesses 3-8 cm, with compositions equivalent to fibroglandular content (breast density) between 0% and 100% and under x-ray beams between 26 kVp and 32 kVp with a Rh/Rh anode/filter. The algorithm was also tested against images from a dedicated breast computed tomography (CT) scanner acquired on 26 volunteers. The CT images were segmented into regions representing adipose, fibroglandular and skin tissues, and then deformed using a finite-element algorithm to simulate the effects of compression in mammography. The mean volume, VBD and thickness of the compressed breast for these deformed images were respectively 558 cm(3), 23.6% and 62 mm. The displaced CT images were then used to generate simulated digital mammograms, considering the effects of the polychromatic x-ray spectrum, the primary and scattered energy transmitted through the breast, the anti-scatter grid and the detector efficiency. The simulated mammograms were analyzed with the VBD algorithm and compared with the deformed CT volumes. With the Rh/Rh anode filter, the root mean square difference between the VBD from CT and from the algorithm was 2.6 PVBD, and a linear regression between the two gave a slope of 0.992 with an intercept of -1.4 PVBD and a correlation with R(2) = 0.963. The results with the Mo/Mo and Mo/Rh anode/filter were similar.

X-ray absorptiometry of the breast using mammographic exposure factors: application to units featuring automatic beam quality selection

A number of studies have identified the relationship between the visual appearance of high breast density at mammography and an increased risk of breast cancer. Approaches to quantify the amount of glandular tissue within the breast from mammography have so far concentrated on image-based methods. Here, it is proposed that the X-ray parameters automatically selected by the mammography unit can be used to estimate the thickness of glandular tissue overlying the automatic exposure sensor area, provided that the unit can be appropriately calibrated. This is a non-trivial task for modern mammography units that feature automatic beam quality selection, as the number of tube potential and X-ray target/filter combinations used to cover the range of breast sizes and compositions can be large, leading to a potentially unworkable number of curve fits and interpolations. Using appropriate models for the attenuation of the glandular breast in conjunction with a constrained set of physical phantom measurements, it is demonstrated that calibration for X-ray absorptiometry can be achieved despite the large number of possible exposure factor combinations employed by modern mammography units. The main source of error on the estimated glandular tissue thickness using this method is shown to be uncertainty in the measured compressed breast thickness. An additional correction for this source of error is investigated and applied. Initial surveys of glandular thickness for a cohort of women undergoing breast screening are presented.

Screen-film mammographic density and breast cancer risk: a comparison of the volumetric standard mammogram form and the interactive threshold measurement methods

BACKGROUND

Mammographic density is a strong risk factor for breast cancer, usually measured by an area-based threshold method that dichotomizes the breast area on a mammogram into dense and nondense regions. Volumetric methods of breast density measurement, such as the fully automated standard mammogram form (SMF) method that estimates the volume of dense and total breast tissue, may provide a more accurate density measurement and improve risk prediction.

METHODS

In 2000-2003, a case-control study was conducted of 367 newly confirmed breast cancer cases and 661 age-matched breast cancer-free controls who underwent screen-film mammography at several centers in Toronto, Canada. Conditional logistic regression was used to estimate odds ratios of breast cancer associated with categories of mammographic density, measured with both the threshold and the SMF (version 2.2beta) methods, adjusting for breast cancer risk factors.

RESULTS

Median percent density was higher in cases than in controls for the threshold method (31% versus 27%) but not for the SMF method. Higher correlations were observed between SMF and threshold measurements for breast volume/area (Spearman correlation coefficient = 0.95) than for percent density (0.68) or for absolute density (0.36). After adjustment for breast cancer risk factors, odds ratios of breast cancer in the highest compared with the lowest quintile of percent density were 2.19 (95% confidence interval, 1.28-3.72; P(t) <0.01) for the threshold method and 1.27 (95% confidence interval, 0.79-2.04; Pt = 0.32) for the SMF method.

CONCLUSION

Threshold percent density is a stronger predictor of breast cancer risk than the SMF version 2.2beta method in digitized images.

Quantification of breast density with dual energy mammography: an experimental feasibility study

PURPOSE

Breast density, the percentage of glandular breast tissue, has been shown to be a strong indicator of breast cancer risk. A quantitative method to measure breast density with dual energy mammography was investigated using physical phantoms.

METHODS

The dual energy mammography system used a tungsten anode x-ray tube with a 50 microm rhodium beam filter for low energy images and a 300 microm copper beam filter for high energy images. Glandular and adipose equivalent phantoms of uniform thickness were used to calibrate a dual energy basis decomposition algorithm. Four different phantom studies were used to evaluate the technique. The first study consisted of phantoms with thicknesses of 2.5-8.5 cm in 0.5 cm steps with variable densities centered at a mean of 28%. The second study consisted of phantoms at a fixed thickness of 4.0 cm, which ranged in densities from 0% to 100% in increments of 12.5%. The third study consisted of 4.0 cm thick phantoms at densities of 25%, 50% and 75% each imaged at three areal sizes, approximately 62.5, 125, and 250 cm2, in order to assess the effect of breast size on density measurement. The fourth study consisted of step phantoms designed to more closely mimic the shape of a female breast with maximal thicknesses from 3.0 to 7.0 cm at a fixed density of 50%. All images were corrected for x-ray scatter.

RESULTS

The RMS errors in breast density measurements were 0.44% for the variable thickness phantoms, 0.64% for the variable density phantoms, 2.87% for the phantoms of different areal sizes, and 4.63% for step phantoms designed to closely resemble the shape of a breast.

CONCLUSIONS

The results of the phantom studies indicate that dual energy mammography can be used to measure breast density with an RMS error of approximately 5%.

The myth of the 50-50 breast

PURPOSE

For dosimetry and for work in optimization of x-ray imaging of the breast, it is commonly assumed that the breast is composed of 50% fibroglandular tissue and 50% fat. The purpose of this study was to assess whether this assumption was realistic.

METHODS

First, data obtained from an experimental breast CT scanner were used to validate an algorithm that measures breast density from digitized film mammograms. Density results obtained from a total of 2831 women, including 191 women receiving CT and from mammograms of 2640 women from three other groups, were then used to estimate breast compositions.

RESULTS

Mean compositions, expressed as percent fibroglandular tissue (including the skin), varied from 13.7% to 25.6% among the groups with an overall mean of 19.3%. The mean compressed breast thickness for the mammograms was 5.9 cm (sigma = 1.6 cm). 80% of the women in our study had volumetric breast density less than 27% and 95% were below 45%.

CONCLUSIONS

Based on the results obtained from the four groups of women in our study, the "50-50" breast is not a representative model of the breast composition.

Cumulative sum quality control for calibrated breast density measurements

PURPOSE

Breast density is a significant breast cancer risk factor. Although various methods are used to estimate breast density, there is no standard measurement for this important factor. The authors are developing a breast density standardization method for use in full field digital mammography (FFDM). The approach calibrates for interpatient acquisition technique differences. The calibration produces a normalized breast density pixel value scale. The method relies on first generating a baseline (BL) calibration dataset, which required extensive phantom imaging. Standardizing prospective mammograms with calibration data generated in the past could introduce unanticipated error in the standardized output if the calibration dataset is no longer valid.

METHODS

Sample points from the BL calibration dataset were imaged approximately biweekly over an extended timeframe. These serial samples were used to evaluate the BL dataset reproducibility and quantify the serial calibration accuracy. The cumulative sum (Cusum) quality control method was used to evaluate the serial sampling.

RESULTS

There is considerable drift in the serial sample points from the BL calibration dataset that is x-ray beam dependent. Systematic deviation from the BL dataset caused significant calibration errors. This system drift was not captured with routine system quality control measures. Cusum analysis indicated that the drift is a sign of system wear and eventual x-ray tube failure.

CONCLUSIONS

The BL calibration dataset must be monitored and periodically updated, when necessary, to account for sustained system variations to maintain the calibration accuracy.