openCV【实践系列】5——使用OpenCV进行Delaunay三角剖分
本帖最后由 fantomas 于 2018-10-19 10:21 编辑使用OpenCV进行Delaunay三角剖分和Voronoi图
图1：左图奥巴马总统使用dlib检测到标志点，中间Delaunay三角剖分的标志点，右图：相应的Voronoi图 在面部标志的众多应用中，首先发现面部地标的三角测量（参见图1），并且这些三角形被扭曲来一些有趣的事情。这篇文章将帮助我们理解Delaunay三角剖分和Voronoi图（又名Voronoi tesselation，Voronoi分解，Voronoi分区和Dirichlet曲面细分），并帮助我们发现OpenCV中几乎没有记录的函数。
[*]什么是Delaunay三角剖分？
图2：有益于小角度的Delaunay三角剖分 给定平面中的一组点，三角剖分指的是将平面细分为三角形，其中点为顶点。在图1中，我们在左图上看到一组标志，在中间图像中看到三角剖分。一组点可以有很多可能的三角形，但Delaunay三角剖分因其具有一些不错的属性而脱颖而出。在Delaunay三角剖分中，选择三角形使得任何三角形的外接圆内都没有点。图2显示Delaunay三角剖分的4个点A，B，C和D。在上图中，为了使三角剖分成为有效的Delaunay三角剖分，点C应该在三角形ABD的外接圆之外，并且点A应该在三角形BCD的外接圆之外。Delaunay三角剖分的一个有趣特性是它不支持“瘦”三角形（即具有一个大角度的三角形）。图2显示了移动点时三角剖分如何变化以选择“胖”三角形。在图2左图中，点B和D的x坐标在x=1.5处，而在图2右图中，它们向右移动到x=1.75。在图2左图中，角度ABC和ABD很大，并且Delaunay三角剖分在B和D之间产生边缘，将两个大角度分成更小的角度ABD，ADB，CDB和CBD。另一方面，在图2右图中，角度BCD太大，并且Delaunay三角剖分产生边缘AC以划分大角度。有许多算法可以找到一组点的Delaunay三角剖分。最明显（但不是最有效）的是从任何三角测量剖分，并检查任何三角形的外接圆是否包含另一个点。如果是这样，翻转边缘（如图2所示）并继续，直到没有三角形的外接圆包含一个点。任何关于Delaunay三角剖分的讨论都必须包括Voronoi图，因为一组点的Voronoi图是其Delaunay三角剖分的数学双重图。
[*]什么是Voronoi图？
给定平面中的一组点，Voronoi图划分空间，使得边界线与相邻点等距。图3显示了用黑点计算的Voronoi图的示例。你会注意到每条边界线都穿过两点的中心。如果连接相邻Voronoi区域中的点，则会得到Delaunay三角剖分！Delaunay三角剖分和Voronoi图以多种方式相关。
[*]OpenCV中的Delaunay三角剖分和Voronoi图
给定一组点，你可以使用Subdiv2D类计算Delaunay三角剖分或Voronoi图。这是操作步骤：1.收集矢量中的所有点Python：points = []
# 添加每一组
points.append((x, y))2. 使用矩形（rect）定义要分区的空间。如果在上一步中定义的点是在图像上定义的，则此矩形可以是（0,0，width，height）。否则，您可以选择一个包含点的矩形。img = cv2.imread("image.jpg");
size = img.shape
rect = (0, 0, size, size)3.使用上一步中获得的矩形创建Subdiv2D的实例subdiv= cv2.Subdiv2D(rect);4.使用bdiv.insert（point）将点插入subdiv。上面的视频显示了三角测量的动画，因为我们将细分添加到细分。5.使用bdiv.getTriangleList获取Delaunay三角形列表。6.使用bdiv.getVoronoiFacetList获取Voronoi方面的列表。
[*]Delaunay三角剖分和Voronoi图的OpenCV示例
这是一个完整的工作示例。我已经从OpenCV附带的示例中复制了一些代码，并对其进行了简化和修改，以满足我们的目的。OpenCV附带的python示例使用旧的（和丑陋的）接口，所以我从头开始编写它。此代码假定图像存储在image.jpg中，点存储在points.txt中。points.txt的每一行包含由空格分隔的点的x和y坐标。import cv2
import numpy as np
import random
# 检查一个点是否在矩形内
def rect_contains(rect, point) :
if point < rect :
return False
elif point < rect :
return False
elif point > rect :
return False
elif point > rect :
return False
return True
# 绘制一个点
def draw_point(img, p, color ) :
cv2.circle( img, p, 2, color, cv2.cv.CV_FILLED, cv2.CV_AA, 0 )
# 绘制 delaunay 三角剖分
def draw_delaunay(img, subdiv, delaunay_color ) :
triangleList = subdiv.getTriangleList();
size = img.shape
r = (0, 0, size, size)
for t in triangleList :
pt1 = (t, t)
pt2 = (t, t)
pt3 = (t, t)
if rect_contains(r, pt1) and rect_contains(r, pt2) and rect_contains(r, pt3) :
cv2.line(img, pt1, pt2, delaunay_color, 1, cv2.CV_AA, 0)
cv2.line(img, pt2, pt3, delaunay_color, 1, cv2.CV_AA, 0)
cv2.line(img, pt3, pt1, delaunay_color, 1, cv2.CV_AA, 0)
# 绘制 voronoi 图
def draw_voronoi(img, subdiv) :
( facets, centers) = subdiv.getVoronoiFacetList([])
for i in xrange(0,len(facets)) :
ifacet_arr = []
for f in facets :
ifacet_arr.append(f)
ifacet = np.array(ifacet_arr, np.int)
color = (random.randint(0, 255), random.randint(0, 255), random.randint(0, 255))
cv2.fillConvexPoly(img, ifacet, color, cv2.CV_AA, 0);
ifacets = np.array()
cv2.polylines(img, ifacets, True, (0, 0, 0), 1, cv2.CV_AA, 0)
cv2.circle(img, (centers, centers), 3, (0, 0, 0), cv2.cv.CV_FILLED, cv2.CV_AA, 0)
if __name__ == '__main__':
win_delaunay = "Delaunay Triangulation"
win_voronoi = "Voronoi Diagram"
# 当绘制三角形剖分时打开动画画板
animate = True
# 定义绘制颜色
delaunay_color = (255,255,255)
points_color = (0, 0, 255)
img = cv2.imread("image.jpg");
img_orig = img.copy();
# 创建用于Subdiv2D 的矩形
size = img.shape
rect = (0, 0, size, size)
# 创建Subdiv2D 实例
subdiv = cv2.Subdiv2D(rect);
points = [];
# 从 text 文件中读取点
with open("points.txt") as file :
for line in file :
x, y = line.split()
points.append((int(x), int(y)))
# 将点依次插入subdiv中
for p in points :
subdiv.insert(p)
# 展示动画画板
if animate :
img_copy = img_orig.copy()
draw_delaunay( img_copy, subdiv, (255, 255, 255) );
cv2.imshow(win_delaunay, img_copy)
cv2.waitKey(100)
# 绘制delaunay 三角剖分
draw_delaunay( img, subdiv, (255, 255, 255) );
for p in points :
draw_point(img, p, (0,0,255))
# 为Voronoi 图分配空间
img_voronoi = np.zeros(img.shape, dtype = img.dtype)
# 绘制 Voronoi 图
draw_voronoi(img_voronoi,subdiv)
cv2.imshow(win_delaunay,img)
cv2.imshow(win_voronoi,img_voronoi)
cv2.waitKey(0)得到的结果就是和图1的中间图和右图一样，如果想要看Dlaunay三角剖分的动态过程可以访问这里，或者自己运行上述的代码网址：https://www.learnopencv.com/delaunay-triangulation-and-voronoi-diagram-using-opencv-c-python/ 藏着很多宝贝啊
这个内容
现在与GAN结合
火爆的一塌糊涂 好久咩有恩关系了 opencv是基础 学会opencv 图像处理会了一大半 再加上deep learning 就如虎添翼
页:
[1]