Complete MTF evaluation of two cone beam CT systems

Objective To evaluate the parallactic unsharpness caused by flat panel detector in digital radiography using different incident angle of central ray. Methods R-1 square-wave phantom was exposed by Kodak DR3000 system with X-ray tube angled 0°, 10°, 20°, 30°, and 40° respectively. Then,presampled modulation transfer function (MTF) was calculated for each case above. Subsequently,experimental data were processed and Wilcoxon signed ranks test were performed by statistic software SPSS 10. 0, in which P ＜ 0. 05 was considered as statistically significant difference. Results The presampled MTF curves of incident angle of 0°-40° degree, were presented orderly from top to bottom, especially the incident angle of 40° was obviously the lowest. The incident angle of 0° was considered as a control group and other groups were compared against it. There was no statistically significant difference for MTF of incident angle of 10°( Z = - 1. 893 ,P = 0. 058 ), while there were significant difference for MTF of incident angle of 20°, 30°, and 40° ( Z = - 2. 547, - 2. 666, - 2. 666, P ＜ 0. 05 ). Conclusions For flat panel detector in digital radiography, the bigger the incident angle of central ray, the larger the parallactic unsharpness. In addition, this effect has less influence on structures of low spatial frequency than those of high spatial frequency. Key words: Radiography; Evaluation studies; Quality control

In indirect conversion flat panel detectors (FPDs) for digital X-ray imaging, scintillating materials such as Terbiumdoped Gadolinium Oxysulfide (Gadox) convert X-ray into visible light, and an amorphous silicon (a-Si) photodiode array converts the light into electrons. It is, however, desired that the detector spatial resolution is improved because the light spreading inside scintillator causes crosstalk to next a-Si photodiode pixels and the resolution is degraded compared with direct conversion FPDs which directly convert X-ray into electrons by scintillating material such as amorphous selenium. In this study, the scintillator was pixelated with same pixel pitch as a-Si photodiode array by barrier rib structure to limit the light spreading, and the detector spatial resolution was improved. The FPD with pixelated scintillator was manufactured as follows. The barrier rib structure with 127μm pitch was fabricated on a substrate by a photosensitive organic-inorganic paste method, and a reflective layer was coated on the surface of the barrier rib, then the structure was filled up with Gadox particles. The pixelated scintillator was aligned with 127μm pixel pitch of a-Si photodiode array and set as a FPD. The FPD with pixelated scintillator showed high modulation transfer function (MTF) and 0.94 at 1cycle/mm and 0.88 at 2cycles/mm were achieved. The MTF values were almost equal to the maximum value that can be theoretically achieved in the FPD with 127μm pixel pitch of a-Si photodiode array. Thus the FPD with pixelated scintillators has great potential to apply for high spatial resolution applications such as mammography and nondestructive testing.

We investigated the effect on image data resampling in an evaluation of the basic imaging properties for a digital radiographic system based on a flat panel detector (FPD). One of the latest digital radiographic systems was used in this study. This system was based on a direct-conversion FPD of amorphous selenium. The basic imaging properties of the system were evaluated by measuring characteristic curve, presampled modulation transfer function (MTF), and Wiener spectrum (WS) using DICOM image with a matrix size of 2048 x 2048. The evaluations were performed under two conditions because matrix size automatically changes according to the selection of imaging size. One of the conditions was a different matrix size between image data acquired on the FPD and the output image (DICOM image for which resampling was performed). The other condition was that these matrices be the same size (DICOM image with no resampling performed). Resampling did not affect the characteristic curves. However, MTF and the WS obtained from the resampled data were different from those of the one not resampled, which is considered to be the "inherent" basic imaging properties, and this phenomenon was remarkable, especially in terms of the MTFs. Our study indicates that the effect on resampling should not be disregarded in evaluating the basic imaging properties of digital radiographic systems. Therefore, it is mandatory to use DICOM images for which no resampling was performed in order to evaluate the inherent basic imaging properties for digital radiographic 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, Philips Digital Diagnost, GE Revolution XQ/i, and Hologic Direct-Ray DR1000 using the IEC-defined RQA5 (approximately 75 kVp, 21 mm Al) and RQA9 (approximately 120 kVp, 40 mm Al) radiographic techniques. The presampled MTF of the systems was measured using an edge method (Samei et al., Med Phys 25:102, 1998). The NPS of the systems was determined for a range of exposure levels by 2D Fourier analysis of uniformly exposed radiographs (Flynn and Samei, Med Phys 26:1612, 1999). The ideal signal-to-noise ratio per exposure for each system was estimated using a semi-empirical x-ray model. The DQE, reported only at the RQA5 technique, was assessed from the measured MTF, NPS, exposure, and the ideal signal-to-noise ratio. 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. The DQE results reflected expected differences based on the absorption efficiency of the different detector materials. Using RQA5 and 0.3 mR exposure, the measured DQE values at spatial frequencies of 0.15 mm-1 and 2.5 mm-1 were 64% and 14% for the XQ/i system and 35% and 19% for DR-1000. Using RQA5 and the averages at all exposures, the corresponding values were 58% and 13% for the XQ/i system and 36% and 19% for DR-1000.

Line-pair resolution phantoms are used to determine the spatial resolution of medical imaging systems. In some cases, these phantoms are used to determine the maximum number of line-pairs per mm that the system can resolve. In other cases, a numerical determination of the modulation transfer function (MTF) is carried out by means of the analysis of the variance of ROIs on the image. In this note, a new procedure is implemented to calculate the presampled MTF of an imaging system. Images of a commercial line-pair phantom are acquired in a flat panel detector. After applying an edge detector and a Radon transform to the image, the direction of the bars in the phantom is calculated. Then, an area of the image that excludes the ends of the bars is determined. Every pixel in this area is used to obtain an oversampled profile on which the MTF is calculated. Every group of line-pairs conforms a periodic wave on this profile. Each of these waves is the output of the system for a square wave input with the same frequency. After extracting a wave, the MTF value for its frequency is calculated as the ratio of its first odd harmonic to the first odd harmonic of the square wave input. The amplitude of this square wave is obtained from two uniform areas on the phantom image. The results obtained are compared to those obtained following the standard edge method, as recommended by the IEC, and show a very good agreement between both methods, in both main directions of the detector and all the dose ranges analyzed. The presented method is shown accurate and can be used to extend the conventional use of line-pair phantoms in conventional radiology.

Planar X-ray mammography is the standard medical imaging modality for the early detection of breast cancer. Based on advancements in digital flat-panel detector technology, dedicated x-ray computed tomography (CT) mammography is a modality under investigation that offers the potential for improved breast tumor imaging. We have implemented a prototype half cone-beam CT breast imaging system that utilizes an indirect flat-panel detector. This prototype can be used to explore and evaluate the effect of varying acquisition and reconstruction parameters on image quality. This report describes our system and characterizes the performance of the system through the analysis of Modulation Transfer Function (MTF) and Noise Power Spectrum (NPS). All CT reconstructions were made using Feldkamp's filtered backprojection algorithm. The 3D MTF was determined by the analysis of the plane spread function (PlSF) derived from the surface spread function (SSF) of reconstructed 6.3mm spheres. 3D NPS characterization was performed through the analysis of a 3D volume extracted from zero-mean CT noise of air reconstructions. The effect of varying locations on MTF and the effect of different Butterworth filter cutoff frequencies on NPS are reported. Finally, we present CT images of mastectomy excised breast tissue. Breast specimen images were acquired on our CTMS using an x-ray technique similar to the one used during performance characterization. Specimen images demonstrate the inherent CT capability to reduce the masking effect of anatomical noise. Both the quantitative system characterization and the breast specimen images continue to reinforce the hope that dedicated flat-panel detector, x-ray cone-beam CT will eventually provide enhanced breast cancer detection capability.

Two main spatial frequency and quantitative calculations and measurements am, used to study the imaging system of three industrial digital radiography (IDR) modules for non-destructive testing (NDT), namely complementary metal oxide semiconductor (CMOS) flat-panel digital detector with 50 mu m pixel pitch, computed radiography (CR) with phosphor imaging plate with 25 mu m pixel pitch and laser-type film digitiser with 50 mu m pixel pitch. The modulation transfer function (MTF) and noise power spectrum (NPS) measurement and calculations were adapted on order to evaluate the image quality of IDR images. The MTF has been used to characterise the resolution properties of analogue and digital X-ray imaging systems and measure image deterioration due to optical factors, such as diffusion of image forming radiation. The NPS is used in determining noise in images produced and the choices of normalisation of the NPS vary slightly from researcher to researcher, but the underlying methodology is consistent. In order to easily compare these two measurements and calculations, the integral of up to the system Nyquist frequency was used as the final image quality quantitative evaluation. For the MTF measurement and calculation, four regions of interest (ROI) were, defined to portray the randomness in performing real testing. The four ROI are then compared with each other. For the NPS measurement and calculation, ten ROI were defined, scattered and averaged to have overall noise estimation. From the measurements and calculations which were done, the averaged MTF at 20% modulation for the laser-type film digitiser system is 6.55 cycles/mm, CR modulated 4.48 and CMOS modulated 2.83 cycles/mm. For NPS measurement and calculation, the result obtained in decreasing order is laser-type film digitiser, CMOS and CR. The laser-type film digitiser system modulates the best transfer function at 20% but does not have the lowest NPS. The lowest and stable NPS is the CR system but it has the lowest modulating capability at 20%. The study shows that in order to perform NDT by using the evaluated modules, the user must know the true capability, of the system and how, it is designed for specific application and discontinuity detection.

To investigate the image quality of a digital radiography system with an amorphous-silicon, large-area, digital flat-panel detector. A flat-panel detector based on a matrix of amorphous silicon was integrated into a projection radiography system. The scintillator consisted of a layer of structured cesium iodide. The active matrix size of 30002 pixels together with a pixel size of 143 microm provided a large image area of 43 x 43 cm2. Basic image quality parameters such as detective quantum efficiency (DQE) and modulation transfer function (MTF) were measured and compared with those obtained with conventional systems. The measurement of DQE yielded a high value of 70% at zero spatial frequency. At a system dose equivalent to 400 speed, the DQE of the digital system was a factor of two larger than the DQE of a storage phosphor or screen-film system within the entire spatial frequency range between zero and the Nyquist limit of 3.5 line pairs per millimeter. The flat-panel detector furthermore has an MTF that is superior to that in regular screen-film systems and also provides a substantially larger dynamic range. This new technology demonstrates its potential to provide equal or superior image quality to conventional screen-film systems and to reduce patient exposure to radiation dose. The advantages of digital radiography systems, based on a flat-panel detector as an instant image display, facilitation of work flow in the radiology department, and digital networking and archiving, are well in sight.

The authors describe a new technique to determine the system presampled modulation transfer function (MTF) in digital radiography using only the detector noise response. A cascaded-linear systems analysis was used to develop an exact relationship between the two-dimensional noise power spectrum (NPS) and the presampled MTF for a generalized detector system. This relationship was then utilized to determine the two-dimensional presampled MTF. For simplicity, aliasing of the correlated noise component of the NPS was assumed to be negligible. Accuracy of this method was investigated using simulated images from a simple detector model in which the "true" MTF was known exactly. Measurements were also performed on three detector technologies (an x-ray image intensifier, an indirect flat panel detector, and a solid state x-ray image intensifier), and the results were compared using the standard edge-response method. Flat-field and edge images were acquired and analyzed according to guidelines set forth by the International Electrotechnical Commission, using the RQA 5 spectrum. The presampled MTF determined using the noise-response method for the simulated detector system was in close agreement with the true MTF with an averaged percent difference of 0.3% and a maximum difference of 1.1% observed at the Nyquist frequency (fN). The edge-response method of the simulated detector system also showed very good agreement at lower spatial frequencies (less than 0.5 fN) with an averaged percent difference of 1.6% but showed significant discrepancies at higher spatial frequencies (greater than 0.5 fN) with an averaged percent difference of 17%. Discrepancies were in part a result of noise in the edge image and phasing errors. For all three detector systems, the MTFs obtained using the two methods were found to be in good agreement at spatial frequencies less than 0.5 fN with an averaged percent difference of 3.4%. Above 0.5 fN, differences increased to an average of 20%. Deviations of the experimental results largely followed the trend seen in the simulation results, suggesting that differences between the two methods could be explained as resulting from the inherent inaccuracies of the edge-response measurement technique used in this study. Aliasing of the correlated noise component was shown to have a minimal effect on the measured MTF for the three detectors studied. Systems with significant aliasing of the correlated noise component (e.g., a-Se based detectors) would likely require a more sophisticated fitting scheme to provide accurate results. Results indicate that the noise-response method, a simple technique, can be used to accurately measure the MTF of digital x-ray detectors, while alleviating the problems and inaccuracies associated with use of precision test objects, such as a slit or an edge.

Quality assurance in computed tomography (CT) is commonly performed with the Fourier-based modulation transfer function (MTF) and the noise variance, while more recently the noise power spectrum (NPS) has increased in popularity. The Fourier-based methods make assumptions such as shift-invariance and cyclostationarity. These assumptions are violated in real clinical systems and consequently are expected to result in systematic errors. A spatial approach, based on the object transfer matrix (T) and the covariance matrix (K) theory, does not require these assumptions and can provide a more general description of the imaging system. In this paper, the authors present an experimental methodology and data treatment for quality assessment of a lab cone-beam CT system by comparing the spatial with the Fourier approach in 2D reconstructed slices. In order to have control over all experimental parameters and image reconstruction, a bench-top flat-panel-based cone-beam CT scanner and a cylindrical water-filled poly(methyl methacrylate) (PMMA) phantom were used for the noise measurements. An aluminum foil inserted in the water phantom enabled the estimation of the line response function (LRF) with a limited number of assumptions. The authors evaluated the spatial blur, the noise and the signal-to-noise ratio (SNR) using the spatial approach as well as the Fourier-based approach. In order to evaluate the degree of noise nonstationarity of their cone-beam CT system, the authors evaluated both the local and global CT noise properties and compared them using both approaches. For the laboratory cone-beam CT, the location-dependent noise evaluation showed that in addition to the noise variance, the NPS and covariance eigenvector symmetry depend on the location in the image. The estimated signal transfer was similar for both approaches. Unlike the Fourier approach which uses the same exponential wave function basis for both MTF and NPS, the eigenvectors of T and K were significantly different. By using the eigenvectors of the noise and object transfer to characterize the system, the spatial approach provides additional information to the Fourier approach and is therefore an important tool for a thorough understanding of a CT system.

Dual Energy (DE) imaging has been widely used in digital radiography and fluoroscopy, as has dual energy CT for various medical applications. In this study, the imaging performance of a dynamic dual-layer a-Si flat panel detector (FPD) prototype was characterized for dual energy imaging tasks. Dual energy cone beam CT (DE CBCT) scans were acquired and used to perform material decomposition in the projection domain, followed by reconstruction to generate material specific and virtual monoenergetic (VM) images. The dual-layer FPD prototype was built on a Varex XRD 4343RF detector by adding a 200 μm thick CsI scintillator and a-Si panel of 150 μm pixel size on top as a low energy detector. A 1 mm copper filter was added as a middle layer to increase energy separation with the bottom layer as a high energy detector. The imaging performance, such as Modulation Transfer Function (MTF), Conversion Factor (CF), and Detector Quantum Efficiency (DQE) of both the top and bottom detector layers were characterized and compared with those of the standard single layer XRD4343 RF detector. Several tissue equivalent cylinders (solid water, liquid water, bone, acrylic, polyethylene, etc.) were placed on a rotating stand, and two separate 450-projection CBCT scans were performed under continuous 120 kV and 80 kV X-ray beams. After an empirical material decomposition calibration, water and bone images were generated for each projection, and their respective volumes were reconstructed using Varex’s CBCT Software Tools (CST 2.0). A VM image, which maximized the contrast-to-noise ratio of water to polyethylene, was generated based on the water and bone images. The MTF at 1.0 lp/mm from the low energy detector was 32% and 22% higher than the high energy detector and the standard detector, respectively; the DQE of both high and low energy detectors is much lower than that of the standard XRD 4343RF detector. The CNR of water to polyethylene from the VM image improved by 50% over that from the low energy image alone at 120 kV, and by 80% at 80 kV. This study demonstrates the feasibility of using a dual-layer FPD in applications such as DE CBCT for contrast enhancement and material decomposition. Further evaluations are underway.

Characterization and application of a GE amorphous silicon flat panel detector in a synchrotron light source

Our work is to investigate and understand the factors affecting the imaging performance of amorphous selenium (a-Se) flat-panel detectors for digital mammography. Both theoretical and experimental methods were developed to investigate the spatial frequency dependent detective quantum efficiency [DQE(f)] of a-Se flat-panel detectors for digital mammography. Since the k-edge of a-Se is 12.66 keV and within the energy range of a mammographic spectrum, a cascaded linear system model was developed which takes into account the effect of k-fluorescence on the modulation transfer function (MTF), noise power spectrum (NPS) and DQE(f) of the detector. This model was used to understand the performance of a prototype detector with 85 mm pixel size. The presampling MTF, NPS and DQE(f) of the prototype were measured, and compared to the theoretical calculation by the model. The calculation showed that k-fluorescence reduces the MTF by 15% at the Nyquist frequency (f<SUB>NY</SUB>) of the prototype detector, and the NPS at f<SUB>NY</SUB> was reduced to 82% of that at zero spatial frequency. Because of the decrease in both MTF and NPS at high spatial frequencies, k-fluorescence only has a small degradation effect on DQE(f) for mammography. The measurement of presampling MTF of the prototype detector revealed an additional source of blurring, which was attributed to the blocking layer at the interface between a-Se and the active matrix. This introduced high frequency drop in both presampling MTF and NPS, and reduced aliasing in the NPS. As a result, the DQE(f) of the prototype detector at f<SUB>NY</SUB> approaches 50% of that at zero spatial frequency.

To compare the results of a new quantitative image quality evaluation method that requires no observers with the results of receiver operating characteristic (ROC) analysis in detecting the mandibular canal (MC) in cone beam CT (CBCT) images. A Teflon (polytetrafluoroethylene) plate phantom with holes of different depths was scanned with two CBCT systems. One CBCT system was equipped with an image intensifier (Experiment 1), and the other was equipped with a flat panel detector (Experiment 2). Holes that were above the threshold gray value (ΔG), calculated using just-noticeable difference (JND), were extracted. The number of extracted holes was used as the index of the image quality, and was compared with the Az values calculated by ROC analysis to detect the MC. The number of extracted holes reflected the influence of different scanning conditions, and showed a strong correlation with the Az values calculated by ROC analysis. Indices of the number of extracted holes corresponding to high Az values for detecting the MC were obtained in both experiments. Our image quality evaluation method applying JND to images of a standardized phantom is a quantitative method that could be useful for evaluating the detectability of the MC in CBCT images.

The purpose of this study was to investigate image quality identification methods among direct- and indirect-type flat panel detectors (FPDs) and a computed radiography (CR) system using two radiation qualities RQA3 and RQA5 defined in the IEC 61267 standard. For each system, the digital characteristic curve, the presampled modulation transfer function (MTF), and the normalized noise power spectrum (NNPS) were measured. Images for a burger phantom and a foot-bone phantom were processed by resolution identification utilizing two-dimensional Fourier transform, and then contrast-to-noise ratio (CNR) for each image was measured. For the RQA3, the direct FPD system indicated the highest DQE value, and for the RQA5 DQE value of indirect FPD system, it was a little higher than that of direct FPD system. The CNR results with the resolution identification displayed good accordance with the DQE results in both phantoms. From the DQE results, dose ratios for image quality identification were determined, and the CNRs of the dose-adjusted images were measured. The results for, the CNRs of all systems showed good coincidence. From these findings, they indicated that the DQE measurement is effective to determine the exposure parameters for equalizing the image quality of different types of radiographic systems.