update RK calculation to be the same as matlab

NYU-DiffusionMRI · Mar 22, 2024 · a0c306a · a0c306a
1 parent d2fcb9c
commit a0c306a
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 1 deletion.
diff --git a/constant/dirs256.mat b/constant/dirs256.mat
diff --git a/lib/tensor.py b/lib/tensor.py
@@ -213,7 +213,11 @@ def extract_parameters(self, dt, b, mask, extract_dti, extract_dki, fit_w=False)
             parameters['trace'] = trace
 
         if extract_dki:
-            dirs = np.array(self.fibonacci_sphere(256, True))
+            #dirs = np.array(self.fibonacci_sphere(256, True))
+            dwd = os.path.dirname(os.path.dirname(os.path.abspath(__file__)))
+            mat = sio.loadmat(os.path.join(dwd,'constant','dirs256.mat'))
+            dirs = mat['dirs']
+
 
             if fit_w:
                 akc = self.w_kurtosis_coeff(dt, dirs)