Discussion :

[Imed Soufi]

bonjour, voila avec ce code je dois calculer les angles de chaque trait, par exemple si j'ai deux traits je dois couper mon image en bloc et calculer le gradient puis estimer mon angle, mais voila j'aimerais que mon résultat soi mis dans un vecteur T, le problème c que les deux trais leurs angle respectivement sont de 37° et 33° alors dans mon vecteur il y a deux valeurs de 73° 73° .
quelqu'un pourrait-il m'éclairer svp.

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
% Initialisations
clear all;
close all;clc ;
% Lecture de l'image originale.
I = imread(' im.bmp');
I=double(I);
[n,m]=size(I);
s3=0;s4=0;
w=20;
f=round(m/w);
for i=1:n-2
    for j=1:w:m-2
        s1=-I(i,j)+I(i,j+2)-2*I(i+1,j)+2*I(i+1,j+2)-I(i+2,j)+I(i+2,j+2);
        s2=I(i,j)+2*I(i,j+1)+I(i,j+2)-I(i+2,j)-2*I(i+2,j+1)-I(i+2,j+2);
        s3=s3-2*s1*s2;
        s4=s4+(s2^2-s1^2);
 
 
        s5=0.5*atand(s3/s4);
        if s4<0 
            s5=s5+90;
        elseif s4>=0 
            if s3<0
                s5=s5+180;
            else
            end
        else
        end
        for p=1:f
                T(p)=s5;
        end
    end 
end
T

[Caro-Line]

Personnellement je ne comprends pas bien ton code.
Tu as 2 boucles imbriquées et donc pour les valeurs de T tu n'auras que celles calculées pour i=n-2 et j=m-2.
De plus tu fais :

Code :

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

T(p)=s5;

Donc quelque soit p, T(p) vaudra toujours s5.
Normal que ton vecteur contienne les mêmes valeurs partout.

conseil : tu as des else sans rien derrière, tu pourrais les enlever ça améliorerait la lecture du code.