试图从 Python 中的很长的一行中删除过时的点

Question

问题描述：

我有大量代表图形边缘的二维点（> 150k）。 点坐标都是integer（点之间的最小距离=1）。 绘制边缘点如下所示：

有的地方看似线条相接，实则不然。 我想将这些点导入到 CAD 程序中并绘制边缘线。 然而，为此，我必须确保每个点恰好有 2 个距离为 1 的邻居。如果一个点有超过 2 个邻居，则需要删除点，直到只有 2 个。删除哪些点并不重要正是因为小的几何偏差不相关。

有没有办法在 Python 中有效地做到这一点？ 不幸的是，我缺乏经验。

方法：

这是我的第一个（非常低效的）方法：

edges_points是输入数组。 它是一个元组列表，每个元组包含一个点的 x 和 y 坐标：

In [20]: len(edges_points)
Out[20]: 156350

In [21]: edges_points[0]
Out[21]: (138, 140)

这是 Python 脚本（edges_points 已经定义）：

def remove_if_more_than_two_neighbors(edges_points,x,y):
    check = 0
    start_over = 0
    for i in [-1,0,1]: #check for all 8 possible neighbors
        for j in [-1,0,1]: 
            if p+i >= len(edges_points) or p+j>= len(edges_points) or p+i <0 or p+j <0: #prevent out-of-range indexing
                continue 
            if (x+i,y+j) in edges_points: #check if there is a neighboring point in the input list
                check = check+1
                if check > 3: #if there are more than 2 neighbors, remove the unwanted point and start over
                    edges_points.remove((x+i,y+j))
                    start_over = 1
                    return [edges_points, start_over]

    return [edges_points, start_over]

for p in range(len(edges_points)):
    stop = 0
    while(stop == 0):
        check = 0
        x = edges_points[p][0]
        y = edges_points[p][1]
        func_return = remove_if_more_than_two_neighbors(edges_points,x,y)
        edges_points=func_return[0]
        start_over=func_return[1]
        if start_over == 0:           
            stop = 1

脚本似乎在删除点，但速度太慢，几乎没有意义。 我意识到关于此任务还有一些其他问题（例如删除断开连接的点或一般的点排序）但我稍后会解决这些问题。

感谢您的帮助！

Answer 1

如果坐标范围在 [0, 10000] x [0, 10000] 或更小，则使用该大小的图像（每像素 1 位或每像素 1 字节）是有意义的。

首先清除图像并绘制每个点。 然后在第二次扫描这些点时，您可以快速测试最近的邻居。

---

如果这对您来说太大了，您可以求助于分箱：将图像视为带有 16 x 16（比方说）瓷砖的瓷砖，并且只表示包含点的瓷砖。 您将需要 hash map 来建立图块坐标和图块图像之间的连接。

设置像素时，首先检查图块是否存在并在必要时创建它； 然后以正确的偏移量写入磁贴。 稍后测试一个像素，如果瓦片不存在，则该点不存在； 否则，您会以正确的偏移量阅读。