Source author record

Philip Richardson

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

2works
1topics
1close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

2 published item(s)

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

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.