Fit T2* values using Python

Fit T2* values using Python

For a project I needed T2*-values for different tissues. These can be estimated by using T2-weighted images with different echo times (TEs). Some free tools exist like the Mipav-based CBS-tools and the web-based MRI Toolbox, but the Mipav toolkit is rather unstable (Java, ...) and needs everything to be in a very specific format. And desptite that I think web-based neuroimaging analysis is very cool, the MRI toolbox is not very flexible yet (you have to manually enter the measurements). Therefore, I set out to find out what this T2(*)-fitting actually does. Basically it fits the following signal decay function: $$S(TE) = S_0 e^{-TE/T2(*)}$$   Continue reading