Change detection by the IR-MAD and kernel MAF methods in Landsat TM data covering a Swedish forest region

• Main Authors: Nielsen, Allan Aasbjerg, Olsson, Håkan
• Other Authors: Miranda, David, Suárez, Juan, Crecente, Rafael
• Format: Conference Object
• Language: English
• Published: 2010
• Subjects: Northern Sweden
• Online Access: https://orbit.dtu.dk/en/publications/389e5947-ede0-449a-878e-ba9b277dfc6c; http://www2.imm.dtu.dk/pubdb/views/edoc_download.php/5930/pdf/imm5930.pdf

Change over time between two 512 by 512 (25 m by 25 m pixels) multispectral Landsat Thematic Mapper images dated 6 June 1986 and 27 June 1988 respectively covering a forested region in northern Sweden, is here detected by means of the iteratively reweighted multivariate alteration detection (IR-MAD) method followed by post-processing by means of kernel maximum autocorrelation factor (kMAF) analysis. The IR-MAD method builds on an iterated version of an established method in multivariate statistics, namely canonical correlation analysis (CCA). It finds orthogonal (i.e., uncorrelated) linear combinations of the multivariate data at two time points that have maximal correlation. These linear combinations are called the canonical variates (CV) and the corresponding correlations are called the canonical correlations. There is one set of CVs for each time point. The difference between the two set of CVs represent the change between the two time points and are called the MAD variates or the MADs for short. The MAD variates are invariant to linear and affine transformations of the original data. The sum of the squared MAD variates (properly normed to unit variance) gives us change variables that will ideally follow a so-called c2 (chi-squared) distribution with p degrees of freedom for the no-change pixels; p is the number of spectral bands in the image data. Here p=6, the thermal band is excluded from the analyses. The c2 image is the basis for calculating an image of probability for no-change, i.e., the probability for finding a higher value of the c2 statistic than the one actually found. This image is the weight image in the iteration scheme mentioned above. Iterations stop when the canonical correlations stop changing. Principal component analysis (PCA) finds orthogonal (i.e., uncorrelated) linear combinations of the multivariate data that have maximal variance. A kernel version of PCA is based on a dual formulation also termed Q-mode analysis in which the data enter into the analysis via inner products in the ...