使用Python，OpenCV对图像进行亚像素点检测，并拟合椭圆进行绘制
2024-06-01 19:38:46
这篇博客将介绍如何使用Python，OpenCV对图像进行亚像素检测，并对亚像素点进行椭圆拟合绘制。
1. 效果图
原始图上绘制拟合椭圆 VS 原始图上绘制拟合椭圆及亚像素点绘制随机半径及颜色的圆 VS 灰度图上绘制亚像素点效果图如下：
我喜欢的颖宝效果图2如下：
圆的随机半径放大一点效果图如下：
2. 源码
需要注意，绘制圆需要圆心都是int值
# 计算图像的亚像素点并绘制（matplotlib以及cv2），并对亚像素点进行3种椭圆方式拟合绘制
# python zernick_detection.pyimport timeimport cv2
import imutils
import matplotlib.pyplot as plt
import numpy as npg_N = 7M00 = np.array([0, 0.0287, 0.0686, 0.0807, 0.0686, 0.0287, 0,0.0287, 0.0815, 0.0816, 0.0816, 0.0816, 0.0815, 0.0287,0.0686, 0.0816, 0.0816, 0.0816, 0.0816, 0.0816, 0.0686,0.0807, 0.0816, 0.0816, 0.0816, 0.0816, 0.0816, 0.0807,0.0686, 0.0816, 0.0816, 0.0816, 0.0816, 0.0816, 0.0686,0.0287, 0.0815, 0.0816, 0.0816, 0.0816, 0.0815, 0.0287,0, 0.0287, 0.0686, 0.0807, 0.0686, 0.0287, 0]).reshape((7, 7))M11R = np.array([0, -0.015, -0.019, 0, 0.019, 0.015, 0,-0.0224, -0.0466, -0.0233, 0, 0.0233, 0.0466, 0.0224,-0.0573, -0.0466, -0.0233, 0, 0.0233, 0.0466, 0.0573,-0.069, -0.0466, -0.0233, 0, 0.0233, 0.0466, 0.069,-0.0573, -0.0466, -0.0233, 0, 0.0233, 0.0466, 0.0573,-0.0224, -0.0466, -0.0233, 0, 0.0233, 0.0466, 0.0224,0, -0.015, -0.019, 0, 0.019, 0.015, 0]).reshape((7, 7))M11I = np.array([0, -0.0224, -0.0573, -0.069, -0.0573, -0.0224, 0,-0.015, -0.0466, -0.0466, -0.0466, -0.0466, -0.0466, -0.015,-0.019, -0.0233, -0.0233, -0.0233, -0.0233, -0.0233, -0.019,0, 0, 0, 0, 0, 0, 0,0.019, 0.0233, 0.0233, 0.0233, 0.0233, 0.0233, 0.019,0.015, 0.0466, 0.0466, 0.0466, 0.0466, 0.0466, 0.015,0, 0.0224, 0.0573, 0.069, 0.0573, 0.0224, 0]).reshape((7, 7))M20 = np.array([0, 0.0225, 0.0394, 0.0396, 0.0394, 0.0225, 0,0.0225, 0.0271, -0.0128, -0.0261, -0.0128, 0.0271, 0.0225,0.0394, -0.0128, -0.0528, -0.0661, -0.0528, -0.0128, 0.0394,0.0396, -0.0261, -0.0661, -0.0794, -0.0661, -0.0261, 0.0396,0.0394, -0.0128, -0.0528, -0.0661, -0.0528, -0.0128, 0.0394,0.0225, 0.0271, -0.0128, -0.0261, -0.0128, 0.0271, 0.0225,0, 0.0225, 0.0394, 0.0396, 0.0394, 0.0225, 0]).reshape((7, 7))M31R = np.array([0, -0.0103, -0.0073, 0, 0.0073, 0.0103, 0,-0.0153, -0.0018, 0.0162, 0, -0.0162, 0.0018, 0.0153,-0.0223, 0.0324, 0.0333, 0, -0.0333, -0.0324, 0.0223,-0.0190, 0.0438, 0.0390, 0, -0.0390, -0.0438, 0.0190,-0.0223, 0.0324, 0.0333, 0, -0.0333, -0.0324, 0.0223,-0.0153, -0.0018, 0.0162, 0, -0.0162, 0.0018, 0.0153,0, -0.0103, -0.0073, 0, 0.0073, 0.0103, 0]).reshape(7, 7)M31I = np.array([0, -0.0153, -0.0223, -0.019, -0.0223, -0.0153, 0,-0.0103, -0.0018, 0.0324, 0.0438, 0.0324, -0.0018, -0.0103,-0.0073, 0.0162, 0.0333, 0.039, 0.0333, 0.0162, -0.0073,0, 0, 0, 0, 0, 0, 0,0.0073, -0.0162, -0.0333, -0.039, -0.0333, -0.0162, 0.0073,0.0103, 0.0018, -0.0324, -0.0438, -0.0324, 0.0018, 0.0103,0, 0.0153, 0.0223, 0.0190, 0.0223, 0.0153, 0]).reshape(7, 7)M40 = np.array([0, 0.013, 0.0056, -0.0018, 0.0056, 0.013, 0,0.0130, -0.0186, -0.0323, -0.0239, -0.0323, -0.0186, 0.0130,0.0056, -0.0323, 0.0125, 0.0406, 0.0125, -0.0323, 0.0056,-0.0018, -0.0239, 0.0406, 0.0751, 0.0406, -0.0239, -0.0018,0.0056, -0.0323, 0.0125, 0.0406, 0.0125, -0.0323, 0.0056,0.0130, -0.0186, -0.0323, -0.0239, -0.0323, -0.0186, 0.0130,0, 0.013, 0.0056, -0.0018, 0.0056, 0.013, 0]).reshape(7, 7)def zernike_detection(path):img = cv2.imread(path)img = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY)blur_img = cv2.medianBlur(img, 13)c_img = cv2.adaptiveThreshold(blur_img, 255, cv2.ADAPTIVE_THRESH_GAUSSIAN_C, cv2.THRESH_BINARY_INV, 7, 4)ZerImgM00 = cv2.filter2D(c_img, cv2.CV_64F, M00)ZerImgM11R = cv2.filter2D(c_img, cv2.CV_64F, M11R)ZerImgM11I = cv2.filter2D(c_img, cv2.CV_64F, M11I)ZerImgM20 = cv2.filter2D(c_img, cv2.CV_64F, M20)ZerImgM31R = cv2.filter2D(c_img, cv2.CV_64F, M31R)ZerImgM31I = cv2.filter2D(c_img, cv2.CV_64F, M31I)ZerImgM40 = cv2.filter2D(c_img, cv2.CV_64F, M40)point_temporary_x = []point_temporary_y = []scatter_arr = cv2.findNonZero(ZerImgM00).reshape(-1, 2)for idx in scatter_arr:j, i = idxtheta_temporary = np.arctan2(ZerImgM31I[i][j], ZerImgM31R[i][j])rotated_z11 = np.sin(theta_temporary) * ZerImgM11I[i][j] + np.cos(theta_temporary) * ZerImgM11R[i][j]rotated_z31 = np.sin(theta_temporary) * ZerImgM31I[i][j] + np.cos(theta_temporary) * ZerImgM31R[i][j]l_method1 = np.sqrt((5 * ZerImgM40[i][j] + 3 * ZerImgM20[i][j]) / (8 * ZerImgM20[i][j]))l_method2 = np.sqrt((5 * rotated_z31 + rotated_z11) / (6 * rotated_z11))l = (l_method1 + l_method2) / 2k = 3 * rotated_z11 / (2 * (1 - l_method2 ** 2) ** 1.5)# h = (ZerImgM00[i][j] - k * np.pi / 2 + k * np.arcsin(l_method2) + k * l_method2 * (1 - l_method2 ** 2) ** 0.5)# / np.pik_value = 20.0l_value = 2 ** 0.5 / g_Nabsl = np.abs(l_method2 - l_method1)if k >= k_value and absl <= l_value:y = i + g_N * l * np.sin(theta_temporary) / 2x = j + g_N * l * np.cos(theta_temporary) / 2point_temporary_x.append(x)point_temporary_y.append(y)else:continuereturn point_temporary_x, point_temporary_yimage = cv2.imread('ml2.jpg')
path = 'ym_600.jpg'
cv2.imwrite(path, imutils.resize(image, height=600))time1 = time.time()
point_temporary_x, point_temporary_y = zernike_detection(path)
time2 = time.time()
print(time2 - time1)# gray : 进行检测的图像
gray = cv2.imread(path, 0)
plt.imshow(gray, cmap="gray")
# point检测出的亚像素点
point = np.array([point_temporary_x, point_temporary_y])
# s:调整显示点的大小
plt.scatter(point[0, :], point[1, :], s=10, marker="*")
plt.show()def cv_fit_ellipse(points, flag=0):# cv2.fitEllipse only can estimate int numpy 2D array datap = np.array(points)p = np.array(p).astype(int)center0, axes0, angle0 = cv2.fitEllipse(p)center1, axes1, angle1 = cv2.fitEllipseAMS(p)center2, axes2, angle2 = cv2.fitEllipseDirect(p)print(type(center0))print("fitEllipse: " + str(np.array(center0)))print("fitEllipseAMS: " + str(np.array(center1)))print("fitEllipseDirect: " + str(np.array(center2)))box = tuple([tuple(np.array(center0)), tuple(np.array(axes0)), angle0])if flag == 0:box = tuple([tuple(np.array(center0)), tuple(np.array(axes0)), angle0])elif flag == 1:box = tuple([tuple(np.array(center1)), tuple(np.array(axes1)), angle1])else:box = tuple([tuple(np.array(center0)), tuple(np.array(axes0)), angle0])return boxpoints = np.array(list(zip(point_temporary_x, point_temporary_y)))
print(type(points), points[0], points.shape)
box1 = cv_fit_ellipse(points, flag=1)image = cv2.imread(path)cv2.ellipse(img=image, box=box1, color=(255, 0, 255), thickness=1)
cv2.imshow("image", image)
for i in points:# 亚像素点绘制随机颜色半径圆radius = np.random.randint(1, high=10)color = np.random.randint(0, high=256, size=(3,)).tolist()# print('radius: ', radius, type(radius))# print('color: ', color, type(color))cv2.circle(image, (int(i[0]), int(i[1])), radius, color, -1)cv2.imshow("image1", image)
cv2.waitKey(0)
cv2.destroyAllWindows()参考
- Python，OpenCV鼠标事件进行矩形、圆形的绘制（随机颜色、随机半径）
- 图像亚像素点计算
- 椭圆拟合
使用Python，OpenCV对图像进行亚像素点检测，并拟合椭圆进行绘制相关推荐
- Python+OpenCV：图像Harris角点检测(Harris Corner Detection)
Python+OpenCV:图像Harris角点检测(Harris Corner Detection) 理论 corners are regions in the image with large v ...
- Python+OpenCV：图像快速角点检测算法(FAST Algorithm for Corner Detection)
Python+OpenCV:图像快速角点检测算法(FAST Algorithm for Corner Detection) 理论 Feature Detection using FAST Select ...
- Python+OpenCV判断图像是黑底还是白底
前言 本篇博客使用Python+OpenCV判断图像是黑底还是白底,利用图像对角线上的黑白像素点个数进行判断,详情见下文. 本篇博客内容包含代码逻辑.说明.依赖.实现,这几部分.代码实现部分包含2种实 ...
- 【Python OpenCV】图像直方图 calcHist方法 equalizeHist方法
[Python OpenCV]图像直方图 calcHist方法 equalizeHist方法 (一)图像直方图 图像的构成是有像素点构成的,每个像素点的值代表着该点的颜色(灰度图或者彩色图).所谓直方 ...
- Python+OpenCV：图像修复（Image Inpainting）
Python+OpenCV:图像修复(Image Inpainting) 理论 Most of you will have some old degraded photos at your home ...
- Python+OpenCV：图像二进制鲁棒独立基本特征(BRIEF, Binary Robust Independent Elementary Features)
Python+OpenCV:图像二进制鲁棒独立基本特征(BRIEF, Binary Robust Independent Elementary Features) 理论 We know SIFT us ...
- Python+OpenCV：图像Shi-Tomasi角点检测器
Python+OpenCV:图像Shi-Tomasi角点检测器 理论 The scoring function in Harris Corner Detector was given by: Inst ...
- Python+OpenCV：图像对比度受限自适应直方图均衡化(CLAHE, Contrast Limited Adaptive Histogram Equalization)
Python+OpenCV:图像对比度受限自适应直方图均衡化(CLAHE, Contrast Limited Adaptive Histogram Equalization) ############ ...
- Python+OpenCV：图像轮廓
Python+OpenCV:图像轮廓 轮廓是什么? 轮廓可以简单地解释为一条连接所有连续点(沿边界)的曲线,具有相同的颜色和强度. 轮廓线是形状分析.目标检测和识别的重要工具. 为了获得更好的精度,可 ...
最新文章
- Visual Studio 2017 版本 15.5.5
- 教授爆笑讲解超级玛丽
- java我的世界1.7.2怎么下载模组_我的世界1.7.2
- 在centos7中静默安装oracle11g
- 链接访问后刷新颜色回到初始_如何使链接可访问（提示：颜色不够）
- JVM—GC垃圾回收器总结
- 关于单点登录的实现方法
- Asp.Net MVC2 控件开发实例(2)
- zabbix监控memcached
- 提交json数据给Java_java 实现通过 post 方式提交json参数操作
- iOS-常用宏定义大全
- 破解wifi密码最简单教程 就三步
- HTML5软件设计大赛,我院成功举行第十七届山东省大学生软件设计大赛 HTML5创意应用命题决赛...
- nginx缓存以及清除缓存
- 麦当劳和肯德鸡的笑话
- 修真院教学模式三大阶段之复盘项目
- 7、【WebGIS实战】专题篇——API key
- 手机安装青龙面板，低功耗，随时随地的薅羊毛（无需服务器）
- mysql interval weekday_Mysql DATE_SUB(NOW(), INTERVAL 1 DAY) 24 hours or weekday?
- 先验分布、后验分布、似然函数