from dipy.viz import fvtk
from dipy.data import fetch_isbi2013_2shell, read_isbi2013_2shell, get_sphere
from dipy.core.gradients import gradient_table
"""
Download and read the data for this tutorial.
fetch_isbi2013_2shell() provides data from the ISBI HARDI contest 2013 acquired
for two shells at b-values 1500 and 2500.
The six parameters of these two functions define the ROI where to reconstruct
the data. They respectively correspond to (xmin,xmax,ymin,ymax,zmin,zmax)
with x, y, z and the three axis defining the spatial positions of the voxels.
"""
fetch_isbi2013_2shell()
img, gtab = read_isbi2013_2shell()
data = img.get_data()
data_small = data[10:40, 22, 10:40]
print('data.shape (%d, %d, %d, %d)' % data.shape)
"""
data contains the voxel data and gtab contains a GradientTable
object (gradient information e.g. b-values). For example, to show the b-values
it is possible to write print(gtab.bvals).
Instantiate the SHORE Model.
radial_order is the radial order of the SHORE basis.
zeta is the scale factor of the SHORE basis.
from dipy.data import fetch_isbi2013_2shell, read_isbi2013_2shell, get_sphere
from dipy.core.gradients import gradient_table
"""
Download and read the data for this tutorial.
two_shells_voxels() provides data from the ISBI HARDI contest 2013 acquired
for two shells at b-values 1500 and 2500.
The six parameters of these two functions define the ROI where to reconstruct
the data. They respectively correspond to (xmin,xmax,ymin,ymax,zmin,zmax)
with x,y and the three axis defining the spatial positions of the voxels.
"""
fetch_isbi2013_2shell()
img, gtab=read_isbi2013_2shell()
data = img.get_data()
data_small=data[10:40,10:40,25]
print('data.shape (%d, %d, %d, %d)' % data.shape)
"""
data contains the voxel data and gtab contains a GradientTable
object (gradient information e.g. b-values). For example to read the b-values
it is possible to write print(gtab.bvals).
Instantiate the SHORE Model.
radial_order is the radial order of the SHORE basis.
zeta is the scale factor of the SHORE basis.