Dual-energy imaging in full-field digital mammography: a phantom study

A dual-energy technique which employs the basis decomposition method is being investigated for application to digital mammography. A three-component phantom, made up of plexiglas, polyethylene, and water, was doubly exposed with the full-field digital mammography system manufactured by General Electric. The 'low' and 'high' energy images were recorded with a Mo/Mo anode-filter combination and a Rh/Rh combination, respectively. The total dose was kept within the acceptable levels of conventional mammography. The first hybrid images obtained with the dual-energy algorithm are presented in comparison with a conventional radiograph of the phantom. Image-quality characteristics at contrast cancellation angles between plexiglas and water are discussed. Preliminary results show that a combination of a standard Mo-anode 28 kV radiograph with a Rh-anode 49 kV radiograph provides the best compromise between image-quality and dose in the hybrid image.     

Comparison of full-field digital mammography to screen-film mammography with respect to contrast and spatial resolution in tissue equivalent breast phantoms

To determine if the improved contrast resolution of full-field digital mammography (FFDM) with reduced spatial resolution allows for superior or equal phantom object detection compared with screen-film mammography (SFM). Tissue equivalent breast phantoms simulating an adipose to glandular ratio of 50/50,30/70, and 20/80 were imaged according to each manufacturers' recommendation with four full-field digital mammography units (Fuji, Sectra, Fischer, and General Electric) and a screen-film mammography unit (MammoMatII 2000, Siemens, Munich, Germany). A total of 20 images were obtained in both hard- and soft-copy formats. For the purpose of soft-copy display, the screen-film hard-copy images were digitized with a 50 microm micron scanner. Six radiologists, experts in breast imaging, and three physicists, experts in scoring mammography phantoms, participated in a reader study where each reader scored each phantom for visibility of line-pairs and for 24 objects (fibers, clusters of specks, and masses). The data were recorded, entered into a database, and analyzed by a mixed-effect model. The limiting spatial resolution in line-pairs per millimeter visible with the digital units was less, regardless of display modality used, than that provided by the screen-film unit. The difference was statistically significant for the General Electric (p < 0.01) and Fuji digital mammography units (p = 0.03). With respect to the number of visible objects, a statistically significant higher number could be detected with the screen-film unit as compared to the Fischer (p < 0.01) and Sectra (p < 0.01) digital mammography units, but there was no significant difference between the other digital units and screen film. Overall, there was significantly better performance on the 50/50 phantom than with the 30/70 and 20/80 phantoms (p = 0.01, p < 0.01) for object visibility. For the digital mammography units, soft-copy display performed better than hard-copy display for the Fischer and Sectra images, but worse for Fuji and General Electric. In addition, soft-copy display of digitized screen-film images was significantly better than hard-copy display (p =0.02) of the original screen films for object visibility, but worse for spatial resolution. The higher contrast resolution of the FFDM units tested did not result in improved detection of line-pair resolution or objects in the phantoms tested versus screen-film mammography. The phantom performance of a digital mammography unit seems to be influenced by the type of detection task (line-pair resolution versus object visibility), the display modality (soft-copy versus hard-copy) chosen to score the phantoms, and the parenchymal pattern composition of the phantom.     

A new test phantom with different breast tissue compositions for image quality assessment in conventional and digital mammography

Our objective is to describe a new test phantom that permits the objective assessment of image quality in conventional and digital mammography for different types of breast tissue. A test phantom, designed to represent a compressed breast, was made from tissue equivalent materials. Three separate regions, with different breast tissue compositions, are used to evaluate low and high contrast resolution, spatial resolution and image noise. The phantom was imaged over a range of kV using a Contour 2000 (Bennett) mammography unit with a Kodak MinR 2190-MinR L screen-film combination and a Senograph 2000D (General Electric) digital mammography unit. Objective image quality assessments for different breast tissue compositions were performed using the phantom for conventional and digital mammography. For a similar mean glandular dose (MGD), the digital system gives a significantly higher contrast-to-noise ratio (CNR) than the screen-film system for 100% glandular tissue. In conclusion, in mammography, a range of exposure conditions is used for imaging because of the different breast tissue compositions encountered clinically. Ideally, the patient dose-image quality relationship should be optimized over the range of exposure conditions. The test phantom presented in this work permits image quality parameters to be evaluated objectively for three different types of breast tissue. Thus, it is a useful tool for optimizing the patient dose-image quality relationship.     

Dual-energy digital mammography for calcification imaging: noise reduction techniques

We have previously developed a dual-energy digital mammography (DEDM) technique for calcification imaging under full-field imaging conditions using a commercially available flat-panel based digital mammography system. Although dual-energy (DE) imaging could suppress the obscuration of calcifications by tissue-structure background, it also increases the intrinsic noise in the DE images. Here we report on the effects of three different noise reduction techniques on DE calcification images: a simple smoothing (boxcar) filter applied to the DE image, a median filter applied to the HE image prior to the computation of the DE image and an adaptation of the Kalender's correlated-noise reduction (KNR) technique for DEDM. We compared the different noise reduction techniques by evaluating their effects on DE calcification images of a 5 cm thick breast-tissue-equivalent slab with continuously varying glandular-tissue ratio superimposed with calcium carbonate crystals of various sizes that simulate calcifications. Evaluations of different noise reducing techniques were performed by comparison of the root-mean-square signal in background regions (no calcifications present) of the DE calcification images and the contrast-to-noise ratios (CNR) of the calcifications in the DE calcification images. Amongst the different noise reduction techniques evaluated in this study, the KNR method was found to be most effective in reducing the image noise and increasing the calcification visibility (or CNR), closely followed by the HE median filter technique. Although the simple smoothing (boxcar) filter reduced the noise, it did not improve calcification visibility. The visible calcification threshold size with DEDM over smoothly varying background at screening mammography doses, assuming a CNR threshold of 4, was estimated to be around 250 microm with both the HE median filter and the KNR techniques. The quality of DE images with noise reduction techniques based on phantom studies were verified with DE images of an animal-tissue phantom that consisted of calcifications superimposed over more realistic tissue structures.     

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.     

The value of scatter removal by a grid in full field digital mammography

Our objective in this study was to investigate the usefulness of an anti-scatter grid in digital mammography using a contrast detail phantom. The mammography system we investigated was a GE Senographe 2000D. We carried out phantom measurements under various conditions with and without using the anti-scatter grid. A new version of the CDMAM phantom (version 3.4) was used. This phantom consists of a matrix of square cells with disks of varying size and contrast. For given exposure conditions detectability of these disks can be determined and used for construction of contrast detail curves. Previously, a computer program was developed at our institute that performs a fully automatic analysis of the phantom recordings using the ideal observer model. Breast thickness was simulated by a homogeneous layer of PMMA in the range of 1 to 7 cm. Series of images were recorded for different KeV and target-filter combinations depending on the simulated thickness. The dose was kept constant for each thickness with and without using a grid. It appeared that image quality improved for simulated breast thickness below 5 cm when the grid was removed. In the range from 5 to 7 cm contrast detail curves obtained with or without a grid were similar. Results suggest that for compressed breast thickness in the range of 1 to 7 cm a grid might not be needed in the digital mammography system we investigated. Below 5 cm, omitting the grid may allow lower dose to the patient without losing image quality.     

Image quality assessment in digital mammography: part II. NPWE as a validated alternative for contrast detail analysis

Assessment of image quality for digital x-ray mammography systems used in European screening programs relies mainly on contrast-detail CDMAM phantom scoring and requires the acquisition and analysis of many images in order to reduce variability in threshold detectability. Part II of this study proposes an alternative method based on the detectability index (d') calculated for a non-prewhitened model observer with an eye filter (NPWE). The detectability index was calculated from the normalized noise power spectrum and image contrast, both measured from an image of a 5 cm poly(methyl methacrylate) phantom containing a 0.2 mm thick aluminium square, and the pre-sampling modulation transfer function. This was performed as a function of air kerma at the detector for 11 different digital mammography systems. These calculated d' values were compared against threshold gold thickness (T) results measured with the CDMAM test object and against derived theoretical relationships. A simple relationship was found between T and d', as a function of detector air kerma; a linear relationship was found between d' and contrast-to-noise ratio. The values of threshold thickness used to specify acceptable performance in the European Guidelines for 0.10 and 0.25 mm diameter discs were equivalent to threshold calculated detectability indices of 1.05 and 6.30, respectively. The NPWE method is a validated alternative to CDMAM scoring for use in the image quality specification, quality control and optimization of digital x-ray systems for screening mammography.     

Improved image quality in digital mammography with image processing

PURPOSE: The effect of image processing, specifically Bayesian image estimation (BIE), on digital mammographic images is studied. BIE is an iterative, nonlinear statistical estimation technique that has previously been used in chest radiography to reduce image scatter content and improve the contrast-to-noise ratio (CNR). We adapt this technique to digital mammography and examine its effect. METHODS/MATERIALS: Images of the American College of Radiologists (ACR) breast phantom were acquired on a calibrated digital mammography system at a normal mammographic exposure both with and without a grid. An iterative Bayesian estimation algorithm was formulated and used to process the images acquired without a grid. Quantitative scatter fractions were measured and compared for the image acquired with the grid, the image acquired without the grid, and the image acquired without the grid and processed by the Bayesian algorithm. CNR values were also computed for the four visible masses in the ACR phantom before and after processing and compared to a grid. RESULTS: Initial images acquired without an antiscatter grid had scatter fractions of 0.46. Processing this image with BIE reduced the scatter content to under 0.04. In comparison, the image acquired with a grid had scatter of 0.19. BIE processing accounted for CNR improvements from 29% to 219% for the masses seen in the ACR phantom as compared to the unprocessed image. Visibility of the four masses in the phantom was improved. CONCLUSIONS: Bayesian image estimation can be used with digital mammography to reduce scatter fractions. This technique is very useful as it can reduce scatter content effectively without introducing any adverse effects, such as grid line aliasing. Bayesian processing can also increase image CNR, which may potentially increase the visualization of subtle masses. Preliminary work shows an improvement in CNR to values greater than that provided by a standard grid.     

Automated analysis of phantom images for the evaluation of long-term reproducibility in digital mammography

The performance of an automatic software package was evaluated with phantom images acquired by a full-field digital mammography unit. After the validation, the software was used, together with a Leeds TORMAS test object, to model the image acquisition process. Process modelling results were used to evaluate the sensitivity of the method in detecting changes of exposure parameters from routine image quality measurements in digital mammography, which is the ultimate purpose of long-term reproducibility tests. Image quality indices measured by the software included the mean pixel value and standard deviation of circular details and surrounding background, contrast-to-noise ratio and relative contrast; detail counts were also collected. The validation procedure demonstrated that the software localizes the phantom details correctly and the difference between automatic and manual measurements was within few grey levels. Quantitative analysis showed sufficient sensitivity to relate fluctuations in exposure parameters (kV(p) or mAs) to variations in image quality indices. In comparison, detail counts were found less sensitive in detecting image quality changes, even when limitations due to observer subjectivity were overcome by automatic analysis. In conclusion, long-term reproducibility tests provided by the Leeds TORMAS phantom with quantitative analysis of multiple IQ indices have been demonstrated to be effective in predicting causes of deviation from standard operating conditions and can be used to monitor stability in full-field digital mammography.     

Étude de la faisabilité du scorage automatique de fantômes mammographiques par traitement d’image

Nowadays, X-ray mammography is one of the most effective methods for early and reliable breast cancer detection and diagnosis. Periodic quality control in mammographic facilities is necessary to provide high quality mammograms. This evaluation is done by visual observation of mammographic phantom films. To prepare Full Field Digital Mammography advent and to get rid of this subjective quality control, digital image and signal processing techniques may be used in order to make this control easier and more objective. This paper presents an automatic method for scoring mammographic phantoms using digital image processing. Phantom films were first digitized and images containing microcalcifications, masses and fibers were extracted. Theses noisy and low contrasted images were preprocessed using an adaptive contrast enhancement method and then segmented in order to extract objects embedded in phantom images. Nine digitized phantom films were studied and results show that a more objective quality control evaluation of mammographic facilities can be done using digital image processing techniques on phantom images.     

Dual-energy digital mammography for calcification imaging: scatter and nonuniformity corrections

Mammographic images of small calcifications, which are often the earliest signs of breast cancer, can be obscured by overlapping fibroglandular tissue. We have developed and implemented a dual-energy digital mammography (DEDM) technique for calcification imaging under full-field imaging conditions using a commercially available aSi:H/CsI:Tl flat-panel based digital mammography system. The low- and high-energy images were combined using a nonlinear mapping function to cancel the tissue structures and generate the dual-energy (DE) calcification images. The total entrance-skin exposure and mean-glandular dose from the low- and high-energy images were constrained so that they were similar to screening-examination levels. To evaluate the DE calcification image, we designed a phantom using calcium carbonate crystals to simulate calcifications of various sizes (212-425 microm) overlaid with breast-tissue-equivalent material 5 cm thick with a continuously varying glandular-tissue ratio from 0% to 100%. We report on the effects of scatter radiation and nonuniformity in x-ray intensity and detector response on the DE calcification images. The nonuniformity was corrected by normalizing the low- and high-energy images with full-field reference images. Correction of scatter in the low- and high-energy images significantly reduced the background signal in the DE calcification image. Under the current implementation of DEDM, utilizing the mammography system and dose level tested, calcifications in the 300-355 microm size range were clearly visible in DE calcification images. Calcification threshold sizes decreased to the 250-280 microm size range when the visibility criteria were lowered to barely visible. Calcifications smaller than approximately 250 microm were usually not visible in most cases. The visibility of calcifications with our DEDM imaging technique was limited by quantum noise, not system noise.     

Contrast Detail Phantom Comparison on a Commercially Available Unit. Digital Breast Tomosynthesis (DBT) versus Full-Field Digital Mammography (FFDM)

The performance of a commercial digital mammographic system working in 2D planar versus tomosynthesis mode was evaluated in terms of the image signal difference to noise ratio (SDNR). A contrast detail phantom was obtained embedding 1 cm Plexiglas, including 49 holes of different diameter and depth, between two layers containing a breast-simulating material. The phantom was exposed with the details plane perpendicular to the X-ray beam using the manufacturer’s standard clinical breast acquisition parameters. SDNR in the digital breast tomosynthesis (DBT) images was higher than that of the full-field digital mammography (FFDM) for 38 out of 49 details in complex background conditions. These differences (p < 0.05) are statistically significant for 19 details out of 38. The relative SDNR results for DBT and FFDM images showed a dependence on the diameter of the details considered. This paper proposes an initial framework for a global image quality evaluation for commercial systems that can operate with different image acquisition modality using the same detector.     

The impact of calibration phantom errors on dual-energy digital mammography

Microcalcification is one of the earliest and main indicators of breast cancer. Because dual-energy digital mammography could suppress the contrast between the adipose and glandular tissues of the breast, it is considered a promising technique that will improve the detection of microcalcification. In dual-energy digital mammography, the imaged object is a human breast, while in calibration measurements only the phantoms of breast tissue equivalent materials are available. Consequently, the differences between phantoms and breast tissues will lead to calibration phantom errors. Based on the dual-energy imaging model, formulae of calibration phantom errors are derived in this paper. Then, this type of error is quantitatively analyzed using publicly available data and compared with other types of error. The results demonstrate that the calibration phantom error is large and dominant in dual-energy mammography, seriously decreasing calculation precision. Further investigations on the physical meaning of calibration phantom error reveal that the imaged objects with the same glandular ratio have identical calibration phantom error. Finally, an error correction method is proposed based on our findings.     

Evaluation of detector dynamic range in the x-ray exposure domain in mammography: a comparison between film-screen and flat panel detector systems

Digital detectors in mammography have wide dynamic range in addition to the benefit of decoupled acquisition and display. How wide the dynamic range is and how it compares to film-screen systems in the clinical x-ray exposure domain are unclear. In this work, we compare the effective dynamic ranges of film-screen and flat panel mammography systems, along with the dynamic ranges of their component image receptors in the clinical x-ray exposure domain. An ACR mammography phantom was imaged using variable mAs (exposure) values for both systems. The dynamic range of the contrast-limited film-screen system was defined as that ratio of mAs (exposure) values for a 26 kVp Mo/Mo (HVL=0.34 mm Al) beam that yielded passing phantom scores. The same approach was done for the noise-limited digital system. Data from three independent observers delineated a useful phantom background optical density range of 1.27 to 2.63, which corresponded to a dynamic range of 2.3 +/- 0.53. The digital system had a dynamic range of 9.9 +/- 1.8, which was wider than the film-screen system (p<0.02). The dynamic range of the film-screen system was limited by the dynamic range of the film. The digital detector, on the other hand, had an estimated dynamic range of 42, which was wider than the dynamic range of the digital system in its entirety by a factor of 4. The generator/tube combination was the limiting factor in determining the digital system's dynamic range.     

Optimization of accelerator target and detector for portal imaging using Monte Carlo simulation and experiment

Megavoltage portal images suffer from poor quality compared to those produced with kilovoltage x-rays. Several authors have shown that the image quality can be improved by modifying the linear accelerator to generate more low-energy photons. This work addresses the problem of using Monte Carlo simulation and experiment to optimize the beam and detector combination to maximize image quality for a given patient thickness. A simple model of the whole imaging chain was developed for investigation of the effect of the target parameters on the quality of the image. The optimum targets (6 mm thick aluminium and 1.6 mm copper) were installed in an Elekta SL25 accelerator. The first beam will be referred to as A16 and the second as Cu1.6. A tissue-equivalent contrast phantom was imaged with the 6 MV standard photon beam and the experimental beams with standard radiotherapy and mammography film/screen systems. The arrangement with a thin Al target/mammography system improved the contrast from 1.4 cm bone in 5 cm water to 19% compared with 2% for the standard arrangement of a thick, high-Z target/radiotherapy verification system. The linac/phantom/detector system was simulated with the BEAM/EGS4 Monte Carlo code. Contrast calculated from the predicted images was in good agreement with the experiment (to within 2.5%). The use of MC techniques to predict images accurately, taking into account the whole imaging system, is a powerful new method for portal imaging system design optimization.     

Quantification of breast arterial calcification using full field digital mammography

Breast arterial calcification is commonly detected on some mammograms. Previous studies indicate that breast arterial calcification is evidence of general atherosclerotic vascular disease and it may be a useful marker of coronary artery disease. It can potentially be a useful tool for assessment of coronary artery disease in women since mammography is widely used as a screening tool for early detection of breast cancer. However, there are currently no available techniques for quantification of calcium mass using mammography. The purpose of this study was to determine whether it is possible to quantify breast arterial calcium mass using standard digital mammography. An anthropomorphic breast phantom along with a vessel calcification phantom was imaged using a full field digital mammography system. Densitometry was used to quantify calcium mass. A calcium calibration measurement was performed at each phantom thickness and beam energy. The known (K) and measured (M) calcium mass on 5 and 9 cm thickness phantoms were related by M=0.964K -0.288 mg (r=0.997 and SEE=0.878 mg) and M=1.004K+0.324 mg (r=0.994 and SEE = 1.32 mg), respectively. The results indicate that accurate calcium mass measurements can be made without correction for scatter glare as long as careful calcium calibration is made for each breast thickness. The results also indicate that composition variations and differences of approximately 1 cm between calibration phantom and breast thickness introduce only minimal error in calcium measurement. The uncertainty in magnification is expected to cause up to 5% and 15% error in calcium mass for 5 and 9 cm breast thicknesses, respectively. In conclusion, a densitometry technique for quantification of breast arterial calcium mass was validated using standard full field digital mammography. The results demonstrated the feasibility and potential utility of the densitometry technique for accurate quantification of breast arterial calcium mass using standard digital mammography.     

Amorphous selenium flat panel detectors for digital mammography: validation of a NPWE model observer with CDMAM observer performance experiments

Model observers have been developed which incorporate a specific imaging task, system performance, and human observer characteristics and can potentially overcome some of the limitations in using detective quantum efficiency for optimization and comparison of detectors. In this paper, a modified nonprewhitening matched filter (NPWE) model observer was developed and validated to predict object detectability for an amorphous selenium (a-Se) direct flat-panel imager (FPI) where aliasing is severe. A preclinical a-Se digital mammography FPI with 85 microm pixel size was used in this investigation. Its physical imaging properties including modulation transfer function (MTF), noise power spectrum, and DQE were fully characterized. An observer performance study was conducted by imaging the CDMAM 3.4 contrast-detail phantom designed specifically for digital mammography and presenting these images to a panel of seven observers. X-ray attenuation and scatter due to the phantom were determined experimentally for use in development of the model observer. The observer study results were analyzed via threshold averaging and signal detection theory (SDT) based techniques to produce contrast-detail curves where threshold contrast is plotted as a function of disk diameter. Validity of the model was established using SDT analysis of the experimental data. The effect of aliasing on the detectability of small diameter disks was determined using the NPWE model observer. The signal spectrum was calculated using the presampling MTF of the detector with and without including the aliased terms. Our results indicate that the NPWE model based on Fourier domain parameters provides reasonable prediction of object detectability for the signal-known-exactly task in uniform image noise for a-Se direct FPI.     

Effects of magnification and zooming on depth perception in digital stereomammography: an observer performance study

We are evaluating the application of stereoscopic imaging to digital mammography. In the current study, we investigated the effects of magnification and zooming on depth perception. A modular phantom was designed which contained six layers of 1-mm-thick Lexan plates, each spaced 1 mm apart. Eight to nine small, thin nylon fibrils were pasted on each plate in horizontal or vertical orientations such that they formed 25 crossing fibril pairs in a projected image. The depth separation between each fibril pair ranged from 2 to 10 mm. A change in the order of the Lexan plates changed the depth separation of the two fibrils in a pair. Stereoscopic image pairs of the phantom were acquired with a GE full-field digital mammography system. Three different phantom configurations were imaged. All images were obtained using a Rh target/Rh filter spectrum at 30 kVp tube potential and a +/- 3 stereo shift angle. Images were acquired in both contact and 1.8X magnification geometry and an exposure range of 4 to 63 mAs was employed. The images were displayed on a Barco monitor driven by a Metheus stereo graphics board and viewed with LCD stereo glasses. Five observers participated in the study. Each observer visually judged whether the vertical fibril was in front of or behind the horizontal fibril in each fibril pair. It was found that the accuracy of depth discrimination increased with increasing fibril depth separation and x-ray exposure. The accuracy was not improved by electronic display zooming of the contact stereo images by 2X. Under conditions of high noise (low mAs) and small depth separation between the fibrils, the observers' depth discrimination ability was significantly better in stereo images acquired with geometric magnification than in images acquired with a contact technique and displayed with or without zooming. Under our experimental conditions, a 2 mm depth discrimination was achieved with over 60% accuracy on contact images with and without zooming, and with over 90% accuracy on magnification images. This study indicates that stereoscopic imaging, especially with magnification, may be useful for visualizing the spatial distribution of microcalcifications in a cluster and for differentiating overlapping tissues from masses on mammograms.     

Scanned-projection digital mammography

The effectiveness of film-screen mammography is limited by tradeoffs between latitude and contrast, film granularity, and the need to increase dose when antiscatter methods are used. We are currently developing a scanned-projection digital mammography (SPDM) system to overcome these limitations. The system consists of a pair of scanning slits, a high-resolution x-ray image intensifier tube, a linear photodiode array, and a digital display. The detective quantum efficiency of the SPDM system at spatial frequencies up to 3 cycles/mm is similar to that of mammographic film-screen combinations, but is lower at high frequencies. For low-contrast objects as small as 0.1 mm in diameter, the signal-to-noise ratio is currently equal to that of optimally exposed mammographic film-screen images for equal dose to the breast and superior for regions which would be underexposed or overexposed on film. This is achieved by the use of a low-noise detector system, geometric magnification, and scatter elimination. Images of a contrast-detail phantom and excised breast tissue illustrate the superior contrast sensitivity of SPDM.     

On the noise variance of a digital mammography system

A recent paper by Cooper et al. [Med. Phys. 30, 2614-2621 (2003)] contains some apparently anomalous results concerning the relationship between pixel variance and x-ray exposure for a digital mammography system. They found an unexpected peak in a display domain pixel variance plot as a function of 1/mAs (their Fig. 5) with a decrease in the range corresponding to high display data values, corresponding to low x-ray exposures. As they pointed out, if the detector response is linear in exposure and the transformation from raw to display data scales is logarithmic, then pixel variance should be a monotonically increasing function in the figure. They concluded that the total system transfer curve, between input exposure and display image data values, is not logarithmic over the full exposure range. They separated data analysis into two regions and plotted the logarithm of display image pixel variance as a function of the logarithm of the mAs used to produce the phantom images. They found a slope of minus one for high mAs values and concluded that the transfer function is logarithmic in this region. They found a slope of 0.6 for the low mAs region and concluded that the transfer curve was neither linear nor logarithmic for low exposure values. It is known that the digital mammography system investigated by Cooper et al. has a linear relationship between exposure and raw data values [Vedantham et al., Med. Phys. 27, 558-567 (2000)]. The purpose of this paper is to show that the variance effect found by Cooper et al. (their Fig. 5) arises because the transformation from the raw data scale (14 bits) to the display scale (12 bits), for the digital mammography system they investigated, is not logarithmic for raw data values less than about 300 (display data values greater than about 3300). At low raw data values the transformation is linear and prevents over-ranging of the display data scale. Parametric models for the two transformations will be presented. Results of pixel variance measurements made on raw data images will be presented. The experimental data are in good agreement with those of Cooper et al. It will be shown that the slope of 0.6 found by Cooper et al. for the log-log plot at low exposure is not due to transfer function nonlinearity, it occurs because of an additive variance term-possibly due to electronic noise. It will also be shown, using population statistics from clinical images, that raw data values below 300 are rare in tissue areas. Those tissue areas with very low raw data values are within about a millimeter of the chest wall or in very dense muscle at comers of images.