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?