Discussion :

[hajlaouirami]

Bonjour à tous

j'ai une problème avec la lecture d'une image sous matlab, qui je ne trouve pas une solution.

le problème est si la longueur est supérieur au largeur mon code fonctionne correctement mais si la longueur d'image est inférieur au largeur mon code donne cette erreur

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
??? Index exceeds matrix dimensions.
 
Error in ==> capte at 3
    X=b(:,1:i,:);
 
Error in ==> test at 26
    [X Y s]=capte(c,a,i);

[Jerome Briot]

Montre nous le code de la fonction capte et explique nous ce qu'elle fait

[hajlaouirami]

voila le code du fonction capte:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
function [X Y s]=capte(b,a,i)
[nl nc z]=size(a);
    X=b(:,1:i,:);
    Y=a(:,end-i+1:end,:);
    s=sum(sum(sum(abs(double(X)-double(Y)))))/(i*nl);
 
end

cette fonction permet de calculer la difference en valeur absolue des pixels de 2 images qui par la suite je l'utilise pour extraire les parties communes entre les images

[Jerome Briot]

Que représente le paramètre i ?

Montre nous le code de test.m

Au passage, la fonction se simplifie en :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
function [X Y s]=capte(b,a,i)
 
    nl = size(a,1);
    X = b(:,1:i,:);
    Y = a(:,end-i+1:end,:);
    s = sum(abs(double(X(:))-double(Y(:))))/(i*nl);

voire, si a et b sont de même classe :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
function [X Y s]=capte(b,a,i)
 
    nl = size(a,1);
    X = b(:,1:i,:);
    Y = a(:,end-i+1:end,:);
    s = sum(abs(double(X(:)-Y(:))))/(i*nl);

[hajlaouirami]

Voici le code du fichier test.m

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
clear all;close all;clc
a=imread('x1.ppm');
c=imread('x2.ppm');
 
imshow(a);figure;
 
[nl nc z]=size(a);
for i=1:nl-1
    [X Y s]=capte(c,a,i);
    S(i)=s;
 
end
figure;plot(S)
R=find (S==min(S))
 
 
[X Y s]=capte(c,a,R);
figure;imshow(X);figure:imshow(Y);

[Jerome Briot]

Ne serait-ce pas plutôt ceci ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
for i=1:nc-1

   ...

end