Guard against division by zero (#1099)

DigitalSlideArchive · Feb 16, 2024 · b0b1242 · b0b1242
1 parent 8e6be49
commit b0b1242
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/histomicstk/features/compute_fsd_features.py b/histomicstk/features/compute_fsd_features.py
@@ -186,7 +186,7 @@ def _FSDs(X, Y, K, Intervals):
  fX = np.fft.fft(Curvature).T
  # spectral energy
  fX = fX * fX.conj()
- fX = fX / fX.sum()
+ fX = (fX / fX.sum()) if fX.sum() else 0
  # calculate 'F' values
  for i in range(L - 1):
  F[i] = np.round(

diff --git a/histomicstk/features/compute_haralick_features.py b/histomicstk/features/compute_haralick_features.py
@@ -370,7 +370,7 @@ def compute_haralick_features(im_label, im_intensity, offsets=None,
  pxy_ijr = np.ravel(pxy_ij)
  HXY1 = -np.dot(nGLCMr, np.log2(pxy_ijr + e))
  HXY2 = -np.dot(pxy_ijr, np.log2(pxy_ijr + e))
- ldata.at[r, 'Haralick.IMC1'] = (HXY - HXY1) / max(HX, HY)
+ ldata.at[r, 'Haralick.IMC1'] = ((HXY - HXY1) / max(HX, HY)) if max(HX, HY) else 0
 
  # computes information measures of correlation
  ldata.at[r, 'Haralick.IMC2'] = \