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)