def show_different_ds(bvals,bvecs):
global CNT
r=fvtk.ren()
r.SetBackground(1.,1.,1.)
Signals=[]
for d in [0.0005,0.0010,0.0015,0.0020,0.0025,0.0030,0.0035,0.0040]:
S1,sticks1=simulations_dipy(bvals,bvecs,d=d,S0=100,angles=[(0,0),(45,0),(90,90)],fractions=[100,0,0],snr=None)
Signals.append(S1[1:]/S1[0])
show_signals(r,Signals,bvecs)
CNT=0
draw_needles(r,sticks1,100,2,off=np.array([0,0,0]))
Signals=[]
for d in [0.0005,0.0010,0.0015,0.0020,0.0025,0.0030,0.0035,0.0040]:
S1,sticks1=simulations_dipy(bvals,bvecs,d=d,S0=100,angles=[(0,0),(45,0),(90,90)],fractions=[50,50,0],snr=None)
Signals.append(S1[1:]/S1[0])
show_signals(r,Signals,bvecs,-200)
CNT=0
Signals=[]
for d in [0.0005,0.0010,0.0015,0.0020,0.0025,0.0030,0.0035,0.0040]:
S1,sticks1=simulations_dipy(bvals,bvecs,d=d,S0=100,angles=[(0,0),(45,0),(90,90)],fractions=[33,33,33],snr=None)
Signals.append(S1[1:]/S1[0])
show_signals(r,Signals,bvecs,-400)
"""
Signals=[]
for d in [0.0005,0.0010,0.0015,0.0025,0.0030,0.0035,0.0040]:
S1,sticks1=simulations_dipy(bvals,bvecs,d=d,S0=100,angles=[(0,0),(90,0),(90,90)],fractions=[20,50,0],snr=None)
Signals.append(S1[1:]/S1[0])
show_signals(r,Signals,bvecs,-600)
Signals=[]
for d in [0.0005,0.0010,0.0015,0.0025,0.0030,0.0035,0.0040]:
S1,sticks1=simulations_dipy(bvals,bvecs,d=d,S0=100,angles=[(0,0),(90,0),(90,90)],fractions=[10,50,0],snr=None)
Signals.append(S1[1:]/S1[0])
show_signals(r,Signals,bvecs,-800)
Signals=[]
for d in [0.0005,0.0010,0.0015,0.0025,0.0030,0.0035,0.0040]:
S1,sticks1=simulations_dipy(bvals,bvecs,d=d,S0=100,angles=[(0,0),(90,0),(90,90)],fractions=[0,50,0],snr=None)
Signals.append(S1[1:]/S1[0])
show_signals(r,Signals,bvecs,-1000)
"""
fvtk.show(r)
Bvecs=np.concatenate([bvecs[1:],-bvecs[1:]])
for (i,S) in enumerate(Signals):
X=np.dot(np.diag(np.concatenate([50*S,50*S])),Bvecs)
mX=np.mean(np.sqrt(X[:,0]**2+X[:,1]**2+X[:,2]**2))
print CNT,mX,np.var(S)
CNT=CNT+1
fvtk.add(r,fvtk.point(X+np.array([(i+1)*200,offset,0]),fvtk.green,1,2,6,6))
#fvtk.show(r)
if __name__=='__main__':
np.set_printoptions(4,suppress=True)
data,bvals,bvecs=get_data(name='118_32',par=0)
Bvecs=np.concatenate([bvecs[1:],-bvecs[1:]])
S0,sticks0=simulations_dipy(bvals,bvecs,d=0.0015,S0=100,angles=[(0,0),(90,0),(90,90)],fractions=[100,0,0],snr=None)
X0=np.dot(np.diag(np.concatenate([S0[1:],S0[1:]])),Bvecs)
#CBK,STK,FRA=needlebook(bvals,bvecs,d=0.0015,subdiv=4,fractions=[.8])
#res0=single_stick_simulations(200,bvals,bvecs,d=0.0015,fractions=[80,0,0],snr=None)
#CBK,STK,FRA=needlebook(bvals,bvecs,d=0.0015,subdiv=3,fractions=[.2])
#res1=single_stick_simulations(200,bvals,bvecs,d=0.0015,fractions=[80,0,0],snr=None)
#res0=single_stick_simulations(1000,bvals,bvecs,d=0.0015,fractions=[100,0,0],snr=None)
#res1=single_stick_simulations(1000,bvals,bvecs,d=0.0015,fractions=[80,0,0],snr=None)
#res2=single_stick_simulations(1000,bvals,bvecs,d=0.0015,fractions=[60,0,0],snr=None)
#res3=single_stick_simulations(1000,bvals,bvecs,d=0.0015,fractions=[40,0,0],snr=None)
#res4=single_stick_simulations(1000,bvals,bvecs,d=0.0015,fractions=[20,0,0],snr=None)
# values=psi(S0[1:]/S0[0],CBK,STK,FRA)
