Source author record

Sean Holman

Sean Holman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.AP math.NA eess.IV math.DG Numerical Analysis

Catalog footprint

What is connected

5works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Simultaneous recovery of attenuation and source density in SPECT

We show that under a certain non-cancellation condition the attenuated Radon transform uniquely determines piecewise constant attenuation $a$ and piecewise $C^2$ source density $f$ with jumps over real analytic boundaries possibly having corners. We also look at numerical examples in which the non-cancellation condition fails and show that unique reconstruction of multi-bang $a$ and $f$ is still appears to be possible although not yet explained by theoretical results.

preprint2020arXiv

Data-consistent neural networks for solving nonlinear inverse problems

Data assisted reconstruction algorithms, incorporating trained neural networks, are a novel paradigm for solving inverse problems. One approach is to first apply a classical reconstruction method and then apply a neural network to improve its solution. Empirical evidence shows that such two-step methods provide high-quality reconstructions, but they lack a convergence analysis. In this paper we formalize the use of such two-step approaches with classical regularization theory. We propose data-consistent neural networks that we combine with classical regularization methods. This yields a data-driven regularization method for which we provide a full convergence analysis with respect to noise. Numerical simulations show that compared to standard two-step deep learning methods, our approach provides better stability with respect to structural changes in the test set, while performing similarly on test data similar to the training set. Our method provides a stable solution of inverse problems that exploits both the known nonlinear forward model as well as the desired solution manifold from data.

preprint2020arXiv

Emission tomography with a multi-bang assumption on attenuation

We consider the problem of joint reconstruction of both attenuation $a$ and source density $f$ in emission tomography in two dimensions. This is sometimes called the Single Photon Emission Computed Tomography (SPECT) identification problem, or referred to as attenuation correction in SPECT. Assuming that $a$ takes only finitely many values and $f \in C_c^1(\mathbb{R}^2)$ we are able to characterise singularities appearing in the Attenuated Radon Transform $R_a f$, which models emission tomography data. Using this characterisation we prove that both $a$ and $f$ can be determined in some circumstances. We also propose a numerical algorithm to jointly compute $a$ and $f$ from $R_af$ based on a weakly convex regularizer when $a$ only takes values from a known finite list, and show that this algorithm performs well on some synthetic examples.

preprint2016arXiv

A Multi-Grid Iterative Method for Photoacoustic Tomography

Inspired by the recent advances on minimizing nonsmooth or bound-constrained convex functions on models using varying degrees of fidelity, we propose a line search multigrid (MG) method for full-wave iterative image reconstruction in photoacoustic tomography (PAT) in heterogeneous media. To compute the search direction at each iteration, we decide between the gradient at the target level, or alternatively an approximate error correction at a coarser level, relying on some predefined criteria. To incorporate absorption and dispersion, we derive the analytical adjoint directly from the first-order acoustic wave system. The effectiveness of the proposed method is tested on a total-variation penalized Iterative Shrinkage Thresholding algorithm (ISTA) and its accelerated variant (FISTA), which have been used in many studies of image reconstruction in PAT. The results show the great potential of the proposed method in improving speed of iterative image reconstruction.

preprint2015arXiv

On the microlocal analysis of the geodesic X-ray transform with conjugate points

We study the microlocal properties of the geodesic X-ray transform $\mathcal{X}$ on a manifold with boundary allowing the presence of conjugate points. Assuming that there are no self-intersecting geodesics and all conjugate pairs are nonsingular we show that the normal operator $\mathcal{N} = \mathcal{X}^t \circ \mathcal{X}$ can be decomposed as the sum of a pseudodifferential operator of order $-1$ and a sum of Fourier integral operators. We also apply this decomposition to prove inversion of $\mathcal{X}$ is only mildly ill-posed in dimension three or higher.