def METHOD_NAME(points, new_width, box_size, width, buf_size, ndim):
    """
    Creates a grid where the voxels are 1 if there is a point there and 0 if
    not.
    """
    d = (new_width,) * ndim
    arr = np.zeros(d)
    img_points = points / (box_size / width)  # points in units of grid spacing
    for pt in img_points:
        shifted_pt = tuple(int(round(pt[i])) for i in range(ndim))
        arr[shifted_pt] = 1
    return arr