New Pair of Primal Dual Algorithms for Bregman Iterated Variational Regularization
For researchers in inverse problems and regularization, this offers a new algorithmic framework with theoretical guarantees, though it is incremental as it extends existing primal-dual and Bregman iteration concepts.
This work introduces a pair of primal-dual algorithms for Bregman iterated variational regularization, providing total error estimation and convergence rates under variational source conditions. The algorithms are validated on GPS tomography and 2-D image reconstruction, demonstrating iterative regularization behavior.
Primal-dual splitting involving proximity operators in order to be able to find some approximation to the minimizer for a general form of Tikhonov type functional is in the focus of this work. This approximation is produced by a pair of iterative variational regularization procedures. Under the assumption of some variational source condition (VSC), total error estimation both in the iterative sense and in the continuous sense has been analysed separately. Rates of convergence will be obtained in terms some concave and positive definite index function. Of the choice of the penalty term, we are interested in Bregman distance penalization associated with the non-smooth total variation (TV) functional. Furthermore, following up the lower and bounds defined for the regularization parameter, some deterministic choice of the regularization parameter is given explicitly. It is in the emphasis of this work that the regularization parameter obeys Morozov`s discrepancy principle (MDP) in order for the stability analysis of regularized solution. In the computerized environment, the algorithms are verified as iterative regularization methods by applying it to an atmospheric tomography problem named as GPS-Tomography. Apart from this 3-D tomographic inverse problem, we also apply the algorithms to some 2-D conventional tomographic image reconstruction problems in order to be able test algorithms` capability of capturing the details and observe that algorithms behave as iterative regularization procedures.