Optimization of x-ray spectra in digital mammography through Monte Carlo simulations

In this work, a Monte Carlo code was used to investigate the performance of different x-ray spectra in digital mammography, through a figure of merit (FOM), defined as FOM = CNR2/(ˉ)D(g), with CNR being the contrast-to-noise ratio in image and [Formula: see text] being the average glandular dose. The FOM was studied for breasts with different thicknesses t (2 cm ≤ t ≤ 8 cm) and glandular contents (25%, 50% and 75% glandularity). The anode/filter combinations evaluated were those traditionally employed in mammography (Mo/Mo, Mo/Rh, Rh/Rh), and a W anode combined with Al or K-edge filters (Zr, Mo, Rh, Pd, Ag, Cd, Sn), for tube potentials between 22 and 34 kVp. Results show that the W anode combined with K-edge filters provides higher values of FOM for all breast thicknesses investigated. Nevertheless, the most suitable filter and tube potential depend on the breast thickness, and for t ≥ 6 cm, they also depend on breast glandularity. Particularly for thick and dense breasts, a W anode combined with K-edge filters can greatly improve the digital technique, with the values of FOM up to 200% greater than that obtained with the anode/filter combinations and tube potentials traditionally employed in mammography. For breasts with t < 4 cm, a general good performance was obtained with the W anode combined with 60 μm of the Mo filter at 24-25 kVp, while 60 μm of the Pd filter provided a general good performance at 24-26 kVp for t = 4 cm, and at 28-30 and 29-31 kVp for t = 6 and 8 cm, respectively.     

Patient dose in digital mammography

In the present investigation, we analyze the dose of 5034 patients (20,137 images) who underwent mammographic examinations with a full-field digital mammography system. Also, we evaluate the system calibration by analyzing the exposure factors as a function of breast thickness. The information relevant to this study has been extracted from the image DICOM header and stored in a database during a 3-year period (March 2001-October 2003). Patient data included age, breast thickness, kVp, mAs, target/filter combination, and nominal dose values. Entrance surface air kerma (ESAK) without backscatter was calculated from the tube output as measured for each voltage used under clinical conditions and from the tube loading (mAs) included in the DICOM header. Mean values for the patient age and compressed breast thickness were 56 years (SD: 11) and 52 mm (SD: 13), respectively. The majority of the images was acquired using the STD (for standard) automatic mode (98%). The most frequent target/filter combination automatically selected for breast smaller than 35 mm was Mo/Mo (75%); for intermediate thicknesses between 35 and 65 mm, the combinations were Mo/Rh (54%) and Rh/Rh (38.5%); Rh/Rh was the combination selected for 91% of the cases for breasts thicker than 65 mm. A wide kVp range was observed for each target/filter combination. The most frequent values were 28 kVp for Mo/Mo, 29 kVp for Mo/Rh, and 29 and 30 kV for Rh/Rh. Exposure times ranged from 0.2 to 4.2 s with a mean value of 1.1 s. Average glandular doses (AGD) per exposure were calculated by multiplying the ESAK values by the conversion factors tabulated by Dance for women in the age groups 50 to 64 and 40 to 49. This approach is based on the dependence of breast glandularity on breast thickness and age. The total mean average glandular dose (AGD(T)) was calculated by summing the values associated with the pre-exposure and with the main exposure. Mean AGD(T) per exposure was 1.88 mGy (CI 0.01) and the mean AGD(T) per examination was 3.8 mGy, with 4 images per examination on average. The mean dose for cranio-caudal view (CC) images was 1.8 mGy, which is lower than that for medio-lateral oblique (MLO) view because the thickness for CC images was on average 10% lower than that for MLO images. Mean AGD(T) for the oldest group of women (1.90) was 3% higher than the AGD(T) for the younger group (1.85) due to the larger compressed breast thickness of women in the older group (10% on average). Differences between the corresponding AGD(T) values of each age group were lowest for breast thicknesses in the range 40-60 mm, being slightly higher for the women in the older group.     

A search for optimal x-ray spectra in iodine contrast media mammography

The aim of this work was to search for the optimal x-ray tube voltage and anode-filter combination in digital iodine contrast media mammography. In the optimization, two entities were of interest: the average glandular dose, AGD, and the signal-to-noise ratio, SNR, for detection of diluted iodine contrast medium. The optimum is defined as the technique maximizing the figure of merit, SNR2/AGD. A Monte Carlo computer program was used which simulates the transport of photons from the x-ray tube through the compression plate, breast, breast support plate, anti-scatter grid and image detector. It computes the AGD and the SNR of an iodine detail inside the compressed breast. The breast thickness was varied between 2 and 8 cm with 10-90% glandularity. The tube voltage was varied between 20 and 55 kV for each anode material (Rh, Mo and W) in combination with either 25 microm Rh or 0.05-0.5 mm Cu added filtration. The x-ray spectra were calculated with MCNP4C (Monte Carlo N-Particle Transport Code System, version 4C). A CsI scintillator was used as the image detector. The results for Rh/0.3 mmCu, Mo/0.3 mmCu and W/0.3 mmCu were similar. For all breast thicknesses, a maximum in the figure of merit was found at approximately 45 kV for the Rh/Cu, Mo/Cu and W/Cu combinations. The corresponding results for the Rh/Rh combination gave a figure of merit that was typically lower and more slowly varying with tube voltage. For a 4 cm breast at 45 kV, the SNR2/AGD was 3.5 times higher for the Rh/0.3 mmCu combination compared with the Rh/Rh combination. The difference is even larger for thicker breasts. The SNR2/AGD increases slowly with increasing Cu-filter thickness. We conclude that tube voltages between 41 and 55 kV and added Cu-filtration will result in significant dose advantage in digital iodine contrast media mammography compared to using the Rh/Rh anode/filter combination at 25-32 kV.     

Experimental investigation of the dose and image quality characteristics of a digital mammography imaging system

Our purpose in this study was to investigate the image quality and absorbed dose characteristics of a digital mammography imaging system with a CsI scintillator, and to identify an optimal x-ray tube voltage for imaging simulated masses in an average size breast with 50% glandularity. Images were taken of an ACR accreditation phantom using a LORAD digital mammography system with a Mo target and a Mo filter. In one experiment, exposures were performed at 80 mAs with x-ray tube voltages varying between 24 and 34 kVp. In a second experiment, the x-ray tube voltage was kept constant at 28 kVp and the technique factor was varied between 5 and 500 mAs. The average glandular dose at each x-ray tube voltage was determined from measurements of entrance skin exposure and x-ray beam half-value layer. Image contrast was measured as the fractional digital signal intensity difference for the image of a 4 mm thick acrylic disk. Image noise was obtained from the standard deviation in a uniformly exposed region of interest expressed as a fraction of the background intensity. The measured digital signal intensity was proportional to the mAs and to the kVp5.8. Image contrast was independent of mAs, and dropped by 21% when the x-ray tube voltage increased from 24 to 34 kVp. At a constant x-ray tube voltage, image noise was shown to be approximately proportional to (mAs)(-05), which permits the image contrast to noise ratio (CNR) to be modified by changing the mAs. At 80 mAs, increasing the x-ray tube voltage from 24 to 34 kVp increased the CNR by 78%, and increased the average glandular dose by 285%. At a constant lesion CNR, the lowest average glandular dose value occurred at 27.3 kVp. Increasing or decreasing the x-ray tube voltage by 2.3 kVp from the optimum kVp increased the average glandular dose values by 5%. These results show that imaging simulated masses in a 4.2 cm compressed breast at approximately 27 kVp with a Mo/Mo target/filter results in the lowest average glandular dose.     

Ambient dose equivalent and effective dose from scattered x-ray spectra in mammography for Mo/Mo, Mo/Rh and W/Rh anode/filter combinations

In this study, scattered x-ray distributions were produced by irradiating a tissue equivalent phantom under clinical mammographic conditions by using Mo/Mo, Mo/Rh and W/Rh anode/filter combinations, for 25 and 30 kV tube voltages. Energy spectra of the scattered x-rays have been measured with a Cd(0.9)Zn(0.1)Te (CZT) detector for scattering angles between 30 degrees and 165 degrees . Measurement and correction processes have been evaluated through the comparison between the values of the half-value layer (HVL) and air kerma calculated from the corrected spectra and measured with an ionization chamber in a nonclinical x-ray system with a W/Mo anode/filter combination. The shape of the corrected x-ray spectra measured in the nonclinical system was also compared with those calculated using semi-empirical models published in the literature. Scattered x-ray spectra measured in the clinical x-ray system have been characterized through the calculation of HVL and mean photon energy. Values of the air kerma, ambient dose equivalent and effective dose have been evaluated through the corrected x-ray spectra. Mean conversion coefficients relating the air kerma to the ambient dose equivalent and to the effective dose from the scattered beams for Mo/Mo, Mo/Rh and W/Rh anode/filter combinations were also evaluated. Results show that for the scattered radiation beams the ambient dose equivalent provides an overestimate of the effective dose by a factor of about 5 in the mammography energy range. These results can be used in the control of the dose limits around a clinical unit and in the calculation of more realistic protective shielding barriers in mammography.     

Calculation of dose and contrast for two mammographic grids.

To aid in selecting optimal conditions for screening mammography practice in Sweden, the performance of two mammographic grids (one moving and one stationary) has been investigated. Monte Carlo techniques were used to simulate photon transport in the breast. Transport through the breast support, grid covers, grid and image receptor (33.9 mg cm-2 Gd2O2S) was treated analytically. The contrast of a 100 microns calcification has been evaluated for three tissue compositions (adipose, glandular, 50:50 fractions by weight of adipose and glandular tissue) as a function of breast thickness (2-8 cm) and potential difference (25-30 kV, Mo anode). Contrast for a 5 cm 'average' breast at 28 kV was improved by 40% using the moving grid and by 30% using the stationary one; the corresponding increases in breast absorbed dose, keeping the energy imparted to the image receptor constant, were 90% and 150%, respectively. The superior properties of the moving grid were due to (i) equal scatter rejection ability and higher transmission of primary photons yielding lower scatter-to-primary ratios behind the grid, and (ii) less attenuation and filtering of the primary photons in the interspace material yielding lower degradation of primary contrast.     

A computer simulation study comparing lesion detection accuracy with digital mammography, breast tomosynthesis, and cone-beam CT breast imaging

Although conventional mammography is currently the best modality to detect early breast cancer, it is limited in that the recorded image represents the superposition of a three-dimensional (3D) object onto a 2D plane. Recently, two promising approaches for 3D volumetric breast imaging have been proposed, breast tomosynthesis (BT) and CT breast imaging (CTBI). To investigate possible improvements in lesion detection accuracy with either breast tomosynthesis or CT breast imaging as compared to digital mammography (DM), a computer simulation study was conducted using simulated lesions embedded into a structured 3D breast model. The computer simulation realistically modeled x-ray transport through a breast model, as well as the signal and noise propagation through a CsI based flat-panel imager. Polyenergetic x-ray spectra of Mo/Mo 28 kVp for digital mammography, Mo/Rh 28 kVp for BT, and W/Ce 50 kVp for CTBI were modeled. For the CTBI simulation, the intensity of the x-ray spectra for each projection view was determined so as to provide a total average glandular dose of 4 mGy, which is approximately equivalent to that given in conventional two-view screening mammography. The same total dose was modeled for both the DM and BT simulations. Irregular lesions were simulated by using a stochastic growth algorithm providing lesions with an effective diameter of 5 mm. Breast tissue was simulated by generating an ensemble of backgrounds with a power law spectrum, with the composition of 50% fibroglandular and 50% adipose tissue. To evaluate lesion detection accuracy, a receiver operating characteristic (ROC) study was performed with five observers reading an ensemble of images for each case. The average area under the ROC curves (Az) was 0.76 for DM, 0.93 for BT, and 0.94 for CTBI. Results indicated that for the same dose, a 5 mm lesion embedded in a structured breast phantom was detected by the two volumetric breast imaging systems, BT and CTBI, with statistically significant higher confidence than with planar digital mammography, while the difference in lesion detection between BT and CTBI was not statistically significant.     

X-ray spectroscopy applied to radiation shielding calculation in mammography

The protective shielding design of a mammography facility requires the knowledge of the scattered radiation by the patient and image receptor components. The shape and intensity of secondary x-ray beams depend on the kVp applied to the x-ray tube, target/filter combination, primary x-ray field size, and scattering angle. Currently, shielding calculations for mammography facilities are performed based on scatter fraction data for Mo/Mo target/filter, even though modern mammography equipment is designed with different anode/filter combinations. In this work we present scatter fraction data evaluated based on the x-ray spectra produced by a Mo/Mo, Mo/Rh and W/Rh target/filter, for 25, 30 and 35 kV tube voltages and scattering angles between 30 and 165 degrees. Three mammography phantoms were irradiated and the scattered radiation was measured with a CdZnTe detector. The primary x-ray spectra were computed with a semiempirical model based on the air kerma and HVL measured with an ionization chamber. The results point out that the scatter fraction values are higher for W/Rh than for Mo/Mo and Mo/Rh, although the primary and scattered air kerma are lower for W/Rh than for Mo/Mo and Mo/Rh target/filter combinations. The scatter fractions computed in this work were applied in a shielding design calculation in order to evaluate shielding requirements for each of these target/filter combinations. Besides, shielding requirements have been evaluated converting the scattered air kerma from mGy/week to mSv/week adopting initially a conversion coefficient from air kerma to effective dose as 1 Sv/Gy and then a mean conversion coefficient specific for the x-ray beam considered. Results show that the thickest barrier should be provided for Mo/Mo target/filter combination. They also point out that the use of the conversion coefficient from air kerma to effective dose as 1 Sv/Gy is conservatively high in the mammography energy range and overestimate the barrier thickness.     

Simulation study of a quasi-monochromatic beam for x-ray computed mammotomography

The purpose of this simulation study was to evaluate the feasibility, benefits, and potential operating parameters of a quasi-monochromatic beam from a tungsten-target x-ray source yielding projection images. The application is intended for newly developed cone beam computed mammotomography (CmT) of an uncompressed breast. The value of a near monochromatic x-ray source for a fully 3D CmT application is the expected improved ability to separate tissues with very small differences in attenuation coefficients. The quasi-monochromatic beam is expected to yield enhanced tomographic image quality along with a low dose, equal to or less than that of dual view x-ray mammography. X-ray spectra were generated with a validated projection x-ray simulation tool (XSpect) for a range of tungsten tube potentials (40-100 kVp), filter materials (Z=51-65), and filter thicknesses (10th to 1000th value layer determined at 60 kVp). The breast was modeled from ICRU-44 breast tissue specifications, and a breast lesion was modeled as a 0.5 cm thick mass. The detector was modeled as a digital flat-panel detector with a 0.06 cm thick CsI x-ray absorption layer. Computed figures of merit (FOMs) included the ratio of mean beam energy post-breast to pre-breast and the ratio of lesion contrasts for edge-located and center-located lesions as indices of breast beam hardening, and SNR2/exposure and SNR2/dose as indices of exposure and dose efficiencies. The impact of optimization of these FOMs on lesion contrast is also examined. For all simulated filter materials at each given attenuation thickness [10th, 100th, 500th, 1000th value layers (VLs)], the mean and standard deviation of the pre-breast spectral full-width at tenth-maximum (FWTM) were 16.1 +/- 2.4, 10.3 +/- 2.2, 7.3 +/- 1.4, and 6.5 +/- 1.5 keV, respectively. The change in beam width at the tenth maximum from pre-breast to post-breast spectra ranged from 4.7 to 1.1 keV, for the thinnest and thickest filters, respectively. The higher Z filters (Z=57-63) produced a quasi-monochromatic beam that allowed the widest tube potential operating range (50-70 kVp) while maintaining minimal beam hardening and maximal SNR2/exposure and SNR2/dose, and providing a contrast greater than that obtained in the unfiltered case. Figures of merit improved with increasing filter thickness, with diminishing returns beyond the 500th value layer attenuation level. Operating parameters required to produce optimal spectra, while keeping exposures equal to that of dual view mammography, are within the capability of the commercial x-ray tube proposed for our experimental study, indicating that use of these highly attenuating filters is viable. Additional simulations comparing Mo/Mo, Mo/Rh, and W/Rh target/filter combinations indicate that they exhibit significantly lower SNR2/exposure than the present approach, precluding them from being used for computed mammotomography, while maintaining dose limitations and obtaining sufficient SNR. Beam hardening was also much higher in the existing techniques (17%-42%) than for our technique (2%). Simulations demonstrate that this quasi-monochromatic x-ray technique may enhance tissue separation for a newly developed cone beam computed mammotomography application for an uncompressed breast.     

A comprehensive analysis of DgN(CT) coefficients for pendant-geometry cone-beam breast computed tomography

The use of a computed tomography (CT) scanner specifically designed for breast imaging has been proposed by several investigators. In this study, the radiation dose due to breast CT was evaluated using Monte Carlo techniques over a range of parameters pertinent to the cone-beam pendant geometry thought to be most appropriate. Monte Carlo dose computations were validated by comparison with physical measurements made on a prototype breast CT scanner under development in our laboratory. The Monte Carlo results were then used to study the influence of cone angle, the use of a beam flattening ("bow-tie") filter, glandular fraction, breast length and source-to-isocenter distance. These parameters were studied over a range of breast diameters from 10 to 18 cm, and for both monoenergetic (8-140 keV by 1 keV intervals) and polyenergetic x-ray beams (30-100 kVp by 5 kVp intervals. Half value layer at 80 kVp = 5.3 mm Al). A parameter referring to the normalized glandular dose in CT (DgN(CT)) was defined which is the ratio of the glandular dose in the breast to the air kerma at isocenter. There was no significant difference (p = 0.743) between physically measured and Monte Carlo derived results. Fan angle, source-to-isocenter distance, and breast length have relatively small influences on the radiation dose in breast CT. Glandular fraction (0% versus 100%) for 10 cm breasts at 80 kVp had approximately a 10% effect on DgN(CT), and a 20% effect was observed for an 18 cm breast diameter. The use of a bow-tie filter had the potential to reduce breast dose by approximately 40%. X-ray beam energy and breast diameter had significant influence on the DgN(CT) parameters, with higher DgN(CT) values for higher energy beams and smaller breast diameters. DgN(CT) values (mGy/mGy) at 80 kVp ranged from 0.95 for an 8 cm diam 50% glandular breast to 0.78 for an 18 cm 50% glandular breast. The results of this investigation should be useful for those interested computing the glandular breast dose for geometries relevant to dedicated breast CT.     

Technique factors and their relationship to radiation dose in pendant geometry breast CT

The use of breast computed tomography (CT) as an alternative to mammography in some patients is being studied at several institutions. However, the radiation dosimetry issues associated with breast CT are markedly different than in the case of mammography. In this study, the spectral properties of an operational breast CT scanner were characterized both by physical measurement and computer modeling of the kVp-dependent spectra, from 40 to 110 kVp (Be window W anode with 0.30 mm added Cu filtration). Previously reported conversion factors, normalized glandular dose for CT-DgN(ct), derived from Monte Carlo methods, were used in concert with the output spectra of the breast scanner to compute the mean glandular dose to the breast based upon different combinations of x-ray technique factors (kVp and mAs). The mean glandular dose (MGD) was measured as a function of the compressed breast thickness (2-8 cm) and three different breast compositions (0%, 50%, and 100% glandular fractions) in four clinical mammography systems in our institution. The average MGD from these four systems was used to compute the technique factors for breast CT systems that would match the two-view mammographic dose levels. For a 14 cm diameter breast (equivalent to a 5 cm thick compressed breast in mammography), air kerma levels at the breast CT scanner's isocenter (468 mm from the source) of 4.4, 6.4, and 9.0 mGy were found to deliver equivalent mammography doses for 0%, 50%, and 100% glandular breasts (respectively) at 80 kVp. At 80 kVp (where air kerma was 11.3 mGy/100 mAs at the isocenter), 57 mAs (integrated over the entire scan) was required to match the mammography dose for a 14 cm 50% glandular breast. At 50 kVp, 360 mAs is required to match mammographic dose levels. Tables are provided for both air kerma at the isocenter and mAs for 0%, 50%, and 100% glandular breasts. Other issues that impact breast CT technique factors are also discussed.     

A survey of clinical factors and patient dose in mammography

A survey was conducted to estimate the mean glandular dose (MGD) for women undergoing mammography and to report the distribution of doses, compressed breast thickness, glandular tissue content, and mammographic technique factors used. From 24,471 mammograms, of 6,006 women, clinical data were collected. The survey data included mammograms from seven modern units using a molybdenum (Mo) anode and either Mo or rhodium (Rh) filter. Exposure factors for each mammogram were recorded automatically onto a floppy disk on each unit. All mammography units were calibrated individually using breast tissue equivalent attenuation slabs of varying glandular content, so the breast glandular content could be estimated on the basis of exposure factors and compressed breast thickness. The MGD was estimated for each mammogram based on the normalized glandular dose and calculated entrance exposure in air. The survey found a median MGD of 2.6 mGy. The median breast glandular tissue content was 28% and the median compressed breast thickness was 5.1 cm. Also, patient attenuation data were converted to equivalent BR-12 and acrylic thickness to help determine appropriate phantom thicknesses required for mammography unit automatic exposure control performance assessment.     

Comparison of three tissue composition measurement techniques using digital mammograms — A signal-to-noise study

Tissue composition measurement may provide a quantitatively and physically meaningful method to objectively determine the “mammographic density” linked to breast cancer risk. A single energy hybrid (SEH) technique is described for measuring the tissue composition on a pixel-by-pixel basis with a single digital mammogram. Theoretical models were derived and used to compute signal-to-noise ratios (SNRs) in tissue composition measurement using the SEH method. The results were compared with those computed for measurements using the dual kVp and dual screen methods. SNRs were theoretically related to the pixel area, total unattenuated detector exposure and fluence spectra of the incident X-rays. SNRs were computed for measurement of the adipose tissue thickness for a 6 cm thick breast, consisting of 50% of adipose tissue and 50% of glandular tissue. Effects of kVp and prepatient filtration were studied by computing the SNRs for various kVps and filters and optimal kVps and filters are determined. The results showed that the SNRs obtained with the SEH method is an order of magnitude better than the dual kVp method, which, in turn, is an order of magnitude better than the dual screen method. When using the optimal kVp’s and no prepatient filtration, the SNRs were computed to be 84.2, 13.2, and 2.0 for the SEH, dual kVp, and dual screen methods, respectively. Prepatient filtration can improve the SNR by as much as 35% for the dual kVp and dual screen techniques with reasonable tube loading factors (8–10).     

High-resolution imager for digital mammography: physical characterization of a prototype sensor

The physical performance characteristics of a high-resolution sensor module for digital mammography were investigated. The signal response of the imager was measured at various detector entrance air kerma and was found to be linear. The spatial resolution was determined by measuring the presampling modulation transfer function, MTF(f), of the system. The noise power spectra, NPS(f), of the system were estimated using 26 kVp: Mo/Mo, 28 kVp: Mo/Rh and 30 kVp: Rh/Rh, with polymethyl methacrylate (PMMA) 'tissue equivalent material' of thickness 20, 45 and 57 mm for each of three x-ray spectra at detector entrance air kerma in the range between approximately 80.2 and 92.3 microGy. The noise equivalent quanta, NEQ(f), and detective quantum efficiencies, DQE(f), for the various spectral conditions were computed. In addition, dose dependence of NPS(f) and DQE(f) was studied at various detector entrance air kerma ranging from 9.4 to 169.7 microGy. A spatial resolution of about 10 cycles mm(-1) was obtained at the 10% MTF(f) level. A small increase in NEQ(f)was observed under higher energy spectral conditions while the DQE(f) decreased marginally. For a given spectrum, increasing PMMA filtration produced negligible change in DQE(f). The estimated DQE values at zero frequency were in the range between 0.45 and 0.55 under the conditions investigated in this study.     

Effect of area x-ray beam equalization on image quality and dose in digital mammography

In mammography, thick or dense breast regions persistently suffer from reduced contrast-to-noise ratio (CNR) because of degraded contrast from large scatter intensities and relatively high noise. Area x-ray beam equalization can improve image quality by increasing the x-ray exposure to under-penetrated regions without increasing the exposure to other breast regions. Optimal equalization parameters with respect to image quality and patient dose were determined through computer simulations and validated with experimental observations on a step phantom and an anthropomorphic breast phantom. Three parameters important in equalization digital mammography were considered: attenuator material (Z = 13-92), beam energy (22-34 kVp) and equalization level. A Mo/Mo digital mammography system was used for image acquisition. A prototype 16 x 16 piston driven equalization system was used for preparing patient-specific equalization masks. Simulation studies showed that a molybdenum attenuator and an equalization level of 20 were optimal for improving contrast, CNR and figure of merit (FOM = CNR2/dose). Experimental measurements using these parameters showed significant improvements in contrast, CNR and FOM. Moreover, equalized images of a breast phantom showed improved image quality. These results indicate that area beam equalization can improve image quality in digital mammography.     

Optimization of technique factors for a silicon diode array full-field digital mammography system and comparison to screen-film mammography with matched average glandular dose

Contrast-detail experiments were performed to optimize technique factors for the detection of low-contrast lesions using a silicon diode array full-field digital mammography (FFDM) system under the conditions of a matched average glandular dose (AGD) for different techniques. Optimization was performed for compressed breast thickness from 2 to 8 cm. FFDM results were compared to screen-film mammography (SFM) at each breast thickness. Four contrast-detail (CD) images were acquired on a SFM unit with optimal techniques at 2, 4, 6, and 8 cm breast thicknesses. The AGD for each breast thickness was calculated based on half-value layer (HVL) and entrance exposure measurements on the SFM unit. A computer algorithm was developed and used to determine FFDM beam current (mAs) that matched AGD between FFDM and SFM at each thickness, while varying target, filter, and peak kilovoltage (kVp) across the full range available for the FFDM unit. CD images were then acquired on FFDM for kVp values from 23-35 for a molybdenum-molybdenum (Mo-Mo), 23-40 for a molybdenum-rhodium (Mo-Rh), and 25-49 for a rhodium-rhodium (Rh-Rh) target-filter under the constraint of matching the AGD from screen-film for each breast thickness (2, 4, 6, and 8 cm). CD images were scored independently for SFM and each FFDM technique by six readers. CD scores were analyzed to assess trends as a function of target-filter and kVp and were compared to SFM at each breast thickness. For 2 cm thick breasts, optimal FFDM CD scores occurred at the lowest possible kVp setting for each target-filter, with significant decreases in FFDM CD scores as kVp was increased under the constraint of matched AGD. For 2 cm breasts, optimal FFDM CD scores were not significantly different from SFM CD scores. For 4-8 cm breasts, optimum FFDM CD scores were superior to SFM CD scores. For 4 cm breasts, FFDM CD scores decreased as kVp increased for each target-filter combination. For 6 cm breasts, CD scores decreased slightly as kVp increased for Mo-Mo, but did not change significantly as a function of kVp for either Mo-Rh or Rh-Rh. For 8 cm breasts, Rh/Rh FFDM CD scores were superior to other target-filter combinations and increased significantly as kVp increased. These results indicate that low-contrast lesion detection was optimized for FFDM by using a softer x-ray beam for thin breasts and a harder x-ray beam for thick breasts, when AGD was kept constant for a given breast thickness. Under this constraint, optimum low-contrast lesion detection with FFDM was superior to that for SFM for all but the thinnest breasts.     

Resolution at oblique incidence angles of a flat panel imager for breast tomosynthesis

Oblique incidence of x rays on an imaging detector causes blurring that reduces spatial resolution. For simple projection imaging this effect is small and often ignored. However, for breast tomosynthesis, the incidence angle can be larger (>20 degrees), leading to increased blur for some of the projections. The modulation transfer function (MTF) is measured for a typical phosphor-coupled flat-panel detector versus angular incidence of the x-ray beam for two x-ray spectra: 26 kV Mo/Mo and 40 kV Rh/Al. At an incidence angle of 40 degrees the MTF at 5 mm(-1) falls by 35% and 40% for each spectrum, respectively (and 65%/80% at 8 mm(-1)). Increasing the detector absorber thickness to achieve improved quantum efficiency will cause the blurring effect due to beam obliquity to become greater. The impact of this blur is likely to cause misregistration and increased relative noise in tomosynthesis reconstructed images.     

Scatter/primary in mammography: comprehensive results

Monte Carlo procedures using the SIERRA code (validated in a companion article) were used to investigate the scatter properties in mammography. The scatter to primary ratio (SPR) was used for quantifying scatter levels as a function of beam spectrum, position in the field, air gap, breast thickness, tissue composition, and the area of the field of view (FOV). The geometry of slot scan mammography was also simulated, and SPR values were calculated as a function of slot width. The influence of large air gaps (to 30 cm) was also studied in the context of magnification mammography. X-ray energy and tissue composition from 100% adipose to 100% glandular demonstrated little effect on the SPR. Air gaps over a range from 0 to 30 mm showed only slight effects. The SPR increased with increased breast thickness and with larger fields of view. Measurements from 82 mammograms provided estimates of the range of compressed breast thickness (median: 5.2 cm, 95% range: 2.4 cm to 7.9 cm) and projected breast area onto the film (left craniocaudal view, median: 146 cm2, 95% range: 58 cm2 to 298 cm2). SPR values for semicircular breast shapes, Mo/Mo spectra, and a 15 mm air gap were parametrized as a function of breast thickness and (semicircular) breast diameter. With the coefficients a = - 2.35452817439093, b = 22.3960980055927, and c = 8.85064260299289, the equation SPR= [a + b x (diameter in cm)--(-1.5) + c x (thickness in cm) --(-0.5)]-- -1 produces SPR data from 2 to 8 cm and from 3 to 30 cm breast diameters with an average error of about 1%.     

Computation of beam quality parameters for Mo/Mo, Mo/Rh, Rh/Rh, and W/Al target/filter combinations in mammography

A computer program was implemented to predict mammography x-ray beam parameters in the range 20-40 kV for Mo/Mo, Mo/Rh, Rh/Rh, and W/Al target/filter combinations. The computation method used to simulate mammography x-ray spectra is based on the Boone et al. model. The beam quality parameters such as the half-value layer (HVL), the homogeneity coefficient (HC), and the average photon energy were computed by simulating the interaction of the spectrum photons with matter. The checking of this computation was done using a comparison of the results with published data and measured values obtained at the Netherlands Metrology Institute Van Swinden Laboratorium, National Institute of Standards and Technology, and International Atomic Energy Agency. The predicted values with a mean deviation of 3.3% of HVL, 3.7% of HC, and 1.5% of average photon energy show acceptable agreement with published data and measurements for all target/filter combinations in the 23-40 kV range. The accuracy of this computation can be considered clinically acceptable and can allow an appreciable estimation for the beam quality parameters.