How to create a T1 Map

Question

0 个投票

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

2 个评论
显示无隐藏无

Lydia Smith 2019-6-13

在 MATLAB Online 中打开

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