Discussion :

[Lost in]

Salut,

Je souhaiterais optimiser le code suivant:

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
39
40
41
42
43
44
45
46
 
 
im = double(im);
B= double(im);
[m n]= size(im);
 
for t=1:10
im=B;
for i=1:m-2
    for j=1:n-2
       nablaNW=abs(im(i,j)-im(i+1,j+1));
       nablaN=abs(im(i,j+1)-im(i+1,j+1));
       nablaNE=abs(im(i,j+2)-im(i+1,j+1));
       nablaW=abs(im(i+1,j)-im(i+1,j+1));
       nablaE=abs(im(i+1,j+2)-im(i+1,j+1));
       nablaSW=abs(im(i+2,j)-im(i+1,j+1));
       nablaS=abs(im(i+2,j+1)-im(i+1,j+1));
       nablaSE=abs(im(i+2,j+2)-im(i+1,j+1));
 
       A=[nablaN nablaS nablaW nablaE nablaNE nablaSE nablaSW nablaNW];
       M=max(A);
 
       lapl=4*im(i+1,j+1)-im(i+1,j+2)-im(i+1,j)-im(i,j+1)-im(i+2,j+1);
       if (M ~= 0)
            cN = max(exp(lapl),1+nablaN);
            cS = max(exp(lapl),1+nablaS);
            cW = max(exp(lapl),1+nablaW);
            cE = max(exp(lapl),1+nablaE);
            cNE =max(exp(lapl),1+nablaNE);
            cSE =max(exp(lapl),1+nablaSE);
            cSW =max(exp(lapl),1+nablaSW);
            cNW =max(exp(lapl),1+nablaNW);
 
 
 
             B(i+1,j+1)= B(i+1,j+1) + ...
                  (...
                  cN.*nablaN + cS.*nablaS + ...
                  cW.*nablaW + cE.*nablaE + ...
                  cNE.*nablaNE + cSE.*nablaSE + ...
                  cSW.*nablaSW + cNW.*nablaNW );
       end;
    end;
 
end;
end;

Avec les 3 boucles imbriquées j'attends un bon bout de temps avant de voir le résultat

Je vous prie de bien vouloir m'aider

[Caro-Line]

Pour avoir une idée de départ et tester c'est quoi la taille de ta matrice im ?

[Lost in]

Pour avoir une idée de départ et tester c'est quoi la taille de ta matrice im ?

480 x 456

[Caro-Line]

Ben écoute, finalement les boucles c'est pas si mal que ça
J'ai tout "condensé", ça fait un code joli et incompréhensible, il n'y a presque plus rien dans les boucles mais c'est 10 fois plus long que ton code initial
Donc les boucles sont bien optimisées sous MATLAB.

Je te donne mon code pour que tu puisses comparer :

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
im = double(im1);
B= im;
[m n]= size(im);
 
for t=1:10
    im=B;
    %initialisation
    nabla=zeros(m-2,n-2,8);
    %Calcul de tous les nabla
    nabla(:,:,1)=abs(im(1:m-2,1:n-2)-im(2:m-1,2:n-1));
    nabla(:,:,2)=abs(im(1:m-2,2:n-1)-im(2:m-1,2:n-1));
    nabla(:,:,3)=abs(im(1:m-2,3:n)-im(2:m-1,2:n-1));
    nabla(:,:,4)=abs(im(2:m-1,1:n-2)-im(2:m-1,2:n-1));
    nabla(:,:,5)=abs(im(2:m-1,3:n)-im(2:m-1,2:n-1));
    nabla(:,:,6)=abs(im(3:m,1:n-2)-im(2:m-1,2:n-1));
    nabla(:,:,7)=abs(im(3:m,2:n-1)-im(2:m-1,2:n-1));
    nabla(:,:,8)=abs(im(3:m,3:n)-im(2:m-1,2:n-1));
    lapl_exp=exp(4*im(2:m-1,2:n-1)-im(2:m-1,3:n)-im(2:m-1,1:n-2)-im(1:m-2,2:n-1)-im(3:m,2:n-1));
    %les maximums
    MAXI=max(nabla,[],3);
    for i=1:m-2
        for j=1:n-2
            if (MAXI(i,j) ~= 0)
                B(i+1,j+1)= B(i+1,j+1) + sum(max(lapl_exp(i,j),1+nabla(i,j,:)).*nabla(i,j,:));
            end
        end
 
    end
end

Peut-être est-il améliorable....Je ne pense pas pouvoir enlever la dernière double boucle.

[Lost in]

Merci caro95470 je vais le tester tout de suite