High-Level Algorithm Prototyping: An Example Extending the TVR-DART Algorithm

• Published in: Discrete Geometry for Computer Imagery Vol. 10502; pp. 109 - 121
• Main Authors: Ringh, Axel, Zhuge, Xiaodong, Palenstijn, Willem Jan, Batenburg, Kees Joost, Öktem, Ozan
• Format: Book Chapter
• Language: English
• Published: Switzerland Springer International Publishing AG 2017; Springer International Publishing
• Series: Lecture Notes in Computer Science
• ISBN: 3319662716; 9783319662718
• ISSN: 0302-9743; 1611-3349
• DOI: 10.1007/978-3-319-66272-5_10

Operator Discretization Library (ODL) is an open-source Python library for prototyping reconstruction methods for inverse problems, and ASTRA is a high-performance Matlab/Python toolbox for large-scale tomographic reconstruction. The paper demonstrates the feasibility of combining ODL with ASTRA to prototype complex reconstruction methods for discrete tomography. As a case in point, we consider the total-variation regularized discrete algebraic reconstruction technique (TVR-DART). TVR-DART assumes that the object to be imaged consists of a limited number of distinct materials. The ODL/ASTRA implementation of this algorithm makes use of standardized building blocks, that can be combined in a plug-and-play manner. Thus, this implementation of TVR-DART can easily be adapted to account for application specific aspects, such as various noise statistics that come with different imaging modalities.