A Regression Mixture Model with Spatial Constraints for Clustering Spatiotemporal Data

We present a new approach for curve clustering designed for analysis of spatiotemporal data. Such data contains both spatial and temporal patterns that we desire to capture. The proposed methodology is based on regression and Gaussian mixture modeling. The novelty of the herein work is the incorporation of spatial smoothness constraints in the form of a prior for the data labels. This allows to take into account the property of spatiotemporal data according to which spatially adjacent data points have higher probability to belong to the same cluster. The proposed model can be formulated as a Maximum a Posteriori (MAP) problem, where the Expectation Maximization (EM) algorithm is used to estimate the model parameters. Several numerical experiments with both simulated data and real cardiac perfusion MRI data are used for evaluating the methodology. The results are promising and demonstrate the value of the proposed approach.

[1]  Wasserman,et al.  Bayesian Model Selection and Model Averaging. , 2000, Journal of mathematical psychology.

[2]  N. Tsekos,et al.  Fast magnetization‐driven preparation for imaging of contrast‐enhanced coronary arteries during intra‐arterial injection of contrast agent , 2006, Journal of magnetic resonance imaging : JMRI.

[3]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[4]  Padhraic Smyth,et al.  Probabilistic curve-aligned clustering and prediction with regression mixture models , 2004 .

[5]  W. DeSarbo,et al.  A maximum likelihood methodology for clusterwise linear regression , 1988 .

[6]  Sally Wood,et al.  Bayesian mixture of splines for spatially adaptive nonparametric regression , 2002 .

[7]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[8]  Nikolas P. Galatsanos,et al.  A spatially constrained mixture model for image segmentation , 2005, IEEE Transactions on Neural Networks.

[9]  Adrian E. Raftery,et al.  Bayesian Regularization for Normal Mixture Estimation and Model-Based Clustering , 2007, J. Classif..

[10]  Donald B. Rubin,et al.  Max-imum Likelihood from Incomplete Data , 1972 .

[11]  Nikolaos V Tsekos,et al.  Dynamic coronary MR angiography and first‐pass perfusion with intracoronary administration of contrast agent , 2002, Journal of magnetic resonance imaging : JMRI.

[12]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[13]  Nikolas P. Galatsanos,et al.  A Class-Adaptive Spatially Variant Mixture Model for Image Segmentation , 2007, IEEE Transactions on Image Processing.

[14]  P. Green Bayesian reconstructions from emission tomography data using a modified EM algorithm. , 1990, IEEE transactions on medical imaging.

[15]  Nikolas P. Galatsanos,et al.  Methods for choosing the regularization parameter and estimating the noise variance in image restoration and their relation , 1992, IEEE Trans. Image Process..

[16]  Nikolas P. Galatsanos,et al.  Segmentation of dynamic PET or fMRI images based on a similarity metric , 2003 .

[17]  Debiao Li,et al.  Magnetic resonance imaging–guided coronary interventions , 2004, Journal of magnetic resonance imaging : JMRI.

[18]  Padhraic Smyth,et al.  Curve Clustering with Random Effects Regression Mixtures , 2003, AISTATS.

[19]  Geoffrey J. McLachlan,et al.  Finite Mixture Models , 2019, Annual Review of Statistics and Its Application.

[20]  Radford M. Neal Pattern Recognition and Machine Learning , 2007, Technometrics.

[21]  Adrian E. Raftery,et al.  How Many Clusters? Which Clustering Method? Answers Via Model-Based Cluster Analysis , 1998, Comput. J..

[22]  Jean-Philippe Thiran,et al.  Counting Pedestrians in Video Sequences Using Trajectory Clustering , 2006, IEEE Transactions on Circuits and Systems for Video Technology.

[23]  Stephen M. Smith,et al.  Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm , 2001, IEEE Transactions on Medical Imaging.