Tomographic image reconstruction can be mapped to a problem of finding solutions to a large system of linear equations which maximize a function that includes extita priori knowledge regarding features of typical images such as smoothness or sharpness. This maximization can be performed with standard local optimization tools when the function is concave, but it is generally intractable for realistic priors, which are non-concave. We introduce a new method to reconstruct images obtained from Radon projections by using Expectation Propagation, which allows us to reframe the problem from an Bayesian inference perspective. We show, by means of extensive simulations, that, compared to state-of-the-art algorithms for this task, Expectation Propagation paired with very simple but non log-concave priors, is often able to reconstruct images up to a smaller error while using a lower amount of information per pixel. We provide estimates for the critical rate of information per pixel above which recovery is error-free by means of simulations on ensembles of phantom and real images.
Non-convex image reconstruction via expectation propagation / Paola Muntoni, Anna; DIAZ HERNANDEZ ROJAS, RAFAEL ALBERTO; Braunstein, Alfredo; PAGNANI, Andrea; PEREZ CASTILLO, Isaac. - In: PHYSICAL REVIEW. E. - ISSN 2470-0053. - 100:23(2018). [10.1103/PhysRevE.100.032134]
Non-convex image reconstruction via expectation propagation
Rafael Díaz Hernández RojasSecondo
;Andrea Pagnani;Isaac Pérez Castillo
2018
Abstract
Tomographic image reconstruction can be mapped to a problem of finding solutions to a large system of linear equations which maximize a function that includes extita priori knowledge regarding features of typical images such as smoothness or sharpness. This maximization can be performed with standard local optimization tools when the function is concave, but it is generally intractable for realistic priors, which are non-concave. We introduce a new method to reconstruct images obtained from Radon projections by using Expectation Propagation, which allows us to reframe the problem from an Bayesian inference perspective. We show, by means of extensive simulations, that, compared to state-of-the-art algorithms for this task, Expectation Propagation paired with very simple but non log-concave priors, is often able to reconstruct images up to a smaller error while using a lower amount of information per pixel. We provide estimates for the critical rate of information per pixel above which recovery is error-free by means of simulations on ensembles of phantom and real images.File | Dimensione | Formato | |
---|---|---|---|
Muntoni_Nonconvex-image-reconstruction_2019.pdf
accesso aperto
Tipologia:
Versione editoriale (versione pubblicata con il layout dell'editore)
Licenza:
Tutti i diritti riservati (All rights reserved)
Dimensione
1.73 MB
Formato
Adobe PDF
|
1.73 MB | Adobe PDF |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.