How to create a T1 Map

Hello,

I am reasonably new to matlab and was wondering if anyone had any ideas on how to create a T1 map for MRI images.

I know I need to create a for loop so that it applys the T1 equation to every pixel in the image, and have a matrix of 0's to put the data into however I am not sure how to put this into action.

I have created the for loop but the equation I have used previously to create an inversion recovery curve does not seem to be compatible. woudl anyone know the equation that I can apply per pixel to create the T1 map?

Thanks

[Lydia Smith]

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

[Nasrin Akter]

Hello

Can you please share where you found the MRI image?

Answer 1

[Prasanth Sikakollu]

Apply the equations in the article attached hereby for each pixel to create the T1 Map using inversion recovery.

For MATLAB script to create a T1 map, have a look at the following link.

http://www.mccauslandcenter.sc.edu/CRNL/teaching/mri

Hope these resources help.

[Walter Roberson]

Wayback Machine archived that second link at https://web.archive.org/web/20180511185434/http://www.mccauslandcenter.sc.edu/CRNL/teaching/mri

The scripts referred to are at https://github.com/rordenlab/spmScripts

The page was part of a course; the current material for the course appears to be at https://crnl.readthedocs.io/psyc589888/index.html