Discussion :

[arni63]

Bonjour,

Code :

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

1
2
3
4
5
6
7
8
9
10
11
cP = ['/d1/data/Petra/rat5/canonical_T2.img';...
      '/d1/data/Petra/LFF4/canonical_T2.img';...
      '/d1/data/Petra/rat3/canonical_T2.img';...
      '/d1/data/Petra/rat4/canonical_T2.img';...
      '/d1/data/Petra/rat6/canonical_T2.img'];
cP = spm_vol(cP);%on recup un vecteur de structures
cdata = spm_read_vols(cP);% cdata de la forme [ Y, XYZ], on aura en sortie(matrice 4D des données d'image, matrice 3*N des X Y Z de retour
for i=1:length(cP)% pour chaque structure de cP
   maxval(i) = max(max(max(cdata(:,:,:,i)))); 
   cdata(:,:,:,i) = cdata(:,:,:,i) / maxval(i);
end

que vaut dans ce cas

Code :

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

maxval(i) = max(max(max(cdata(:,:,:,i))));

Merci

[Jean Dumoncel]

Bonjour,

Code :

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

help max

[arni63]

si je comprends bien alors:
je prends à chaque fois le max sur ma matrice 4*4et je remplace les valeurs minimales successivement par la valeur de i

c'est bien ça ou pas??

[Jean Dumoncel]

Non, tu prends la valeur maximum de cdata(:,:,:,i) et tu normalises cdata(:,:,:,i) par cette valeur (la plus grande valeur de cdata(:,:,:,i) deviendra 1), et tu le fait pour chaque i.

[arni63]

Merci.
j'avais comrpis qu'il s'agissait de normalisation.
en fait ce qui m'echape c'est le fonctionnement du max(vect(:,:,1)).
comment réussir à le faire toute seule sur une matrice.

je vais essayer encore de chercher, ça m'embete de comprendre des choses simple comme ça, pour mieux comprendre mon code!!!

[Jean Dumoncel]

Un exemple :

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
>> g = rand(4,4,3)

g(:,:,1) =

    0.4634    0.8760    0.9790    0.7575
    0.3362    0.4877    0.8633    0.8974
    0.3091    0.3051    0.3307    0.2394
    0.4017    0.4940    0.4704    0.7895


g(:,:,2) =

    0.4674    0.4389    0.1414    0.5822
    0.8565    0.4621    0.5368    0.2969
    0.4597    0.0311    0.3823    0.3852
    0.6724    0.1569    0.4124    0.6749


g(:,:,3) =

    0.1552    0.2664    0.8727    0.5646
    0.3005    0.2764    0.6728    0.0440
    0.3689    0.1433    0.2698    0.6422
    0.6436    0.3400    0.1716    0.7154

Si on recherche le max de g(:,:,2) (en rouge):

Code :

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

max(max(g(:,:,2)))

On applique le max plusieurs fois car la fonction max renvoie le maximum de chaque colonne lorsqu'on lui met une matrice en entrée. (donc dans ce cas, il faut à nouveau prendre le max des max de chaque colonne)

Dans ton cas tu as une matrice 4-dimensions, donc chaque cdata(:,:,:,i) est une matrice 3-D, c'est pourquoi on applique 3 fois la fonction max pour obtenir le maximum de cette matrice 3-D.