A posteriori error estimation in biomedical imaging
Paper i proceeding, 2007

We present an adaptive hybrid FEM/FDM method for an inverse scattering problem in scanning acoustic microscopy with a special focus on new application in medical imaging. The problem takes the form of reconstructing an unknown sound velocity c(x) from boundary displacement data measured in acoustic microscopy in order to obtain the pathological defects in bone. The inverse problem is formulated as an optimal control problem, where we solve the equations of optimality expressing stationarity of an associated Lagrangian by a quasi-Newton method. We present a posteriori error estimate for the error in the Lagrangian which couples residuals of the computed solution to weights of the reconstruction obtained by solving and associated linearized problem for the Hessian of the Lagrangian. The performance of the adaptive hybrid method and usefulness of the a posteriori error estimator are illustrated in numerical examples. © 2007 IEEE.

Biomedical imaging

Adaptive finite element methods

Inverse scattering

A posteriori error estimation

Transient wave equation

Hybrid finite element/difference method


4th IEEE International Symposium on Biomedical Imaging: From Nano to Macro - Proceedings

1945-7928 (ISSN)