Total Generalized Variation for Inverse Problems Involving Piecewise Linear Finite Elements

Abstract

Higher-order regularization has shown improved results over first-order methods such as total variation. In this paper, convergence of a piecewise linear finite element discretization is shown for the second-order total generalized variation functional. Additionally, this convergence result is extended to a spatiotemporal setting after the spatiotemporal function space is comprehensively introduced. In both cases, discrete minimizers converge to the continuous minimizer with rate $h^{1/4}$, matching the best known rate for finite element discretizations of total variation in the same general setting. Numerical experiments confirm the convergence rates and indicate faster convergence in practice. Comparisons with related spatial image reconstruction algorithms show comparable reconstruction results. The spatiotemporal case is discussed using the example of the inverse problem in electrocardiographic imaging, where the proposed functional improves the state of the art.