Inversion of large-support ill-conditioned linear operators using a Markov model with a line process
We propose a method for the reconstruction of an image, only partially observed through a linear integral operator. As such an inverse problem is ill-posed, prior information must be introduced. We consider the case of a compound Markov random field with a non-interacting line process. In order to maximise the posterior likelihood function, we propose an extension of the graduated non convexity principle pioneered by Blake and Zisserman (1987) which allows its use for ill-posed linear inverse problems. We discuss the role of the observation scale and some aspects of the implemented algorithm. Finally, we present an application of the method to a diffraction tomography imaging problem.<<ETX>>