sh demo for glsl

sh_demo.py

import numpy as np 
from dipy.sims.voxel import multi_tensor_odf 
from dipy.data import get_sphere
from dipy.reconst.shm import sf_to_sh


mevals = np.array([[0.0015, 0.0003, 0.0003],
                   [0.0015, 0.0003, 0.0003]])
angles = [(0, 0), (60, 0)]
fractions = [50, 50]
sphere = get_sphere('repulsion724')
sphere = sphere.subdivide(2)

# Get ODF
odf = multi_tensor_odf(sphere.vertices, mevals, angles, fractions)

#Get spherical harmonics coefficient
odf_sh = sf_to_sh(odf, sphere)
print(odf_sh) # contains 15 coefficients to pass to the shader


# to get some informations about the basis k01, K02, etc..... look at 
# - https://github.com/dipy/dipy/blob/master/dipy/reconst/shm.py#L335
# - https://github.com/dipy/dipy/blob/master/dipy/reconst/shm.py#L233
# - https://github.com/dipy/dipy/blob/master/dipy/reconst/shm.py#L857