Graph- and finite element-based total variation models for the inverse problem in diffuse optical tomography

Wenqi Lu, Jinming Duan, David Orive Miguel, Lionel Herve, Iain Styles

Research output: Contribution to journalArticlepeer-review

5 Citations (Scopus)
204 Downloads (Pure)


Total variation (TV) is a powerful regularization method that has been widely applied in different imaging applications, but is difficult to apply to diffuse optical tomography (DOT) image reconstruction (inverse problem) due to unstructured discretization of complex geometries, non-linearity of the data fitting and regularization terms, and non-differentiability of the regularization term. We develop several approaches to overcome these difficulties by: i) defining discrete differential operators for TV regularization using both finite element and graph representations; ii) developing an optimization algorithm based on the alternating direction method of multipliers (ADMM) for the non-differentiable and non-linear minimization problem; iii) investigating isotropic and anisotropic variants of TV regularization, and comparing their finite element (FEM)- and graph-based implementations. These approaches are evaluated on experiments on simulated data and real data acquired from a tissue phantom. Our results show that both FEM and graph-based TV regularization is able to accurately reconstruct both sparse and non-sparse distributions without the over-smoothing effect of Tikhonov regularization and the over-sparsifying effect of L1 regularization. The graph representation was found to out-perform the FEM method for low-resolution meshes, and the FEM method was found to be more accurate for high-resolution meshes.
Original languageEnglish
Pages (from-to)2684-2707
JournalBiomedical Optics Express
Issue number6
Publication statusPublished - 2 May 2019

ASJC Scopus subject areas

  • Biotechnology
  • Atomic and Molecular Physics, and Optics


Dive into the research topics of 'Graph- and finite element-based total variation models for the inverse problem in diffuse optical tomography'. Together they form a unique fingerprint.

Cite this