def dsi_voxels():
fimg, fbvals, fbvecs = get_fnames('small_101D')
bvals = np.loadtxt(fbvals)
bvecs = np.loadtxt(fbvecs).T
data, _ = load_nifti(fimg)
gtab = gradient_table(bvals, bvecs)
return data, gtab
def dsi_deconv_voxels():
from dipy.sims.voxel import sticks_and_ball
gtab = gradient_table(np.loadtxt(get_fnames('dsi515btable')))
data = np.zeros((2, 2, 2, 515))
for ix in range(2):
for iy in range(2):
for iz in range(2):
data[ix, iy, iz], _ = sticks_and_ball(gtab,
d=0.0015,
S0=1.,
angles=[(0, 0), (90, 0)],
fractions=[50, 50],
snr=None)
return data, gtab
def two_cingulum_bundles():
fname = get_fnames('cb_2')
res = np.load(fname)
cb1 = relist_streamlines(res['points'], res['offsets'])
cb2 = relist_streamlines(res['points2'], res['offsets2'])
return cb1, cb2