f_th = np.zeros((a, b))