Abstract

The aim of the present study was to examine the image quality performance of a CMOS digital imaging optical sensor coupled to custom made gadolinium oxysulfide powder scintillators, doped with praseodymium, cerium, and fluorine (Gd2O2S:Pr,Ce,F). The screens, with coating thicknesses 35.7 and 71.2 mg/cm2, were prepared in our laboratory from Gd2O2S:Pr,Ce,F powder (Phosphor Technology, Ltd.) by sedimentation on silica substrates and were placed in direct contact with the optical sensor. Image quality was determined through single index (information capacity, IC) and spatial frequency dependent parameters, by assessing the Modulation Transfer Function (MTF) and the Normalized Noise Power Spectrum (NNPS). The MTF was measured using the slanted-edge method. The CMOS sensor/Gd2O2S:Pr,Ce,F screens combinations were irradiated under the RQA-5 (IEC 62220-1) beam quality. The detector response function was linear for the exposure range under investigation. Under the general radiography conditions, both Gd2O2S:Pr,Ce,F screen/CMOS combinations exhibited moderate imaging properties, in terms of IC, with previously published scintillators, such as CsI:Tl, Gd2O2S:Tb, and Gd2O2S:Eu.

1. Introduction

Indirect medical imaging detectors incorporate a scintillator detector coupled to an optical sensor (CCD, CMOS, etc.). The scintillator is used in the form of either screens or needle phosphor layers [1]. Optical ceramic scintillators have been also developed in order to replace single crystals in some applications [14]. Ceramic gadolinium oxysulfide doped with praseodymium, cerium, and fluorine (Gd2O2S:Pr,Ce,F) (ultrafast ceramics, UFC) has been proposed for Computed Tomography (CT) applications [4]. Gd2O2S:Pr,Ce,F is an efficient and fast scintillator (decay time of the order of 3-4 μs) based on the well-known traditional Gd2O2S host material [5].

A semiconductor technology that has been used widely recently is complementary metal oxide (CMOS) semiconductors [69]. The initial CMOS pixel structure was in the form of passive pixel sensors (PPSs) and was followed by the active pixel sensors (APS), which have been developed with the goal of improving the image quality [10]. APS CMOS sensors provide high resolution even at high framing rates and, in association with scintillating screens, have been increasingly investigated for medical imaging applications [11]. In the present study, image quality of a CMOS digital imaging sensor, coupled to Gd2O2S:Pr,Ce,F powder phosphor screens, was investigated by single index and spatial frequency dependent parameters.

The single index information capacity (IC) was considered as an overall image quality metric, including the effects of sharpness, contrast, and noise of the CMOS sensor coupled to Gd2O2S:Pr,Ce,F powder phosphor screens. IC was estimated by the experimental measurement of the Signal Transfer Property (STP), the Modulation Transfer Function (MTF), and the Normalized Noise Power Spectrum (NNPS).

2. Materials and Methods

2.1. Phosphor Screens

Gd2O2S:Pr,Ce,F was purchased in powder form (Phosphor Technology Ltd., England, code: UKL59CF/N-R1) with a mean grain size (estimated by ultrasonic dispersion with a coulter counter having 100 μm aperture) of approximately 9.7 μm (Phosphor Technology Ltd., datasheet) [5]. Gd2O2S:Pr,Ce,F has effective atomic number , density of 7.34 g/cm3, and a decay time of the order of a 3 × 10−6 s [3].

The phosphor was used in the form of thin layers to simulate the intensifying screens employed in X-ray imaging. Two screens with coating thicknesses 35.7 and 71.2 mg/cm2 were prepared by sedimentation of Gd2O2S:Pr,Ce,F powder on fused silica substrates (spectrosil B). Sodium orthosilicate (Na2SiO3) was used as binding material between the powder grains.

Gd2O2S:Pr,Ce,F powder phosphor particles morphology was verified via scanning electron microscope (SEM) micrographs using the Jeol JSM 5310 SEM combined with the INCA software (Figure 1).

2.2. CMOS Sensor

The Gd2O2S:Pr,Ce,F scintillating screens were manually coupled to an optical readout device including a CMOS Remote RadEye HR photodiode pixel array [9, 12]. The CMOS photodiode array consists of 1200 × 1600 pixels with 22.5 μm pixel spacing. The Gd2O2S:Pr,Ce,F screens were directly overlaid onto the active area of the CMOS photodiode array, consisting of an N-well diffusion on p-type epitaxial silicon, and held by using a thin polyurethane foam layer for compression between the screen and a 1 mm thick graphite cover. A 70 kV (RQA-5) X-ray beam was used, following the IEC standards [13]. IEC standard X-ray spectrum was achieved by adding 21 mm Al filtration in the beam to simulate beam quality alternation by a human body. A BMΙ General Medical Merate tube with rotating tungsten anode and inherent filtration equivalent to 2 mm Al was used for the RQA-5 beam quality. The source-to-detector distance (SDD) between the X-ray focal spot and the surface of the detector was set to 176 cm. The added filtration was placed as close as possible to the source.

2.3. Image Quality
2.3.1. Signal Transfer Property (STP)

The Signal Transfer Property (STP) provides the relationship between mean pixel value (MPV) and Entrance Surface Air Kerma (ESAK) at the detector surface. This relationship was obtained by plotting pixel values (PV) versus ESAK at the detector, as described in the IEC method [13]. A sequence of uniform images was acquired at different exposure levels. MPV was evaluated in a 1 × 1 cm region of interest (ROI). The system’s response curve was fitted using a linear equation of the following form:where and are fit parameters and is the incident air kerma (ESAK). From the slope of the system’s response curve, the value of the gain factor () was obtained. This value is relating MPV to the incident exposure at the detector (in digital units per ) [14]. The magnitude of the pixel offset at zero air kerma was also estimated [15].

2.3.2. Modulation Transfer Function (MTF)

The spatial frequency transfer characteristics of linear and spatially invariant imaging systems can be characterized by the system’s Modulation Transfer Function (MTF) [16, 17]. MTF describes the variation of contrast with spatial frequency and hence expresses spatial resolution. MTF of an imaging system is defined as the absolute value of its optical transfer function, normalized to unity at spatial frequency zero [14]:

In this study, the MTF was measured using the slanted-edge technique, following the procedures described in IEC standard [13, 15]. A PTW Freiburg tungsten edge test device was used to obtain the slanted-edge images. The edge test device consists of a 1 mm thick W edge plate (100 × 75 mm2) fixed on a 3 mm thick lead plate. Images of the edge, placed at a slight angle in order to avoid aliasing effects, were obtained. The edge spread function (ESF) was calculated by the extraction of a 1 × 1 cm2 ROI with the edge roughly at the center. The restricted ROI area could not accurately measure the low frequency drop, in powder screens; however, the size is limited from the size of the screen which is 2 × 2 cm2 and the active area of the CMOS sensor (2.7 × 3.6 cm). Neglecting the low frequency drop could result in overestimation of the MTF and thus also of the DQE values.

2.3.3. Normalized Noise Power Spectrum (NNPS)

The NPS was calculated according to IEC 62220-1-1:2015 [13]. For each ROI, the pixel values were converted into air-kerma units. The slowly varying spatial background effects including the heel effect were corrected by fitting and subtracting a two-dimensional second-order polynomial to the original acquired image data [18]. All images were acquired with offset and gain correction. The area of analysis was 1024 × 1024 pixels. Half-overlapping ROIs with a size of 128 × 128 pixels were then taken from the subimages [13]. A total of 128 ROIs were taken from each flood image. For all the ROIs taken from each image 2D fast Fourier transform (FFT) of each ROI was calculated and added to the NPS ensemble. NNPS was obtained by dividing NPS by the square of the corresponding and afterwards the ensemble average was obtained.

2.3.4. Detective Quantum Efficiency (DQE) and Image Quality Figure of Merit (IQFOM)

The Detective Quantum Efficiency has been defined as the efficiency of a system to transfer the Signal-to-Noise Ratio (SNR) through the imaging chain [19]. The DQE of the system can be calculated by the following equation:where is the number of photons per unit kerma (μGy) per mm2 [9]. A value of 29563 photons × μGy−1 × mm−2 is required, for the RQA-5 beam quality, according to [13]. However, in the present study, a value of 21627 photons × μGy−1 × mm−2 was calculated from direct X-ray spectra measurements with a portable Amptek XR-100T spectrometer and the extracted ESAK values were validated against direct measurements with a Victoreen 4000M+ dosimeter [12]. Accordingly, (3) was used for calculations; however, an alternative term was employed, namely, Image Quality Figure of Merit (IQFOM), due to the difference in the values proposed by the IEC protocol [13].

2.3.5. Information Capacity (IC)

The concept of image information capacity (IC) has been introduced within the context of Shannon’s information theory, in order to assess image information content [20]. In digital imaging, the continuous spatial distribution of an optically generated image is sampled by the discrete sensitive pixels on a photodiode array, whose outputs are converted into digitized signals and stored in an image processing system for numerical evaluation. According to Shannon’s information theory, the image information capacity, per unit of image area, may be defined as follows [20]:where are the different signal intensity levels, is the duration of the signal, and is the number of pixels of the imaging system. In this study, information capacity was calculated according to [21]where is the signal power spectrum which can be defined as , where is the gain factor [22].

3. Results and Discussion

Figure 2 shows the detector response curves (STP) of the CMOS sensor combined with the 35.7 and 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screens, respectively, under the RQA-5 (70 kVp) beam quality. The detector was found to have a linear response, covering the whole exposure range, with a pixel value offset of −8.193 and 15.19 for the 35.7 and 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screens, respectively.

The linear no-threshold fits gave correlation coefficients () greater than 0.9992 and 0.9994 for the 35.7 and 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screens. The gain factor was determined as the slope of the characteristic curve, relating the mean pixel value to the incident exposure. Using flat-field images, the gain factors were determined by linear regression to be and 7.902 digital units per for the 35.7 and 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screens under the RQA-5 beam quality. However, the values of are restricted by the bit depth of the CMOS sensor under investigation, compared to large flat panel detectors, dedicated for medical imaging [23].

Figure 3 shows MTF curves, for the 35.7 and 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screens, coupled to the CMOS sensor. Under the RQA-5 conditions, the 35.7 mg/cm2 screen/CMOS combinations show higher MTF values, compared to the thicker 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screen, across the spatial frequency range, except for some frequencies, where the thick screen shows systematic discrepancies in the MTF values, which may be due to screen nonuniformities. For comparison purposes, data for Gd2O2S:Tb [12] and CsI:Tl [24] scintillators detectors are shown. Both Gd2O2S:Tb and CsI:Tl screens were produced industrially (not custom made). Furthermore CsI:Tl is in columnar form; thus, light can be easily transmitted to the output in order to produce narrow light pulses, which in turn lead to improved resolution. The Gd2O2S:Pr,Ce,F screen combinations showed lower MTF values compared to previously published scintillators. However, both Gd2O2S:Pr,Ce,F screens show comparable MTF results with the 65.1 mg/cm2 Europium doped Gd2O2S (Gd2O2S:Eu) screen [21] in all the spatial frequency range, showing that the resolution properties of this material could be considered for general radiography applications.

Figure 4 shows the extracted 1D Normalized Noise Power Spectrum (NNPS), for the direction, obtained from uniformly exposed images, at an exposure level of 11.1 , which is close to the usual clinical dose ranges for medical imaging systems [23]. The NNPS curves are falling with increasing spatial frequency. A slight decrease from 4.56 × 10−6 mm2 (at 2.08 cycles/mm2) to 1.1 × 10−6 (at 9 cycles/mm2) was observed in the NNPS values of the 35.7 mg/cm2 Gd2O2S:Pr,Ce,F screen coupled to the CMOS sensor. This is attributed to the combined noise effects of the screen and the CMOS sensor, in which a 2D correction is applied reducing the structured noise in each image. Noise levels in the 35.7 mg/cm2 Gd2O2S:Pr,Ce,F screen are in acceptable levels and can be compared with the CsI:Tl screen, for spatial frequencies up to 4 cycles/mm. In the thicker 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screen coupled to the CMOS sensor, NNPS levels are rather higher due to the presence of low frequency components (spatial frequencies up to 4.5 cycles/mm) that can be attributed to screen inhomogeneities, due to the sedimentation production process. A decrease from 4.92 × 10−5 mm2 (at 2.08 cycles/mm2) to 8.26 × 10−7 (at 9 cycles/mm2) was observed in the NNPS values of the 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screen coupled to the CMOS sensor. Noise levels of the 65.1 mg/cm2 Gd2O2S:Eu screen are in between both Gd2O2S:Pr,Ce,F screens, being in accordance with the noise properties following screen thickness. As was expected, Gd2O2S:Tb screen shows the lowest noise levels, combining screen thickness, internal properties of the material, and the industrial production process.

Figure 5 shows IQFOM curves obtained for the 35.7 and 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screens, coupled to the CMOS sensor. DQE and IQFOM values are principally affected by the variation of the ratio of the squared MTF over NNPS. IQFOM of the 35.7 mg/cm2 screen/CMOS combination was found to maximize in the medium spatial frequency range and drop thereafter, whereas the corresponding IQFOM values of 71.2 mg/cm2 screen/CMOS combination maximize in the low spatial frequency range.

Table 1 shows information capacity values for the combination of the CMOS sensor with the 35.7 and 71.2 mg/cm2 Gd2O2S:Pr,Ce,F screens under investigation and previously published IC values for a CMOS sensor coupled to Gd2O2S powder scintillators, activated either with terbium (Tb) or Europium (Eu). The comparison, under the RQA-5 beam quality, was obtained at the same exposure level for all screen/sensor combinations. The thinner Gd2O2S:Tb screen (33.91 mg/cm2) showed the highest IC values due to the screen thickness and the higher MTF values of this screen, compared to the 35.7 mg/cm2 Gd2O2S:Pr,Ce,F. The IC values of the thicker 71.2 mg/cm2 Gd2O2S:Pr,Ce,F were lower than all the other screens, even from the IC value of the 65.1 mg/cm2 Gd2O2S:Eu screen due to the fact of the higher thickness and screen uniformity, leading to lower MTF values and higher noise. These data show that, for a given level of incident X-ray fluence, information capacity is mainly determined by the intrinsic phosphor material properties and by the screen thickness of the imaging system. In thick screens, the lateral light trajectories are very long causing a large fraction of the laterally directed photons to be absorbed before reaching the screen output.

4. Conclusion

In the present study, image quality of two custom made Gd2O2S:Pr,Ce,F powder scintillator screens, coupled to a CMOS digital imaging sensor, was investigated under X-ray radiography imaging conditions. Image quality was investigated in terms of single index and spatial frequency dependent parameters. The detector response function was linear for the exposure range under investigation. The overall imaging properties, in terms of the single index IC, of both Gd2O2S:Pr,Ce,F screen/CMOS combinations were found moderate compared to previously published scintillators, such as Gd2O2S:Tb, CsI:Tl, and Gd2O2S:Eu.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research has been cofunded by the European Union (European Social Fund) and Greek National Resources under the framework of the “ARISTEIA” project MISCIRLU Code 1476 of the “Education and Lifelong Learning” Operational Programme.