问题描述
在为Delaunay三角剖分开发代码之后,我列出了节点坐标,并发现有少量重复节点。
所以我避免了重复的节点:
distort = tri.Triangulation(mesh_x,mesh_y) #triangulation
#making list of nodes coordinates
data = np.array([mesh_x,mesh_y])
data = np.transpose(data)
#sorting avoids duplicated nodes
unique_data = np.unique(data,axis = 0)
现在,连接节点的边线存在问题。删除重复的节点后,将重新分配边缘,而不是使用平滑网格,而不是像这样:image
如何检查顶点匹配并替换边缘的重复索引并获得平滑的网格? (例如image)
解决方法
每次删除一个节点时,都需要向所删除的节点索引较大的三角形顶点的所有索引中减去一个。