BF: KFA is zero for small MKT values (MKT < tol, tol=1e-8)

dipy · Jun 14, 2021 · 1094a97 · 1094a97
1 parent f6f36ed
commit 1094a97
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/dipy/reconst/dki.py b/dipy/reconst/dki.py
@@ -1457,7 +1457,8 @@ def kurtosis_fractional_anisotropy(dki_params):
     # Compute KFA
     KFA = np.zeros(A.shape)
     cond1 = B > 0  # Avoiding Singularity (if B = 0, KFA = 0)
-    cond2 = W > 0  # Avoiding overestimating KFA for MKT = 0
+    # Avoiding overestimating KFA for small MKT values (KFA=0, MKT < tol)
+    cond2 = W > 1e-8
     cond = np.logical_and(cond1, cond2)
     KFA[cond] = np.sqrt(A[cond]/B[cond])