Discussion :

[coolzy]

Bonjour,

j'ai une séquence d'images IRM dans le temps, et j'étudie les courbes d'intensité de chaque pixel pour effectuer certains calculs. Mais avant de faire ces calculs, ces courbes doivent être approximées avec l'équation d'une courbe gamma de la forme:

C(t) = K.t^alpha.exp(-t/beta)

Avec C(t) l'intensité lumineuse du pixel au temps t.

Je dois donc à partir des valeurs de mon pixel, trouver les trois paramètres K, alpha et beta pour faire coller cette courbe à cette equation.
J'utilise donc nlinfit, je l'avais deja utilisé pour une autre équation et cela marchait bien, mais ici j'ai sans arrêt l'erreur suivante:

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
??? Error using ==> nlinfit>checkFunVals at 356
MODELFUN has returned Inf or NaN values.
 
Error in ==> nlinfit>LMfit at 305
    if funValCheck && ~isfinite(sse), checkFunVals(r); end
 
Error in ==> nlinfit at 168
    [beta,J,lsiter,cause] = LMfit(X,y, model,beta,options,verbose,maxiter);
 
Error in ==> project>buttonperfusion_Callback at 3226
            [BETAt] = nlinfit(timeline(1:tr-t0+1)',C_t(i,1:tr-t0+1)',@gamma_fit, BETA0);
 
Error in ==> gui_mainfcn at 95
        feval(varargin{:});
 
Error in ==> project at 64
    gui_mainfcn(gui_State, varargin{:});
 
??? Error using ==> project('buttonperfusion_Callback',gcbo,[],guidata(gcbo))
Error using ==> nlinfit>checkFunVals at 356
MODELFUN has returned Inf or NaN values.
 
??? Error while evaluating uicontrol Callback

comment puis-je m'en sortir?

[Giansolo]

Salut,

Même problème ici sous Octave, avec la fonction leasqr (fonctions identiques ou très proches sans doute) et des retours NaN...
Grrr...

[coolzy]

toujours pas de solution? J'aimerai juste gerer l'erreure, du style dire 'si nlinfit renvoie cette erreure, continuez sur les autres pixels comme si de rien n'etait', pour que la boucle ne soit pas interrompue.

[tug83]

Pour gerer cette erreur, vovii une idée:

Code :

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

1
2
3
4
5
6
try
   do_something
catch
   do_cleanup
   rethrow(lasterror)
end

[coolzy]

ok oui j'utilisais deja un try, mais pas tres bien place, jvais arranger tout ca. merci!

[coolzy]

desolee mais je ne vois pas trop comment m'en sortir.
j'ai le code general:

Code :

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

1
2
3
4
5
6
7
8
9
10
 
for i = 1:n
 
           try
               [BETA(:,i),R,J] = nlinfit(timeline, Cbis(i,:), @tofts_tissue_conc_func, BETA0);
                catch
                   rethrow(lasterror)    
                end
 
end

mais ca ne change pas grand chose. je ne sais pas comment lui dire: si il y a une erreure, tant pis mets BETA(:,i) a zero et continue la boucle.