Discussion :

[Alucard9800XT]

Bonjour ou Bonsoir
Bon pour mon programme qui est censé calculé le nombre de pixel pour chaque niveau de gris d'une images ' sur 256 niveaux de gris ' mon probleme c'est ce message d'erreur

Code :

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

1
2
3
4
??? Attempted to access I(1,234); index out of bounds because size(I)=[238,233].
 
Error in ==> histogram at 9
            if I(i,j) == v(k)

et voici mon programme

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
v=(0:1:255)
z=v*0
I=X;
[n,m,t]=size(I);
for i=1:m
    for j=1:n
        for l=1:t
        for k=1:256
            if I(i,j) == v(k)
                z(k)=z(k)+1;
            end
        end
        end
    end
end

sachant que X est une image couleur pour mon cas charger par matlab

dois je prendre en compte l'image en matrice 3dimensions RVB?
et merci

[Alucard9800XT]

mais toujour erroné du poin de vue Matlab

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
 
v=(0:1:255)
z=v*0
I=X;
[n,m,t]=size(I);
for i=1:m
    for j=1:n
        for l=1:t
        for k=1:256
            if I(i,j,t) == v(k)
                z(k)=z(k)+1;
            end
        end
        end
    end
end

[rostomus]

Bonjour,


Lis la FAQ: probleme d'indixage

Bon courage

[ol9245]

Tu peux supprimer au moins deux boucles facilement.
peut-être la troisième avec plus d'astuce que j'en ai.

Code :

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

1
2
3
4
z = [] ;
for i = 1:255
z(i) = sum(sum(I==i)) ;
end

pour comprendre mon petit programme tape en ligne de commande

Code :

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

magic(4)==3 ;

OL

[Jerome Briot]

Il faut utiliser la fonction HISTC ici:

Code :

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

1
2
3
IM=round(rand(10)*255);
N=histc(IM(:),0:255);
[[0:255].' N]

-----------------------------------------

Pour ol9245,

Code :

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

1
2
3
z = [] ;
for i = 1:255
...

Lire la FAQ à propos de la préallocation de matrice

Code :

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

z(i) = sum(sum(I==i)) ;

Ou simplement

Code :

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

z(i)=sum(I(:)==i);

[ol9245]

Cool. je me doutais bien que la troisième boucle pouvait sauter !
OL

[larimoise]

On peut aussi trouver l'histogramme d'une image avec la commande imhist dans l'Image Processing toolbox