Tomographic imaging of dynamic objects with the ensemble Kalman filter.
We address the image formation of a dynamic object from projections by formulating it as a state estimation problem. The problem is solved with the ensemble Kalman filter (EnKF), a Monte Carlo algorithm that is computationally tractable when the state dimension is large. In this paper, we first rigorously address the convergence of the EnKF. Then, the effectiveness of the EnKF is demonstrated in a numerical experiment where a highly variable object is reconstructed from its projections, an imaging modality not yet explored with the EnKF. The results show that the EnKF can yield estimates of almost equal quality as the optimal Kalman filter but at a fraction of the computational effort. Further experiments explore the rate of convergence of the EnKF, its performance relative to an idealized particle filter, and implications of modeling the system dynamics as a random walk.
Butala
Mark D
MD
University of Illinois at Urbana-Champaign, Urbana, IL 61801 USA. butala@uiuc.edu
Frazin
Richard A
RA
Chen
Yuguo
Y
Kamalabadi
Farzad
F
D000465:Algorithms; D001921:Brain / Q000033:anatomy & histology; D006801:Humans; D007091:Image Processing, Computer-Assisted / Q000379:methods*; D008962:Models, Theoretical; D009010:Monte Carlo Method; D013269:Stochastic Processes; D014054:Tomography / Q000379:methods*