Dabeen Park/Matlab/MRI
From Biolecture.org
x=TR;
S0=0.0001;
T1=0.0002;
resultS0=zeros(128,128);
resultT1=zeros(128,128);
for i=1:128
for j=1:128
y=squeeze(img(i,j,1:16));
result=fit(x',y,'S0*(1-exp(-x/T1))', 'StartPoint', [S0,T1]);
resultS0(i,j)=result.S0;
resultT1(i,j)=result.T1;
end
end
surf(resultS0)
view(2)
