Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 308520, 6 pages

http://dx.doi.org/10.1155/2013/308520

Research Article

Material Discrimination Based on K-edge Characteristics

Peng He,1,2 Biao Wei,¹ Peng Feng,¹ Mianyi Chen,¹ and Deling Mi¹

¹The Key Lab of Optoelectronic Technology and Systems of the Education Ministry of China, Chongqing University, Chongqing 400044, China

²The Key Lab of Biorheological Science and Technology of the Education Ministry of China, Chongqing University, Chongqing 400044, China

Received 25 July 2013; Revised 28 September 2013; Accepted 2 October 2013

Academic Editor: Liang Li

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.

Abstract.

Spectral/multienergy CT employing the state-of-the-art energy-discriminative photon-counting detector can identify absorption features in the multiple ranges of photon energies and has the potential to distinguish different materials based on K-edge characteristics. K-edge characteristics involve the sudden attenuation increase in the attenuation profile of a relatively high atomic number material. Hence, spectral CT can utilize material K-edge 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 K-edge 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 contrast-to-noise ratio (CNR) criterion to optimize the energy bin width after the K-edge 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

X-ray 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 X-ray spectrum contains much information; the conventional CT system collects photons over the whole X-ray 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 K-edge characteristics of contrast agents to distinguish different materials. K-edge 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 K-edges 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 K-edge 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 contrast-to-noise 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 K-edge imaging model to determine two energy bins on both sides of the K-edge and analyzed the effect of K-edge 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 K-edge jump to distinguishing contrast enhancement region and background region.

For current spectral CT system, its spectral detector (e.g., Medipix-3 [4–6]) is a photon-counting 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 X-ray 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 K-edge 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 K-edge 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 K-edge 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 K-edge jump; we have

Figure 1: Attenuation profiles of a typical contrast agent and a soft tissue (background material).

For a given energy bin after the K-edge jump, from (2) we have the received photon number

The reconstructed images can be evaluated as contrast-to-noise 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.uni-erlangen.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 K-edge 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).

Figure 2: Thorax phantom.

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 X-ray 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.

Table 1: Biologically relevant densities of phantom materials.

Materials
Blood
Bone
Lung
Tissue

Density (g/cm³)
1.05
1.9
0.26
1.0

Then, we used a free-of-charge software program (SpekCalc) [27] to calculate X-ray spectra from tungsten anode tubes. The X-ray 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 parallel-beam geometry, and the scanning range was from 0° to 180° with a 1° angular increment in the given energy bins. According to the Beer-Lambert law, we can capture the thorax phantom projection data.

Figure 3: Source photon emission spectra.

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 X-ray 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 K-edge 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 K-edge 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.

Table 2: Summary of maximum CNR for different imaging models.

Imaging models
CNR

Optimal energy bin (5079 keV) imaging
642.1

Broad energy spectrum (25100 keV) imaging
438.0

Figure 4: Relationship between the energy bin width () and CNR. The curve for Gadolinium solution (0.5%) in the thorax phantom.

(a)

(b)

Figure 5: Two model imaging results. (a) is the reconstructed thorax phantom image in the optimal energy bin, and (b) is the reconstructed thorax phantom image in a broad energy spectrum. The display window for the two images is .

Figure 6: Profiles corresponding to the broken lines in Figure 5.

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.

Table 3: Summary of the optimal width and maximum CNR for different concentration contrast agents.

Contrast agents
Optimal (keV)
Maximum CNR

Gadolinium solution (0.5%)
29
642.1

Gadolinium solution (1%)
30
1084.9

Gadolinium solution (5%)
32
1519.0

Figure 7: Relationships between the energy bin width () and CNR. The curves for Gadolinium solutions (0.5%, 1%, and 5%) in the thorax phantom.

5. Discussions and Conclusion

This paper is a follow-up 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 X-ray emission spectrum in our study; the X-ray emission spectra with 1 keV energy bins are obtained by the free-of-charge 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 X-ray emission spectra. The unsmooth reason is that the X-ray 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 application-specific. Nevertheless, our optimization theory is rigorous and can be applied once the application context or the class of images is known. In a follow-up 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 X-ray CT research and development,” Medical Physics, vol. 35, no. 3, pp. 1051–1064, 2008.

M. Campbell, E. H. M. Heijne, G. Meddeler, E. Pernigotti, and W. Snoeys, “A Readout Chip for a 64 × 64 Pixel Matrix with 15-bit Single Photon Counting,” IEEE Transactions on Nuclear Science, vol. 45, no. 3, pp. 751–753, 1998.

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.

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.

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.

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.

J. Giersch, M. Firsching, D. Niederlöhner, and G. Anton, “Material reconstruction with spectroscopic pixel X-ray detectors,” Nuclear Instruments and Methods in Physics Research A, vol. 546, no. 1-2, pp. 125–130, 2005.

M. Firsching, P. T. Talla, T. Michel, and G. Anton, “Material resolving X-ray imaging using spectrum reconstruction with Medipix2,” Nuclear Instruments and Methods in Physics Research A, vol. 591, no. 1, pp. 19–23, 2008.

E. C. Frey, X. Wang, Y. Du, et al., “Investigation of the use of photon counting X-ray detectors with energy discrimination capability for material decomposition in micro-computed tomography,” Physics of Medical Imaging, vol. 6510, pp. 93–103, 2007.

S. J. Nik, J. Meyer, and R. Watts, “Optimal material discrimination using spectral X-ray imaging,” Physics in Medicine and Biology, vol. 56, no. 18, pp. 5969–5983, 2011.

B. Norlin, A. Manuilskiy, H.-E. Nilsson, and C. Fröjdh, “Material recognition with the Medipix photon counting colour X-ray system,” Nuclear Instruments and Methods in Physics Research A, vol. 531, no. 1-2, pp. 265–269, 2004.

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.

M. Firsching, D. Niederlöhner, T. Michel, and G. Anton, “Quantitative material reconstruction in CT with spectroscopic X-ray pixel detectors—A Simulation Study,” in Proceedings of the IEEE Nuclear Science Symposium Conference Record, vol. 4, pp. 2257–2259, November 2006.

N. G. Anderson, A. P. Butler, N. J. A. Scott et al., “Spectroscopic (multi-energy) CT distinguishes iodine and barium contrast material in MICE,” European Radiology, vol. 20, no. 9, pp. 2126–2134, 2010.

E. Roessl, B. Brendel, K.-J. Engel, J.-P. Schlomka, A. Thran, and R. Proksa, “Sensitivity of photon-counting based K-Edge Imaging in X-ray computed tomography,” IEEE Transactions on Medical Imaging, vol. 30, no. 9, pp. 1678–1690, 2011.

A. Abudurexiti, M. Kameda, E. Sato et al., “Demonstration of iodine K-edge imaging by use of an energy-discrimination X-ray computed tomography system with a cadmium telluride detector,” Radiological Physics and Technology, vol. 3, no. 2, pp. 127–135, 2010.

Q. Xu, H. Y. Yu, J. Bennett, et al., “Image reconstruction for hybrid true-color micro-CT,” IEEE Transactions on Biomedical Engineering, vol. 59, pp. 1711–1719, 2012.

P. He, H. Y. Yu, P. Thayer, et al., “Preliminary experimental results from a MARS micro-CT system,” Journal of X-Ray Science and Technology, vol. 20, pp. 199–211, 2012.

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.

S. J. Riederer and C. A. Mistretta, “Selective iodine imaging using K-edge energies in computerized X-ray tomography,” Medical Physics, vol. 4, no. 6, pp. 474–481, 1977.

J. P. Schlomka, E. Roessl, R. Dorscheid et al., “Experimental feasibility of multi-energy photon-counting K-edge imaging in pre-clinical computed tomography,” Physics in Medicine and Biology, vol. 53, no. 15, pp. 4031–4047, 2008.

P. He, B. Wei, W. X. Cong, and G. Wang, “Optimization of K-edge Imaging with Spectral CT,” Medical Physics, vol. 39, pp. 6572–6579, 2012.

H. Hinghofer-Szalkay and J. E. Greenleaf, “Continuous monitoring of blood volume changes in humans,” Journal of Applied Physiology, vol. 63, no. 3, pp. 1003–1007, 1987.

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 high-resolution CT densities,” Journal of Applied Clinical Medical Physics, vol. 10, no. 2, pp. 92–103, 2009.

N. H. Clinthorne, “A constrained dual-energy reconstruction method for material-selective transmission tomography,” Nuclear Instruments and Methods in Physics Research A, vol. 353, no. 1–3, pp. 347–348, 1994.

G. Poludniowski, G. Landry, F. Deblois, P. M. Evans, and F. Verhaegen, “SpekCalc: a program to calculate photon spectra from tungsten anode X-ray tubes,” Physics in Medicine and Biology, vol. 54, no. 19, pp. N433–N438, 2009.

A. C. Kak and M. Slaney, Principles of Computerized Tomographic Imaging, IEEE Press, 1988.