scaling: ensure mad is at least 1.0

nanoporetech · Feb 4, 2022 · 6e91a9d · 6e91a9d
1 parent 6f23467
commit 6e91a9d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/bonito/fast5.py b/bonito/fast5.py
@@ -82,7 +82,7 @@ def __init__(self, read, filename, meta=False):
 
         if len(scaled) > 8000:
             med, mad = med_mad(scaled)
-            self.signal = (scaled - med) / mad
+            self.signal = (scaled - med) / max(1.0, mad)
         else:
             self.signal = norm_by_noisiest_section(scaled)