How can I optimize the DCE-MRI fitting to obtain a parameter map?

Hi everybody! I'm working on a script for my master thesis. I have to extract parameters from DCE-MRI data. The aim is to fit data with known functions for every pixel of the MRI image.

I give you one example: let us take one compartimental model, as the Tofts one (attached image).

From the DCE-MRI acquisitions I know the terms Cp(t) and Ct(t) for every t; the unknow parameters are Ktrans and ve. They must be fitted with the known function which include the operation of convolution between Cp(t) and Ht(t).

If we can imagine to have a single pixel, we have, for example, t = 65 and so we have 65 data for Cp and for Ct to do the fitting.

The fitting must be done for every pixel of the image; in my case I have 16 slices, each with 256x256 pixels. As you can imagine, this process takes a lot of time.

Can you suggest a way to do the same thing faster? Thank you in advance!

Answer 1

[Simone Lorenzi]

up

[Carlos_D]

Hi Simone Lorenzi

I am also working with DCE-MRI images although I am starting. Could you share the script to adjust that data since there is none in Matlab or the information to do so?

My best wishes