Harden haralick computation. (#1122)

Should resolve #1116.
DigitalSlideArchive · Aug 1, 2024 · 4e2e2eb · 4e2e2eb
1 parent d0a8fd9
commit 4e2e2eb
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/histomicstk/features/compute_haralick_features.py b/histomicstk/features/compute_haralick_features.py
@@ -364,7 +364,7 @@ def compute_haralick_features(im_label, im_intensity, offsets=None,
             IMC1 = (HXY - HXY1) / max(HX, HY)
 
             # computes information measures of correlation
-            IMC2 = np.sqrt(1 - np.exp(-2.0 * (HXY2 - HXY)))
+            IMC2 = np.sqrt(np.maximum(0, 1 - np.exp(-2.0 * (HXY2 - HXY))))
 
             features_per_offset[r] = [
                 ASM,