CN109784121B - Dot-peen DPM码的识别方法及装置 - Google Patents

Dot-peen DPM码的识别方法及装置 Download PDF

Info

Publication number
CN109784121B
CN109784121B CN201910031689.7A CN201910031689A CN109784121B CN 109784121 B CN109784121 B CN 109784121B CN 201910031689 A CN201910031689 A CN 201910031689A CN 109784121 B CN109784121 B CN 109784121B
Authority
CN
China
Prior art keywords
point
points
line
coordinates
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.)
Active
Application number
CN201910031689.7A
Other languages
English (en)
Other versions
CN109784121A (zh
Inventor
何学智
黄自力
郑艺祥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Newland Digital Technology Co ltd
Original Assignee
Newland Digital Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Newland Digital Technology Co ltd filed Critical Newland Digital Technology Co ltd
Priority to CN201910031689.7A priority Critical patent/CN109784121B/zh
Publication of CN109784121A publication Critical patent/CN109784121A/zh
Application granted granted Critical
Publication of CN109784121B publication Critical patent/CN109784121B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明公开了Dot‑peen DPM码的识别方法及装置,首先通过极值检测获取图像中的点;其次通过检测出的点的位置进行主方向确认,主方向为二维码边缘L线的两个方向,检测定位主方向上的L线;而后通过迭代透视变换生成采样网络;最后对透视变换后的所有点坐标进行采样及二值化处理。采用本方案,可以有效应对DPM码中的Dot‑peen类型的点,并对金属反光、光照不均匀、噪点复杂等情况的DPM码有很好的解码效果;此外,本方案对Dot‑peen DPM码中存在部分点污损的情况有很好的容忍度,可以精确识别带有污渍的Dot‑peen DPM码。

Description

Dot-peen DPM码的识别方法及装置
技术领域
本发明涉及图像识别技术领域,特别涉及Dot-peen DPM码的识别方法及装置。
背景技术
DPM码(direct part marking,直接零件标识)技术是一种直接刻印在产品表面的永久条码标识,而Dot-peenDPM码,是通过针打(Dot-peen)将二维码永久标在产品表面的方法。不同于纸质条码,DPM条码由于其附着的金属面粗糙、存在刮痕、高反光、背景噪点多等原因,使得拍摄所得的条码图像一般都属于低品质条码。对于条码的读取,在目前的定位和解码算法中,普遍采用边缘检测方法,对二维码的特征部分进行搜索、定位、识别。而对于Dot-peen DPM,由于其模块之间不连续,边缘检测无法使用,导致目前的解码算法无法对其解码,因此需要对该情况的条码研发其特有的解码算法。
发明内容
本发明要解决的技术问题是如何提供一种针对Dot-peenDPM码边缘不连续、图像质量低的特点的Dot-peen DPM码的识别方法及装置。
为了解决上述技术问题,本发明的技术方案为:
Dot-peen DPM码的识别方法,包括步骤:
通过极值检测获取图像中的点;
通过检测出的点的位置进行主方向确认,所述主方向为二维码边缘L线的两个方向;
检测定位主方向上的L线;
通过迭代透视变换生成采样网络;
对透视变换后的所有点坐标进行采样及二值化处理。
优选地,步骤通过检测出的点的位置进行主方向确认,包括:
对于每个点,找到与其曼哈顿距离最小的两个相邻点;
将所有相邻点连接起来,形成连通分量;
统计所有相邻点连线的方向,将数量最多的两个方向作为条码的两个主方向θ1,θ2
优选地,步骤检测定位主方向上的L线,包括:
剔除与主方向之间的夹角大于预设值的连线,获得θ1方向的连线集合L′1、θ2方向的连线集合L′2
将L′1的平均长度作为θ1方向的模块宽度s1,将L′2的平均长度作为θ2方向的模块宽度s2
采用Hough变换得到主方向的L线。
优选地,步骤通过迭代透视变换生成采样网络,包括:
由L线交点向两个主方向延伸,获取L线上所有点的模板坐标;
建立L线上点的实际坐标及模板坐标之间的透视变换公式;
将L线上的点的实际坐标及模板坐标代入透视变换公式,求得透视变换系数;
从L线向条码区域内部延伸,获取最靠近L的线段Lin-1所有模板坐标,通过透视变换计算可以得到Lin-1的所有模板坐标对应的伪实际坐标;
设置搜索半径,若在伪实际坐标的搜索半径内存在实际的图像检测点,则建立Lin-1上的实际坐标与模板坐标的映射对;
将L、Lin-1上所有的实际坐标、模板坐标映射对代入透视变换公式,更新透视变换系数;
从Lin-1向条码区域内部延伸,重复直至遍历DPM上所有检测点,生成采样网络。
优选地,通过极值检测获取图像中的点具体为:
设I、J是定义于同一离散域D上的两张灰度图像,I为掩模图像,J为标识图像,且
Figure BDA0001944455160000021
其中p表示图像上的任意像素点,J(p)、I(p)是p点在图像J、I中的灰度值,则用J对I做形态学重建表示为:
Figure BDA0001944455160000022
其中,
Figure BDA0001944455160000023
Figure BDA0001944455160000024
表示形态学膨胀操作,B为尺度为1的结构元素,∧表示逐点比较取最小值;
采用形态学重建,标识图像J满足:
Figure BDA0001944455160000025
形态学重建表示为ρI(J)=ρI(I-h),其中h为预设常数;
用于点检测的公式为:Sh(I)=I-ρI(I-h),S表示得到的标注区域极值点的图像。
本发明还提出一种Dot-peen DPM码的识别装置,包括:
点检测模块:通过极值检测获取图像中的点;
主方向模块:通过检测出的点的位置进行主方向确认,所述主方向为二维码边缘L线的两个方向;
L线模块:检测定位主方向上的L线;
迭代模块:通过迭代透视变换生成采样网络;
二值化模块:对透视变换后的所有点坐标进行采样及二值化处理。
优选地,主方向模块包括:
距离单元:对于每个点,找到与其曼哈顿距离最小的两个相邻点;
连通单元:将所有相邻点连接起来,形成连通分量;
统计单元:统计所有相邻点连线的方向,将数量最多的两个方向作为条码的两个主方向θ12
优选地,L线模块包括:
剔除单元:剔除与主方向之间的夹角大于预设值的连线,获得θ1方向的连线集合L′1、θ2方向的连线集合L′2
宽度单元:将L′1的平均长度作为θ1方向的模块宽度s1,将L′2的平均长度作为θ2方向的模块宽度s2
Hough变换单元:采用Hough变换得到主方向的L线。
优选地,迭代模块包括:
延伸单元:由L线交点向两个主方向延伸,获取L线上所有点的模板坐标;
公式单元:建立L线上点的实际坐标及模板坐标之间的透视变换公式;
系数单元:将L线上的点的实际坐标及模板坐标代入透视变换公式,求得透视变换系数;
伪实际坐标单元:从L线向条码区域内部延伸,获取最靠近L的线段Lin-1所有模板坐标,通过透视变换计算可以得到Lin-1的所有模板坐标对应的伪实际坐标;
映射对单元:设置搜索半径,若在伪实际坐标的搜索半径内存在实际的图像检测点,则建立Lin-1上的实际坐标与模板坐标的映射对;
系数更新单元:将L、Lin-1上所有的实际坐标、模板坐标映射对代入透视变换公式,更新透视变换系数;
遍历单元:从Lin-1向条码区域内部延伸,重复直至遍历DPM上所有检测点,生成采样网络。
优选地,点检测模块中,通过极值检测获取图像中的点具体为:
设I、J是定义于同一离散域D上的两张灰度图像,I为掩模图像,J为标识图像,且
Figure BDA0001944455160000031
其中p表示图像上的任意像素点,J(p)、I(p)是p点在图像J、I中的灰度值,则用J对I做形态学重建表示为:
Figure BDA0001944455160000032
其中,
Figure BDA0001944455160000033
Figure BDA0001944455160000034
表示形态学膨胀操作,B为尺度为1的结构元素,∧表示逐点比较取最小值;
采用形态学重建,标识图像J满足:
Figure BDA0001944455160000041
形态学重建表示为ρI(J)=ρI(I-h),其中h为预设常数;
用于点检测的公式为:Sh(I)=I-ρI(I-h),S表示得到的标注区域极值点的图像。
采用上述技术方案,通过对针打点进行检测,对点进行连线进而找出各个点形成的采样网络,生成二维码图像。通过由点组成线,再由线构成面的形式,对Dot-peenDPM进行图像识别。采用本方案,可以有效应对DPM码中的Dot-peen类型,并对金属反光、光照不均匀、噪点复杂等情况有很好的解码效果;此外,本方案对Dot-peenDPM码中存在部分点污损的情况有很好的容忍度,可以精确识别带有污渍的Dot-peen DPM码。
附图说明
图1为本发明Dot-peen DPM码的识别方法一实施例的流程图;
图2为图1中步骤S20的流程图;
图3为图1中步骤S30的流程图;
图4为图1中步骤S40的流程图;
图5为本发明Dot-peen DPM码的识别装置一实施例的原理模块图;
图6为本发明Dot-peen DPM码的识别方法未处理之前的图像;
图7为图6中极值点的图像;
图8为图6中标注区域极值点后的图像;
图9为步骤S23中相邻点连线的方向的统计直方图;
图10为找到L线后的图像;
图11为步骤S40中的第一张迭代点的示意图;
图12为步骤S40中的第二张迭代点的示意图。
图13为步骤S40中的第三张迭代点的示意图。
图14为步骤S40中的第四张迭代点的示意图。
图15为步骤S40中的第五张迭代点的示意图。
图16为步骤S40中的第六张迭代点的示意图。
图17为步骤S40中的第七张迭代点的示意图。
图18为图6通过Dot-peen DPM码的识别方法处理后提取的二维码图像。
图中,10-点检测模块,20-L线模块,30-主方向模块,40-迭代模块,50-二值化模块。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明。在此需要说明的是,对于这些实施方式的说明用于帮助理解本发明,但并不构成对本发明的限定。此外,下面所描述的本发明各个实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互组合。
参照图1,本发明的技术方案为Dot-peen DPM码的识别方法,包括步骤:
S10:通过极值检测获取图像中的点;
通过极值检测获取图像中的点具体为:
设I、J是定义于同一离散域D上的两张灰度图像,I为掩模图像,J为标识图像,且
Figure BDA0001944455160000051
其中p表示图像上的任意像素点,J(p)、I(p)是p点在图像J、I中的灰度值,则用J对I做形态学重建表示为:
Figure BDA0001944455160000052
其中,
Figure BDA0001944455160000053
Figure BDA0001944455160000054
表示形态学膨胀操作,B为尺度为1的结构元素,∧表示逐点比较取最小值;
采用形态学重建,标识图像J满足:
Figure BDA0001944455160000055
形态学重建表示为ρI(J)=ρI(I-h),其中h为预设常数;
用于点检测的公式为:Sh(I)=I-ρI(I-h),S表示得到的标注区域极值点的图像。
通过极值检测获取图像中的点具体为:
设I、J是定义于同一离散域D上的两张灰度图像,I为掩模图像,J为标识图像,且
Figure BDA0001944455160000056
其中p表示图像上的任意像素点,J(p)、I(p)是p点在图像J、I中的灰度值,则用J对I做形态学重建表示为:
Figure BDA0001944455160000057
其中,
Figure BDA0001944455160000058
Figure BDA0001944455160000059
表示形态学膨胀操作,B为尺度为1的结构元素,∧表示逐点比较取最小值;
采用形态学重建,标识图像J满足:
Figure BDA00019444551600000510
形态学重建表示为ρI(J)=ρI(I-h),其中h为预设常数;
用于点检测的公式为:Sh(I)=I-ρI(I-h),S表示得到的标注区域极值点的图像。
S20:通过检测出的点的位置进行主方向确认,主方向为二维码边缘L线的两个方向;
参照图2,步骤S20具体的包括:
S21:对于每个点,找到与其曼哈顿距离最小的两个相邻点;
S22:将所有相邻点连接起来,形成连通分量;
S23:统计所有相邻点连线的方向,将数量最多的两个方向作为条码的两个主方向θ12
参照图3,S30:检测定位主方向上的L线;
步骤S30具体的包括
S31:剔除与主方向之间的夹角大于预设值的连线,获得θ1方向的连线集合L′1、θ2方向的连线集合L′2
S32:将L′1的平均长度作为θ1方向的模块宽度s1,将L′2的平均长度作为θ2方向的模块宽度s2
S33:采用Hough变换得到主方向的L线。本实施例中,Hough变换的具体方法为:
采用Hough变换得到主方向的L线。Hough变换的原理是将点变换到一组参数空间(λ,θ)上,根据参数空间点的累计结果找到一个极大值对应的解,那么这个解就对应着要寻找的直线参数。由于这里只需要寻找两个主方向θ12,所以采用Hough变换得到主方向的L线的步骤如下:
先计算L线的θ1方向的线段:
1、建立参数空间(λ,θ),这里设θ=θ1
2、顺序搜索U′={Pi,i=1,2...,n′}中所有点,根据λ=xcosθ1+ysinθ1计算λ={λi,i=1,2...,n′};
3、对数组λ进行直方图统计,以2为步进尺寸,求得
Figure BDA0001944455160000061
其中λj为直方图上落入最大统计量区间的所有M个λ值。
4、求得L线的θ1方向的线段L1方程为:λ1=xcosθ1+ysinθ1
同理可计算得L线的θ2方向的线段L2方程为:λ2=xcosθ2+ysinθ2
参照图4,S40:通过迭代透视变换生成采样网络;
步骤S40具体的包括:
S41:由L线交点向两个主方向延伸,获取L线上所有点的模板坐标;
S42:建立L线上点的实际坐标及模板坐标之间的透视变换公式;
S43:将L线上的点的实际坐标及模板坐标代入透视变换公式,求得透视变换系数;
S44:从L线向条码区域内部延伸,先将最靠近L的线段Lin-1,可获取Lin-1的所有模板坐标,通过透视变换计算可以得到Lin-1的所有模板坐标对应的伪实际坐标;
S45:设置搜索半径,若在搜索半径内存在实际的图像检测点,则建立Lin-1上的实际坐标与模板坐标的映射对;
S46:将L、Lin-1、上所有的实际坐标、模板坐标映射对代入透视变换公式,更新透视变换系数;
从Lin-1向条码区域内部延伸,重复直至遍历DPM上所有检测点,生成采样网络。
S50:对透视变换后的所有点坐标进行采样及二值化处理。
在本发明的另一实施例中,以图6的图像为例,该方法通过以下步骤实现:
1、点检测;
这里我们采用基于H-dome形态学重建的区域极值检测方法。由于针打的点,相对于其领域是极大值或极小值,因此我们可以采用极值检测的方法来得到各个针打的点分布。形态学重建是基于测地膨胀的形态学运算,主要用来提取区域极值。
对于灰度形态学重建的定义如下:设I、J是定义于同一离散域D上的两张灰度图像,I为掩模图像,J为标识图像,且
Figure BDA0001944455160000071
其中p表示图像上的任意像素点,J(p),I(p)是p点在图像J、I中的灰度值。用J对I做形态学重建可表示为:
Figure BDA0001944455160000072
其中,
Figure BDA0001944455160000073
Figure BDA0001944455160000074
表示形态学膨胀操作,B为尺度为1的结构元素,∧表示逐点比较取最小值。由于I的限制,n次迭代最终会趋于收敛。
这里,我们采用H-dome形态学重建,标识图像J满足:
Figure BDA0001944455160000075
其中h为预设常数。即将图像I的每个像素的灰度值降低h,用J对I进行重建,得到的图像可以保持图像I的形态,消除幅度小于h的微小差异,有效的从灰度图中提取出各极大区域。对于Dot-peen DPM的点,在图像上会呈现其亮度较高(或较暗),而目标四周的灰度较低(或较高),考虑到Dot-peen DPM码存在反光、不均匀光照、背景噪点多的问题,这里我们采用自适应计算h的方法,即对不同图像区域,计算不同的h,选取h小于Dot与背景的灰度均值的差,而大于一般的局部极大值点的高度,就能很好的保留Dot点而消除周围的复杂背景。
最终用于点检测的公式如下:Sh(I)=I-ρI(I-h),S表示得到的标注区域极值点的图像,如图7及8所示。
2、条码主方向估计;
标准的二维DPM码两个主方向是互为垂直的,以DM为例,两个主方向为边缘L线的两个方向。步骤1已经获得了点的空间分布,由于Dot-peenDPM码存在背景噪点、对比度不均匀等问题,因此检测出来的点不一定就是DPM码的信息点模块,而DPM码的信息点模块也不一定会完全被检测出来。因此通过条码主方向估计这个步骤对已检测的点空间分布进行DPM码主方向估计,为后续的条码信息提取做准备,也可以排除掉其他噪点的干扰。
本步骤的主要过程如下:
2.1对于检测出来的每个点,找到其最近相邻距离的两个点,定义为相邻点;
这里距离的定义如下:对于空间中任意两个点,其坐标分别为(xi,yi),(xj,yj),采用曼哈顿距离(ManhattanDistance),Disij=|xi-xj|+yi-yj|;
采用K-近邻(KNN)算法,这里k设为2,可以获得每个点的最近两个距离的点。K-近邻算法搜索的复杂度为O(N2),这里可以引入KD-Tree的优化,将复杂度下降为O(NlogN),加快速度。
2.2将所有相邻点连接起来,形成连通分量。
将空间中所有的相邻点连接起来,会形成很多个连通分量。由于DPM码相对于背景比较明显,因此虽存在一定的噪点,但是连通分量中包含点最多的,就是DPM区域。
2.3统计DPM区域所有相邻点连线的方向,得到如图9所示的方向直方图;
所有相邻点连线的集合记为L={Lij},所有相邻点连线的方向的集合记为θ={θij},下标ij表示相邻点(xi,yi),(xj,yj)。θij方向定义为:
Figure BDA0001944455160000081
设定eps=0.1,
Figure BDA0001944455160000082
表示向下取整。对θ={θij},可以绘制其直方图,得到方向直方图。
2.4方向直方图最大的两个方向即为条码的两个主方向,两个主方向表示为θ12。主方向确认后的图像如图10所示。
3、L线搜索;
采用Hough变换得到主方向的L线。Hough变换的原理是将点变换到一组参数空间(λ,θ)上,根据参数空间点的累计结果找到一个极大值对应的解,那么这个解就对应着要寻找的直线参数。由于这里只需要寻找两个主方向θ12,所以采用Hough变换得到主方向的L线的步骤如下:
先计算L线的θ1方向的线段:
(1)、建立参数空间(λ,θ),这里设θ=θ1
(2)、顺序搜索U′={Pi,i=1,2...,n′}中所有点,根据λ=xcosθ1+ysinθ1计算λ={λi,i=1,2...,n′};
(3)、对数组λ进行直方图统计,以2为步进尺寸,求得
Figure BDA0001944455160000091
其中λj为直方图上落入最大统计量区间的所有M个λ值。
(4)、求得L线的θ1方向的线段L1方程为:λ1=xcosθ1+ysinθ1
同理可计算得L线的θ2方向的线段L2方程为:λ2=xcosθ2+ysinθ2
4、迭代透视变换生成采样网络;
(1)、求得L1、L2两直线的交点;
(2)、通过求得的两个方向模块宽度s1,s2,从交点分别向两个方向延伸,可获得DPML线上所有检测点的模板坐标。
(3)、根据L线上所有点的模板坐标和实际坐标映射对做迭代透视变换,生成采样网络。迭代的具体步骤如下:
建立透视变换公式,U=(ax+by+c)/(gx+hy+1),V=(dx+ey+f)/(gx+hy+1),其中,(U,V)为点的实际坐标,(x,y)为点的模板坐标。
将L线上的点的实际坐标(UL,VL),模板坐标(xL,yL)代入透视变换公式,求得透视变换系数a、b、c、d、e、f、g、h;
如图11至图17所示,本实施例每次迭代递进两行,向二维码内部位置迭代,图案中方形点为处理后的点,圆形点为待处理的点。从L线往条码区域内部延伸,考虑精度和速度的折中,先考虑最靠近L的两条线段Lin-1,Lin-2,可获取Lin-1,Lin-2的所有模块中心的模板坐标,通过透视变换计算可以得到Lin-1,Lin-2的所有模块中心模板坐标对应的伪实际坐标,对这些伪实际坐标搜索其距离最近的实际点坐标。搜索方式为:设置搜索半径γ,若在搜索半径γ内存在实际的图像检测点,则建立Lin-1线上的实际坐标(Uin-1,Vin-1)与模板坐标(xin-1,yin-1)的映射对以及Lin-2线上的实际坐标(Uin-2,Vin-2)与模板坐标(xin-2,yin-2)的映射对。
将L、Lin-1、Lin-2上所有的实际坐标、模板坐标映射对代入透视变换公式,更新透视变换系数a、b、c、d、e、f、g、h;
从Lin-2再往内延伸,重复进行步骤iii、iv,直到遍历DPM上所有检测点,即可得到最终的透视变换系数a、b、c、d、e、f、g、h;
根据透视变换,生成采样网络。
需要说明的是,本实施例中,一次迭代为递进两层线Lin-1、Lin-2。在实际实施时,一次迭代的递进线数可以为一层或者两层等情况,不仅仅限制为递进一层。
通过上述迭代透视变换求解,可以有效地减少远离L线的检测点的透视变换的误差,使得采样网络更准确。
5、二值化;
对透视变换完的所有点坐标进行采样,进行二值化,判断其是0模块还是1模块。这里可以采用Niblack二值化方法。得到如图18的二值化后的编制二维码图像结果。
6、纠错、译码。
对上述二值化的结果进行纠错、译码等常规条码识别操作,即可获取该Dot-peenDPM码的信息。
通过以上步骤,将边缘不连续的Dot-peenDPM码进行点检测、主方向确认、边缘L线确认、迭代获取采样网络、二值化、纠错、译码等步骤,实现了对Dot-peenDPM码的识别,抗噪能力强。
如图5,本发明还提出了Dot-peen DPM码的识别装置,包括:
点检测模块10:通过极值检测获取图像中的点;
主方向模块20:通过检测出的点的位置进行主方向确认,主方向为二维码边缘L线的两个方向;
L线模块30:检测定位主方向上的L线;
迭代模块40:通过迭代透视变换生成采样网络;
二值化模块50:对透视变换后的所有点坐标进行采样及二值化处理。
具体地,主方向模块20包括:
距离单元:对于每个点,找到与其曼哈顿距离最小的两个相邻点;
连通单元:将所有相邻点连接起来,形成连通分量;
统计单元:统计所有相邻点连线的方向,将数量最多的两个方向作为条码的两个主方向θ12
具体地,L线模块30包括:
剔除单元:剔除与主方向之间的夹角大于预设值的连线,获得θ1方向的连线集合L′1、θ2方向的连线集合L′2
宽度单元:将L′1的平均长度作为θ1方向的模块宽度s1,将L′2的平均长度作为θ2方向的模块宽度s2
Hough变换单元:采用Hough变换得到主方向的L线。
具体地,迭代模块40包括:
延伸单元:由L线交点向两个主方向延伸,获取L线上所有点的模板坐标;
公式单元:建立L线上点的实际坐标及模板坐标之间的透视变换公式;
系数单元:将L线上的点的实际坐标及模板坐标代入透视变换公式,求得透视变换系数;
伪实际坐标单元:从L线向条码区域内部延伸,获取最靠近L的线段Lin-1所有模板坐标,通过透视变换计算可以得到Lin-1的所有模板坐标对应的伪实际坐标;
映射对单元:设置搜索半径,若在伪实际坐标的搜索半径内存在实际的图像检测点,则建立Lin-1上的实际坐标与模板坐标的映射对;
系数更新单元:将L、Lin-1上所有的实际坐标、模板坐标映射对代入透视变换公式,更新透视变换系数;
遍历单元:从Lin-1向条码区域内部延伸,重复直至遍历DPM上所有检测点,生成采样网络。
具体地,点检测模块中10,通过极值检测获取图像中的点具体为:
设I、J是定义于同一离散域D上的两张灰度图像,I为掩模图像,J为标识图像,且
Figure BDA0001944455160000111
其中p表示图像上的任意像素点,J(p)、I(p)是p点在图像J、I中的灰度值,则用J对I做形态学重建表示为:
Figure BDA0001944455160000112
其中,
Figure BDA0001944455160000113
Figure BDA0001944455160000114
表示形态学膨胀操作,B为尺度为1的结构元素,∧表示逐点比较取最小值;
采用形态学重建,标识图像J满足:
Figure BDA0001944455160000115
形态学重建表示为ρI(J)=ρI(I-h),其中h为预设常数;
用于点检测的公式为:Sh(I)=I-ρI(I-h),S表示得到的标注区域极值点的图像。
采用上述技术方案,通过对针打点进行检测,对点进行连线进而找出各个点形成的采样网络,生成二维码图像。通过由点组成线,再由线构成面的形式,对Dot-peenDPM进行图像识别。采用本方案,可以有效应对DPM码中的Dot-peen类型,并对金属反光、光照不均匀、噪点复杂等情况有很好的解码效果;此外,本方案对Dot-peen DPM码中存在部分点污损的情况有很好的容忍度,可以精确识别带有污渍的Dot-peen DPM码。
以上结合附图对本发明的实施方式作了详细说明,但本发明不限于所描述的实施方式。对于本领域的技术人员而言,在不脱离本发明原理和精神的情况下,对这些实施方式进行多种变化、修改、替换和变型,仍落入本发明的保护范围内。

Claims (6)

1.一种Dot-peen DPM码的识别方法,其特征在于,包括步骤:
通过极值检测获取图像中的点;
通过检测出的点的位置进行主方向确认,所述主方向为二维码边缘L线的两个方向;
检测定位主方向上的L线;
通过迭代透视变换生成采样网络;
对透视变换后的所有点坐标进行采样及二值化处理;
其中,通过极值检测获取图像中的点具体为:
设I、J是定义于同一离散域D上的两张灰度图像,I为掩模图像,J为标识图像,且
Figure FDA0003592576160000011
其中p表示图像上的任意像素点,J(p)、I(p)是p点在图像J、I中的灰度值,则用J对I做形态学重建表示为:
Figure FDA0003592576160000012
其中,
Figure FDA0003592576160000013
Figure FDA0003592576160000014
表示形态学膨胀操作,B为尺度为1的结构元素,∧表示逐点比较取最小值;
采用形态学重建,标识图像J满足:
Figure FDA0003592576160000015
形态学重建表示为ρI(J)=ρI(I-h),其中h为预设常数;
用于点检测的公式为:Sh(I)=I-ρI(I-h),S表示得到的标注区域极值点的图像;
其中,步骤检测定位主方向上的L线,包括:
剔除与主方向之间的夹角大于预设值的连线,获得θ1方向的连线集合L′1、θ2方向的连线集合L′2
将L′1的平均长度作为θ1方向的模块宽度s1,将L′2的平均长度作为θ2方向的模块宽度s2
采用Hough变换得到主方向的L线。
2.根据权利要求1所述的Dot-peenDPM码的识别方法,其特征在于,步骤通过检测出的点的位置进行主方向确认,包括:
对于每个点,找到与其曼哈顿距离最小的两个相邻点;
将所有相邻点连接起来,形成连通分量;
统计所有相邻点连线的方向,将数量最多的两个方向作为条码的两个主方向θ12
3.根据权利要求2所述的Dot-peenDPM码的识别方法,其特征在于,步骤通过迭代透视变换生成采样网络,包括:
由L线交点向两个主方向延伸,获取L线上所有点的模板坐标;
建立L线上点的实际坐标及模板坐标之间的透视变换公式;
将L线上的点的实际坐标及模板坐标代入透视变换公式,求得透视变换系数;
从L线向条码区域内部延伸,获取最靠近L的线段Lin-1所有模板坐标,通过透视变换计算可以得到Lin-1的所有模板坐标对应的伪实际坐标;
设置搜索半径,若在伪实际坐标的搜索半径内存在实际的图像检测点,则建立Lin-1上的实际坐标与模板坐标的映射对;
将L、Lin-1上所有的实际坐标、模板坐标映射对代入透视变换公式,更新透视变换系数;
从Lin-1向条码区域内部延伸,重复直至遍历DPM上所有检测点,生成采样网络。
4.一种Dot-peen DPM码的识别装置,其特征在于,包括:
点检测模块:通过极值检测获取图像中的点;
主方向模块:通过检测出的点的位置进行主方向确认,所述主方向为二维码边缘L线的两个方向;
L线模块:检测定位主方向上的L线;
迭代模块:通过迭代透视变换生成采样网络;
二值化模块:对透视变换后的所有点坐标进行采样及二值化处理;
其中,点检测模块中,通过极值检测获取图像中的点具体为:
设I、J是定义于同一离散域D上的两张灰度图像,I为掩模图像,J为标识图像,且
Figure FDA0003592576160000021
其中p表示图像上的任意像素点,J(p)、I(p)是p点在图像J、I中的灰度值,则用J对I做形态学重建表示为:
Figure FDA0003592576160000022
其中,
Figure FDA0003592576160000023
Figure FDA0003592576160000024
表示形态学膨胀操作,B为尺度为1的结构元素,∧表示逐点比较取最小值;
采用形态学重建,标识图像J满足:
Figure FDA0003592576160000031
形态学重建表示为ρI(J)=ρI(I-h),其中h为预设常数;
用于点检测的公式为:Sh(I)=I-ρI(I-h),S表示得到的标注区域极值点的图像;
其中,L线模块包括:
剔除单元:剔除与主方向之间的夹角大于预设值的连线,获得θ1方向的连线集合L′1、θ2方向的连线集合L′2
宽度单元:将L′1的平均长度作为θ1方向的模块宽度s1,将L′2的平均长度作为θ2方向的模块宽度s2
Hough变换单元:采用Hough变换得到主方向的L线。
5.根据权利要求4所述的Dot-peenDPM码的识别装置,其特征在于,主方向模块包括:
距离单元:对于每个点,找到与其曼哈顿距离最小的两个相邻点;
连通单元:将所有相邻点连接起来,形成连通分量;
统计单元:统计所有相邻点连线的方向,将数量最多的两个方向作为条码的两个主方向θ12
6.根据权利要求4所述的Dot-peenDPM码的识别装置,其特征在于,迭代模块包括:
延伸单元:由L线交点向两个主方向延伸,获取L线上所有点的模板坐标;
公式单元:建立L线上点的实际坐标及模板坐标之间的透视变换公式;
系数单元:将L线上的点的实际坐标及模板坐标代入透视变换公式,求得透视变换系数;
伪实际坐标单元:从L线向条码区域内部延伸,获取最靠近L的线段Lin-1所有模板坐标,通过透视变换计算可以得到Lin-1的所有模板坐标对应的伪实际坐标;
映射对单元:设置搜索半径,若在伪实际坐标的搜索半径内存在实际的图像检测点,则建立Lin-1上的实际坐标与模板坐标的映射对;
系数更新单元:将L、Lin-1上所有的实际坐标、模板坐标映射对代入透视变换公式,更新透视变换系数;
遍历单元:从Lin-1向条码区域内部延伸,重复直至遍历DPM上所有检测点,生成采样网络。
CN201910031689.7A 2019-01-14 2019-01-14 Dot-peen DPM码的识别方法及装置 Active CN109784121B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910031689.7A CN109784121B (zh) 2019-01-14 2019-01-14 Dot-peen DPM码的识别方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910031689.7A CN109784121B (zh) 2019-01-14 2019-01-14 Dot-peen DPM码的识别方法及装置

Publications (2)

Publication Number Publication Date
CN109784121A CN109784121A (zh) 2019-05-21
CN109784121B true CN109784121B (zh) 2022-09-23

Family

ID=66500544

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910031689.7A Active CN109784121B (zh) 2019-01-14 2019-01-14 Dot-peen DPM码的识别方法及装置

Country Status (1)

Country Link
CN (1) CN109784121B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116468060B (zh) * 2023-04-23 2024-01-19 创新奇智(上海)科技有限公司 一种dpm码处理方法、装置、存储介质及设备

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101765861A (zh) * 2007-03-08 2010-06-30 微扫描系统公司 二维码的检测和分割
CN102096795A (zh) * 2010-11-25 2011-06-15 西北工业大学 磨损二维条码图像识别方法
CN102930268A (zh) * 2012-08-31 2013-02-13 西北工业大学 一种用于污染及多视角情况下DataMatrix码的精确定位方法
KR101552202B1 (ko) * 2015-06-25 2015-09-10 주식회사 아이엔솔루션 포터블 dpm 바코드 스캐너
CN104915951A (zh) * 2015-05-08 2015-09-16 天津大学 一种点刻式dpm二维码区域定位方法
CN105654015A (zh) * 2014-11-18 2016-06-08 康耐视公司 解码二维矩阵符号的方法和系统
US9792477B1 (en) * 2016-05-16 2017-10-17 Symbol Technologies, Llc System for, and method of, controlling illumination of direct part marking (DPM) targets to be read by image capture
CN107301368A (zh) * 2017-06-28 2017-10-27 昂纳自动化技术(深圳)有限公司 一种DataMatrix二维码的识别方法
CN107609451A (zh) * 2017-09-14 2018-01-19 斯坦德机器人(深圳)有限公司 一种基于二维码的高精度视觉定位方法及系统
CN109086644A (zh) * 2018-07-27 2018-12-25 广东奥普特科技股份有限公司 亚像素精度的DataMatrix二维码的定位方法
CN109190434A (zh) * 2018-06-25 2019-01-11 广东工业大学 一种基于亚像素级角点检测的条码识别算法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4708186B2 (ja) * 2003-05-02 2011-06-22 豊 木内 2次元コード解読プログラム
CN103530590B (zh) * 2013-10-19 2016-02-24 高韬 Dpm二维码识别系统
CN103593695B (zh) * 2013-11-15 2016-08-10 天津大学 一种dpm二维码区域定位的方法
US10534970B2 (en) * 2014-12-24 2020-01-14 Datalogic Ip Tech S.R.L. System and method for reading direct part marking (DPM) codes on objects

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101765861A (zh) * 2007-03-08 2010-06-30 微扫描系统公司 二维码的检测和分割
CN102096795A (zh) * 2010-11-25 2011-06-15 西北工业大学 磨损二维条码图像识别方法
CN102930268A (zh) * 2012-08-31 2013-02-13 西北工业大学 一种用于污染及多视角情况下DataMatrix码的精确定位方法
CN105654015A (zh) * 2014-11-18 2016-06-08 康耐视公司 解码二维矩阵符号的方法和系统
CN104915951A (zh) * 2015-05-08 2015-09-16 天津大学 一种点刻式dpm二维码区域定位方法
KR101552202B1 (ko) * 2015-06-25 2015-09-10 주식회사 아이엔솔루션 포터블 dpm 바코드 스캐너
US9792477B1 (en) * 2016-05-16 2017-10-17 Symbol Technologies, Llc System for, and method of, controlling illumination of direct part marking (DPM) targets to be read by image capture
CN107301368A (zh) * 2017-06-28 2017-10-27 昂纳自动化技术(深圳)有限公司 一种DataMatrix二维码的识别方法
CN107609451A (zh) * 2017-09-14 2018-01-19 斯坦德机器人(深圳)有限公司 一种基于二维码的高精度视觉定位方法及系统
CN109190434A (zh) * 2018-06-25 2019-01-11 广东工业大学 一种基于亚像素级角点检测的条码识别算法
CN109086644A (zh) * 2018-07-27 2018-12-25 广东奥普特科技股份有限公司 亚像素精度的DataMatrix二维码的定位方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DPM码快速定位关键技术研究;刘志洪梅芳;《浙江工业大学学报》;20160225;第44卷(第1期);第28-33页 *
Yang Liu ; Zihuai Lin ; Ping Hao ; Ning Hong.Rate Flexible Coded Digital Phase Modulation (DPM) over Rings.《2014 8th International Conference on Future Generation Communication and Networking》.2015,第21-24页. *
复杂光照下DPM图像自适应多阈值分割方法研究;王娟 王萍 刘敏;《计算机工程与应用》;20170510;第194-200页 *

Also Published As

Publication number Publication date
CN109784121A (zh) 2019-05-21

Similar Documents

Publication Publication Date Title
CN108388896B (zh) 一种基于动态时序卷积神经网络的车牌识别方法
CN107507208B (zh) 一种基于轮廓上曲率估计的图像特征点提取方法
CN110866924A (zh) 一种线结构光中心线提取方法及存储介质
CN105654015A (zh) 解码二维矩阵符号的方法和系统
CN110647795A (zh) 一种表格识别方法
CN112560538B (zh) 一种根据图像冗余信息快速定位受损qr码的方法
CN102169581A (zh) 一种基于特征向量的快速高精度鲁棒性匹配方法
CN105894521A (zh) 基于高斯拟合的亚像素边缘检测方法
CN115170669A (zh) 基于边缘特征点集配准的识别定位方法及系统、存储介质
CN108960280B (zh) 一种图片相似度检测方法及系统
US11030478B1 (en) System and method for correspondence map determination
CN111582000B (zh) 一种条形码定位方法、装置及相关设备
CN109190625A (zh) 一种大角度透视变形的集装箱箱号识别方法
CN113903024A (zh) 一种手写票据数值信息识别方法、系统、介质及装置
CN113592923A (zh) 一种基于深度局部特征匹配的批图像配准方法
CN111368573A (zh) 一种基于几何特征约束的定位方法
CN109784328B (zh) 定位条码的方法、终端及计算机可读存储介质
CN109784121B (zh) Dot-peen DPM码的识别方法及装置
CN115713556A (zh) 二维码定位方法、装置及计算机可读存储介质
CN113793309B (zh) 一种基于形态学特征的亚像素级椭圆检测方法
CN114612412A (zh) 三维点云数据的处理方法及其用途、电子设备、存储介质
CN111931537B (zh) 一种颗粒状qr二维码定位方法
CN111160142B (zh) 一种基于数值预测回归模型的证件票据定位检测方法
CN110378167B (zh) 一种基于深度学习的条码图像补正方法
CN116740062A (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