Dynamic Cell Imaging in PET with Optimal Transport Regularization
arXiv:1902.07521 · doi:10.1109/TMI.2019.2953773
Abstract
We propose a novel dynamic image reconstruction method from PET listmode data that could be particularly suited to tracking single or small numbers of cells. In contrast to conventional PET reconstruction our method combines the information from all detected events not only to reconstruct the dynamic evolution of the radionuclide distribution, but also to improve the reconstruction at each single time point by enforcing temporal consistency. This is achieved via optimal transport regularization where in principle, among all possible temporally evolving radionuclide distributions consistent with the PET measurement, the one is chosen with least kinetic motion energy. The reconstruction is found by convex optimization so that there is no dependence on the initialization of the method. We study its behaviour on simulated data of a human PET system and demonstrate its robustness even in settings with very low radioactivity. In contrast to previously reported cell tracking algorithms, our technique is oblivious to the number of tracked cells. Without any additional complexity one or multiple cells can be reconstructed, and the model automatically determines the number of particles. For instance, four radiolabelled cells moving at a velocity of 3.1 mm/s and a PET recorded count rate of 1.1 cps (for each cell) could be simultaneously tracked with a tracking accuracy of 5.3 mm inside a simulated human body.
arXiv commentRevised version, to appear in IEEE Trans Med Imaging. Supplementary material attached as last page
References in corpus (2)
Cited by in corpus (9)
- A generalized conditional gradient method for dynamic inverse problems with optimal transport regularization
- On the extremal points of the ball of the Benamou-Brenier energy
- An optimal transport approach for solving dynamic inverse problems in spaces of measures
- Asymptotic linear convergence of fully-corrective generalized conditional gradient methods
- Towards a mathematical theory of trajectory inference
- Regularization with optimal space-time priors
- Complex extension of optical flow and its practical evaluation for undersampled dynamic MRI
- Motion simulation of radio-labeled cells in whole-body positron emission tomography
- Sparsity for dynamic inverse problems on Wasserstein curves with bounded variation