OpenCV计算机图像处理 —— 凸性缺陷 + 点多边形测试 + 形状匹配 + 轮廓分层与cv.findContours()

OpenCV计算机图像处理 —— 凸性缺陷 + 点多边形测试 + 形状匹配 + 轮廓分层与cv.findContours()

1. 凸性缺陷

一般来说凸曲线都是凸出或平坦的曲线，如果在内部凸出了（凹进去了）我们就称其为凸性缺陷，OpenCV提供了一个方法cv.convexityDefects()

这个函数返回一个数组，其中每行包含这些值-【起点，终点，最远点，到最远点的近似距离】，我们可以用图像把它形象化，我们画一条连接起点和终点的线，然后在最远处画一个圆

import cv2 as cv import numpy as np img = cv.imread(r'E:\image\test14.png') img_gray = cv.cvtColor(img, cv.COLOR_BGR2GRAY) ret, thresh = cv.threshold(img_gray, 127, 255, 0) contours, hierarchy = cv.findContours(thresh, 2, 1) cnt = contours[0] hull = cv.convexHull(cnt, returnPoints=False) defects = cv.convexityDefects(cnt, hull) for i in range(defects.shape[0]): s, e, f, d = defects[i, 0] start = tuple(cnt[s][0]) end = tuple(cnt[e][0]) far = tuple(cnt[f][0]) cv.line(img, start, end, [0, 255, 0], 2) cv.circle(img, far, 5, [0, 0, 255], -1) cv.imshow('img', img) cv.waitKey(0) cv.destroyAllWindows()

在以前轮廓的内容中我们提到了了一个函数cv.convexHull()，当这个函数的属性returnPoints = False时，我们就可以找到凸性缺陷

代码解释：又到了激动人心的解析代码的时间啦！第4-5行即是老生常谈的读取图片并转化为灰度图，第6行通过threshold()函数设置阈值，第7行通过findContours()方法找到轮廓，将第一个轮廓拿出来后用convexHull()函数检查轮廓的凸度缺陷，然后使用convexityDefects()函数得到这些凸度缺陷对应点的相关信息，包括起点终点等

2. 点多边形测试

点多边形测试中涉及了一个新的函数：cv.pointPolygonTest()，顾名思义它是一个测试函数，第一个参数是轮廓，第二个参数就是我们的测试点了，第三个参数指的是measureDist，如果其取值为True，这个函数会找到有符号的距离，反之它只会反映出这个点是在轮廓线内部、上面还是外部（分别返回1、-1和0）

import cv2 as cv import numpy as np img = cv.imread(r'E:\image\test12.png') img_gray = cv.cvtColor(img, cv.COLOR_BGR2GRAY) ret, thresh = cv.threshold(img_gray, 127, 255, 0) contours, hierarchy = cv.findContours(thresh, 2, 1) cnt = contours[0] img02 = cv.drawContours(img, [cnt], 0, (0, 0,255), 3) dist = cv.pointPolygonTest(cnt, (20, 20), True) dist2 = cv.pointPolygonTest(cnt, (192.5, 156), True) print("point:(20,20) = " + str(dist)) print("point:(192.5,156) = " + str(dist2)) cv.imshow('test', img02) cv.waitKey(0) cv.destroyWindow()

3. 形状匹配

OpenCV中的形状匹配与一个函数有关：cv.matchShapes()，这个函数使我们能够比较两个形状（或两个轮廓），并返回一个显示相似性的度量，这个值越低说明匹配越好，在实质上它是根据矩值计算出来的

值得注意的是，即使是图像旋转也不会对结果产生很大的影响

import cv2 as cv import numpy as np from matplotlib import pyplot as pyt img1 = cv.imread(r'E:/image/star01.png', 0) img2 = cv.imread(r'E:/image/star02.png', 0) img3 = cv.imread(r'E:/image/rectangle01.png', 0) ret, thresh = cv.threshold(img1, 127, 255, 0) ret, thresh2 = cv.threshold(img2, 127, 255, 0) ret, thresh3 = cv.threshold(img3, 127, 255, 0) contours, hierarchy = cv.findContours(thresh, 2, 1) cnt1 = contours[0] contours, hierarchy = cv.findContours(thresh2, 2, 1) cnt2 = contours[0] contours, hierarchy = cv.findContours(thresh3, 2, 1) cnt3 = contours[0] ret = cv.matchShapes(cnt1, cnt2, 1, 0.0) ret2 = cv.matchShapes(cnt1, cnt3, 1, 0.0) print("星图1与星图2的匹配度量 = " + str(ret)) print("星图1与方图1的匹配度量 = " + str(ret2)) pyt.subplot(1, 3, 1), pyt.imshow(img1, cmap="gray") pyt.title("star_img01"), pyt.xticks([]), pyt.yticks([]) pyt.subplot(1, 3, 2), pyt.imshow(img2, cmap="gray") pyt.title("star_img02"), pyt.xticks([]), pyt.yticks([]) pyt.subplot(1, 3, 3), pyt.imshow(img3, cmap="gray") pyt.title("rectangle_img01"), pyt.xticks([]), pyt.yticks([]) pyt.show()

4. 轮廓分层与cv.findContours()

4.1 轮廓分层与cv.findcontours()的关系

在前面的文章我们聊了怎么使用cv.findcontours()函数找到轮廓，并且提到了它的第三个参数轮廓检索模式，对于这个参数我们通常使用cv.RETR_LIST或cv.RETR_TREE就能得到很好的效果，而且对应这个函数的三个输出我们也很有必要了解一下，它的输出包括三个数组，包括图像（image）、轮廓（contours）和hierarchy，它其实就是我们在这一节要聊的”轮廓层次结构”

在某些情况下，某些形状位于其他形状中，就像嵌套的图形一样，在这种情况下，我们把外部的称为父类，把内部的称为子类，这样，图像中的轮廓就有了一定的相互关系。我们可以指定一个轮廓是如何相互连接的，比如，它是另一个轮廓的子轮廓，还是父轮廓等等。这种关系的表示称为层次结构

所以每个轮廓都有它自己的信息关于它是什么层次，谁是它的孩子，谁是它的父母等等，OpenCV将它表示为一个包含四个值的数组: [Next, Previous, First_Child, Parent]，这就是我们在findcontours()函数中得到的第三个参数

4.2 轮廓检索模式（四种参数）

RETR_LIST

这是四个标志中最简单的一个，它只是检索所有的轮廓，但不创建任何亲子关系，在这个规则下，父轮廓和子轮廓是平等的，他们只是轮廓，并且都属于同一层级，如果我们在开发中没有使用任何层次结构特性，它就是我们最佳的选择

RETR_EXTERNAL

如果使用此标志，它只返回极端外部标志，所有孩子的轮廓都被留下了，更形象地来说，根据这项规则，每个家庭只有长子得到关注，它不关心家庭的其他成员，所以当我们只关心外部轮廓时可以使用它

RETR_CCOMP

此标志检索所有轮廓并将其排列为2级层次结构，物体的外部轮廓(即物体的边界)放在层次结构-1中，对象内部孔洞的轮廓(如果有)放在层次结构-2中，如果其中有任何对象，则其轮廓仅在层次结构1中重新放置，以及它在层级2中的漏洞等等

RETR_TREE

它就是最完美的一个家伙了，它检索所有的轮廓并创建一个完整的家族层次结构列表，当我们需要检索所有轮廓时会使用它

接下来通过代码看看吧，实践永远是最好的老师

import cv2 as cv import numpy as np from matplotlib import pyplot as pyt img = cv.imread(r'E:/image/the_first.png') img_gray = cv.cvtColor(img, cv.COLOR_BGR2GRAY) ret, thresh = cv.threshold(img_gray, 127, 255, 0) contours, hierarchy = cv.findContours(thresh, cv.RETR_LIST, cv.CHAIN_APPROX_NONE) cnt = contours[0] contours, hierarchy = cv.findContours(thresh, cv.RETR_EXTERNAL, cv.CHAIN_APPROX_NONE) cnt2 = contours[0] contours, hierarchy = cv.findContours(thresh, cv.RETR_CCOMP, cv.CHAIN_APPROX_NONE) cnt3 = contours[0] contours, hierarchy = cv.findContours(thresh, cv.RETR_TREE, cv.CHAIN_APPROX_NONE) cnt4 = contours[0] img1 = cv.drawContours(thresh, [cnt], 0, (255, 0, 0), 2) img2 = cv.drawContours(thresh, [cnt2], 0, (0, 255, 0), 2) img3 = cv.drawContours(thresh, [cnt3], 0, (0, 0, 255), 2) img4 = cv.drawContours(thresh, [cnt4], 0, (0, 0, 0), 2) pyt.subplot(2, 2, 1), pyt.imshow(img1, cmap="gray") pyt.title("RETR_LIST"), pyt.xticks([]), pyt.yticks([]) pyt.subplot(2, 2, 2), pyt.imshow(img2, cmap="gray") pyt.title("RETR_EXTERNAL"), pyt.xticks([]), pyt.yticks([]) pyt.subplot(2, 2, 3), pyt.imshow(img3, cmap="gray") pyt.title("RETR_CCOMP"), pyt.xticks([]), pyt.yticks([]) pyt.subplot(2, 2, 4), pyt.imshow(img3, cmap="gray") pyt.title("RETR_TREE"), pyt.xticks([]), pyt.yticks([]) pyt.show()

（注：文章内容参考OpenCV4.1中文官方文档）

如果文章对您有所帮助，记得一键三连支持一下哦

AI OpenCV 图像处理 机器学习 机器视觉

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。