Past AAPM President: 2021

The following article was published on Wiley Online Library on 07 December 2022 before its intended inclusion in the 50th Anniversary of the Medical Physics Journal special issue. The article has been temporarily removed and will be republished as part of the special issue. Wiley would like to apologize to the author(s) and the academic community for this mistake.

Population of 224 realistic human subject-based computational breast phantoms

PURPOSE

To create a database of highly realistic and anatomically variable 3D virtual breast phantoms based on dedicated breast computed tomography (bCT) data.

METHODS

A tissue classification and segmentation algorithm was used to create realistic and detailed 3D computational breast phantoms based on 230 + dedicated bCT datasets from normal human subjects. The breast volume was identified using a coarse three-class fuzzy C-means segmentation algorithm which accounted for and removed motion blur at the breast periphery. Noise in the bCT data was reduced through application of a postreconstruction 3D bilateral filter. A 3D adipose nonuniformity (bias field) correction was then applied followed by glandular segmentation using a 3D bias-corrected fuzzy C-means algorithm. Multiple tissue classes were defined including skin, adipose, and several fractional glandular densities. Following segmentation, a skin mask was produced which preserved the interdigitated skin, adipose, and glandular boundaries of the skin interior. Finally, surface modeling was used to produce digital phantoms with methods complementary to the XCAT suite of digital human phantoms.

RESULTS

After rejecting some datasets due to artifacts, 224 virtual breast phantoms were created which emulate the complex breast parenchyma of actual human subjects. The volume breast density (with skin) ranged from 5.5% to 66.3% with a mean value of 25.3% ± 13.2%. Breast volumes ranged from 25.0 to 2099.6 ml with a mean value of 716.3 ± 386.5 ml. Three breast phantoms were selected for imaging with digital compression (using finite element modeling) and simple ray-tracing, and the results show promise in their potential to produce realistic simulated mammograms.

CONCLUSIONS

This work provides a new population of 224 breast phantoms based on in vivo bCT data for imaging research. Compared to previous studies based on only a few prototype cases, this dataset provides a rich source of new cases spanning a wide range of breast types, volumes, densities, and parenchymal patterns.

Multi-Institutional Evaluation of Digital Tomosynthesis, Dual-Energy Radiography, and Conventional Chest Radiography for the Detection and Management of Pulmonary Nodules

Purpose To conduct a multi-institutional, multireader study to compare the performance of digital tomosynthesis, dual-energy (DE) imaging, and conventional chest radiography for pulmonary nodule detection and management. Materials and Methods In this binational, institutional review board-approved, HIPAA-compliant prospective study, 158 subjects (43 subjects with normal findings) were enrolled at four institutions. Informed consent was obtained prior to enrollment. Subjects underwent chest computed tomography (CT) and imaging with conventional chest radiography (posteroanterior and lateral), DE imaging, and tomosynthesis with a flat-panel imaging device. Three experienced thoracic radiologists identified true locations of nodules (n = 516, 3-20-mm diameters) with CT and recommended case management by using Fleischner Society guidelines. Five other radiologists marked nodules and indicated case management by using images from conventional chest radiography, conventional chest radiography plus DE imaging, tomosynthesis, and tomosynthesis plus DE imaging. Sensitivity, specificity, and overall accuracy were measured by using the free-response receiver operating characteristic method and the receiver operating characteristic method for nodule detection and case management, respectively. Results were further analyzed according to nodule diameter categories (3-4 mm, >4 mm to 6 mm, >6 mm to 8 mm, and >8 mm to 20 mm). Results Maximum lesion localization fraction was higher for tomosynthesis than for conventional chest radiography in all nodule size categories (3.55-fold for all nodules, P < .001; 95% confidence interval [CI]: 2.96, 4.15). Case-level sensitivity was higher with tomosynthesis than with conventional chest radiography for all nodules (1.49-fold, P < .001; 95% CI: 1.25, 1.73). Case management decisions showed better overall accuracy with tomosynthesis than with conventional chest radiography, as given by the area under the receiver operating characteristic curve (1.23-fold, P < .001; 95% CI: 1.15, 1.32). There were no differences in any specificity measures. DE imaging did not significantly affect nodule detection when paired with either conventional chest radiography or tomosynthesis. Conclusion Tomosynthesis outperformed conventional chest radiography for lung nodule detection and determination of case management; DE imaging did not show significant differences over conventional chest radiography or tomosynthesis alone. These findings indicate performance likely achievable with a range of reader expertise. ^© RSNA, 2016 Online supplemental material is available for this article.

Development and application of a suite of 4-D virtual breast phantoms for optimization and evaluation of breast imaging systems

Mammography is currently the most widely utilized tool for detection and diagnosis of breast cancer. However, in women with dense breast tissue, tissue overlap may obscure lesions. Digital breast tomosynthesis can reduce tissue overlap. Furthermore, imaging with contrast enhancement can provide additional functional information about lesions, such as morphology and kinetics, which in turn may improve lesion identification and characterization. The performance of these imaging techniques is strongly dependent on the structural composition of the breast, which varies significantly among patients. Therefore, imaging system and imaging technique optimization should take patient variability into consideration. Furthermore, optimization of imaging techniques that employ contrast agents should include the temporally varying breast composition with respect to the contrast agent uptake kinetics. To these ends, we have developed a suite of 4-D virtual breast phantoms, which are incorporated with the kinetics of contrast agent propagation in different tissues and can realistically model normal breast parenchyma as well as benign and malignant lesions. This development presents a new approach in performing simulation studies using truly anthropomorphic models. To demonstrate the utility of the proposed 4-D phantoms, we present a simplified example study to compare the performance of 14 imaging paradigms qualitatively and quantitatively.

The effect of averaging adjacent planes for artifact reduction in matrix inversion tomosynthesis

PURPOSE

Matrix inversion tomosynthesis (MITS) uses linear systems theory and knowledge of the imaging geometry to remove tomographic blur that is present in conventional backprojection tomosynthesis reconstructions, leaving in-plane detail rendered clearly. The use of partial-pixel interpolation during the backprojection process introduces imprecision in the MITS modeling of tomographic blur, and creates low-contrast artifacts in some MITS planes. This paper examines the use of MITS slabs, created by averaging several adjacent MITS planes, as a method for suppressing partial-pixel artifacts.

METHODS

Human chest tomosynthesis projection data, acquired as part of an IRB-approved pilot study, were used to generate MITS planes, three-plane MITS slabs (MITSa3), five-plane MITS slabs (MITSa5), and seven-plane MITS slabs (MITSa7). These were qualitatively examined for partial-pixel artifacts and the visibility of normal and abnormal anatomy. Additionally, small (5 mm) subtle pulmonary nodules were simulated and digitally superimposed upon human chest tomosynthesis projection images, and their visibility was qualitatively assessed in the different reconstruction techniques. Simulated images of a thin wire were used to generate modulation transfer function (MTF) and slice-sensitivity profile curves for the different MITS and MITS slab techniques, and these were examined for indications of partial-pixel artifacts and frequency response uniformity. Finally, mean-subtracted, exposure-normalized noise power spectra (ENNPS) estimates were computed and compared for MITS and MITS slab reconstructions, generated from 10 sets of tomosynthesis projection data of an acrylic slab. The simulated in-plane MTF response of each technique was also combined with the square root of the ENNPS estimate to yield stochastic signal-to-noise ratio (SNR) information about the different reconstruction techniques.

RESULTS

For scan angles of 20° and 5 mm plane separation, seven MITS planes must be averaged to sufficiently remove partial-pixel artifacts. MITSa7 does appear to subtly reduce the contrast of high-frequency "edge" information, but the removal of partial-pixel artifacts makes the appearance of low-contrast, fine-detail anatomy even more conspicuous in MITSa7 slices. MITSa7 also appears to render simulated subtle 5 mm pulmonary nodules with greater visibility than MITS alone, in both the open lung and regions overlying the mediastinum. Finally, the MITSa7 technique reduces stochastic image variance, though the in-plane stochastic SNR (for very thin objects which do not span multiple MITS planes) is only improved at spatial frequencies between 0.05 and 0.20 cycles∕mm.

CONCLUSIONS

The MITSa7 method is an improvement over traditional single-plane MITS for thoracic imaging and the pulmonary nodule detection task, and thus the authors plan to use the MITSa7 approach for all future MITS research at the authors' institution.

Generation of a suite of 3D computer-generated breast phantoms from a limited set of human subject data

PURPOSE

The authors previously reported on a three-dimensional computer-generated breast phantom, based on empirical human image data, including a realistic finite-element based compression model that was capable of simulating multimodality imaging data. The computerized breast phantoms are a hybrid of two phantom generation techniques, combining empirical breast CT (bCT) data with flexible computer graphics techniques. However, to date, these phantoms have been based on single human subjects. In this paper, the authors report on a new method to generate multiple phantoms, simulating additional subjects from the limited set of original dedicated breast CT data. The authors developed an image morphing technique to construct new phantoms by gradually transitioning between two human subject datasets, with the potential to generate hundreds of additional pseudoindependent phantoms from the limited bCT cases. The authors conducted a preliminary subjective assessment with a limited number of observers (n = 4) to illustrate how realistic the simulated images generated with the pseudoindependent phantoms appeared.

METHODS

Several mesh-based geometric transformations were developed to generate distorted breast datasets from the original human subject data. Segmented bCT data from two different human subjects were used as the "base" and "target" for morphing. Several combinations of transformations were applied to morph between the "base' and "target" datasets such as changing the breast shape, rotating the glandular data, and changing the distribution of the glandular tissue. Following the morphing, regions of skin and fat were assigned to the morphed dataset in order to appropriately assign mechanical properties during the compression simulation. The resulting morphed breast was compressed using a finite element algorithm and simulated mammograms were generated using techniques described previously. Sixty-two simulated mammograms, generated from morphing three human subject datasets, were used in a preliminary observer evaluation where four board certified breast radiologists with varying amounts of experience ranked the level of realism (from 1 = "fake" to 10 = "real") of the simulated images.

RESULTS

The morphing technique was able to successfully generate new and unique morphed datasets from the original human subject data. The radiologists evaluated the realism of simulated mammograms generated from the morphed and unmorphed human subject datasets and scored the realism with an average ranking of 5.87 ± 1.99, confirming that overall the phantom image datasets appeared more "real" than "fake." Moreover, there was not a significant difference (p > 0.1) between the realism of the unmorphed datasets (6.0 ± 1.95) compared to the morphed datasets (5.86 ± 1.99). Three of the four observers had overall average rankings of 6.89 ± 0.89, 6.9 ± 1.24, 6.76 ± 1.22, whereas the fourth observer ranked them noticeably lower at 2.94 ± 0.7.

CONCLUSIONS

This work presents a technique that can be used to generate a suite of realistic computerized breast phantoms from a limited number of human subjects. This suite of flexible breast phantoms can be used for multimodality imaging research to provide a known truth while concurrently producing realistic simulated imaging data.

Estimation of the two-dimensional presampled modulation transfer function of digital radiography devices using one-dimensional test objects

PURPOSE

The modulation transfer function (MTF) of medical imaging devices is commonly reported in the form of orthogonal one-dimensional (1D) measurements made near the vertical and horizontal axes with a slit or edge test device. A more complete description is found by measuring the two-dimensional (2D) MTF. Some 2D test devices have been proposed, but there are some issues associated with their use: (1) they are not generally available; (2) they may require many images; (3) the results may have diminished accuracy; and (4) their implementation may be particularly cumbersome. This current work proposes the application of commonly available 1D test devices for practical and accurate estimation of the 2D presampled MTF of digital imaging systems.

METHODS

Theory was developed and applied to ensure adequate fine sampling of the system line spread function for 1D test devices at orientations other than approximately vertical and horizontal. Methods were also derived and tested for slit nonuniformity correction at arbitrary angle. Techniques were validated with experimental measurements at ten angles using an edge test object and three angles using a slit test device on an indirect-detection flat-panel system [GE Revolution XQ∕i (GE Healthcare, Waukesha, WI)]. The 2D MTF was estimated through a simple surface fit with interpolation based on Delaunay triangulation of the 1D edge-based MTF measurements. Validation by synthesis was also performed with simulated images from a hypothetical direct-detection flat-panel device.

RESULTS

The 2D MTF derived from physical measurements yielded an average relative precision error of 0.26% for frequencies below the cutoff (2.5 mm(-1)) and approximate circular symmetry at frequencies below 4 mm(-1). While slit analysis generally agreed with the results of edge analysis, the two showed subtle differences at frequencies above 4 mm(-1). Slit measurement near 45° revealed radial asymmetry in the MTF resulting from the square pixel aperture (0.2 mm × 0.2 mm), a characteristic which was not necessarily appreciated with the orthogonal 1D MTF measurements. In simulation experiments, both slit- and edge-based measurements resolved the radial asymmetries in the 2D MTF. The average absolute relative accuracy error in the 2D MTF between the DC and cutoff (2.5 mm(-1)) frequencies was 0.13% with average relative precision error of 0.11%. Other simulation results were similar to those derived from physical data.

CONCLUSIONS

Overall, the general availability, acceptance, accuracy, and ease of implementation of 1D test devices for MTF assessment make this a valuable technique for 2D MTF estimation.

Automated lung segmentation in digital chest tomosynthesis

PURPOSE

The purpose of this study was to develop an automated lung segmentation method for computerized detection of lung nodules in digital chest tomosynthesis.

METHODS

The authors collected 45 digital tomosynthesis scans and manually segmented reference lung regions in each scan to assess the performance of the method. The authors automated the technique by calculating the edge gradient in an original image for enhancing lung outline and transforming the edge gradient image to polar coordinate space. The authors then employed a dynamic programming technique to delineate outlines of the unobscured lungs in the transformed edge gradient image. The lung outlines were converted back to the original image to provide the final segmentation result. The above lung segmentation algorithm was first applied to the central reconstructed tomosynthesis slice because of the absence of ribs overlapping lung structures. The segmented lung in the central slice was then used to guide lung segmentation in noncentral slices. The authors evaluated the segmentation method by using (1) an overlap rate of lung regions, (2) a mean absolute distance (MAD) of lung borders, (3) a Hausdorff distance of lung borders between the automatically segmented lungs and manually segmented reference lungs, and (4) the fraction of nodules included in the automatically segmented lungs.

RESULTS

The segmentation method achieved mean overlap rates of 85.7%, 88.3%, and 87.0% for left lungs, right lungs, and entire lungs, respectively; mean MAD of 4.8, 3.9, and 4.4 mm for left lungs, right lungs, and entire lungs, respectively; and mean Hausdorrf distance of 25.0 mm, 25.5 mm, and 30.1 mm for left lungs, right lungs, and entire lungs, respectively. All of the nodules inside the reference lungs were correctly included in the segmented lungs obtained with the lung segmentation method.

CONCLUSIONS

The method achieved relatively high accuracy for lung segmentation and will be useful for computer-aided detection of lung nodules in digital tomosynthesis.

An analysis of the mechanical parameters used for finite element compression of a high-resolution 3D breast phantom

PURPOSE

The authors previously introduced a methodology to generate a realistic three-dimensional (3D), high-resolution, computer-simulated breast phantom based on empirical data. One of the key components of such a phantom is that it provides a means to produce a realistic simulation of clinical breast compression. In the current study, they have evaluated a finite element (FE) model of compression and have demonstrated the effect of a variety of mechanical properties on the model using a dense mesh generated from empirical breast data. While several groups have demonstrated an effective compression simulation with lower density finite element meshes, the presented study offers a mesh density that is able to model the morphology of the inner breast structures more realistically than lower density meshes. This approach may prove beneficial for multimodality breast imaging research, since it provides a high level of anatomical detail throughout the simulation study.

METHODS

In this paper, the authors describe methods to improve the high-resolution performance of a FE compression model. In order to create the compressible breast phantom, dedicated breast CT data was segmented and a mesh was generated with 4-noded tetrahedral elements. Using an explicit FE solver to simulate breast compression, several properties were analyzed to evaluate their effect on the compression model including: mesh density, element type, density, and stiffness of various tissue types, friction between the skin and the compression plates, and breast density. Following compression, a simulated projection was generated to demonstrate the ability of the compressible breast phantom to produce realistic simulated mammographic images.

RESULTS

Small alterations in the properties of the breast model can change the final distribution of the tissue under compression by more than 1 cm; which ultimately results in different representations of the breast model in the simulated images. The model properties that impact displacement the most are mesh density, friction between the skin and the plates, and the relative stiffness of the different tissue types.

CONCLUSIONS

The authors have developed a 3D, FE breast model that can yield high spatial resolution breast deformations under uniaxial compression for imaging research purposes and demonstrated that small changes in the mechanical properties can affect images generated using the phantom.

Chest tomosynthesis: technical principles and clinical update

Digital tomosynthesis is a radiographic technique that can produce an arbitrary number of section images of a patient from a single pass of the X-ray tube. It utilizes a conventional X-ray tube, a flat-panel detector, a computer-controlled tube mover, and special reconstruction algorithms to produce section images. While it does not have the depth resolution of computed tomography (CT), tomosynthesis provides some of the tomographic benefits of CT but at lower cost and radiation dose than CT. Compared to conventional chest radiography, chest tomosynthesis results in improved visibility of normal structures such as vessels, airway and spine. By reducing visual clutter from overlying normal anatomy, it also enhances detection of small lung nodules. This review article outlines the components of a tomosynthesis system, discusses results regarding improved lung nodule detection from the recent literature, and presents examples of nodule detection from a clinical trial in human subjects. Possible implementation strategies for use in clinical chest imaging are discussed.

Effective DQE (eDQE) and speed of digital radiographic systems: an experimental methodology

Prior studies on performance evaluation of digital radiographic systems have primarily focused on the assessment of the detector performance alone. However, the clinical performance of such systems is also substantially impacted by magnification, focal spot blur, the presence of scattered radiation, and the presence of an antiscatter grid. The purpose of this study is to evaluate an experimental methodology to assess the performance of a digital radiographic system, including those attributes, and to propose a new metric, effective detective quantum efficiency (eDQE), a candidate for defining the efficiency or speed of digital radiographic imaging systems. The study employed a geometric phantom simulating the attenuation and scatter properties of the adult human thorax and a representative indirect flat-panel-based clinical digital radiographic imaging system. The noise power spectrum (NPS) was derived from images of the phantom acquired at three exposure levels spanning the operating range of the clinical system. The modulation transfer function (MTF) was measured using an edge device positioned at the surface of the phantom, facing the x-ray source. Scatter measurements were made using a beam stop technique. The eDQE was then computed from these measurements, along with measures of phantom attenuation and x-ray flux. The MTF results showed notable impact from the focal spot blur, while the NPS depicted a large component of structured noise resulting from use of an antiscatter grid. The eDQE was found to be an order of magnitude lower than the conventional DQE. At 120 kVp, eDQE(0) was in the 8%-9% range, fivefold lower than DQE(0) at the same technique. The eDQE method yielded reproducible estimates of the system performance in a clinically relevant context by quantifying the inherent speed of the system, that is, the actual signal to noise ratio that would be measured under clinical operating conditions.

Methodology for generating a 3D computerized breast phantom from empirical data

The initial process for creating a flexible three-dimensional computer-generated breast phantom based on empirical data is described. Dedicated breast computed-tomography data were processed to suppress noise and scatter artifacts in the reconstructed image set. An automated algorithm was developed to classify the breast into its primary components. A preliminary phantom defined using subdivision surfaces was generated from the segmented data. To demonstrate potential applications of the phantom, simulated mammographic image data were acquired of the phantom using a simplistic compression model and an analytic projection algorithm directly on the surface model. The simulated image was generated using a model for a polyenergetic cone-beam projection of the compressed phantom. The methods used to create the breast phantom generate resulting images that have a high level of tissue structure detail available and appear similar to actual mammograms. Fractal dimension measurements of simulated images of the phantom are comparatively similar to measurements from images of real human subjects. A realistic and geometrically defined breast phantom that can accurately simulate imaging data may have many applications in breast imaging research.

Tomosynthesis imaging: at a translational crossroads

Tomosynthesis is a decades-old technique for section imaging that has seen a recent upsurge in interest due to its promise to provide three-dimensional information at lower dose and potentially lower cost than CT in certain clinical imaging situations. This renewed interest in tomosynthesis began in the late 1990s as a new generation of flat-panel detectors became available; these detectors were the one missing piece of the picture that had kept tomosynthesis from enjoying significant utilization earlier. In the past decade, tomosynthesis imaging has been investigated in a variety of clinical imaging situations, but the two most prominent have been in breast and chest imaging. Tomosynthesis has the potential to substantially change the way in which breast cancer and pulmonary nodules are detected and managed. Commercial tomosynthesis devices are now available or on the horizon. Many of the remaining research activities with tomosynthesis will be translational in nature and will involve physicist and clinician alike. This overview article provides a forward-looking assessment of the translational questions facing tomosynthesis imaging and anticipates some of the likely research and clinical activities in the next five years.

Detector or system? Extending the concept of detective quantum efficiency to characterize the performance of digital radiographic imaging systems

PURPOSE

To develop an experimental method for measuring the effective detective quantum efficiency (eDQE) of digital radiographic imaging systems and evaluate its use in select imaging systems.

MATERIALS AND METHODS

A geometric phantom emulating the attenuation and scatter properties of the adult human thorax was employed to assess eight imaging systems in a total of nine configurations. The noise power spectrum (NPS) was derived from images of the phantom acquired at three exposure levels spanning the operating range of the system. The modulation transfer function (MTF) was measured by using an edge device positioned at the anterior surface of the phantom. Scatter measurements were made by using a beam-stop technique. All measurements, including those of phantom attenuation and estimates of x-ray flux, were used to compute the eDQE.

RESULTS

The MTF results showed notable degradation owing to focal spot blur. Scatter fractions ranged between 11% and 56%, depending on the system. The eDQE(0) results ranged from 1%-17%, indicating a reduction of up to one order of magnitude and different rank ordering and performance among systems, compared with that implied in reported conventional detective quantum efficiency results from the same systems.

CONCLUSION

The eDQE method was easy to implement, yielded reproducible results, and provided a meaningful reflection of system performance by quantifying image quality in a clinically relevant context. The difference in the magnitude of the measured eDQE and the ideal eDQE of 100% provides a great opportunity for improving the image quality of radiographic and mammographic systems while reducing patient dose.

Impulse response and Modulation Transfer Function analysis for Shift-And-Add and Back Projection image reconstruction algorithms in Digital Breast Tomosynthesis (DBT)

Breast cancer is second only to lung cancer as the leading cause of non-preventable cancer death in women. Digital Breast Tomosynthesis (DBT) is a promising technique to improve early breast cancer detection. In this paper, we present the impulse response and Modulation Transfer Function (MTF) analysis to quantitatively compare Shift-And-Add (SAA) and point-by-point Back Projection (BP) three-dimensional image reconstruction algorithms in DBT. A Filtered Back Projection (FBP) deblurring algorithm based on point-by-point BP was used to demonstrate deblurred tomosynthesis images.

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.

Assessment of detective quantum efficiency: intercomparison of a recently introduced international standard with prior methods

PURPOSE

To prospectively evaluate the recently introduced international standard method for measurement of the detective quantum efficiency (DQE) of digital radiography systems, in comparison with representative prior methods.

MATERIALS AND METHODS

A recently introduced international standard method (International Electrotechnical Commission [IEC] 62220-1, 2003) for DQE measurement and two previously described DQE evaluation methods were considered. In addition to an overall comparison, evaluations of the following method factors were performed: beam quality, beam-limiting devices (apertures or collimators), noise power spectrum (NPS) analysis algorithms and parameters (area, region of interest size, background detrending), and modulation transfer function (MTF) test devices and methods.

RESULTS

Overall, at low to middle frequencies, the IEC method yielded DQE estimates that were 3.3% and 6.5% lower than the values yielded by the two previous methods. Averaged over the frequency range of 1.5-2.5 mm(-1), the DQE estimate derived by using the IEC method was 7.1% lower and 12.4% higher than the estimates derived by using the other two methods. Results obtained with the two previous DQE evaluation methods agreed well (within 2.0%) in the low- to middle-frequency range but diverged by up to 10% at higher frequencies. When the DQE method factors were evaluated separately, the largest percentage deviations in DQE were associated with (in order of decreasing influence) the MTF analysis method ( approximately 11%), the beam limitation (about 7%-10%), the beam quality ( approximately 9%), and the NPS analysis method ( approximately 3%).

CONCLUSION

Comparison of DQE estimates obtained by using the recently introduced international standard technique with those obtained by using prior methods revealed that the overall measurement method can affect the DQE estimate by as much as 12%. Findings further suggest that both beam limitation achieved by means of internal collimation (rather than external apertures) and use of a radio-opaque edge MTF device yield a more accurate estimation of the DQE.

Multiprojection correlation imaging for improved detection of pulmonary nodules

OBJECTIVE

The purpose of this study was the development and preliminary evaluation of multiprojection correlation imaging with 3D computer-aided detection (CAD) on chest radiographs for cost- and dose-effective improvement of early detection of pulmonary nodules.

SUBJECTS AND METHODS

Digital chest radiographs of 10 configurations of a chest phantom and of seven human subjects were acquired in multiple angular projections with an acquisition time of 11 seconds (single breath-hold) and total exposure comparable with that of a posteroanterior chest radiograph. An initial 2D CAD algorithm with two difference-of-gaussians filters and multilevel thresholds was developed with an independent database of 44 single-view chest radiographs with confirmed lesions. This 2D CAD algorithm was used on each projection image to find likely suspect nodules. The CAD outputs were reconstructed in 3D, reinforcing signals associated with true nodules while simultaneously decreasing false-positive findings produced by overlapping anatomic features. The performance of correlation imaging was tested on two to 15 projection images.

RESULTS

Optimum performance of correlation imaging was attained when nine projection images were used. Compared with conventional, single-view CAD, correlation imaging decreased as much as 79% the frequency of false-positive findings in phantom cases at a sensitivity level of 65%. The corresponding reduction in false-positive findings in the cases of human subjects was 78%.

CONCLUSION

Although limited by a relatively simple CAD implementation and a small number of cases, the findings suggest that correlation imaging performs substantially better than single-view CAD and may greatly enhance identification of subtle solitary pulmonary nodules on chest radiographs.

Intercomparison of methods for image quality characterization. II. Noise power spectruma)

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 functiona)

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.

Optimization of the matrix inversion tomosynthesis (MITS) impulse response and modulation transfer function characteristics for chest imaging

Matrix inversion tomosynthesis (MITS) uses linear systems theory, along with a priori knowledge of the imaging geometry, to deterministically distinguish between true structure and overlying tomographic blur in a set of conventional tomosynthesis planes. In this paper we examine the effect of total scan angle (ANG), number of input projections (N), and plane separation/number of reconstructed planes (NP) on the MITS impulse response (IR) and modulation transfer function (MTF), with the purpose of optimizing MITS imaging of the chest. MITS IR and MTF data were generated by simulating the imaging of a very thin wire, using various combinations of ANG, N, and NP. Actual tomosynthesis data of an anthropomorphic chest phantom were acquired with a prototype experimental system, using the same imaging parameter combinations as those in the simulations. Thoracic projection data from two human subjects were collected for corroboration of the system response analysis in vivo. Results suggest that ANG=20 degrees, N=71, NP=69 is the optimal combination for MITS chest imaging given the inherent constraints of our prototype system. MITS chest data from human subjects demonstrates that the selected imaging strategy can effectively produce high-quality MITS thoracic images in vivo.

Measurement of the detective quantum efficiency in digital detectors consistent with the IEC 62220-1 standard: Practical considerations regarding the choice of filter material

As part of a larger evaluation we attempted to measure the detective quantum efficiency (DQE) of an amorphous silicon flat-panel detector using the method described in the International Electrotechnical Commission standard 62220-1 published in October 2003. To achieve the radiographic beam conditions specified in the standard, we purchased scientific-grade ultrahigh purity aluminum (99.999% purity, type-11999 alloy) filters in thicknesses ranging from 0.1 through 10.0 mm from a well-known, specialty metals supplier. Qualitative evaluation of flat field images acquired at 71 kV (RQA5 beam quality) with 21 mm of ultrahigh purity aluminum filtration demonstrated a low frequency mottle that was reproducible and was not observed when the measurement was repeated at 74 kV (RQA5 beam quality) with 21 mm of lower-purity aluminum (99.0% purity, type-1100 alloy) filtration. This finding was ultimately attributed to the larger grain size (approximately 1-2 mm) of high purity aluminum metal, which is a well-known characteristic, particularly in thicknesses greater than 1 mm. The impact of this low frequency mottle is to significantly overestimate the noise power spectrum (NPS) at spatial frequencies < or = 0.2 mm(-1), which in turn would cause an underestimation of the DQE in this range. A subsequent evaluation of ultrahigh purity aluminum, purchased from a second source, suggests, that reduced grain size can be achieved by the process of annealing. Images acquired with this sample demonstrated vertical striated nonuniformities that are attributed to the manufacturing method and which do not appear to appreciably impact the NPS at spatial frequencies > or = 0.5 mm(-1), but do result in an asymmetry in the x- and y-NPS at spatial frequencies < or = 0.2 mm(-1). Our observations of markedly visible nonuniformities in images acquired with high purity aluminum filtration suggest that the uniformity of filter materials should be carefully evaluated and taken into consideration when measuring the DQE.

Comparative scatter and dose performance of slot-scan and full-field digital chest radiography systems

PURPOSE

To evaluate the scatter, dose, and effective detective quantum efficiency (DQE) performance of a slot-scan digital chest radiography system compared with that of a full-field digital radiography system.

MATERIALS AND METHODS

Scatter fraction of a slot-scan system was measured for an anthropomorphic and a geometric phantom by using a posterior beam-stop technique at 117 and 140 kVp. Measurements were repeated with a full-field digital radiography system with and without a 13:1 antiscatter grid at 120 and 140 kVp. For both systems, the effective dose was measured on posteroanterior and lateral views for standard clinical techniques by using dosimeters embedded in a female phantom. The effective DQEs of the two systems were assessed by taking into account the scatter performance and the DQE of each system. The statistical significance of all the comparative differences was ascertained by means of t test analysis.

RESULTS

The slot-scan system and the full-field system with grid yielded scatter fractions of 0.13-0.14 and 0.42-0.48 in the lungs and 0.30-0.43 and 0.69-0.78 in the mediastinum, respectively. The sum of the effective doses for posteroanterior and lateral views for the slot-scan system (0.057 mSv +/- 0.003 [+/- standard deviation]) was 34% lower than that for the full-field system (0.086 mSv +/- 0.001, P < .05) at their respective clinical peak voltages (140 and 120 kVp, respectively). The effective DQE of the slot-scan system was equivalent to that of the full-field system in the lung region but was 37% higher in the dense regions (P < .05).

CONCLUSION

The slot-scan design leads to marked scatter reduction compared with the more conventional full-field geometries with a grid. The improved scatter performance of a slot-scan geometry can effectively compensate for low DQE and lead to improved image quality.

A framework for optimising the radiographic technique in digital X-ray imaging

The transition to digital radiology has provided new opportunities for improved image quality, made possible by the superior detective quantum efficiency and post-processing capabilities of new imaging systems, and advanced imaging applications, made possible by rapid digital image acquisition. However, this transition has taken place largely without optimising the radiographic technique used to acquire the images. This paper proposes a framework for optimising the acquisition of digital X-ray images. The proposed approach is based on the signal and noise characteristics of the digital images and the applied exposure. Signal is defined, based on the clinical task involved in an imaging application, as the difference between the detector signal with and without a target present against a representative background. Noise is determined from the noise properties of uniformly acquired images of the background, taking into consideration the absorption properties of the detector. Incident exposure is estimated or otherwise measured free in air, and converted to dose. The main figure of merit (FOM) for optimisation is defined as the signal-difference-to-noise ratio (SdNR) squared per unit exposure or (more preferably) dose. This paper highlights three specific technique optimisation studies that used this approach to optimise the radiographic technique for digital chest and breast applications. In the first study, which was focused on chest radiography with a CsI flat-panel detector, a range of kV(p) (50-150) and filtration (Z = 13-82) were examined in terms of their associated FOM as well as soft tissue to bone contrast, a factor of importance in digital chest radiography. The results indicated that additive Cu filtration can improve image quality. A second study in digital mammography using a selenium direct flat-panel detector indicated improved SdNR per unit exposure with the use of a tungsten target and a rhodium filter than conventional molybdenum target/molybdenum filter techniques. Finally, a third study focusing on cone-beam computed tomography of the breast using a CsI flat-panel detector indicated that high Z filtration of a tungsten target X-ray beam can notably improve the signal and noise characteristics of the image. The general findings highlight the fact that the techniques that are conventionally assumed to be optimum may need to be revisited for digital radiography.

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.

Digital x-ray tomosynthesis: current state of the art and clinical potential

Digital x-ray tomosynthesis is a technique for producing slice images using conventional x-ray systems. It is a refinement of conventional geometric tomography, which has been known since the 1930s. In conventional geometric tomography, the x-ray tube and image receptor move in synchrony on opposite sides of the patient to produce a plane of structures in sharp focus at the plane containing the fulcrum of the motion; all other structures above and below the fulcrum plane are blurred and thus less visible in the resulting image. Tomosynthesis improves upon conventional geometric tomography in that it allows an arbitrary number of in-focus planes to be generated retrospectively from a sequence of projection radiographs that are acquired during a single motion of the x-ray tube. By shifting and adding these projection radiographs, specific planes may be reconstructed. This topical review describes the various reconstruction algorithms used to produce tomosynthesis images, as well as approaches used to minimize the residual blur from out-of-plane structures. Historical background and mathematical details are given for the various approaches described. Approaches for optimizing the tomosynthesis image are given. Applications of tomosynthesis to various clinical tasks, including angiography, chest imaging, mammography, dental imaging and orthopaedic imaging, are also described.

Digital tomosynthesis of hand joints for arthritis assessment

The two principal forms of hand arthritis, rheumatoid arthritis (RA) and osteoarthritis (OA) have large clinical and economic costs. Radiography has been shown to be a useful tool to assess the condition of the disease. A hand radiograph, however, is a two-dimensional projection of a three-dimensional object. In this report we present the results of a study that applied digital tomosynthesis to hand radiography in order to extract three-dimensional outcome measures that should be more sensitive to arthritis progression. The study was performed using simulated projection radiographs created using micro computed tomography (microCT) and a set of five dry-bone hand skeletons. These simulated projection images were then reconstructed into tomographic slices using the matrix inversion tomosynthesis (MITS) algorithm. The accuracy of the tomosynthesis reconstruction was evaluated by comparing the reconstructed images to a gold standard created using the microCT data. A parameter from image registration science, normalized mutual information, provided a quantifiable figure of merit. This study examined the effects of source displacement, number of reconstructed planes, number of acquisitions, noise added to the gray scale images, and errors in the location of a fiducial marker. We also optimized the reconstruction as a function of two variables k and alpha, that controlled the mixing of MITS with conventional shift-and-add tomosynthesis. A study using hand delineated joint margins demonstrated that MITS images provided a better measurement of average joint space width. We found good agreement between the MITS slices and the true planes. Both joint margins and trabecular structure were visible and the reconstructed slices showed additional structures not visible with the standard projection image. Using hand-delineated joint margins we compared the average joint space width of the gold standard slices to the MITS and projection images. A root-mean square deviation (RMSD), calculated for this comparison, gave RMSDproj = 0.18 mm and RMSDMITS = 0.14 mm for the projection and MITS images, respectively. We have demonstrated the potential of digital tomosynthesis for imaging of the hand to assess arthritic changes. We have also developed a methodology that can be used to optimize the technique and have studied the issues that will control the feasibility of clinical implementation.

Quantitative evaluation of noise reduction strategies in dual-energy imaging

In this paper we describe a quantitative evaluation of the performance of three dual-energy noise reduction algorithms: Kalender's correlated noise reduction (KCNR), noise clipping (NOC), and edge-predictive adaptive smoothing (EPAS). These algorithms were compared to a simple smoothing filter approach, using the variance and noise power spectrum measurements of the residual noise in dual-energy images acquired with an a-Si TFT flat-panel x-ray detector. An estimate of the true noise was made through a new method with subpixel accuracy by subtracting an individual image from an ensemble average image. The results indicate that in the lung regions of the tissue image, all three algorithms reduced the noise by similar percentages at high spatial frequencies (KCNR=88%, NOC=88%, EPAS=84%, NOC/KCNR=88%) and somewhat less at low spatial frequencies (KCNR=45%, NOC=54%, EPAS=52%, NOC/KCNR=55%). At low frequencies, the presence of edge artifacts from KCNR made the performance worse, thus NOC or NOC combined with KCNR performed best. At high frequencies, KCNR performed best in the bone image, yet NOC performed best in the tissue image. Noise reduction strategies in dual-energy imaging can be effective and should focus on blending various algorithms depending on anatomical locations.

Chest radiography: optimization of X-ray spectrum for cesium iodide-amorphous silicon flat-panel detector

PURPOSE

To ascertain the optimum x-ray spectrum for chest radiography with a cesium iodide-amorphous silicon flat-panel detector.

MATERIALS AND METHODS

End points for optimization included the ratio of tissue contrast to bone contrast and a figure of merit (FOM) equal to the square of the signal-to-noise ratio of tissue divided by incident exposure to the patient. Studies were conducted with both computer spectrum modeling and experimental measurement in narrow-beam and full-field exposure conditions for four tissue thicknesses (8-32 cm). Three parameters that affect spectra were considered: the atomic number (Z) of filter material (Z = 13, 26, 29, 42, 50, 56, 64, 74, and 82), kilovoltage (from 50 to 150 kVp), and filter thickness (from 0.25 to 2.00 half-value layer [HVL]).

RESULTS

Computer modeling and narrow-beam experimental data showed similar trends for the full range of parameters evaluated. Spectrum model results showed that copper filtration at 120 kVp or more was optimum for FOM. The ratio of contrasts showed a trend to be higher with higher kilovoltage and only a minor variation with filter material. Full-field experimental results, which reflect the added contribution of x-ray scatter, differed in magnitude but not trends from the narrow-beam data in all cases except the ratio of contrasts in the mediastinum.

CONCLUSION

The best performance overall, including both FOM and ratio of contrasts, was at 120 kVp with 1-HVL copper filtration (0.2 mm). With this beam spectrum and an increase in tube output (ie, milliampere seconds) of about 50%, a chest radiograph can be obtained with image quality approximately equal to that with a conventional spectrum but with about 25% less patient exposure.

Imaging characteristics of an amorphous silicon flat-panel detector for digital chest radiography

PURPOSE

To evaluate the imaging characteristics of an amorphous silicon flat-panel detector (FPD) for digital chest radiography.

MATERIALS AND METHODS

The 41 x 41-cm digital FPD is constructed on a single monolithic glass substrate with a structured cesium iodide scintillator layer and an amorphous silicon thin-film transistor array for image readout. Basic imaging characteristics of the FPD and associated image processing system were assessed on acquired images, including linearity, repeatability, uniformity of response, modulation transfer function (MTF), noise power spectrum, detective quantum efficiency (DQE), contrast sensitivity, and scatter content. Results with the FPD system were compared to those with a storage phosphor computed radiography (CR) system.

RESULTS

Images obtained with the FPD demonstrated excellent uniformity, repeatability, and linearity, as well as MTF and DQE that were superior to those with the storage phosphor CR system. The contrast and scatter content of images acquired with the FPD were equivalent to those acquired with the storage phosphor system.

CONCLUSION

The FPD provides radiographic images with excellent inherent physical image quality.

Performance characteristics of a Kodak computed radiography system

The performance characteristics of a photostimulable phosphor based computed radiographic (CR) system were studied. The modulation transfer function (MTF), noise power spectra (NPS), and detective quantum efficiency (DQE) of the Kodak Digital Science computed radiography (CR) system (Eastman Kodak Co.-model 400) were measured and compared to previously published results of a Fuji based CR system (Philips Medical Systems-PCR model 7000). To maximize comparability, the same measurement techniques and analysis methods were used. The DQE at four exposure levels (30, 3, 0.3, 0.03 mR) and two plate types (standard and high resolution) were calculated from the NPS and MTF measurements. The NPS was determined from two-dimensional Fourier analysis of uniformly exposed plates. The presampling MTF was determined from the Fourier transform (FT) of the system's finely sampled line spread function (LSF) as produced by a narrow slit. A comparison of the slit type ("beveled edge" versus "straight edge") and its effect on the resulting MTF measurements was also performed. The results show that both systems are comparable in resolution performance. The noise power studies indicated a higher level of noise for the Kodak images (approximately 20% at the low exposure levels and 40%-70% at higher exposure levels). Within the clinically relevant exposure range (0.3-3 mR), the resulting DQE for the Kodak plates ranged between 20%-50% lower than for the corresponding Fuji plates. Measurements of the presampling MTF with the two slit types have shown that a correction factor can be applied to compensate for transmission through the relief edges.

A rapid and automatic image registration algorithm with subpixel accuracy

A number of digital imaging techniques in medicine require the combination of multiple images. Using these techniques, it is essential that the images be adequately aligned and registered prior to addition, subtraction, or any other combination of the images. This paper describes an alignment routine developed to register an image of a fixed object containing a global offset error, rotation error, and magnification error relative to a second image. The described routine uses sparsely sampled regional correlation in a novel way to reduce computation time and avoid the use of markers and human interaction. The result is a fast, robust, and automatic alignment algorithm, with accuracy better than about 0.2 pixel in a test with clinical computed radiography images.

Plate scatter correction for improved performance in dual-energy imaging

We have developed a method for correcting for the scatter produced by the computed radiography plates in single-shot dual-energy imaging. In our technique, which uses a series of four computed radiography plates (the middle two plates act as a filter), the back plate can have scatter fractions as high as 30% due to the scatter from the first three plates. We use the deconvolution method to correct for this scatter, as the scatter medium is uniform and isotropic with a constant air gap. The ratio of the modulation transfer function (MTF) of a standard computed radiography plate to that of a dual-energy plate gives us the correction factor we need to scale the Fourier transforms of our images, thus compensating for the reduced response of the system due to scatter. We also investigated the effect of having an air gap between the back dual-energy plate and the first three plates. Our measurements show that this degrades the signal-to-noise response of the back plate slightly at low spatial frequencies. We also used independent scatter fraction measurements to verify the validity of our MTF measurements. Applying the correction to a back plate image improves the fine detail contrast by approximately 30%.

DQE(f) of four generations of computed radiography acquisition devices

Measurements were made of the MTF(f), NPS(f), and DQE(f) of four generations of computed radiography (CR) imaging plates and three generations of CR image readers. The MTF generally showed only a minor change between generations of plates and readers, but the DQE(f) has improved substantially from a very early plate/reader combination to a more recent one. The DQE in the more recent plate/reader combination is 1.3X greater at low frequencies and about 3X greater at high frequencies than the much earlier versions. Thus there has been substantial improvement in the imaging performance obtainable with CR since some of the early observer studies which indicated poorer performance with CR than with screen-film.

Effects of undersampling on the proper interpretation of modulation transfer function, noise power spectra, and noise equivalent quanta of digital imaging systems

The proper understanding of modulation transfer function (MTF), noise power spectra (NPS), and noise equivalent quanta (NEQ) in digital systems is significantly hampered when the systems are undersampled. Undersampling leads to three significant complications: (1) MTF and NPS do not behave as transfer amplitude and variance, respectively, of a single sinusoid, (2) the response of a digital system to a delta function is not spatially invariant and therefore does not fulfill certain technical requirements of classical analysis, and (3) NEQ loses its common meaning as maximum available SNR2 (signal-to-noise) at a particular frequency. These three complications cause the comparisons of MTF and NEQ between undersampled digital systems to depend on the frequency content of the images being evaluated. A tutorial of MTF, NPS, and NEQ concepts for digital systems is presented, along with a complete theoretical treatment of the above-mentioned complications from undersampling.

Variable compensation chest radiography performed with a computed radiography system: design considerations and initial clinical experience

The authors describe a variable compensation (VC) technique in which an x-ray equalizer and a computed radiography system are used. The VC technique allows retrospective alteration of equalized chest appearance with maintenance of improved signal-to-noise ratio in dense regions. Two imaging plates are used: one upstream of the patient to record the incident beam profile and one down-stream to record the equalized image. Subtraction of a weighted version of the upstream image from the down-stream image permits alteration of the appearance of the VC image, from the extremes of stimulated-unequalized to highly equalized. VC image appearance was optimized with a real-time workstation. The quality of VC images obtained in 33 patients was evaluated by three chest radiologists. Mediastinal appearance was better on VC equalized images than on conventional screen-film images. The stimulation of the appearance of a conventional radiograph with VC proved useful in interpretation of lung appearances on equalized radiographs.

Digital chest radiography with photostimulable storage phosphors: signal-to-noise ratio as a function of kilovoltage with matched exposure risk

A photostimulable storage phosphor (PSP) digital radiography system was evaluated regarding the signal-to-noise ratio (S/N) on chest images acquired at differing peak kilovoltage settings but with matched risk from radiation exposure. Images of two chest phantoms were acquired by using bedside (portable) imaging equipment at tube voltages ranging from 60 to 120 kV. Exposure factors were set so that the effective dose equivalent, a risk estimator weighted for various organs, was approximately equal in all exposures. The S/N in the lung-equivalent regions was found to be slightly higher (maximum, 15%) in the low-energy exposures, while the S/N values in the mediastinum- and subdiaphragm-equivalent regions were approximately equal at all kilovoltage settings. The absence of a high sensitivity of S/N to kilovoltage in risk-matched PSP images should enable institutions to select x-ray beam quality on the basis of other imaging criteria.

Direct digitization of optical images using a photostimulable phosphor system

The authors describe a method for directly digitizing optical images with a photostimulable phosphor (PSP) system. A PSP plate is initially charged with an exposure to a uniform x-ray field, and is then exposed to an optical image which discharges the plate in relation to the amount of incident light. Two applications were investigated: a contact-print technique for digitizing film radiographs, and a projection technique for digitizing transparent objects such as histology slides. Spatial uniformity was found to be adequate, and linearity of optical density response was excellent from 0.0-2.9 o.d. after look-up table correction. Spatial frequency response was degraded with the optical technique relative to the x-ray imaging properties of the plates, but was restorable by Fourier filtering. Image noise following spatial enhancement was satisfactory at intermediate to high optical densities using a high-resolution PSP plate, but was somewhat degraded at low densities.

Threshold perception performance with computed and screen-film radiography: implications for chest radiography

Images of a phantom obtained with computed radiography and standard screen-film imaging were compared to evaluate observer threshold perception performance with a modified contrast-detail technique. Optimum exposure necessary for performance with the imaging plate technique to match that with screen-film techniques was determined, as was comparative performance with variation in kilovoltages, plate type, spatial enhancement, and hard-copy interpolation method. It was found that computed radiography necessitates about 75%-100% more exposure than screen-film radiography to optimally match performance with Ortho-C film with Lanex regular or medium screens (Eastman Kodak, Rochester, NY) for detection of objects 0.05-2.0 cm in diameter. However, only minimal loss of detection performance (approximately 10% overall) was experienced if standard screen-film exposures were used with computed radiography. Little change in observer performance was found with variation in plate type, spatial enhancement, or method of hard-copy interpolation. However, perception performance with computed radiographic images was better at lower kilovoltages.

Single-exposure conventional and computed radiography image acquisition

A technique for simultaneously acquiring a conventional film-screen radiographic image and a digital computed radiography (CR) image with a single x-ray exposure is described. Measurements of image contrast, spatial resolution, and signal-to-noise ratios demonstrate that a modified film cassette in which the first intensifier screen has been replaced with a CR imaging plate permits dual-image, single-exposure imaging with only nominal degradation in film and CR image quality relative to the two standard image counterparts. This technique may be used to acquire matched image pairs for research or as a way to provide full-size conventional film images in the clinical environment, while retaining the advantages offered by computed radiography systems.

Scatter fractions in AMBER imaging

Images of two phantoms were obtained with use of an advanced multiple-beam equalization radiography system, and scatter fractions were estimated with use of a photostimulable phosphor imaging system. Scatter fractions in the equalized images were lower in the mediastinum-equivalent areas and higher in the lung-equivalent areas, relative to images that were conventionally acquired with use of an antiscatter grid. The differences are attributed to a reduction in incident exposure in the lungs and the presence of cross-scatter between lung and mediastinal regions.