Table of Contents Author Guidelines Submit a Manuscript
Computational and Mathematical Methods in Medicine
Volume 2018, Article ID 2527516, 9 pages
Research Article

Image Decomposition Algorithm for Dual-Energy Computed Tomography via Fully Convolutional Network

1National Digital Switching System Engineering & Technological R&D Centre, Zhengzhou 450002, China
2153 Central Hospital of Henan Province, Zhengzhou 450002, China

Correspondence should be addressed to Bin Yan; moc.liamtoh@ecapsby

Received 13 April 2018; Revised 17 July 2018; Accepted 30 July 2018; Published 5 September 2018

Academic Editor: Maria E. Fantacci

Copyright © 2018 Yifu Xu 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.


Background. Dual-energy computed tomography (DECT) has been widely used due to improved substances identification from additional spectral information. The quality of material-specific image produced by DECT attaches great importance to the elaborated design of the basis material decomposition method. Objective. The aim of this work is to develop and validate a data-driven algorithm for the image-based decomposition problem. Methods. A deep neural net, consisting of a fully convolutional net (FCN) and a fully connected net, is proposed to solve the material decomposition problem. The former net extracts the feature representation of input reconstructed images, and the latter net calculates the decomposed basic material coefficients from the joint feature vector. The whole model was trained and tested using a modified clinical dataset. Results. The proposed FCN delivers image with about 60% smaller bias and 70% lower standard deviation than the competing algorithms, suggesting its better material separation capability. Moreover, FCN still yields excellent performance in case of photon noise. Conclusions. Our deep cascaded network features high decomposition accuracies and noise robust property. The experimental results have shown the strong function fitting ability of the deep neural network. Deep learning paradigm could be a promising way to solve the nonlinear problem in DECT.

1. Introduction

Conventional single-energy X-ray technique provides information about the examined object which is not sufficient to characterize it precisely. Dual-energy computed tomography (DECT) provides additional information by using two different energy spectra to scan the object, which has been presented as a valid alternative to conventional single-energy X-ray imaging. In recent years, the adoption of DECT has gained increased attention in public security [1] and medical field [2, 3]. The advantage of DECT is the ability for material characterization and differentiation [4]. This decomposition of mixture into two basic materials depends on the principle that the attenuation coefficient is material and energy dependent. Thus, measurements at two distinct energies should permit the separation of the attenuation into its basic components.

The quality of material-specific image produced by DECT attaches great importance to the elaborated design of the basis material decomposition method. The existing decomposition methods can be divided into two main categories: projection-based [57] and image-based [810]. Projection-based methods pass the projection data through a decomposition function, followed by image reconstruction such as filtered backprojection (FBP). It commonly provides better accuracy and reconstructed image with reduced beam-hardening artifacts in comparison with image-based methods. However, projection-based methods need matched projection datasets. This means that physically the same lines need to be measured for each spectrum, which is usually not the case in today’s CT scanners. Image-based methods use linear combinations of reconstructed images to get an image that contains material-selective DECT information. It is an approximative technique, and the resulting images are less quantitative than with projection-based methods. But image-based methods can handle mismatched projection datasets and are applicable to the decomposition of three or more constituent materials, which is more expedient in practice. Thus, they have been employed more frequently in modern DECT implementations.

The material decomposition problem in image domain can be described by the following equation:where and are the pixels in reconstructed images from low- and high-energy projections, respectively, and and are the corresponding points in decomposed basic materials images. The subscripts 1 and 2 indicate two specific materials. and are the average attenuation coefficients of the two basic materials under low/high-energy spectra. These attenuation coefficients are usually obtained by manually selecting two uniform regions of interest (ROIs) on the CT images that contain the basic materials [9, 11, 12]. Direct material decomposition via matrix inversion is a way of calculating the points and in the decomposed image, which is written as follows:

Equation (2) can be easily solved as long as the value of is not null. However, values of the two terms in do not differ significantly from each other. Therefore, the decomposition result is very sensitive to the noise in the input reconstructed images. Various methods have been proposed to solve this noise suppression problem. Precorrection [13, 14] methods reconstruct two water-precorrected images, followed by a linear combination, to yield images that are free from cupping artifacts usually in water-equivalent materials. The noise reduction techniques after image decomposition include Kalender’s correlated noise reduction (KCNR) [15, 16], noise forcing (NOF) [17], and noise clipping (NOC) [18], whose most fundamental strategy is the application of a smoothing filter. Recent advanced iterative methods [9, 10] consider the statistical properties of the decomposition process, producing high-quality edge-preserving images. These methods have shown great success on the decomposition problem. Their well performances rely on the well-handcrafted design of the algorithm.

In recent years, deep learning techniques, which use neural networks having a deep structure with three or more layers, have attracted widespread attention, mainly by outperforming alternative machine learning methods in numerous important applications. The current most popular deep model is the convolutional neural network (CNN) which has emerged as a powerful class of models for image classification [19, 20] and object detection [21]. In the field of computed tomography, some of the recent studies have already attempted to use deep neural networks to solve the problems such as low-dose image denoising [22] and artifact reduction [23]. Wang [24] provides an analytical and global perspective to the combination of tomographic imaging and deep learning. For the material decomposition problem in DECT, several neural network-based methods have also been proposed, but they all decompose the material in the projection domain [7, 25, 26].

Inspired by the recent learning-based methods [27, 28], in this paper, we propose an end-to-end image decomposition algorithm via deep learning techniques. A modified fully convolutional network is applied to extract the feature of reconstructed images and suppress the image noise at the same time. The last layer of the model is a fully connected layer to calculate the decomposed images from the extracted features. We demonstrate the effectiveness of our algorithm by the experiment on a clinical dataset. Two conventional algorithms are implemented and compared to the proposed FCN.

2. Methods

2.1. Fully Convolutional Network

Fully convolutional network (FCN) is one kind of CNN, which is firstly proposed and used for semantic segmentation [29]. The standard CNN generally is composed of a pooling layer and a convolutional layer which are alternately connected. The convolutional layers learn the features of the input. The pooling layers guarantee that the deeper layers can extract higher scale-level features through downsampling. In order to map the feature to the class labels, a fully connected layer will be added to the last output layer, which has fixed dimensions and throws away spatial coordinates. Due to this structural design, the naive CNN requires fixed-sized inputs and produces no-spatial outputs.

The main idea of FCN is transforming the last fully connected layer into a convolution layer with kernels that cover its entire input region. This replacement policy brings about several advantages for FCN. First, the input of the net can be the images of arbitrary sizes, which means that the net can be trained on image patches and then tested on the full-sized images. Second, it can efficiently learn to make dense predictions for per-pixel tasks such as semantic segmentation. Lastly, per-pixel tasks for naive CNN generate a huge amount of redundant convolution computations at adjacent patches. FCN avoids such problems by computing all convolutions in the first layer on the entire input image, leading to significant speedup in the forward-propagation process.

Because of these advantages, FCN is especially suitable for solving the image-based material decomposition problem which can also be regarded as a per-pixel prediction task. In addition, convolution operation to image is interpretable, since it can be seen as a kind of image filtering.

2.2. Image Decomposition Model

For image decomposition, we designed an end-to-end decomposition model based on FCN. The proposed model takes reconstructed images as inputs and predicts the basic material coefficients pixel by pixel in the decomposed image, completing image decomposition and noise suppression at one time.

An overview of our model is illustrated in Figure 1. It is composed of two types of layers: convolutional and fully connected layers. Since the pooling layer may discard important structural details in feature maps, we omit it from the model to avoid losing the quality of result images. But no downsampling process by the pooling layer will lead to the same size of the feature maps at different layers. We hope the model can still catch the multiscale features of the image at different layers, so the strides of the convolutional layers are set to 2 to finish the downsampling operation. The input of the model is the image patch of size in reconstructed images. There are two independent fully convolutional nets corresponding to the reconstructed images from low- and high-energy projections. The two nets have the same layer structure and are called the L-FCN and H-FCN in short in this study. They are composed of four convolutional layers. The output of layer can be formulated as follows:where is the input feature map or images and and represent the convolutional kernel weights and bias parameter, respectively. ∗ is the convolutional operation. is the nonlinear active function of the neuron. The outputs of L-FCN or H-FCN () are a vector which represents the feature of the current input patch. The two feature vectors from L-FCN and H-FCN are merged into a joint vector. Then, a fully connected layer is used to calculate the decomposed basic material coefficients from the joint vector, which follows the following equation:where is the predicted material coefficients vector, and are the unsolved parameter matrixes, and represents the merged vector from L-FCN and H-FCN.

Figure 1: Overall architecture of the proposed network.

The whole decomposed images can be obtained by traversing all the patches in the input images. The specific information about each layer of the proposed FCN is listed in Table 1.

Table 1: Detailed configuration of L-FCN/H-FCN.
2.3. The Training Detail

The proposed FCN is implemented via the TensorFlow [30] framework on a computer platform containing two Titan X GPUs (a total of 24 GB video memory). The base learning rate of the model is , which decays by an exponential power of 0.9. There are 1200 training samples in one batch. The mean squared error (MSE) is utilized as the loss function:where is the true value of the decomposed image. We used Adam [31] to optimize the loss function in this study. The entire model contains about 64k unsolved parameters and is trained for 40 epochs in 37 hours. The loss curve for training is plotted in Figure S1 in the Supplementary Materials.

3. Experimental Design

3.1. Experimental Dataset

The experimental data are acquired from a real clinical dataset which contains 5987 pleural and cranial cavity images from 12 patients. These raw data are obtained by one single-energy scan. The tissue and bone regions in the images are all manually sketched out. The images from 10 patients were selected to generate training samples, and the images from the rest of the patients were used for testing. All the images are split up into two partitions. Each partition includes regions of bone or tissue only and is used as the ground truth of the decomposed images. In order to generate dual-energy images, we processed the original raw data and simulated the imaging system. The original image is inconvenient to process for its small value. So, firstly, we amplified the value of the raw data to a proper range via a linear transform that follows the following equation:where and are the pixel values of tissue and bone regions in original images and and are the corresponding pixel values in transformed images. The values and in the experiment are set to 50 and 15, respectively. Here, the different setting of and is for the purpose of better visual contrast in the transformed images. Secondly, we applied a BM3D [32] algorithm for attenuation of additive white Gaussian noise from the image. Thirdly, we used SpekCalc [33] software to generate 80 kVp and 140 kVp energy spectra, calculated the projection under the simulated scan of dual energy, and obtained the reconstructed images via filtered backprojection (FBP). Lastly, for each patient in the training set, we selected one slice every 10 images. Then, for each image, we extracted patches with the sliding interval of 5 pixels. The patch size was set to , the same as the input layer of the proposed FCN, getting totally 2,454,300 training patches.

3.2. Evaluation Metrics

The proposed FCN is compared with two other algorithms, direct decomposition (matrix inversion) and iterative decomposition [9]. We choose the bias and standard deviation to evaluate the performance of these methods. Bias shows the difference between the measured value and expected value, which can be a measure of the precision of the result. Standard deviation (SD) reflects the degree of dispersion of the result. They are calculated as follows:where and are the predicted value and true value at point of the image, respectively, is the mean value of the material, and is the number of points in ROI.

In order to further investigate the robustness of the proposed FCN, before reconstruction via FBP, photon noise is introduced into the dual-energy projections. There are two major types of noise in X-ray projection images [34]. One type is due to the electrical and roundoff error, which is image independent and can be considered as the Gaussian noise; the other type is due to the statistical fluctuation of the X-ray photons, which is image dependent and can be considered as the Poisson noise. The noise of the first type is small and is omitted in this study. The noise of the second type can be calculated as follows:where and are the noise-corrupted low- and high-energy projections, is a random process according to Poisson’s distribution with mean , and and are the number of photons of low- and high-energy incident X-rays. We set and in the experiments, respectively.

4. Results

We test our model on a cranial image and a pleural image which are excluded from the training dataset. Figure 2 shows the decomposition results by using three algorithms. The first column is the ground truth. Bone and tissue are chosen as the basis materials. Matrix inversion achieves similar results in vision as iterative decomposition. Loss of details and noticeable blocky artifacts are observed for the tissue and bone images from both algorithms. Figure 3 shows the zoom-in images whose areas are indicated in Figure 2 with a dashed rectangle. The iterative decomposition delivers smooth image due to its smoothness regularization term in loss function. It is noticeable that the proposed FCN suppresses most artifacts while preserving the structural features better than the competing algorithms. But there are not distinct improvements in edge-preserving. We guess this is mainly caused by the convolution kernel in the model. The convolutional operation of image can be seen as a kind of filtering.

Figure 2: The decomposed images by using three methods.
Figure 3: Result comparisons in the zoom-in area which is indicated in Figure 2 with a dashed rectangle.

For quantitative evaluation, the bias and SD are calculated on the images generated by using different algorithms inside material’s ROI and summarized in Table 2. Generally, the estimate of bone is more accurate than that of tissue. The proposed FCN achieves results closest in values to the ground truth, with about 60% smaller bias and 70% lower standard deviation than the competing algorithms, suggesting its better material separation capability.

Table 2: A list of bias and SD on the images generated by using different algorithms.

To evaluate the potential improvement by FCN, we investigate the effects of photon noise on the material decomposition algorithms. The reconstructed image is generated from noise-corrupted projections as described in Section 3.2. Figure 4 presents the decomposition results on same testing images. It can be seen that direct matrix inversion magnifies the noises both in ROI and background. Iterative decomposition also suffers from serious artifacts. This indicates that both algorithms are more sensitive to the noise. The proposed FCN yields the decomposed images that have not much noticeable change in comparison with the results in Figure 2.

Figure 4: The decomposition results on data with photon noise.

Figure 5 illustrates the absolute value of the difference between images in Figures 2 and 4, providing a visual comparison of the performance of noise suppression. For matrix inversion, the noise is statistically independent and evenly distributed in the images because the value of each pixel in decomposed images is calculated by using the corresponding pixel in projections. For iterative decomposition, the noise demonstrates a regional distribution characteristic. The region of tissue and background contain larger amount of noises than bone. In contrast, there are not much obvious differences in the result produced by the proposed FCN. Clearly, it outperforms the other two algorithms, more effectively suppressing image noise while keeping subtle structures.

Figure 5: The absolute value of the difference between images in Figures 2 and 4.

The quantitative results are listed in Table 3. In the case of photon noise, the bias and SD of the competing algorithms have increased in varying degrees. FCN still demonstrates good agreement to the true value, indicating its advantages on the antinoise capability.

Table 3: A list of bias and SD on the images in case of photon noise.

5. Discussion

We have designed a cascaded neural network for the material decomposition problem. The reconstructed images are pixel wisely mapped to decomposed images via several convolutional layers and a fully connected layer. The size of the input layer is , based on the hypothesis that the value of the material coefficient depends largely on the local region in reconstructed images. The proposed FCN processes data in an end-to-end way, without any needs of precorrected images or other prior knowledge. The experimental results show its strong performance in capturing the localized structural information and suppressing image noise. The decomposed images generated by matrix inversion and iterative decomposition contain relatively a large amount of artifacts. In the robustness-testing experiment, the noise-corrupted inputs will have a negative impact on the performance of the other competing algorithms, but not much on the FCN. The proposed FCN still achieves excellent results which have low bias and standard deviation. Data augmentation was used in the training process. It brought no boost in performance but costs more training time. We guess the main reason for this issue is that the material decomposition is a regression problem. The value of the label is in a continuous space. Data augmentation assumes that the examples in vicinity share the same class. This hypothesis is usually plausible to the classification problem in which the label is a discrete variable, but unnecessary for the regression problem. The main drawback of our algorithm is the requirement of the specific type of material. Tissue and bone are selected as the basic material in the experiment. The whole model needs to be retrained if one of the materials was changed. So, we hope the proposed algorithm can be used in some applications such as medical diagnosis where the selection of the material is relatively fixed. The amount of training samples is another main factor contributing to the effectiveness of our model. Normally, more data bring better performance of the model. But it may be difficult to collect enough data in some conditions.

6. Conclusions and Further Work

In this study, we present a deep learning approach for the image decomposition problem in DECT. According to the preliminary decomposition results, we successfully prove the feasibility of the proposed algorithm which delivers image with 70% smaller bias and 60% lower standard deviation than the competing algorithms. A deep learning paradigm promises to improve the ability of solving the nonlinear problem in DECT.

We think there are two directions of work that are worth further researching. One is to extend our model to make it applicable to the three-materials decomposition problem. The other is the attempt of using the deconvolutional network which will output the whole decomposed images in a forward-propagation calculation rather than pixel wisely prediction.

Data Availability

The code and data used in the research can be obtained from

Conflicts of Interest

The authors have no conflicts of interest to declare.


This work was supported by the National Key R&D Program of China under grant no. 2017YFB1002502, the National Natural Science Foundation of China (nos. 61601518 and 61372172) and the National Natural Science Foundation of Henan Province of China (no. 162300410333).

Supplementary Materials

Figure S1: the proposed model contains about 64k unsolved parameters and is trained for 40 epochs in 37 hours. The training batch size is 1200 reconstructed images from the noise-corrupted low- and high-energy projections. Figures S2 and S3: more testing results to show the superiority of the proposed method. All the testing images are reconstructed from the noise-corrupted low- and high-energy projections. (Supplementary Materials)


  1. Z. Ying, R. Naidu, and C. R. Crawford, “Dual energy computed tomography for explosive detection,” Journal of X-Ray Science and Technology, vol. 14, no. 4, pp. 235–256, 2006. View at Google Scholar
  2. H. W. Goo and J. M. Goo, “Dual-energy CT: new horizon in medical imaging,” Korean Journal of Radiology, vol. 18, no. 4, pp. 555–569, 2017. View at Publisher · View at Google Scholar · View at Scopus
  3. C. Hong, T. Y. Chin, and W. C. G. Peh, “Dual-energy CT in gout-a review of current concepts and applications,” Journal of Medical Radiation Sciences, vol. 64, no. 1, pp. 41–51, 2017. View at Publisher · View at Google Scholar · View at Scopus
  4. A. A. Postma, M. Das, A. A. R. Stadler, and J. E. Wildberger, “Dual-energy CT: what the neuroradiologist should know,” Current Radiology Reports, vol. 3, no. 5, 2015. View at Publisher · View at Google Scholar
  5. B. Li and Y. Zhang, “Projection decomposition algorithm of x-ray dual-energy computed tomography based on projection matching,” Acta Optica Sinica, vol. 31, no. 3, Article ID 311002, 2011. View at Publisher · View at Google Scholar
  6. B. Brendel and J. P. Schlomka, “Empirical projection-based basis-component decomposition method,” in Proceedings of Medical Imaging 2009: Physics of Medical Imaging, vol. 7258, p. 72583Y, International Society for Optics and Photonics, Bellingham, WA, USA, March 2009. View at Publisher · View at Google Scholar
  7. T. G. Schmidt and K. C. Zimmerman, “Material decomposition of multi-spectral x-ray projections using neural networks,” 2015, US Patent US20150371378. View at Google Scholar
  8. J. W. Lambert, Y. Sun, R. G. Gould, M. A. Ohliger, Z. Li, and B. M. Yeh, “An image-domain contrast material extraction method for dual-energy computed tomography,” Investigative Radiology, vol. 52, no. 4, pp. 245–254, 2017. View at Publisher · View at Google Scholar · View at Scopus
  9. T. Niu, X. Dong, M. Petrongolo, and L. Zhu, “Iterative image-domain decomposition for dual-energy CT,” Medical Physics, vol. 41, no. 4, Article ID 041901, 2014. View at Publisher · View at Google Scholar · View at Scopus
  10. X. Dong, T. Niu, and L. Zhu, “Combined iterative reconstruction and image-domain decomposition for dual energy CT using total-variation regularization,” Medical Physics, vol. 41, no. 5, Article ID 051909, 2014. View at Publisher · View at Google Scholar · View at Scopus
  11. T. P. Szczykutowicz and G. H. Chen, “Dual energy CT using slow kVp switching acquisition and prior image constrained compressed sensing,” Physics in Medicine and Biology, vol. 55, no. 21, pp. 6411–6429, 2010. View at Publisher · View at Google Scholar · View at Scopus
  12. P. V. Granton, S. I. Pollmann, N. L. Ford, M. Drangova, and D. W. Holdsworth, “Implementation of dual- and triple-energy cone-beam micro-CT for postreconstruction material decomposition,” Medical Physics, vol. 35, no. 11, pp. 5030–5042, 2008. View at Publisher · View at Google Scholar · View at Scopus
  13. M. Baer, C. Maaß, W. A. Kalender, and M. Kachelrieß, “Image-based dual energy CT using optimized precorrection functions: a practical new approach to material decomposition in the image domain,” in World Congress on Medical Physics and Biomedical Engineering, pp. 205–208, Springer Berlin Heidelberg, Munich, Germany, September 2009.
  14. R. A. Brooks, “A quantitative theory of the Hounsfield unit and its application to dual energy scanning,” Journal of Computer Assisted Tomography, vol. 1, no. 4, pp. 487–493, 1977. View at Publisher · View at Google Scholar · View at Scopus
  15. W. A. Kalender, E. Klotz, and L. Kostaridou, “An algorithm for noise suppression in dual energy CT material density images,” IEEE Transactions on Medical Imaging, vol. 7, no. 3, pp. 218–224, 1988. View at Publisher · View at Google Scholar · View at Scopus
  16. D. L. Ergun, C. A. Mistretta, D .E. Brown et al., “Single-exposure dual-energy computed radiography: improved detection and processing,” Radiology, vol. 174, no. 1, pp. 243–249, 1990. View at Publisher · View at Google Scholar · View at Scopus
  17. J. T. Dobbins, “Recent progress in noise reduction and scatter correction in dual-energy imaging,” in Proceedings of SPIE-the International Society for Optical Engineering, pp. 134–142, San Diego, CA, USA, May 1995.
  18. W. W. Peppler, J. T. Dobbins, E. B. Bellers et al., “Dual-energy computed radiography: improvements in processing,” in Proceedings of SPIE-the International Society for Optical Engineering, vol. 2167, pp. 663–671, Newport Beach, CA, USA, May 1994.
  19. L. Kaiser, A. N. Gomez, N. Shazeer et al., “One model to learn them all,”, 2017.
  20. K. He, X. Zhang, S. Ren, and J. Sun, “Deep residual learning for image recognition,”, 2015.
  21. T. Kong, A. Yao, Y. Chen, and F. Sun, “HyperNet: towards accurate region proposal generation and joint object detection,” in Proceedings of Computer Vision and Patter Recognition IEEE, pp. 845–853, Piscataway, NJ, USA, April 2016.
  22. D. Wu, K. Kim, G. E. Fakhri et al., “A cascaded convolutional neural network for x-ray low-dose CT image denoising,”, 2017.
  23. H. Zhang, L. Li, K. Qiao et al., “Image prediction for limited-angle tomography via deep learning with convolutional neural network,”, 2016.
  24. G. Wang, “A perspective on deep imaging,” IEEE Access, vol. 4, no. 99, pp. 8914–8924, 2016. View at Publisher · View at Google Scholar · View at Scopus
  25. K. C. Zimmerman and T. G. Schmidt, “Experimental comparison of empirical material decomposition methods for spectral CT,” Physics in Medicine & Biology, vol. 60, no. 8, pp. 3175–3191, 2015. View at Publisher · View at Google Scholar · View at Scopus
  26. W. J. Lee, D. S. Kim, S. W. Kang, and W. J. Yi, “Material depth reconstruction method of multi-energy X-ray images using neural network,” in Proceedings of International Conference of the IEEE Engineering in Medicine and Biology Society, pp. 1514–1517, Piscataway, NJ, USA, August–September 2012.
  27. Y. Lu, M. Berger, M. Manhart et al., “Bridge to real data: empirical multiple material calibration for learning-based material decomposition,” in Proceedings of IEEE International Symposium on Biomedical Imaging, pp. 457–460, Prague, Czech Republic, April 2016.
  28. Y. Lu, J. Geret, M. Unberath et al., “Projection-based material decomposition by machine learning using image-based features for computed tomography,” in Proceedings of the 13th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, pp. 448–452, Newport, RI, USA, July 2015.
  29. E. Shelhamer, J. Long, and T. Darrell, “Fully convolutional networks for semantic segmentation,” in Proceedings of IEEE Transactions on Patter Analysis & Machine Intelligence, vol. 39, no. 4, pp. 640–651, Washington, DC, USA, 2014.
  30. M. Abadi, A. Agarwal, P. Barham et al., “TensorFlow: large-scale machine learning on heterogeneous distributed systems,”, 2016.
  31. D. P. Kingma and J. Ba, “Adam: a method for stochastic optimization,”, 2014.
  32. K. Dabov, A. Foi, and V. Katkovnik, “Image denoising by sparse 3-D transform-domain collaborative filtering,” IEEE Transactions on Image Processing, vol. 16, no. 8, pp. 2080–2095, 2007. View at Publisher · View at Google Scholar · View at Scopus
  33. 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 & Biology, vol. 54, no. 19, pp. N433–N438, 2009. View at Publisher · View at Google Scholar · View at Scopus
  34. L. Zhu, J. Wang, and L. Xing, “Noise suppression in scatter correction for cone-beam CT,” Medical Physics, vol. 36, no. 3, pp. 741–752, 2009. View at Publisher · View at Google Scholar · View at Scopus