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.     

Computation of the glandular radiation dose in digital tomosynthesis of the breast

Tomosynthesis of the breast is currently a topic of intense interest as a logical next step in the evolution of digital mammography. This study reports on the computation of glandular radiation dose in digital tomosynthesis of the breast. Previously, glandular dose estimations in tomosynthesis have been performed using data from studies of radiation dose in conventional planar mammography. This study evaluates, using Monte Carlo methods, the normalized glandular dose (DgN) to the breast during a tomosynthesis study, and characterizes its dependence on breast size, tissue composition, and x-ray spectrum. The conditions during digital tomosynthesis imaging of the breast were simulated using a computer program based on the Geant4 toolkit. With the use of simulated breasts of varying size, thickness and tissue composition, the DgN to the breast tissue was computed for varying x-ray spectra and tomosynthesis projection angle. Tomosynthesis projections centered about both the cranio-caudal (CC) and medio-lateral oblique (MLO) views were simulated. For each projection angle, the ratio of the glandular dose for that projection to the glandular dose for the zero degree projection was computed. This ratio was denoted the relative glandular dose (RGD) coefficient, and its variation under different imaging parameters was analyzed. Within mammographic energies, the RGD was found to have a weak dependence on glandular fraction and x-ray spectrum for both views. A substantial dependence on breast size and thickness was found for the MLO view, and to a lesser extent for the CC view. Although RGD values deviate substantially from unity as a function of projection angle, the RGD averaged over all projections in a complete tomosynthesis study varies from 0.91 to 1.01. The RGD results were fit to mathematical functions and the resulting equations are provided.     

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.     

Dual-energy contrast-enhanced breast tomosynthesis: optimization of beam quality for dose and image quality

Dual-energy contrast-enhanced breast tomosynthesis is a promising technique to obtain three-dimensional functional information from the breast with high resolution and speed. To optimize this new method, this study searched for the beam quality that maximized image quality in terms of mass detection performance. A digital tomosynthesis system was modeled using a fast ray-tracing algorithm, which created simulated projection images by tracking photons through a voxelized anatomical breast phantom containing iodinated lesions. The single-energy images were combined into dual-energy images through a weighted log subtraction process. The weighting factor was optimized to minimize anatomical noise, while the dose distribution was chosen to minimize quantum noise. The dual-energy images were analyzed for the signal difference to noise ratio (SdNR) of iodinated masses. The fast ray-tracing explored 523?776 dual-energy combinations to identify which yields optimum mass SdNR. The ray-tracing results were verified using a Monte Carlo model for a breast tomosynthesis system with a selenium-based flat-panel detector. The projection images from our voxelized breast phantom were obtained at a constant total glandular dose. The projections were combined using weighted log subtraction and reconstructed using commercial reconstruction software. The lesion SdNR was measured in the central reconstructed slice. The SdNR performance varied markedly across the kVp and filtration space. Ray-tracing results indicated that the mass SdNR was maximized with a high-energy tungsten beam at 49 kVp with 92.5 μm of copper filtration and a low-energy tungsten beam at 49 kVp with 95 μm of tin filtration. This result was consistent with Monte Carlo findings. This mammographic technique led to a mass SdNR of 0.92 ± 0.03 in the projections and 3.68 ± 0.19 in the reconstructed slices. These values were markedly higher than those for non-optimized techniques. Our findings indicate that dual-energy breast tomosynthesis can be performed optimally at 49 kVp with alternative copper and tin filters, with reconstruction following weighted subtraction. The optimum technique provides best visibility of iodine against structured breast background in dual-energy contrast-enhanced breast tomosynthesis.     

Estimation of mean glandular dose for breast tomosynthesis: factors for use with the UK, European and IAEA breast dosimetry protocols

A formalism is proposed for the estimation of mean glandular dose for breast tomosynthesis, which is a simple extension of the UK, European and IAEA protocols for dosimetry in conventional projection mammography. The formalism introduces t-factors for the calculation of breast dose from a single projection and T-factors for a complete exposure series. Monte Carlo calculations of t-factors have been made for an imaging geometry with full-field irradiation of the breast for a wide range of x-ray spectra, breast sizes and glandularities. The t-factors show little dependence on breast glandularity and tables are provided as a function of projection angle and breast thickness, which may be used for all x-ray spectra simulated. The T-factors for this geometry depend upon the choice of projection angles and weights per projection, but various example calculations gave values in the range 0.93-1.00. T-factors are also provided for the Sectra tomosynthesis system, which employs a scanned narrow-beam imaging geometry. In this quite different configuration, the factor (denoted T(S)) shows an important dependence on breast thickness, varying between 0.98 and 0.76 for 20 and 110 mm thick breasts, respectively. Additional data are given to extend the current tabulations of g-, c- and s-factors used for dosimetry of conventional 2D mammography.     

Scatter radiation in digital tomosynthesis of the breast

Digital tomosynthesis of the breast is being investigated as one possible solution to the problem of tissue superposition present in planar mammography. This imaging technique presents various advantages that would make it a feasible replacement for planar mammography, among them similar, if not lower, radiation glandular dose to the breast; implementation on conventional digital mammography technology via relatively simple modifications; and fast acquisition time. One significant problem that tomosynthesis of the breast must overcome, however, is the reduction of x-ray scatter inclusion in the projection images. In tomosynthesis, due to the projection geometry and radiation dose considerations, the use of an antiscatter grid presents several challenges. Therefore, the use of postacquisition software-based scatter reduction algorithms seems well justified, requiring a comprehensive evaluation of x-ray scatter content in the tomosynthesis projections. This study aims to gain insight into the behavior of x-ray scatter in tomosynthesis by characterizing the scatter point spread functions (PSFs) and the scatter to primary ratio (SPR) maps found in tomosynthesis of the breast. This characterization was performed using Monte Carlo simulations, based on the Geant4 toolkit, that simulate the conditions present in a digital tomosynthesis system, including the simulation of the compressed breast in both the cranio-caudal (CC) and the medio-lateral oblique (MLO) views. The variation of the scatter PSF with varying tomosynthesis projection angle, as well as the effects of varying breast glandular fraction and x-ray spectrum, was analyzed. The behavior of the SPR for different projection angle, breast size, thickness, glandular fraction, and x-ray spectrum was also analyzed, and computer fit equations for the magnitude of the SPR at the center of mass for both the CC and the MLO views were found. Within mammographic energies, the x-ray spectrum was found to have no appreciable effect on the scatter PSF and on the SPR. Glandular fraction and compressed breast size were found to have a small effect, while compressed breast thickness and projection angle, as expected, introduced large variations in both the scatter PSF and SPR. The presence of the breast support plate and the detector cover plate in the simulations introduced important effects on the SPR, which are also relevant to the scatter content in planar mammography.     

Imaging performance of an amorphous selenium digital mammography detector in a breast tomosynthesis system

In breast tomosynthesis a rapid sequence of N images is acquired when the x-ray tube sweeps through different angular views with respect to the breast. Since the total dose to the breast is kept the same as that in regular mammography, the exposure used for each image of tomosynthesis is 1/N. The low dose and high frame rate pose a tremendous challenge to the imaging performance of digital mammography detectors. The purpose of the present work is to investigate the detector performance in different operational modes designed for tomosynthesis acquisition, e.g., binning or full resolution readout, the range of view angles, and the number of views N. A prototype breast tomosynthesis system with a nominal angular range of +/-25 degrees was used in our investigation. The system was equipped with an amorphous selenium (a-Se) full field digital mammography detector with pixel size of 85 microm. The detector can be read out in full resolution or 2 x 1 binning (binning in the tube travel direction). The focal spot blur due to continuous tube travel was measured for different acquisition geometries, and it was found that pixel binning, instead of focal spot blur, dominates the detector modulation transfer function (MTF). The noise power spectrum (NPS) and detective quantum efficiency (DQE) of the detector were measured with the exposure range of 0.4-6 mR, which is relevant to the low dose used in tomosynthesis. It was found that DQE at 0.4 mR is only 20% less than that at highest exposure for both detector readout modes. The detector temporal performance was categorized as lag and ghosting, both of which were measured as a function of x-ray exposure. The first frame lags were 8% and 4%, respectively, for binning and full resolution mode. Ghosting is negligible and independent of the frame rate. The results showed that the detector performance is x-ray quantum noise limited at the low exposures used in each view of tomosynthesis, and the temporal performance at high frame rate (up to 2 frames per second) is adequate for tomosynthesis.     

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.     

Cone-beam volume CT breast imaging: feasibility study

X-ray projection mammography, using a film/screen combination, or digital techniques, has proven to be the most effective imaging modality currently available for early detection of breast cancer. However, the inherent superimposition of structures makes a small carcinoma (a few millimeters in size) difficult to detect when it is occult or in dense breasts, leading to a high false-positive biopsy rate. Cone-beam x-ray-projection-based volume imaging using flat panel detectors (FPDs) may allow obtaining three-dimensional breast images, resulting in more accurate diagnosis of structures and patterns of lesions while eliminating the hard compression of breasts. This article presents a novel cone-beam volume computed tomographic breast imaging (CBVCTBI) technique based on the above techniques. Through a variety of computer simulations, the key issues of the system and imaging techniques were addressed, including the x-ray imaging geometry and corresponding reconstruction algorithms, x-ray characteristics of breast tissue and lesions, x-ray setting techniques, the absorbed dose estimation, and the quantitative effect of x-ray scattering on image quality. The preliminary simulation results support the proposed CVBCTBI modality for breast imaging in respect to its feasibility and practicability. The absorbed dose level is comparable to that of current mammography and will not be a prominent problem for this imaging technique. Compared to conventional mammography, the proposed imaging technique with isotropic spatial resolution will potentially provide significantly better low-contrast detectability of breast tumors and more accurate location of breast lesions.     

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.     

Normalized glandular dose (DgN) coefficients for flat-panel CT breast imaging

The development of new digital mammography techniques such as dual-energy imaging, tomosynthesis and CT breast imaging will require investigation of optimal camera design parameters and optimal imaging acquisition parameters. In optimizing these acquisition protocols and imaging systems it is important to have knowledge of the radiation dose to the breast. This study presents a methodology for estimating the normalized glandular dose to the uncompressed breast using the geometry proposed for flat-panel CT breast imaging. The simulation uses the GEANT 3 Monte Carlo code to model x-ray transport and absorption within the breast phantom. The Monte Carlo software was validated for breast dosimetry by comparing results of the normalized glandular dose (DgN) values of the compressed breast to those reported in the literature. The normalized glandular dose was then estimated for a range of breast diameters from 10 cm to 18 cm using an uncompressed breast model with a homogeneous composition of adipose and glandular tissue, and for monoenergetic x-rays from 10 keV to 120 keV. These data were fit providing expressions for the normalized glandular dose. Using these expressions for the DgN coefficients and input variables such as the diameter, height and composition of the breast phantom, the mean glandular dose for any spectra can be estimated. A computer program to provide normalized glandular dose values has been made available online. In addition, figures displaying energy deposition maps are presented to better understand the spatial distribution of dose in CT breast imaging.     

Evaluating the impact of X-ray spectral shape on image quality in flat-panel CT breast imaging

In recent years, there has been an increasing interest in exploring the feasibility of dedicated computed tomography (CT) breast imaging using a flat-panel digital detector in a truncated cone-beam imaging geometry. Preliminary results are promising and it appears as if three-dimensional tomographic imaging of the breast has great potential for reducing the masking effect of superimposed parenchymal structure typically observed with conventional mammography. In this study, a mathematical framework used for determining optimal design and acquisition parameters for such a CT breast imaging system is described. The ideal observer signal-to-noise ratio (SNR) is used as a figure of merit, under the assumptions that the imaging system is linear and shift invariant. Computation of the ideal observer SNR used a parallel-cascade model to predict signal and noise propagation through the detector, as well as a realistic model of the lesion detection task in breast imaging. For all evaluations, the total mean glandular dose for a CT breast imaging study was constrained to be approximately equivalent to that of a two-view conventional mammography study. The framework presented was used to explore the effect of x-ray spectral shape across an extensive range of kVp settings, filter material types, and filter thicknesses. The results give an indication of how spectral shape can affect image quality in flat-panel CT breast imaging.     

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.     

Quantification of breast density with spectral mammography based on a scanned multi-slit photon-counting detector: a feasibility study

A simple and accurate measurement of breast density is crucial for the understanding of its impact in breast cancer risk models. The feasibility to quantify volumetric breast density with a photon-counting spectral mammography system has been investigated using both computer simulations and physical phantom studies. A computer simulation model involved polyenergetic spectra from a tungsten anode x-ray tube and a Si-based photon-counting detector has been evaluated for breast density quantification. The figure-of-merit (FOM), which was defined as the signal-to-noise ratio of the dual energy image with respect to the square root of mean glandular dose, was chosen to optimize the imaging protocols, in terms of tube voltage and splitting energy. A scanning multi-slit photon-counting spectral mammography system has been employed in the experimental study to quantitatively measure breast density using dual energy decomposition with glandular and adipose equivalent phantoms of uniform thickness. Four different phantom studies were designed to evaluate the accuracy of the technique, each of which addressed one specific variable in the phantom configurations, including thickness, density, area and shape. In addition to the standard calibration fitting function used for dual energy decomposition, a modified fitting function has been proposed, which brought the tube voltages used in the imaging tasks as the third variable in dual energy decomposition. For an average sized 4.5 cm thick breast, the FOM was maximized with a tube voltage of 46 kVp and a splitting energy of 24 keV. To be consistent with the tube voltage used in current clinical screening exam (～32 kVp), the optimal splitting energy was proposed to be 22 keV, which offered a FOM greater than 90% of the optimal value. In the experimental investigation, the root-mean-square (RMS) error in breast density quantification for all four phantom studies was estimated to be approximately 1.54% using standard calibration function. The results from the modified fitting function, which integrated the tube voltage as a variable in the calibration, indicated a RMS error of approximately 1.35% for all four studies. The results of the current study suggest that photon-counting spectral mammography systems may potentially be implemented for an accurate quantification of volumetric breast density, with an RMS error of less than 2%, using the proposed dual energy imaging technique.     

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.     

The effect of breast compression on mass conspicuity in digital mammography

This study analyzed how the inherent quality of diagnostic information in digital mammography could be affected by breast compression. A digital mammography system was modeled using a Monte Carlo algorithm based on the Penelope program, which has been successfully used to model several medical imaging systems. First, the Monte Carlo program was validated against previous measurements and simulations. Once validated, the Monte Carlo software modeled a digital mammography system by tracking photons through a voxelized software breast phantom, containing anatomical structures and breast masses, and following photons until they were absorbed by a selenium-based flat-panel detector. Simulations were performed for two compression conditions (standard compression and 12.5% reduced compression) and three photon flux conditions (constant flux, constant detector signal, and constant glandular dose). The results showed that reduced compression led to higher scatter fractions, as expected. For the constant photon flux condition, decreased compression also reduced glandular dose. For constant glandular dose, the SdNR for a 4 cm breast was 0.60 +/- 0.11 and 0.62 +/- 0.11 under standard and reduced compressions, respectively. For the 6 cm case with constant glandular dose, the SdNR was 0.50 +/- 0.11 and 0.49 +/- 0.10 under standard and reduced compressions, respectively. The results suggest that if a particular imaging system can handle an approximately 10% increase in total tube output and 10% decrease in detector signal, breast compression can be reduced by about 12% in terms of breast thickness with little impact on image quality or dose.     

Why should breast tumour detection go three dimensional?

Although x-ray mammography is widely developed for breast tumour detection, it suffers from spatial superposition in its two-dimensional (2D) representation of a three-dimensional (3D) breast structure. Accordingly, 3D breast imaging, such as cone-beam computed tomography (CT), arises at the historic moment. In this paper, we theoretically elucidate the spatial superposition effect associated with x-ray mammography on breast tumour detection. This explanation is based on the line integral of x-ray traversing a composite breast model. As a result, we can characterize the difficulty of detecting small tumours in terms of local intensity contrast in x-ray images. In comparison, we also introduce cone-beam CT breast imaging for 3D breast volume representation, which offers advantages for breast mass segmentation and measurement. The discussion is demonstrated with an experiment with a breast surgical specimen. In conclusion, we strongly believe that 3D volumetric representation allows for more accurate breast tumour detection.     

Optimized radiographic spectra for small animal digital subtraction angiography

The increasing use of small animals in basic research has spurred interest in new imaging methodologies. Digital subtraction angiography (DSA) offers a particularly appealing approach to functional imaging in the small animal. This study examines the optimal x-ray, molybdenum (Mo) or tungsten (W) target sources, and technique to produce the highest quality small animal functional subtraction angiograms in terms of contrast and signal-difference-to-noise ratio squared (SdNR2). Two limiting conditions were considered-normalization with respect to dose and normalization against tube loading. Image contrast and SdNR2 were simulated using an established x-ray model. DSA images of live rats were taken at two representative tube potentials for the W and Mo sources. Results show that for small animal DSA, the Mo source provides better contrast. However, with digital detectors, SdNR2 is the more relevant figure of merit. The W source operated at kVps >60 achieved a higher SdNR2. The highest SdNR2 was obtained at voltages above 90 kVp. However, operation at the higher potential results in significantly greater dose and tube load and reduced contrast quantization. A reasonable tradeoff can be achieved at tube potentials at the beginning of the performance plateau, around 70 kVp, where the relative gain in SdNR2 is the greatest.     

Development and validation of a simulation procedure to study the visibility of micro calcifications in digital mammograms

The visibility of micro calcifications is a determining factor for digital mammography. To address the problem of quantification, we developed a procedure to simulate micro calcifications into real mammograms. First, the shapes, sizes and x-ray transmission coefficients of real micro calcifications were derived from the appearance of biopsy specimens in the raw data of magnified, digital images acquired at 27 kVp and Mo/Mo anode-filter combination. This allowed us to create "ideal templates" of micro calcifications. The x-ray transmissions of the real micro calcifications values were expressed in Al-equivalent thickness. This made it possible to recalculate the x-ray transmission characteristics of a particular ideal template for other x-ray beam qualities. Extra corrections for differences in spatial resolution were based on the presampled two-dimensional modulation transfer functions and on the difference in pixel size. Three radiologists compared the appearance of real and simulated micro calcifications in a two-alternative forced choice (2AFC) evaluation. They perceived no differences between real and simulated lesions. Preliminary results show that it is possible to simulate micro calcifications with well defined characteristics that are indistinguishable from real ones. It should be noted, however, that the full potential of the approach has not been proven. In future work, these templates may be useful to evaluate particular aspects of digital mammograms, such as the effects of processing and of viewing conditions on the visibility of micro calcifications.