Research Article  Open Access
Peng He, Biao Wei, Peng Feng, Mianyi Chen, Deling Mi, "Material Discrimination Based on Kedge Characteristics", Computational and Mathematical Methods in Medicine, vol. 2013, Article ID 308520, 6 pages, 2013. https://doi.org/10.1155/2013/308520
Material Discrimination Based on Kedge Characteristics
Abstract
Spectral/multienergy CT employing the stateoftheart energydiscriminative photoncounting detector can identify absorption features in the multiple ranges of photon energies and has the potential to distinguish different materials based on Kedge characteristics. Kedge characteristics involve the sudden attenuation increase in the attenuation profile of a relatively high atomic number material. Hence, spectral CT can utilize material Kedge characteristics (sudden attenuation increase) to capture images in available energy bins (levels/windows) to distinguish different material components. In this paper, we propose an imaging model based on Kedge characteristics for maximum material discrimination with spectral CT. The wider the energy bin width is, the lower the noise level is, but the poorer the reconstructed image contrast is. Here, we introduce the contrasttonoise ratio (CNR) criterion to optimize the energy bin width after the Kedge jump for the maximum CNR. In the simulation, we analyze the reconstructed image quality in different energy bins and demonstrate that our proposed optimization approach can maximize CNR between target region and background region in reconstructed image.
1. Introduction
Xray computed tomography (CT) has been widely applied in clinical and preclinical applications, since Hounsfield’s Nobel Prize winning breakthrough. A typical conventional CT system employs a broad energy spectrum source and a digital integrating sensor whose output is proportional to the energy fluence integrated over the entire incidence spectrum. Physically, the Xray spectrum contains much information; the conventional CT system collects photons over the whole Xray spectrum to ignore spectral responses of materials. Hence, the conventional CT often does not have sufficiently high contrast resolution for biological soft tissues [1].
With the development of spectral detectors and novel contrast agents, CT image contrast resolution could be significantly improved. Recent advances in spectral/multienergy detector technology have allowed for spectral CT systems to identify absorption features in the multiple ranges of photon energies [2–6]. Spectral CT has a stronger capability to distinguish different materials because it can capture images in available energy bins [7–13]. Meanwhile, contrast agent has been widely applied in biomedical imaging to enhance tissue contrast [14–19]. Spectral CT imaging utilizes not only density characteristics of contrast agents but also Kedge characteristics of contrast agents to distinguish different materials. Kedge characteristics involve the sudden attenuation increase in the attenuation profile of some contrast agents, which could be captured by spectral CT in available energy bins. Hence, different materials can be easily distinguished according to their Kedges characteristics [20, 21], while their Hounsfield numbers may be very similar in conventional CT images. This opens a door for spectral CT to support functional, cellular, and molecular imaging studies.
For contrast agent imaging by spectral CT, threshold settings for available energy bins have a major impact on spectral image quality in terms of image contrast and noise level. Hence, it is important to partition the energy bin optimally based on Kedge characteristics. In this paper, we propose a contrast agent imaging model to optimize energy bin for maximum material discrimination with spectral CT. Based on this model, we investigate how to select one energy bin for optimal contrast agent imaging to distinguish different materials, introducing a contrasttonoise ratio (CNR) where the signal difference is defined between contrast enhancement region (CER) and background region values.
This paper is organized as follows. Section 2 introduces our proposed contrast agent imaging model. Section 3 describes our simulation experiment. Section 4 demonstrates our experimental results. Section 5 discusses relevant issues and concludes the paper.
2. Materials and Methods
An earlier paper [22] has introduced the theoretical formalism of Kedge imaging model to determine two energy bins on both sides of the Kedge and analyzed the effect of Kedge energy bins on the resultant image quality. Here, we only briefly reproduce the optimization scheme and analyze how to optimize one energy bin after the Kedge jump to distinguishing contrast enhancement region and background region.
For current spectral CT system, its spectral detector (e.g., Medipix3 [4–6]) is a photoncounting system with selectable thresholds, which depends on a threshold equalization mask to adjust each pixel to record different energy photons. We assume that the energy distribution function of an Xray source is ; in a given energy threshold , we have the photon number received by the spectral detector: where is the detector efficiency.
For a given energy bin defined by two energy thresholds , the received photon number can be expressed as This paper focuses on how to set the thresholds in the energy bin imaging to distinguish contrast enhancement region and background region based on Kedge characteristics. First, we need to study the linear attenuation characteristics of background materials and contrast agents. Let be the linear attenuation coefficient function of a background material at an energy ; we have we assume that is the linear attenuation coefficient function of a contrast agent at an energy , and can be expressed as where and are the mass attenuation coefficients of background material and contrast agent, and are the densities of background material and contrast agent, respectively, and is the concentration of the contrast agent.
For a given contrast agent concentration , we plot two linear attenuation profiles of a typical contrast agent and background material (i.e., tissue) according to their mass attenuation coefficients and densities, which are shown as in Figure 1. In Figure 1, the attenuation coefficients of the contrast agent have a sudden increment at an energy , which reflects the Kedge characteristics. Theoretically, if we perform an energy bin imaging at the point with spectral CT, there is a maximum material discrimination for background region and contrast agent region in reconstructed image. However, the narrower the energy bin width is, the higher the noise level is, and there are few photons to carry the information. Hence, we select one energy bin of finite width to study contrast agent imaging in this paper. Let be the average attenuation coefficient of the background material within the energy bin after the Kedge jump; we have where is the energy bin width, and let be the average attenuation coefficient of the contrast agent within the energy bin after the Kedge jump; we have
For a given energy bin after the Kedge jump, from (2) we have the received photon number The reconstructed images can be evaluated as contrasttonoise ratio (CNR), and the CNR can be defined as where and are corresponding variances of the contrast agent region and background region in reconstructed image.
The difference between the mean value of contrast agent region and background region relies on the energy bin width and contrast agent concentration . For a given reconstructed object, the variances of reconstructed image rely on the photon number determined by the energy bin width . Hence, we can search for the optimal value to maximize the CNR for maximum material discrimination. In the following, we will make numerical simulation to test the proposed imaging model, including phantom design and image reconstruction protocols.
3. Numerical Simulation
In the simulation, a thorax phantom (Figure 2) was designed to be more preclinically relevant, which is defined on http://www.imp.unierlangen.de/forbild/. The phantom contains a heart region, a tissue region, a lung region, a vertebra region, and a contrast enhancement region (CER). The phantom was made 25 cm × 25 cm in size and discretized into a 500 × 500 matrix. We used Gadolinium solution whose Kedge is 50 keV, as a testing contrast agent in the CER inside the heart region, and the whole heart region is considered as the region of interest (ROI).
To investigate the proposed imaging theory, we study how to search for the optimal energy bin to maximize the CNR for maximum material discrimination. First, we study the characteristics of the thorax phantom materials. For tomographic imaging, the linear attenuation coefficient represents the gray value of reconstructed image. Here, we can obtain the mass attenuation coefficient according to the Xray attenuation databases reported by the National Institute of Standards and Technology (NIST). To calculate the linear attenuation coefficients of the phantom materials, the densities of these materials were selected in reference to the biomedical literature [23–26] and summarized in Table 1. In our simulation, we used blood attenuation characteristics to substitute heart attenuation characteristics.

Then, we used a freeofcharge software program (SpekCalc) [27] to calculate Xray spectra from tungsten anode tubes. The Xray tube voltage is assumed as 120 kVp with a 2.5 mm Al filter, and its emission spectra are shown in Figure 3. In our study, it was assumed that the detector efficiency was 90%, the spectral CT system was viewed as in a typical parallelbeam geometry, and the scanning range was from 0° to 180° with a 1° angular increment in the given energy bins. According to the BeerLambert law, we can capture the thorax phantom projection data.
To perform the energy bin imaging with spectral CT, we used a typical analytical reconstruction protocol: filtered backprojection (FBP), and a reconstructed image using FBP formula can be expressed as where is the integral of the linear attenuation coefficient distribution along an Xray path.
Then, we can calculate the expected image and the variance of reconstructed image [22, 28] From (10) and (11), CNR of ROI in reconstructed image can be written as follows: where is the reconstructed sinogram of the contrast agent region and is the reconstructed sinogram of the background material region. From (12), we can find that the CNR will depend on the energy bin width . In Section 4, we will calculate the best energy bin width value after the Kedge jump to maximize the CNR of ROI.
4. Results
We used our proposed approach to analyze the thorax phantom and determined the best energy bin for imaging based on Kedge characteristics and plotted the relationships between and CNR, as shown in Figure 4. Then, we used the optimal value (29 keV) to perform energy bin imaging for the thorax phantom, and the reconstructed image is shown in Figure 5(a). Meanwhile, we chose a broad energy spectrum (25~100 keV) to reconstruct the thorax phantom which can be considered as the conventional CT imaging, and the reconstructed result is shown in Figure 5(b). Compared to the broad energy spectrum imaging result, it is easier to distinguish Gadolinium solution region and heart region in optimal energy bin imaging result. Finally, we calculated the CNR of ROI in Figures 5(a) and 5(b), which is summarized in Table 2. From Table 2, we can see that CNR of ROI in Figure 5(a) is smaller than that in Figure 5(b). To compare the reconstructed results, we plotted the profiles along the broken lines in Figure 5, as shown in Figure 6.

(a)
(b)
The CNR of ROI also relies on the contrast agent concentration , and then we analyze the relationships between CNR and the concentration of contrast agent. We chose different concentrations of Gadolinium solution (0.5%, 1%, and 5%) as the testing contrast agents and plotted the relationships between and CNR with different concentrations, as shown in Figure 7. Finally, we calculated the optimal and maximum CNR for different concentration contrast agents, and the results are summarized in Table 3. From Table 3, we can see that the higher the concentration of Gadolinium solution is, the wider the optimal energy bin width is, and the bigger the CNR of ROI is.

5. Discussions and Conclusion
This paper is a followup study for an earlier paper [22]. Although some relevant theories are similar, this study focuses on how to distinguish contrast agents and background materials (i.e., tissue) in biomedical imaging with spectral CT, which can be readily generalized to deal with more general settings and able to determine the best energy bin for maximum material discrimination.
There are several issues worth further discussion in the simulation. First, we apply a mimetic Xray emission spectrum in our study; the Xray emission spectra with 1 keV energy bins are obtained by the freeofcharge software program (SpekCalc). To analyze relationships between energy bin and CNR of ROI, we calculate the photon number with the given energy bin width value of an integer. As a result, it is inevitable to introduce some errors, compromising the estimation of the optimal energy bin width . Second, the curves of CNR in Figure 7 are not so smooth in some energy bins, which reflects the real CNR characteristics in the given Xray emission spectra. The unsmooth reason is that the Xray emission spectra have some drastic jumps in some energies. If the imaging energy bin contains these drastic jumps, the photon number will drastically increase, and the variance of ROI in the reconstructed image will drastically decrease, so the curves of CNR in Figure 7 have some jumps in corresponding energy bins. Additionally, the proposed approach depends on specific phantom configurations, and the optimal energy bin width is applicationspecific. Nevertheless, our optimization theory is rigorous and can be applied once the application context or the class of images is known. In a followup study, we will study the biomedical samples with the spectral CT based on our proposed imaging theory.
In conclusion, contrast agent imaging with spectral CT has a great potential for clinical applications including, but not limited to, tissue characterization and contrast studies. We proposed a contrast agent imaging model to optimize energy bin for maximum material discrimination; it established guidelines for optimization of energy thresholds and could be readily generalized for biomedical imaging.
Acknowledgment
This work was partially supported by the National Natural Science Foundation of China Grants 61171157 and 61201346.
References
 G. Wang, H. Yu, and B. De Man, “An outlook on Xray CT research and development,” Medical Physics, vol. 35, no. 3, pp. 1051–1064, 2008. View at: Publisher Site  Google Scholar
 M. Campbell, E. H. M. Heijne, G. Meddeler, E. Pernigotti, and W. Snoeys, “A Readout Chip for a 64 × 64 Pixel Matrix with 15bit Single Photon Counting,” IEEE Transactions on Nuclear Science, vol. 45, no. 3, pp. 751–753, 1998. View at: Google Scholar
 X. Llopart, M. Campbell, R. Dinapoli, D. San Segundo, and E. Pernigotti, “Medipix2: a 64 k pixel readout chip with 55μm square elements working in single photon counting mode,” IEEE Transactions on Nuclear Science, vol. 49, no. 5, pp. 2279–2283, 2002. View at: Publisher Site  Google Scholar
 R. Ballabriga, M. Campbell, E. Heijne, X. Llopart, L. Tlustos, and W. Wong, “Medipix3: a 64 k pixel detector readout chip working in single photon counting mode with improved spectrometric performance,” Nuclear Instruments and Methods in Physics Research A, vol. 633, no. 1, pp. S15–S18, 2011. View at: Publisher Site  Google Scholar
 R. Ballabriga, G. Blaj, M. Campbell et al., “Characterization of the Medipix3 pixel readout chip,” Journal of Instrumentation, vol. 6, no. 1, pp. 1052–1058, 2011. View at: Publisher Site  Google Scholar
 R. Ballabriga, M. Campbell, E. H. M. Heijne, X. Llopart, and L. Tlustos, “The Medipix3 prototype, a pixel readout chip working in single photon counting mode with improved spectrometric performance,” IEEE Transactions on Nuclear Science, vol. 54, no. 5, pp. 1824–1829, 2007. View at: Publisher Site  Google Scholar
 J. Giersch, M. Firsching, D. Niederlöhner, and G. Anton, “Material reconstruction with spectroscopic pixel Xray detectors,” Nuclear Instruments and Methods in Physics Research A, vol. 546, no. 12, pp. 125–130, 2005. View at: Publisher Site  Google Scholar
 M. Firsching, P. T. Talla, T. Michel, and G. Anton, “Material resolving Xray imaging using spectrum reconstruction with Medipix2,” Nuclear Instruments and Methods in Physics Research A, vol. 591, no. 1, pp. 19–23, 2008. View at: Publisher Site  Google Scholar
 E. C. Frey, X. Wang, Y. Du et al., “Investigation of the use of photon counting Xray detectors with energy discrimination capability for material decomposition in microcomputed tomography,” Physics of Medical Imaging, vol. 6510, pp. 93–103, 2007. View at: Google Scholar
 S. J. Nik, J. Meyer, and R. Watts, “Optimal material discrimination using spectral Xray imaging,” Physics in Medicine and Biology, vol. 56, no. 18, pp. 5969–5983, 2011. View at: Publisher Site  Google Scholar
 B. Norlin, A. Manuilskiy, H.E. Nilsson, and C. Fröjdh, “Material recognition with the Medipix photon counting colour Xray system,” Nuclear Instruments and Methods in Physics Research A, vol. 531, no. 12, pp. 265–269, 2004. View at: Publisher Site  Google Scholar
 H. Q. Le and S. Molloi, “Least squares parameter estimation methods for material decomposition with energy discriminating detectors,” Medical Physics, vol. 38, no. 1, pp. 245–255, 2011. View at: Publisher Site  Google Scholar
 M. Firsching, D. Niederlöhner, T. Michel, and G. Anton, “Quantitative material reconstruction in CT with spectroscopic Xray pixel detectors—A Simulation Study,” in Proceedings of the IEEE Nuclear Science Symposium Conference Record, vol. 4, pp. 2257–2259, November 2006. View at: Publisher Site  Google Scholar
 N. G. Anderson, A. P. Butler, N. J. A. Scott et al., “Spectroscopic (multienergy) CT distinguishes iodine and barium contrast material in MICE,” European Radiology, vol. 20, no. 9, pp. 2126–2134, 2010. View at: Publisher Site  Google Scholar
 E. Roessl, B. Brendel, K.J. Engel, J.P. Schlomka, A. Thran, and R. Proksa, “Sensitivity of photoncounting based KEdge Imaging in Xray computed tomography,” IEEE Transactions on Medical Imaging, vol. 30, no. 9, pp. 1678–1690, 2011. View at: Publisher Site  Google Scholar
 A. Abudurexiti, M. Kameda, E. Sato et al., “Demonstration of iodine Kedge imaging by use of an energydiscrimination Xray computed tomography system with a cadmium telluride detector,” Radiological Physics and Technology, vol. 3, no. 2, pp. 127–135, 2010. View at: Publisher Site  Google Scholar
 Q. Xu, H. Y. Yu, J. Bennett et al., “Image reconstruction for hybrid truecolor microCT,” IEEE Transactions on Biomedical Engineering, vol. 59, pp. 1711–1719, 2012. View at: Google Scholar
 P. He, H. Y. Yu, P. Thayer et al., “Preliminary experimental results from a MARS microCT system,” Journal of XRay Science and Technology, vol. 20, pp. 199–211, 2012. View at: Google Scholar
 M. Firsching, A. P. Butler, N. Scott, N. G. Anderson, T. Michel, and G. Anton, “Contrast agent recognition in small animal CT using the Medipix2 detector,” Nuclear Instruments and Methods in Physics Research A, vol. 607, no. 1, pp. 179–182, 2009. View at: Publisher Site  Google Scholar
 S. J. Riederer and C. A. Mistretta, “Selective iodine imaging using Kedge energies in computerized Xray tomography,” Medical Physics, vol. 4, no. 6, pp. 474–481, 1977. View at: Publisher Site  Google Scholar
 J. P. Schlomka, E. Roessl, R. Dorscheid et al., “Experimental feasibility of multienergy photoncounting Kedge imaging in preclinical computed tomography,” Physics in Medicine and Biology, vol. 53, no. 15, pp. 4031–4047, 2008. View at: Publisher Site  Google Scholar
 P. He, B. Wei, W. X. Cong, and G. Wang, “Optimization of Kedge Imaging with Spectral CT,” Medical Physics, vol. 39, pp. 6572–6579, 2012. View at: Google Scholar
 H. HinghoferSzalkay and J. E. Greenleaf, “Continuous monitoring of blood volume changes in humans,” Journal of Applied Physiology, vol. 63, no. 3, pp. 1003–1007, 1987. View at: Google Scholar
 J. R. Cameron, J. G. Skofronick, and R. M. Grant, Physics of the Body, Medical Physics, 1999.
 A. I. Saito, J. G. Li, C. Liu, K. R. Olivier, and J. F. Dempsey, “Accurate heterogeneous dose calculation for lung cancer patients without highresolution CT densities,” Journal of Applied Clinical Medical Physics, vol. 10, no. 2, pp. 92–103, 2009. View at: Google Scholar
 N. H. Clinthorne, “A constrained dualenergy reconstruction method for materialselective transmission tomography,” Nuclear Instruments and Methods in Physics Research A, vol. 353, no. 1–3, pp. 347–348, 1994. View at: Google Scholar
 G. Poludniowski, G. Landry, F. Deblois, P. M. Evans, and F. Verhaegen, “SpekCalc: a program to calculate photon spectra from tungsten anode Xray tubes,” Physics in Medicine and Biology, vol. 54, no. 19, pp. N433–N438, 2009. View at: Publisher Site  Google Scholar
 A. C. Kak and M. Slaney, Principles of Computerized Tomographic Imaging, IEEE Press, 1988.
Copyright
Copyright © 2013 Peng He et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.