Robert S. Saunders

A method for modifying the image quality parameters of digital radiographic images

A new computer simulation approach is presented that is capable of modeling several varieties of digital radiographic systems by their image quality characteristics. In this approach, the resolution and noise characteristics of ideal supersampled input images are modified according to input modulation transfer functions (MTFs) and noise power spectra (NPS). The modification process is separated into two routines-one for modification of the resolution and another for modification of the noise characteristics of the input image. The resolution modification routine blurs the input image by applying a frequency filter described by the input MTF. The resulting blurred image is then reduced to its final size to account for the sampling process of the digital system. The noise modification routine creates colored noise by filtering the frequency components of a white noise spectrum according to the input noise power. This noise is then applied to the image by a moving region of interest to account for variations in noise due to differences in attenuation. In order to evaluate the efficacy of the modification routines, additional routines were developed to assess the resolution and noise of digital images. The MTFs measured from the output images of the resolution modification routine were within 3% of the input MTF The NPS measured from the output images of the noise modification routine were within 2% of the input NPS. The findings indicate that the developed modification routines provide a good means of simulating the resolution and noise characteristics of digital radiographic systems for optimization or processing purposes.

Dose dependence of mass and microcalcification detection in digital mammography: free response human observer studies.

The purpose of this study was to evaluate the effect of dose reduction in digital mammography on the detection of two lesion types-malignant masses and clusters of microcalcifications. Two free-response observer studies were performed-one for each lesion type. Ninety screening images were retrospectively selected; each image was originally acquired under automatic exposure conditions, corresponding to an average glandular dose of 1.3 mGy for a standard breast (50 mm compressed breast thickness with 50% glandularity). For each study, one to three simulated lesions were added to each of 40 images (abnormals) while 50 were kept without lesions (normals). Two levels of simulated system noise were added to the images yielding two new image sets, corresponding to simulated dose levels of 50% and 30% of the original images (100%). The manufacturers standard display processing was subsequently applied to all images. Four radiologists experienced in mammography evaluated the images by searching for lesions and marking and assigning confidence levels to suspicious regions. The search data were analyzed using jackknife free-response (JA-FROC) methodology. For the detection of masses, the mean figure-of-merit (FOM) averaged over all readers was 0.74, 0.71, and 0.68 corresponding to dose levels of 100%, 50%, and 30%, respectively. These values were not statistically different from each other (F= 1.67, p=0.19) but showed a decreasing trend. In contrast, in the microcalcification study the mean FOM was 0.93, 0.67, and 0.38 for the same dose levels and these values were all significantly different from each other (F = 109.84, p < 0.0001). The results indicate that lowering the present dose level by a factor of two compromised the detection of microcalcifications but had a weaker effect on mass detection.

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.

Fundamental imaging characteristics of a slot-scan digital chest radiographic system

Our purpose in this study was to evaluate the fundamental image quality characteristics of a new slot-scan digital chest radiography system (ThoraScan, Delft Imaging Systems/Nucletron, Veenendaal, The Netherlands). The linearity of the system was measured over a wide exposure range at 90, 117, and 140 kVp with added Al filtration. System uniformity and reproducibility were established with an analysis of images from repeated exposures. The modulation transfer function (MTF) was evaluated using an established edge method. The noise power spectrum (NPS) and the detective quantum efficiency (DQE) of the system were evaluated at the three kilo-voltages over a range of exposures. Scatter fraction (SF) measurements were made using a posterior beam stop method and a geometrical chest phantom. The system demonstrated excellent linearity, but some structured nonuniformities. The 0.1 MTF values occurred between 3.3-3.5 mm(-1). The DQE(0.15) and DQE(2.5) were 0.21 and 0.07 at 90 kVp, 0.18 and 0.05 at 117 kVp, and 0.16 and 0.03 at 140 kVp, respectively. The system exhibited remarkably lower SFs compared to conventional full-field systems with anti-scatter grid, measuring 0.13 in the lungs and 0.43 in the mediastinum. The findings indicated that the slot-scan design provides marked scatter reduction leading to high effective DQE (DQEeff) of the system and reduced patient dose required to achieve high image quality.

An Anthropomorphic Breast Model for Breast Imaging Simulation and Optimization

RATIONALE AND OBJECTIVES Optimization studies for x-ray-based breast imaging systems using computer simulation can greatly benefit from a phantom capable of modeling varying anatomical variability across different patients. This study aimed to develop a three-dimensional phantom model with realistic and randomizable anatomical features. MATERIALS AND METHODS A voxelized breast model was developed consisting of an outer layer of skin and subcutaneous fat, a mixture of glandular and adipose, stochastically generated ductal trees, masses, and microcalcifications. Randomized realization of the breast morphology provided a range of patient models. Compression models were included to represent the breast under various compression levels along different orientations. A Monte Carlo (MC) simulation code was adapted to simulate x-ray based imaging systems for the breast phantom. Simulated projections of the phantom at different angles were generated and reconstructed with iterative methods, simulating mammography, breast tomosynthesis, and computed tomography (CT) systems. Phantom dose maps were further generated for dosimetric evaluation. RESULTS Region of interest comparisons of simulated and real mammograms showed strong similarities in terms of appearance and features. Noise-power spectra of simulated mammographic images demonstrated that the phantom provided target properties for anatomical backgrounds. Reconstructed tomosynthesis and CT images and dose maps provided corresponding data from a single breast enabling optimization studies. Dosimetry result provided insight into the dose distribution difference between modalities and compression levels. CONCLUSION The anthropomorphic breast phantom, combined with the MC simulation platform, generated a realistic model for a breast imaging system. The developed platform is expected to provide a versatile and powerful framework for optimizing volumetric breast imaging systems.

A mathematical model platform for optimizing a multiprojection breast imaging system

Multiprojection imaging is a technique in which a plurality of digital radiographic images of the same patient are acquired within a short interval of time from slightly different angles. Information from each image is combined to determine the final diagnosis. Projection data are either reconstructed into slices as in the case of tomosynthesis or analyzed directly as in the case of multiprojection correlation imaging technique, thereby avoiding reconstruction artifacts. In this study, the authors investigated the optimum geometry of acquisitions of a multiprojection breast correlation imaging system in terms of the number of projections and their total angular span that yield maximum performance in a task that models clinical decision. Twenty-five angular projections of each breast from 82 human subjects in our breast tomosynthesis database were each supplemented with a simulated 3 mm mass. An approach based on Laguerre-Gauss channelized Hotelling observer was developed to assess the detectability of the mass in terms of receiver operating characteristic (ROC) curves. Two methodologies were developed to integrate results from individual projections into one combined ROC curve as the overall figure of merit. To optimize the acquisition geometry, different components of acquisitions were changed to investigate which one of the many possible configurations maximized the area under the combined ROC curve. Optimization was investigated under two acquisition dose conditions corresponding to a fixed total dose delivered to the patient and a variable dose condition, based on the number of projections used. In either case, the detectability was dependent on the number of projections used, the total angular span of those projections, and the acquisition dose level. In the first case, the detectability approximately followed a bell curve as a function of the number of projections with the maximum between 8 and 16 projections spanning angular arcs of about 23 degrees-45 degrees, respectively. In the second case, the detectability increased with the number of projections approaching an asymptote at 11-17 projections for an angular span of about 45 degrees. These results indicate the inherent information content of the multi-projection image data reflecting the relative role of quantum and anatomical noise in multiprojection breast imaging. The optimization scheme presented here may be applied to any multiprojection imaging modalities and may be extended by including reconstruction in the case of digital breast tomosynthesis and breast computed tomography.

Does image quality matter? Impact of resolution and noise on mammographic task performance

The purpose of this study was to examine the effects of different resolution and noise levels on task performance in digital mammography. This study created an image set with images at three different resolution levels, corresponding to three digital display devices, and three different noise levels, with noise magnitudes similar to full clinical dose, half clinical dose, and quarter clinical dose. The images were read by five experienced breast imaging radiologists. The data were then analyzed to compute two accuracy statistics (overall classification accuracy and lesion detection accuracy) and performance at four diagnostic tasks (detection of microcalcifications, benign masses, malignant masses, and discrimination of benign and malignant masses). Human observer results showed decreasing display resolution had little effect on overall classification accuracy and individual diagnostic task performance, but increasing noise caused overall classification accuracy to decrease by a statistically significant 21% as the breast dose went to one quarter of its normal clinical value. The noise effects were most prominent for the tasks of microcalcification detection and mass discrimination. When the noise changed from full clinical dose to quarter clinical dose, the microcalcification detection performance fell from 89% to 67% and the mass discrimination performance decreased from 93% to 79%, while malignant mass detection performance remained relatively constant with values of 88% and 84%, respectively. As a secondary aim, the image set was also analyzed by two observer models to examine whether their performance was similar to humans. Observer models differed from human observers and each other in their sensitivity to resolution degradation and noise. The primary conclusions of this study suggest that quantum noise appears to be the dominant image quality factor in digital mammography, affecting radiologist performance much more profoundly than display resolution.

Resolution and noise measurements of five CRT and LCD medical displays.

The performance of soft-copy displays plays a significant role in the overall image quality of a digital radiographic system. In this work, we discuss methods to characterize the resolution and noise of both cathode ray tube (CRT) and liquid crystal display (LCD) devices. We measured the image quality of five different commercial display devices, representing both CRT and LCD technologies, using a high-quality charge-coupled device (CCD) camera. The modulation transfer function (MTF) was calculated using the line technique, correcting for the MTF of the CCD camera and the display pixel size. The normalized noise power spectrum (NPS) was computed from two-dimensional Fourier analysis of uniform images. To separate the effects of pixel structure from interpixel luminance variations, we created structure-free images by eliminating the pixel structures of the display device. The NPS was then computed from these structure-free images to isolate interpixel luminance variations. We found that the MTF of LCDs remained close to the theoretical limit dictated by their inherent pixel size (0.85 +/- 0.08 at Nyquist frequency), in contrast to the MTF for the two CRT displays, which dropped to 0.15 +/- 0.08 at the Nyquist frequency. However, the NPS of LCDs showed significant peaks due to the subpixel structure, while the NPS of CRT displays exhibited a nearly flat power spectrum. After removing the pixel structure, the structured noise peaks for LCDs were eliminated and the overall noise magnitude was significantly reduced. The average total noise-to-signal ratio for CRT displays was 6.55% +/- 0.59%, of which 6.03% +/- 0.24% was due to interpixel luminance variations, while LCD displays had total noise to signal ratios of 46.1% +/- 5.1% of which 1.50% +/- 0.41% were due to interpixel luminance variations. Depending on the extent of the blurring and prewhitening processes of the human visual system, the magnitude of the display noise (including pixel structure) potentially perceived by the observer was reduced to 0.43% +/- 0.01% (accounting for blurring only) and 0.40 +/- 0.01% (accounting for blurring and prewhitening) for CRTs, and 1.02% +/- 0.22% (accounting for blurring only) and 0.36% +/- 0.08% (accounting for blurring and prewhitening) for LCDs.

Effect of dose reduction on the detection of mammographic lesions: A mathematical observer model analysis

The effect of reduction in dose levels normally used in mammographic screening procedures on the detection of breast lesions were analyzed. Four types of breast lesions were simulated and inserted into clinically-acquired digital mammograms. Dose reduction by 50% and 75% of the original clinically-relevant exposure levels were simulated by adding corresponding simulated noise into the original mammograms. The mammograms were converted into luminance values corresponding to those displayed on a clinical soft-copy display station and subsequently analyzed by Laguerre-Gauss and Gabor channelized Hotelling observer models for differences in detectability performance with reduction in radiation dose. Performance was measured under a signal known exactly but variable detection task paradigm in terms of receiver operating characteristics (ROC) curves and area under the ROC curves. The results suggested that luminance mapping of digital mammograms affects performance of model observers. Reduction in dose levels by 50% lowered the detectability of masses with borderline statistical significance. Dose reduction did not have a statistically significant effect on detection of microcalcifications. The model results indicate that there is room for optimization of dose level in mammographic screening procedures.

Micro-CT imaging of breast tumors in rodents using a liposomal, nanoparticle contrast agent.

A long circulating liposomal, nanoscale blood pool agent encapsulating traditional iodinated contrast agent (65 mg I/mL) was used for micro-computed tomography (CT) imaging of rats implanted with R3230AC mammary carcinoma. Three-dimensional vascular architecture of tumors was imaged at 100-micron isotropic resolution. The image data showed good qualitative correlation with pathologic findings. The approach holds promise for studying tumor angiogenesis and for evaluating anti-angiogenesis therapies.