Пример 11.13. Главные компоненты.

f1 = imread('C: \Users\Alex\Desktop\Washington1.tif');

f2 = imread('C: \Users\Alex\Desktop\Washington2.tif');

f3 = imread('C: \Users\Alex\Desktop\Washington3.tif');

f4 = imread('C: \Users\Alex\Desktop\Washington4.tif');

f5 = imread('C: \Users\Alex\Desktop\Washington5.tif');

f6 = imread('C: \Users\Alex\Desktop\Washington6.tif');

S = cat(3, f1, f2, f3, f4, f5, f6);

[X, R] = imstack2vectors(S);

P = princomp(X, 6);

g1 = P.Y(:, 1);

g1 = reshape(g1, 512, 512);

imshow(g1, [ ])

d = diag(P.Cy);

P = princomp(X, 2);

h1 = P.X(:, 1);

h1 = reshape(h1, 512, 512);

D1 = double(f1), double(h1);

D1 = gscale(D1);

imshow(D1)

P.ems

ans =

1.7311e+003

Пример 12.1. Использование корреляции для сопоставления изображений.

f = imread('C: \Users\Alex\Desktop\Hurricane_Andrew.tif');

w = imread('C: \Users\Alex\Desktop\hurricane_mask.tif');

g = dftcorr(f, w);

gs = gscale(g);

imshow(gs)

[I, J] = find(g == max(g(:)))

I =

J =

imshow(gs > 254)