Ehsan Samei

A method for measuring the presampled MTF of digital radiographic systems using an edge test device

The modulation transfer function (MTF) of radiographic systems is frequently evaluated by measuring the systems line spread function (LSF) using narrow slits. The slit method requires precise fabrication and alignment of a slit and high radiation exposure. An alternative method for determining the MTF uses a sharp, attenuating edge device. We have constructed an edge device from a 250-microm-thick lead foil laminated between two thin slabs of acrylic. The device is placed near the detector and aligned with the aid of a laser beam and a holder such that a polished edge is parallel to the x-ray beam. A digital image of the edge is processed to obtain the presampled MTF. The image processing includes automated determination of the edge angle, reprojection, sub-binning, smoothing of the edge spread function (ESF), and spectral estimation. This edge method has been compared to the slit method using measurements on standard and high-resolution imaging plates of a digital storage phosphor (DSP) radiography system. The experimental results for both methods agree with a mean MTF difference of 0.008. The edge method provides a convenient measurement of the presampled MTF for digital radiographic systems with good response at low frequencies.

Low-tube-voltage, high-tube-current multidetector abdominal CT: improved image quality and decreased radiation dose with adaptive statistical iterative reconstruction algorithm--initial clinical experience.

PURPOSE To investigate whether an adaptive statistical iterative reconstruction (ASIR) algorithm improves the image quality at low-tube-voltage (80-kVp), high-tube-current (675-mA) multidetector abdominal computed tomography (CT) during the late hepatic arterial phase. MATERIALS AND METHODS This prospective, single-center HIPAA-compliant study was institutional review board approved. Informed patient consent was obtained. Ten patients (six men, four women; mean age, 63 years; age range, 51-77 years) known or suspected to have hypervascular liver tumors underwent dual-energy 64-section multidetector CT. High- and low-tube-voltage CT images were acquired sequentially during the late hepatic arterial phase of contrast enhancement. Standard convolution FBP was used to reconstruct 140-kVp (protocol A) and 80-kVp (protocol B) image sets, and ASIR (protocol C) was used to reconstruct 80-kVp image sets. The mean image noise; contrast-to-noise ratio (CNR) relative to muscle for the aorta, liver, and pancreas; and effective dose with each protocol were assessed. A figure of merit (FOM) was computed to normalize the image noise and CNR for each protocol to effective dose. Repeated-measures analysis of variance with Bonferroni adjustment for multiple comparisons was used to compare differences in mean CNR, image noise, and corresponding FOM among the three protocols. The noise power spectra generated from a custom phantom with each protocol were also compared. RESULTS When image noise was normalized to effective dose, protocol C, as compared with protocols A (P = .0002) and B (P = .0001), yielded an approximately twofold reduction in noise. When the CNR was normalized to effective dose, protocol C yielded significantly higher CNRs for the aorta, liver, and pancreas than did protocol A (P = .0001 for all comparisons) and a significantly higher CNR for the liver than did protocol B (P = .003). Mean effective doses were 17.5 mSv +/- 0.6 (standard error) with protocol A and 5.1 mSv +/- 0.3 with protocols B and C. Compared with protocols A and B, protocol C yielded a small but quantifiable noise reduction across the entire spectrum of spatial frequencies. CONCLUSION Compared with standard FBP reconstruction, an ASIR algorithm improves image quality and has the potential to decrease radiation dose at low-tube-voltage, high-tube-current multidetector abdominal CT during the late hepatic arterial phase.

An experimental comparison of detector performance for direct and indirect digital radiography systems

Current flat-panel detectors either directly convert x-ray energy to electronic charge or use indirect conversion with an intermediate optical process. The purpose of this work was to compare direct and indirect detectors in terms of their modulation transfer function (MTF), noise power spectrum (NPS), and detective quantum efficiency (DQE). Measurements were made on three flat-panel detectors, Hologic Direct-Ray DR-1000 (DRC), GE Revolution XQ/i (XQ/i), and Philips Digital Diagnost (DiDi) using the IEC-defined RQA5 (approximately 74 kVp, 21 mm Al) and RQA9 (approximately 120 kVp, 40 mm Al) radiographic techniques. The presampled MTFs of the systems were measured using an edge method [Samei et al., Med. Phys. 25, 102 (1998)]. The NPS of the systems were determined for a range of exposure levels by two-dimensional (2D) Fourier analysis of uniformly exposed radiographs [Flynn and Samei, Med. Phys. 26, 1612 (1999)]. The DQEs were assessed from the measured MTF, NPS, exposure, and estimated ideal signal-to-noise ratios. For the direct system, the MTF was found to be significantly higher than that for the indirect systems and very close to an ideal function associated with the detector pixel size. The NPS for the direct system was found to be constant in relation to frequency. For the XQ/i and DRC systems, the DQE results reflected expected differences based on the absorption efficiency of the different detector materials. Using RQA5, the measured DQE values in the diagonal (and axial) direction(s) at spatial frequencies of 0.15 mm(-1) and 2.5 mm(-1) were 64% (64%) and 20% (15%) for the XQ/i system, and 38% (38%) and 20% (20%) for the DRC, respectively. The DQE results of the DiDi system were difficult to interpret due to additional preprocessing steps in that system.

Towards task-based assessment of CT performance: System and object MTF across different reconstruction algorithms

PURPOSE To investigate a measurement method for evaluating the resolution properties of CT imaging systems across reconstruction algorithms, dose, and contrast. METHODS An algorithm was developed to extract the task-based modulation transfer function (MTF) from disk images generated from the rod inserts in the ACR phantom (model 464 Gammex, WI). These inserts are conventionally employed for HU accuracy assessment. The edge of the disk objects was analyzed to determine the edge-spread function, which was differentiated to yield the line-spread function and Fourier-transformed to generate the object-specific MTF for task-based assessment, denoted MTFTask . The proposed MTF measurement method was validated against the conventional wire technique and further applied to measure the MTF of CT images reconstructed with an adaptive statistical iterative algorithm (ASIR) and a model-based iterative (MBIR) algorithm. Results were further compared to the standard filtered back projection (FBP) algorithm. Measurements were performed and compared across different doses and contrast levels to ascertain the MTFTask dependencies on those factors. RESULTS For the FBP reconstructed images, the MTFTask measured with the inserts were the same as the MTF measured from the wire-based method. For the ASIR and MBIR data, the MTFTask using the high contrast insert was similar to the wire-based MTF and equal or superior to that of FBP. However, results for the MTFTask measured using the low-contrast inserts, the MTFTask for ASIR and MBIR data was lower than for the FBP, which was constant throughout all measurements. Similarly, as a function of mA, the MTFTask for ASIR and MBIR varied as a function of noise--with MTFTask being proportional to mA. Overall greater variability of MTFTask across dose and contrast was observed for MBIR than for ASIR. CONCLUSIONS This approach provides a method for assessing the task-based MTF of a CT system using conventional and iterative reconstructions. Results demonstrated that the object-specific MTF can vary as a function of dose and contrast. The analysis highlighted the paradigm shift for iterative reconstructions when compared to FBP, where iterative reconstructions generally offer superior noise performance but with varying resolution as a function of dose and contrast. The MTFTask generated by this method is expected to provide a more comprehensive assessment of image resolution across different reconstruction algorithms and imaging tasks.PURPOSE To investigate a measurement method for evaluating the resolution properties of CT imaging systems across reconstruction algorithms, dose, and contrast. METHODS An algorithm was developed to extract the task-based modulation transfer function (MTF) from disk images generated from the rod inserts in the ACR phantom (model 464 Gammex, WI). These inserts are conventionally employed for HU accuracy assessment. The edge of the disk objects was analyzed to determine the edge-spread function, which was differentiated to yield the line-spread function and Fourier-transformed to generate the object-specific MTF for task-based assessment, denoted MTF(Task). The proposed MTF measurement method was validated against the conventional wire technique and further applied to measure the MTF of CT images reconstructed with an adaptive statistical iterative algorithm (ASIR) and a model-based iterative (MBIR) algorithm. Results were further compared to the standard filtered back projection (FBP) algorithm. Measurements were performed and compared across different doses and contrast levels to ascertain the MTF(Task) dependencies on those factors. RESULTS For the FBP reconstructed images, the MTF(Task) measured with the inserts were the same as the MTF measured from the wire-based method. For the ASIR and MBIR data, the MTF(Task) using the high contrast insert was similar to the wire-based MTF and equal or superior to that of FBP. However, results for the MTF(Task) measured using the low-contrast inserts, the MTF(Task) for ASIR and MBIR data was lower than for the FBP, which was constant throughout all measurements. Similarly, as a function of mA, the MTF(Task) for ASIR and MBIR varied as a function of noise--with MTF(Task) being proportional to mA. Overall greater variability of MTF(Task) across dose and contrast was observed for MBIR than for ASIR. CONCLUSIONS This approach provides a method for assessing the task-based MTF of a CT system using conventional and iterative reconstructions. Results demonstrated that the object-specific MTF can vary as a function of dose and contrast. The analysis highlighted the paradigm shift for iterative reconstructions when compared to FBP, where iterative reconstructions generally offer superior noise performance but with varying resolution as a function of dose and contrast. The MTF(Task) generated by this method is expected to provide a more comprehensive assessment of image resolution across different reconstruction algorithms and imaging tasks.

Intercomparison of methods for image quality characterization. II. Noise power spectrum

Second in a two-part series comparing measurement techniques for the assessment of basic image quality metrics in digital radiography, in this paper we focus on the measurement of the image noise power spectrum (NPS). Three methods were considered: (1) a method published by Dobbins et al. [Med. Phys. 22, 1581-1593 (1995)], (2) a method published by Samei et al. [Med. Phys. 30, 608-622 (2003)], and (3) a new method sanctioned by the International Electrotechnical Commission (IEC 62220-1, 2003), developed as part of an international standard for the measurement of detective quantum efficiency. In addition to an overall comparison of the estimated NPS between the three techniques, the following factors were also evaluated for their effect on the measured NPS: horizontal versus vertical directional dependence, the use of beam-limiting apertures, beam spectrum, and computational methods of NPS analysis, including the region-of-interest (ROI) size and the method of ROI normalization. Of these factors, none was found to demonstrate a substantial impact on the amplitude of the NPS estimates (< or = 3.1% relative difference in NPS averaged over frequency, for each factor considered separately). Overall, the three methods agreed to within 1.6% +/- 0.8% when averaged over frequencies > 0.15 mm(-1).

Intercomparison of methods for image quality characterization. I. Modulation transfer function.

The modulation transfer function (MTF) and the noise power spectrum (NPS) are widely recognized as the most relevant metrics of resolution and noise performance in radiographic imaging. These quantities have commonly been measured using various techniques, the specifics of which can have a bearing on the accuracy of the results. As a part of a study aimed at comparing the relative performance of different techniques, in this paper we report on a comparison of two established MTF measurement techniques: one using a slit test device [Dobbins et al., Med. Phys. 22, 1581-1593 (1995)] and another using a translucent edge test device [Samei et al., Med. Phys. 25, 102-113 (1998)], with one another and with a third technique using an opaque edge test device recommended by a new international standard (IEC 62220-1, 2003). The study further aimed to substantiate the influence of various acquisition and processing parameters on the estimated MTF. The slit test device was made of 2 mm thick Pb slabs with a 12.5 microm opening. The translucent edge test device was made of a laminated and polished Pt(0.9)Ir(0.1). alloy foil of 0.1 mm thickness. The opaque edge test device was made of a 2 mm thick W slab. All test devices were imaged on a representative indirect flat-panel digital radiographic system using three published beam qualities: 70 kV with 0.5 mm Cu filtration, 70 kV with 19 mm Al filtration, and 74 kV with 21 mm Al filtration (IEC-RQA5). The latter technique was also evaluated in conjunction with two external beam-limiting apertures (per IEC 62220-1), and with the tube collimator limiting the beam to the same area achieved with the apertures. The presampled MTFs were deduced from the acquired images by Fourier analysis techniques, and the results analyzed for relative values and the influence of impacting parameters. The findings indicated that the measurement technique has a notable impact on the resulting MTF estimate, with estimates from the overall IEC method 4.0% +/- 0.2% lower than that of Dobbins et al. and 0.7% +/- 0.4% higher than that of Samei et al. averaged over the zero to cutoff frequency range. Over the same frequency range, keeping beam quality and limitation constant, the average MTF estimate obtained with the edge techniques differed by up to 5.2% +/- 0.2% from that of the slit, with the opaque edge providing lower MTF estimates at lower frequencies than those obtained with the translucent edge or slit. The beam quality impacted the average estimated MTF by as much as 3.7% +/- 0.9% while the use of beam limiting devices alone increased the average estimated MTF by as much as 7.0% +/- 0.9%. While the slit method is inherently very sensitive to misalignment, both edge techniques were found to tolerate misalignments by as much as 6 cm. The results suggest the use of the opaque edge test device and the tube internal collimator for beam limitation in order to achieve an MTF result most reflective of the overall performance of the imaging system and least susceptible to misalignment and scattered radiation. Careful attention to influencing factors is warranted to achieve accurate results.

An experimental comparison of detector performance for computed radiography systems

The intrinsic resolution, noise, and signal-to-noise transfer characteristics of five commercial digital computed radiography (CR) systems were compared using identical experimental methods. The reader/screen combinations evaluated were Agfa ADC-Compact/MD-10, Agfa ADC-Compact/MD-30, Agfa ADC-Solo/MD-10, Agfa ADC-Solo/MD-30, Lumisys CR-2000/MD-10, Fuji FCR-9501 (HQ)/ST-Va, Kodak CR-400/GP-25, and Kodak CR-400/HR. Measurements were made at 70 and 115 kVp with 19 mm added aluminum filtration. The presampled modulation transfer functions (MTFs) of the systems were measured using an edge method. The noise power spectra (NPS) were determined by 2D Fourier analysis of uniformly exposed radiographs. The frequency-dependent detective quantum efficiencies (DQEs) were computed from the MTF, NPS, exposure measurements, and computational estimates of the ideal signal-to-noise ratios. Using 70 kVp and 0.1-0.12 mm pixel sizes, spatial frequencies of 2.1, 2.0, 2.2, 1.9, 2.0, 2.0, 2.3, 2.3, and 3.5 cycles/mm were measured at 0.2 MTF for the eight reader/screen combinations, respectively. Using 70 kVp, 7.74 x 10(-8) C/kg (0.3 mR), and 0.1-0.12 mm pixel sizes, DQE(0.15) values of 20.3%, 22.9%, 24.6%, 28.6%, 22.2%, 30.0%, 29.5%, and 17.3% were obtained for the eight combinations, respectively. The corresponding values at 115 kVp were 15.9%, 18.5%, 21.5%, 21.8%, 15.3%, 23.1%, 22.3%, and 13.8%, respectively. The findings of the study demonstrate the pixel size, orientation, beam quality, screen, and reader dependencies of image quality in CR systems. The physical performance of the systems having standard-resolution screens demonstrated similar resolution performance but more notable variations in DQE. The one high-resolution screen tested had reduced DQE and increased MTF at high frequencies.

Experimental comparison of noise and resolution for 2k and 4k storage phosphor radiography systems.

The purpose of this study was to compare the image quality for a digital storage phosphor system using 1760 x 2140 (2k) and 3520 x 4280 (4k) image arrays. Measurements were made on a chest radiography system (Fuji FCR-9501) with special provisions to be operated in both 2k (standard) and 4k (HQ) modes. Presampled modulation transfer functions (MTF) were measured using an edge method. Noise power spectra (NPS) were determined for different input exposures by two-dimensional Fourier analysis. These measures along with exposure measurements and an x-ray spectral model were used to determine the frequency-dependent detective quantum efficiency DQE (f) of the system for the 4k and the 2k modes. The magnitude of the NPS for the 4k mode was about 1/2 that of the 2k mode. A MTF value of 0.5 was found at 1.25 cycles/mm for the 4k system and 1.50 cycles/mm for the 2k system. The 4k images had an extended MTF of 0.1 at 4.5 cycles/mm in the plate-scan direction. Overall, the DQE (f) of the 4k mode was slightly better than that for the 2k mode by about 0.02 due primarily to its better noise characteristics.

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.

Detection of Pancreatic Tumors, Image Quality, and Radiation Dose during the Pancreatic Parenchymal Phase: Effect of a Low-Tube-Voltage, High-Tube-Current CT Technique—Preliminary Results

PURPOSE To intraindividually compare a low-tube-voltage (80 kVp), high-tube-current (675 mA) computed tomographic (CT) technique with a high-tube-voltage (140 kVp) CT protocol for the detection of pancreatic tumors, image quality, and radiation dose during the pancreatic parenchymal phase. MATERIALS AND METHODS This prospective, single-center, HIPAA-compliant study was approved by the institutional review board, and written informed consent was obtained. Twenty-seven patients (nine men, 18 women; mean age, 64 years) with 23 solitary pancreatic tumors underwent dual-energy CT. Two imaging protocols were used: 140 kVp and 385 mA (protocol A) and 80 kVp and 675 mA (protocol B). For both protocols, the following variables were compared during the pancreatic parenchymal phase: contrast enhancement for the aorta, the pancreas, and the portal vein; pancreas-to-tumor contrast-to-noise ratio (CNR); noise; and effective dose. Two blinded, independent readers qualitatively scored the two data sets for tumor detection and image quality. Random-effect analysis of variance tests were used to compare differences between the two protocols. RESULTS Compared with protocol A, protocol B yielded significantly higher contrast enhancement for the aorta (508.6 HU vs 221.5 HU, respectively), pancreas (151.2 HU vs 67.0 HU), and portal vein (189.7 HU vs 87.3 HU), along with a greater pancreas-to-tumor CNR (8.1 vs 5.9) (P < .001 for all comparisons). No statistically significant difference in tumor detection was observed between the two protocols. Although standard deviation of image noise increased with protocol B (11.5 HU vs 18.6 HU), this protocol significantly reduced the effective dose (from 18.5 to 5.1 mSv; P < .001). CONCLUSION A low-tube-voltage, high-tube-current CT technique has the potential to improve the enhancement of the pancreas and peripancreatic vasculature, improve tumor conspicuity, and reduce patient radiation dose during the pancreatic parenchymal phase.