This study focuses on the use of coarse spatial resolution (CR, pixel size about 1 km2) remote sensing data for land cover change detection and qualification. Assuming the linear mixing model for CR pixels, the problem is that both the multitemporal class features and the pixel composition in terms of classes are unknown. The proposed algorithm is then based on the iterative alternate estimation of each unknown variable. At each iteration, the class features are estimated, thanks to the knowledge of the composition of some pixels, and then the pixel composition is re-estimated knowing the class features. The subset of known composition pixels is the subset of pixels where no change has occurred, i.e. the previous land cover map is still valid. It is derived automatically by removing at each iteration the pixels where the new composition estimation disagrees with the former one. Finally, for the final estimation of the pixel composition, a Markovian chain model is used to guide the solution, i.e. the previous land cover map is used as a ‘reminder’ or ‘memory’ term.This approach has been first validated using simulated data with different spatial resolution ratios. Then, the detection of forest change with SPOT/VGT-S10 has been considered as an actual application case. Finally, the method has been applied to change detection on the Val de Saône watershed between the 1980s and 2000. The results obtained from three coarse resolution series, NOAA/AVHRR, SPOT/VGT-S10 and SPOT/VGT-P, have been compared.