Having read the previous answer I have come up with the code I put underneath (256x256 matrix, 11 inversion times, 9 slices)
this is the error I have been getting :
Error using fit>iFit (line 135)
X must be a matrix with one or two columns.
Error in fit (line 116)
[fitobj, goodness, output, convmsg] = iFit( xdatain, ydatain, fittypeobj, ...
Error in T1map (line 9)
T1f=fit(TI,signal,eq);
anyone have an idea where I am going wrong?
re_im=reshape(im,256,256,11,9);
TI=[0.01 0.03 0.4 0.5 1 1.25 1.5 1.75 3 5 9]
eq=fittype('a*(1-2*exp(-x/T1))','coefficients',{'a','T1'});
x=zeros(128,128,9);
for row=[1:256]
for col=[1:256]
for im =[1:15]
signal=re_im(row,col,:,im);
T1f=fit(TI,signal,eq);
x(row,col,im)=T1f;
end
end
end