Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Multishell #97

Open
wants to merge 19 commits into
base: master
Choose a base branch
from
Open

Multishell #97

wants to merge 19 commits into from

Conversation

bcmcpher
Copy link
Collaborator

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

Brent McPherson added 19 commits November 11, 2020 17:34
…rtosis coefficient (akc); modifications for multishell data - bookmark, not for release
…ensor model still - this is still locked on even when the kurtosis model isnt valid. Split one long print statement into 3 lines in feConnectomeEncoding. comments changed in feGet
…parently the original orientation of the data for linear indexing
…ned data, updated comments to match / describe
…he functions need to be updated to correctly parse inputs for these options - but that has to happen when the options work in the first place
…ed feConnectomeInit to pass warnings if a mismatch in file orientations is detected
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

1 participant