Discussion :

[kmaniche]

Bonsoir,

J'ai une image 3D Im3D dont les dimensions sont :
/x = sx
/y = sy
/z = sz

je cherche à construire une image 2D tel que chaque element de cette image soit égale au maximum selon Z (profondeur).

Particulièrement, utilisant la méthode de vectorisation :

Code :

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

1
2
3
 
Im2D = zeros(sx, sy) ;
Im2D(:, :) = max(Im3D(:, :, 1:sz)) ; % <<== Là ça ne va pas ??

Votre aide svp.

[Jerome Briot]

Essaie ceci :

Code :

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

Im2D = max(Im3D,3)

[kmaniche]

Bonjour,


Je vous remercie pour la réponse, mais cela ne semble le pas marcher. Pourquoi, parceque ici il va comparer chaque élement de la matrice 3D (Im3D) par rapport à la valeur 3, c'est à dir tous les élements de cette matrice qui sont inferieurs à 3 (0:2), le résultat de retour sera 3.

Pour bien vous expliquer, voici une façon iterative :

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
 
clear all ;
close all ;
clc ;
 
load mri ; % pour l'exemple ;)
im3D = squeeze(D) ;
 
im2Dmx = zeros(siz(1), siz(2)) ;
 
for i=1:siz(1)
    for j=1:siz(2)
        mx = im3D(j, i, 1) ;
        for k=1:siz(3)
            if mx<im3D(j, i, k)
                mx = im3D(j, i, k) ;
            end
        end
        im2Dmx(j, i) = mx ;
    end
end

je cherche à vectoriser cela, genre :

Code :

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

1
2
 
im2Dmx = max(im3D) ;

Il est logique que im2Dmx est une matrice 2D (128x128).

Que pensez-vous ?

[Jerome Briot]

Je suis allé un peu vite dans ma réponse

Mais une petite lecture de la documentation donnait :

[Y,I] = MAX(X,[],DIM) operates along the dimension DIM.

Donc :

Code :

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

Im2D = max(Im3D, [] ,3)

[kmaniche]

J'ai déjà essayé cette solution (doc max), mais avrait dire que le résultat n'est pas le même que celui des boucles ?

Etrange non !!

[kmaniche]

Zut, une erreur dans le code du post (la boucle for de k) j'ai corrigé .


ça ne marche toujours pas ! ! !

[Jerome Briot]

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
clear all ;
close all ;
clc ;
 
load mri ; % pour l'exemple ;)
im3D = squeeze(D) ;
 
im2Dmx = zeros(siz(1), siz(2)) ;
 
for i=1:siz(1)
    for j=1:siz(2)
        mx = im3D(j, i, 1) ;
        for k=1:siz(3)
            if mx<im3D(j, i, k)
                mx = im3D(j, i, k) ;
            end
        end
        im2Dmx(j, i) = mx ;
    end
end
 
im2Dmx2 = max(im3D,[],3);
 
all(im2Dmx(:)==im2Dmx2(:))

Code :

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

1
2
3
4
 
ans =
 
     1

[kmaniche]

Interéssnte cette fonction "all",

OK cela fonctionne, moi j'ai eu à utiliser ce code pour voir si les deux matrice sont égales :

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
 
clear all ;
close all ;
clc ;
 
load mri ; % pour l'exemple ;)
im3D = squeeze(D) ;
 
im2Dmx = zeros(siz(1), siz(2)) ;
 
for i=1:siz(1)
    for j=1:siz(2)
        mx = im3D(j, i, 1) ;
        for k=1:siz(3)
            if mx<im3D(j, i, k)
                mx = im3D(j, i, k) ;
            end
        end
        im2Dmx(j, i) = mx ;
    end
end
 
sum im2Dmx2
sum im2Dmx

le résultat (les deux sum) n'étant pas le même , je sais pas pourquoi ?

[kmaniche]

Bonjour,

Merci pour la solution, ça marche parfaitement, j'ai fait le test avec deux matrices de 3x3x3. J'ai eu le resultat escompté.

Merci encore une fois.