58 lines
2.4 KiB
Matlab
58 lines
2.4 KiB
Matlab
% Calcul d'une mosaique d'image à partir de 2 images : I1 et I2
|
|
% connaissant l'homographie H entre les 2.
|
|
% L'image resultat est stockee dans Res.
|
|
% On choisit de projeter I2 dans I1 pour construire la mosaique.
|
|
% Attention !!!
|
|
% On suppose un axe de rotation parallèle aux colonnes.
|
|
% C'est la raison pour laquelle on inverse les lignes et les colonnes
|
|
% dans la reconstruction de la mosaique.
|
|
|
|
function [Imos] = mosaiqueN(TailleFenetre, NbPoints, k, seuil, varargin)
|
|
Imos = im2double(varargin{1});
|
|
for i = 2:nargin - 4
|
|
Im1 = Imos;
|
|
Im2 = im2double(varargin{i});
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
% Detection des points d'interet avec Harris %
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
Im1_gray = rgb2gray(Im1);
|
|
Im2_gray = rgb2gray(Im2);
|
|
XY_1 = harris(Im1_gray, TailleFenetre, NbPoints, k);
|
|
XY_2 = harris(Im2_gray, TailleFenetre, NbPoints, k);
|
|
|
|
% affichage temporaire
|
|
figure;
|
|
affichage_POI(Im1,XY_1,'POI Image 1',1,2,1);
|
|
affichage_POI(Im2,XY_2,'POI Image 2',1,2,2);
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
% Appariement des points d'interet %
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
[XY_C1, XY_C2] = apparier_POI(Im1_gray, XY_1, Im2_gray, XY_2, TailleFenetre, seuil);
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
% Estimation (et verification) de l'homographie %
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
s = 4;
|
|
t = 0.999;
|
|
fitFuncHomo = @(XY) homographie(XY(:, 1:2), XY(:, 3:4));
|
|
disFuncHomo = @(H, XY) distance_homographie(XY(:, 1:2), XY(:, 3:4), H);
|
|
[H, idx] = ransac([XY_C1, XY_C2], fitFuncHomo, disFuncHomo, s, t);
|
|
|
|
figure('Name', 'Appariement des points d''interet RANSAC');
|
|
affichage_appariement(Im1, XY_C1(idx, 1), XY_C1(idx, 2), 'Points d''interet Image 1', 1, 2, 1);
|
|
affichage_appariement(Im2, XY_C2(idx, 1), XY_C2(idx, 2), 'Points d''interet correspondants Image 2', 1, 2, 2);
|
|
|
|
% https://github.com/hero9968/Multiple-View-Geometry-in-Computer-Vision/blob/master/vgg_examples/ransacfithomography_vgg.m
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
% Collage des deux images %
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
Imos = mosaiquecoul(Im1, Im2, H);
|
|
figure;
|
|
affichage_image(Imos,'Mosaique obtenue a partir des 2 images initiales',1,1,1);
|
|
end
|
|
|
|
end
|