Introduction
Computed tomography (CT) is one of the most commonly utilized methods for acquiring three-dimensional images of anatomy in vivo. CT generates reconstructed images as a distribution of the linear attenuation coefficient, which is averaged by the X-ray energy spectrum.
Dual-energy CT (DECT) was first developed in the 1970s [1], and the clinical application of DECT has recently been realized as a result of robust improvements in performance. The principle of DECT is the acquisition of two images of the same anatomy with different kilo-voltage (kV) settings (usually 80/140 kV). Clinically, different DECT systems have been developed by the commercial industry, which include a rapid kV switching system with a single-layer detector [2], dual-source system with dual detectors [3], and single-source system with a dual-layer detector [4]. The DECT technique is utilized for artifact reduction [5], contrast enhancement [6], and noise reduction by mixing or synthesizing between two image sets with different spectral information [7].
An approach to generate a single set of images is to synthesize virtual monochromatic images using projection data [5, 6]. This approach was described by previous work [8, 9], which involves basis-material decomposition in the projection domain and a linear combination of the density map the basis materials. In principle, such a process can provide quantitative information and remove beam-hardening artifacts on the imaged anatomy. Beam-hardening artifact reduction is considered one of the main benefits of virtual monochromatic DECT images.
Another approach is to generate processed images from low- and high-energy images based on the image domain. The image-based virtual monochromatic images contain beam-hardening artifacts propagated from the low- and high-energy images. However, the objective of image-based monochromatic images is primarily to generate a single optimized set of images for routine diagnosis [10]. The image quality of DECT was previously evaluated in comparison to that of single-energy CT (SECT) images in routine diagnostic interpretation by Yu et al. and Alvarez and Macovski [7, 10]. They evaluated the image quality of linearly mixed images from low- (80 kV) and high-energy (140 kV) images after image reconstruction at the same total radiation dose, taking into account the effect on patient size. The results indicated that the iodine contrast-to-noise ratio (CNR) in DECT images was similar to or better than that in SECT images.
Thus, DECT techniques are useful for gaining material information by using different linear attenuation coefficients. Usually, the tube voltages are set as 80/140 kV when the DE mode is applied for a patient. Although the low energy is typically set as 80 kV because of the resulting improvement in contrast of a material, the noise level is increased at this voltage because of the increased absorption coefficients of the material. On the other hand, the high-energy scan (140 kV) provides images with lower material contrast with high-atomic-number (Z) materials, but the noise is reduced. Therefore, it is necessary to investigate different kV combinations to determine their effect on the image quality of virtual monochromatic images acquired by DE scans. Previous research on different kV combinations of 80/140 kV and 100/140 kV focused on spectral separation by means of additional filtration [11]. In their work, additional filtration for the high-energy spectrum was modified to maximize the separation between the two spectra. The results proved that the effect of filtration could allow improved DE imaging with the 100/140 kV protocol.
Prior studies focused on image quality for DECT with the 80/140 kV protocol or on spectral separation for two spectra. The reason for using 80 kV as the low energy in the DECT scan is to avoid spectral overlapping. A wider separation between spectra should improve material separation while allowing radiation dose to the patient. To achieve this, an additional filter is applied to the high energy when using 100/140 kV scanning. Therefore, it is necessary to investigate the image quality with DE scan modes such as 80/140 kV and 100/140 kV without an additional filter.
In the present study, we investigated the image quality of virtual monochromatic images synthesized from DECT scans at voltages of 80/140 kV and 100/140 kV without an additional filter at the same radiation dose. Then, the virtual monochromatic images obtained with the two scan modes were compared to SE images in terms of noise and CNR.
Materials and Methods
1. Image-based Virtual Monochromatic Imaging
Image-based virtual monochromatic imaging was reported in a previous study [10]. Virtual monochromatic images are created from reconstructed low- and high-energy images. Assuming the mass attenuation coefficients of the two basis materials at low and high energy are
( μ ρ ) i j , where j=L, H and i=1, 2, the virtual monochromatic image at energy E is given by
where L and H represent the low and high energy, respectively, and 1 and 2 represent the two basis materials. By solving the two linear equations, one obtains the mass density of the two basis materials.
where
Δ = ( μ ρ ) 1 L · ( μ ρ ) 2 H - ( μ ρ ) 1 H · ( μ ρ ) 2 L .
The virtual monochromatic image at energy E is given by
Rewriting the linear attenuation coefficients in Equation 4 in terms of the CT number and assuming that one of the basis materials is water and the other is cortical bone, one can show that the virtual monochromatic image at energy E can be expressed as a weighted average of the images at low- and high-energy scans, which is given by
where the weighting factor is given by
Thus, the virtual monochromatic image generated from image-domain data is a linear combination of the two CT images at low and high energies, where the sum of the two weighting factors equals 1.
An optimal monochromatic energy exists that yields either the highest CNR or the lowest noise in the virtual monochromatic image. We define the standard deviations of CT numbers in the background and signal regions for low- and high-energy images as σj,b and σj,s, respectively, where j=L, H. The contrast (the absolute difference in CT numbers between the signal region and background region) for the low- and high-energy images is denoted by Cj (j=L, H). According to the definition in a previous paper [11], the weighting factor for the highest CNR in the linearly mixed image and the corresponding maximum CNR are given by
and
respectively.
The optimal weighting factor for either the lowest noise (Equation 9) or the highest CNR (Equation 10) corresponds to an optimal monochromatic energy, which can be obtained by solving either
for the minimum noise or
for the maximum CNR.
2. Phantom Information and Image Scan Setup
As shown in Figure 1, a tissue-characterization phantom (Gammex 467, Gammex Inc., Middleton, WI) was scanned by a CT scanner (Brilliance 64, Philips, Eindhoven, Holland) in both the SE and DE mode. Gammex 467 phantom has a diameter of 33 cm and consists of tissue substitutes such as lung exhale, adipose, breast, solid water, brain, liver, inner bone, and cortical bone. The specifications of the inserts are listed in Table 1. The phantom is scanned in the SE and DE mode in the step and shoot mode. The gantry rotation time was 1 second, and the slice thickness was 10 mm for all scans. SE scans were performed at 80, 100, 120, and 140 kV. We utilized DE scans of the 80/140 kV protocol and 100/140 kV protocol to obtain virtual monochromatic images.
3. Dose Partitioning
Computed tomography dose index (CTDI) is used for dose measurement. A CTDI obtained using a 100-mm chamber is referred to as CTDI100, which is calculated using the following equation:
where M is the measured charge in roentgens, F is the exposure-to-dose conversion factor (0.78 cGy/R), TP is the correction factor for both temperature and pressure, CF is the chamber calibration factor, N is the number of slices, and T is the slice thickness (mm). The weighted CTDI (CTDIw) can be obtained with a weighted CTDI100 by using the following equation:
where CTDI100,C is the CTDI100 measured at the center and CTDI100,P is the mean of the CTDI100 measured at the 3, 6, 9, and 12 o’clock positions. Since the pitch is 1 in one scan, CTDIw is equal to CTDIvol. A detailed description of the derivation of CTDI was reported in a previous publication [12].
4. Image Metrics
Physical metrics that influence image quality are evaluated with respect to CNR. For each DE image (80/140 kV and 100/140 kV), the CT numbers of the materials (lung exhale, adipose, breast, solid water, brain, liver, inner bone, and cortical bone) and the acryl background were measured. Virtual monochromatic images were generated from 40 to 70 keV by using Equation 2. For DE and SE images, noise is estimated by measuring the standard deviation of the region-of-interest (ROI) of the signal and background in the images of the materials. CNR is given by
where S is the signal of the materials and B is the signal of the background. σs and σb represent the noise expressed by standard deviation. CNR is investigated for virtual monochromatic images from 40 to 70 keV and for SE images at 80, 100, 120, and 140 kV. The maximum CNR in the virtual monochromatic images is calculated using Equation 6, and the results were compared to that in the SE images.
Results and Discussion
1. Noise Evaluation
Noise is measured as a function of energy through the standard deviation of the ROI on various materials in virtual monochromatic CT images and SECT images. The trend of noise values with respect to energy is illustrated in Figure 2.
As shown in Figure 2A–2D, the noise in DE images with the 80/140 kV and 100/140 kV protocol is evaluated with respect to the virtual monochromatic energy range. Here, the virtual monochromatic images are obtained using Equation 2. For all inserts (lung exhale, adipose, breast, solid water, brain, liver, inner bone, and cortical bone), Figure 2A and 2B show plots of the noise as a function of virtual monochromatic energy with the 80/140 kV protocol. Image noise is minimized for virtual monochromatic energy ranging from 47 to 51 keV for various materials, as shown in Figure 2A and 2B. In Figure 2C and 2D, image noise is plotted for the 100/140 kV protocol. The minimum noise is observed within the energy range from 51 to 54 keV. In the SE scans, the minimum noise in the images for all inserts is observed at 120 and 140 kV. Generally, CT images are obtained with scans at a tube voltage of 120 kV. Therefore, the quality of DE images is compared to that of SE images at 120 kV. In Figure 2E and 2F, image noise is plotted for the SE scans at 80, 100, 120, and 140 kV. The mean energies of the SE scan are 41.8, 47.7, 52.5, and 56.5 keV for 80, 100, 120, and 140 kV, respectively.
The minimum noise level in virtual monochromatic images and the image noise in SE images at 120 kV are illustrated in Figure 3. Noise in DE images with the 100/140 kV protocol is less than that with the 80/140 kV protocol for all materials. Consequently, the low-energy images at 80 kV had greater noise [14–16]. When the low energy is set to 100 kV, the image noise is reduced. Therefore, the 100/140 kV scan is better than the 80/140 kV scan for achieving a reduced noise. Although the noise in DE images with the 100/140 kV protocol was reduced up to 2.63%, 1.25%, and 3.27% for lung exhale, brain, and cortical bone, the noise in SE images at 120 kV is low for adipose, breast, solid water, liver, and inner bone.
2. CNR Trend According to Energy for DE and SE images
Figure 4 shows CNR values as a function of virtual monochromatic energy for DE and SE images. DE scans with the 80/140 kV and 100/140 kV protocol are used to plot CNRs from 40 to 70 keV. SE scans are performed at 80, 100, 120, and 140 kV. In DE images, CNR is maximized in the energy range of 47–53 keV for all materials. These trends are thought to be at the minimum noise level at the optimal monochromatic energies, and this phenomenon has already been reported in previous work [10]. Although a maximum CNR exists, the CNR in SE images is higher than that in DE images. However, the maximum CNR of breast and brain in DE images is higher than that in SE images. The CNR trend of brain is different between DE and SE images because the Hounsfield unit (HU) of brain is similar to the HU of the background.
3. Maximum CNR
The maximum CNR for various materials in DE images is compared to the CNR in SE scans at 120 kV, as illustrated in Figure 5, which is derived from Equations 5 and 6. The CNR of SE images at 120 kV is measured and compared to the maximum CNR of DE images with both the 80/140 kV and 100/140 kV scans. For comparison of CNR between DE images, the results of a scan at 100/140 kV are higher than that at 80/140 kV, except for adipose and brain. As shown in Figure 5, CNR improvement in the 100/140 kV scan is 6.75%, 1.23%, 5.33%, 8.84%, 3.38%, and 2.00% for lung exhale, breast, solid water, liver, inner bone, and cortical bone. Therefore, the 100/140 kV scan is better than the 80/140 kV scan for achieving an improved CNR. The reason for this effect is that noise level is reduced when using 100 kV as the low energy.
In DECT, mixed images are created from low- and high-energy images to provide a single set of images for routine diagnosis [17]. The virtual monochromatic image is synthesized from the sum of low- and high-energy images using linear weightings. In this process, noise could be increased by summation. For comparison between the CNRs of DE and SE, we found that the CNR of SE images for various materials is slightly higher than that of DE images owing to the synthesis process of virtual monochromatic images. In addition, the effect of linear attenuation coefficients on the X-ray energy is caused by the CNR values. If 100 kV is applied as the low-energy tube voltage, image noise is reduced, and contrast between the signal and background is low. Despite this result, the CNR improvement with the 100/140 kV protocol is 7.89% and 13.33% for the breast and brain compared to the corresponding SE images.
On comparing the virtual monochromatic images obtained using the 80/140 kV and 100/140 kV protocols, the 100/140 kV protocol could be considering for application to lung exhale, solid water, liver, inner bone, and cortical bone CT in accordance with the CNR results because the CNR of DE images scanned with the 100/140 kV protocol is higher than the CNR of 80/140 scan images.
The present work demonstrates that kV combinations affect the DE virtual monochromatic image. The effect of technical factors including combinations of low and high energy was considered in this evaluation, which provides essential information for physicists, technologists, and imaging physicians regarding the image quality at an optimal polychromatic tube potential for the same radiation dose.
There are some limitations of this study. We only assessed the image quality with respect to kV combinations in virtual monochromatic CT. Image quality depends on the radiation dose to the patient as well as the patient size. Further studies are needed to understand how body size and dose partitioning can influence image quality.
Conclusion
The image quality of virtual monochromatic images from DECT scan was evaluated and compared to that of SECT at the same radiation dose. The image quality of virtual monochromatic images mixed from different DECT scans was evaluated and compared to that of SECT scans at the same radiation dose. Overall, at optimal monochromatic energies, the maximum CNR with the 100/140 kV protocol was better than that with the 80/140 kV protocol. Furthermore, the maximum CNR with the 100/140 kV protocol is similar or higher than that of the SE scan at 120 kV.