The statistical techniques of multivariate alteration detection, maximum autocorrelation factor transformation, expectation maximization, fuzzy maximum likelihood estimation and probabilistic label relaxation are combined in a unified scheme to classify changes in multispectral satellite data. An example involving bitemporal LANDSAT TM imagery is given.