lucasperso.asv

Lukas Kanade Multi resolution program in matlab language.

% Gaetan Boehringer, Bruno Jacquot & Fr閐閞ic Piegay
%% Algorithme de Lucas & Kanade

%% Chargement des images

clear all;
close all;
clc;

% Lecture des images
im1o=double(imread('LKtest3im1.bmp'));
im2o=double(imread('LKtest3im2.bmp'));
[Lim,Cim]=size(im1o);

% Affichage des images
figure;
subplot(1,2,1); imagesc(im1o);axis equal; colormap(gray);title('Image initiale');
subplot(1,2,2); imagesc(im2o);axis equal; colormap(gray); title ('Image d閏al閑');

%% Filtrage

h = fspecial('gaussian',10,3);
figure; imshow(h,[]);title('gradient')
im1 = imfilter(im1o,h);
im2 = imfilter(im2o,h);

% Affichage des images filtr閑s
figure;
subplot(1,2,1); imagesc(im1);axis equal; colormap(gray);title('Image initiale filtr閑');
subplot(1,2,2); imagesc(im2);axis equal; colormap(gray);title('Image d閏al閑 filtr閑');

%% Calcul en 1D
% Affichage d'un profil

% Extraction de la ligne centrale de chacune des images
y1=im1(:,Cim/2)';
y2=im2(:,Cim/2)';
[L,C]=size(y1);

% Affichage
figure;
plot(y1,'b');
hold on;
plot(y2,'m'); title('Affichage du profil de la ligne des 2 images');

grady=zeros(1,C);
ypf=conv(y1,(1/12)*[-1 8 0 -8 1]); % On applique un filtre