Discussion :

[Lost in]

Salut !!

J'ai un autre problème

J'essaye d'écrire un programme de binarisation par la méthode d'otsu:

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
 
for i=0:255
    hist(i)=0;
end;
for i=1:Hauteur
    for j=1:Largeur
        hist(img(i,j))=hist(img(i,j))+1;
    end;
end;
 
for i=0:255
    prob(i)=hist/(Hauteur*Largeur);
end;
omega(0)=prob(0);
myu(0)=0;
for i=1:255
    omega(i)=omega(i-1)+prob(i);
    myu(i)=myu(i-1)+ i*prob(i);
end;
seuil=0;
max_sigma=0;
for i=0:254
    if(omega(i)~= 0 & omega (i)~=1)
        sigma(i)=power(myu(255)*omega(i)-myu(i),2)/(omega(i)*(1- omega(i)));
    else sigma(i)=0;
    end;
    if(sigma(i)>max_sigma)
        max_sigma=sigma(i);
        seuil=i;
    end;
end;

ça ne marche pas à cause des indices qui commencent par 0 dans ce programme la valeur 0 "noir" est indispensable alors comment faire ?

A+

[Caro-Line]

Je ne vois pas bien le problème, tu n'as qu'à décaler tes indices de 1.

Code :

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

1
2
3
4
5
6
hist = zeros(1,256);%la fonction ZEROS est bien pratique
for i=1:Hauteur
    for j=1:Largeur
        hist(img(i,j)+1)=hist(img(i,j)+1)+1;
    end;
end;

En supposant que img n'a que des valeurs entières comprises entre 0 et 255 si je suis ton code (mais ça serait pas plutôt 254 ?).

Et pour la suite :

Code :

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

prob=ones(1,256)*hist/(Hauteur*Largeur);%la fonction ONES est bien aussi

[Lost in]

Merci caro95470

[Lost in]

Salut !!

L'autre jour lorsque j'ai posé le problème de binarisation je travaillais avec des images au format bmp, à présent je souhaiterais réutiliser le même programme pour des images dcm (16bits)

J'ai remplacé la valeur 255 du programme par 65536 mais ça n'a rien donné, bien au contraire, mon programme a bloqué

Pouvez vous m'aider SVP

[Jerome Briot]

J'ai remplacé la valeur 255 du programme par 65536 mais ça n'a rien donné, bien au contraire, mon programme a bloqué

Je pense que c'est simplement le temps d'exécution qui est plus long, non ?

Sinon, change le nom de la variable hist dans ton code. C'est déjà le nom d'une fonction MATLAB

[Jerome Briot]

J'essaye d'écrire un programme de binarisation par la méthode d'otsu:

La fonction GRAYTHRESH le fait très bien

[Lost in]

Je ne pense pas que ce soit une lenteur d'execution du programme

Le voici tu peux le tester si tu veux:

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
35
36
37
38
 
histo = zeros(1,65537);
for i=1:Hauteur
    for j=1:Largeur
        histo(img(i,j)+1)=histo(img(i,j)+1)+1;
    end;
end;
 
for i=1:65537 
    prob(i)=histo(i)/(Hauteur*Largeur);
end;
omega(1)=prob(1);
myu(1)=0;
for i=2:65537 
    omega(i)=omega(i-1)+prob(i);
    myu(i)=myu(i-1)+ i*prob(i);
end;
seuil=0;
max_sigma=0;
for i=1:65536 
    if(omega(i)~= 0 & omega (i)~=1)
        sigma(i)=power(myu(65535 )*omega(i)-myu(i),2)/(omega(i)*(1- omega(i)));
    else sigma(i)=0;
    end;
    if(sigma(i)>max_sigma)
        max_sigma=sigma(i);
        seuil=i;
    end;
end;
for i=1:Hauteur
    for j=1:Largeur
        if (img(i,j)> seuil)
            img(i,j)=0;
            else img(i,j)=65536 ;
                end;
    end;
end;
imshow(img);

A+

[Lost in]

Salut !!

Tu as raison Dut c'est le temps d'execution qui est très très lent

Il n'y a pas un moyen pour rendre l'execution plus rapide ?

A+

[Jerome Briot]

Pourquoi n'utilises-tu pas la fonction GRAYTHRESH ?

[Lost in]

Pourquoi n'utilises-tu pas la fonction GRAYTHRESH ?

Quand je l'ai testé pour des images bmp elle a donné de bons résultats, mais dès que je suis passé au format dicom j'obtiens toujours une image totalement blanche

A+