A probabilistic framework for uncertainty quantification in positron emission particle tracking. (1st May 2023)

• Record Type: Journal Article
• Authors: Offner, Avshalom; Manger, Sam; Vanneste, Jacques
• Is Part Of: Inverse problems. Volume 39:Number 5(2023)
• Publication Date: 2023-05-01
• Subjects: positron emission particle tracking -- Bayesian inference -- uncertainty quantification; Inverse problems (Differential equations) -- Periodicals; 515.357
• DOI: 10.1088/1361-6420/acc47d
• Languages: English
• ISSNs: 0266-5611

Abstract

Abstract: Positron emission particle tracking (PEPT) is an imaging method for the visualization of fluid motion, capable of reconstructing three-dimensional trajectories of small tracer particles suspended in nearly any medium, including fluids that are opaque or contained within opaque vessels. The particles are labeled radioactively, and their positions are reconstructed from the detection of pairs of back-to-back photons emitted by positron annihilation. Current reconstruction algorithms are heuristic and typically based on minimizing the distance between the particles and the so-called lines of response (LoRs) joining the detection points, while accounting for spurious LoRs generated by scattering. Here we develop a probabilistic framework for the Bayesian inference and uncertainty quantification of particle positions from PEPT data. We formulate a likelihood by describing the emission of photons and their noisy detection as a Poisson process in the space of LoRs. We derive formulas for the corresponding Poisson rate in the case of cylindrical detectors, accounting for both undetected and scattered photons. We illustrate the formulation by quantifying the uncertainty in the reconstruction of the position of a single particle on a circular path from data generated by state-of-the-art Monte Carlo simulations. The results show how the observation time Δ t can be chosen optimally to balance the need for a large number of LoRs with the requirement of small particle displacement imposed by the assumption that the particle is static over Δ t . We further show how this assumption can be relaxed by inferring jointly the position and velocity of the particle, with clear benefits for the accuracy of the reconstruction.