EM Based Approximation of Empirical Distributions with Linear Combinations of Discrete Gaussians

We propose novel expectation maximization (EM) based algorithms for accurate approximation of an empirical probability distribution of discrete scalar data. The algorithms refine our previous ones in that they approximate the empirical distribution with a linear combination of discrete Gaussians (LCDG). The use of the DGs results in closer approximation and considerably better convergence to a local likelihood maximum compared to previously involved conventional continuous Gaussian densities. Experiments in segmenting multimodal medical images show the proposed algorithms produce more adequate region borders.