Discussion :

[slaima15]

Je débute en matlab, j'ai programmé des lignes de code cité ci dessous. Mais j'ai constaté que mon code et gourmand en temps de calcul:

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
 
 K=[0.25 0.5 1];
 INIT=K(1);
 I0=imresize(ImageReel,INIT,'nearest');
 s=size(I0);
  Speau=0.77;
  Sfond=0.22;
%M0 et P0 se sont des matrices binaires
for i=1:s(1)%parcour row
       for j=1:s(2)%parcour col
           if(i==1|j==1|i==s(1)|j==s(2))
            if (P0(i,j)==1)
                 M0(i,j)=1;
            end
           else
             E=energie(P0,i,j);%fonction qui calcule l'énergie de chaque pixel dans un bloc de taille 3*3 
          if(E>Speau|E<Sfond)
              M0(i,j)=0;
          else
              M0(i,j)=1;
          end
     end 
       end
 end

Pouvez–vous m’aidez à réecrire le code pour évitez le parcourt matricielle dans le but d'optimiser dans le temps ??

[Merel]

Salut,

Avant de faire des modifications, tu peux lancer un profil de ton code.

Code :

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

help profile

En gros, tu mets profile on en début de ton programme et une fois fini, tu tapes profview.

Code :

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

help profview

Ca te permettra de cibler les passages qui sont gourmand en temps de calcul

[slaima15]

oui j'ai essayé avec ça, toute la fonction consomme 0,6s

Donc mon but c'est ya-il une autre façon pour réecrire le code tout en évitant les deux boucle For(c-a-d le parcour de l'image pixel par pixel) dans le but d'optimiser le temps?

[Invité]

Bonsoir,

Une solution, si tu possèdes l'Image Processing Toolbox: utliser blkproc

Code :

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

E = blkproc(P0,[1 1],[1 1],@energie);

Par contre, la fonction energie devra être codée de façon à traiter chaque bloc 3x3 directement, tu n'auras plus besoin de i et j qui te positionneront dans la matrice.

Si tu ne possèdes pas l'Image Processing Toolbox,
une solution serai déjà de faire partir tes boucles à 2 et les faire terminer à s(1) -1 et s(2) -1, ce qui t'éconimiserai un if/else, traitant les bords en dehors des boucles.
De plus la fonction energie, ne faisant que quelques lignes, tu pourrais l'écrire directement dans la boucle, économisant ainsi l'appel d'une fonction.

Je ne vois pas d'autres modifications significatives, quelqu'un d'autre pourra peut-être

[slaima15]

merci Winjerome pour ton aide,
j'ai essayé avec ceci:

Code :

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

1
2
 
E = blkproc(P0,[1 1],[1 1],@energie);

mais j'ai pas compris comment écrire le code de la fonction energie
j'ai essayé avec ça:

Code :

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

1
2
3
 
function [E]=energie(P0)
E=(P0(1,1)^2+P0(1,2)^2+P0(1,3)^2+P0(2,1)^2+P0(2,2)^2+P0(2,3)^2+P0(3,1)^2+P0(3,2)^2+P0(3,3)^2)/9;

mais je sens que c'est faux.
En plus de ça avec ce blkproc, comment je peux faire ce test:

Code :

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

1
2
3
4
5
6
 
if(E>Speau|E<Sfond)
              M0(i,j)=0;
          else
              M0(i,j)=1;
          end

[Invité]

En plus de ça avec ce blkproc, comment je peux faire ce test:

Code :

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

1
2
3
4
5
if(E>Speau|E<Sfond)
    M0(i,j)=0;
else
    M0(i,j)=1;
end

Code :

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

1
2
E = blkproc(P0,[1 1],[1 1],@energie); 
M0 = (E<=Speau|E>=Sfond);

Ensuite:

Code :

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

1
2
function [E]=energie(P0)
E=(P0(1,1)^2+P0(1,2)^2+P0(1,3)^2+P0(2,1)^2+P0(2,2)^2+P0(2,3)^2+P0(3,1)^2+P0(3,2)^2+P0(3,3)^2)/9;

Est correct, mais j'écrirais plutôt:

Code :

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

E = sum(P0(:).^2)/9;

[slaima15]

Avec ceci

Code :

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

1
2
 
M0 = (E<=Speau|E>=Sfond);

Je peux pas inverser les conditions car ils me donne pas les même résultats, est ce que je peux exprimer ça autrement en conservant les conditions?

[Invité]

Avec ceci

Code :

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

M0 = (E<=Speau|E>=Sfond);

Je peux pas inverser les conditions car ils me donne pas les même résultats, est ce que je peux exprimer ça autrement en conservant les conditions?

Oui petite étourderie (il se fait tard ):

Code :

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

1
2
3
M0 = (E<=Speau & E>=Sfond); 
% Ou
M0 = ~(E>Speau|E<Sfond);