Haase V, Stierstorfer K, Hahn K, Schöndube H, Maier A, Noo F (2019)
Publication Type: Conference contribution
Publication year: 2019
Publisher: SPIE
Book Volume: 10948
Conference Proceedings Title: Progress in Biomedical Optics and Imaging - Proceedings of SPIE
ISBN: 9781510625433
DOI: 10.1117/12.2513375
Because total variation (TV) is non-differentiable, iterative reconstruction using a TV penalty comes with technical difficulties. To avoid these, it is popular to use a smooth approximation of TV instead, which is defined using a single parameter, herein called δ, with the convention that the approximation is better when δ is smaller. To our knowledge, it is not known how well image reconstruction with this approximation can approach a converged non-smooth TV-regularized result. In this work, we study this particular question in the context of X-ray computed tomography (CT). Experimental results are reported with real CT data of a head phantom and supplemented with a theoretical analysis. To address our question, we make proficient use of a generalized iterative soft-thresholding algorithm that allows us to handle TV and its smooth approximation in the same framework. Our results support the following conclusions. First, images reconstructed using the smooth approximation of TV appears to smoothly converge towards the TV result as δ tends to zero. Second, the value of δdoes not need to be overly small to obtain a result that is essentially equivalent to TV, implying that numerical instabilities can be avoided. Last, though it is smooth, the convergence with δ is not particularly fast, as the mean absolute pixel difference decreases only as √δ in our experiments. Altogether, we conclude that the approximation is a theoretically valid way to approximate the non-smooth TV penalty for CT, opening the door to safe utilization of a wide variety of optimization algorithms.
APA:
Haase, V., Stierstorfer, K., Hahn, K., Schöndube, H., Maier, A., & Noo, F. (2019). Exploring the space between smoothed and non-smooth total variation for 3D iterative CT reconstruction. In Hilde Bosmans, Guang-Hong Chen, Taly Gilat Schmidt (Eds.), Progress in Biomedical Optics and Imaging - Proceedings of SPIE. San Diego, CA, US: SPIE.
MLA:
Haase, Viktor, et al. "Exploring the space between smoothed and non-smooth total variation for 3D iterative CT reconstruction." Proceedings of the Medical Imaging 2019: Physics of Medical Imaging, San Diego, CA Ed. Hilde Bosmans, Guang-Hong Chen, Taly Gilat Schmidt, SPIE, 2019.
BibTeX: Download