python + opencv图像处理——圆检测

本文详细介绍了霍夫圆检测的基本原理,包括从平面坐标到极坐标的转换过程,以及如何利用霍夫变换实现圆的检测。文章还强调了霍夫圆检测对噪声的敏感性,并提出使用中值滤波作为预处理步骤。最后,通过OpenCV的实现展示了如何在图像中检测圆,包括边缘保留滤波、圆心检测和半径计算。

霍夫圆检测原理
1、从平面坐标到极坐标转换三个参数C(x0,y0,r),其中x0,y0是圆心;
2、假设平面坐标的任意一个圆上的点,转换到极坐标中,C(x0,y0,r)处有最大值,霍夫变换正是利用这个原理实现圆检测;
3、因为霍夫圆检测对噪声比较敏感,所以首先要对图像做中值滤波;
4、基于效率考虑,opencv中实现的霍夫变换圆检测是基于图像梯度的实现,分为两步:
①检测边缘,发现可能的圆心;
②基于第一步的基础上从候选圆心开始计算最佳半径大小;

x = a + Rcos(θ)

y = b + Rsin(θ)

from matplotlib import pyplot as plt 
from cv2 import cv2 as cv
import numpy as np 

def detect_circles_demo(image):
    #边缘保留滤波   #也可用高斯滤波
    dst = cv.pyrMeanShiftFiltering(image,10,100)
    '''
    scr:输入的8比特,3信道图像
    dst:和源图像相同大小、格式的输出图像
    sp:空间窗的半径
    sr:色彩窗的半径
    '''
    cimage = cv.cvtColor(dst,cv.COLOR_BGR2GRAY)
    circles = cv.HoughCircles(cimage,cv.HOUGH_GRADIENT,1,20,param1=50,param2=28,minRadius=0,maxRadius=0)
    circles1 = np.uint16(np.around(circles))#转换成整数
    for i in circles1[0,:]:
        cv.circle(image,(i[0],i[1]),i[2],(0,0,255),2)#画圆
        cv.circle(image,(i[0],i[1]),2,(255,0,0),2)#圆心
    cv.imshow("circles",image)

if __name__ == "__main__":
    filepath = "C:\\pictures\\Test paper\\unqualified\\unqualifed-1.jpg"
    img = cv.imread(filepath)       # blue green red
    cv.namedWindow("input image",cv.WINDOW_AUTOSIZE)
    cv.imshow("input image",img)

    detect_circles_demo(img)

    cv.waitKey(0)
    cv.destroyAllWindows()

结果如下:
在这里插入图片描述

PythonOpenCV中实现直线检测,可以使用Hough变换来检测直线。Hough变换是一种常用的图像处理方法,可用于检测直线、等几何形状。 以下是一个简单的示例代码,使用Hough变换来检测直线并计算交点: ```python import cv2 import numpy as np # 读取图像 img = cv2.imread('test.jpg') gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) # 边缘检测 edges = cv2.Canny(gray, 50, 150, apertureSize=3) # Hough变换检测直线 lines = cv2.HoughLines(edges, 1, np.pi/180, 200) # 计算交点 points = [] for i in range(len(lines)): for j in range(i+1, len(lines)): rho1, theta1 = lines[i][0] rho2, theta2 = lines[j][0] if abs(theta1 - theta2) < np.pi/4: continue a = np.array([[np.cos(theta1), np.sin(theta1)], [np.cos(theta2), np.sin(theta2)]]) b = np.array([rho1, rho2]) x0, y0 = np.linalg.solve(a, b) points.append((int(x0), int(y0))) # 绘制直线和交点 for line in lines: rho, theta = line[0] a = np.cos(theta) b = np.sin(theta) x0 = a*rho y0 = b*rho x1 = int(x0 + 1000*(-b)) y1 = int(y0 + 1000*(a)) x2 = int(x0 - 1000*(-b)) y2 = int(y0 - 1000*(a)) cv2.line(img, (x1,y1), (x2,y2), (0,0,255), 2) for point in points: cv2.circle(img, point, 5, (0,255,0), -1) # 显示图像 cv2.imshow('image', img) cv2.waitKey(0) cv2.destroyAllWindows() ``` 在代码中,首先读取图像并进行灰度转换和边缘检测。然后使用Hough变换检测直线,并计算交点。最后绘制直线和交点,并显示图像。 需要注意的是,在计算交点时,需要将两条直线的极坐标表示转换为直角坐标表示,并使用线性方程组求解。 希望这个例子能够帮助到你实现直线检测并计算交点。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

xiao黄

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值