Discussion :

[davidmorice]

Bonjour,
Je suis nouveau sur le forum et je suis coincé sur MATLAB
Je travaille sur image binarisé qui représente des porosité dans un matériaux
Pour illustrer, l'image que j'ai est du type :

Code :

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

1
2
3
4
I = 0 0 1 0 1 0 1 
    0 1 0 1 0 1 0
    0 0 0 0 0 0 1
    1 0 1 0 1 1 0

Je souhaiterai créer un programme qui calcule la densité de porosité c'est à dire pour chaque point de l'image, calculer le nombre de 1 et le diviser par le nombre total.

Exemple avec une "boite d'étude de 3*3"

J(2,2) =  (0+0+1+
           0+1+0+
           0+0+0)/9 = 0,22 soit 22%
J= 0,11 0,22 0,33 ...
   0,11 0,22 0,33 ...
   0,22 0,33 0,33 ...
   ...

Voila ce que j'ai fait actuellement:

Code :

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

1
2
3
4
5
6
7
TX=3;
TY=TX;
for X = 1+TX:size(I,1)-TX;
    for Y = 1+TY:size(I,2)-TY;
        J(X-TX,Y-TY)= 100 - (sum(sum(I(X-TX:X+TX,Y-TY:Y+TY)))/((2*TX+1)*(2*TY+1))*100);
    end
end

Le résultat que j'obtiens n'est pas bon,
merci d'avance de bien vouloir m'aider!

PS: j'espère que mon problème est clair, n'hésiter pas à me dire si vous ne m'avaez pas compris

[Jerome Briot]

Que devrait selon toi retourner ceci ?

Code :

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

X = 1+TX:size(I,1)-TX

Pour ma part :

Code :

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

1
2
3
4
5
>> X = 1+TX:size(I,1)-TX
 
X =
 
   Empty matrix: 1-by-0

car

Code :

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

1
2
3
4
5
6
7
8
9
10
11
>> 1+TX
 
ans =
 
     4
 
>> size(I,1)-TX
 
ans =
 
     1

Sinon ton problème se résout avec :

Code :

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

1
2
3
4
5
6
I = [0 0 1 0 1 0 1
    0 1 0 1 0 1 0
    0 0 0 0 0 0 1
    1 0 1 0 1 1 0]
 
J = filter2(ones(3),I,'same')/9;

ce qui donne :

J =

    0.1111    0.2222    0.3333    0.3333    0.3333    0.3333    0.2222
    0.1111    0.2222    0.3333    0.3333    0.3333    0.4444    0.3333
    0.2222    0.3333    0.3333    0.3333    0.4444    0.4444    0.3333
    0.1111    0.2222    0.1111    0.2222    0.2222    0.3333    0.2222

[davidmorice]

en écrivant ceci, je voulais "balayer" entre 1+TX et size(I,1) avec un pas de 1

Merci en tout cas pour ta réponse rapide! Ca marche nikel!

[Jerome Briot]

en écrivant ceci, je voulais "balayer" entre 1+TX et size(I,1) avec un pas de 1

Dans ce cas, le -TX à la fin de l'expression ne sert à rien :

Code :

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

X = 1+TX:size(I,1)