programmation matlab (boucles et sélection)

Version imprimable

Bonjour
J'essaye de simuler ces 2 étapes ss Matlab pour une utilisation ultérieure du résultat.

•1ère étape : Calcul de l’image moyenne M

Pour chaque Pixel G(i,j) de l’image G faire

Code:

1
2
3
4
5
 
           Si  G(i,j) est un Pixel de bordure alors 
                M(i,j)= G(i,j)
           Sinon   
            M(i,j)= [som de (G(m,n))]/9; avec m= i-1:i+1; n=j-1:j+1;)

• 2ème étape : calcul de l’image variance locale VL

Pour chaque Pixel G(i,j) de l’image G faire

Code:

1
2
3
 Si G(i,j) est un Pixel de bordure alors 
              VL(i,j) = G(i,j)
           Sinon VL(i,j)= [som de (carré (G(m,n)-M(i,j))]/9; avec m= i-1:i+1; n=j-1:j+1;)

• 3ème étape : calcul de l’image variance V

Code:

1
2
3
 Pour chaque Pixel VL(i,j) de l’image VL faire
 
                        V (i,j) =  30*log10(1+sqrt(VL(i,j)))

Exemple numérique :

Afin de simuler l'algorithme, Il prennent G une matrice 3*3,

Code:

G=[200 230 0;150 90 97; 160 40 255]

et ils obtiennent

Code:

1
2
3
4
	G(i,j)  =90 ;
	M(i,j)= (200+230+0+150+90+97+160+40+255)/9  = 135 (on prend le nombre entier)
	VL(i,j)= [(200-135)²+ (230-135)²+ (0-135)²+ (150-135)²+ (90-135)²+ (97-135)²+ (160-135)²+ (40-135)²+ (255-135)²] /9  =   7557
	      =  58

Voici le programme que j'ai écris sous MATLAB

Code:

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
clear
clc
close all
G=[200 230 0;150 90 97;160 40 255]
%Calcul de M, VL, V, moyenne, variance locale et variance  de l'image G
[L,C]=size(G);
for i=1:L
 
    for j=1:C
    % traitement des pixels de bord de G
    if i==1 | i==L |j==1 | j==C
        M(i,j)=G(i,j)
        VL(i,j)=G(i,j)
        V(i,j)=30*log10(1+sqrt(VL(i,j)))
    else
        % traitement des autres pixels de G 
        som1=0;
        som2=0;
        for k=i-1:i+1
            for l=j-1:j+1
        %pavet de 9 pixels
        som1=som1+ G(k,l)
        % calcul de la matrice de variance locale sur des pavets de 3*3
        % pixel centré sur le pixel P(i,j)
        som2=som2+ (G(k,l)-M(i,j))^2
            end
        end
        M(i,j)=floor(som1/9)
        VL(i,j)=floor(som2/9)
        V(i,j)=30*log10(1+sqrt(VL(i,j)))
   end
  %keyboard
   end 
end

Ceci me donne le bon résultat pour M(i,j) mais pas la bonne variance pour L(i,j)=VL(2,2). Je ne trouve pas 7557.
Où est le problème? Je ne fais que recommencer la compilation mais vainement, voudriez vous me dire où est l’erreur dans mon code source MATLAB?

Mes remerciements

Essaie ceci :

Code:

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
 
% Initilaisation
    clc ;
    clear all ;
    G = [200 230   0  ; ...
         150  90  97  ; ...
         160  40 255] ;
 
    [nL nC] = size(G) ;
 
% Traitement
    M  = zeros(size(G)) ;
    VL = zeros(size(G)) ;
    for i=1:nL 
        for j=1:nC
            % Vérification si on est sur les bords de l'image
            if i==1 || i==nL || j==1 || j==nC
                M(i, j)  = G(i, j) ;
                VL(i, j) = G(i, j) ;
            else
                M(i, j)  = sum(sum(G(i-1:i+1, j-1:j+1)))/9 ;
                VL(i, j) = sum(sum( (G(i-1:i+1, j-1:j+1)-M(i, j)).^2)) ;
            end
        end
    end
 
    VL(2,2)

Bonjour
Tout d'abord grand merci pour toujours donné suite à mes questions
dans la formule de VL il manquait /9 que j'ai rajouté. Les résultats doivent être entiers et dc j'ai pensé utiliser la fonction floor comme vous pourriez le voir ci dessous. J'ai aussi recalculé la quantité VL(2,2) selon l'algorithme qu'ils ont donné. le résultat est affecté à la variable r.
Code:

1 2 3 4 5 6 7 8 9 10 M(i, j) = sum(sum(G(i-1:i+1, j-1:j+1)))/9 ; M(i,j)=floor(M(i,j)) VL(i, j) = sum(sum( (G(i-1:i+1, j-1:j+1)-M(i, j)).^2))/9 ; VL(i,j)=floor(VL(i,j)) end end end VL(2,2) r= ((200-135)^2+ (230-135)^2+ (0-135)^2+ (150-135)^2+ (90-135)^2+ (97-135)^2+ (160-135)^2+ (40-135)^2+ (255-135)^2)/9
en utilisant la fonction floor on trouve exactement le mm résultat puisque la différence est ds la partie décimale, (mais c différent de ce qu'il disent avoir trouvé 7557)
sans la fonction floor, je trouve respectivement
VL(2,2)=ans =

6.5793e+003

r =

6.5799e+003

que pourrait on dire??
Mes remerciements

21/02/2010, 15h31
kmaniche

Et maintenant, tu as les bons résultats ou il existe toujours une erreur de résultat ?
21/02/2010, 16h07
maraval

boucle et sélection

justement voici ce que je trouve

d'après le code Matlab
VL(2,2)=ans =

6.5793e+003

d'après leur algorithme
r =

6.5799e+003

il y a y une petite différence dans la partie décimale (en gras). Comment est jugé le résultat du code Matlab (juste ou faux) , puisque la fonction sum n'utilise pas du calcul approché. Les étapes de calcul sont similaire à une somme ordinaire de termes.
En fait floor donne la partie entière d'un nombre. Elle est la seule???
21/02/2010, 16h13
maraval

en fait :-)! un help de floor donne ceil, round, et fix.
fix et floor posent une certaine cnfusion. Pourrais je avoir un exemple numérique s'il vous plait
21/02/2010, 16h24
kmaniche

FAQ MATLAB >>

http://matlab.developpez.com/faq/?page=mat#mat_decimale
21/02/2010, 16h54
maraval

boucle et sélection

Merci pour cette réorientation. Je suis servie avec des exemples et tout. Parfait!
Pour ce qui est du code Matlab, comment sont jugés les résultats obtenus avec sum et directement avec l'expression mathsématique
21/02/2010, 19h13
kmaniche
Soit une matrice a définit entre aléatoirement entre [0 et 10] :
Code:

1 2 3 4 5 6 7 8 >> a=fix(rand(4,5)*10) a = 8 6 9 9 4 9 0 9 4 9 1 2 1 8 7 9 5 9 1 9
Ci-après te donne la somme des cellules colonnes par colonnes :
Code:

1 2 3 4 5 >> sum(a) ans = 27 13 28 22 29
Ci-après te donne la somme des cellules lignes par lignes :
Code:

1 2 3 4 5 6 7 8 9 10 >> sum(a,2) ans = 36 31 19 33 >>
Ci-après te donne la somme de toutes les cellules
Code:

1 2 3 4 5 6 7 8 9 10 11 >> sum(sum(a)) % est la même que sum(a(:)) ans = 119 >> sum(a(:)) ans = 119
Je penses que c'est simple, sinon si tu veux plus, fais moi signe :wink:

Pourquoi tu appelles floor à chaque itération et sans ";" :

Code:

1
2
3
4
5
6
7
8
9
10
 
                M(i, j)  = sum(sum(G(i-1:i+1, j-1:j+1)))/9 ;
                M(i,j)=floor(M(i,j))
                VL(i, j) = sum(sum( (G(i-1:i+1, j-1:j+1)-M(i, j)).^2))/9 ;
                VL(i,j)=floor(VL(i,j))
            end
        end
    end
 
    VL(2,2)

Plus simple

Code:

1
2
3
4
5
6
7
8
9
10
11
 
                M(i, j)  = sum(sum(G(i-1:i+1, j-1:j+1)))/9 ;
                VL(i, j) = sum(sum( (G(i-1:i+1, j-1:j+1)-M(i, j)).^2))/9  ;
            end
        end
    end
 
M  = floor(M) ;
 VL = floor(VL) ;
 
sprintf('La valeur de VL(%d, %d) = %d', 2, 2, VL(2,2))

21/02/2010, 21h42
maraval

boucle et sélection

bonsoir
Oui en effet, grâce à vos consignes et çà devient très clair.
Pour floor c'est plutôt fix pour n'en extraire que la partie entière :-)
J'ai pas mis de ; pour visualiser le résultat. Vous m'avez donné mieux 'fprintf' et c parfait.
Mes remerciements les plus sincères!
21/02/2010, 21h52
maraval

affichage des résultats

En fait, comme pour les fonctions floor, ceil et autres y aurait il une FAQ qui renvoie vers des exemples d'illustration pour les fonctions d'affichage fprintf, sprintf, scanf ...., afin de les utiliser avec le meilleur des usages?
21/02/2010, 22h14
kmaniche
Avec Matlab, tu as une documentation très très riche. profitez-on alors via :
Code:

1 2 >> doc sprintf
21/02/2010, 22h27
maraval

affichage des résultats

j'y suis déjà :-)
Merci @ vous et bonne soirée!