A bilinear inverse problem with forward operator inaccuracy applied to neonatal atlas-based diffuse optical tomography

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Hakula, Aada, Hirvi, Pauliina, Hyvönen, Nuutti
Format: Preprint
Publié: 2026
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866911530419748864
author Hakula, Aada
Hirvi, Pauliina
Hyvönen, Nuutti
author_facet Hakula, Aada
Hirvi, Pauliina
Hyvönen, Nuutti
contents Linear inverse problems are highly common in practical real-world applications from industry to medical imaging. The forward operator is often built on some approximations of the studied system. Handling inaccuracies in the forward operator in the context of inverse problems is a relatively unstudied problem. In this work, we assume that we have a set of candidate forward operator matrices and suggest principal component analysis (PCA) for modeling their variation from the mean. We combine the original linear problem with the included forward operator inaccuracy into a bilinear tensor inverse problem and present two optimization algorithms and Gibbs sampling for approximately solving the problem. As a real-world test case, we apply the algorithms to account for the inaccuracy that is present in the sensitivity profiles or Jacobian matrices in diffuse optical tomography when an atlas-based model of the head anatomy is used instead of the subject's own anatomical model in neonates over a wide range of gestational ages (29--44 weeks). We report visual and numerical improvements in the spatial localization and contrast-to-noise-ratio in reconstructions of simulated hemodynamic activity.
format Preprint
id arxiv_https___arxiv_org_abs_2603_18980
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle A bilinear inverse problem with forward operator inaccuracy applied to neonatal atlas-based diffuse optical tomography
Hakula, Aada
Hirvi, Pauliina
Hyvönen, Nuutti
Numerical Analysis
65F10, 78A46, 62H25, 92C55
Linear inverse problems are highly common in practical real-world applications from industry to medical imaging. The forward operator is often built on some approximations of the studied system. Handling inaccuracies in the forward operator in the context of inverse problems is a relatively unstudied problem. In this work, we assume that we have a set of candidate forward operator matrices and suggest principal component analysis (PCA) for modeling their variation from the mean. We combine the original linear problem with the included forward operator inaccuracy into a bilinear tensor inverse problem and present two optimization algorithms and Gibbs sampling for approximately solving the problem. As a real-world test case, we apply the algorithms to account for the inaccuracy that is present in the sensitivity profiles or Jacobian matrices in diffuse optical tomography when an atlas-based model of the head anatomy is used instead of the subject's own anatomical model in neonates over a wide range of gestational ages (29--44 weeks). We report visual and numerical improvements in the spatial localization and contrast-to-noise-ratio in reconstructions of simulated hemodynamic activity.
title A bilinear inverse problem with forward operator inaccuracy applied to neonatal atlas-based diffuse optical tomography
topic Numerical Analysis
65F10, 78A46, 62H25, 92C55
url https://arxiv.org/abs/2603.18980