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.     

Anisotropic imaging performance in indirect x-ray imaging detectors

We report on the variability in imaging system performance due to oblique x-ray incidence, and the associated transport of quanta (both x rays and optical photons) through the phosphor, in columnar indirect digital detectors. The analysis uses MANTIS, a combined x-ray, electron, and optical Monte Carlo transport code freely available. We describe the main features of the simulation method and provide some validation of the phosphor screen models considered in this work. We report x-ray and electron three-dimensional energy deposition distributions and point-response functions (PRFs), including optical spread in columnar phosphor screens of thickness 100 and 500 microm, for 19, 39, 59, and 79 keV monoenergetic x-ray beams incident at 0 degrees, 10 degrees, and 15 degrees. In addition, we present pulse-height spectra for the same phosphor thickness, x-ray energies, and angles of incidence. Our results suggest that the PRF due to the phosphor blur is highly nonsymmetrical, and that the resolution properties of a columnar screen in a tomographic, or tomosynthetic imaging system varies significantly with the angle of x-ray incidence. Moreover, we find that the noise due to the variability in the number of light photons detected per primary x-ray interaction, summarized in the information or Swank factor, is somewhat independent of thickness and incidence angle of the x-ray beam. Our results also suggest that the anisotropy in the PRF is not less in screens with absorptive backings, while the noise introduced by variations in the gain and optical transport is larger. Predictions from MANTIS, after additional validation, can provide the needed understanding of the extent of such variations, and eventually, lead to the incorporation of the changes in imaging performance with incidence angle into the reconstruction algorithms for volumetric x-ray imaging systems.     

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.     

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.     

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.     

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.     

A computer simulation platform for the optimization of a breast tomosynthesis system

In breast tomosynthesis there is a compromise between resolution, noise, and acquisition speed for a given glandular dose. The purpose of the present work is to develop a simulation platform to investigate the potential imaging performance for the many possible tomosynthesis system configurations. The simulation platform was used to investigate the dependence of image blur and signal difference to noise ratio (SDNR) for several different tomosynthesis acquisition configurations. Simulated projections of a slanted thin tungsten wire placed in different object planes were modified according to the detector's modulation transfer function (MTF), with or without pixel binning. In addition, the focal spot blur (FSB), which depends on the location of the wire, the system geometry, the source-detector movement speed, and the exposure time, was also incorporated into the projections. Both expectation maximization (EM) and filtered back projection (FBP) based algorithms were used for 3D image reconstruction. The in-plane MTF was calculated from the reconstructed image of the wire. To evaluate the noise performance, simulated noiseless projections of calcification and tumor in uniform breast tissue were modified with the noise power spectrum (NPS) calculated from a cascaded linear system model for the detector for a given x-ray dose. The SDNR of the reconstructed images was calculated with different tomosynthesis configurations, e.g., pixel binning, view number, and angular range. Our results showed that for a source-to-imager distance (SID) of 66 cm, pixel binning (2 x 2) caused more degradation to the in-plane MTF than the blur caused by the moving focal spot and reconstruction. The in-depth resolution can be improved by increasing the angular range.     

Importance of point-by-point back projection correction for isocentric motion in digital breast tomosynthesis: relevance to morphology of structures such as microcalcifications

Digital breast tomosynthesis is a three-dimensional imaging technique that provides an arbitrary set of reconstruction planes in the breast from a limited-angle series of projection images acquired while the x-ray tube moves. Traditional shift-and-add (SAA) tomosynthesis reconstruction is a common mathematical method to line up each projection image based on its shifting amount to generate reconstruction slices. With parallel-path geometry of tube motion, the path of the tube lies in a plane parallel to the plane of the detector. The traditional SAA algorithm gives shift amounts for each projection image calculated only along the direction of x-ray tube movement. However, with the partial isocentric motion of the x-ray tube in breast tomosynthesis, small objects such as microcalcifications appear blurred (for instance, about 1-4 pixels in blur for a microcalcification in a human breast) in traditional SAA images in the direction perpendicular to the direction of tube motion. Some digital breast tomosynthesis algorithms reported in the literature utilize a traditional one-dimensional SAA method that is not wholly suitable for isocentric motion. In this paper, a point-by-point back projection (BP) method is described and compared with traditional SAA for the important clinical task of evaluating morphology of small objects such as microcalcifications. Impulse responses at different three-dimensional locations with five different combinations of imaging acquisition parameters were investigated. Reconstruction images of microcalcifications in a human subject were also evaluated. Results showed that with traditional SAA and 45 degrees view angle of tube movement with respect to the detector, at the same height above the detector, the in-plane blur artifacts were obvious for objects farther away from x-ray source. In a human subject, the appearance of calcifications was blurred in the direction orthogonal to the tube motion with traditional SAA. With point-by-point BP, the appearance of calcifications was sharper. The point-by-point BP method demonstrated improved rendition of microcalcifications in the direction perpendicular to the tube motion direction. With wide angles or for imaging of larger breasts, this point-by-point BP rather than the traditional SAA should also be considered as the basis of further deblurring algorithms that work in conjunction with the BP method.     

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.     

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.     

Evaluation of dual-energy subtraction of digital mammography images under conditions found in a commercial unit

Radiological contrast-to-noise ratio (CNR) is evaluated in subtracted images of microcalcifications in breast tissue. CNR is calculated for dual-kVp subtraction combining beams available in a Senographe 2000D, assuming single breast compression. Spectra were obtained from Boone et al (1997 Med. Phys. 24 1863-73), and the study was limited to lowest 25 kV Mo/Mo and highest 40 kV Rh/Rh beams, for 2.58 x 10(-4) C kg(-1) (1R) total exposure. For a standard case combining 25 kVp Mo/Mo and 40 kVp Rh/Rh beams, predicted maximum CNR for 300 microm calcification in 5 cm thick, 50% glandular, breast is about 1.2, below Rose's criterion for visualization. Total mean glandular doses are about 2.5 cGy for a standard case. The effect that input factors might have on predictions has been evaluated. Choice between alternative spectra can affect CNR by 50%. Assumed calcification composition leads to differences of 67% in calculated CNR, and assumed breast tissue composition can alter CNR by 45%; these results are weakly dependent on calcification or breast thickness, or on the assumed fraction of glandular tissue. CNR values are related to detected spectra effective energy. Calculations predict that above 37 kVp Mo/Mo beams are more energetic than Rh/Rh at the same kVp, due to beam hardening.     

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.     

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.     

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.     

Physical characteristics of GE Senographe Essential and DS digital mammography detectors

The purpose of this study was to investigate physical characteristics of two full field digital mammography (FFDM) systems (GE Senographe Essential and DS). Both are indirect conversion (x ray to light) alpha-Si flat panels coupled with a CsI(Tl) scintillator. The examined systems have the same pixel size (100 microm) but a different field of view: a conventional size 23 x 19.2 cm2 and a large field 24 X 30.7 cm2, specifically designed to image large breasts. In the GE Senographe Essential model relevant improvements in flat panel design were implemented and new deposition tools for metal, alpha-Si, and CsI(Tl) were introduced by GE. These changes in detector design are expected to be beneficial for advanced applications such as breast tomosynthesis. The presampling modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE) were measured for a wide range of exposure (25-240 microGy) with a RQA-M2 technique (28 kVp with a Mo/Mo target/filter combination and 2 mm of additional aluminum filtration). At 1, 2, and at 4 lp/mm MTF is equal to 0.9, 0.76, and 0.46 for the conventional field detector and to 0.85, 0.59, and 0.24 for the large field detector. The latter detector exhibits an improved NNPS due to a lower electronic noise and a better DQE that reaches 60%. In addition a contrast-detail analysis was performed with CDMAM 3.4 phantom and CDCOM software: GE Senographe DS showed statistically significant poorer detection ability in comparison with the GE Senographe Essential. These results could have been expected, at least qualitatively, considering the relative DQE of the two systems.     

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.     

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.     

Physical characterization of a prototype selenium-based full field digital mammography detector

The purpose of this study was to measure experimentally the physical performance of a prototype mammographic imager based on a direct detection, flat-panel array design employing an amorphous selenium converter with 70 microm pixels. The system was characterized for two different anode types, a molybdenum target with molybdenum filtration (Mo/Mo) and a tungsten target with rhodium filtration (W/Rh), at two different energies, 28 and 35 kVp, with approximately 2 mm added aluminum filtration. To measure the resolution, the presampled modulation transfer function (MTF) was measured using an edge method. The normalized noise power spectrum (NNPS) was measured by two-dimensional Fourier analysis of uniformly exposed mammograms. The detective quantum efficiencies (DQEs) were computed from the MTFs, the NNPSs, and theoretical ideal signal to noise ratios. The MTF was found to be close to its ideal limit and reached 0.2 at 11.8 mm(-1) and 0.1 at 14.1 mm(-1) for images acquired at an RQA-M2 technique (Mo/Mo anode, 28 kVp, 2 mm Al). Using a tungsten technique (MW2; W/Rh anode, 28 kVp, 2 mm Al), the MTF went to 0.2 at 11.2 mm(-1) and to 0.1 at 13.3 mm(-1). The DQE reached a maximum value of 54% at 1.35 mm(-1) for the RQA-M2 technique at 1.6 microC/kg and achieved a peak value of 64% at 1.75 mm(-1) for the tungsten technique (MW2) at 1.9 microC/kg. Nevertheless, the DQE showed strong exposure and frequency dependencies. The results indicated that the detector offered high MTFs and DQEs, but structured noise effects may require improved calibration before clinical implementation.     

Analysis of the detective quantum efficiency of a developmental detector for digital mammography

We are developing a modular detector for applications in full field digital mammography and for diagnostic breast imaging. The detector is based on a design that has been refined over the past decade for applications in x-ray crystallography [Kalata et al., Proc. SPIE 1345, 270-279 (1990); Phillips et al. ibid. 2009, 133-138 (1993), Phillips et al., Nucl. Instrum. Methods Phys. Rev. A 334, 621-630 (1993)]. The full field mammographic detector, currently undergoing clinical evaluation, is formed from a 19 cm x 28 cm phosphor screen, read out by a 2 x 3 array of butted charge-coupled device (CCD) modules. Each 2k x 2k CCD is optically coupled to the phosphor via a fiber optic taper with dimensions of 9.4 cm x 9.4cm at the phosphor. This paper describes the imaging performance of a two-module prototype, built using a similar design. In this paper we use cascaded linear systems analysis to develop a model for calculating the spatial frequency dependent noise power spectrum (NPS) and detective quantum efficiency (DQE) of the detector using the measured modulation transfer function (MTF). We compare results of the calculation with the measured NPS and DQE of the prototype. Calculated and measured DQEs are compared over a range of clinically relevant x-ray exposures and kVps. We find that for x-ray photon energies between 10 and 28 keV, the detector gain ranges between 2.5 and 3.7 CCD electrons per incident x-ray, or approximately 5-8 electrons per absorbed x ray. Using a Mo/Mo beam and acrylic phantom, over a detector entrance exposure range of approximately 10 to 80 mR, the volume under the measured 2-d NPS of the prototype detector is proportional to the x-ray exposure, indicating quantum limited performance. Substantial agreement between the calculated and measured values was obtained for the frequency and exposure dependent NPS and DQE over a range of tube voltage from 25 to 30 kVp.