Practical limitations during the image acquisition like obese patients often lead to the fact that the raw data is recorded incompletely. Thus, the quality of the reconstructed image is decreases by streaking artifacts which often prohibits usage in clinical routines.
As Rudin et al. showed in , the minimization of the total variation can be used to reduce artifacts and noise while edges in the image domain are preserved. Recently Chen et al.  extended this approach by introducing anisotropic filtering with directional derivatives to enhance the image quality and sharpen the edges.
Rudin, Osher and Fatemi outlined the benefits of Total variation (TV) or L1-norm of derivatives for denoising compared with other methods like L2-norm related minimization algorithms or maximum entropy methods. The total variation of an image f in the image domain is defined as(1)
where ||·||1 describes the L1-norm of the gradient of f. The fact that limited angle artifacts appear perpendicular to the direction of the missing projections motivates the idea of anisotropic regularization to reduce these streaks and their impact of disturbances. To minimize such an expression, one has to find the roots of the partial derivatives:(2)
Since discontinuities may arise using (2), Persson et al. suggest to use the approximation defined by
where [·]x, [·]y and [·]z denote the x, y or z-component of the given vector. Consequently, (2) can be written as
Thereby, ||·||ε defines the modified norm to stabilize the derivatives.
2.1 Anisotropic regularization
The investigated and extended algorithm in this paper is the algorithm for iterative anisotropic total variation minimization (aTVM) of Chen et al. . It was designed to remove artifacts that strongly depend on the direction of projection (by using directional derivatives). The fact that limited angle artifacts appear perpendicular to the direction of the missing projections encourage the use of anisotropic regularization to reduce these streaks. The anisotropic TV minimization can be described as an optimization problem constrained by the fact that the projection of the recovered image χ(f) should match the projections p.(3)
Here i denotes the ith direction under the angle αi used for the directional derivative while Nα denotes the total number of directions used. The directional derivatives are used to minimize the TV of the image f, while preserving fine details. A weighted and normalized edge image is calculated by(4)
and minimized by a gradient descent with the update parameter(5)
The starting and final regularization parameter sstart and sfinal, respectively, cause a cooling down effect of the strength of regularization to prevent oversmoothing. The weights ωi are defined by(6)
where Nv denotes the number of projections and rj(θ) the angle between the ray j and the rotated axis. In polar coordinates the weight for a given angle is defined by the distance towards the origin.
Finally, in an iterative reconstruction the weighted correction term
in the nth iteration is added to the previous image and the step parameter updated afterwards:(7)
For our results we used a scanning range of 160°, sstart = 1, sfinal = 10−4, NTV = 100, N = 3 and ε = 10−8 with equiangular distributed i.We used this scheme in an iterative reconstruction approach, the regularization is applied once an iteration is finished. Afterwards, the smoothed image is provided as a prior for the next iteration. Figure 1 shows the resulting images for all different methods under investigation. The corresponding similarity values are given in Table 1. The similarity is expressed in terms of the structural similarity index (SSIM) , which considers not only the gray values but also the structural information present in the image with floating windows or patches. Table 1 clearly states that the anisotropic method (SART + aTVM) outperforms the other algorithms. Our algorithm provides comparable results to the standard denoising approach, the ROF-denoising model, but also uses directional derivatives. Figure 2 shows the trend depending on the number of iterations. It can clearly be seen that the ROF-denoising and the anisotropic method give comparable results regarding the SSIM values. Moreover, both trend to even better results and clearly outperform the isotropic method.
This paper features a new method to preserve fine details in three dimensional limited angle tomography image reconstruction by using an anisotropic regularization scheme. The results demonstrate clear advantages over the isotropic method and can compete with standard approaches. Moreover, this algorithm handles high noise levels with low photon counts exceptionally well.
L. I. Rudin, S. Osher, and E. Fatemi. Nonlinear total variation based noise removal algorithms. Physica D: Nonlinear Phenomena, 60 (1-4): 259-268, 1992. Google Scholar
Z. Chen, X. Jin, L. Li, and G. Wang. A limited-angle ct reconstruction method based on anisotropic tv minimization. Physics in Medicine and Biology, 58 (7): 2119, 2013. Google Scholar
M. Persson, D. Bone, and H. Elmqvist. Total variation norm for three-dimensional iterative reconstruction in limited view angle tomography. Physics in Medicine and Biology, 46 (3): 853, 2001Google Scholar
Z. Wang, A. C. Bovik, H. R. Sheikh, and E. P. Simoncelli. Image quality assessment: From error visibility to structural similarity. Image Processing, IEEE Transactions on, 13 (4): 600-612, 2004. Google Scholar
About the article
Published Online: 2015-09-12
Published in Print: 2015-09-01
Conflict of interest: Authors state no conflict of interest. Material and Methods: Informed consent: Informed consent has been obtained from all individuals included in this study. Ethical approval: The research related to human use has been complied with all the relevant national regulations, institutional policies and in accordance the tenets of the Helsinki Declaration, and has been approved by the authors’ institutional review board or equivalent committee.