传统直方图灰度级映射的py代码

    # def cal_new_gray(self, y_mat, cdf_obj, cdf_org, org_idx):
    #     """
    #     计算映射后的新灰度级
    #     """
    #     new = []
    #     dis = []
    #     size = len(cdf_org)
    #     for i in range(size):
    #         t = np.abs(cdf_org[i] - cdf_obj)
    #         index = t.argmin()
    #         temp = index - i
    #         new.append(index)
    #         dis.append(temp)
    #     new_gray = new
    #     assert isinstance(y_mat, np.ndarray)  # 判断类型
    #     h, w = y_mat.shape
    #     new_y = np.zeros_like(y_mat)
    #     new_gray = np.array(new_gray)
    #     idx = np.array(org_idx)
    #     for i in range(h):
    #         for j in range(w):
    #             new_y[i, j] = new_gray[int(np.where(idx == y_mat[i, j])[0])]
    #     return new_y