def METHOD_NAME(original_points: npt.NDArray, voxel_size: float):
    """
    Voxel Point Sampling function.
    This function sample N-dimensional points with voxel grid.
    Points in a same voxel grid will be merged by mean operation for sampling.
    Parameters
    ----------
    original_points :  (M, N) N-dimensional points for sampling.
                        The number of points is M.
    voxel_size : voxel grid size
    Returns
    -------
    sampled points (M', N)
    """
    voxel_dict = defaultdict(list)
    for i in range(original_points.shape[0]):
        xyz = original_points[i, :]
        xyz_index = tuple(xyz // voxel_size)
        voxel_dict[xyz_index].append(xyz)
    points = np.vstack([np.mean(v, axis=0) for v in voxel_dict.values()])
    return points