Abstract:Single Particle Tracking (SPT) is a well known class of tools for studying the dynamics of biological macromolecules moving inside living cells. In this paper, we focus on the problem of localization and parameter estimation given a sequence of segmented images. In the standard paradigm, the location of the emitter inside each frame of a sequence of camera images is estimated using, for example, Gaussian fitting (GF), and these locations are linked to provide an estimate of the trajectory. Trajectories are then analyzed by using Mean Square Displacement (MSD) or Maximum Likelihood Estimation (MLE) techniques to determine motion parameters such as diffusion coefficients. However, the problems of localization and parameter estimation are clearly coupled. Motivated by this, we have created an Expectation Maximization (EM) based framework for simultaneous localization and parameter estimation. We demonstrate this framework through two representative methods, namely, Sequential Monte Carlo combined with Expectation Maximization (SMC-EM) and Unscented Kalman Filter combined with Expectation Maximization (U-EM). Using diffusion in two-dimensions as a prototypical example, we conduct quantitative investigations on localization and parameter estimation performance across a wide range of signal to background ratios and diffusion coefficients and compare our methods to the standard techniques based on GF-MSD/MLE. To demonstrate the flexibility of the EM based framework, we do comparisons using two different camera models, an ideal camera with Poisson distributed shot noise but no readout noise, and a camera with both shot noise and the pixel-dependent readout noise that is common to scientific complementary metal-oxide semiconductor (sCMOS) camera. Our results indicate our EM based methods outperform the standard techniques, especially at low signal levels. While U-EM and SMC-EM have similar accuracy, U-EM is significantly more computationally efficient, though the use of the Unscented Kalman Filter limits U-EM to lower diffusion rates.

EM-based algorithms for single particle tracking of Ornstein-Uhlenbeck motion from sCMOS camera data

Simultaneous Localization and Parameter Estimation for Single Particle Tracking via Sigma Points based EM

Expectation maximization based framework for joint localization and parameter estimation in single particle tracking from segmented images

A time-varying approach to single particle tracking with a nonlinear observation model

Joint Estimation of Trajectory and Model Parameters for Single Particle Tracking of 3D Confined Diffusion Using the Double-Helix Point Spread Function

An EM Algorithm for Nonlinear State Estimation With Model Uncertainties

Multiple Dense Particle Tracking In Fluorescence Microscopy Images Based On Multidimensional Assignment

Embedding Ensemble Tracking in a Stochastic Framework for Robust Object Tracking

Method for simultaneous localization and parameter estimation in particle tracking experiments

Visual Tracking Via Subspace Motion Model

Estimating mixed-effects state-space models via particle filters and the EM algorithm

An expectation-maximization algorithm for positron emission particle tracking

Tracking A Large Number of Closely Spaced Objects Based on the Particle Probability Hypothesis Density Filter Via Optical Sensor

Subcellular particles tracking in time-lapse confocal microscopy images

An efficient particle-based method for maximum likelihood estimation in nonlinear state-space models

A Novel EM Implementation for Initial Alignment of SINS Based on Particle Filter and Particle Swarm Optimization

UmUTracker: A versatile MATLAB program for automated particle tracking of 2D light microscopy or 3D digital holography data

Tracking of Extended Object Using Random Matrix With Non-Uniformly Distributed Measurements

3D Human Motion Tracking Combining Part Segmentation and Particle Filter