Multishell

Open bcmcpher opened this issue 5 years ago • 0 comments

A working estimate for the LiFE model with multishell data. This pull request is to share a working development version.

This branch is still being actively tested. Forward model may not work. Single shell data will have problems being passed with the functions as they are. Change between tensor forward model and kurtosis forward model by commenting line 132/133 in feConnectomeBuildDictionaries.m

This is for testing changes required for multishell analysis / forward model.

% this will work for a test fg = fgRead('cc_crop.tck'); nv = [ 0.0017 0.0003 0.0003; 0.0014 0.0002 0.0002 ]; k = [ 88.2725 -129.2530 -159.8982 -24.8088 32.5064 -0.3677 3.9417 -5.1648 2.4475 -32.0594 -24.1665 -45.5590 -12.9820 9.5266 -12.6976 ]; fe = feConnectomeInit('dwi.nii.gz', fg, [], [], [], 't1.nii.gz', 720, nv, k); fit = feFitModel(feGet(fe, 'model'), feGet(fe, 'dsigdemeaned'), 'bbnnls', 500, 'preconditioner');

THIS IS A DEV BRANCH - NOT FULLY FUNCTIONAL

Nov 16 '20 23:11 bcmcpher