最近在处理多b值MRI曲线拟合的数据,每组的数据结构如下
b:[0 20 50 80 100 150 200 400 600 800 1000];
S:[297 283.8 265.2 257.2 256.1 225.8 215.2 169.9 138.5 109.8 101.5];
DWI-MRI成像中b值和信号强度的比值关系如下:
- 单指数模型:S(b)/S(0) = exp(-b*ADC);
- 双指数模型:
- 拉伸指数:
- DKI
- …
有了b值和S值,怎么拟合得到这些参数呢?利用MATLAB的曲线拟合工具箱可以很好的解决这些问题;
举个例子DKI拟合,对应的Matlab函数
ft = fittype( 'exp(-a*x+x^2*a^2*b/6)', 'independent', 'x', 'dependent', 'y' );
opts = fitoptions( ft );
opts.Algorithm = 'Levenberg-Marquardt';
opts.Display = 'Off';