CN103996191B - 一种基于最小二乘优化的黑白棋盘格图像角点的检测方法 - Google Patents
一种基于最小二乘优化的黑白棋盘格图像角点的检测方法 Download PDFInfo
- Publication number
- CN103996191B CN103996191B CN201410195996.6A CN201410195996A CN103996191B CN 103996191 B CN103996191 B CN 103996191B CN 201410195996 A CN201410195996 A CN 201410195996A CN 103996191 B CN103996191 B CN 103996191B
- Authority
- CN
- China
- Prior art keywords
- black
- angle point
- image
- gridiron pattern
- pattern image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 31
- 238000006243 chemical reaction Methods 0.000 claims abstract description 4
- 238000001914 filtration Methods 0.000 claims description 6
- 238000002187 spin decoupling employing ultra-broadband-inversion sequences generated via simulated annealing Methods 0.000 description 7
- 238000007689 inspection Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 230000002452 interceptive effect Effects 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种基于最小二乘优化的黑白棋盘格图像角点的检测方法,将彩色棋盘格图像转化为灰度棋盘格图像,将棋盘格图像进行二值化处理得到二值化图像,根据二值化图像和灰度棋盘图像确定出角点区域,根据角点区域内的像素点坐标用最小二乘拟合圆方法拟合出圆,该圆心即为该棋盘格图像角点,从而生成角点坐标集。该方法无需标定相机镜头的参数,适合检测任意畸变的黑白棋盘格角点,在摄像机标定、全景泊车、视觉检测和单目测距等领域有广泛的应用。
Description
技术领域
本发明涉及机器视觉领域,具体涉及一种基于最小二乘优化的黑白棋盘格图像角点的检测方法。
背景技术
黑白棋盘格模板已经广泛应用于摄像机标定获取镜头内参和外参的过程,前期的摄像头标定得到的参数的准确性决定了后期各种图像应用的精度,而黑白棋盘格标定板角点位置提取的精度很大情况下决定了摄像机标定的精度。
角点检测的方法多种多样,但大致上可以分为4类:基于边缘特征的角点检测、基于灰度图像的角点检测、基于二值图像的角点检测和数学形态学。大家常用的采用鼠标点击式的人机交互式方法是基于边缘特征的角点检测方法,也有人根据棋盘格的矩形或正方形的大小设计出的矩形模板进行比对提取方法是基于二值图像的角点检测和数学形态学,还有人设计出的Harris角点提取或SUSAN方法等专门的方法是基于灰度图像的角点检测方法,但结果都差强人意,满足不了嵌入式设备和移动式设备的方法快速性的需求。
其中最为简单的鼠标点击式的人机交互式方法通过人为选取棋盘角点位置用计算机显示,所提取到的角点坐标位置准确度高,但是算法周期长,无法在计算机设备上实现,满足不了实时性的要求。根据棋盘格的矩形或正方形的大小设计矩形模板进行比对提取的方法比较适用于图像畸变没有或者畸变较小的情况,然而摄像机标定就是为了通过畸变棋盘格角点提取来确定摄像机镜头畸变角度,实际情况下黑白棋盘格畸变较大,并且可能由于光照不均,边缘信息丢失较严重,用矩形模板提取到的角点在边缘处误差较大。Harris角点检测方法是由Chris Harris和MikeStephens在1988年提出,该方法是在Moravec方法的基础上发展起来的。Harris角点检测方法是在一个局部窗口在不同方向进行少量的偏移后,考察窗口内图像亮度值的平均变化使用灰度的高斯梯度作为响应值判断棋盘格图像中的角点,然而,在棋盘格图像中,由于角点处图像的模糊现象使得实际角点附近的一个或者多个点的Harris响应值较高,很难准确确定角点所在位置。SUSAN方法是关于特征点检测的一种直观的新方法,其思想简单,易于实现,由于不需要微分运算,所以SUSAN方法抗噪能力强、运算速度快。SUSAN该方法用一个预定义的掩板对图像中的每个象素点进行局部的运算操作,给出相应的边缘响应,这个响应就作为输出得到图像的边缘。但SUSAN方法也有不足:由于其在检测角点时只考虑SUSAN区域的大小,而不考虑SUSAN区域的形状,所以对一些特殊类型的角点会出现漏检的情况。
发明内容
针对现有技术的不足,本发明提供了一种基于最小二乘优化的黑白棋盘格图像角点的检测方法。
本发明的技术方案是:
一种基于最小二乘优化的黑白棋盘格图像角点的检测方法,包括以下步骤:
步骤1、将彩色棋盘格图像转化为灰度棋盘格图像。
步骤2、调整二值化阈值大小,使灰度棋盘格图像中每一个黑色格子与周围四个黑色格子分开,从而确定二值化阈值,对灰度棋盘格图像进行二值化处理,得到二值化图像。
步骤3、运用邻域滤波法去掉二值化图像的噪声,根据棋盘形态学边界屏蔽边界外的信息。
步骤4、将二值化图像中的白色格子区域中与灰度棋盘格图像中两个黑色格子邻接区域重合的区域作为角点区域。
步骤5、根据角点区域内的像素点坐标用最小二乘拟合圆方法拟合出圆,该圆的圆心即为黑白棋盘格图像角点。
步骤6:得到的黑白棋盘格图像的各角点坐标,生成角点坐标集,即完成检测。
本发明的有益效果是:采用提取角点所在区域,利用最小二乘拟合圆方法得到圆心作为角点的思路,解决了传统方法对角点界定难的问题,提高了角点位置提取的精度和准确度,算法简洁、实时性强,利于嵌入式设备和移动式设备的方法移植,能够比较好的提取出任意畸变的黑白棋盘格图像角点。
附图说明
图1为本发明具体实施方式的基于最小二乘优化的黑白棋盘格图像角点的检测方法的流程图;
图2为本发明具体实施方式的基于最小二乘优化的黑白棋盘格图像角点的检测方法中二值化图像部分区域放大图;
图3为本发明具体实施方式的基于最小二乘优化的黑白棋盘格图像角点的检测方法中二值化后的图像;
图4为本发明具体实施方式的基于最小二乘优化的黑白棋盘格图像角点的检测方法中消除噪声和棋盘外界环境影响后的二值化图像;
图5为本发明具体实施方式的基于最小二乘优化的黑白棋盘格图像角点的检测方法中所有提取出的角点图。
具体实施方式
下面结合附图对本发明一种实施例做进一步说明。
本实施方式的软件环境为WINDOWS 7***,仿真环境为MATLAB R2014a,处理器为第二代酷睿i5-2410。
本实施方式的基于最小二乘优化的黑白棋盘格图像角点的检测方法的流程图,如图1所示,包括以下步骤:
步骤1、将彩色棋盘格图像转化为灰度棋盘格图像。
利用MATLAB R2014a***灰度转化函数将彩色棋盘图像转化为灰度棋盘格图像。
步骤2、调整二值化阈值大小,使灰度棋盘格图像中每一个黑色格子与周围四个黑色格子分开,从而确定二值化阈值,对灰度棋盘格图像进行二值化处理,得到二值化图像。
MATLAB***二值化方法无法满足本发明对二值化后图像的要求,因此本发明根据灰度棋盘格图像全局灰度大小调节二值化阈值获得二值化图像,令灰度棋盘格图像灰度f(x,y),二值化图像灰度为g(x,y)。根据图像全局灰度大小调节二值化阈值Threshold的大小,使得棋盘每一个黑色的格子与周围四个黑色格子分隔开。
如图2所示,本实施方式中调节二值化阈值Threshold的大小获得的二值化图像部分区域放大图,黑色格子与周围黑色格子分离2~3个像素点,本次确定的二值化阈值为Threshold=122。
对灰度棋盘格图像二值化处理后得到二值化图像如图3所示。
步骤3、运用邻域滤波法去掉二值化图像的噪声,根据棋盘形态学边界屏蔽边界外的信息。
本实施方式利用邻域滤波法对单个噪点和双噪点进行消除,在消除噪点的时候先进行横向邻域滤波,再进行纵向邻域滤波。针对棋盘外界环境影响,利用棋盘形态学边界,将边界外的无用信息全部消除。如图4所示,为消除噪声和棋盘外界环境影响后的二值化图像。
步骤4、将二值化图像中的白色格子区域中与灰度棋盘格图像中两个黑色格子邻接区域重合的区域作为角点区域。
由于步骤二中调整二值化阈值Threshold的大小已将棋盘每一个黑色的格子与周围四个黑色格子分离开,即角点已经与黑色的格子分离,角点被包含在两黑色格子相邻的区域,需要将在灰度棋盘格图像中黑色格子相邻区域同时在二值化图像中不属于黑色格子的区域确定为角点区域。
步骤5、根据角点区域内的像素点坐标用最小二乘拟合圆方法拟合出圆,该圆的圆心即为黑白棋盘格图像角点。
设待拟合的圆方程为
x2+y2+mx+ny+c=0 (2)
其中,m、n、c分别为待拟合圆的参数。
点(xi,yi)到圆心的平方和与半径平方的差为
令Q(m,n,c)为di的平方和:
其中N为边界点个数。
利用偏导数为零求得m、n、c,即
解线性方程组(4),令
解得m、n、c为
那么该鱼眼图像的畸变中心u0、v0和畸变半径R分别为:
将步骤4得到的角点区域点集代入公式(7),得到m=-534,n=-1292,c=488589,代入公式(8),那么角点为(u0,v0)=(267,646)。
如图5所示,为所有提取出的角点图。
步骤6:得到的黑白棋盘格图像的各角点坐标,生成角点坐标集,即完成检测。
本实施方式提取的320个角点的坐标(坐标顺序为从左到右)如下:
Claims (1)
1.一种基于最小二乘优化的黑白棋盘格图像角点的检测方法,其特征在于,包括以下步骤:
步骤1、将彩色棋盘格图像转化为灰度棋盘格图像;
步骤2、调整二值化阈值大小,使灰度棋盘格图像中每一个黑色格子与周围四个黑色格子分开,从而确定二值化阈值,对灰度棋盘格图像进行二值化处理,得到二值化图像;
步骤3、运用邻域滤波法去掉二值化图像的噪声,根据棋盘形态学边界屏蔽边界外的信息;
步骤4、将二值化图像中的白色格子区域中与灰度棋盘格图像中两个黑色格子邻接区域重合的区域作为角点区域;
步骤5、根据角点区域内的像素点坐标用最小二乘拟合圆方法拟合出圆,该圆的圆心即为黑白棋盘格图像角点;
步骤6:得到的黑白棋盘格图像的各角点坐标,生成角点坐标集,即完成检测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410195996.6A CN103996191B (zh) | 2014-05-09 | 2014-05-09 | 一种基于最小二乘优化的黑白棋盘格图像角点的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410195996.6A CN103996191B (zh) | 2014-05-09 | 2014-05-09 | 一种基于最小二乘优化的黑白棋盘格图像角点的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103996191A CN103996191A (zh) | 2014-08-20 |
CN103996191B true CN103996191B (zh) | 2016-10-26 |
Family
ID=51310346
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410195996.6A Expired - Fee Related CN103996191B (zh) | 2014-05-09 | 2014-05-09 | 一种基于最小二乘优化的黑白棋盘格图像角点的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103996191B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105818746A (zh) * | 2015-01-05 | 2016-08-03 | 上海纵目科技有限公司 | 环视高级辅助驾驶***的标定方法及*** |
CN106201062B (zh) * | 2016-06-23 | 2019-05-07 | 南京仁光电子科技有限公司 | 一种基于递归算法的棋盘格角点定位方法 |
CN107231553A (zh) * | 2017-05-31 | 2017-10-03 | 信利光电股份有限公司 | 角点位置获取方法及装置 |
CN107845101B (zh) * | 2017-09-20 | 2020-04-10 | 深圳市易成自动驾驶技术有限公司 | 车载环视图像的特征点标定方法、装置和可读存储介质 |
CN107480678A (zh) * | 2017-09-29 | 2017-12-15 | 北京深度奇点科技有限公司 | 一种棋盘识别方法和识别*** |
CN108399412A (zh) * | 2018-03-07 | 2018-08-14 | 鞍钢集团矿业有限公司 | 一种x型角点亚像素提取方法 |
CN109003311B (zh) * | 2018-08-22 | 2022-03-18 | 上海庄生晓梦信息科技有限公司 | 一种鱼眼镜头的标定方法 |
CN111047639B (zh) * | 2018-10-12 | 2023-08-18 | 深圳富联富桂精密工业有限公司 | 单目测距装置 |
CN109509200B (zh) * | 2018-12-26 | 2023-09-29 | 深圳市繁维医疗科技有限公司 | 基于轮廓提取的棋盘格角点检测方法以及计算机可读存储介质 |
CN112446895B (zh) * | 2019-08-28 | 2024-07-09 | 深圳市光鉴科技有限公司 | 棋盘格角点自动提取方法、***、设备及介质 |
CN112284330A (zh) * | 2020-11-18 | 2021-01-29 | 青岛科美创视智能科技有限公司 | 一种基于鱼眼摄像机的距离测量***及方法 |
CN117157668A (zh) * | 2021-03-29 | 2023-12-01 | 哈曼国际工业有限公司 | 对dvs与相机之间的统一标定板进行角点检测的方法 |
CN114529616B (zh) * | 2022-04-22 | 2022-07-26 | 武汉精视遥测科技有限公司 | 基于内壁刻度的广角镜头参数标定方法、***及计算机 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103093451A (zh) * | 2011-11-03 | 2013-05-08 | 北京理工大学 | 一种棋盘格交叉点识别算法 |
CN103177439A (zh) * | 2012-11-26 | 2013-06-26 | 惠州华阳通用电子有限公司 | 一种基于黑白格角点匹配的自动标定方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2465793A (en) * | 2008-11-28 | 2010-06-02 | Sony Corp | Estimating camera angle using extrapolated corner locations from a calibration pattern |
-
2014
- 2014-05-09 CN CN201410195996.6A patent/CN103996191B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103093451A (zh) * | 2011-11-03 | 2013-05-08 | 北京理工大学 | 一种棋盘格交叉点识别算法 |
CN103177439A (zh) * | 2012-11-26 | 2013-06-26 | 惠州华阳通用电子有限公司 | 一种基于黑白格角点匹配的自动标定方法 |
Non-Patent Citations (3)
Title |
---|
An Automatic and Robust Chessboard Corner Extraction;Xifan Shi et al;《Lecture Notes in Computer Science》;20111231;第6758卷;49-56 * |
Automatic Distortion Correction of Endoscopic Images Captured With Wide-Angle Zoom Lens;Tung-Ying Lee et al;《IEEE Transactions on Biomedical Engineering》;20130930;第60卷(第9期);2603 - 2613 * |
一种自动检测棋盘角点的新算法;*** 等;《计算机工程》;20040731;第30卷(第14期);19-21,35 * |
Also Published As
Publication number | Publication date |
---|---|
CN103996191A (zh) | 2014-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103996191B (zh) | 一种基于最小二乘优化的黑白棋盘格图像角点的检测方法 | |
WO2021138995A1 (zh) | 一种棋盘格角点全自动检测方法 | |
WO2018024030A1 (zh) | 一种夜视红外图像中的基于显著性的道路目标提取方法 | |
WO2017084204A1 (zh) | 一种二维视频流中的人体骨骼点追踪方法及*** | |
WO2018023916A1 (zh) | 一种彩色图像去阴影方法和应用 | |
WO2017128605A1 (zh) | 一种教育玩具套件及其基于hsv的魔方颜色识别方法 | |
CN103679677B (zh) | 一种基于模型互更新的双模图像决策级融合跟踪方法 | |
CN104077759A (zh) | 一种基于色觉感知及全局质量因子的多曝光度图像融合方法 | |
CN102930534A (zh) | 人体背部穴位自动定位方法 | |
CN103198319B (zh) | 用于矿山井筒环境下的模糊图像角点提取方法 | |
CN104268853A (zh) | 一种红外图像与可见光图像配准方法 | |
CN103455813A (zh) | 一种ccd图像测量***光斑中心定位的方法 | |
CN101916370A (zh) | 人脸检测中非特征区域图像处理的方法 | |
CN105844655A (zh) | 一种激光线条纹边缘提取方法 | |
CN104008542B (zh) | 一种针对特定平面图形的快速角点匹配方法 | |
CN108921881A (zh) | 一种基于单应性约束的跨摄像头目标跟踪方法 | |
CN103218605A (zh) | 一种基于积分投影与边缘检测的快速人眼定位方法 | |
CN106162124A (zh) | 投影机输出图像的校准方法 | |
CN103093451A (zh) | 一种棋盘格交叉点识别算法 | |
CN105147311A (zh) | 用于ct***中的可视化设备辅助扫描定位方法和*** | |
CN106408580A (zh) | 一种基于模糊c均值和均值漂移的肝脏区域提取方法 | |
CN104537688A (zh) | 一种基于背景差分和hog特征的运动目标检测方法 | |
CN106686347B (zh) | 一种基于视频的判定地铁摄像机移位的方法 | |
CN106327464A (zh) | 一种边缘检测方法 | |
CN105550981A (zh) | 一种基于Lucas-Kanade算法的图像配准和拼接方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161026 |
|
CF01 | Termination of patent right due to non-payment of annual fee |