Discussion :

[matt67]

Bonjour,

Pourriez-vous m'indiquer une méthode matricielle plus rapide de mon code ci-dessous utilisant par défaut des boucles, qui prend donc un temps infini... :

Code :

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

1
2
3
4
5
6
7
8
9
FiberFA{MaxSizeFibers} = [];
 
for i = 1 : MaxSizeFibers
    for j = 1 : NfibersROI 
        if length(ROIfibersFA.fiber(j).FA) >= MaxSizeFibers
            FiberFA{i} = [ FiberFA{i} ; ROIfibersFA.fiber(j).FA(i) ];
        end
    end
end

"MaxSizeFibers" et "NfibersROI" sont deux entiers. "ROIfibersFA.fiber" est un vecteur de structure de taille "NfibersROI" et "ROIfibersFA.fiber(j).FA" est un vecteur d'entiers de taille au maximum égale à "MaxSizeFibers".

Merci d'avance pour votre aide !

Matthieu

[Gooby]

Bonjour

Tu peux essayer d'utiliser ceci:

Code :

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

1
2
3
A=struct2cell(ROIfibersFA.fiber);
longueur = cellfun(@length,A);
idx = find(longueur>=MaxSizeFibers);

Il existe probablement une multitude d'autres méthodes.

Bon ensuite comment se passer de la boucle sur i, dans l'immédiat je ne vois pas bien parce que j'ai dû mal à me représenter le contenu des variables. ça doit être possible également.

[matt67]

Bonjour,

J'ai essayé le code suivant mais la boucle sur i est interminable lorsque je teste "MaxSizeFibers" = 36000

Code :

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

1
2
3
4
5
6
7
8
9
10
11
[MaxSizeFibers,index] = max(cat(1,RightFCS_FA.fiber.nFiberLength));
 
length_FA = arrayfun(@(x) length(x.FA),RightFCS_FA.fiber);
 
FiberFA{MaxSizeFibers} = [];
 
idx = arrayfun(@(x) find(x<=length_FA), 1:MaxSizeFibers, 'UniformOutput', false);
 
for i = 1 : MaxSizeFibers
    FiberFA{i} = arrayfun(@(x) x.FA(i), RightFCS_FA.fiber(idx{i}));
end

Avez-vous une idée d'optimisation pour la boucle i ?

Merci !

Matthieu