In this paper, we study the problem of detecting and estimating change-points in a time series of multivariate images. We extend existent works to take into account the heterogeneity of the dataset on a spatial neighbourhood. The classic complex Gaussian assumption of the data is replaced by a complex elliptically symmetric assumption. Then a robust statistics are derived using Generalised Likelihood Ratio Test (GLRT). These statistics are coupled to an estimation strategy for one or several changes. Performance of these robust statistics have been analyzed in simulation and compared to the one associated with standard multivariate normal assumption. When the data is heterogeneous, the detection and estimation strategy yields better results with the new statistics.