Abstract:The single-scatter approximation is fundamental in many tomographic imaging problems including x-ray scatter imaging and optical scatter imaging for certain media. In all cases, noisy measurements are affected by both local scatter events and nonlocal attenuation. Prior works focus on reconstructing one of two images: scatter density or total attenuation. However, both images are media specific and useful for object identification. Nonlocal effects of the attenuation image on the data are summarized by the broken ray transform (BRT). While analytic inversion formulas exist, poor conditioning of the inverse problem is only exacerbated by noisy measurements and sampling errors. This has motivated interest in the related star transforms incorporating BRT measurements from multiple source-detector pairs. However, all analytic methods operate on the log of the data. For media comprising regions with no scatter a new approach is required. We are the first to present a joint estimation algorithm based on Poisson data models for a single-scatter measurement geometry. Monotonic reduction of the log-likelihood function is guaranteed for our iterative algorithm while alternating image updates. We also present a fast algorithm for computing the discrete BRT forward operator. Our generalized approach can incorporate both transmission and scatter measurements from multiple source-detector pairs. Transmission measurements resolve low-frequency ambiguity in the joint image estimation problem, while multiple scatter measurements resolve the attenuation image. The benefits of joint estimation, over single-image estimation, vary with problem scaling. Our results quantify these benefits and should inform design of future acquisition systems.

Inversion formulas for the broken-ray Radon transform

Iterative Algorithms for Joint Scatter and Attenuation Estimation From Broken Ray Transform Data

Inversion formulae for ray transforms in vector and tensor tomography

Exact reconstruction formulas for a Radon transform over cones

Revisit to the Inverse Exponential Radon Transform

Inversion of the star transform

Spline-based Inverse Radon Transform in Two and Three Dimensions

On an Analytical Inversion Formula for the Modulo Radon Transform

Generalized Abel equations and applications to translation invariant Radon transforms

Reconstructions from Parallel- and Fan-Beam Data with Truncations

Inversion of the Radon Transform Associated with the Classical Domain of Type One

Weighted Radon transforms of vector fields, with applications to magnetoacoustoelectric tomography

An Alternative Derivation for Bronnikov’s Reconstruction Formula in X-Ray Phase Contrast Tomography

Inversion of generalized V-line transforms of vector fields in $\mathbb{R}^2$

A simple shearlet-based reconstruction for computer tomography

An improved exact inversion formula for cone beam vector tomography

Plancherel formula for the attenuated Radon transform

Inversion of a restricted transverse ray transform with sources on a curve

Image Reconstruction for Source Trajectories with Kinks

Survey on Broken Ray Transforms

Inversion of the pair of weighted and classical circular Radon transforms in C ( R 2 )