CN107167169B - 基于机器视觉系统的指针式仪表读数识别测量方法 - Google Patents

基于机器视觉系统的指针式仪表读数识别测量方法 Download PDF

Info

Publication number
CN107167169B
CN107167169B CN201710532391.5A CN201710532391A CN107167169B CN 107167169 B CN107167169 B CN 107167169B CN 201710532391 A CN201710532391 A CN 201710532391A CN 107167169 B CN107167169 B CN 107167169B
Authority
CN
China
Prior art keywords
pointer
instrument
image
dial plate
camera
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
CN201710532391.5A
Other languages
English (en)
Other versions
CN107167169A (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.)
Jilin University
Original Assignee
Jilin University
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 Jilin University filed Critical Jilin University
Priority to CN201710532391.5A priority Critical patent/CN107167169B/zh
Publication of CN107167169A publication Critical patent/CN107167169A/zh
Application granted granted Critical
Publication of CN107167169B publication Critical patent/CN107167169B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D9/00Recording measured values
    • G01D9/42Recording indications of measuring instruments by photographic means, e.g. of counters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D13/00Component parts of indicators for measuring arrangements not specially adapted for a specific variable

Abstract

本发明涉及一种基于机器视觉系统的指针式仪表读数识别测量方法,属于仪表测量技术领域。建立测量多个指针式仪表的视觉成像系统模型;使用CCD摄像机平面标定法,标定摄像机内部参数、畸变系数和指针式仪表表盘的外参;利用指针式仪表表盘的外参,旋转指针式仪表图像到平行于指针式仪表表盘的位置;平移各个指针式仪表图像到光学中心附近,再利用图像处理技术提取仪表图像上的刻度、指针及其中心的位置;运用角度法通过计算获得指针式仪表的读数。优点在于:采用视觉测量技术,通过图像旋转消除透视投影畸变,提高测量精度,适合用一台摄像机同时测量多个仪表的读数。实现了高精度智能检测,提高了仪表读数效率,提高了摄像机标定效率。

Description

基于机器视觉系统的指针式仪表读数识别测量方法
技术领域
本发明涉及仪表测量技术领域,特别涉及一种基于机器视觉系统的指针式仪表读数识别测量方法。
背景技术
通常仪表可分为指针式仪表和数字式仪表两大类,其中数字式仪表显示的测量值为数字, 因此具有读数简单、直观、准确的特点;而指针式仪表中包含指针和与该指针对应的刻度盘, 内部结构简单、维护方便、制造成本低、无需运行电源、抗干扰能力强,有着数字式仪表不可替代的特点,例如家用水表、汽车的速度指示表、指针式欧姆表等。因此,从科学实验和工业生产对指针式仪表的需求看,指针式仪表具有庞大的市场需求。
目前,大部分的指针式仪表并没有通信接口,指针式仪表的读数不能够自动获取并且完成数据传输,大部分读数仍然采用肉眼识别的读数方式。由于人为因素的影响或视觉偏差, 容易造成人为错误或不可避免的误差。因此,需要一种指针式仪表自动读数技术来取代人工读数方式。若一旦同时检测多个指针式仪表的读数,这对指针式仪表读数的精确度提出了更高的要求。
近些年,提出了一些基于数字图像处理的指针式仪表自动读数算法,大多数是在理想环境下进行的,保证摄像机与指针式仪表表盘水平放置。目前,机器视觉测量方法大部分是利用光学标定板进行摄像机标定获取摄像机的内部参数和外部参数。首先,通过图像采集设备获取指针式仪表的数字图像,利用连通域标记法提取出仪表表盘的轮廓和刻度线区域。其次,基于质心法提取出起始刻度线中心点和终止刻度线中心点的亚像素坐标。然后,采用霍夫变换检测出指针两端的亚像素坐标和指针式仪表表盘旋转中心的亚像素坐标。最后,运用角度法测量指针式仪表的读数。但在实际应用中,由于光照不均匀或者指针式仪表无法与摄像机保持平行等因素,导致视觉误差对读数的精度有严重影响,直接应用上述算法无法得到准确的测量结果。亟待改进。
发明内容
本发明的目的在于提供一种基于机器视觉系统的指针式仪表读数识别测量方法,解决了现有图像识别检测方法存在的上述问题。保证了较低的硬件成本、较高的检测效率和较好的测量精度。是一种更有效的多个指针式仪表读数的同时检测方法,是采用机器视觉成像的方法,考虑到CCD摄像机与指针式仪表之间的位置关系,综合使用摄像机标定和图像识别技术实现测量。这种方法利用CCD摄像机检测指针式仪表表盘的刻度线,计算指针式仪表表盘在该位置下的外部参数,取代传统的光学标定板。通过旋转指针式仪表图像到平行于指针式仪表表盘的位置,消除透视投影畸变,最后同时计算多个指针式仪表的读数。测量过程中,不要求被测指针式仪表与测量系统有特定的位置关系,放松了对测量系统的摆放要求,更有利于实际应用。本发明在已知摄像机内部参数和镜头的畸变系数的前提下,可以直接利用指针式仪表表盘的刻度线信息进行摄像机标定获取摄像机的外部参数,提高摄像机标定效率和准确性。在确定摄像机的外部参数之后,使仪表图像旋转到平行于指针式仪表表盘的位置,并通过双线性插值法使得旋转之后图像的像素坐标是整数,就可以通过对其他仪表的平移,直接读取同一平面上多个指针式仪表的读数。
本发明的上述目的通过以下技术方案实现:
基于机器视觉系统的指针式仪表读数识别测量方法,基于机器视觉成像系统,建立多个指针式仪表的测量模型,结合数字图像处理技术和摄像机标定技术实现多个指针式仪表读数的测量,其包括以下步骤:
1)标定CCD摄像机的内部参数和镜头的畸变系数;
2)提取指针式仪表图像上的刻度线中心,计算指针式仪表表盘在该刻度线中心位置下的外部参数;
3)利用指针式仪表表盘的外部参数,旋转仪表图像到平行于指针式仪表表盘的位置;
4)平移各个指针式仪表图像到光学中心附近,计算指针式仪表的读数。
上述步骤1)中所述的标定CCD摄像机的内部参数和镜头的畸变系数,是基于摄像机平面标定算法,使用高精度标定板在不同姿态下图像的角点坐标,标定出摄像机内部参数和镜头的畸变系数,具体过程包括下列步骤:
1.1)利用CCD摄像机采集16幅标定板在不同姿态下的图像;
1.2)使用Matlab软件中的Bouguet工具箱检测出图像中角点的亚像素坐标;
1.3)基于摄像机平面标定算法,利用检测得到的角点的像素坐标和世界坐标求解摄像机内参、畸变系数和外参的初值;
1.4)利用所有不同姿态的标定板图像中提取到的角点的像素坐标和世界坐标,对摄像机内部参数矩阵、畸变系数、标定板外部参数进行优化求解。
上述步骤2)中所述的提取指针式仪表图像上的刻度线中心,计算指针式仪表表盘在该刻度线中心位置下的外部参数,具体过程包括以下步骤:
2.1)保持CCD摄像机位置不变,用CCD摄像机采集多个指针式仪表图像;
2.2)利用连通域标记法提取出其中一个仪表表盘轮廓和刻度线区域;
2.3)根据仪表表盘的量程,利用质心法对提取出来的刻度线连通区域的中心进行定位,依次标记起始刻度线到终止刻度线中心点的亚像素坐标,以及长、短刻度线各自的个数;
2.4)根据指针式仪表本身刻度线的分布规律,分别测量仪表表盘圆心到长短刻度线中心点的实际距离,依据仪表表盘刻度线的个数,获得仪表表盘刻度线中心点的世界坐标;
2.5)基于张正友的平面标定法,利用检测得到的仪表表盘刻度线中心点的像素坐标和世界坐标求解仪表表盘的外部参数。
上述步骤3)中所述的利用指针式仪表表盘的外部参数,旋转指针式仪表图像到平行于指针式仪表表盘的位置,从而消除透视投影畸变,具体过程包括以下步骤:
3.1)使用通过步骤2.3)得到的刻度线中心点的亚像素坐标,利用步骤2.5)获得的指针式仪表表盘的外部参数,计算旋转后仪表表盘上的刻度线中心点在世界坐标系下的坐标;
3.2)通过步骤3.1)计算得到的旋转后仪表盘上刻度线中心点对应的目标点,利用双线性插值法使得旋转之后仪表表盘图像的像素坐标是整数,得到旋转后平行于指针式仪表表盘的仪表表盘图像;
3.3)提取出旋转后的仪表表盘图像中每一个表盘的轮廓和刻度线区域;
3.4)提取出旋转后的仪表表盘图像中每一个表盘的刻度线中心点的亚像素坐标;
3.5)采用霍夫变换方法检测出旋转后的仪表表盘图像中每一个表盘的指针直线,并确定每一个仪表表盘指针针尖的亚像素坐标和表盘旋转中心亚像素坐标。
上述步骤4)中所述的平移各个指针式仪表图像到光学中心附近,计算指针式仪表的读数,具体过程包括以下步骤:
4.1)利用步骤3.5)获得的每一个表盘旋转中心亚像素坐标,计算每个指针式仪表到距离光学中心最近的指针式仪表的平移距离;
4.2)运用角度法测量距离光学中心最近的指针式仪表的读数;
4.3)根据步骤4.1)得到的平移距离,平移其他指针式仪表表盘到离光学中心最近的指针式仪表的位置;
4.4)运用步骤4.2)测量同一平面上的其他指针式仪表的读数。
本发明的有益效果在于:引入机器视觉测量方法,基于张正友的平面标定法,使用高精度标定板在不同姿态下图像的角点坐标,标定出摄像机内部参数和镜头的畸变系数。然后利用指针式仪表表盘进行摄像机外部参数的提取,获取指针式仪表表盘与摄像机的空间位置关系之后,旋转仪表图像到平行于指针式仪表表盘的位置,消除透视投影畸变。可以任意角度检测指针式仪表读数,实现高精度智能检测;可以同时检测同一平面上的多个指针式仪表,提高仪表读数效率;可以直接利用指针式仪表表盘的刻度线信息进行摄像机标定,提高摄像机标定效率。
基于经典的平面标定法对仪表表盘进行标定,旋转指针式仪表图像到平行于指针式仪表表盘的位置,有效地修正了多个指针式仪表表盘中被测点的透视投影畸变;利用仪表表盘刻度线的分布规律,提取刻度线中心点作为标定板的角点,简化了标定过程,提高了测量精度;在对仪表表盘视觉成像系统完成标定之后,测量时,只需要保证CCD摄像机的位置不变,放松了测量系统与被测仪表表盘之间摆放位置的要求,更利于实际应用;本发明采取了视觉测量技术,因此能实现同时检测多个指针式仪表,同时采用了精确的图像识别技术,能保证较高的测量精度。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实例及其说明用于解释本发明,并不构成对本发明的不当限定。
图1为本发明的指针式仪表读数机器视觉检测方法的CCD摄像机标定内部参数硬件示意图;
图2为本发明的指针式仪表读数机器视觉检测方法的指针式仪表表盘的摄像机外部参数标定的硬件示意图;
图3为本发明的指针式仪表的读数示意图;
图4为本发明的平移多个指针式仪表的示意图;
图5为本发明的多个指针式仪表测量系统模型。
图中:1、计算机;2、CCD摄像机;3、镜头;4、指针式仪表;5、标定板;6、背光源。
具体实施方式
下面结合附图进一步说明本发明的详细内容及其具体实施方式。
参见图1至图5所示,本发明的基于机器视觉系统的指针式仪表读数识别测量方法,基于机器视觉成像系统,建立多个指针式仪表的测量模型,结合数字图像处理技术和摄像机标定技术实现多个指针式仪表读数的测量。能够在线同时测量多个指针式仪表的读数。首先,建立测量多个指针式仪表的视觉成像系统模型;其次,使用CCD摄像机平面标定法,标定摄像机内部参数、畸变系数和指针式仪表表盘的外参;再次,利用指针式仪表表盘的外参,旋转指针式仪表图像到平行于指针式仪表表盘的位置,消除透视投影畸变;然后,平移各个指针式仪表图像到光学中心附近,再利用图像处理技术提取仪表图像上的刻度、指针及其中心的位置;最后,运用角度法通过计算获得指针式仪表的读数。本发明采用视觉测量技术,通过图像旋转消除透视投影畸变,提高测量精度,适合用一台摄像机同时测量多个仪表的读数。具体步骤如下:
1、标定CCD摄像机的内部参数及镜头的畸变系数,是基于摄像机平面标定算法,使用高精度标定板在不同姿态下图像的角点坐标,标定出摄像机内部参数和镜头的畸变系数,具体过程包括下列步骤:
1.1)利用CCD摄像机采集16幅标定板在不同姿态下的图像。采集图像时,如图1所示,标定板5置于背光源6和CCD摄像机2之间。
1.2)使用Matlab软件中的Bouguet工具箱检测出图像中角点的亚像素坐标;
1.3)基于摄像机平面标定算法,利用检测得到的角点的像素坐标和世界坐标求解摄像机内参、畸变系数和外参的初值;
在标定摄像机内部参数、外部参数以及畸变系数的过程中,涉及的四个坐标系如下:PCS 像素坐标系O0uv;ICS图像坐标系Oxy;CCS摄像机坐标系OCXCYCZC;WCS世界坐标系OWXWYWZW。摄像机的成像过程即是空间物点在这四个坐标系中的一系列变换过程,本实施例采用的标定模型如下:
其中H=A[r1 r2 t]为3×3矩阵,称为单应性矩阵,表征了 WCS和PCS之间的线性几何关系。在[r1 r2 t]中的6个参数θ,ψ,Tx,Ty,Tz被称为摄像机外部参数,它们决定了CCS与WCS之间的相对位置关系,如公式(2):
其中t=(Tx,Ty,Tz)T是一个3×1的平移向量,是一个3×3的单位正交旋转矩阵,这个旋转矩阵由三个旋转角度确定,分别为绕CCS的ZC轴旋转角度θ,绕YC轴旋转角度ψ, XC轴旋转角度由三个欧拉角表示的旋转矩阵R的形式如下:
对于针孔摄像机模型,从CCS到ICS的变换是透视投影,则图像投影点在ICS中的坐标可表示为:
由于实际摄像机成像时,光学镜头几何畸变的存在,需要考虑到镜头的径向畸变系数k1, k2和切向畸变系数p1,p2
其中,(xu,yu)和(xd,yd)分别为理想图像坐标和实际图像坐标。由数字图像的产生原理可得,从ICS到PCS的变换是数字化图像过程,
矩阵A中的参数α,β,γ,u0,v0称为摄像机线性成像模型的内部参数,表达式如下:
1.4)利用所有不同姿态的标定板图像中提取到的角点的像素坐标和世界坐标,对摄像机内部参数矩阵、畸变系数、标定板外部参数进行优化求解。
2、提取指针式仪表图像上的刻度线中心,计算指针式仪表表盘在该刻度线中心的外部参数,具体过程包括以下步骤:
2.1)保持CCD摄像机位置不变,将多个指针式仪表固定在同一墙面上,使得多个指针式仪表是在同一平面上。用CCD摄像机采集多个指针式仪表图像时,如图2所示,CCD摄像机2与镜头3的组合以及同一平面上的多个指针式仪表4呈非平行放置。
2.2)将采集到的多个指针式仪表图像送到计算机1系统里,采用3×3中值滤波器进行降噪处理,对降噪后的灰度图像利用阈值分割法转化为二值图像,有效地划分开仪表表盘与背景。针对二值图像进行连通域标记,计算仪表表盘中各个连通区域的面积。根据连通区域面积的大小保留仪表表盘轮廓、刻度线和指针,去除仪表表盘背景。结合数学形态学的相关理论,通过膨胀腐蚀运算和骨架细化提取出一个仪表表盘轮廓以及刻度线区域。
2.3)根据仪表表盘的量程,利用质心法对提取出来的刻度线连通区域的中心进行定位,从起始刻度到终止刻度依次标记刻度线中心点的亚像素坐标(xp,yp),以及长、短刻度线各自的个数np。若目标像素点分布在一个矩形窗口内,窗口左上角像素坐标为(p1,q1),右下角像素坐标为(p2,q2),则它的中心位置计算公式如下:
式(8)中,(xn,yn)为目标像素点质心坐标估计值;I(i,j)为(i,j)处像素点的灰度值。
2.4)根据指针式仪表本身刻度线的分布规律,用游标卡尺分别测量仪表表盘圆心到长、短刻度线中心点的实际距离r1和r2,依据仪表表盘刻度线的个数np,分别获得仪表的长、短刻度线中心点的世界坐标(Xn (1),Yn (1),Zn (1))和(Xn (2),Yn (2),Zn (2))。
由于指针式仪表表盘中的刻度线都是均匀分布的,式(9)中η=(θ×π)/(180×np)是指相邻刻度线之间的弧度,θ是指起始刻度线逆时针到终止刻度线的角度,并将长、短刻度线的世界坐标重新排序整合成世界坐标(Xn,Yn,Zn)。
2.5)基于张正友的平面标定法,利用检测得到的仪表表盘刻度线中心点的像素坐标和世界坐标求解仪表表盘的外部参数,建立该位置仪表表盘平面上被测点的世界坐标与摄像机坐标系之间的变换关系(即两坐标系之间的旋转矩阵和平移矩阵)。在已知摄像机内部参数和畸变系数的条件下,根据奇异值分解(SVD)可以解得单应性矩阵,即仪表表盘的摄像机外部参数。
3、利用指针式仪表表盘的外部参数,旋转指针式仪表图像到平行于指针式仪表表盘的位置,从而消除透视投影畸变,具体过程包括以下步骤:
3.1)使用通过步骤2.3)计算得到的刻度线中心点的亚像素坐标(xp,yp),利用步骤2.5) 获得的指针式仪表表盘的外部参数,计算旋转后仪表表盘上的刻度线中心点在世界坐标系下的坐标。将刻度线中心点坐标(xp,yp)带入公式(6),并通过公式(5)、(4)和(2)便可以确定该目标点在世界坐标系下的坐标。
3.2)通过步骤3.1)计算得到的旋转后仪表盘上刻度线中心点对应的目标点,利用双线性插值法使得旋转之后仪表表盘图像的像素坐标是整数,得到旋转后平行于指针式仪表表盘的仪表表盘图像。
3.3)提取出旋转后的仪表表盘图像中每一个表盘的轮廓和刻度线区域;
3.4)提取出旋转后的仪表表盘图像中每一个表盘的起始刻度线中心点和终止刻度线中心点的亚像素坐标;
3.5)采用霍夫变换方法检测出旋转后的仪表表盘图像中每一个表盘的指针直线,并确定每一个仪表表盘指针针尖的亚像素坐标和表盘旋转中心亚像素坐标。
4、平移各个指针式仪表图像到光学中心附近,计算指针式仪表的读数,具体过程包括以下步骤:
4.1)利用步骤3.5)获得的每一个表盘旋转中心亚像素坐标,如图3所示,计算每个指针式仪表到距离光学中心最近的指针式仪表的平移距离d。
其中Tx=xc2′-xc1′是两个指针式仪表的水平距离;
Ty=yc2′-yc1′是两个指针式仪表的垂直距离。
4.2)运用角度法测量距离光学中心最近的指针式仪表的读数,如图4所示。
其中:α=arccos[(xc'-x1',yc'-y1'),(xc'-xn',yc'-yn')]零刻度线和满刻度线之间的夹角β=arccos[(xc'-x1',yc'-y1'),(xc'-xh',yc'-yh')]指针和零刻度线之间的夹角 M即指针式仪表的总量程。
4.3)根据步骤4.1)得到的平移距离,平移其他指针式仪表表盘到离光学中心最近的指针式仪表的位置;
4.4)运用步骤4.2)测量同一平面上的其他指针式仪表的读数。
以上所述仅为本发明的优选实例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡对本发明所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于机器视觉系统的指针式仪表读数识别测量方法,其特征在于:基于机器视觉成像系统,建立多个指针式仪表的测量模型,结合数字图像处理技术和摄像机标定技术实现多个指针式仪表读数的测量,其包括以下步骤:
1)标定CCD摄像机的内部参数和镜头的畸变系数;
2)提取指针式仪表图像上的刻度线中心,计算指针式仪表表盘在该刻度线中心位置下的外部参数;
3)利用指针式仪表表盘的外部参数,旋转仪表图像到平行于指针式仪表表盘的位置;
4)平移各个指针式仪表图像到光学中心附近,计算指针式仪表的读数;
步骤2)中所述的提取指针式仪表图像上的刻度线中心,计算指针式仪表表盘在该刻度线中心位置下的外部参数,具体过程包括以下步骤:
2.1)保持CCD摄像机位置不变,将多个指针式仪表固定在同一墙面上,使得多个指针式仪表是在同一平面上;用CCD摄像机采集多个指针式仪表图像时,CCD摄像机(2)与镜头(3)的组合以及同一平面上的多个指针式仪表(4)呈非平行放置;
2.2)将采集到的多个指针式仪表图像送到计算机(1)系统里,采用3×3中值滤波器进行降噪处理,对降噪后的灰度图像利用阈值分割法转化为二值图像,有效地划分开仪表表盘与背景;针对二值图像进行连通域标记,计算仪表表盘中各个连通区域的面积;根据连通区域面积的大小保留仪表表盘轮廓、刻度线和指针,去除仪表表盘背景;结合数学形态学的理论,通过膨胀腐蚀运算和骨架细化提取出一个仪表表盘轮廓以及刻度线区域;
2.3)根据仪表表盘的量程,利用质心法对提取出来的刻度线连通区域的中心进行定位,从起始刻度到终止刻度依次标记刻度线中心点的亚像素坐标(xp,yp),以及仪表表盘刻度线的个数np;若目标像素点分布在一个矩形窗口内,窗口左上角像素坐标为(p1,q1),右下角像素坐标为(p2,q2),则它的中心位置计算公式如下:
式(8)中,(xn,yn)为目标像素点质心坐标估计值;I(i,j)为(i,j)处像素点的灰度值;
2.4)根据指针式仪表本身刻度线的分布规律,用游标卡尺分别测量仪表表盘圆心到长、短刻度线中心点的实际距离r1和r2,依据仪表表盘刻度线的个数np,分别获得仪表的长、短刻度线中心点的世界坐标(Xn (1),Yn (1),Zn (1))和(Xn (2),Yn (2),Zn (2));
由于指针式仪表表盘中的刻度线都是均匀分布的,式(9)中η=(θ×π)/(180×np)是指相邻刻度线之间的弧度,θ是指起始刻度线逆时针到终止刻度线的角度,并将长、短刻度线的世界坐标重新排序整合成世界坐标(Xn,Yn,Zn);
2.5)基于张正友的平面标定法,利用检测得到的仪表表盘刻度线中心点的像素坐标和世界坐标求解仪表表盘的外部参数,建立仪表表盘平面上被测点的世界坐标与摄像机坐标系之间的变换关系,即两坐标系之间的旋转矩阵和平移矩阵;在已知摄像机内部参数和畸变系数的条件下,根据奇异值分解SVD解得单应性矩阵,即仪表表盘的摄像机外部参数;
步骤3)中所述的利用指针式仪表表盘的外部参数,旋转指针式仪表图像到平行于指针式仪表表盘的位置,从而消除透视投影畸变,具体过程包括以下步骤:
3.1)使用通过步骤2.3)计算得到的刻度线中心点的亚像素坐标(xp,yp),利用步骤2.5)获得的指针式仪表表盘的外部参数,计算旋转后仪表表盘上的刻度线中心点在世界坐标系下的坐标;将刻度线中心点坐标(xp,yp)带入公式(6),并通过公式(5)、(4)和(2)确定目标点在世界坐标系下的坐标;
在标定摄像机内部参数、外部参数以及畸变系数的过程中,涉及的四个坐标系如下:PCS像素坐标系O0uv;ICS图像坐标系Oxy;CCS摄像机坐标系OCXCYCZC;WCS世界坐标系OWXWYWZW
其中t=(Tx,Ty,Tz)T是一个3×1的平移向量,是一个3×3的单位正交旋转矩阵;对于针孔摄像机模型,从CCS到ICS的变换是透视投影,则图像投影点在ICS中的坐标可表示为:
由于实际摄像机成像时,光学镜头几何畸变的存在,需要考虑到镜头的径向畸变系数k1,k2和切向畸变系数p1,p2
其中,(xu,yu)和(xd,yd)分别为理想图像坐标和实际图像坐标;由数字图像的产生原理可得,从ICS到PCS的变换是数字化图像过程,
矩阵A中的参数α,β,γ,u0,v0称为摄像机线性成像模型的内部参数;
3.2)通过步骤3.1)计算得到的旋转后仪表盘上刻度线中心点对应的目标点,利用双线性插值法使得旋转之后仪表表盘图像的像素坐标是整数,得到旋转后平行于指针式仪表表盘的仪表表盘图像;
3.3)提取出旋转后的仪表表盘图像中每一个表盘的轮廓和刻度线区域;
3.4)提取出旋转后的仪表表盘图像中每一个表盘的刻度线中心点的亚像素坐标;
3.5)采用霍夫变换方法检测出旋转后的仪表表盘图像中每一个表盘的指针直线,并确定每一个仪表表盘指针针尖的亚像素坐标和表盘旋转中心亚像素坐标。
2.根据权利要求1所述的基于机器视觉系统的指针式仪表读数识别测量方法,其特征在于:步骤1)中所述的标定CCD摄像机的内部参数和镜头的畸变系数,是基于摄像机平面标定算法,使用高精度标定板在不同姿态下图像的角点坐标,标定出摄像机内部参数和镜头的畸变系数,具体过程包括下列步骤:
1.1)利用CCD摄像机采集16幅标定板在不同姿态下的图像;
1.2)使用Matlab软件中的Bouguet工具箱检测出图像中角点的亚像素坐标;
1.3)基于摄像机平面标定算法,利用检测得到的角点的像素坐标和世界坐标求解摄像机内参、畸变系数和外参的初值;
1.4)利用所有不同姿态的标定板图像中提取到的角点的像素坐标和世界坐标,对摄像机内部参数矩阵、畸变系数、标定板外部参数进行优化求解。
3.根据权利要求1所述的基于机器视觉系统的指针式仪表读数识别测量方法,其特征在于:步骤4)中所述的平移各个指针式仪表图像到光学中心附近,计算指针式仪表的读数,具体过程包括以下步骤:
4.1)利用步骤3.5)获得的每一个表盘旋转中心亚像素坐标,计算每个指针式仪表到距离光学中心最近的指针式仪表的平移距离;
4.2)运用角度法测量距离光学中心最近的指针式仪表的读数;
4.3)根据步骤4.1)得到的平移距离,平移其他指针式仪表表盘到离光学中心最近的指针式仪表的位置;
4.4)运用步骤4.2)测量同一平面上的其他指针式仪表的读数。
CN201710532391.5A 2017-07-03 2017-07-03 基于机器视觉系统的指针式仪表读数识别测量方法 Active CN107167169B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710532391.5A CN107167169B (zh) 2017-07-03 2017-07-03 基于机器视觉系统的指针式仪表读数识别测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710532391.5A CN107167169B (zh) 2017-07-03 2017-07-03 基于机器视觉系统的指针式仪表读数识别测量方法

Publications (2)

Publication Number Publication Date
CN107167169A CN107167169A (zh) 2017-09-15
CN107167169B true CN107167169B (zh) 2019-04-26

Family

ID=59828112

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710532391.5A Active CN107167169B (zh) 2017-07-03 2017-07-03 基于机器视觉系统的指针式仪表读数识别测量方法

Country Status (1)

Country Link
CN (1) CN107167169B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107729853B (zh) * 2017-10-24 2021-05-25 武汉科技大学 一种适用于变电站窄刻度指针式仪表的自动识别方法
CN108038489A (zh) * 2017-12-20 2018-05-15 杭州讯能科技有限公司 一种指针式仪表识别方法、系统及嵌入式系统
CN108764257B (zh) * 2018-05-23 2022-03-04 郑州金惠计算机系统工程有限公司 一种多视角的指针式仪表识别方法
CN108458733B (zh) * 2018-06-20 2023-09-22 安徽炬视科技有限公司 基于指针检测的指针式仪表自动读数装置系统及识别方法
WO2020019287A1 (en) * 2018-07-27 2020-01-30 Linde Aktiengesellschaft Attachable gauge reading device
CN109063717B (zh) * 2018-07-30 2021-06-15 安徽慧视金瞳科技有限公司 一种获取仪表中心点方法
CN109974623B (zh) * 2018-09-13 2021-01-05 北京伟景智能科技有限公司 基于线激光和双目视觉的三维信息获取方法和装置
CN109544628B (zh) * 2018-11-23 2021-08-31 福州大学 一种指针式仪表的准确读数识别系统及方法
CN111458337A (zh) * 2019-03-18 2020-07-28 珠海经济特区南森科技有限公司 一种钟表检测装置及其检测方法
CN110189314B (zh) * 2019-05-28 2023-04-28 长春大学 基于机器视觉的汽车仪表盘图像定位方法
CN110189375B (zh) * 2019-06-26 2022-08-23 中国科学院光电技术研究所 一种基于单目视觉测量的图像目标识别方法
CN112749598A (zh) * 2019-10-31 2021-05-04 中移物联网有限公司 一种指针式表计的信息处理方法、终端及数据平台
CN110789238A (zh) * 2019-11-29 2020-02-14 无锡星月科技有限公司 一种用于仪表刻度盘自动打印设备
CN111476787A (zh) * 2020-04-23 2020-07-31 中科开创(广州)智能科技发展有限公司 一种指针表自适应畸变的自动读数方法及装置
CN112001379A (zh) * 2020-05-19 2020-11-27 西安工程大学 基于机器视觉的汽车仪表固定视点读取仪表的修正算法
CN112285625A (zh) * 2020-09-21 2021-01-29 广东电网有限责任公司广州供电局 一种通用型电力仪表自动检定系统及其检定方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6935566B1 (en) * 1997-02-03 2005-08-30 Symbol Technologies, Inc. Portable instrument for electro-optically reading indicia and for projecting a bit-mapped image
CN101498592B (zh) * 2009-02-26 2013-08-21 北京中星微电子有限公司 指针式仪表的读数方法及装置
CN102176228B (zh) * 2011-01-28 2013-02-06 河海大学常州校区 一种用于识别多个指针式仪表表盘信息的机器视觉方法
CN103499302B (zh) * 2013-09-27 2015-11-25 吉林大学 基于结构光视觉成像系统的凸轮轴直径尺寸在线测量方法
CN103471531B (zh) * 2013-09-27 2016-01-20 吉林大学 轴类零件直线度在线非接触测量方法
CN105758439A (zh) * 2014-12-15 2016-07-13 沈阳晨讯希姆通科技有限公司 机器人检测装置及系统
CN106599897B (zh) * 2016-12-09 2019-10-18 广州供电局有限公司 基于机器视觉的指针式仪表读数识别方法及装置

Also Published As

Publication number Publication date
CN107167169A (zh) 2017-09-15

Similar Documents

Publication Publication Date Title
CN107167169B (zh) 基于机器视觉系统的指针式仪表读数识别测量方法
CN106989669B (zh) 基于虚拟立体靶标的大视场高精度视觉系统标定方法
CN102376089B (zh) 一种标靶校正方法及系统
CN112818988B (zh) 一种指针式仪表自动识别读数方法及系统
CN109558871B (zh) 一种指针式仪表读数识别方法及装置
CN102622747B (zh) 一种用于视觉测量的摄像机参数优化方法
CN104634248B (zh) 一种双目视觉下的转轴标定方法
CN109579695B (zh) 一种基于异构立体视觉的零件测量方法
CN110223355B (zh) 一种基于双重极线约束的特征标志点匹配方法
CN102800096B (zh) 一种摄像机参数的鲁棒性估计算法
CN105783773A (zh) 一种线结构光视觉传感器的数值标定方法
Yang et al. An image-based intelligent system for pointer instrument reading
CN102183214A (zh) 一种大口径非球面镜结构光检测方法
CN105278260B (zh) 一种pcb曝光图形正确性验证方法
CN109341668A (zh) 基于折射投影模型和光束追踪法的多相机测量方法
CN107941471A (zh) 一种自由曲面的检测方法
CN110455222A (zh) 一种高精度旋转角测量方法、装置及设备
CN109360794A (zh) 一种晶硅光伏太阳能电池电极二次印刷精度视觉检测方法及装置
CN106403838A (zh) 一种手持式线结构光视觉三维扫描仪的现场标定方法
WO2024011981A1 (zh) 一种圆柱形靶光斑校验方法
CN109712157A (zh) 一种基于单目视觉的重力场法加速度计校准方法
CN109211268A (zh) 一种仪表指针校准方法、装置及计算机可读存储介质
CN105423975A (zh) 一种大型工件的标定系统及方法
CN106447729A (zh) 一种基于坐标变换的二维数字图像相关补偿方法及二维光学数字图像相关引伸计
CN112488030A (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