Discussion :

[arni63]

Bonjour,
Est ce que quelqu'un peux m'aider à comprendre ce qu'on fait dans ces lignes de code? ça renvoie quoi? vecteur?matrice???

Merci

Code :

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

1
2
3
4
5
6
7
8
9
10
 
 
% VG        - template handle(s)
% VF        - handle of image to estimate params from
vx1 = sqrt(sum(VG(1).mat(1:3,1:3).^2));
vx2 = vx1;
 
kx = (pi*((1:k(1))'-1)/VG(1).dim(1)/vx1(1)).^2; ox=ones(k(1),1);
ky = (pi*((1:k(2))'-1)/VG(1).dim(2)/vx1(2)).^2; oy=ones(k(2),1);
kz = (pi*((1:k(3))'-1)/VG(1).dim(3)/vx1(3)).^2; oz=ones(k(3),1);

[arni63]

Code :

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

k    = max(round((VG(1).dim(1:3).*tmp)/cutoff),[1 1 1]);

un élément qui manquait pour comprendre le code

[vampirella]

Bonjour,

A première vue, ça m'a l'air d'un traitement d'image, mais impossible d'en dire plus sans le contexte (d'autant plus que les variables sont peu explicites).

Pourquoi ne pas travailler avec le débuggueur ? Tu pourras ainsi suivre ligne par ligne les exécutions des commandes, et comprendre sur un exemple concret ...

[arni63]

je ne peux pas utiliser le débogueur, vu que je ne peux pas faire tourner mon code, on m'a juste donné le code pour le détailler et donner une signification pour chaque résultat!!!
c'est une méthode de SPM pour faire du recalage d'images.
c'est juste un bout de code, sachant que les lignes d'avant on y a fait une Transformée en Cosinus Discrète.

[arni63]

Ma première partie du code:

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
% Number of basis functions for x, y & z
%-----------------------------------------------------------------------
tmp  = sqrt(sum(VG(1).mat(1:3,1:3).^2));
k    = max(round((VG(1).dim(1:3).*tmp)/cutoff),[1 1 1]);
 
% Scaling is to improve stability.
%-----------------------------------------------------------------------
stabilise = 8;  %on fait par bloc de 8*8
% on fait une DCT(transformée en cosinus discrète), ça nous donne la
% transformation??
basX = spm_dctmtx(VG(1).dim(1),k(1))*stabilise; 
basY = spm_dctmtx(VG(1).dim(2),k(2))*stabilise;
basZ = spm_dctmtx(VG(1).dim(3),k(3))*stabilise;
                  % j'obtiens les dérivées de la DCT, ça nous donne
                  % l'intensité?
dbasX = spm_dctmtx(VG(1).dim(1),k(1),'diff')*stabilise;
dbasY = spm_dctmtx(VG(1).dim(2),k(2),'diff')*stabilise;
dbasZ = spm_dctmtx(VG(1).dim(3),k(3),'diff')*stabilise;
 
vx1 = sqrt(sum(VG(1).mat(1:3,1:3).^2));
vx2 = vx1;
%   %définition des trois directions orthogonales???
kx = (pi*((1:k(1))'-1)/VG(1).dim(1)/vx1(1)).^2; ox=ones(k(1),1);
ky = (pi*((1:k(2))'-1)/VG(1).dim(2)/vx1(2)).^2; oy=ones(k(2),1);
kz = (pi*((1:k(3))'-1)/VG(1).dim(3)/vx1(3)).^2; oz=ones(k(3),1);

[arni63]

Code :

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

1
2
3
4
  %définition des trois directions orthogonales???
kx = (pi*((1:k(1))'-1)/VG(1).dim(1)/vx1(1)).^2; ox=ones(k(1),1);
ky = (pi*((1:k(2))'-1)/VG(1).dim(2)/vx1(2)).^2; oy=ones(k(2),1);
kz = (pi*((1:k(3))'-1)/VG(1).dim(3)/vx1(3)).^2; oz=ones(k(3),1);

est ce que à votre avis, ce qu'on fait là est une décomposition de l'image de base VG en sous matrice??