Discussion :

[transistor]

bonjour a vous
je voudrait détecté des point d'intérêt dans une image pour cela j ai utilise un balayage de la photo par une matrice de 3*3 suivant les lignes et les colonnes et je calcule une fonction qui définie le point d'intérêt.
le balayage de l'image et fait comme suit

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
I=imread('I.bmp');
Es=2; %le valeur sueil liée a la fonction
[N M]=size(I);
 
for j=1 : N
    for i=1 : M
        for k=0 : 2
            for r=0 : 2
                E=E+sqr(abs((I(i+r,j+k-1))-(I(i+r-1,j+k-1)))) ;
            end ;
        end ;
    end ;
    if E>Es
        text (i,j,'p','color','b') ;
    end ;
end ;

mais j'obtien le msg d'erreur suivant

Code :

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

??? Subscript indices must either be real positive integers or logicals.

vous pouvez m'aider s'il vpus plais.

[Jerome Briot]

Pour j=1 et k=0 => j+k-1=0

Tu essaies donc d'indexer I(1,0) ce qui est impossible sous MATLAB

=> : Subscript indices must either be real positive integers or logicals.

[transistor]

donc je doit commencé mes variable pour la valeur de K et R a partir de 1

[Jean Dumoncel]

Bonjour,

Quel est la méthode utilisée pour détecter les points d'intérêts?

[transistor]

c'est le détecteur de MORAVEC :on doit balayer l'image par une matrice de taille 3*3 et mesuré la différence entre les intensité des pixels est si elle est supérieur a un seuil donc c'est un point d'intérêt.

[transistor]

sil vous plait ca fait déja presque une semaine et toujours pas de réponce merci de m'aidé

[Jean Dumoncel]

c'est le détecteur de MORAVEC

[totalement HS]me demande si il a inspiré le nom des moravecs d'ILLIUM de Dan Simmons...[\totalement HS]

Bon après rapide lecture de cet algorithme, tu dois calculer les différences par rapport à un pixel de référence (celui que tu considères pour une fenêtre donnée) (comme dans ce fichier, diapo 41). Donc j'écrirai plutôt ceci :

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
I=imread('I.bmp');
Es=2; %le valeur sueil liée a la fonction
[N M]=size(I);

for j=1 : N-2 % Attention à ne pas dépasser la taille de ton image
    for i=1 : M-2
        for k=0 : 2
            for r=0 : 2
                E=E+sqr(abs((I(i+r,j+k))-(I(i,j)))) ;
            end ;
        end ;
        if E>Es
            text (i,j,'p','color','b') ;
        end ;
    end ;
end ;

Une première optimisation :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
I=imread('I.bmp');
Es=2; %le valeur sueil liée a la fonction
[N M]=size(I);
 
for j=1 : N-2
    for i=1 : M-2
        fenetre = ones(3,3).*I(i,j);
        E = sqrt(sum(abs((I(i:i+2,j:j+2)-fenetre))));
        if E>Es
            text (i,j,'p','color','b') ;
        end ;
    end ;
end ;

Non testé.

[Jerome Briot]

Une première optimisation :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
I=imread('I.bmp');
Es=2; %le valeur sueil liée a la fonction
[N M]=size(I);
 
for j=1 : N-2
    for i=1 : M-2
        fenetre = ones(3,3).*I(i,j);
        E = sqrt(sum(abs((I(i:i+2,j:j+2)-fenetre))));
        if E>Es
            text (i,j,'p','color','b') ;
        end ;
    end ;
end ;

Ou :

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
I=imread('I.bmp');
Es=2; %le valeur sueil liée a la fonction
[N M]=size(I);
 
fenetre = zeros(3);
 
for j=1 : N-2
    for i=1 : M-2
        fenetre(:) = I(i,j);
        E = sqrt(sum(abs((I(i:i+2,j:j+2)-fenetre))));
        if E>Es
            text (i,j,'p','color','b') ;
        end 
    end 
end

Mais tu n'aurais pas inverser la somme et la racine carrée par rapport au premier code ?

[Jean Dumoncel]

Ouais exact et même pas besoin de créer une matrice fenêtre, et puis la somme n'était pas faite correctement... (en même temps, c'est vendredi soir, j'ai une excuse).


Sauf nouvelle erreur :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
I=imread('I.bmp');
Es=2; %le valeur sueil liée a la fonction
[N M]=size(I);
 
for i=1 : N-2
    for j=1 : M-2
        E = sum(sum(sqrt(abs((I(i:i+2,j:j+2)-I(i,j))))));
        if E>Es
            text (j,i,'p','color','b') ;
        end 
    end 
end

[transistor]

merci a vous tous je test