# Effective x-ray attenuation coefficient measurements from two full field digital mammography systems for data calibration applications

- John J Heine
^{1}Email author and - Jerry A Thomas
^{2}

**7**:13

https://doi.org/10.1186/1475-925X-7-13

© Heine and Thomas; licensee BioMed Central Ltd. 2008

**Received: **06 November 2007

**Accepted: **28 March 2008

**Published: **28 March 2008

## Abstract

### Background

Breast density is a significant breast cancer risk factor. Currently, there is no standard method for measuring this important factor. Work presented here represents an essential component of an ongoing project that seeks to determine the appropriate method for calibrating (standardizing) mammography image data to account for the x-ray image acquisition influences. Longer term goals of this project are to make accurate breast density measurements in support of risk studies.

### Methods

Logarithmic response calibration curves and effective x-ray attenuation coefficients were measured from two full field digital mammography (FFDM) systems with breast tissue equivalent phantom imaging and compared. Normalization methods were studied to assess the possibility of reducing the amount of calibration data collection. The percent glandular calibration map functional form was investigated. Spatial variations in the calibration data were used to assess the uncertainty in the calibration application by applying error propagation analyses.

### Results

Logarithmic response curves are well approximated as linear. Measured effective x-ray attenuation coefficients are characteristic quantities independent of the imaging system and are in agreement with those predicted numerically. Calibration data collection can be reduced by applying a simple normalization technique. The calibration map is well approximated as linear. Intrasystem calibration variation was on the order of four percent, which was approximately half of the intersystem variation.

### Conclusion

FFDM systems provide a quantitative output, and the calibration quantities presented here may be used for data acquired on similar FFDM systems.

## Keywords

## 1. Background

Early detection is a key element in reducing breast cancer mortality [1]. Mammography screening is an essential surveillance component for early detection [2]. Similarly, there is interest in developing total cancer care methods in clinical practice so that disease screening and treatment can be tailored to the patient [3]. The development of accurate breast cancer risk models may play an important role in designing risk based cancer control strategies. Because breast density is a significant breast cancer risk factor [4], it may be useful to include it in the clinical setting for risk assessment. The Gail breast cancer risk model is used for intervention studies and counseling [5] but does not include breast density beyond research purposes. There is a critical need to incorporate all available information on breast cancer risk to ensure that risk models are useful for clinical decision making [6].

Mammographic density and breast density are terms used synonymously to describe the degree of bright areas in mammograms, which is related to dense breast tissue. Substantial evidence indicates that women with greater quantities of breast density have a significantly increased breast cancer risk [7]. Breast density may be assessed from either film mammograms or from the newer full field digital mammography (FFDM) systems. To date, most published breast density related research is based on film analysis without applying standardization [4, 8], which may be due to the availability of archived film data relative to FFDM data. The overall equivalence of film mammography and FFDM for breast cancer detection [9] may indicate that the FFDM system use will increase, which is the case at this center.

Methods used for assessing breast density may be loosely grouped into (1) techniques that do not consider the acquisition influences, and (2) radiometric standardization techniques that compensate for the acquisition influences. Standardization (calibration) methods are to correct for interpatient variations in the x-ray exposure, beam type, compression height, and the detector response. Standardization techniques are under development for both film [10–12] and FFDM applications [13]. To the best of our knowledge, there are no published reports on calibrated techniques with risk assessments using breast cancer as the endpoint. A non-calibration user-assisted approach has been shown consistently to produce a measure that correlates with breast cancer [4, 8], which produces a binary labeled image defining dense and fatty tissue. However, there is no universal standard used for measuring breast density [14].

Our previous work showed that it is possible to make effective attenuation coefficient measurements using the General Electric (GE) FFDM system with phantom imaging [8]. In this prior work, theoretical arguments were developed without considering the detector interaction or specific x-ray spectral form, which predicted linear logarithmic response curves. Idealized x-ray equations were then used with simulated x-ray spectra, tabulated attenuation data, and numerical integration to generate logarithmic response curves. These numerically generated curves were well approximated as linear and were used to make effective attenuation coefficient estimates, which were then compared with measured coefficient values obtained experimentally with phantom imaging. This work suggested that the effective attenuation coefficients may be used for calibration purposes as known quantities. This previous work did not provide measured effective attenuation coefficient values for glandular or fatty breast tissue.

Findings presented here represent essential components of an ongoing project for developing calibration techniques for FFDM breast cancer risk applications [8] that extends earlier research [13]. This project includes measuring the intra/inter system concordance of calibration data generated with phantom imaging using two similar FFDM systems. Serial stability comparisons will be analyzed after collecting sufficient data over an extended period. The overall goals of this multiyear project are to (1) determine the necessary amount of calibration data that must be collected (sampled) to calibrate the prospective data with a given accuracy using breast cancer as the endpoint, (2) analyze the serial stability of the systems by sampling the calibration data, and (3) study the wider applicability of calibration data generated on a given system by assessing the necessary requirements to merge data from multiple systems. The serial stability analysis is in support of a multiyear case control study that involves collecting prospective patient image/risk data. By hypothesis, this data will be used to validate the calibration procedure by showing that it produces a stronger risk relation when compared with methods that do not consider the acquisition influences when applied to the same dataset.

Specifically in this report, effective attenuation coefficients and breast tissue equivalent calibration curves generated with two similar FFDM systems were compared for the three filter\target combinations and multiple voltage settings. This may be considered the baseline (BL) data collection and analysis for the ongoing serial study. Effective attenuation coefficients were also compared with values generated from numerical integration by using idealized x-ray attenuation equations. Methods for generating these quantities were described previously [8] and will not be repeated in detail here. A technique for reducing the amount of calibration data collection was investigated.

The percent glandular (*PG*) equivalent transform (map) construction and functional form were investigated with a specific acquisition example and numerical methods. In practice, the *PG* map will be used to calibrate arbitrary image data. Tissue equivalent logarithmic response curve analyses enabled the map construction. Uncertainty in the *PG* map was estimated with error propagation techniques.

## 2. Methods

*y*direction (from top to bottom), and the distance from left hand border along the

*x*direction to the furthest distance to the semicircular perimeter is 13 cm.

Baseline tissue equivalent logarithmic response curves were measured for both locations. Response curves were generated by measuring the response for a particular acquisition technique (for a given voltage setting, target\filter combination, and exposure) as a function of phantom configuration height above the detector. These response curves permit the calibration of data collected with other acquisition techniques. Molybdenum/Molybdenum (Mo/Mo), Molybdenum/Rhodium (Mo/Rh), and Rhodium/Rhodium (Rh/Rh) target\filter combinations were used to generate the response curves. The voltage setting abbreviation used here corresponds with the peak voltage (kVp) system console indicator. Each tissue type response curve was generated for these settings: Mo/Mo – [26,27] kVp; Mo/Rh – [28, 29] kVp; and Rh/Rh – [30 kVp]. All BL response curves were generated with the 160 milliamp second beam-current system setting, which corresponds to the system console *mAs* indicator. The system *mAs* value was used as surrogate for the incident exposure. Phantoms were imaged either with [2–6] cm, [2–7] cm, or with [3–7] cm heights depending on the target\filter combination using one cm height increments. The 160 *mAs* setting was used because it provided signal when imaging the more attenuating [6–7] cm glandular configuration while not saturating the intensity values for less attenuating [2–3] cm fatty phantom arrangements. All image data was acquired as left cranial caudal views. Florida BL data was acquired in October 2006, and KS data was acquired in November 2006. The compression paddle was in place and in contact with the top phantom for all configurations while imaging to mimic patient imaging. The flat field quality assurance procedure (assessed by the built-in system diagnostics) was implemented before imaging to ensure that the system was operating within the manufacturer's tolerances.

An 800 × 1400 pixel rectangular region was used for the analysis, which defined the field of view (FOV). The aim was to restrict the comparison analysis to the same detector regions for both systems. The FOV, which is shown in Fig. 1 and Fig. 2, was dictated by the KS phantom geometry. The FOV is approximately the largest rectangle that fits within the semicircular geometry with the left hand margin.

Logarithmic response (*LR*) curves were determined for each 50 × 50 pixel sector within the FOV as a function of phantom height and averaged for a given tissue type. Each sector was fitted with this linear model
*LR*(*t*) = -*μt* + *l*,

where *t* is a given phantom configuration height measured in cm, μ is the effective x-ray attenuation coefficient for a given tissue type measured in cm^{-1}, and *l* is the logarithmic intercept value (unitless) for a given tissue type. Within a given sector, the pixel values (*V*) were averaged and natural logarithm transformed. This procedure was applied to each phantom for given height, given sector, and then for every sector and height. Plotting ln [mean(*V*)/*mAs* ], where *mAs* is the system acquisition readout, as a function of phantom height above the breast support surface (bucky) gives the empirical *LR* for a given square. The *mAs* normalization is discussed below in more detail. The sectored analysis resulted in a distribution of 16 × 28 values for each of the two model parameters for each tissue type *LR* curve. Findings will be presented as the mean and standard deviation for the various distributions. Confidence intervals (CIs) were generated for the regression parameters. For one example, the *LR* plot will be provided for one individual sector centered on the detector in the *y* direction positioned on the left-hand side of FOV (approximately under the focal point). The R^{2} statistic was used to assess the linear model agreement, because it is the fraction of variability explained by the model (unity implies perfect agreement) and therefore represents an important goodness of fit measure derived from the residual analysis. The distribution of R^{2} coefficients was generated and summarized to assess the linear model agreement. Effective attenuation coefficients were derived with numerical methods discussed previously [8] and compared with the corresponding measured values.

*V*=

*m*(

*mAs*) +

*b*, where

*m*and

*b*are constants for a given beam-type (fixed kVp and target\filter settings) [8]. This was achieved by acquiring exposures of the detector (no attenuation) for

*mAs*settings up to the point of detector saturation (not including saturation). Response curves acquired with phantom imaging were analyzed by first transforming

*V*to the equivalent

*mAs*value for a arbitrary acquisition by inverting the detector response:

*mAs*equivalent = (

*V*-

*b*)

*m*

^{-1}. The equivalent

*mAs*value was then normalized by the acquisition

*mAs*and logarithmically transformed giving the logarithmic relative exposure relation:

*LRE*= ln[(

*V*-

*b*) (

*mAs*×

*m*)

^{-1}]. The response generated with the reference

*mAs*was then equated with the response generated for arbitrary

*mAs*value giving

*V*>

*b*, this relation reduces to

where *V* is the average sector value. When the Eq. (3) approximation holds, the need to characterize the detector response is eliminated. Moreover, Eq. (3) suggests that for fixed voltage setting and target\filter combination, a response curve generated with one *mAs* value is equivalent to all other curves acquired with varying *mAs* values. In this work, Eq. (4) was measured as a function of phantom height and fitted with the linear model defined in Eq. (1) for both tissue types.

For the *mAs* normalization analysis, additional phantom imaging was performed at the FL location (only) to generate the Mo/Mo – 26 kVp glandular tissue *LR* curves for 160 *mAs* and 200 *mAs* acquisitions to assess Eq. (3). Additional data was also acquired to generate the fatty tissue *LR* curves for 160 *mAs* and 110 *mAs* for the same voltage setting and target\filter. In this analysis, *mAs*
_{r} = 160. This additional data collection was used to eliminate the possibility of serial drift influence in the normalization analysis. Mo/Mo – 26 kVp acquisition parameters are approximately the average settings used in clinical practice [16].

To assess the linear calibration model, the calibration procedure derived from the *LR* regression parameters was compared with the measured calibration curve for a fixed height of 6 cm (*T* = 6 cm total height). For a given acquisition technique, two calibration points are required [13] to calibrate arbitrary image data. These points correspond to the *LR* value for each tissue type for a given height, which can be estimated from Eq. (1) with *t* = *T*. In practice, the detected *LR* value is likely to fall between the endpoints. For fixed height and arbitrary logarithmic response, *LR*
_{a} = *z*, the *PG* map is given by
*PG* = *Mz* + *B*.

(5)

The validity of the linear map form is addressed in the Appendix. The *M* and *B* parameters can be expressed as functions of the *LR* regression parameters by using *PG* = 100% and *PG* = 0% as boundary conditions with Eq. (1), which gives
*M* = 100 × [(*μ*
_{
g
}- *μ*
_{
f
})*T* + (*l*
_{
g
}- *l*
_{
f
})]^{-1}

where the *g* and *f* subscripts denote glandular and fatty tissue values, respectively. Equation 1 was modified to derive Eq. (6) with this form: *LR* = *μt* + *l*. A sequence of five phantom images was acquired to verify Eq. (6). This sequence represents a surrogate composite fatty/glandular tissue mixture for *T* = 6 cm corresponding to the known [0, 16.6, 50, 83.3, 100] *PG* mixtures. The known *PG* is determined from this relation: (height of the 100% glandular tissue component)/*T* × 100% for a given phantom image in the sequence. In [glandular, fatty] cm component phantom heights, these proportions correspond to composite heights of [0,6], [1,5], [3,3], [5,1], and [6,0], where the two component heights within a given bracket add to 6 cm. The sequence represents a surrogate because the various mixtures are not present in one image but in five separate images acquired consecutively. As a specific example, the Mo/Mo – 26 kVp acquisition settings were used to generate the sequence. Plotting *z* (see Eq. (5)) along the *x* axis and plotting the known *PG* values along the *y* axis shows the map. The agreement between *PG* = *Mz* + *B* using the *LR* regression distribution quantities substituted into Eq. (6) was compared with the known and measured *PG* determined from the surrogate composite phantom images.

For the map analysis, *LR* s were generated again for [2–6] cm heights using the same phantom imaging procedure outlined above by acquiring another dataset not included in the BL acquisition. This additional data collection will show the best scenario because the data (the phantom images for the tissue *LR* s and *PG* curves) were collected on the same day in January 2007 to eliminate the possibility of serial drift uncertainty. A sector example will also (the same sector as above) be used for demonstration purposes.

*PG*=

*f*(μ

_{g}, μ

_{f},

*l*

_{g},

*l*

_{f}, z) giving

*T*, is considered as exact and therefore does not contribute. The quadratic form used in Eq. (7) assumes each component contribution is independent. Equation (7) was used to estimate both intersystem variability (external variation) and intrasystem variability (internal variation). The total variation is given by

where each term in Eq. (7) has been relabeled, respectively. Intersystem analysis was used to derive the Δ terms for the external variation, which provided an estimate of the uncertainty when applying the map generated at one facility for data collected at another location. The internal variation was estimated with the FL spatial *LR* analysis. For both internal and external analyses, z and *Δz* were estimated for each *PG* surrogate component (6 cm height) as the average and standard deviation of log-transformed and *mAs* normalized FOV pixel distribution, respectively. Although *z* and Δz were estimated with FL data, they support reasonable estimates of the anticipated uncertainty due to the variations in the other *PG* map parameters for both uncertainty forms. The 26 kVp – Mo/Mo acquisition was used for the uncertainty analysis.

## 3. Results

### 3. 1 Inter and intra system comparisons

*l*values vary somewhat relative to attenuation coefficient variation. Agreement may be gauged from the parenthetical entries and CIs. Parenthetical quantities listed in the attenuation coefficient and intercept columns are absolute values of the percent difference between the location parameters. For example, the glandular quantity was calculated as |(μ

_{gFL}- μ

_{gKS})/μ

_{gKS}| × 100%, which was used for all calculations of this type. For either location, the intrasystem variation is relatively small for the regression parameters, which follows from the σ and σ

_{l}column entries. Glandular tissue attenuation coefficient variation is larger than the corresponding fatty tissue variation as indicated in the Δ column. Signal is attenuated greater at a given height for the glandular tissue relative to the fatty tissue indicating there may be greater uncertainty due to the decreased signal. The R

^{2}statistic for all cases is close to unity indicating agreement with the linear model.

Molybdenum/Molybdenum regression distribution quantities. Quantities are summaries of the field of view: μ is the mean effective x-ray attenuation coefficient, σ is the standard deviation of the μ distribution, Δ is the mean variation in μ derived from averaging the individual regression variation, *l* is the mean logarithmic intercept, σ_{
l
}is the standard deviation of the *l* distribution, and R^{2} is the mean of the associated distribution. Two sided 95% mean confidence interval is cited under the respective μ and *l* quantities. Location is indicated by the subscripts of the left hand column entries. Parenthetical FL entries indicate the absolute value percent difference between respective location quantities.

Mo/Mo 26 kVp | μ | σ | Δ |
| σ | R |
---|---|---|---|---|---|---|

Fat | 0.572 | 0.001 | 0.005 | 5.093 | 0.010 | 0.998 |

0.557–0.586 (1.9%) | 5.03–5.15 (2.6%) | |||||

Fat | 0.583 | 0.001 | 0.005 | 4.964 | 0.010 | 0.998 |

0.568–0.599 | 4.89–5.03 | |||||

Gland | 0.833 | 0.002 | 0.015 | 4.893 | 0.001 | 0.998 |

0.784–0.882 (0.60%) | 4.68–5.10 (2.6%) | |||||

Gland | 0.838 | 0.010 | 0.012 | 4.767 | 0.020 | 0.998 |

0.799–0.877 | 4.60–4.93 | |||||

Mo/Mo 27 kVp | μ | σ | Δ |
| σ | R |

Fat | 0.559 | 0.001 | 0.005 | 5.262 | 0.010 | 0.999 |

0.544–0.575 (2.1%) | 5.19–5.32 (1.8%) | |||||

Fat | 0.571 | 0.001 | 0.005 | 5.165 | 0.010 | 0.999 |

0.554–0.586 | 5.09–5.16 | |||||

Gland | 0.805 | 0.003 | 0.020 | 5.043 | 0.010 | 0.998 |

0.750–0.860 (0.40%) | 4.81–5.28 (1.9%) | |||||

Gland | 0.808 | 0.005 | 0.014 | 4.950 | 0.020 | 0.999 |

0.763–0.854 | 4.75–5.14 |

^{2}is not provided because it was very close to zero in all cases. Individual sector

*LR*example plots for the Mo/Mo – 26 kVp settings are shown in Fig. 3 (glandular) and Fig. 4 (fatty). Plots for a given tissue type are shown in the same figure for intersystem comparisons. In these plots, measured points from the FL location are denoted by diamonds and asterisks for the KS measured points. Fitted lines from the regression analysis are solid for the FL plots and dashed for KS plots. Regression parameters are listed in Table 2. Regression parameter distribution quantities for the other Mo/Rh and Rh/Rh settings are shown in Table 3. The corresponding individual sector plot examples for these other settings are not shown because they were similar to the prior examples. For all distribution examples, variation within given location is relatively small, and the effective attenuation coefficients are in close agreement between the locations. Moreover, the linear model approximation holds when assessing the R

^{2}statistic for all situations. To show the intersystem agreement, a paired t-test was applied by pairing all of the 10 average attenuation coefficient quantities for the FL location with the corresponding values from the KS location (second column from Table 1 and Table 3), which gave

*p*= 0.56 (linear correlation of = 0.97). Following the same procedure for the 10 average intercept quantities (using the fifth columns from the same tables), gave

*p*= 0.101 (linear correlation = 0.99). Using the

*p*value as a gauge, effective attenuation coefficients as an ensemble are similar across the two systems. Intersystem intercept comparisons indicate a weak system dependence.

Molybdenum/Molybdenum individual sector example regression parameters. Labeling is analogous to Table 1 as applied to a single realization.

Mo/Mo 26 kVp | μ | Δ |
| R |
---|---|---|---|---|

Fat | 0.573 | 0.004 | 5.103 | 0.999 |

0.559–0.586 (2.4%) | 5.05–5.16 (2.4%) | |||

Fat | 0.587 | 0.005 | 4.99 | 0.999 |

0.573–0.600 | 4.92–5.04 | |||

Gland | 0.837 | 0.014 | 4.904 | 0.999 |

0.791–0.883 (0.12%) | 4.70–5.09 (2.4%) | |||

Gland | 0.845 | 0.011 | 4.788 | 0.999 |

0.809–0.882 | 4.59–4.98 |

Molybdenum/Rhodium and Rhodium/Rhodium regression distribution quantities (using the same format as in Table 1).

Mo/Rh 28 kVp | μ | σ | Δ | l | σ | R |
---|---|---|---|---|---|---|

Fat | 0.503 | 0.002 | 0.005 | 5.343 | 0.011 | 0.999 |

0.490–0.516 (1.9%) | 5.28–5.40 (0.1%) | |||||

Fat | 0.511 | 0.001 | 0.005 | 5.34 | 0.010 | 0.999 |

0.498–0.525 | 5.27–5.40 | |||||

Gland | 0.718 | 0.002 | 0.013 | 5.155 | 0.013 | 0.998 |

0.683–0.753 (1.0%) | 4.98–5.32 (1.0%) | |||||

Gland | 0.711 | 0.010 | 0.012 | 5.118 | 0.020 | 0.999 |

0.677–0.745 | 4.95–5.28 | |||||

Mo/Rh 29 kVp | ||||||

Fat | 0.495 | 0.002 | 0.004 | 5.47 | 0.011 | 0.999 |

0.482–0.507 (1.2%) | 5.41–5.53 (0.20%) | |||||

Fat | 0.501 | 0.001 | 0.004 | 5.48 | 0.010 | 0.999 |

0.489–0.513 | 5.42–5.54 | |||||

Gland | 0.700 | 0.002 | 0.014 | 5.274 | 0.014 | 0.998 |

0.661–0.739 (1.2%) | 5.08–5.46 (0.20%) | |||||

Gland | 0.692 | 0.010 | 0.013 | 5.263 | 0.020 | 0.998 |

0.655–0.729 | 5.08–5.44 | |||||

Rh/Rh 30 kVp | ||||||

Fat | 0.442 | 0.001 | 0.002 | 5.648 | 0.010 | 0.999 |

0.435–0.449 (0.21%) | 5.61–5.68 (1.1%) | |||||

Fat | 0.443 | 0.001 | 0.002 | 5.712 | 0.010 | 0.999 |

0.434–0.451 | 5.66–5.75 | |||||

Gland | 0.644 | 0.002 | 0.010 | 5.557 | 0.010 | 0.999 |

0.621–0.666 (2.4%) | 5.45–5.66 (1.0%) | |||||

Gland | 0.629 | 0.005 | 0.010 | 5.593 | 0.014 | 0.999 |

0.609–0.648 | 5.49–5.69 |

^{2}values are also provided for this modeling, which indicate the linear model approximation holds in the theoretical developments as well. Numerical effective attenuation values are listed in the μ

_{ T }column. The associated CIs indicate approximate agreement with the measured values. Measured value agreement is also indicated by the parenthetical percent differences cited with the respective

*μ*

_{T}values. Intercept quantities were irrelevant.

Effective attenuation coefficient comparisons. This shows the numerically derived attenuation coefficients, μ_{T}, the associated regression variation, Δ, and R^{2} values for the Mo/Mo, Mo/Rh, and Rh/Rh configurations for the same voltage settings (kVp) used for the phantom imaging. Attenuation coefficient confidence intervals are cited beneath the respective values. The <*μ* > column gives the averaged FL + KS measured attenuation coefficients for the respective tissue types. Percent differences between < μ > and *μ*
_{T} are listed parenthetically next to the respective *μ*
_{T} values. Other formatting is the same used above.

Mo/Mo |
| Δ | R | < |
---|---|---|---|---|

Fat | 0.595 (2.8) | 0.010 | 0.999 | 0.578 |

0.575–0.613 | ||||

Gland | 0.871 (4.1) | 0.020 | 0.998 | 0.835 |

0.818–0.923 | ||||

Fat | 0.585 (2.9) | 0.010 | 0.999 | 0.565 |

0.562–0.604 | ||||

Gland | 0.836 (3.4) | 0.021 | 0.997 | 0.807 |

0.777–0.894 | ||||

Mo/Rh | ||||

Fat | 0.529 (4.0) | 0.010 | 0.999 | 0.508 |

0.511–0.547 | ||||

Gland | 0.764 (6.8) | 0.014 | 0.998 | 0.713 |

0.772–0.807 | ||||

Fat | 0.516 (3.8) | 0.010 | 0.999 | 0.498 |

0.498–0.535 | ||||

Gland | 0.735 (5.3) | 0.020 | 0.998 | 0.696 |

0.690–0.779 | ||||

Rh/Rh | ||||

Fat | 0.462 (4.3) | 0.010 | 0.999 | 0.442 |

0.448–0.476 | ||||

Gland | 0.675 (5.7) | 0.010 | 0.999 | 0.636 |

0.647–0.703 |

### 3. 2 Response normalization

*LR*distribution parameters generated for the

*mAs*response normalization analysis using [160, 110]

*mAs*for the fatty tissue type and [160, 200]

*mAs*for the glandular tissue type. For each tissue type,

*mAs*

_{r}= 160. The quantities in Table 5 show that the Eq. (3) approximation holds. A comparison of the respective attenuation coefficients for the varying

*mAs*values indicates close agreement. This suggests the

*LRs*acquired for a given acquisition

*mAs*determine the

*LR*s for an arbitrary

*mAs*readout value.

Normalization regression parameters. This gives the logarithmic response curve fitted parameter Molybdenum/Molybdenum (Mo/Mo) distribution quantities for both tissue types using the beam-current readout normalization (*mAs*). Acquisition *mAs* system values are indicated by the subscript of the left column entries.

Mo/Mo 26 kVp | μ | σ | Δ | l | σ | R |
---|---|---|---|---|---|---|

Fat | 0.568 | 0.002 | 0.005 | 5.054 | 0.01 | 0.999 |

0.550–0.586 | 4.97–5.13 | |||||

Fat | 0.569 | 0.002 | 0.006 | 5.058 | 0.010 | 0.999 |

0.549–0.590 | 4.97–5.14 | |||||

Gland | 0.827 | 0.002 | 0.016 | 4.851 | 0.010 | 0.998 |

0.777–0.878 | 4.63–5.07 | |||||

Gland | 0.826 | 0.002 | 0.015 | 4.853 | 0.010 | 0.998 |

0.777–0.875 | 4.64–5.06 |

### 3.3 Percent glandular map

*PG*row of Table 6. Internal variation in

*M*is small relative to the best estimate of

*M*, spatial variations in both

*M*and

*B*are relatively small, and R

^{2}indicates the linear model approximation holds. The bottom row of Table 6 shows

*M*and

*B*estimated with Eq. (6) using the quantities from Table 5. Deviation of the measured

*PG*value from the known phantom configuration is discussed below. Quantities in Table 6 indicate that the measured map agrees with the model and is well approximated as linear, as predicted by the work in the Appendix; a related numerically generated map is shown in Fig. 5. As further demonstration, the measured map determined from the individual sector example is shown in Fig. 6. For this case,

*M*= -57.23 ± 1.23,

*B*= 93.93 with R

^{2}= 0.999. Both

*M*and

*B*values for the sector agree with the respective entries in Table 6.

Percent glandular map parameters. The percent glandular row shows distribution quantities for the Molybdenum/Molybdenum (Mo/Mo) map determined with surrogate mixture phantom imaging. The PG_{160} row gives *M* and *B* calculated with Eq. (7) using the 160 *mAs* regression parameters from Table 5 with Eq. (6).

Mo/Mo 26 kVp | M | σ | Δ | B | σ | R |
---|---|---|---|---|---|---|

PG | -57.86 | 0.31 | 1.15 | 94.45 | 0.77 | 0.998 |

-61.5 – -54.2 | 90.7–98.2 | |||||

PG | -57.07 | / | / | 93.76 | / | / |

### 3.4 Uncertainty analysis

*PG*components expressed in Eqs. (7–8) were generated for the intersystem variation (external) and the intrasystem (internal) variation for a Mo/Mo – 26 kVp example. External variations in the log intercepts and attenuation values were estimated by taking the absolute value difference between respective location quantities from Table 1. The 160

*mAs*generated effective attenuation coefficients and the log intercepts listed in Table 5 were used as the known values for both the external and internal variation calculations. Internal variation was estimated by using the 160

*mAs*spatial distribution values cited in Table 5. Internal variations in the attenuation coefficients and log intercepts were taken from the σ and σ

_{l}columns in Table 5. Estimated uncertainties for the regression analysis parameters are summarized in Table 7. Uncertainties are broken down into their respective components along with the totals in Table 8, which follow from Eqs. (7–8);

*z*distribution quantities are cited in the bottom rows. External variation is approximately twice that of the internal. External logarithmic intercept differences represent major contributions, which are cited in the n

_{3}and n

_{4}columns; as one component increases the other decreases and vice versa, which maintains the variation throughout the

*PG*composition. Related internal components show a similar trend but with smaller contributions. If a bias exists with this specific example, it will tend to over estimate the external variation because this dataset has at least equivalent if not more intersystem variability in comparison with the other datasets.

Estimated uncertainties for the Molybdenum/Molybdenum example.

Mo/Mo 26 kVp | Δμ | Δ | Δμ | Δ | ΔLR |
---|---|---|---|---|---|

external | 0.01 | 0.130 | 0.01 | 0.140 | 0.15 |

internal | 0.002 | 0.01 | .002 | 0.01 | 0.012 |

Percent glandular map absolute % variation. This shows the individual external error (top) and internal error (bottom) components for the Molybdenum/Molybdenum – 26 kVp map output for each mixture. Totals follow from Eq. (8). The last two rows list the field of view average and standard deviation for each composition.

PG external variation | 0 | 16.6 | 50 | 83.3 | 100 |
---|---|---|---|---|---|

n | 0.02 | 0.30 | 0.89 | 1.42 | 1.68 |

n | 3.8 | 3.1 | 1.8 | 0.64 | 0.07 |

n | 0.07 | 1.3 | 3.7 | 6.0 | 7.0 |

n | 7.4 | 6.0 | 3.5 | 1.3 | 0.14 |

n | 0.78 | 0.95 | 1.42 | 1.9 | 2.3 |

total external variation | 8.4 | 7.0 | 5.7 | 6.6 | 7.6 |

Internal variation | |||||

n | 0.02 | 0.03 | 0.89 | 1.42 | 1.68 |

n | 3.80 | 3.10 | 1.81 | 0.64 | 0.08 |

n | 0.0 | 0.10 | 0.30 | 0.47 | 0.56 |

n | 0.58 | 0.47 | 0.30 | 0.10 | 0.01 |

n | 0.78 | 0.95 | 1.42 | 1.90 | 2.3 |

total internal variation | 3.9 | 3.3 | 2.5 | 2.5 | 2.9 |

<z> | 1.659 | 1.329 | 0.735 | 0.188 | -0.738 |

Δz | 0.013 | 0.016 | 0.024 | 0.033 | 0.040 |

## 4. Discussion

*V*response for the KS system is about 1.6 (ratio of the slopes) that of the FL system. The KS detector saturates near 60

*mAs*, whereas the FL system saturates slightly below 100

*mAs*for this example. Intercept confidence intervals both span zero but the KS width is much wider than the FL width (about 114 in comparison with 451).

Detector response parameter comparison. This shows the detector response regression parameters (distribution quantities) for both detectors for the Molybdenum/Molybdenum – 26 kVp acquisition settings fitted to *V* = *m*(*mAs*)+*b*. Parenthetical entries are distribution standard deviations.

Slope ( | Location | Intercept ( | Intercept CI | R | |
---|---|---|---|---|---|

FL | 183.2 (1.22) | 180.9–185.4 | 58.8 (12.1) | -55.3–173.1 | 0.999 |

KS | 300.0 (3.36) | 286.5–313.5 | 256.1 (14.2) | -195.6–707.8 | 0.999 |

The work shows that there is relatively little spatial variation in the two regression parameters for a given location. The choice of sector size will require further analysis to determine the optimal size because it influences the variation. Without a clinical endpoint comparison (cancer/no cancer patients), the appropriate spatial resolution of the mapping cannot be determined. The optimal resolution will be determined by its ability to assign breast cancer risk after accruing sufficient patient data. Likewise, the normalization work further demonstrates that when the voltage setting and target/filter are held constant, only two calibration curves are necessary, which permits a considerable reduction of *mAs* calibration data collection (sampling). Numerical evidence developed in the Appendix shows that the map is well approximated as linear, which is supported by experimental *PG* analysis. Moreover, the calibration map generated by the calibration curve parameters is in close agreement with the measured values obtained from the surrogate composite phantom imaging, which is another indication that these FFDM systems provide a quantitative output. As the *PG* external variation analysis showed, intercepts are source of uncertainty in the map. A similar argument would apply to a given system if the system drifts in time relative to the calibration curves. If the regression parameters are stable over time and the intercept variation proves acceptable, the work suggests that calibration parameters determined on a given GE-FFDM system will apply to other similar systems without the need for additional phantom imaging if the variation estimated here is acceptable. To the best of our knowledge, the calibration values presented here are not yet available in the public domain. This analysis was limited to two like systems due to availability. However, the calibration quantities presented here may serve as reference for comparisons of related work in the field implemented on other FFDM systems or can be used in calibration research performed at other sites without additional imaging. These findings extend the previous work [8, 13] by demonstrating the relationship between the calibration map, the measured calibration points, and the associated uncertainty.

## 5. Conclusion

Data used for this study is insufficient to assess the natural serial variation in the regression parameters. Longer term goals are to compare the inter/intra serial uncertainty with this BL data. Control measures from the BL data have been sampled on a biweekly basis since October 2007. Future work will use this serial data to determine the serial variation and its impact on the calibration procedure when assigning breast cancer risk to patient data. If the calibration approach is serially stable in the forward direction, it implies that data collected previously can be calibrated. It is often thought that FFDM data does not suffer from the same technical difficulties associated with film data due to varying response curves for example. This assumption may be generally true but will require validation for quantitative calibration type measurements.

Ideally, generating calibration curves for all possible acquisition techniques should produce the most accurate data normalization, which is obviously not practical. Thus, some form of interpolation is necessary. In part, the longer term aim of this work is to minimize the amount of a priori collected calibration data while maintaining accurate standardization. Showing that the effective attenuation coefficients are characteristic quantities and demonstrating the *mAs* normalization represent incremental steps in this direction. The non-specialty center calibration application acceptance (or usage) may be enhanced if the calibration data collection is minimized. In analogous fashion, future work includes minimizing the number of sample points in the voltage setting space. The guiding principle is that if the calibration requires minimal experimental effort, it may gain wider usage in the future.

## Appendix

### A. 1 Error propagation

Explicit calculation for Eq. (8) is derived using *f* = *M*(*μ*
_{
g
}, *μ*
_{
f
}, *l*
_{
g
}, *l*
_{
f
})*z* + *B*(*μ*
_{
g
}, *μ*
_{
f
}, *l*
_{
g
}, *l*
_{
f
}), where *z* is an arbitrary *LR* measurement with *f* = *PG* and *LR* = *μt* + *l*, which is a modified form of Eq. (1). The following substitutions are also used
*G* = [(*μ*
_{
g
}+ *μ*
_{
f
})*T* + *l*
_{
g
}+ *l*
_{
f
}]

and
*H* = [(*μ*
_{
g
}- *μ*
_{
f
})*T* + (*l*
_{
g
}- *l*
_{
f
})].

*f*with respect to each of the four parameters and

*z*gives

These five partial derivative expressions were used to evaluate Eqs. (7–8).

### A. 2 Calibration map analysis

Numerical methods, idealizations, abbreviations, and quantities used previously [8] will be used to show the map is approximately linear. The normalized energy photon distribution is given by *p*(*ε*). The response for composite structure of *T* = total height corresponding to *t*
_{1} cm of glandular tissue with Δ*t* = *T*-*t*
_{1} cm of fatty tissue is given by
*z'* = ∫*p*(*ε*)exp[-*μ*
_{
g
}(*ε*)*t*
_{1} - *μ*
_{
f
}(*ε*)Δ*t*]*dε*.

Attenuation coefficients are implicit functions of the x-ray energy (see previous work for details). The known *PG* = 100 × *t*
_{1}/*T*. Performing numerical integration using the 26 kVp – Mo/Mo simulated spectra for *t*
_{1} = [1,2,3,4,5,6] cm with the known *PG* = [16.6, 33.3, 50.0, 66.6, 83.3, 100.0], gives R^{2} = 0.999 when applying regression analysis to the ordered pairs: [ln(*z'*), *PG*]. The fitted line is shown in Fig. 5 for this example (regression line parameters as irrelevant here). Linear agreement is important because it validates the map form. The same analysis was carried out for *T* = [2,3,6] cm for all filter/target combinations and kVp values used for the phantom imaging based *LR* analysis with same the *t*
_{1} values and *PG* sample points used to generate Fig. 6. For all combinations, the results were similar with R^{2} = 0.999, which validates the use of Eq. (5).

## Declarations

### Acknowledgements

The work was supported by NCI grant # R01CA114491. The authors would like to thank Christy Smallwood and all of the staff at the Lifetime Cancer Screening Center at the Moffitt Cancer Center for their support of the work. The Authors also thank Dr. Thomas Sellers and Madhusmita Behera for their review of the work.

## Authors’ Affiliations

## References

- Tabar L, Dean PB:
**Mammography and breast cancer: the new era.***Int J Gynaecol Obstet*2003,**82:**319–326. 10.1016/S0020-7292(03)00262-5View ArticleGoogle Scholar - Duffy SW, Smith RA, Gabe R, Tabar L, Yen AM, Chen TH:
**Screening for breast cancer.***Surg Oncol Clin N Am*2005,**14:**671–697. 10.1016/j.soc.2005.06.001View ArticleGoogle Scholar - Dalton WS:
**The "total cancer care" concept: linking technology and health care.***Cancer Control*2005,**12:**140–141.Google Scholar - McCormack VA, dos Santos Silva I:
**Breast density and parenchymal patterns as markers of breast cancer risk: a meta-analysis.***Cancer Epidemiol Biomarkers Prev*2006,**15:**1159–1169. 10.1158/1055-9965.EPI-06-0034View ArticleGoogle Scholar - Decarli A, Calza S, Masala G, Specchia C, Palli D, Gail MH:
**Gail model for prediction of absolute risk of invasive breast cancer: independent evaluation in the Florence-European Prospective Investigation Into Cancer and Nutrition cohort.***J Natl Cancer Inst*2006,**98:**1686–1693.View ArticleGoogle Scholar - Freedman AN, Seminara D, Gail MH, Hartge P, Colditz GA, Ballard-Barbash R, Pfeiffer RM:
**Cancer risk prediction models: a workshop on development, evaluation, and application.***J Natl Cancer Inst*2005,**97:**715–723.View ArticleGoogle Scholar - Boyd NF, Rommens JM, Vogt K, Lee V, Hopper JL, Yaffe MJ, Paterson AD:
**Mammographic breast density as an intermediate phenotype for breast cancer.***Lancet Oncol*2005,**6:**798–808. 10.1016/S1470-2045(05)70390-9View ArticleGoogle Scholar - Heine JJ, Behera M:
**Effective x-ray attenuation measurements with full field digital mammography.***Med Phys*2006,**33:**4350–4366. 10.1118/1.2356648View ArticleGoogle Scholar - Pisano ED, Gatsonis C, Hendrick E, Yaffe M, Baum JK, Acharyya S, Conant EF, Fajardo LL, Bassett L, D'Orsi C,
*et al*.:**Diagnostic performance of digital versus film mammography for breast-cancer screening.***N Engl J Med*2005,**353:**1773–1783. 10.1056/NEJMoa052911View ArticleGoogle Scholar - Highnam R, Brady M:
*Mammographic image analysis*. Dordrecht; Boston: Kluwer Academic Publishers; 1999.View ArticleGoogle Scholar - Pawluczyk O, Augustine BJ, Yaffe MJ, Rico D, Yang J, Mawdsley GE, Boyd NF:
**A volumetric method for estimation of breast density on digitized screen-film mammograms.***Med Phys*2003,**30:**352–364. 10.1118/1.1539038View ArticleGoogle Scholar - Shepherd JA, Herve L, Landau J, Fan B, Kerlikowske K, Cummings SR:
**Novel use of Single X-Ray Absoptiometry for Measuring Breast Density.***Technology in Cancer Research & Treatment*2005,**4:**173–182.View ArticleGoogle Scholar - Kaufhold J, Thomas JA, Eberhard JW, Galbo CE, Trotter DE:
**A calibration approach to glandular tissue composition estimation in digital mammography.***Med Phys*2002,**29:**1867–1880. 10.1118/1.1493215View ArticleGoogle Scholar - Couzin J:
**Breast cancer. Fine-tuning breast density measures.***Science*2005,**309:**1665. 10.1126/science.309.5741.1665View ArticleGoogle Scholar - Vedantham S, Karellas A, Suryanarayanan S, Albagli D, Han S, Tkaczyk EJ, Landberg CE, Opsahl-Ong B, Granfors PR, Levis I,
*et al*.:**Full breast digital mammography with an amorphous silicon-based flat panel detector: physical characteristics of a clinical prototype.***Med Phys*2000,**27:**558–567. 10.1118/1.598895View ArticleGoogle Scholar - Cooper VN 3rd, Oshiro T, Cagnon CH, Bassett LW, McLeod-Stockmann TM, Bezrukiy NV:
**Evaluation of detector dynamic range in the x-ray exposure domain in mammography: a comparison between film-screen and flat panel detector systems.***Med Phys*2003,**30:**2614–2621. 10.1118/1.1606450View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.