Merge pull request #1228 from nipreps/fix/cc-snr-usemad

ENH: Use MAD for robust estimation of sigma in the CC mask
nipreps · Apr 3, 2024 · aff0d63 · aff0d63
2 parents 3fa5bd9 + b176260
commit aff0d63
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/mriqc/qc/diffusion.py b/mriqc/qc/diffusion.py
@@ -92,6 +92,7 @@
 from __future__ import annotations
 
 import numpy as np
+from statsmodels.robust.scale import mad
 
 
 def noise_b0(
@@ -231,7 +232,7 @@ def cc_snr(
     """
 
     cc_mask = cc_mask > 0  # Ensure it's a boolean mask
-    std_signal = in_b0[cc_mask].std()
+    std_signal = mad(in_b0[cc_mask])
 
     cc_snr_estimates = {}