I get numpy warnings on blank float arrays if not use masked arrays.

    arr = np.empty((3,3), dtype=np.float64)
    arr.fill(100)
    threshold_isodata(arr)

    RuntimeWarning: invalid value encountered in divide
      l = np.cumsum(pmf * binnums, dtype=float) / cpmfl