@@ -22,22 +22,22 @@ def ivim_model(b, S0=1, f=0.1, Dstar=0.01, D=0.001):
22
22
return S0 * (f * np .exp (- b * Dstar ) + (1 - f )* np .exp (- b * D ))
23
23
24
24
signals = ivim_model (bvalues )
25
- data = np .array ([signals , signals , signals ])
25
+ data = np .array ([[ signals , signals , signals ], [ signals , signals , signals ] ])
26
26
#print(data)
27
27
signals = data
28
28
29
29
#model = ETP_SRI_LinearFitting(thresholds=[200])
30
30
if kwargs :
31
- results = model .osipi_fit (signals , bvalues , ** kwargs )
31
+ results = model .osipi_fit_full_volume (signals , bvalues , ** kwargs )
32
32
else :
33
- results = model .osipi_fit (signals , bvalues )
33
+ results = model .osipi_fit_full_volume (signals , bvalues )
34
34
print (results )
35
35
#test = model.osipi_simple_bias_and_RMSE_test(SNR=20, bvalues=bvalues, f=0.1, Dstar=0.03, D=0.001, noise_realizations=10)
36
36
37
37
#model1 = ETP_SRI_LinearFitting(thresholds=[200])
38
- # model2 = IAR_LU_biexp(bounds=([0,0,0,0], [1,1,1,1]))
38
+ model2 = IAR_LU_biexp (bounds = ([0 ,0 ,0 ,0 ], [1 ,1 ,1 ,1 ]))
39
39
#model2 = IAR_LU_modified_mix()
40
- model2 = OGC_AmsterdamUMC_biexp ()
40
+ # model2 = OGC_AmsterdamUMC_biexp()
41
41
42
42
#dev_test_run(model1)
43
43
dev_test_run (model2 )
0 commit comments