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 neighborhood. The classic complex Gaussian assumption of the data is replaced by a complex elliptically symmetric assumption. Then robust statistics are derived using Generalized 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.