def METHOD_NAME(zz_lst, data, z_spaccing):
    data_interpol = np.copy(data)
    for z_hole_start, z_hole_end in list(_list2range(zz_lst)):
        z_ref_start, z_ref_end = z_hole_start - 1, z_hole_end
        slice_ref_start, slice_ref_end = data[:, :, z_ref_start], data[:, :, z_ref_end]
        hole_cur_lst = list(range(z_hole_start, z_hole_end))
        lenght_hole = len(hole_cur_lst) + 1
        phys_lenght_hole = lenght_hole * z_spaccing
        denom_interpolation = (lenght_hole + 1)
        if phys_lenght_hole < 10:
            logger.warning('Filling a hole in the segmentation around z_slice #:' + str(z_ref_start))
            for idx_z, z_hole_cur in enumerate(hole_cur_lst):
                num_interpolation = (lenght_hole - idx_z - 1) * slice_ref_start  # Contribution of the bottom ref slice
                num_interpolation += (idx_z + 1) * slice_ref_end  # Contribution of the top ref slice
                slice_interpolation = num_interpolation * 1. / denom_interpolation
                slice_interpolation = (slice_interpolation > 0).astype(int)
                data_interpol[:, :, z_hole_cur] = slice_interpolation
    return data_interpol