CN105335973A - 运用于带钢加工生产线的视觉处理方法 - Google Patents
运用于带钢加工生产线的视觉处理方法 Download PDFInfo
- Publication number
- CN105335973A CN105335973A CN201510684833.9A CN201510684833A CN105335973A CN 105335973 A CN105335973 A CN 105335973A CN 201510684833 A CN201510684833 A CN 201510684833A CN 105335973 A CN105335973 A CN 105335973A
- Authority
- CN
- China
- Prior art keywords
- point
- target image
- template
- profile
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
- G06T7/001—Industrial image inspection using an image reference approach
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30136—Metal
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明提供了一种运用于带钢加工生产线的视觉处理方法,其针对搬运堆垛过程的带钢粗略定位,提取带钢上标志物轮廓信息并转化为极坐标下的表示形式,制定匹配规则完成匹配和定位。针对落料加工过程中,带钢待切割前的精确定位,建立带旋转的图像金字塔,并充分利用已知模板的形状信息,计算区域内满足条件点的梯度幅值与方向信息,逐层搜索各金字塔图层特征点,最终返回精确位置平移以及旋转坐标信息;针对匹配过程中遇到的由于相机视野有限带来的目标图像残缺问题,自定义残缺轮廓的类曲率特征,实现了残缺目标的匹配和定位。针对带钢切割后分离效果的检测,对图像进行自适应阈值区域分割,完成了分离检测。
Description
技术领域
本发明涉及计算机视觉领域,具体地,涉及一种运用于带钢加工生产线的视觉处理方法。
背景技术
目前国内大多数钢板切割过程的钢板定位主要靠以下方式进行:人工目测(精度需求不高时)、电火花碰边、辅助工装固定、标尺测量等方式,这些方式存在精度不高、拆装工作量大、自动化程度低、时效性不高的问题。
目前基于形状的图像匹配领域有了很多研究,主要分为两个大类。
一是基于边缘的形状描述匹配:将形状用周长、主轴、形状数、傅立叶描述、小波描述、尺度空间等一定的特征重新数学化表达,定义相关的匹配函数完成匹配。目前在基于形状的图像匹配中,傅里叶描述方法较为出色且易于实现,但轮廓上较尖锐的变化被映射到全部系数中而无法侦测到,且只能实现定性匹配,而不利于完成模板的精确定位。
二是基于区域的形状描述匹配:利用图像的区域面积、重心、形变参数、偏心率、各种类型的不变矩等来完成匹配。简单的区域形状描述符实现简单,有很明确的物理意义,但有时候在区分不同形状的时候能力较弱。将图像灰度转化为数学空间的各种不变矩也是常用的方法,它能很好的识别形状特征,但一些高阶矩阵的运算量太大而不足以满足时效性要求,而且对噪声很敏感。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种运用于带钢加工生产线的视觉处理方法。
根据本发明提供的一种运用于带钢加工生产线的视觉处理方法,其特征在于,包括如下任一个或任多个步骤:
搬运堆垛定位步骤,记为步骤S1,具体地:根据在带钢加工生产的搬运堆垛过程中获取的目标图像,对带钢进行定位,其中,在带钢加工生产的搬运堆垛过程中获取的目标图像记为第一目标图像;
落料加工定位步骤,记为步骤S2,具体地:根据在带钢加工生产的落料加工过程中获取的目标图像,对带钢进行定位,其中,在带钢加工生产的落料加工过程中获取的目标图像记为第二目标图像;
匹配定位步骤,记为步骤S3,具体地:针对目标图像为残缺图像的情况,进行残缺图像的匹配与定位。
优选地,所述步骤S1包括如下步骤:
步骤S1.1、提取带钢原始图像上目标图形的轮廓特征;
步骤S1.2、进行目标图形的轮廓特征量的极坐标变换,包括如下步骤:
步骤S1.2.1、利用目标图形的轮廓边缘上每一个像素点的直角坐标(x,y),得到轮廓几何中心点的直角坐标(x0,y0);
步骤S1.2.2、将像素点直角坐标(x,y)转换为极坐标(r,θ),其中,r为极径,θ为极角;
步骤S1.3、将基于轮廓点极坐标信息的模板与第一目标图像进行匹配,包括步骤:
步骤S1.3.1、分别寻找对应极径r为最大值时,模板轮廓上的轮廓点P、第一目标图像中目标轮廓上的轮廓点Pd,计算尺度比例L,其中,尺度比例L的计算方式为:假设模板轮廓上的轮廓点P对应的最大极径值为R1,目标轮廓上的轮廓点Pd对应的最大极径值为R2,则L=R2/R1;
步骤S1.3.2、用尺度比例L对目标轮廓极坐标信息进行归一化操作,从极径r最大时对应的轮廓点开始,以对应模板轮廓与目标轮廓各自尺寸的步长逆时针遍历取点,得到各组对应点,然后计算各组对应点极径差值的和值,将该和值除以步数作为平均误差;
步骤S1.3.3、若平均误差小于设置阈值,则认为模板轮廓与目标轮廓匹配,进入步骤S1.3.4继续执行,否则,则认为不匹配,将当前的模板轮廓替换为不同的模板轮廓,进入步骤S1.3.1继续执行;
步骤S1.3.4、将目标轮廓的轮廓中心点与模板轮廓的轮廓中心点之间的差值作为偏移量,利用各自对应极径区域极值点的极角差值平均值作为旋转角度,利用各自对应极径区域极值点的比值平均值作为尺度变化量。
优选地,所述步骤S2包括如下步骤:
步骤S2.1、利用高斯内核卷积带钢原始图像获得高斯金字塔;
步骤S2.2、利用Sobel算子卷积高斯金字塔各层图像获得每层的梯度幅值图,将梯度方向保存成灰度图,并计算灰度图中每一个像素点的归一化方向向量;
步骤S2.3、基于形状梯度信息的模板与第二目标图像的匹配,包括步骤:
步骤S2.3.1、在对应于最小分辨率的高斯金字塔顶层图层,将模板子图在第二目标图像中不断平移直至遍历第二目标图像,根据如下公式计算匹配分数S1:
其中,m、n分别为第二目标图像的长、宽,Th,k为模板子图中像素点的归一化方向向量,Dh,k为第二目标图像中对应像素点的归一化方向向量;下标h,k表示像素点在当前模板子图内的相对位置;
将第二目标图像以角度ψ为步长,以第二目标图像的几何中心为旋转中心不断旋转第二目标图像,利用二次线性插值生成新的图像矩阵,每次旋转过后都计算一次当前S1取最大值时模板子图在第二目标图像中的平移量(Sx,Sy),旋转第二目标图像直至360度,返回S1最大时对应的旋转角度ψM,记录当前平移量(SxM,SyM);Sx、SxM表示沿图像宽度方向的平移量,Sy、SyM表示沿图像高度方向的平移量;
步骤S2.3.2、对于下一层高斯金字塔图层,以小于角度ψ的步长在ψM邻域内旋转第二目标图像,获得新的图像矩阵,每次旋转结束后,在(SxM,SyM)为中心的模板子图邻域平移模板子图,利用自适应阈值分割当前层梯度幅值图像提取边缘点,具体根据如下公式计算匹配分数S2:
其中,(x,y)为边缘点坐标,T为模板子图中的归一化方向向量,D为第二目标图像中对应点的归一化方向向量,N表示第二目标图像中模板子图大小的区域中像素点梯度幅值大于设定阈值的点的数目
计算并返回S2最大时对应的旋转角度和平移量,更新ψM和(SxM,SyM),进入高斯金字塔下一层;
其中,不断重复执行步骤S2.3.2直至高斯金字塔底层,返回最终旋转角度和平移量完成匹配。
优选地,所述步骤S3包括如下步骤:
步骤S3.1、获取目标图像中图形的轮廓;
步骤S3.2、针对目标图像为残缺图像缺失信息而只能获得部分信息的情况,将模板轮廓与目标残缺轮廓的类曲率特征保存在一维数组中,寻找目标残缺轮廓的类曲率数组在模板轮廓类曲率数组中最相似的段落,完成残缺图像匹配;
步骤S3.3、计算模板轮廓中心点,建立模板轮廓上像素点与由模板轮廓所有轮廓点的坐标加权平均获得的中心点之间的对应关系,以该对应关系计算出目标残缺轮廓虚拟中心点,完成残缺图像定位。
优选地,所述步骤S1还包括如下步骤:
步骤S1.2.3、处理轮廓中心点在轮廓外部情况,具体地:
步骤S1.2.3通过以下步骤完成:
-在没有对应极径的极角范围,r取值为0;
-在没有对应极径的极角范围之外的其他角度范围,如果存在多个r值,则进行区别保存。
优选地,步骤S2.1中所述的高斯内核为:
优选地,步骤S2.2中所述归一化方向向量通过以下步骤获得:
a.利用Sobel算子对灰度图进行滤波,分别获得灰度图的高度方向Y的梯度Dy、宽度方向X的梯度Dx;
b.根据如下公式建立归一化方向向量:
EdgeDir=atan-1(Dy/Dx)
e=(COS(EdgeDir),SIN(EdgeDir))
其中,EdgeAmp为像素点的梯度方向,EdgeDir为像素点的梯度幅值,e为像素点的单位方向向量。
优选地,步骤S3.2中轮廓的类曲率特征通过以下步骤完成:
a.目标残缺轮廓视为连续轮廓,预处理目标残缺轮廓,截除首尾部分轮廓点;
b.从目标残缺轮廓上某一点开始逆时针遍历轮廓,取3个点,根据如下公式计算类曲率K作为类曲率特征:
其中,(xi-1,yi-1)、(xi,yi)、(xi+1,yi+1)分别为目标残缺轮廓的首端点、中间点、末端点的坐标,(xi,o,yi,o)为首端点与中间点的中垂线相交于中间点与末端点的中垂线的交点坐标。
优选地,还包括如下步骤:
分离检测步骤,记为步骤S4,具体为:在带钢切割后采集第三目标图像,进行分离检测。
优选地,所述步骤SS4包括如下步骤:
步骤S4.1、利用自适应阈值分割分离图像钢板和切割后产生的空洞部分,将空洞部分作为感兴趣区域;
步骤S4.2、根据如下公式计算匹配指标M的值,判断与设定的阈值是否一致,若一致则判定切割完成,否则判定为未完成;
其中,M为匹配指标,L为感兴趣区域周长,A为感兴趣区域面积。
与现有技术相比,本发明具有如下的有益效果:
1、本发明能适应带钢加工生产线高亮、多粉尘的工业环境,能够解决实际生产中由于相机视野有限导致的目标图像残缺问题,完成完整的带钢加工生产任务,定位误差限定在2个像素以内,并且能满足实时性要求。
2、带钢加工生产中需要运用的视觉处理的部分主要为以下几个:一、带钢搬运到传送带上的搬运堆垛过程中的粗略定位;二、带钢搬运到传送带上后,仍在x和y方向存在一定的误差和旋转角度误差,在切割前,需要对钢板进行精确定位;三、带钢在切割完成后的分离检测过程,对此,本发明提出了具体视觉处理方法,能够应对噪声,局部形变,光照变化等鲁邦性问题,定位效果良好。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为搬运堆垛定位步骤(即带钢粗略定位)的流程示意图。
图2为落料加工定位步骤(即带钢精确定位)的流程示意图。
图3为匹配定位步骤(即残缺图像匹配定位)的流程示意图。
图4为分离检测步骤的流程示意图。
图5为完整图像的举例示意图。
图6为残缺图像的举例示意图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
根据本发明提供的一种运用于带钢加工生产线的视觉处理方法,包括如下任一个或任多个步骤:
第一次定位步骤,记为步骤S1,具体地:在带钢加工生产的搬运堆垛过程中获取第一目标图像,对带钢进行第一次定位,记为粗定位;
第二次定位步骤,记为步骤S2,具体地:在带钢加工生产的落料加工过程中获取第二目标图像,对带钢进行第二次定位,记为细定位;
匹配定位步骤,记为步骤S3,具体地:针对第一目标图像和/或第二目标图像为残缺图像的情况,进行残缺图像的匹配与定位;由于相机拍摄视野有限,拍摄到的图片中待切割形状不能完整获取,例如,如图5所示图为完整图像,图6所示为残缺图像。
所述步骤S1包括如下步骤:
步骤S1.1、通过canny算子提取带钢原始图像上目标图形的轮廓特征;
步骤S1.2、进行目标图形的轮廓特征量的极坐标变换,包括如下步骤:
步骤S1.2.1、利用目标图形的轮廓边缘上每一个像素点的直角坐标(x,y),得到轮廓几何中心点的直角坐标(x0,y0);
步骤S1.2.2、将像素点直角坐标(x,y)转换为极坐标(r,θ),其中,r为极径,θ为极角;公式如下:
步骤S1.3、进行基于轮廓点极坐标信息的模板与第一目标图像的匹配,其中:模板的获取有两种可能情况,一是直接采用目标切割图形的标准CAD图作为模板,二是用相机拍摄切割好的标准带钢图像,截取相机拍摄图像中完整包括了目标切割图形的矩形区域作为模板图片;第一目标图像即是待匹配图像,也就是实际运作过程中相机拍摄的图像,目标图形即第一目标图像中待切割部分图形。所述进行基于轮廓点极坐标信息的模板与第一目标图像的匹配,包括步骤:
步骤S1.3.1、分别寻找对应极径r为最大值时,模板轮廓上的轮廓点P、目标轮廓上的轮廓点Pd,计算尺度比例L,其中,尺度比例L的计算方式为:假设模板轮廓上的轮廓点P对应的最大极径值为R1,目标轮廓上的轮廓点Pd对应的最大极径值为R2,则L=R2/R1;其中,目标轮廓的获取方法为:第一目标图像通过提取轮廓算子获得了轮廓,去除其中由于杂质或者带钢上水雾导致的部分干扰轮廓后获得的图形轮廓边缘作为目标轮廓,即待匹配轮廓;
步骤S1.3.2、用尺度比例L对目标轮廓极坐标信息进行归一化操作,从极径r最大时对应的轮廓点开始,以对应模板轮廓与目标轮廓各自尺寸的步长逆时针遍历取点,得到各组对应点,然后计算各组对应点极径差值的和值,将该和值除以步数作为平均误差;其中,如果相机在拍摄图像时在竖直方向上有平移,则拍摄到的图像相对大小不同,图像中待切割图形的相对大小也会发生变化,所获得的轮廓点数多少也发生了变化,这里的尺寸即为这种相对大小,具体计算中可采用轮廓点的数目作为尺寸大小的判定指标;
步骤S1.3.3、若平均误差小于设置阈值,则认为模板轮廓与目标轮廓匹配,进入步骤S1.3.4继续执行,否则,则认为不匹配,将当前的模板轮廓替换为不同的模板轮廓,进入步骤S1.3.1继续执行;
步骤S1.3.4、将目标轮廓的轮廓中心点与模板轮廓的轮廓中心点之间的差值作为偏移量,利用各自对应极径区域极值点的极角差值平均值作为旋转角度,利用各自对应极径区域极值点的比值平均值作为尺度变化量;其中,将轮廓点坐标表示为极径和极角时,有时候轮廓中心点会处于整个轮廓外部,在某一个极角范围不存在对应的极径值,所以把极角存在对应极径值的极角范围叫做极径区域;所述比值为目标轮廓某上某一点的极径值除以模板轮廓上某一点的极径值获得。
所述步骤S2包括如下步骤:
步骤S2.1、利用高斯内核卷积带钢原始图像获得高斯金字塔;
步骤S2.2、利用Sobel算子卷积高斯金字塔各层图像获得每层的梯度幅值图,将梯度方向保存成灰度图,并计算灰度图中每一个像素点的归一化方向向量;
步骤S2.3、进行基于形状梯度信息的模板与第二目标图像的匹配,其中:基于形状梯度信息的模板的获取有两种可能情况,一是直接采用目标切割图形的标准CAD图作为模板,二是用相机拍摄切割好的标准带钢图像,截取相机拍摄图像中完整包括了目标切割图形的矩形区域作为模板图片,其中,所述进行基于形状梯度信息的模板与第二目标图像的匹配,包括步骤:
步骤S2.3.1、在对应于最小分辨率的高斯金字塔顶层图层,利用归一化交叉相关的基于灰度值的匹配方法,将模板子图在第二目标图像中不断平移直至遍历第二目标图像,根据如下公式计算匹配分数S1:
其中,m、n分别为第二目标图像的长、宽,Th,k为模板子图中像素点的归一化方向向量,Dh,k为第二目标图像中对应像素点的归一化方向向量;下标h,k表示像素点在当前模板子图内的相对位置;其中,模板子图指的是:在目标图像中选取模板大小的图形区域,称为模板子图;
将第二目标图像以一定的角度ψ为步长,以第二目标图像的几何中心为旋转中心不断旋转第二目标图像,利用二次线性插值生成新的图像矩阵,每次旋转过后都计算一次当前S1取最大值时模板子图在第二目标图像中的平移量(Sx,Sy),旋转第二目标图像直至360度,返回S1最大时对应的旋转角度ψM,记录当前平移量(SxM,SyM);Sx、SxM表示沿图像宽度方向的平移量,Sy、SyM表示沿图像高度方向的平移量;
步骤S2.3.2、对于下一层高斯金字塔图层,以小于角度ψ的步长在ψM邻域内旋转第二目标图像,获得新的图像矩阵,每次旋转结束后,在(SxM,SyM)为中心的模板子图邻域平移模板子图,利用自适应阈值分割当前层梯度幅值图像提取边缘点,具体根据如下公式计算匹配分数S2:
其中,(x,y)为边缘点坐标,T为模板子图中的归一化方向向量,D为第二目标图像中对应点的归一化方向向量,N表示第二目标图像中模板子图大小的区域中像素点梯度幅值大于设定阈值的点的数目
计算并返回S2最大时对应的旋转角度和平移量,更新ψM和(SxM,SyM),进入高斯金字塔下一层;
其中,不断重复执行步骤S2.3.2直至高斯金字塔底层(分辨率最大),返回最终旋转角度和平移量完成匹配。
所述步骤S3包括如下步骤:
步骤S3.1、通过canny边缘检测算法获取第一目标图像和/或第二目标图像中图形的轮廓;其中,整个带钢加工生产线分为3个加工部分:搬运堆垛过程、落料加工过程、分离检测过程;各个加工部分分别有不同精度的相机采集目标图像,各相机拍摄到的图像都称为各自部分的目标图像,其中,搬运堆垛过程中相机采集到的目标图像记为第一图标图像,落料加工过程中相机采集到的目标图像记为第二图标图像,分离检测过程中相机采集到的目标图像记为第三图标图像;
步骤S3.2、针对目标图像为残缺图像缺失信息而只能获得部分信息的情况,将模板轮廓与目标残缺轮廓的类曲率特征保存在一维数组中,寻找目标残缺轮廓的类曲率数组在模板轮廓类曲率数组中最相似的段落,完成残缺图像匹配;
步骤S3.3、计算模板轮廓中心点,建立模板轮廓上像素点与由模板轮廓所有轮廓点的坐标加权平均获得的中心点之间的对应关系,以该对应关系计算出目标残缺轮廓虚拟中心点,完成残缺图像定位。
5.根据权利要求2所述的运用于带钢加工生产线的视觉处理方法,其特征在于,所述步骤S1还包括如下步骤:
步骤S1.2.3、处理轮廓中心点在轮廓外部情况,具体地:
步骤S1.2.3通过以下步骤完成:
-在没有对应极径的极角范围,r取值为0;
-在没有对应极径的极角范围之外的其他角度范围,如果存在多个r值,则进行区别保存。
步骤S2.1中所述的高斯内核为:
高斯内核与当前层图像卷积完成后,删除偶数行列,获得上一层图像。
所述步骤S2.2中所述归一化方向向量通过以下步骤获得:
a.利用Sobel算子对灰度图进行滤波,分别获得灰度图的高度方向Y的梯度Dy、宽度方向X的梯度Dx;
b.根据如下公式建立归一化方向向量:
EdgeDir=atan-1(Dy/Dx)
e=(COS(EdgeDir),SIN(EdgeDir))
其中,EdgeAmp为像素点的梯度方向,EdgeDir为像素点的梯度幅值,e为像素点的单位方向向量。
9.根据权利要求4所述的运用于带钢加工生产线的视觉处理方法,其特征在于,步骤S3.2中轮廓的类曲率特征通过以下步骤完成:
a.带钢加工过程中获得的残缺轮廓是由于相机视野有限导致的,目标残缺轮廓可视为连续轮廓,预处理目标残缺轮廓,截除首尾部分轮廓点以减小误差;
b.从目标残缺轮廓上某一点开始逆时针遍历轮廓,以适应轮廓点总数目的间隔取3个点,根据如下公式计算类曲率K作为类曲率特征:
其中,(xi-1,yi-1)、(xi,yi)、(xi+1,yi+1)分别为目标残缺轮廓的首端点、中间点、末端点的坐标,(xi,o,yi,o)为首端点与中间点的中垂线相交于中间点与末端点的中垂线的交点坐标。
所述运用于带钢加工生产线的视觉处理方法,还包括如下步骤:
分离检测步骤,记为步骤S4,具体为:在带钢切割后采集第三目标图像,进行分离检测;
所述步骤SS4包括如下步骤:
步骤S4.1、利用自适应阈值分割分离图像钢板和切割后产生的空洞部分,将空洞部分作为感兴趣区域;
步骤S4.2、根据如下公式计算匹配指标M的值,判断与设定的阈值是否一致,若一致则判定切割完成,否则判定为未完成;
其中,M为匹配指标,L为感兴趣区域周长,A为感兴趣区域面积。
在优选例中,所述步骤S3.2包括如下步骤:
a.设定初始比较位置为模板轮廓类曲率数组第一个下标;
b.计算目标残缺轮廓的类曲率数组与模板轮廓类曲率数组对应长度段的欧式距离之和,移动初始比较位置至下一个下标,直至初始比较位置遍历整个模板轮廓类曲率数组,将欧式距离之和最小时对应的下标作为最佳匹配起始点,欧式距离除以预处理后的目标残缺轮廓点数目作为匹配分数,该值越小,认为匹配程度越高;
c.当匹配分数小于设定阈值时,认为目标残缺轮廓属于模板轮廓,完成匹配。
在优选例中,所述步骤S3.3包括如下步骤:
a.计算模板轮廓中心点O1;
b.从模板轮廓上某一点开始逆时针遍历轮廓,以一定的间隔获取模板轮廓上3个点A1、B1、C1和目标残缺轮廓上3个点A2、B2、C2,分别计算点A1、B1、C1到模板轮廓中心的距离D1,D2、D3;
c.分别以A2、B2为圆心,D1,D2为半径画圆,保留两交点中的合理交点O2,计算O2到以C2为圆心,D3为半径的圆的距离,若该距离小于设定阈值,则将点O2作为计算目标残缺轮廓虚拟中心点的可用点,计算并保存向量O2-A2与向量01-A1的夹角,为计算旋转角做准备;
d.A1取模板轮廓上下一个点,重复整个过程直至A1遍历模板轮廓上所有点,计算步骤c中所有可用点的平均值,作为目标残缺轮廓虚拟中心点,计算步骤c中所有向量夹角的平均值,作为目标残缺轮廓相对于模板轮廓的旋转角度,完成定位。
所述切割是否完成,包含以下任一种情况:
情况a.正常切割;
情况b.切断部分钢板掉落时某一边缘与带钢仍有连结,切断钢板区域大部分掉落,第三目标图像中呈现小部分与钢板背景色相同区域,即空洞区域小部分凹陷;
情况c.切断部分钢板掉落时多边缘与带钢仍有连结,第三目标图像中呈现大部分与钢板背景色相同区域,即空洞区域大部分凹陷,成弯月形状态;
情况d.钢板未被切断,仅切出几条缝隙,空洞区域为多条细长型区域;
情况e.钢板完全掉落,但检测过程中拍摄第三目标图像时,空洞区域存在同带钢灰度相近的色块干扰(传送带底部有已掉落钢板或者其他杂物反光导致),这种情况会被误检为未切割完成,需要加入进一步的判定;具体地,传送带底部钢板或者其他杂物反光部分与完整切割边缘存在一条细长的暗色区域,直接采用canny算子检测整个区域将不会把该暗色区域作为边缘,而通过以下方法检测此区域是否存在,即可判定切割是否完成,步骤如下:
步骤I.确定完整切割情况的空洞边缘;
步骤II.对于待检测图像,将检测区域缩小到空洞边缘位置比边缘略大的区域,同样采用canny算子检测边缘,若检测到细长暗色区域,则切割是正常的;若不能检测到细长暗色区域,则认为切割没有完成,切断部分仍与带钢相连。
下面对本发明的一个优选的具体实施方式进行说明。
原始图像分辨率为6576*4384(对于第一部分带钢的粗略定位过程,相机精度不需要很高,但此实例仍采用高分辨率的图片来检测算法的精度);
将原始图像中创建模板的形状部分以(1571,1683)为旋转中心顺时针旋转5度后,平移(200,100)个像素作为目标图像。
搬运堆垛过程中带钢的粗略定位过程:
通过canny算子提取带钢上目标图形的轮廓特征,低阈值10,高阈值25,轮廓点数为3740;
其中Dx为X方向卷积邻域,Dy为Y方向卷积邻域;
利用轮廓边缘上每一个像素点的坐标(x,y),计算轮廓中心点(x0,y0);
模板轮廓中心点为:(1570.840,1683.259);
目标轮廓中心点为:(1771.067,1783.394);
根据公式:
其中r为极径,θ为极角;
将xy平面的像素点坐标转换到极坐标平面;
所用实例不存在r取值为0的极角范围,也不存在极角对应多个r值的情况;
模板轮廓r最大值为675.766,对应轮廓点坐标为(1377,1032);
目标轮廓r最大值为676.549,对应轮廓点坐标为(1637,1120);
尺度比例L=1.00115,可视为1,小数点部分视为误差,因此不进行归一化操作;
从r最大时对应轮廓点开始,以对应模板轮廓与目标轮廓各自尺寸的步长逆时针遍历取点,计算对应点极径差值的和,除以步数作为平均误差;
平均误差阈值设定为2,本例平均误差为0.594,认为模板轮廓与目标轮廓匹配;
目标轮廓的轮廓中心点与模板轮廓中心点差值为(200.277,100.135),取整(200,100)作为偏移量,旋转角度平均值为-5.19522度,即目标轮廓相对于模板轮廓顺时针旋转了5.19522度,判定不存在尺度变化。
落料加工过程中的带钢精确定位过程;
利用高斯内核卷积原始图像获得图像高斯金字塔:
高斯内核与当前层图像卷积完成后,删除偶数行列,可获得上一层图像;
实际构建7层金字塔图像,金字塔顶层图像分辨率为102*68
利用Sobel算子卷积高斯金字塔各层图像获得每层的梯度幅值图,以一定规则将梯度方向保存成灰度图,并计算每一个像素点的归一化方向向量;
利用Sobel算子对图像进行滤波即可分别获得图像的Y,X方向梯度Dy,Dx;
根据公式:
EdgeDir=atan-1(Dy/Dx)
e=(COS(EdgeDir),SIN(EdgeDir))
其中EdgeAmp为像素点的梯度方向,EdgeDir为像素点的梯度幅值,e为像素点的单位方向向量;
可建立归一化方向向量。
模板子图在金字塔底层的分辨率大小为1272*1272的矩形区域,在金子塔顶层分辨率大小为20*20
将模板子图在目标图像中不断平移,根据公式
其中m,n分别为模板子图的长和宽,Th,k为模板子图中点的归一化方向向量,Dh,k为目标图像中对应点的归一化方向向量,S1为匹配分数;
将目标图像以一定的角度Ψ=1度为步长,以图像中心为旋转中心不断逆时针旋转目标图像,利用二次线性插值生成新的图像矩阵,每次旋转过后都计算一次当前S1取最大值时模板子图在目标图像中的平移量(Sx,Sy),旋转目标图像直至360度,返回S1最大时的旋转角度ψM=5度,当前平移量(SxM,SyM);
对于下一层高斯金字塔图层,以10分之1于上一层图像的旋转步长在ψM正负5个步长范围内内旋转目标图像,获得新的图像矩阵,每次旋转结束后,在(SxM,SyM)为中心的模板子图上下左右各大于模板子图当前层5个像素的区域平移模板子图,利用自适应阈值分割当前层梯度幅值图像提取边缘点,取得阈值为46,根据公式
其中(x,y)为边缘点坐标,T为模板子图中的归一化方向向量,D为目标图像中对应点的归一化方向向量;
计算并返回S2最大时对应的旋转角度和平移量,更新ψM和(SxM,SyM),进入高斯金字塔下一层;
不断重复步骤进入下一层直至高斯金字塔底层(分别率最大)
最终匹配分数S2为0.984111
返回最终旋转角度为顺时针旋转了5.032131度
模板子图从目标图像左上角平移至目标位置平移量为(1135,1147);
残缺图像的匹配与定位;
人为将原始图像右半部分截断,将构建的新区域作为目标图像区域。
通过canny边缘检测算法获取目标图像中图形的轮廓;
针对残缺图像缺失信息而只能获得部分信息的情况,设计基于类曲率特征量的匹配方法,进行必要的轮廓细化处理,以便计算轮廓的类曲率特征;
带钢加工过程中获得的残缺轮廓是由于相机视野有限导致的,目标残缺轮廓可视为连续轮廓,预处理目标残缺轮廓,截除首尾部分轮廓点以减小误差,本例去除首尾各20个轮廓点;
从轮廓上某一点开始逆时针遍历轮廓,以轮廓总数的100分之1为间隔取3个点,根据公式
其中(xi-1,yi-1),(xi,yi),(xi+1,yi+1)分别为首端点,中间点和末端点的坐标,(xi,o,yi,o)为首端点与中间点的中垂线和中间点与末端点的中垂线的交点坐标;
可计算类曲率K;
计算模板轮廓与目标残缺轮廓的类曲率特征,保存在一维数组cur[]中,寻找残缺轮廓的类曲率数组在模板轮廓类曲率数组中最相似的段落;
设定初始比较位置为模板轮廓类曲率数组第一个下标;
计算目标残缺轮廓的类曲率数组与模板轮廓类曲率数组对应长度段的欧式距离之和;
移动初始比较位置至下一个下标;
重复过程直至初始比较位置遍历整个模板轮廓类曲率数组,将欧式距离之和最小时对应的下标作为最佳匹配起始点,欧式距离除以预处理后的目标残缺轮廓点数目作为匹配分数,该值越小,认为匹配程度越高;
最终匹配分数为6.16892x10^-5,小于设定阈值3x10^-4,认为目标残缺轮廓属于模板轮廓,
该时刻匹配起始点为模板轮廓类曲率数组下标1341,记录该点,为下一步定位工作做准备;
计算模板轮廓中心点,建立模板轮廓上像素点与中心点对应关系,以相同的对应关系计算出残缺轮廓虚拟中心点;
模板轮廓中心点O1=(1570.840,1683.259);
从模板轮廓上某一点开始逆时针遍历轮廓,以150个像素点为间隔获取模板轮廓上3个点A1、B1、C1和目标残缺轮廓上3个点A2、B2、C2,分别计算点A1、B1、C1到模板轮廓中心的距离D1,D2、D3;
分别以A2、B2为圆心,D1,D2为半径画圆,保留两交点中的合理交点O2,计算O2到以C2为圆心,D3为半径的圆的距离,若该距离小于设定阈值,则将点O2作为计算目标残缺轮廓虚拟中心点的可用点,计算并保存向量O2-A2与向量01-A1的夹角,为计算旋转角做准备;
A1取模板轮廓上下一个点,重复整个过程直至A1遍历模板轮廓上所有点,计算所有可用点的平均值,作为目标残缺轮廓虚拟中心点,计算所有向量夹角的平均值,作为目标残缺轮廓相对于模板轮廓的旋转角度
最终轮廓虚拟中心点为(1771.299,1783.24)
模板轮廓旋转角度为-4.85117.即顺时针旋转4.85117度
带钢切割后的分离检测;
利用自适应阈值分割分离图像钢板和切割后产生的空洞部分,将空洞部分作为感兴趣区域;
根据公式
其中M为匹配指标,L为感兴趣区域周长,A为感兴趣区域面积;
计算各种可能的不同切割情况下的指标M的值,与完美切割情况下的指标M的值比较,判定切割是否完成。
1.正常切割情况下指标M的值为:1.60555
对于任意平移,旋转,尺度变化情况下的切割过程,认为M值在1.5至1.7范围内或者M的值为1.60555时为正常切割。
2.切断部分钢板掉落时某一边缘与带钢仍有连结,切断钢板区域大部分掉落,图片中呈现小部分与钢板背景色相同区域,即空洞区域小部分凹陷;
该情况下指标M的值为2.8868;
3.切断部分钢板掉落时多边缘与带钢仍有连结,图片中呈现大部分与钢板背景色相同区域,即空洞区域大部分凹陷,成弯月形状态;
该情况下指标M的值为8.08486;
4.钢板未被切断,仅切出几条缝隙,空洞区域为多条细长型区域;
该情况下指标M的值为227.002
5.钢板完全掉落,但检测过程中拍摄照片时,空洞区域存在同带钢灰度相近的色块干扰(传送带底部有已掉落钢板或者其他杂物反光导致),这种情况会被误检为未切割完成,需要加入进一步的判定;
传送带底部钢板或者其他杂物反光部分与完美切割边缘存在一条细长的暗色区域,直接采用canny算子检测整个区域将不会把该暗色区域作为边缘,而通过以下方法检测此区域是否存在,即可判定切割是否完成;
确定完整切割情况的空洞边缘;
对于待检测图像,将检测区域缩小到空洞边缘位置比边缘略大的区域,同样采用canny算子检测边缘,低阈值10,高阈值25,此时可检测到细长暗色区域,则切割是正常的。若不能检测到细长暗色区域,则认为切割不完成,切断部分仍与带钢相连。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (10)
1.一种运用于带钢加工生产线的视觉处理方法,其特征在于,包括如下任一个或任多个步骤:
搬运堆垛定位步骤,记为步骤S1,具体地:根据在带钢加工生产的搬运堆垛过程中获取的目标图像,对带钢进行定位,其中,在带钢加工生产的搬运堆垛过程中获取的目标图像记为第一目标图像;
落料加工定位步骤,记为步骤S2,具体地:根据在带钢加工生产的落料加工过程中获取的目标图像,对带钢进行定位,其中,在带钢加工生产的落料加工过程中获取的目标图像记为第二目标图像;
匹配定位步骤,记为步骤S3,具体地:针对目标图像为残缺图像的情况,进行残缺图像的匹配与定位。
2.根据权利要求1所述的运用于带钢加工生产线的视觉处理方法,其特征在于,所述步骤S1包括如下步骤:
步骤S1.1、提取带钢原始图像上目标图形的轮廓特征;
步骤S1.2、进行目标图形的轮廓特征量的极坐标变换,包括如下步骤:
步骤S1.2.1、利用目标图形的轮廓边缘上每一个像素点的直角坐标(x,y),得到轮廓几何中心点的直角坐标(x0,y0);
步骤S1.2.2、将像素点直角坐标(x,y)转换为极坐标(r,θ),其中,r为极径,θ为极角;
步骤S1.3、将基于轮廓点极坐标信息的模板与第一目标图像进行匹配,包括步骤:
步骤S1.3.1、分别寻找对应极径r为最大值时,模板轮廓上的轮廓点P、第一目标图像中目标轮廓上的轮廓点Pd,计算尺度比例L,其中,尺度比例L的计算方式为:假设模板轮廓上的轮廓点P对应的最大极径值为R1,目标轮廓上的轮廓点Pd对应的最大极径值为R2,则L=R2/R1;
步骤S1.3.2、用尺度比例L对目标轮廓极坐标信息进行归一化操作,从极径r最大时对应的轮廓点开始,以对应模板轮廓与目标轮廓各自尺寸的步长逆时针遍历取点,得到各组对应点,然后计算各组对应点极径差值的和值,将该和值除以步数作为平均误差;
步骤S1.3.3、若平均误差小于设置阈值,则认为模板轮廓与目标轮廓匹配,进入步骤S1.3.4继续执行,否则,则认为不匹配,将当前的模板轮廓替换为不同的模板轮廓,进入步骤S1.3.1继续执行;
步骤S1.3.4、将目标轮廓的轮廓中心点与模板轮廓的轮廓中心点之间的差值作为偏移量,利用各自对应极径区域极值点的极角差值平均值作为旋转角度,利用各自对应极径区域极值点的比值平均值作为尺度变化量。
3.根据权利要求1所述的运用于带钢加工生产线的视觉处理方法,其特征在于,所述步骤S2包括如下步骤:
步骤S2.1、利用高斯内核卷积带钢原始图像获得高斯金字塔;
步骤S2.2、利用Sobel算子卷积高斯金字塔各层图像获得每层的梯度幅值图,将梯度方向保存成灰度图,并计算灰度图中每一个像素点的归一化方向向量;
步骤S2.3、基于形状梯度信息的模板与第二目标图像的匹配,包括步骤:
步骤S2.3.1、在对应于最小分辨率的高斯金字塔顶层图层,将基于形状梯度信息的模板子图在第二目标图像中不断平移直至遍历第二目标图像,根据如下公式计算匹配分数S1:
其中,m、n分别为第二目标图像的长、宽,Th,k为模板子图中像素点的归一化方向向量,Dh,k为第二目标图像中对应像素点的归一化方向向量;下标h,k表示像素点在当前模板子图内的相对位置;
将第二目标图像以角度ψ为步长,以第二目标图像的几何中心为旋转中心不断旋转第二目标图像,利用二次线性插值生成新的图像矩阵,每次旋转过后都计算一次当前S1取最大值时模板子图在第二目标图像中的平移量(Sx,Sy),旋转第二目标图像直至360度,返回S1最大时对应的旋转角度ψM,记录当前平移量(SxM,SyM);Sx、SxM表示沿图像宽度方向的平移量,Sy、SyM表示沿图像高度方向的平移量;
步骤S2.3.2、对于下一层高斯金字塔图层,以小于角度ψ的步长在ψM邻域内旋转第二目标图像,获得新的图像矩阵,每次旋转结束后,在(SxM,SyM)为中心的模板子图邻域平移模板子图,利用自适应阈值分割当前层梯度幅值图像提取边缘点,具体根据如下公式计算匹配分数S2:
其中,(x,y)为边缘点坐标,T为模板子图中的归一化方向向量,D为第二目标图像中对应点的归一化方向向量,N表示第二目标图像中模板子图大小的区域中像素点梯度幅值大于设定阈值的点的数目
计算并返回S2最大时对应的旋转角度和平移量,更新ψM和(SxM,SyM),进入高斯金字塔下一层;
其中,不断重复执行步骤S2.3.2直至高斯金字塔底层,返回最终旋转角度和平移量完成匹配。
4.根据权利要求1所述的运用于带钢加工生产线的视觉处理方法,其特征在于,所述步骤S3包括如下步骤:
步骤S3.1、获取目标图像中图形的轮廓;
步骤S3.2、针对目标图像为残缺图像缺失信息而只能获得部分信息的情况,将模板轮廓与目标残缺轮廓的类曲率特征保存在一维数组中,寻找目标残缺轮廓的类曲率数组在模板轮廓类曲率数组中最相似的段落,完成残缺图像匹配;
步骤S3.3、计算模板轮廓中心点,建立模板轮廓上像素点与由模板轮廓所有轮廓点的坐标加权平均获得的中心点之间的对应关系,以该对应关系计算出目标残缺轮廓虚拟中心点,完成残缺图像定位。
5.根据权利要求2所述的运用于带钢加工生产线的视觉处理方法,其特征在于,所述步骤S1还包括如下步骤:
步骤S1.2.3、处理轮廓中心点在轮廓外部情况,具体地:
步骤S1.2.3通过以下步骤完成:
-在没有对应极径的极角范围,r取值为0;
-在没有对应极径的极角范围之外的其他角度范围,如果存在多个r值,则进行区别保存。
6.根据权利要求3所述的运用于带钢加工生产线的视觉处理方法,其特征在于,步骤S2.1中所述的高斯内核为:
7.根据权利要求3所述的运用于带钢加工生产线的视觉处理方法,其特征在于,步骤S2.2中所述归一化方向向量通过以下步骤获得:
a.利用Sobel算子对灰度图进行滤波,分别获得灰度图的高度方向Y的梯度Dy、宽度方向X的梯度Dx;
b.根据如下公式建立归一化方向向量:
EdgeDir=atan-1(Dy/Dx)
e=(COS(EdgeDir),SIN(EdgeDir))
其中,EdgeAmp为像素点的梯度方向,EdgeDir为像素点的梯度幅值,e为像素点的单位方向向量。
8.根据权利要求4所述的运用于带钢加工生产线的视觉处理方法,其特征在于,步骤S3.2中轮廓的类曲率特征通过以下步骤完成:
a.目标残缺轮廓视为连续轮廓,预处理目标残缺轮廓,截除首尾部分轮廓点;
b.从目标残缺轮廓上某一点开始逆时针遍历轮廓,取3个点,根据如下公式计算类曲率K作为类曲率特征:
其中,(xi-1,yi-1)、(xi,yi)、(xi+1,yi+1)分别为目标残缺轮廓的首端点、中间点、末端点的坐标,(xi,o,yi,o)为首端点与中间点的中垂线相交于中间点与末端点的中垂线的交点坐标。
9.根据权利要求1所述的运用于带钢加工生产线的视觉处理方法,其特征在于,还包括如下步骤:
分离检测步骤,记为步骤S4,具体为:在带钢切割后采集第三目标图像,进行分离检测。
10.根据权利要求1所述的运用于带钢加工生产线的视觉处理方法,其特征在于,所述步骤SS4包括如下步骤:
步骤S4.1、利用自适应阈值分割分离图像钢板和切割后产生的空洞部分,将空洞部分作为感兴趣区域;
步骤S4.2、根据如下公式计算匹配指标M的值,判断与设定的阈值是否一致,若一致则判定切割完成,否则判定为未完成;
其中,M为匹配指标,L为感兴趣区域周长,A为感兴趣区域面积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510684833.9A CN105335973B (zh) | 2015-10-20 | 2015-10-20 | 运用于带钢加工生产线的视觉处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510684833.9A CN105335973B (zh) | 2015-10-20 | 2015-10-20 | 运用于带钢加工生产线的视觉处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105335973A true CN105335973A (zh) | 2016-02-17 |
CN105335973B CN105335973B (zh) | 2018-08-03 |
Family
ID=55286478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510684833.9A Active CN105335973B (zh) | 2015-10-20 | 2015-10-20 | 运用于带钢加工生产线的视觉处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105335973B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106097312A (zh) * | 2016-06-01 | 2016-11-09 | 同济大学 | 一种基于机器视觉的手套撕破和油污检测方法 |
CN106504231A (zh) * | 2016-10-11 | 2017-03-15 | 广州视源电子科技股份有限公司 | 元件缺陷检测方法和系统 |
CN107220964A (zh) * | 2017-05-03 | 2017-09-29 | 长安大学 | 一种线状特征提取用于地质危岩体稳定性评估方法 |
CN107452028A (zh) * | 2017-07-28 | 2017-12-08 | 浙江华睿科技有限公司 | 一种确定目标图像位置信息的方法及装置 |
CN107671008A (zh) * | 2017-11-13 | 2018-02-09 | 中国科学院合肥物质科学研究院 | 一种基于视觉的零部件流水线自动分拣装箱设备 |
CN108492327A (zh) * | 2018-02-07 | 2018-09-04 | 广州视源电子科技股份有限公司 | Aoi图像匹配方法、系统、可读存储介质及智能设备 |
CN109205327A (zh) * | 2018-06-11 | 2019-01-15 | 杭州灵西机器人智能科技有限公司 | 一种基于多尺度轮廓直角检测的拆垛方法 |
CN112308072A (zh) * | 2020-11-06 | 2021-02-02 | 中冶赛迪重庆信息技术有限公司 | 一种废钢料场散落料识别方法、系统、电子设备及介质 |
CN113177959A (zh) * | 2021-05-21 | 2021-07-27 | 广州普华灵动机器人技术有限公司 | 一种快速运动过程中qr码实时提取算法 |
CN113378886A (zh) * | 2021-05-14 | 2021-09-10 | 珞石(山东)智能科技有限公司 | 一种自动训练形状匹配模型的方法 |
CN113421265A (zh) * | 2021-08-25 | 2021-09-21 | 苏州维嘉科技股份有限公司 | 线路板切割路径的生成方法和生成装置、分板机 |
CN114202467A (zh) * | 2020-09-16 | 2022-03-18 | 中国船舶重工集团公司第七六0研究所 | 一种水中低信噪比图像中标准月牙形残缺气泡修补方法 |
CN115170669A (zh) * | 2022-09-05 | 2022-10-11 | 合肥安迅精密技术有限公司 | 基于边缘特征点集配准的识别定位方法及系统、存储介质 |
CN117422714A (zh) * | 2023-12-18 | 2024-01-19 | 大陆汽车电子(济南)有限公司 | 装配检测方法、设备和存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441769A (zh) * | 2008-12-11 | 2009-05-27 | 上海交通大学 | 单目摄像机实时视觉定位方法 |
CN101651784A (zh) * | 2009-09-24 | 2010-02-17 | 上海交通大学 | 云台摄像机全景的视频跟踪系统 |
CN103218834A (zh) * | 2013-04-25 | 2013-07-24 | 重庆大学 | 一种基于点扩展函数的工业ct图像重建中心定位方法 |
CN103544709A (zh) * | 2013-11-07 | 2014-01-29 | 东北大学 | 基于显著活动轮廓模型的硅钢板表面微小缺陷检测方法 |
CN103612860A (zh) * | 2013-11-23 | 2014-03-05 | 冶金自动化研究设计院 | 基于机器视觉的宽厚板成品库入库定位与位置识别系统 |
CN104390591A (zh) * | 2014-11-27 | 2015-03-04 | 上海江南长兴造船有限责任公司 | 大型曲面板测量中圆形标记物的精确定位方法 |
-
2015
- 2015-10-20 CN CN201510684833.9A patent/CN105335973B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441769A (zh) * | 2008-12-11 | 2009-05-27 | 上海交通大学 | 单目摄像机实时视觉定位方法 |
CN101651784A (zh) * | 2009-09-24 | 2010-02-17 | 上海交通大学 | 云台摄像机全景的视频跟踪系统 |
CN103218834A (zh) * | 2013-04-25 | 2013-07-24 | 重庆大学 | 一种基于点扩展函数的工业ct图像重建中心定位方法 |
CN103544709A (zh) * | 2013-11-07 | 2014-01-29 | 东北大学 | 基于显著活动轮廓模型的硅钢板表面微小缺陷检测方法 |
CN103612860A (zh) * | 2013-11-23 | 2014-03-05 | 冶金自动化研究设计院 | 基于机器视觉的宽厚板成品库入库定位与位置识别系统 |
CN104390591A (zh) * | 2014-11-27 | 2015-03-04 | 上海江南长兴造船有限责任公司 | 大型曲面板测量中圆形标记物的精确定位方法 |
Non-Patent Citations (5)
Title |
---|
CONG JIA-HUI 等: "Real-time surface defects inspection of steel strip based on difference image", 《PHOTOELECTRONIC DETECTION AND IMAGING》 * |
王英: "BeckHoff 控制器在堆垛控制系统中的应用", 《科技情报开发与经济》 * |
王路: "板带钢表面质量智能检测系统研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
程万胜 等: "一种带钢表面缺陷实时检测系统的研制及应用", 《检测与仪表》 * |
赵薇: "带钢表面缺陷检测及分割技术研究", 《机械设计与制造》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106097312A (zh) * | 2016-06-01 | 2016-11-09 | 同济大学 | 一种基于机器视觉的手套撕破和油污检测方法 |
CN106097312B (zh) * | 2016-06-01 | 2019-10-01 | 同济大学 | 一种基于机器视觉的手套撕破和油污检测方法 |
CN106504231A (zh) * | 2016-10-11 | 2017-03-15 | 广州视源电子科技股份有限公司 | 元件缺陷检测方法和系统 |
CN107220964A (zh) * | 2017-05-03 | 2017-09-29 | 长安大学 | 一种线状特征提取用于地质危岩体稳定性评估方法 |
CN107452028A (zh) * | 2017-07-28 | 2017-12-08 | 浙江华睿科技有限公司 | 一种确定目标图像位置信息的方法及装置 |
CN107671008A (zh) * | 2017-11-13 | 2018-02-09 | 中国科学院合肥物质科学研究院 | 一种基于视觉的零部件流水线自动分拣装箱设备 |
CN108492327A (zh) * | 2018-02-07 | 2018-09-04 | 广州视源电子科技股份有限公司 | Aoi图像匹配方法、系统、可读存储介质及智能设备 |
CN108492327B (zh) * | 2018-02-07 | 2021-01-26 | 广州视源电子科技股份有限公司 | Aoi图像匹配方法、系统、可读存储介质及智能设备 |
CN109205327B (zh) * | 2018-06-11 | 2020-08-28 | 杭州灵西机器人智能科技有限公司 | 一种基于多尺度轮廓直角检测的拆垛方法 |
CN109205327A (zh) * | 2018-06-11 | 2019-01-15 | 杭州灵西机器人智能科技有限公司 | 一种基于多尺度轮廓直角检测的拆垛方法 |
CN114202467A (zh) * | 2020-09-16 | 2022-03-18 | 中国船舶重工集团公司第七六0研究所 | 一种水中低信噪比图像中标准月牙形残缺气泡修补方法 |
CN112308072A (zh) * | 2020-11-06 | 2021-02-02 | 中冶赛迪重庆信息技术有限公司 | 一种废钢料场散落料识别方法、系统、电子设备及介质 |
CN113378886A (zh) * | 2021-05-14 | 2021-09-10 | 珞石(山东)智能科技有限公司 | 一种自动训练形状匹配模型的方法 |
CN113177959A (zh) * | 2021-05-21 | 2021-07-27 | 广州普华灵动机器人技术有限公司 | 一种快速运动过程中qr码实时提取算法 |
CN113421265A (zh) * | 2021-08-25 | 2021-09-21 | 苏州维嘉科技股份有限公司 | 线路板切割路径的生成方法和生成装置、分板机 |
CN115170669A (zh) * | 2022-09-05 | 2022-10-11 | 合肥安迅精密技术有限公司 | 基于边缘特征点集配准的识别定位方法及系统、存储介质 |
CN117422714A (zh) * | 2023-12-18 | 2024-01-19 | 大陆汽车电子(济南)有限公司 | 装配检测方法、设备和存储介质 |
CN117422714B (zh) * | 2023-12-18 | 2024-03-29 | 大陆汽车电子(济南)有限公司 | 装配检测方法、设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105335973B (zh) | 2018-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105335973A (zh) | 运用于带钢加工生产线的视觉处理方法 | |
US20200364849A1 (en) | Method and device for automatically drawing structural cracks and precisely measuring widths thereof | |
CN102865859B (zh) | 一种基于surf特征的航空序列图像位置估计方法 | |
CN109100741A (zh) | 一种基于3d激光雷达及图像数据的目标检测方法 | |
CN103727930B (zh) | 一种基于边缘匹配的激光测距仪与相机相对位姿标定方法 | |
CN108007388A (zh) | 一种基于机器视觉的转盘角度高精度在线测量方法 | |
CN103020632B (zh) | 一种室内环境中移动机器人定位标志点的快速识别方法 | |
CN103824080B (zh) | 动态稀疏环境下机器人slam物体状态检测方法 | |
CN107392929B (zh) | 一种基于人眼视觉模型的智能化目标检测及尺寸测量方法 | |
CN103136525B (zh) | 一种利用广义Hough变换的异型扩展目标高精度定位方法 | |
US8666170B2 (en) | Computer system and method of matching for images and graphs | |
CN105096317A (zh) | 一种复杂背景中的高性能相机全自动标定方法 | |
CN105865344A (zh) | 一种基于机器视觉的工件尺寸测量方法和装置 | |
CN101650784B (zh) | 一种利用结构上下文特征进行图像匹配的方法 | |
CN113379712B (zh) | 一种基于计算机视觉的钢桥螺栓病害检测方法及系统 | |
CN108710909B (zh) | 一种可变形旋转不变装箱物体清点方法 | |
CN102903109B (zh) | 一种光学影像和sar影像一体化分割配准方法 | |
CN112308916A (zh) | 一种基于图像靶标的目标位姿识别方法 | |
CN104700395A (zh) | 一种构造物外观裂缝检测方法及系统 | |
Yuan et al. | Combining maps and street level images for building height and facade estimation | |
CN105488541A (zh) | 增强现实系统中基于机器学习的自然特征点识别方法 | |
CN116844147A (zh) | 一种基于深度学习的指针式仪表识别和异常警报方法 | |
Ma et al. | An intelligent object detection and measurement system based on trinocular vision | |
CN110044358B (zh) | 基于现场场线特征的移动机器人定位方法 | |
CN113705564B (zh) | 一种指针式仪表识别读数方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |