CN113129396A - 一种基于区域分割的平行四边形编码标志的解码方法 - Google Patents
一种基于区域分割的平行四边形编码标志的解码方法 Download PDFInfo
- Publication number
- CN113129396A CN113129396A CN202110392326.3A CN202110392326A CN113129396A CN 113129396 A CN113129396 A CN 113129396A CN 202110392326 A CN202110392326 A CN 202110392326A CN 113129396 A CN113129396 A CN 113129396A
- Authority
- CN
- China
- Prior art keywords
- coding
- parallelogram
- zeta
- mark
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000011218 segmentation Effects 0.000 title claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims description 94
- 108091026890 Coding region Proteins 0.000 claims description 57
- 238000004422 calculation algorithm Methods 0.000 claims description 14
- 108700026244 Open Reading Frames Proteins 0.000 claims description 11
- 239000003550 marker Substances 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 8
- 238000003491 array Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 230000000007 visual effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000007797 corrosion Effects 0.000 claims description 2
- 238000005260 corrosion Methods 0.000 claims description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 abstract description 2
- 238000003672 processing method Methods 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 10
- 230000008901 benefit Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000003086 colorant Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T9/00—Image coding
- G06T9/20—Contour coding, e.g. using detection of edges
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Multimedia (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于区域分割的平行四边形编码标志的解码方法,该平行四边形编码标志由平行四边形背景图案和编码图案构成,编码图案在平行四边形背景图案的内部,包括定位图案、定向图案和编码标志图案,定向图案和定位图案用于平行四边形编码标志方向的判断;编码标志图案用于编码标志四个顶点进行编码。本发明将平行四边形编码标志分割为六个区域进行编码,因此对于仿射变换的图像也有较高的处理精度;该平行四边形编码标志具有较大的编码容量,在不失准确度的情况下,能应对大尺寸目标或大视场场景下的标定或三维点云拼接;所述的解码方法通过数字图像处理的方法解析出平行四边形编码标志所包含的方向信息和编码信息,从而完成对每个角点的解码。
Description
技术领域
本发明涉及计算机视觉中的视觉测量领域,具体为一种基于区域分割的平行四边形编码标志的解码方法,适用于相机标定、目标特征提取、立体匹配和三维数据拼接等领域。
背景技术
随着现代工业的发展,为提高产品质量和产品效益,解决大尺寸和小尺寸及复杂背景的测量问题,视觉测量技术顺应而生。视觉测量技术具有非接触、不易受温度变化和振动的影响、高精度、便携性好等特点。在视觉测量领域中,大视场、复杂背景下的相机标定困难,以及大尺寸目标如何实现立体匹配和三维数据拼接等问题仍是研究热点。
为有效解决上述问题,编码标记点技术飞速发展,出现的编码标志点类型很多,但大多都为环形、圆形或扇形等,而这些编码标记点在面对仿射变换图像,或大尺寸目标三维拼接时,不能保证结果的准确度和精确度。因此,本发明提出的一种基于区域分割编码的平行四边形编码标志的解码方法,有效弥补了以上缺陷,在实现复杂条件的相机标定的同时,仍保证了较高的精度和鲁棒性,同时由于编码容量较大,面对大尺寸目标的三维数据拼接仍能保持较高的准确性,且操作简便。
发明内容
本发明旨在克服现有编码标志的缺陷,提出了一种基于区域分割的平行四边形编码标志的解码方法,利用分割区域进行编码和解码工作的方法,在保持精确度和准确度的情况下,仍能应对大视场、复杂背景下的相机标定,以及大尺寸目标的立体匹配和三维数据拼接。
为实现上述效果,本发明采用的技术方案为:
一种基于区域分割的平行四边形编码标志的解码方法,所述编码标志为正方形的编码方格,以平行四边形编码标志的四个顶点作为特征角点,在视觉结构光三维测量过程中,用于目标曲面的三维点云拼接;编码方格的表面设置有平行四边形背景图案和编码图案,所述编码图案位于平行四边形背景图案的内部,所述编码图案包括定位图案、定向图案和编码标志图案,所述定向图案和定位图案用于平行四边形编码标志方向的判断,所述编码标志图案用于平行四边形编码标志各个角点的编码。利用数字图像处理的方法获得摄像机拍摄的包含N个平行四边形编码标志的图像中每个特征角点的编码序号和亚像素坐标,进而完成平行四边形编码标志解码的步骤如下:
步骤1.1、利用摄像机拍摄空间中放置的N个平行四边形编码标志,获得编码标志图像,所述编码标志图像当中包含N个平行四边形编码标志;
步骤1.2、建立角点的像素坐标系;
在上述拍摄的包含N个平行四边形编码标志的图像中,取图像左上角作为角点像素坐标系的原点o,自左向右为角点像素坐标系的x轴方向,自上向下作为角点像素坐标系的y轴方向,由此建立角点像素坐标系o-xy;
步骤2.1、对编码标志图像进行灰度处理,得到编码标志灰度图像P1;其中,编码标志灰度图像P1为8位灰度图;
步骤2.2、对编码标志灰度图像P1进行复制备份,得到编码标志灰度备份图像P1′;
步骤2.3、对编码标志灰度备份图像P1′进行二值化处理,使得每个平行四边形编码标志的背景颜色均为黑色(即灰度值为0),每个平行四边形编码标志内的定向图案和定位图案均为白色(即灰度值为255),编码单元图案的颜色根据平行四边形编码标志的编码规则可以为白色也可以为黑色,进而得到编码标志二值化图像P2;
步骤3、设定黑色连通域圆度阀值λ′,像素点个数阀值λ″;
步骤4.1、计算编码标志二值化图像P2中所有黑色连通域的圆度值λn和像素点个数cn,其中每个黑色连通域的圆度值由式(1)得出,
其中,ln代表对应黑色连通域轮廓的周长,sn代表对应黑色连通域的面积,n=1,2,3,4…;
步骤4.2、选择合适的黑色连通域圆度阀值λ′和像素点个数阀值λ″可以在编码标志二值化图像P2中获得N个平行四边形编码标志中各自所对应的N个定向圆环中心的黑色连通域;将N个定向圆环中心的黑色连通域分别记为编码标志二值化图像P2中的圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域ΜN,并将圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域ΜN依次放入圆环黑色连通域集合A′中,即圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域ΜN分别为圆环黑色连通域集合A′中的第1个元素、第2个元素、…、第N个元素;由此获得编码标志二值化图像P2中的所有N个平行四边形编码标志中定向圆环内部的黑色连通域;其中,编码标志二值化图像P2中的圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域ΜN的圆度值均小于圆度阀值λ′,并且圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域ΜN的像素点个数均大于像素点个数阀值λ″;
步骤4.3、取整数变量i并赋值i=1;
步骤4.4、计算圆环黑色连通域集合A′中的第i个圆环黑色连通域Μi的质心像素坐标并记为o″d,i(x″d,i,y″d,i),将得到的第i个圆环黑色连通域Μi的质心像素坐标o″d,i(x″d,i,y″d,i)作为圆环黑色连通域质心集合A的第i元素;
步骤4.5、判断i是否小于N,若i<N,则将i+1赋值给i,返回步骤4.4顺序执行;否则,即获得了编码标志二值化图像P2中各个圆环黑色连通域的质心像素坐标o″d,1(x″d,1,y″d,1)、o″d,2(x″d,2,y″d,2)、…、o″d,N(x″d,N,y″d,N),并依次放入圆环黑色连通域质心集合A中;
步骤5、取整数变量ζ并赋值ζ=1;
步骤6.1、将编码标志二值化图像P2进行两次复制备份,分别得到第ζ组第1个备份二值化图像P′ζ,1和第ζ组第2个备份二值化图像P′ζ,2;
步骤6.2、利用去除编码标志复杂背景算法对第ζ组第1个备份二值化图像P′ζ,1进行数字图像处理,进而得到第ζ个无复杂背景平行四边形二值化图像P″ζ,1;在所述第ζ个无复杂背景平行四边形二值化图像P″ζ,1中,包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志内部的所有像素点的灰度值被置为255(即白色),包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志外部的所有像素点的灰度值值被置为0(即黑色);
步骤7.1、对编码标志灰度图像P1再次进行复制备份,获得第ζ个备份灰度图像Pζ″′;
步骤7.2、利用角点检测算法在第ζ个无复杂背景平行四边形二值化图像P″ζ,1中提取4个角点的亚像素坐标,并记为第ζ个无复杂背景平行四边形二值化图像P″ζ,1上的4个备选角点C″ζ,1(x″ζ,1,y″ζ,1)、C″ζ,2(x″ζ,2,y″ζ,2)、C″ζ,3(x″ζ,3,y″ζ,3)、C″ζ,4(x″ζ,4,y″ζ,4);
步骤7.3、在第ζ个无复杂背景平行四边形二值化图像P″ζ,1中做黑色连通域腐蚀,使得第ζ个无复杂背景平行四边形二值化图像P″ζ,1中的白色连通域扩大,在后续进行图像乘法运算时,能够保证运算结果图像中提取出完整的包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志的灰度图案且其4个特征角点保存完整;并将第ζ个无复杂背景平行四边形二值化图像P″ζ,1经本步骤处理所得的图像记为第ζ个运算二值化图像P″ζ,2;
步骤7.4、将第ζ个运算二值化图像P″ζ,2中灰度值为255的部分的值置为1,将第ζ个运算二值化图像P″ζ,2中灰度值为0的部分的值置为0,而后利用第ζ个运算二值化图像P″ζ,2和第ζ个备份灰度图像Pζ″′进行乘法运算,将得到的结果图像记为第ζ个覆盖结果图像Tζ;所述第ζ个覆盖结果图像Tζ中只含有包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志的图像信息,去除了其它平行四边形编码标志的图像信息;
步骤7.5、利用角点检测算法在第ζ个覆盖结果图像Tζ中提取所有角点的亚像素坐标并存入第ζ个初选角点亚像素坐标集合Bζ中;其中,第ζ个初选角点亚像素坐标集合Bζ中存放了第ζ个覆盖结果图像Tζ上包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志的4个特征角点亚像素坐标值和干扰角点的亚像素坐标值,所述干扰角点是除特征角点以外的干扰角点;
步骤7.6、取整数变量i并赋值i=1;
步骤7.7、在所述第ζ个初选角点亚像素坐标集合Bζ中寻找距离亚像素坐标值(x″ζ,i,y″ζ,i)最近的角点亚像素坐标并记为C′ζ,i(x′ζ,i,y′ζ,i);
步骤7.8、判断i是否小于4,若i<4,则将i+1赋值给i并返回执行步骤7.7;否则,此时已经获得了第ζ个覆盖结果图像Tζ上包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志的4个特征角点的亚像素坐标C′ζ,1(x′ζ,1,y′ζ,1)、C′ζ,2(x′ζ,2,y′ζ,2)、C′ζ,3(x′ζ,3,y′ζ,3)、C′ζ,4(x′ζ,4,y′ζ,4),且执行步骤8.1;
步骤8.1、在第ζ组第2个备份二值化图像P′ζ,2上,将包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志记为第ζ个平行四边形编码标志;
步骤8.2、在第ζ组第2个备份二值化图像P′ζ,2上,选取亚像素坐标值分别为(x′ζ,1,y′ζ,1)、(x′ζ,2,y′ζ,2)、(x′ζ,3,y′ζ,3)、(x′ζ,4,y′ζ,4)的4个像素点作为第ζ个平行四边形编码标志的角点四边形Sζ的4个顶点,将4个顶点相连进而获得第ζ个平行四边形编码标志的角点四边形Sζ;
步骤9、在第ζ组第2个备份二值化图像P′ζ,2上,寻找距离质心像素坐标值(x″d,ζ,y″d,ζ)最近的白色连通域,并将此白色连通域的灰度值赋值为0;
步骤10、在第ζ组第2个备份二值化图像P′ζ,2上,第ζ个平行四边形编码标志的角点四边形Sζ内部的所有像素点的灰度值保持不变,将第ζ个平行四边形编码标志的角点四边形Sζ以外的所有像素点的灰度值赋值为255;
步骤11.1、在第ζ组第2个备份二值化图像P′ζ,2上,提取第ζ个平行四边形编码标志的角点四边形Sζ中所有的内、外轮廓,记为第ζ个平行四边形编码标志内的轮廓集合Dζ;
步骤11.2、统计第ζ个平行四边形编码标志内的轮廓集合Dζ内每个轮廓包含的像素点数,将包含像素点数第二多的轮廓记为第ζ组第2个备份二值化图像P′ζ,2上第ζ个平行四边形编码标志内的定位圆轮廓Gζ,计算其质心像素坐标并记为第ζ组第2个备份二值化图像P′ζ,2上第ζ个平行四边形编码标志内的定位圆质心像素坐标o′l,ζ(x′l,ζ,y′l,ζ);
步骤12.1、根据第ζ个平行四边形编码标志内的轮廓集合Dζ中,除去包含像素点数最多的两个轮廓,判断剩余轮廓个数κζ的大小,分为如下情况:若κζ==0,则执行步骤13;否则这κζ个轮廓即为第ζ个平行四边形编码标志内编码标志图案的轮廓,分别记为编码标志圆编码标志圆…、编码标志圆
步骤12.2、给整数变量i并赋予初值i=1;
步骤12.3、在第ζ组第2个备份二值化图像P′ζ,2上,计算编码标志圆轮廓的质心将i+1重新赋值给i后继续执行此步骤,直到i>κζ结束;由此可以得到对应编码标志圆编码标志圆…、编码标志圆的质心像素坐标
步骤13.1、在第ζ组第2个备份二值化图像P′ζ,2上,记在第ζ个平行四边形编码标志中位于第1编码区域、第3编码区域、第4编码区域和第6编码区域的特征角点的亚像素坐标分别为Cζ,1(xζ,1,yζ,1)、Cζ,2(xζ,2,yζ,2)、Cζ,3(xζ,3,yζ,3)、Cζ,4(xζ,4,yζ,4);并在第ζ组第2个备份二值化图像P′ζ,2上,将像素坐标值为的(x″d,ζ,y″d,ζ)的像素点记为第ζ组第2个备份二值化图像P′ζ,2上第ζ个平行四边形编码标志的定向圆环质心o′d,ζ(x′d,ζ,y′d,ζ),并将亚像素坐标值为(x′ζ,1,y′ζ,1)、(x′ζ,2,y′ζ,2)、(x′ζ,3,y′ζ,3)、(x′ζ,4,y′ζ,4)的4个像素点分别记为C″′ζ,1(x″′ζ,1,y″′ζ,1)、C″′ζ,2(x″′ζ,2,y″′ζ,2)、C″′ζ,3(x″′ζ,3,y″′ζ,3)、C″′ζ,4(x″′ζ,4,y″′ζ,4);
步骤14、在第ζ组第2个备份二值化图像P′ζ,2上,找出4个像素点C″′ζ,1(x″′ζ,1,y″′ζ,1)、C″′ζ,2(x″′ζ,2,y″′ζ,2)、C″′ζ,3(x″′ζ,3,y″′ζ,3)、C″′ζ,4(x″′ζ,4,y″′ζ,4)中距离定位圆质心的像素坐标o′l,ζ(x′l,ζ,y′l,ζ)最近的2个像素点,分别记为Cζ,1min(xζ,1min,yζ,1min)和Cζ,2min(xζ,2min,yζ,2min);
分以下两种情况执行:
(1)若sinαζ<0,sinβζ>0,则Cζ,1min(xζ,1min,yζ,1min)为第ζ个平行四边形编码标志中第1编码区域的特征角点,把Cζ,1min(xζ,1min,yζ,1min)的亚像素坐标值赋值给Cζ,1(xζ,1,yζ,1);Cζ,2min(xζ,2min,yζ,2min)为第ζ个平行四边形编码标志中第6编码区域的特征角点,把Cζ,2min(xζ,2min,yζ,2min)的亚像素坐标值赋值给Cζ,4(xζ,4,yζ,4);
(2)若sinαζ>0,sinβζ<0,则Cζ,2min(xζ,2min,yζ,2min)为第ζ个平行四边形编码标志中第1编码区域的特征角点,把Cζ,2min(xζ,2min,yζ,2min)的像素坐标值赋值给Cζ,1(xζ,1,yζ,1);Cζ,1min(xζ,1min,yζ,1min)为第ζ个平行四边形编码标志中第6编码区域的特征角点,把Cζ,1min(xζ,1min,yζ,1min)的像素坐标值赋值给Cζ,4(xζ,4,yζ,4);
步骤15、在第ζ组第2个备份二值化图像P′ζ,2上,通过步骤14已经在4个像素点C″′ζ,1(x″′ζ,1,y″′ζ,1)、C″′ζ,2(x″′ζ,2,y″′ζ,2)、C″′ζ,3(x″′ζ,3,y″′ζ,3)、C″′ζ,4(x″′ζ,4,y″′ζ,4)中找到第ζ个平行四边形编码标志中第1编码区域和第6编码区域的特征角点的亚像素坐标值并分别赋值给了Cζ,1(xζ,1,yζ,1)和Cζ,4(xζ,4,yζ,4),将余下的2个像素点的亚像素坐标值分别赋值给第ζ个平行四边形编码标志的临时坐标值1,记为C′ζ,5(x′ζ,5,y′ζ,5),以及临时坐标值2,记为C′ζ,6(x′ζ,6,y′ζ,6);
根据sinωζ和sinξζ的取值情况,将C′ζ,5(x′ζ,5,y′ζ,5)或C′ζ,6(x′ζ,6,y′ζ,6)的亚像素坐标值赋值给Cζ,2(xζ,2,yζ,2)或Cζ,3(xζ,3,yζ,3);具体分为以下两种情况执行:
(1)若sinωζ==0,sinξζ≠0,则第1临时坐标C′ζ,5(x′ζ,5,y′ζ,5)即为第ζ个平行四边形编码标志中第3编码区域的特征角点的亚像素坐标,把C′ζ,5(x′ζ,5,y′ζ,5)的亚像素坐标值赋值给Cζ,2(xζ,2,yζ,2);第2临时坐标C′ζ,6(x′ζ,6,y′ζ,6)即为第ζ个平行四边形编码标志中第4编码区域的特征角点的亚像素坐标,把C′ζ,6(x′ζ,6,y′ζ,6)的亚像素坐标值赋值给Cζ,3(xζ,3,yζ,3);
(2)若sinωζ≠0,sinξζ==0,则第2临时坐标C′ζ,6(x′ζ,6,y′ζ,6)即为第ζ个平行四边形编码标志中第3编码区域的特征角点的亚像素坐标,把C′ζ,6(x′ζ,6,y′ζ,6)的亚像素坐标值赋值给Cζ,2(xζ,2,yζ,2);第1临时坐标C′ζ,5(x′ζ,5,y′ζ,5)即为第ζ个平行四边形编码标志中第4编码区域的特征角点的亚像素坐标,把C′ζ,5(x′ζ,5,y′ζ,5)的亚像素坐标值赋值给Cζ,3(xζ,3,yζ,3);
至此,在第ζ组第2个备份二值化图像P′ζ,2上找到了第ζ个平行四边形编码标志中第1编码区域的特征角点的亚像素坐标Cζ,1(xζ,1,yζ,1)、第3编码区域的特征角点的亚像素坐标Cζ,2(xζ,2,yζ,2)、第4编码区域的特征角点的亚像素坐标Cζ,3(xζ,3,yζ,3)和第6编码区域的特征角点的亚像素坐标Cζ,4(xζ,4,yζ,4);
在第ζ组第2个备份二值化图像P′ζ,2上,以第ζ个平行四边形编码标志中的定位圆质心o′l,ξ(x′l,ζ,y′l,ζ)为起点做与第5判断向量平行且同向的单位向量,记为以第ζ个平行四边形编码标志中定向圆环质心o′d,ζ(x′d,ζ,y′d,ζ)为起点做与第5判断向量平行且同向的单位向量,记为
步骤18、定义6个浮点型二维数组 用于存放第ζ个平行四边形编码标志中分别位于第1编码区域、第2编码区域、第3编码区域、第4编码区域、第5编码区域和第6编码区域的编码标志圆轮廓质心在第ζ组第2个备份二值化图像P′ζ,2上的像素坐标,初始化这6个二维数组中的所有元素,且赋值为-1;另外取6个整数变量并将其初始化,取整数变量i并赋值i=1;
步骤19、在第ζ组第2个备份二值化图像P′ζ,2上,根据式(12)、式(13)计算出第ζ个平行四边形编码标志中编码标志圆轮廓的质心像素坐标分别与定位圆质心o′l,ζ和定向环质心o′d,ζ所形成的第i组第1个象限向量和第i组第2个象限向量
步骤21.1:
步骤21.2、本步骤分为以下情况执行:
步骤21.3、本步骤分为以下情况执行:
Wζ=Vζ T·U (18)
至此获得了第ζ个平行四边形编码标志中4个特征角点各自的编码号与其相应的亚像素坐标值:
步骤24、判断ζ是否小于N,若满足,则将ζ+1重新赋值给ζ,并返回步骤6.1顺序执行;若不满足,则已经找到每个平行四边形编码标志上4个特征角点的亚像素坐标值和对应的编码号。
进一步的,步骤6.1中,以获得的第ζ组第1个备份二值化图像P′ζ,1以及圆环黑色连通域质心集合A中第ζ个圆环黑色连通域质心像素坐标值(x″′d,ζ,y″′d,ζ)作为输入条件,利用去除编码标志复杂背景算法去除第ζ组第1个备份二值化图像P′ζ,1的复杂背景,并最终获得第ζ个无复杂背景平行四边形二值化图像P″ζ的步骤如下:
步骤6.1.1、取整数变量i并赋值i=1;
步骤6.1.2、在第ζ组第1个备份二值化图像P′ζ,1上寻找距离像素坐标值(x″′d,ζ,y″′d,ζ)最近的白色连通域并记为第ζ组第1个备份二值化图像P′ζ,1上的第i个圆环白色连通域Μ′i,并将第i个圆环白色连通域Μ′i的灰度值均赋值为0;
步骤6.1.3、判断i是否小于N,若i<N,则将i+1赋值给i,返回步骤1.2顺序执行;否则,将第ζ组第1个备份二值化图像P′ζ,1中所有平行四边形编码标志中的定向圆环置为黑色,并将第ζ组第1个备份二值化图像P′ζ,1经过本步骤处理后的图像记为第ζ个黑色背景提取图像P′ζ,b,并执行步骤6.1.4;
步骤6.1.4、在第ζ个黑色背景提取图像P′ζ,b上的所有黑色连通域中,寻找距离像素坐标值(x″′d,ζ,y″′d,ζ)最近的黑色连通域记为第ζ个黑色背景提取图像P′ζ,b上的第ζ个平行四边形背景黑色连通域Ωζ;
步骤6.1.5、在第ζ个黑色背景提取图像P′ζ,b上提取第ζ个平行四边形背景黑色连通域Ωζ的所有内、外轮廓,获得第ζ个平行四边形背景黑色连通域Ωζ的内、外轮廓集合Vζ;
步骤6.1.6、比较轮廓集合Vζ中所有轮廓的周长,将周长最长的轮廓记为第ζ个平行四边形背景黑色连通域Ωζ的外轮廓Εζ;
步骤6.1.7、在第ζ个黑色背景提取图像P′ζ,b上,将外轮廓Εζ以内的所有像素点的灰度值赋值为255,外轮廓Εζ以外的所有像素点的灰度值赋值为0,从而获得第ζ个无复杂背景平行四边形二值化图像P″ζ;使得在第ζ个无复杂背景平行四边形二值化图像P″ζ中,包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志内部的所有像素点灰度值被置为255(即为白色),包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志外部的所有像素点灰度值被置为0(即为黑色)。
还提供了一种计算机可读存储介质,包括与具有图像处理功能的电子设备结合使用的计算机程序,所述计算机程序可被处理器执行以所述的解码方法。
与现有技术相比较,本发明的有益效果如下:
1、本发明提出的平行四边形编码标志中的编码图案,由定位图案、定向图案和编码标志图案构成,其中,定位图案和定向图案能够实现在相机标定、三维数据拼接等过程中,自动判断所测目标与摄像机相对旋转方向,提高了灵活性;
2、本发明提出的平行四边形编码标志中,分为6个编码区域,每个编码区域包含2个编码标志圆,因此本发明提出的平行四边形编码标志共可产生4096种不同的编码,容量得以扩充,因此,针对尺寸较大的目标也能完成三维数据拼接;
3、本发明提出的一种基于区域分割的平行四边形编码标志的解码方法,能够有效地去除噪声的干扰,提出的去除复杂背景的方法极大地提高了解码的准确性;
4、本发明提出的一种基于区域分割的平行四边形编码标志的解码方法具有高效性、解码速度快的优点,且可以实现实时解码。
附图说明
图1为本发明的平行四边形编码标志的示意图;
图3编码标志灰度图像P1的展示图;
图4为编码标志二值化图像P2的展示图;
图5为在编码标志二值化图像P2上提取圆环黑色连通域质心集合A的展示图;
图6为第1个无复杂背景平行四边形二值化图像P″1,1的展示图;
图7为利用Harris检点检测算法提取第1个无复杂背景平行四边形二值化图像P1″中角点的亚像素坐标结果展示图;
图8为第1个覆盖结果图像T1的展示图;
图9为利用Harris检点检测算法提取第1个覆盖结果图像T1中角点的亚像素坐标结果展示图;
图10为在备份二值化图像P′1,2上将角点四边形S1以外的所有像素点的灰度值赋值为255后所得图像的展示图;
图11为本发明所述解码方法的流程示意图。
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述,以使本发明的优点和特征能更易于被本领域技术人员理解,从而对本发明的保护范围做出更为清楚明确的界定。
请参阅图1,一种基于区域分割的平行四边形编码标志,以平行四边形编码标志的4个顶点作为特征角点,在视觉结构光三维测量过程中,用于目标曲面的三维点云拼接;所述平行四边形编码标志为长度为a且宽度为b的平行四边形,a和b均大于零;本实施例中,平行四边形编码标志为长度和宽度均为24毫米的正方形。
在平行四边形编码标志背景内包括定位图案、定向图案和编码标志图案,其中编码标志图案又由多个编码单元图案组成。由定向图案和定位图案实现平行四边形编码标志旋转方向的判断,编码标志图案用于对平行四边形编码标志中的每一个角点进行编码。
所述的一种基于区域分割编码的平行四边形编码标志,所述平行四边形编码标志中的定位图案、定向图案和编码标志图案均不重叠且不连通;
所述平行四边形编码标志中的定位图案、定向图案和编码标志图案均不重叠且不连通;定位图案为圆形连通域,定向图案为圆环连通域,所有编码单元图案为圆形连通域。本实施例中,每一个平行四边形编码单元中的定位图案为直径6毫米的实心圆,定向图案均为外直径6毫米且内直径3毫米的圆环,编码单元图案均为直径3毫米的实心圆。在所述平行四边形编码标志中,每个编码单元图案的轮廓长度均小于定位图案的轮廓长度,本实施例中,定位图案的周长小于96毫米,每个编码单元图案的轮廓长度9.42毫米小于定位图案的轮廓长度18.85毫米;
在所述平行四边形编码标志中,将所有编码单元图案中面积最大的编码单元图案的面积记为编码单元图案最大面积,并且定向图案、定位图案和所有编码单元图案的面积大小满足如下关系:定向图案的面积大于定位图案和编码单元图案最大面积,各个编码单元图案的面积可以互不相同。
本实施例中,记平行四边形编码标志的背景颜色为黑色(即颜色1),定位图案和定向图案的颜色为白色(即颜色2),而所有编码单元图案的颜色与黑色(颜色1)相同或与白色(颜色2)相同,并且颜色1与颜色2具有明显差异;编码单元图案的颜色根据平行四边形编码标志的编码规则可以为白色(即颜色2)也可以为黑色(即颜色1)。
所述平行四边形编码标志中,定位图案的质心和定向图案的质心的连线的中点在平行四边形编码标志的中心位置附近处;
如图2所示,任取平行四边形编码标志的一个特征角点记为向量确定第一顶点o″1,在平行四边形编码标志上将相交形成向量确定第一顶点o″1的任意一条边记为向量确定第一边Νv1,在向量确定第一边Νv1上取平行四边形编码标志的顶点记为向量确定第二点o″2,其中向量确定第二点o″2与向量确定第一顶点o″1是互不重合的2个点,记向量为规定向量并且平行四边形编码标志内的定位图案和定向图案的位置关系如下:平行四边形编码标志内由定向图案质心指向定位图案质心的方向与规定向量的方向相同;
将平行四边形编码标志在的平面记为标志平面Pt,以向量确定第一顶点o″1为起点做一个与规定向量同向的单位向量记为第1个规定单位向量当人正视看向平行四边形编码标志时,以向量确定第一顶点o″1为旋转中心,在标志平面Pt内将第1个规定单位向量逆时针旋转β′角度(0°<β′<90°)得到第2个规定单位向量在空间中以向量确定第一顶点o″1为起点做一个与所得向量的方向相同的单位向量,并记为正向向量将距离平行四边形编码标志中定向图案最近的两个特征角点分别记为第1临时顶点o″3和第2临时顶点o″4;若所得向量的方向与正向向量的方向相同,则将记为向量辅助向量若所得向量的方向与正向向量的方向相同,则将向量记为辅助向量
参阅图11,利用数字图像处理的方法获得摄像机拍摄的包含21个平行四边形编码标志的图像中每个特征角点的编码序号和亚像素坐标,完成平行四边形编码标志解码的步骤如下:
步骤1.1、利用摄像机拍摄空间中放置的21个平行四边形编码标志,进而获得编码标志图像,所述编码标志图像包含21个平行四边形编码标志;
步骤1.2、建立角点的像素坐标系:
在上述拍摄的包含21个平行四边形编码标志的图像中,取图像左上角作为角点像素坐标系的原点o,自左向右为角点像素坐标系的x轴方向,自上向下作为角点像素坐标系的y轴方向,由此建立角点像素坐标系o-xy;
步骤2.1、对编码标志图像进行灰度处理,得到编码标志灰度图像P1;其中,编码标志灰度图像P1为8位灰度图,如图3所示;
步骤2.2、对编码标志灰度图像P1进行复制备份,得到编码标志灰度备份图像P1′;
步骤2.3、对编码标志灰度备份图像P1′进行二值化处理,使得每个平行四边形编码标志的背景颜色均为黑色(即灰度值为0),每个平行四边形编码标志内的定向图案和定位图案均为白色(即灰度值为255),编码单元图案的颜色根据平行四边形编码标志的编码规则可以为白色也可以为黑色,进而得到编码标志二值化图像P2,如图4所示;
步骤3、设定黑色连通域圆度阀值λ′=1,像素点个数阀值λ″=30;
步骤4.1、计算编码标志二值化图像P2中所有黑色连通域的圆度值λn(其中n=1,2,3,4...)和像素点个数cn(其中n=1,2,3,4...),其中每个黑色连通域的圆度值可由式(1)得出,
其中ln代表对应黑色连通域轮廓的周长,sn代表对应黑色连通域的面积;
步骤4.2、本实施例中选择的黑色连通域圆度阀值λ′=1和像素点个数阀值λ″=30可以在编码标志二值化图像P2中获得21个平行四边形编码标志中各自所对应的21个定向圆环中心的黑色连通域;将21个定向圆环中心的黑色连通域分别记为编码标志二值化图像P2中的圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域Μ21,并将圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域Μ21依次放入圆环黑色连通域集合A′中,即圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域Μ21分别为圆环黑色连通域集合A′中的第1个元素、第2个元素、…、第21个元素;由此获得编码标志二值化图像P2中的每个平行四边形编码标志中定向圆环内部的黑色连通域;其中编码标志二值化图像P2中的圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域Μ21的圆度值均小于1,并且圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域Μ21的像素点个数均大于30;
步骤4.3、取整数变量i并赋值i=1;
步骤4.4、计算圆环黑色连通域集合A′中第i个圆环黑色连通域Μi的质心像素坐标并记为o″d,i(x″d,i,y″d,i),并且将得到第i个圆环黑色连通域Μi的质心像素坐标o′d,i(x′d,i,y′d,i)放入圆环黑色连通域质心集合A中,作为圆环黑色连通域质心集合A的第i元素;
步骤4.5、判断i是否小于21,若i<21,则将i+1赋值给i,返回步骤4.4顺序执行;否则,如图3所示,即获得了编码标志二值化图像P2中的圆环黑色连通域Μ1、圆环黑色连通域Μ2、圆环黑色连通域Μ3、…、圆环黑色连通域Μ21的质心像素坐标o″d,1(437,213)、o″d,2(x″d,2,y″d,2)、…、o″d,21(x″d,21,y″d,21),并依次放入圆环黑色连通域质心集合A中,如图5所示;
步骤5、取整数变量ζ并赋值ζ=1;
步骤6.1、将编码标志二值化图像P2进行两次复制备份,分别得到第1组第1个备份二值化图像P′1,1和第1组第2个备份二值化图像P′1,2;
步骤6.2、利用去除编码标志复杂背景算法对第1组第1个备份二值化图像P′1,1进行数字图像处理,进而得到第1个无复杂背景平行四边形二值化图像P″1,1;在所述第1个无复杂背景平行四边形二值化图像P′1,1中,包含质心像素坐标值(437,213)的平行四边形编码标志内部的所有像素点的灰度值被置为255(即白色),包含质心像素坐标值(437,213)的平行四边形编码标志外部的所有像素点的灰度值值被置为0(即黑色),如图6所示;
步骤7.1、对编码标志灰度图像P1再次进行复制备份,获得第1个备份灰度图像P1″′;
步骤7.2、利用Harris角点检测算法在第1个无复杂背景平行四边形二值化图像P″1,1中提取4个角点的亚像素坐标,如图7所示,并记为第1个无复杂背景平行四边形二值化图像P″1,1上的4个备选角点C″1,1(480.1343,153.2671)、C″1,2(480.2826,236.5667)、C″1,3(392.6782,236.2432)、C″1,4(392.2600,152.8951);
步骤7.3、在第1个无复杂背景平行四边形二值化图像P″1,1中做黑色连通域腐蚀,使得第1个无复杂背景平行四边形二值化图像P″1,1中的白色连通域扩大,在后续进行图像乘法运算时,能够保证运算结果图像中提取出完整的包含质心像素坐标值(437,213)的平行四边形编码标志的灰度图案且其4个特征角点保存完整;并将第1个无复杂背景平行四边形二值化图像P″1,1经本步骤处理所得的图像记为第1个运算二值化图像P″1,2;
步骤7.4、将第1个运算二值化图像P″1,2中灰度值为255的部分的值置为1,将第1个运算二值化图像P″1,2中灰度值为0的部分的值置为0,而后第1个运算二值化图像P″1,2和第1个备份灰度图像P1″′进行图像乘法运算,将得到的结果图像记为第1个覆盖结果图像T1,如图8所示;所述第1个覆盖结果图像T1中只含有包含质心像素坐标值(437,213)的平行四边形编码标志的图像信息,去除了其它平行四边形编码标志的图像信息;
步骤7.5、利用Harris角点检测算法在第1个覆盖结果图像T1中提取所有角点的亚像素坐标并存入第1个初选角点亚像素坐标集合B1中,如图9所示;其中,第1个初选角点亚像素坐标集合B1中存放了第1个覆盖结果图像T1上包含质心像素坐标值(437,213)的平行四边形编码标志的4个特征角点亚像素坐标值和干扰角点的亚像素坐标值,所述干扰角点是除特征角点以外的干扰角点;
步骤7.6、取整数变量i并赋值i=1;
步骤7.7、在所述第1个初选角点亚像素坐标集合B1中寻找距离亚像素坐标C″1,i(x″1,i,y″1,i)最近的角点亚像素坐标并记为C′1,i(x′1,i,y′1,i);
步骤7.8、判断i是否小于4,若i<4,则将i+1赋值给i并返回执行步骤7.7;否则此时已经获得了第1个覆盖结果图像T1上包含质心像素坐标值(437,213)的平行四边形编码标志的4个特征角点的亚像素坐标C′1,1(477.4545,155.4545)、C′1,2(477.3571,238.3571)、C′1,3(394.4545,238.4545)、C′1,4(394.3000,155.3000),且执行步骤8.1;
步骤8.1、在第1组第2个备份二值化图像P′1,2上,将包含质心像素坐标值(437,213)的平行四边形编码标志记为第ζ个平行四边形编码标志;
步骤8.2、在第1组第2个备份二值化图像P′1,2上,选取亚像素坐标值分别为(480.1343,153.2671)、(480.2826,236.5667)、(392.6782,236.2432)、(392.2600,152.8951)的4个像素点作为第1个平行四边形编码标志的角点四边形S1的4个顶点,将4个顶点相连进而获得第1个平行四边形编码标志的角点四边形S1;
步骤9、在第1组第2个备份二值化图像P′1,2上,寻找距离上述圆环黑色连通域质心集合A中第1个定向圆环质心像素坐标值(437,213)最近的白色连通域,并将此白色连通域的灰度值赋值为0;
步骤10、在第1组第2个备份二值化图像P′1,2上,保持第1个平行四边形编码标志的角点四边形S1内部的所有像素点灰度值不变,第1个平行四边形编码标志的角点四边形S1以外的所有像素点灰度值赋值为255,如图10所示;
步骤11.1、在第1组第2个备份二值化图像P′1,2上,提取第1个平行四边形编码标志的角点四边形S1中所有的内、外轮廓,记为第1个平行四边形编码标志内的轮廓集合D1;
步骤11.2统计第1个平行四边形编码标志内的轮廓集合D1内每个轮廓包含的像素点数,将包含像素点数第二多的轮廓记为第1组第2个备份二值化图像P′1,2上第1个平行四边形编码标志内的定位圆轮廓G1,计算其质心像素坐标并记为第1组第2个备份二值化图像P′1,2上第1个平行四边形编码标志内的定位圆质心像素坐标o′l,1(437,181);
步骤12.1、根据第1个平行四边形编码标志内的轮廓集合D1中,除去包含像素点数最多的两个轮廓。本实施例中剩余轮廓个数κ1==2,即κ1≠0,则这2个轮廓即为第1个平行四边形编码标志内编码标志圆的轮廓,分别记为编码标志圆编码标志圆
步骤12.2、给整数变量i并赋予初值i=1;
步骤12.3、在第1组第2个备份二值化图像P′1,2上,计算编码标志圆轮廓的质心像素坐标o′1 i(x′1 i,y′1 i),将i+1重新赋值给i后继续执行此步骤,直到i>2结束;由此可以得到对应编码标志圆轮廓编码标志圆轮廓的质心像素坐标o′1 1(459,176)、o′1 2(459,190);
步骤13.1、在第1组第2个备份二值化图像P′1,2上,记在第1个平行四边形编码标志中位于第1编码区域、第3编码区域、第4编码区域和第6编码区域的特征角点的亚像素坐标分别为C1,1(x1,1,y1,1)、C1,2(x1,2,y1,2)、C1,3(x1,3,y1,3)、C1,4(x1,4,y1,4);并在第1组第2个备份二值化图像P′1,2上,将像素坐标值为(437,213)的像素点记为第1组第2个备份二值化图像P′1,2上第1个平行四边形编码标志的定向圆环质心o′d,1(437,213),将亚像素坐标值为(477.4545,155.4545)、(477.3571,238.3571)、(394.4545,238.4545)、(394.3000,155.3000)的4个像素点分别记为C″′1,1(477.4545,155.4545)、C″′1,2(477.4571,238.3571)、C″′1,3(394.4545,238.4545)、C″′1,4(394.3000,155.3000);
步骤14、在第1组第2个备份二值化图像P′1,2上,找出4个像素点C″′1,1(477.4545,155.4545)、C″′1,2(477.4571,238.3571)、C″′1,3(394.4545,238.4545)、C″′1,4(394.3000,155.3000)中距离定位圆质心像素坐标o′l,1(437,181)最近的2个像素点,并分别记为C1,1min(x1,1min,y1,1min)和C1,2min(x1,2min,y1,2min)。本实施例中,即为C1,1min(477.4545,155.4545)和C1,2min(394,3000,155.3000);通过公式(3)和(4)计算出在第1个平行四边形编码标志中的第1判断向量和第2判断向量并通过式(5)和式(6)计算出第1区域划分正弦值sinα1和第2区域划分正弦值sinβ1;
本实施例中,由于sinα1<0,sinβ1>0,则C1,1min(477.4545,155.4545)为第1个平行四边形编码标志中第1编码区域的特征角点,把C1,1min(477.4545,155.4545)的亚像素坐标值赋值给C1,1(x1,1,y1,1);C1,2min(394,3000,155.3000)为第1个平行四边形编码标志中第6编码区域的特征角点,把C1,2min(394,3000,155.3000)的亚像素坐标值赋值给C1,4(x1,4,y1,4);
步骤15、在上述步骤14中,已在4个像素点C″′1,1(477.4545,155.4545)、C″′1,2(477.4571,238.3571)、C″′1,3(394.4545,238.4545)、C″′1,4(394.3000,155.3000)中找到第1个平行四边形编码标志中第1编码区域和第6编码区域的特征角点C1,1(477.4545,155.4545)和C1,4(394.3000,155.3000);将余下的2个像素点的亚像素坐标分别赋值给第1个平行四边形编码标志的第1临时坐标值,记为C′1,5(477.3571,238.3571),以及第2临时坐标值,记为C′1,6(394.4545,238.4545);根据公式(7)和(8)可以求出在第1个平行四边形编码标志中的第3判断向量和第4判断向量
预先设定平行向量正弦阀值g=0.01,若sinω1<g则判定认为第3判断向量和方向向量平行,在本实施例中,sinω1≈0,sinξ1≠0,则C′1,5(477.3571,238.3571)即为第1个平行四边形编码标志中第3编码区域的特征角点的亚像素坐标,把C′1,5(477.3571,238.3571)的亚像素坐标值赋值给C1,2(x1,2,y1,2);C′1,6(394.4545,238.4545)为第1个平行四边形编码标志中第4编码区域的特征角点的亚像素坐标,把C′1,6(394.4545,238.4545)的亚像素坐标值赋值给C1,3(x1,3,y1,3);
至此,在第1组第2个备份二值化图像P′1,2上找到了第1个平行四边形编码标志中第1编码区域的特征角点的亚像素坐标C1,1(477.4545,155.4545)、第3编码区域的特征角点的亚像素坐标C1,2(477.3571,238.3571)、第4编码区域的特征角点的亚像素坐标C1,3(394.4545,238.4545)和第6编码区域的特征角点的亚像素坐标C1,4(394.3000,155.3000);
步骤17、在第1组第2个备份二值化图像P1′,2上,根据得出的第1个平行四边形编码标志中第1编码区域的角点C1,1(477.4545,155.4545),第6编码区域的角点C1,4(394.3000,155.3000),可由公式(11)得出第1个平行四边形编码标志中的第5判断向量
过第1个平行四边形编码标志中的定位圆质心o′l,1(437,181)做与第5判断向量平行且同向的单位向量,记为过第1个平行四边形编码标志中定向圆环质心o′d,1(437,213)做与第5判断向量平行且同向的单位向量,记为将整数变量i重新赋值i=1;
步骤18、定义6个浮点型二维数组 用于存放第1个平行四边形编码标志中分别位于第1编码区域、第2编码区域、第3编码区域、第4编码区域、第5编码区域和第6编码区域的编码标志圆轮廓质心在第1组第2个备份二值化图像P′1,2上的像素坐标,初始化这6个二维数组中的所有元素,且赋值为-1;另外取6个整数变量并将其初始化,
步骤19、计算步骤15.2中,第1个平行四边形编码标志中编码标志圆轮廓的质心像素坐标o′1 i(x′1 i,y′1 i)分别与定位圆中心o′l,1和定向环中心o′d,1所形成的第i组第1判断向量和第i组第2判断向量本实施例中可分别求得结果如下:
本实施例中,编码标志圆所属编码区域的判断结果如下:
步骤21.1、在本实施例中,根据本步骤可得:
步骤21.2、在本实施例中,根据本步骤可得:
步骤21.3、在本实施例中,根据步骤可得:
W1=V1 T·U=10 (18)
其中,列向量U=(20,21,22,...211)T,列向量V1=(0,1,0,1,0,0,0,0,0,0,0,0)T;
在本实施例中可得:
步骤24、判断ζ是否小于21,若满足,则将ζ+1重新赋值给ζ,并返回步骤6顺序执行;若不满足,则已经找到每个平行四边形编码标志上4个特征角点的亚像素坐标值和对应的编码号。
将步骤6.1获得第1组第1个备份二值化图像P′1,1以及圆环黑色连通域质心集合A中第1个圆环黑色连通域质心像素坐标值(437,213)作为输入,利用去除编码标志复杂背景算法去除第1组第1个备份二值化图像P′1,1中的复杂背景,并最终获得第1个无复杂背景平行四边形二值化图像P1″的步骤如下:
步骤6.1.1、取整数变量i并赋值i=1;
步骤6.1.2、在第1组第1个备份二值化图像P′1,1上寻找距离像素坐标值(437,213)最近的白色连通域并记为第1组第1个备份二值化图像P′1,1上第i个圆环白色连通域Μ′i,并将第1组第1个备份二值化图像P′1,1上第i个圆环白色连通域Μ′i的灰度值均赋值为0;
步骤6.1.3、判断i是否小于21,若i<21,则将i+1赋值给i,返回步骤1.2顺序执行;否则,已将第1组第1个备份二值化图像P′1,1中所有平行四边形编码标志中的定向圆环置为黑色,并将第1组第1个备份二值化图像P′1,1经过本步骤1.3处理后所得的图像记为第1个黑色背景提取图像P′1,b,并执行步骤6.1.4;
步骤6.1.4、在第1个黑色背景提取图像P′1,b上的所有黑色连通域中,寻找距离像素坐标值(437,213)最近的黑色连通域记为第1个黑色背景提取图像P′1,b上的第1个平行四边形背景黑色连通域Ω1;
步骤6.1.5、在第1个黑色背景提取图像P′1,b上提取第1个平行四边形背景黑色连通域Ω1的所有内、外轮廓,获得第1个平行四边形背景黑色连通域Ω1的内、外轮廓集合V1;
步骤6.1.6、比较轮廓集合V1中所有轮廓的周长,将周长最长的轮廓记为第1个黑色背景提取图像P′1,b上第1个平行四边形背景黑色连通域Ω1的外轮廓Ε1;
步骤6.1.7、在第1个黑色背景提取图像P′1,b上利用射线法区分平行四边形背景黑色连通域Ω1的外轮廓Ε1以外的所有像素点和外轮廓Ε1以内的所有像素点,并将外轮廓Ε1以内的所有像素点灰度值赋值为255,外轮廓Ε1以外的所有像素点灰度值赋值为0,进而获得第1个无复杂背景平行四边形二值化图像P1″;使得第1个无复杂背景平行四边形二值化图像P1″中,包含质心像素坐标值(437,213)的平行四边形编码标志内部的所有像素点灰度值被置为255,包含质心像素坐标值(437,213)的个平行四边形编码标志外部的所有像素点灰度值被置为0,如图10所示。
本发明提供的基于区域分割的平行四边形编码标志的解码方法,需要编制相应的计算机程序,并在计算机上执行程序以实现相应的运算处理及逻辑控制功能,因而本发明也提供一种计算机可读存储介质,包括与具有图像处理功能的电子设备结合使用的计算机程序,所述计算机程序可被处理器执行以所述的解码方法。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (7)
1.一种基于区域分割的平行四边形编码标志的解码方法,以所述平行四边形编码标志的四个顶点作为特征角点,在视觉结构光三维测量过程中,用于目标曲面的三维点云拼接;编码方格的表面设置有平行四边形背景图案和编码图案,所述编码图案位于平行四边形背景图案的内部,所述编码图案包括定位图案、定向图案和编码标志图案,所述定向图案和定位图案用于平行四边形编码标志方向的判断,所述编码标志图案用于平行四边形编码标志各个角点的编码,其特征在于:利用数字图像处理的方法获得摄像机拍摄的包含N个平行四边形编码标志的图像中每个特征角点的编码序号和亚像素坐标,进而完成平行四边形编码标志解码的步骤如下:
步骤1.1、利用摄像机拍摄空间中放置的N个平行四边形编码标志,获得编码标志图像,所述编码标志图像当中包含N个平行四边形编码标志;
步骤1.2、建立角点的像素坐标系;
步骤2.1、对编码标志图像进行8位灰度处理,得到编码标志灰度图像P1;
步骤2.2、对编码标志灰度图像P1进行复制备份,得到编码标志灰度备份图像P1′;
步骤2.3、对编码标志灰度备份图像P1′进行二值化处理,得到编码标志二值化图像P2;
步骤3、设定黑色连通域圆度阀值λ′,像素点个数阀值λ″;
步骤4.1、计算编码标志二值化图像P2中所有黑色连通域的圆度值λn和像素点个数cn,其中每个黑色连通域的圆度值由式(1)得出,
其中,ln代表对应黑色连通域轮廓的周长,sn代表对应黑色连通域的面积,n=1,2,3,4…;
步骤4.2、选择合适的黑色连通域圆度阀值λ′和像素点个数阀值λ″以在编码标志二值化图像P2中获得N个平行四边形编码标志中各自所对应的N个定向圆环中心的黑色连通域,并放入圆环黑色连通域集合A′中;
其中,编码标志二值化图像P2中的各个圆环黑色连通域的圆度值均小于圆度阀值λ′,且各个圆环黑色连通域的像素点个数均大于像素点个数阀值λ″;
步骤4.3、取整数变量i并赋值i=1;
步骤4.4、计算圆环黑色连通域集合A′中的第i个圆环黑色连通域Μi的质心像素坐标并记为o″d,i(x″d,i,y″d,i),将得到的第i个圆环黑色连通域Μi的质心像素坐标o″d,i(x″d,i,y″d,i)作为圆环黑色连通域质心集合A的第i元素;
步骤4.5、判断i是否小于N,若i<N,则将i+1赋值给i,返回步骤4.4顺序执行;否则,即获得了编码标志二值化图像P2中各个圆环黑色连通域的质心像素坐标o″d,1(x″d,1,y″d,1)、o″d,2(x″d,2,y″d,2)、…、o″d,N(x″d,N,y″d,N),并依次放入圆环黑色连通域质心集合A中;
步骤5、取整数变量ζ并赋值ζ=1;
步骤6.1、将编码标志二值化图像P2进行两次复制备份,分别得到第ζ组第1个备份二值化图像P′ζ,1和第ζ组第2个备份二值化图像P′ζ,2;
步骤6.2、利用去除编码标志复杂背景算法对第ζ组第1个备份二值化图像P′ζ,1进行数字图像处理,得到第ζ个无复杂背景平行四边形二值化图像P″ζ,1;
步骤7.1、对编码标志灰度图像P1再次进行复制备份,获得第ζ个备份灰度图像P″′ζ;
步骤7.2、利用角点检测算法在第ζ个无复杂背景平行四边形二值化图像P″ζ,1中提取4个角点的亚像素坐标,并记为第ζ个无复杂背景平行四边形二值化图像P″ζ,1上的4个备选角点C″ζ,1(x″ζ,1,y″ζ,1)、C″ζ,2(x″ζ,2,y″ζ,2)、C″ζ,3(x″ζ,3,y″ζ,3)、C″ζ,4(x″ζ,4,y″ζ,4);
步骤7.3、在第ζ个无复杂背景平行四边形二值化图像P″ζ,1中做黑色连通域腐蚀,使得第ζ个无复杂背景平行四边形二值化图像P″ζ,1中的白色连通域扩大;并将处理所得的图像记为第ζ个运算二值化图像P″ζ,2;
步骤7.4、将第ζ个运算二值化图像P″ζ,2中灰度值为255的部分的值置为1,将第ζ个运算二值化图像P″ζ,2中灰度值为0的部分的值置为0,而后利用第ζ个运算二值化图像P″ζ,2和第ζ个备份灰度图像P″′ζ进行乘法运算,将得到的结果图像记为第ζ个覆盖结果图像Tζ;
步骤7.5、利用角点检测算法在第ζ个覆盖结果图像Tζ中提取所有角点的亚像素坐标并存入第ζ个初选角点亚像素坐标集合Bζ中;
步骤7.6、取整数变量i并赋值i=1;
步骤7.7、在所述第ζ个初选角点亚像素坐标集合Bζ中寻找距离亚像素坐标值(x″ζ,i,y″ζ,i)最近的角点亚像素坐标并记为C′ζ,i(x′ζ,i,y′ζ,i);
步骤7.8、判断i是否小于4,若i<4,则将i+1赋值给i并返回执行步骤7.7;否则,此时已经获得了第ζ个覆盖结果图像Tζ上包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志的4个特征角点的亚像素坐标C′ζ,1(x′ζ,1,y′ζ,1)、C′ζ,2(x′ζ,2,y′ζ,2)、C′ζ,3(x′ζ,3,y′ζ,3)、C′ζ,4(x′ζ,4,y′ζ,4),且执行步骤8.1;
步骤8.1、在第ζ组第2个备份二值化图像P′ζ,2上,将包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志记为第ζ个平行四边形编码标志;
步骤8.2、在第ζ组第2个备份二值化图像P′ζ,2上,选取亚像素坐标值分别为(x′ζ,1,y′ζ,1)、(x′ζ,2,y′ζ,2)、(x′ζ,3,y′ζ,3)、(x′ζ,4,y′ζ,4)的4个像素点作为第ζ个平行四边形编码标志的角点平行四边形Sζ的4个顶点,将4个顶点相连进而获得第ζ个平行四边形编码标志的角点四边形Sζ;
步骤9、在第ζ组第2个备份二值化图像P′ζ,2上,寻找距离质心像素坐标值(x″d,ζ,y″d,ζ)最近的白色连通域,并将此白色连通域的灰度值赋值为0;
步骤10、在第ζ组第2个备份二值化图像P′ζ,2上,第ζ个平行四边形编码标志的角点四边形Sζ内部的所有像素点的灰度值保持不变,将第ζ个平行四边形编码标志的角点四边形Sζ以外的所有像素点的灰度值赋值为255;
步骤11.1、在第ζ组第2个备份二值化图像P′ζ,2上,提取第ζ个平行四边形编码标志的角点四边形Sζ中所有的内、外轮廓,记为第ζ个平行四边形编码标志内的轮廓集合Dζ;
步骤11.2、统计第ζ个平行四边形编码标志内的轮廓集合Dζ内每个轮廓包含的像素点数,将包含像素点数第二多的轮廓记为第ζ组第2个备份二值化图像P′ζ,2上第ζ个平行四边形编码标志内的定位圆轮廓Gζ,计算其质心像素坐标并记为第ζ组第2个备份二值化图像P′ζ,2上第ζ个平行四边形编码标志内的定位圆质心像素坐标o′l,ζ(x′l,ζ,y′l,ζ);
步骤12.1、根据第ζ个平行四边形编码标志内的轮廓集合Dζ中,除去包含像素点数最多的两个轮廓,判断剩余轮廓个数κζ的大小,分为如下情况:若κζ==0,则执行步骤13;否则这κζ个轮廓即为第ζ个平行四边形编码标志内编码标志圆的轮廓,分别记为编码标志圆编码标志圆…、编码标志圆步骤12.2、给整数变量i并赋予初值i=1;
步骤12.3、在第ζ组第2个备份二值化图像P′ζ,2上,计算编码标志圆轮廓的质心将i+1重新赋值给i后继续执行此步骤,直到i>κζ结束;由此可以得到对应编码标志圆编码标志圆…、编码标志圆的质心像素坐标
步骤13.1、在第ζ组第2个备份二值化图像P′ζ,2上,记在第ζ个平行四边形编码标志中位于第1编码区域、第3编码区域、第4编码区域和第6编码区域的特征角点的亚像素坐标分别为Cζ,1(xζ,1,yζ,1)、Cζ,2(xζ,2,yζ,2)、Cζ,3(xζ,3,yζ,3)、Cζ,4(xζ,4,yζ,4);并在第ζ组第2个备份二值化图像P′ζ,2上,将像素坐标值为的(x″d,ζ,y″d,ζ)的像素点记为第ζ组第2个备份二值化图像P′ζ,2上第ζ个平行四边形编码标志的定向圆环质心o′d,ζ(x′d,ζ,y′d,ζ),并将亚像素坐标值为(x′ζ,1,y′ζ,1)、(x′ζ,2,y′ζ,2)、(x′ζ,3,y′ζ,3)、(x′ζ,4,y′ζ,4)的4个像素点分别记为C″′ζ,1(x″′ζ,1,y″′ζ,1)、C″′ζ,2(x″′ζ,2,y″′ζ,2)、C″′ζ,3(x″′ζ,3,y″′ζ,3)、C″′ζ,4(x″′ζ,4,y″′ζ,4);
步骤14、在第ζ组第2个备份二值化图像P′ζ,2上,找出4个像素点C″′ζ,1(x″′ζ,1,y″′ζ,1)、C″′ζ,2(x″′ζ,2,y″′ζ,2)、C″′ζ,3(x″′ζ,3,y″′ζ,3)、C″′ζ,4(x″′ζ,4,y″′ζ,4)中距离定位圆质心的像素坐标o′l,ζ(x′l,ζ,y′l,ζ)最近的2个像素点,分别记为Cζ,1min(xζ,1min,yζ,1min)和Cζ,2min(xζ,2min,yζ,2min);
步骤15、在第ζ组第2个备份二值化图像P′ζ,2上,将余下的2个像素点的亚像素坐标值分别赋值给第ζ个平行四边形编码标志的临时坐标值1,记为C′ζ,5(x′ζ,5,y′ζ,5),以及临时坐标值2,记为C′ζ,6(x′ζ,6,y′ζ,6);
根据sinωζ和sinξζ的取值情况,将C′ζ,5(x′ζ,5,y′ζ,5)或C′ζ,6(x′ζ,6,y′ζ,6)的亚像素坐标值赋值给Cζ,2(xζ,2,yζ,2)或Cζ,3(xζ,3,yζ,3);
在第ζ组第2个备份二值化图像P′ζ,2上,以第ζ个平行四边形编码标志中的定位圆质心o′l,ξ(x′l,ζ,y′l,ζ)为起点做与第5判断向量平行且同向的单位向量,记为以第ζ个平行四边形编码标志中定向圆环质心o′d,ζ(x′d,ζ,y′d,ζ)为起点做与第5判断向量平行且同向的单位向量,记为
步骤18、定义6个浮点型二维数组 用于存放第ζ个平行四边形编码标志中分别位于第1编码区域、第2编码区域、第3编码区域、第4编码区域、第5编码区域和第6编码区域的编码标志圆轮廓质心在第ζ组第2个备份二值化图像P′ζ,2上的像素坐标,初始化这6个二维数组中的所有元素,且赋值为-1;另外取6个整数变量并将其初始化,取整数变量i并赋值i=1;
步骤19、在第ζ组第2个备份二值化图像P′ζ,2上,根据式(12)、式(13)计算出第ζ个平行四边形编码标志中编码标志圆轮廓的质心像素坐标分别与定位圆质心o′l,ζ和定向环质心o′d,ζ所形成的第i组第1个象限向量和第i组第2个象限向量
Wζ=Vζ T·U (18)
至此获得了第ζ个平行四边形编码标志中4个特征角点各自的编码号与其相应的亚像素坐标值:
步骤24、判断ζ是否小于N,若满足,则将ζ+1重新赋值给ζ,并返回步骤6.1顺序执行;若不满足,则已经找到每个平行四边形编码标志上4个特征角点的亚像素坐标值和对应的编码号。
2.根据权利要求1所述的解码方法,其特征在于:步骤6.1中,以获得的第ζ组第1个备份二值化图像P′ζ,1以及圆环黑色连通域质心集合A中第ζ个圆环黑色连通域质心像素坐标值(x″d,ζ,y″d,ζ)作为输入条件,利用去除编码标志复杂背景算法去除第ζ组第1个备份二值化图像P′ζ,1的复杂背景,并最终获得第ζ个无复杂背景平行四边形二值化图像P″ζ的步骤如下:
步骤6.1.1、取整数变量i并赋值i=1;
步骤6.1.2、在第ζ组第1个备份二值化图像P′ζ,1上寻找距离像素坐标值(x″d,ζ,y″d,ζ)最近的白色连通域并记为第ζ组第1个备份二值化图像P′ζ,1上的第i个圆环白色连通域Μ′i,并将第i个圆环白色连通域Μ′i的灰度值均赋值为0;
步骤6.1.3、判断i是否小于N,若i<N,则将i+1赋值给i,返回步骤1.2顺序执行;否则,已将第ζ组第1个备份二值化图像P′ζ,1中所有平行四边形编码标志中的定向圆环置为黑色,并将第ζ组第1个备份二值化图像P′ζ,1经过本步骤处理后的图像记为第ζ个黑色背景提取图像P′ζ,b,并执行步骤6.1.4;
步骤6.1.4、在第ζ个黑色背景提取图像P′ζ,b上的所有黑色连通域中,寻找距离像素坐标值(x″d,ζ,y″d,ζ)最近的黑色连通域记为第ζ个黑色背景提取图像P′ζ,b上的第ζ个平行四边形背景黑色连通域Ωζ;
步骤6.1.5、在第ζ个黑色背景提取图像P′ζ,b上提取第ζ个平行四边形背景黑色连通域Ωζ的所有内、外轮廓,获得第ζ个平行四边形背景黑色连通域Ωζ的内、外轮廓集合Vζ;
步骤6.1.6、比较轮廓集合Vζ中所有轮廓的周长,将周长最长的轮廓记为第ζ个平行四边形背景黑色连通域Ωζ的外轮廓Εζ;
步骤6.1.7、在第ζ个黑色背景提取图像P′ζ,b上,将外轮廓Εζ以内的所有像素点的灰度值赋值为255,外轮廓Εζ以外的所有像素点的灰度值赋值为0,从而获得第ζ个无复杂背景平行四边形二值化图像P″ζ;使得在第ζ个无复杂背景平行四边形二值化图像P″ζ中,包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志内部的所有像素点灰度值被置为255,包含质心像素坐标值(x″d,ζ,y″d,ζ)的平行四边形编码标志外部的所有像素点灰度值被置为0。
3.根据权利要求1所述的解码方法,其特征在于:步骤14中,分以下两种情况执行:
(1)若sinαζ<0,sinβζ>0,则Cζ,1min(xζ,1min,yζ,1min)为第ζ个平行四边形编码标志中第1编码区域的特征角点,把Cζ,1min(xζ,1min,yζ,1min)的亚像素坐标值赋值给Cζ,1(xζ,1,yζ,1);Cζ,2min(xζ,2min,yζ,2min)为第ζ个平行四边形编码标志中第6编码区域的特征角点,把Cζ,2min(xζ,2min,yζ,2min)的亚像素坐标值赋值给Cζ,4(xζ,4,yζ,4);
(2)若sinαζ>0,sinβζ<0,则Cζ,2min(xζ,2min,yζ,2min)为第ζ个平行四边形编码标志中第1编码区域的特征角点,把Cζ,2min(xζ,2min,yζ,2min)的像素坐标值赋值给Cζ,1(xζ,1,yζ,1);Cζ,1min(xζ,1min,yζ,1min)为第ζ个平行四边形编码标志中第6编码区域的特征角点,把Cζ,1min(xζ,1min,yζ,1min)的像素坐标值赋值给Cζ,4(xζ,4,yζ,4)。
4.根据权利要求1所述的解码方法,其特征在于:步骤16中,分以下两种情况执行:
(1)若sinωζ==0,sinξζ≠0,则第1临时坐标C′ζ,5(x′ζ,5,y′ζ,5)即为第ζ个平行四边形编码标志中第3编码区域的特征角点的亚像素坐标,把C′ζ,5(x′ζ,5,y′ζ,5)的亚像素坐标值赋值给Cζ,2(xζ,2,yζ,2);第2临时坐标C′ζ,6(x′ζ,6,y′ζ,6)即为第ζ个平行四边形编码标志中第4编码区域的特征角点的亚像素坐标,把C′ζ,6(x′ζ,6,y′ζ,6)的亚像素坐标值赋值给Cζ,3(xζ,3,yζ,3);
(2)若sinωζ≠0,sinξζ==0,则第2临时坐标C′ζ,6(x′ζ,6,y′ζ,6)即为第ζ个平行四边形编码标志中第3编码区域的特征角点的亚像素坐标,把C′ζ,6(x′ζ,6,y′ζ,6)的亚像素坐标值赋值给Cζ,2(xζ,2,yζ,2);第1临时坐标C′ζ,5(x′ζ,5,y′ζ,5)即为第ζ个平行四边形编码标志中第4编码区域的特征角点的亚像素坐标,把C′ζ,5(x′ζ,5,y′ζ,5)的亚像素坐标值赋值给Cζ,3(xζ,3,yζ,3)。
5.根据权利要求1所述的解码方法,其特征在于:步骤19中,分为以下情况执行:
步骤21.1:
步骤21.2、本步骤分为以下情况执行:
步骤21.3、本步骤分为以下情况执行:
7.一种计算机可读存储介质,包括与具有图像处理功能的电子设备结合使用的计算机程序,所述计算机程序可被处理器执行以如权利要求1所述的解码方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011546292 | 2020-12-23 | ||
CN2020115462920 | 2020-12-23 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113129396A true CN113129396A (zh) | 2021-07-16 |
CN113129396B CN113129396B (zh) | 2022-10-14 |
Family
ID=76775872
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110392326.3A Active CN113129396B (zh) | 2020-12-23 | 2021-04-13 | 一种基于区域分割的平行四边形编码标志的解码方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113129396B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116652440A (zh) * | 2023-07-31 | 2023-08-29 | 中建安装集团西安建设投资有限公司 | 钢结构数字化焊接方法及控制*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050286741A1 (en) * | 2004-06-29 | 2005-12-29 | Sanyo Electric Co., Ltd. | Method and apparatus for coding images with different image qualities for each region thereof, and method and apparatus capable of decoding the images by adjusting the image quality |
CN101398895A (zh) * | 2007-09-26 | 2009-04-01 | 杨高波 | 一种数据矩阵(dm)二维条形码识别的图像预处理方法 |
CN108537217A (zh) * | 2018-04-04 | 2018-09-14 | 湖南科技大学 | 基于字符编码标志的识别及定位方法 |
-
2021
- 2021-04-13 CN CN202110392326.3A patent/CN113129396B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050286741A1 (en) * | 2004-06-29 | 2005-12-29 | Sanyo Electric Co., Ltd. | Method and apparatus for coding images with different image qualities for each region thereof, and method and apparatus capable of decoding the images by adjusting the image quality |
CN101398895A (zh) * | 2007-09-26 | 2009-04-01 | 杨高波 | 一种数据矩阵(dm)二维条形码识别的图像预处理方法 |
CN108537217A (zh) * | 2018-04-04 | 2018-09-14 | 湖南科技大学 | 基于字符编码标志的识别及定位方法 |
Non-Patent Citations (2)
Title |
---|
D. OLIVIER等: "Supervised segmentation at low bit rates for region representation and color image compression", 《IEEE XPLORE》 * |
金滔等: "一种环状编码标记点的设计及解码算法", 《计算机应用研究》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116652440A (zh) * | 2023-07-31 | 2023-08-29 | 中建安装集团西安建设投资有限公司 | 钢结构数字化焊接方法及控制*** |
CN116652440B (zh) * | 2023-07-31 | 2023-10-20 | 中建安装集团西安建设投资有限公司 | 钢结构数字化焊接方法及控制*** |
Also Published As
Publication number | Publication date |
---|---|
CN113129396B (zh) | 2022-10-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111243032B (zh) | 一种棋盘格角点全自动检测方法 | |
CN109215016B (zh) | 一种编码标志的识别定位方法 | |
Ha et al. | Deltille grids for geometric camera calibration | |
CN110246168A (zh) | 一种移动巡检机器人双目图像拼接的特征匹配方法 | |
CN108573511B (zh) | 点状分布合作编码标志及其识别定位方法 | |
CN110738273B (zh) | 图像特征点的匹配方法、装置、设备及存储介质 | |
CN113129397B (zh) | 一种基于图形几何关系的平行四边形编码标志的解码方法 | |
CN109559346B (zh) | 一种3d点云测量中的待测部位的定位及分割方法、扫描仪 | |
CN113129385B (zh) | 基于空间中多编码平面靶标的双目摄像机内外参数标定方法 | |
CN111444948B (zh) | 一种图像特征提取与匹配方法 | |
CN113096191B (zh) | 基于编码平面靶标的单目摄像机智能标定方法 | |
CN113129396B (zh) | 一种基于区域分割的平行四边形编码标志的解码方法 | |
CN113160329B (zh) | 用于摄像机标定的编码平面靶标及其解码方法 | |
CN112257721A (zh) | 一种基于Fast ICP的图像目标区域匹配方法 | |
CN113129386B (zh) | 基于编码平面靶标的双目摄像机内外参数的智能标定方法 | |
CN115222828A (zh) | 用于视觉定位的正六边形编码标志解码方法 | |
CN112784761A (zh) | 一种特殊纹理背景遥感图像匹配方法 | |
CN115272459A (zh) | 一种四角二维码定向平面靶标装置及其标定方法 | |
CN113112550B (zh) | 一种用于摄像机内外参数标定的编码平面靶标及其编码方法 | |
CN113112549B (zh) | 一种基于编码立体靶标的单目摄像机快速标定方法 | |
CN110197509B (zh) | 一种基于彩色人工标识的相机位姿求解法 | |
Guo et al. | Automatic shape-based target extraction for close-range photogrammetry | |
Wang et al. | A novel cross-circular coded target for photogrammetry | |
CN113192143B (zh) | 一种用于摄像机快速标定的编码立体靶标及其解码方法 | |
CN116485882A (zh) | 一种编码容量可变的编码标志的解码方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |