def refine_positions(npa, atol=1e-05):
    """
    给定一个精度，让靠近1的元素变成零。则使得坐标在0～1之间可以直接比较。
    """
    flatten_npa = npa.flatten()
    for i, val in enumerate(flatten_npa):
        if val < atol or val > 1 - atol:
            flatten_npa[i] = 0
        else:
            flatten_npa[i] = val
    return flatten_npa.reshape(npa.shape)
