CN105300482A - 基于图像处理的水表检定方法、装置及系统 - Google Patents

Abstract

本发明公开基于图像处理的水表检定方法、装置及系统，该水表检定方法包括：从指针二值化图中提取图像指针轮廓线，根据指针轮廓线计算所有轮廓点到指针圆圆心的第一距离，定义最大值为longest_diatance；计算每个指针轮廓点到指针圆圆心的第二距离；搜索第二距离中范围在0.8倍longest_diatance至1倍longest_diatance之间的轮廓点；求取轮廓点的平均坐标值，并依据平均坐标值确定指针针尖位置；根据指针针尖位置与指针圆圆心确定指针方向；根据指针方向对预设水表进行示数判读，得到判读结果；计算判读结果并确定预设水表的示值状态。本发明在图像识别技术的基础上，对水表图像进行处理，使得特征区域提取出来，然后进行判读，最后利用获取的多幅图像对判读准确性进行验证。

Description

基于图像处理的水表检定方法、装置及系统

技术领域

本发明涉及仪表检测领域，更具体的，涉及一种基于图像处理的水表检定方法、装置及系统。

背景技术

冷水水表检定过程是一个确定水表读数误差大小的过程，在检定时可以适当采用计算机识别图像读数代替人工方式进行检测。现阶段，指针式水表检定系统的相关研究日趋增多，众多的国内外研究者将研究的方向设定在机器视觉以及图像处理上，即采集指针式水表图像，再利用计算机处理图像矩阵识别水表指针位置，进而读取水表示值。图像处理检测技术具有极强的对象性，不同的检测对象的图像结构信息相差很大，而且识别的精度速度等要求也会不尽相同。近几年，国内外图像处理专家在指针式仪表识别方面不断探索，也取得了一定的研究成果。

关于指针式水表只是众多指针式仪表检测中的一种，而基于图像处理的指针式仪表的研究在国外起步较早，上世纪八九十年代开始，欧美许多发达国家的研究人员就已经对此逐渐探索，最早是于1989年J.P.Jones选用机械手对指针式仪表图像采集处理，进而进行指针读数判读；1994年，美国韦恩州立大学的Sablatning等人就已在自己文章中详细叙述了以圆形表盘且刻度分布均匀的水表作为识别对象，基于Hough变换识别直线即指针位置，再求取指针指向读取示数，成为仪表识别领域代表性突破；2000年，葡萄牙里斯本技术大学的CorreaAlegria等创新性的提出针对多幅表盘图像，联合减影法以及Hough变换在将圆形表盘变换成水平刻度线样式后再进行指针位置的识别；2004年，S.L.Pang等人提出基于最大灰度相减法利用含指针图像减无指针模板图像实现角度判读；2006年出现以抑制光晕为目的的基于中心环绕模型；2007年，美国科研人员利用指针对环状光栅盘的遮挡所得到的明暗光信号来读取指针示数。

国外的很多学者对于基于图像处理的指针式仪表识别算法理论研究大大推动了实现仪表自动识别的发展，尤其是对于欧美等发达国家来讲，使用的图像处理逐渐增多，技术的发展也就相对较快，经验比较丰富，应用方面日渐普及到各个行业如汽车、制药、食品、饮料和包装等。即便如此，国外的图像处理系统由于知识产权等原因价格昂贵、维修困难，在国内无法被广大用户所接受。因此，我国研制出具备自主知识产权的指针式仪表识别系统地必要性很强。

同国外相比，我国针对指针式仪表图像识别的研究较晚，研究方向主要集中在图像分割和特征识别算法上。最早进行此方面研究的是哈尔滨工业大学的李铁桥教授，他指出以表盘圆心为极点在极坐标系中根据识别出的指针的极坐标角度确定读数；接着王三武等人对于多子表盘的水表图像识别系统进行了研究，主要采用的是模版匹配；华北电力大学的李宝树指出识别指针刻度线然后根据预先设定的待测检点与实际检测点之间的误差读数进行判读的方法；丁庆生、李世平等人提出以指针点、表盘圆心和零点这三点为顶点构成三角模型进行指针读取的方法；香港理工大学W.L.Chan等人设计出一套指针仪表的自动检定系统；华南理工大学的陈彬将中心投影同Hough变换相结合确定指针位置；南京航空航天大学的吴欢欢、游有鹏等将ActiveX控件技术引入指针仪表监测系统，提高的检测系统稳定性；北京邮电大学的曲仁军等对基于嵌入式环境下仪表快速识别理论进行了探索，结合圆环提取的改进以及边沿梯度圆心查找法对指针读数进行判读；武汉大学方彦军以基于机器视觉的仪表自动判读项目为依托对自动校验装置对准控制进行了深入钻研，我国学者对实现指针式仪表示数识别一直都在不断探索。

基于以上对前人指针式仪表识别算法的研究不难发现这些指针式仪表识别系统以指针作为识别判读的关键，对指针方向识别的方法大体分为四类：步长法、圆周灰度检测、模版匹配、圆周弦长法，对指针方向进行判读的主要方法基本就是两种，即角度法和距离法。

步长法根本原理是以指针旋转轴为种子原点分别沿周围八个邻域方向中某一方向进行搜索，直到搜索到邻域灰度值是0的位置为止，再以此位置为原点重复搜索过程，在震荡中搜索与原点距离最远的位置，此点即为最终点，即指针针尖位置。然后连接原点位置与针尖位置所成的线即为指针指向，但此方法是以指针轮廓线为基础，如果获得的轮廓线有任何断裂情况就会使得震荡搜索停止，从而造成指针针尖位置提取错误，而且读错的概率很大；圆周灰度检测法的基本思想是在确定了指针长度L的基础上以轴心为圆心以L为半径做圆并确定与指针边缘的两个交点，再确定两交点间的圆弧的中点，将轴心与中点连线，此连线即为指针方向，但此方法必须在指针长度确定而且有一定宽度的情况下才较为合适的使用，如果指针相对较细则造成两交点距离过近，而有的厂商生产的指针式水表也有指针相对较细的情况，应用此种方法效果并不是很好；而对于模版匹配的方法需要提前设定好一个固定模版，然后将采集的图像与模版相比较，确定指针的方向，但在本实施例想要设计系统是手持式图像采集终端，模版并不固定，由此对于模版匹配的方法更加不适合。

针对传统方法各个方面的局限性，很多学者提出了运用圆周弦长检测法来进行水表指针的方向识别，这种方法的基本原理是将指针以旋转轴心为原点将指针分为四个象限，然后从原点出发向两个相反的方向搜索指针轮廓，记录轮廓点间距离，每搜索完一条线上的两个轮廓点就以1°为步长顺时针或逆时针继续搜索直到旋转180°完成搜索，此搜索过程中所得的最长的轮廓距离即对应指针指示方向。虽然此种方法在指针式仪表识别领域具有十分广泛的用途，但如果指针较尖锐，与旋转指针圆没有平滑过度的情况下会浪费许多资源，而且如果指针在提取过程中出现平针的现象会造成指针方向识别不准确甚至是错误。

因此，现有技术中水表检定还存在效率及准确率均不高的问题。

发明内容

本发明公开一种基于图像处理的水表检定方法、装置及系统，用于解决现有技术中水表检定还存在效率及准确率均不高的问题。

为实现上述目的，本发明提供一种基于图像处理的水表检定方法，并采用如下技术方案：

基于图像处理的水表检定方法，包括：获取预设表盘图像的指针二值化图；从所述指针二值化图中提取图像指针轮廓线，根据所述指针轮廓线计算所有轮廓点到指针圆圆心的第一距离；确定所述第一距离中的最大值，定义所述最大值为longest_diatance；计算每个指针轮廓点到所述指针圆圆心的第二距离；搜索所述第二距离中范围在0.8倍longest_diatance至1倍longest_diatance之间的轮廓点；求取所述轮廓点的平均坐标值，并依据所述平均坐标值确定指针针尖位置；根据所述指针针尖位置与所述指针圆圆心确定指针方向；根据所述指针方向对所述预设水表进行示数判读，得到判读结果；计算所述判读结果并确定所述预设水表的示值状态。

进一步地，所述获取预设表盘图像的指针二值化图包括：对所述表盘图像进行灰度化处理，得到第一处理结果；在所述第一处理结果的基础上对所述表盘图像进行表盘图像增强，得到第二处理结果；在所述第二处理结果的基础上对所述表盘图像进行滤波去噪，得到第三处理结果；在所述第三处理结果的基础上对所述表盘图像进行二值化处理，得到第四处理结果；在所述第四处理结果的基础上对所述表盘图像进行特征区域提取，得到第五处理结果；在所述第五处理结果的基础上对所述表盘图像进行指针提取，得到第六处理结果；在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图。

进一步地，所述在所述第三处理结果的基础上对所述表盘图像进行二值化处理包括：计算经过滤波去噪的所述表盘图像的灰度均值，计算公式如下：

$\stackrel{\‾}{m}=\frac{1}{height*width}{\mathrm{\Σ}}_{x=0}^{height}{\mathrm{\Σ}}_{y=0}^{width}f(x,y)$

以0为阈值将图像分为大于等于均值的第一部分和小于均值的第二部分；计算所述第一部分与所述第二部分内像素数目n1、n2，再计算所述第一部分的平均值与所述第二部分的平均值计算公式如下：

$\left\{\begin{array}{cc}\stackrel{\&OverBar;}{m_1}=\frac{{\mathrm{\&Sigma;}}_{n=0}^{n_1}f{(x,y)}_n}{n_1}& f(x,y)<\stackrel{\&OverBar;}{m}\\ \stackrel{\&OverBar;}{m_2}=\frac{{\mathrm{\&Sigma;}}_{n=0}^{n_2}f{(x,y)}_n}{n_2}& f(x,y)\&GreaterEqual;\stackrel{\&OverBar;}{m}\end{array}\right.$

根据所述第一部分的平均值与所述第二部分的平均值计算类间方差，计算公式如下：

$\mathrm{\&sigma;}=\frac{n_1}{height*width}*{(\stackrel{\&OverBar;}{m_1}-\stackrel{\&OverBar;}{m})}^2+\frac{n_2}{height*width}*{(\stackrel{\&OverBar;}{m_2}-\stackrel{\&OverBar;}{m})}^2$

将阈值范围在[0，255]处的类间方差做比较，获取类方差最大时的阈值，并根据所述类方差最大时的阈值处理所述表盘图像；其中，width为宽度、height为高度、(x，y)位置处灰度大小用f(x，y)表示，表盘为图像的灰度均值σ表示类间方差。

进一步地，所述在所述第四处理结果的基础上对所述表盘图像进行特征区域提取：在所述表盘图像二值化的基础上对所述表盘图像进行霍夫变换算法，识别表盘，再识别子表盘；依据LRCD算法对所述子表盘特征区域提取红色指针，具体方法包括：采用预设参数模型对所述表盘图像进行灰度变换，所述预设参数模型为：

g(x，y)＝R-Y＝0.701R-0.587G-0.114B

式中g(x，y)为处理后灰度大小，R、G、B分别代表图像某位置的三通道灰度值大小，Y为加权平均值法灰度化所得灰度值；利用红色信道的灰度与蓝色或绿色信道相应灰度做差，获取灰度差即：

g(x，y)＝R-G或g(x，y)＝R-B

根据所述灰度差对所述表盘图像进行所述第四处理。

进一步地，所述根据所述指针方向对所述预设水表进行示数判读包括：确定所述预设水表的表盘圆心位置以及所述指针圆圆心的位置；根据预设计算公式以及所述表盘圆心位置、所述指针圆圆心的位置、识别的所述指针方向计算所述预设水表的指针读数；所述预设计算公式为

$\tan \mathrm{\&theta;}=\frac{(q-n)*(x-m_2)-(y-n_2)*(p-m)}{(x-m_2)*(p-m)+(y-n_2)*(q-n)};$

其中，所述表盘圆心位置为(x，y)，所述指针圆圆心的位置按照从左到右的顺序依次为(m1,n1)、(m2,n2)、(m3,n3)、(m4,n4),对应的针尖的位置依次为(p1,q1)、(p2,q2)、(p3,q3)、(p4,q4)，以向量(x-m2,y-n2)为零向量，计算各个指针方向向量(p-m，q-n)与零向量所成的角度。

进一步地，在所述根据所述指针方向对所述预设水表进行示数判读之前，所述水表检定方法还包括：判断所述指针方向向量与所述零向量的数量积是否为零；在所述指针方向向量与所述零向量的数量积为零时，确定所述角度为90度或270度；在所述指针方向向量与所述零向量的数量积不为零时，则利用所述预设计算公式求取所述角度的向量夹角θ。

进一步地，所述在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图包括：基于霍夫变换算法提取所述子表盘，并对所述子表盘进行区域分割；依据所述LRCD色差变换模型对所述子表盘中的红色指针进行提取，得到所述指针二值化图。

根据本发明的另外一个方面，提供一种基于图像处理的水表检定装置，并采用如下技术方案：

该基于图像处理的水表检定装置包括：获取模块，用于获取预设表盘图像的指针二值化图；提取模块，用于从所述指针二值化图中提取图像指针轮廓线，根据所述指针轮廓线计算所有轮廓点到指针圆圆心的第一距离；第一确定模块，用于所述第一距离中的最大值，定义所述最大值为longest_diatance；第一计算模块，用于计算每个指针轮廓点到所述指针圆圆心的第二距离；搜索模块，用于搜索所述第二距离中范围在0.8倍longest_diatance至1倍longest_diatance之间的轮廓点；求取模块，用于求取所述轮廓点的平均坐标值，并依据所述平均坐标值确定指针针尖位置；第二确定模块，用于根据所述指针针尖位置与所述指针圆圆心确定指针方向；判读模块，用于根据所述指针方向对所述预设水表进行示数判读，得到判读结果；第二计算模块，用于计算所述判读结果并确定所述预设水表的示值状态。

进一步所述获取模块包括：灰度化处理模块，用于对所述表盘图像进行灰度化处理，得到第一处理结果；增强模块，用于在所述第一处理结果的基础上对所述表盘图像进行表盘图像增强，得到第二处理结果；去噪模块，用于在所述第二处理结果的基础上对所述表盘图像进行滤波去噪，得到第三处理结果；二值化处理模块，用于在所述第三处理结果的基础上对所述表盘图像进行二值化处理，得到第四处理结果；特征区域提取模块，用于在所述第四处理结果的基础上对所述表盘图像进行特征区域提取，得到第五处理结果；指针提取模块，用于在所述第五处理结果的基础上对所述表盘图像进行指针提取，得到第六处理结果；指针二值化处理模块，用于在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图。

根据本发明的又一个方面，提供一种基于图像处理的水表检定系统，并采用如下技术方案：

基于图像处理的水表检定系统包括上述的水表检定装置。

本发明提高的水表读数识别方法，具有如下技术效果：

第一、对指针式水表图像预处理过程中的各种方法进行了分析，图像预处理尤其对相对清晰度不高干扰较大的表盘图像显得更为重要。

第二、实现了对指针式水表图像特征区域子表盘的提取以及指针提取。在指针子表盘提取过程中并非针对识别子表盘圆内信息提取而是采用扩充提取以保证指针针尖部分的完整性。对于指针的提取过程并非直接对灰度图像进行处理，而是充分利用图像彩色信息，利用LRCD彩色差分原理提取指针。

第三、实现了指针圆完整准确检定。利用提取指针信息准确识别指针转轴即指针圆圆心位置。

附图说明

附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1表示本发明实施例所述的基于图像处理的水表检定方法的流程图；

图2表示本发明实施例所述的基于水表图像进行处理后所得到的指针二值化图；

图3表示本发明实施例所述的确定指针轮廓示意图；

图4表示本发明实施例所述的确定指针方向的原理图；

图5表示本发明实施例所述的确定指针方向的又一原理图；

图6表示本发明实施例所述的确定指针方向示意图；

图7表示本发明实施例所述的确定表盘圆心的示意图；

图8表示本发明实施例所述的确定表盘圆心的又一示意图；

图9表示本发明实施例所述的确定表盘圆心的原理图；

图10表示本发明实施例所述的水表图像的灰度图；

图11表示图10所述的灰度图的直方图；

图12表示本发明实施所述的圆形空间参数示意图；

图13表示本发明实施例所述的特征区域提取中的图像表盘识别示意图；

图14表示本发明实施例所述的特征区域提取中的图像子表盘识别示意图；

图15表示本发明实施例所述的原始识别子表盘提取特征区域图；

图16表示本发明实施例所述的圆的半径扩大后的子表盘特征区域提取图；

图17表示本发明实施例所述的水表图像的红色通道灰度图；

图18表示本发明实施例所述的水表图像的绿色通道灰度图；

图19表示本发明实施例所述的水表图像的蓝色通道灰度图；

图20表示本发明实施例所述的水表图像的红色指针提取灰度图；

图21表示本发明实施例所述的R-G红色指针灰度图；

图22表示本发明实施例所述的R-B红色指针灰度图；

图23表示本发明实施例所述的指针直方图；

图24表示本发明实施例所述的指针二值化效果图；

图25表示本发明实施例所述的基于图像处理的水表检定装置的结构示意图。

具体实施方式

在通过实施例详细阐述本发明技术方案及技术效果之前，先介绍一下水表检测的相关步骤及方法。

对于指针式水表检测项目主要包括三类：外观和功能检测、密封性排查、示值误差计算。尤以最后一项示值误差作为检定关键指标。

示值误差检定是水表处于水温低于30℃环境中水流速度分别在最小流量点Q₁、分界流量点Q₂、常用流量点Q₃三个流量点处的相对误差，是检测结果对仪表准确度等级判别的根本体现。流量是水表实际流过的水流体积与所用时间之商，当然在实际实验时，不可能将流量点设置的如此精确，所以一般的实际流量分别控制在Q₁～1.1Q₁、Q₂～1.1Q₂、0.9Q₃～Q₃范围内即可。指针式水表的相对误差E是用百分数表示的一个重要检测指标，计算如下：

$E=\frac{V_i-V_a}{V_a}\&times;100\%$

式中V_i表示水表指示体积差，V_a表示标准设备实际体积差。当进行检定时，水表的最大允许误差在检定低区(Q₁≤Q＜Q₂)是±5％，检定高区是±2％。

检定步骤

1、首先观察被检水表的外表是否存在问题，尤其是察看流量点、水表流向等各项标注是否完整，若不完整直接返厂或淘汰。

2、将外观符合标准的水表按照水流流向与水表流向指示箭头一致的方向正确装到检定试验台，打开出水阀门再缓慢打开进水阀门，对整组水表以及管道进行排气，根据玻璃转子流量计显示的水表排气情况确定是否终止排气，如果玻璃转子流量计内没有气泡冒出时并不关闭进水阀门和出水阀门，而是继续流水过程，目的是稳定流量和检查水表是否发生水泄漏情况。水流流过一定合适时间(本人在检定过程中的合适时间采取的是大于等于一分钟)关闭进水阀门，开始检定过程的初始状态读数，首先读取的是标准超声波流量计的数据V₀，然后读取各个指针式水表表盘数据V₀’。

3、再次打开进水阀门开始水表检定。检定的终止条件为流经水表的水流体积相当于水表检定流量下排放一分钟的体积。关闭试验台水表上游进水阀门同时关闭下游出水阀门。

4、记录标准超声波流量计检定终止末态时数据V₁，以及各个指针式水表的检定末态数据V₁’。

5、根据误差计算公式，计算被检水表示值误差，分析误差统计参数[25]，若处于最大允许误差范围内即可证明此水表为合格水表，符合标注精度等级，此时才能打印合格证明以及检定证书。

以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。

图1表示本发明实施例一所述的基于图像处理的水表检定方法的流程图。

参见图1所示，基于图像处理的水表检定方法包括：

S101:获取预设表盘图像的指针二值化图；

S103:从所述指针二值化图中提取图像指针轮廓线，根据所述指针轮廓线计算所有轮廓点到指针圆圆心的第一距离；

S105:确定所述第一距离中的最大值，定义所述最大值为longest_diatance；

S107:计算每个指针轮廓点到所述指针圆圆心的第二距离；

S109:搜索所述第二距离中范围在0.8倍longest_diatance至1倍longest_diatance之间的轮廓点；

S111:求取所述轮廓点的平均坐标值，并依据所述平均坐标值确定指针针尖位置；

S113:根据所述指针针尖位置与所述指针圆圆心确定指针方向；

S115:根据所述指针方向对所述预设水表进行示数判读，得到判读结果；

S117:计算所述判读结果并确定所述预设水表的示值状态。

在步骤S101中，预设表盘图像，目前实现数字图像采集设备的构成主要包括工业相机、图像采集卡、传输设备、上位机。目前按照成像原理将工业相机分为模拟和数字两类。本实施例所用设备内嵌的数字工业相机是基于逐行曝光的CMOS图像传感器，内置心跳功能可自动诊断工作状态从异常中恢复，支持多种采集方式即连续采集和触发采集。其自动配备通用性1394接口，符合DCAM(discriminatingcontentaddressablememory)规范标准，在向主机传输图像信息过程中占用资源较少，同时此设备还具有较好的环境适应性，十分适合图像采集传输等应用场合。

在图像摄取过程中，尽量以表盘直径为基准，使得表盘尽量接近图像宽度长度。将图像通过图像采集设备成功采集后需要将其传入计算机，然后开始进行处理分析。

在获取步骤S101中表盘图像后，对表盘图像进行处理，以期获取表盘图像的指针二值化图，参见图2。

在步骤S103，从所述指针二值化图中提取图像指针轮廓线，具体参见图3所示，根据所述指针轮廓线计算所有轮廓点到指针圆圆心的第一距离。

在得到指针二值化图后就需要对指针的方向开始进行识别，对于指针式水表指针方向的识别方法多种多样。针对现有技术在水表检测方面存在的缺陷，本实施例指针二值化图为基础，识别出指针提取图像中的指针圆。

结合指针本质特征，由于无论是指针尖端与指针圆过渡平滑的情况还是指针指示部分与指针圆剧烈过渡，通过观察不难发现指针针尖所对应的位置均是距离指针圆圆心最大的，由此通过确定指针圆圆心位置以及指针针尖位置来识别指针方向,具体可参见图4至图5所示。

算法的具体描述如下：

计算所有轮廓点pt到指针圆圆心的第一距离distance。

$dis\tan ce=\sqrt{{(pt\&RightArrow;x-cvRound(p\&lsqb;0\&rsqb;\left)\right)}^2+{(pt\&RightArrow;y-cvRound(p\&lsqb;1\&rsqb;\left)\right)}^2}---(5-1)$

其中，pt→x代表所述轮廓点的横坐标，pt→y代表所述轮廓点的纵坐标；

S105:确定所述第一距离中的最大值，定义所述最大值为longest_diatance；

从上述的第一距离distance中获取最大距离longest_diatance。

S107:计算每个指针轮廓点到所述指针圆圆心的第二距离；

计算每个指针轮廓点到指针圆圆心位置的距离，称为第二距离。

S109:搜索所述第二距离中范围在0.8倍longest_diatance至1倍longest_diatance之间的轮廓点；

在第二距离中搜索距离范围在[0.8*longest_diatance，longest_diatance]的轮廓点，统计上述轮廓点。

S111:求取所述轮廓点的平均坐标值，并依据所述平均坐标值确定指针针尖位置；

求取上述轮廓点的平均坐标值[average_x，average_y]。

S113:根据所述指针针尖位置与所述指针圆圆心确定指针方向。

指针圆圆心同轮廓点平均坐标值的连线是指针指向，参见图6所示。

S115:根据所述指针方向对所述预设水表进行示数判读，得到判读结果；

在获得指针的方向后需要对指针读数进行判别，指针读数的判别方法可以采用距离法或角度法。

角度法是利用识别出来的指针方向与表盘指针指向零刻度时所成的角度与整个子表盘量程之比计算指针读数。这种方法最重要的是对指针角度的判断[65]，可以在未知刻度线的情况下读取示值，而且水表子表盘的刻度线与刻度值呈线性关系，刻度分布均匀，因此本实施例采用角度法。

观察水表子表盘发现各个子表盘的圆心共圆，自然指针圆圆心也在此圆上，指针圆圆心成为此圆圆周上的点。所以在检测出指针圆后连接各圆心组成六条圆周弦，参见图7。由于每天圆周弦的垂直平分线必过表盘圆心，根据两两垂直平分线的交点计算交点均值从而确定出表盘圆心位置。其连接各圆心的圆周弦效果图以及表盘圆心位置确定，参见图8。

一般根据确定的表盘圆心位置、指针圆圆心位置以及识别的指针方向就可以利用角度法进行指针读数的确定，本实施例是基于向量求取读数。

设在指针式水表表盘中检测到的表盘圆心位置为(x，y)，指针圆的位置按照从左到右的顺序依次为(m1,n1)、(m2,n2)、(m3,n3)、(m4,n4),对应的针尖的位置依次为(p1,q1)、(p2,q2)、(p3,q3)、(p4,q4)，以向量(x-m2,y-n2)为零向量，计算各个指针方向向量(p-m，q-n)与零向量所成的角度，原理图如图9所示。

基于公式： $\tan \mathrm{\&theta;}=\frac{(q-n)*(x-m_2)-(y-n_2)*(p-m)}{(x-m_2)*(p-m)+(y-n_2)*(q-n)}$

在确定所述预设水表的表盘圆心位置以及所述指针圆圆心的位置；根据预设计算公式以及所述表盘圆心位置、所述指针圆圆心的位置、识别的所述指针方向计算所述预设水表的指针读数。

在步骤S117中，计算所述判读结果并确定所述预设水表的示值状态。

根据所求角度θ，以及整个子表盘均匀分布十个刻度即每36度刻度加一，所以角度与表盘读数关系如下表所示：

虽然按照对应关系可以求得指针对数，但当指针处于36度倍数的临界角等敏感刻度点时往往很小的误差就会导致读数变换差别很大，尤其是随着以立方米为单位的子表盘×0.0001、×0.001、×0.01、×0.1所代表的计数单位增大所带来的判别误差逐渐增大甚至不可估量。

为了消除计算所带来的误差，可以根据各子表盘之间十进制关系利用小一位的计数表盘对进制子表盘进行读数纠正，如最小子表盘×0.0001的指向为9，×0.001指示为比2稍大一点则判其为2而非3。

本实施例在基于指针二值化水表图像的基础上，实现了指针圆完整准确检定，利用提取指针信息准确识别指针转轴即指针圆圆心位置。然后结合改进的角度判读方法，实现指针的准确判读。

优选地，在所述根据所述指针方向对所述预设水表进行示数判读之前，所述水表检定方法还包括：判断所述指针方向向量与所述零向量的数量积是否为零；在所述指针方向向量与所述零向量的数量积为零时，确定所述角度为90度或270度；在所述指针方向向量与所述零向量的数量积不为零时，则利用所述预设计算公式求取所述角度的向量夹角θ。

由于在计算过程中存在指针方向向量与零向量垂直的情况，因此以上计算公式针对非垂直情况下进行计算，所以计算之前首先需要判断指针方向向量与零向量的数量积是否为零，即(x-m2)*(p-m)+(y-n2)*(q-n)是否为零，若为零则判断(x-m2)*(p-m)的值，若不小于零则为90°，否则为270°。若(x-m2)*(p-m)+(y-n2)*(q-n)不为零则利用公式求取向量夹角θ，然后判断夹角所属类别，由于求的θ的范围为[0，180]范围内，所以若(x-m2)*(p-m)大于零则所成角度为θ，否则为θ＝180+θ。

作为优选的实施例，本实施例给出指针二值化图的获取方式，具体包括：

对所述表盘图像进行灰度化处理，得到第一处理结果；在所述第一处理结果的基础上对所述表盘图像进行表盘图像增强，得到第二处理结果；在所述第二处理结果的基础上对所述表盘图像进行滤波去噪，得到第三处理结果；在所述第三处理结果的基础上对所述表盘图像进行二值化处理，得到第四处理结果；在所述第四处理结果的基础上对所述表盘图像进行特征区域提取，得到第五处理结果；在所述第五处理结果的基础上对所述表盘图像进行指针提取，得到第六处理结果；在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图。

针对对所述表盘图像进行灰度化处理，得到第一处理结果，具体的实施方式可以为：

将所述表盘图像的各像素点的R、G、B三个分量分别乘以一个加权系数即权值，数学表达式如下：

g(x，y)＝(W_R*R+W_G*G+W_B*B)/3

采用如下加权参数模型：

g(x，y)＝0.299R+0.587G+0.114B

对所述各像素点进行加权平均，获取所述表盘图像的灰度图，所得到的灰度图即为第一处理结果。

在所述第一处理结果的基础上对所述表盘图像进行表盘图像增强，得到第二处理结果，即图像增强，具体如下：

图像增强是针对采集图像灰度分布集中问题通过一系列图像变换手段扩大灰度直方图灰度级覆盖领域，增强目标区域与背景信息之间对比度差异，使得图像转换成一种更加适于计算机分析和处理的形式。

为了进行指针式水表图像增强算法的研究，需要首先对采集的图像灰度化后的灰度直方图进行研究分析，具体参见图10至图11，通过水表灰度图像的统计直方图的灰度级分布是否充分利用了[0,255]灰度空间。图11即为图10的直方统计图。

从水表图像的统计直方图不难看出，摄取的表盘存在灰度级分布过于集中，分布范围狭窄，对比度较低的缺点。为了增强边缘轮廓等细节清晰度需要对其进行增强。

灰度范围变换是表盘图像增强的主要方法，用数学表达式的方式表示如下：

g(x，y)＝T[f(x，y)]

f(x,y)代表原水表图像中(x,y)位置处像素值，g(x,y)代表增强后表盘相应位置像素值，T为水表图像变换函数。根据灰度变换函数的不同，下面对几种灰度变换的常用灰度变换方法分析研究。

线性比例灰度拉伸变换原理为对原图的某段灰度范围基于线性函数进行变换的过程，灰度分布过于集中的水表图像可以直接映射到大于此范围的灰度级空间，对于灰度级范围较大但仅某狭窄范围内同一灰度值的数量统计较多，分布概率较大的图像进行灰度变换时可以将有限个关键的灰度级大于某一阈值的拓展到另一合适灰度范围，提高所要识别的目标的亮度，达到增强有效像素的目的。

所述在所述第三处理结果的基础上对所述表盘图像进行二值化处理包括：计算经过滤波去噪的所述表盘图像的灰度均值，计算公式如下：

$\stackrel{\&OverBar;}{m}=\frac{1}{height*width}{\mathrm{\&Sigma;}}_{x=0}^{height}{\mathrm{\&Sigma;}}_{y=0}^{width}f(x,y)$

以0为阈值将图像分为大于等于均值的第一部分和小于均值的第二部分；计算所述第一部分与所述第二部分内像素数目n1、n2，再计算所述第一部分的平均值与所述第二部分的平均值计算公式如下：

$\left\{\begin{array}{cc}\stackrel{\&OverBar;}{m_1}=\frac{{\mathrm{\&Sigma;}}_{n=0}^{n_1}f{(x,y)}_n}{n_1}& f(x,y)<\stackrel{\&OverBar;}{m}\\ \stackrel{\&OverBar;}{m_2}=\frac{{\mathrm{\&Sigma;}}_{n=0}^{n_2}f{(x,y)}_n}{n_2}& f(x,y)\&GreaterEqual;\stackrel{\&OverBar;}{m}\end{array}\right.$

根据所述第一部分的平均值与所述第二部分的平均值计算类间方差，计算公式如下：

$\mathrm{\&sigma;}=\frac{n_1}{height*width}*{(\stackrel{\&OverBar;}{m_1}-\stackrel{\&OverBar;}{m})}^2+\frac{n_2}{height*width}*{(\stackrel{\&OverBar;}{m_2}-\stackrel{\&OverBar;}{m})}^2$

将阈值范围在[0，255]处的类间方差做比较，获取类方差最大时的阈值，并根据所述类方差最大时的阈值处理所述表盘图像；其中，width为宽度、height为高度、(x，y)位置处灰度大小用f(x，y)表示，表盘为图像的灰度均值σ表示类间方差。

所述在所述第四处理结果的基础上对所述表盘图像进行特征区域提取：在所述表盘图像二值化的基础上对所述表盘图像进行霍夫变换算法，识别表盘，再识别子表盘；依据LRCD算法对所述子表盘特征区域提取红色指针，具体方法包括：采用预设参数模型对所述表盘图像进行灰度变换，所述预设参数模型为：

g(x，y)＝R-Y＝0.701R-0.587G-0.114B

式中g(x，y)为处理后灰度大小，R、G、B分别代表图像某位置的三通道灰度值大小，Y为加权平均值法灰度化所得灰度值；利用红色信道的灰度与蓝色或绿色信道相应灰度做差，获取灰度差即：

g(x，y)＝R-G或g(x，y)＝R-B

根据所述灰度差对所述表盘图像进行所述第四处理。

在所述第四处理结果的基础上对所述表盘图像进行特征区域提取，得到第五处理结果；

在完成了指针式水表图像一系列图像预处理及二值化后，紧接着就是需要对特征区域进行提取。指针式水表表盘图像的表盘、子表盘均为圆形，而Hough变换在线、圆、椭圆的识别上具有很大优势，所以本实施例结合Hough变换算法识别表盘，再识别子表盘，进而基于子表盘特征区域引入LRCD算法提取红色指针。

Hough变换是一种基于点-线的对偶性使用了投票表决机制将图像复杂边缘特征信息映射为所要检测的图像形状的参数估计技术，这种映射是根据所要检测的形状的参数定义的函数，这种变换将所有符合参数方程的像素点映射为统计为参数空间中参数峰值特征。Hough变化在直线、圆、椭圆等检测方面具有很大优势，能够不受几何形状中间断点影响，鲁棒性较强，被广泛应用于图像处理领域。

Hough变换为对图像进行两种不同的坐标系数据变换，将X-y二维坐标系中数据变换至以参数为基底的多维坐标系中，维数的数量就是由参数方程中自由参数的数量。其总体变换过程是首先将检测到的连续的图形边缘变换到参数空间，然后将参数空间按照提前设定的步长单位离散化，形成离散参数矩阵，将所有二值图像中的有效像素点带入参数表达式，查看满足参数表达式的参数矩阵元素位置并进行累加，位置不同得到的累加值不同，但峰值点所在的位置是原函数的参数具体值的概率最大，所以去峰值点为参数值。Hough变换最早应用于直线检测，在指针式仪表识别时只要是针对线状指针如电流表、高精度压力表等，因此在这方面的应用研究起初集中通过减少变换区域提高Hough变换直线检定速度，但本实施例是用Hough变换检测圆形，所以下面以圆的检测过程为例进行详细介绍。

在图像处理中基于Hough变换对圆进行检测可以分为两种情况，一种是已知圆的半径情况下对图像中的圆形的检测，另一种是未知圆的半径情况下对图像中的圆形的检测，但一般图像中圆形的半径是未知量，因此在对圆形进行检测过程中需要以圆心位置和半径范围为基底建立三维坐标系进行圆的识别。以(a,b)为圆心，以r为半径的圆的方程为：

(x-a)²+(x-b)²＝r²

方程中的a,b,r即为映射空间的维数基底，将图像空间中的任意一点映射到此空间为一个圆锥曲面，原图像中的圆心坐标即为所有在此圆上的像素映射的圆锥曲面的交点，如图12所示。

以上这种Hough变换检测适合于圆上像素点位置未知，需要以整幅图像的所有提取轮廓为变换对象进行变换，但如果已知图像中某几个点共圆，则采用快速Hough更为合理简便。快速Hough变换的本质思想即为共圆的点的连线即为此圆的弦，而弦的垂直平分线必过圆心，已知A、B、C、D四点共圆，通过连接A、B和C、D得到两条圆周弦，再分别做圆周弦的垂直平分线，交点就是共圆圆心位置，也就是E点。

特征提取的实质是删除原始数据中的无用信息，选择最为有效最能反映事物本质的特征，并数学的方式进行描述的过程。在进行识别过程中，针对对象不同，提取的图像的特征不具有广泛性，所以需要具体问题具体分析，具体情况选择合适的方法进行特征提取。

对于图像的特征区域提取主要是根据采集的指针式水表图像特征提取水表子表盘区域部分。为了较为准确的提取，首先需要识别水表表盘，然后识别水表子表盘，提取子表盘特征区域。

在对图像进行采集的时，将表盘的轮廓尽量铺满所采集到的图像，然后对图像进行Hough变换识别水表表盘，识别时将圆的半径R范围设置为[width/3,width/2](其中width为图像的宽度)，阈值设置小于100，其实在识别表盘过程中有了圆的半径的限制使得阈值的设置只要不是太大均能检测出来。然后根据表盘的结构信息，每个表盘的半径大约为子表盘的4-6倍，根据识别表盘半径的识别子表盘，利用Hough算法对子表盘识别过程中将子表盘圆的半径r范围设置为[R/4,R/6]，对图像表盘以及子表盘的识别分别如图13至图14所示。

在图13至图14中，黑色加粗部分代表红色。

由于实施例对红色指针的提取是基于LRCD方法，需要用到颜色信息，所以为了下一步的指针提取需要以原始彩色图像为基准，然后对识别的子表盘进行区域保留，在对子表盘识别过程中已将子表盘的圆心位置、半径等信息分别存储于提前构建的数组当中，因此在提取彩色特征区域时只需构建一幅背景图像，再将子表盘区域内像素保留即可。

在图14中，从图14子表盘识别效果图可以看出子表盘的识别并不是十分精准，如果在彩色特征区域提取时直接仅仅将圆内像素保留会很容易导致红色指针针尖提取不全，出现平针、斜针的情况。但针尖的位置对于后续的指针方向的判读具有绝对性作用，由此为了完全提取红色指针需要将提取区域的半径设置为比检测半径稍大，因而设置为[(R+10)/4,(R+10)/6]，即将半径扩大十个像素左右就不会影响指针的提取。如图15、图16分别为原始识别子表盘提取特征区域图和圆的半径扩大后的子表盘特征区域提取图，其中黑色指针代表红色指针，因为附图不能有颜色。

在所述第五处理结果的基础上对所述表盘图像进行指针提取，得到第六处理结果；

第五处理结果，即为上述对原始识别子表盘提取特征区域和圆的半径扩大后的子表盘特征区域提取，在上述结果上，对表盘图像进行指针提取，具体方法如下：

LRCD是一种色差变换方法，目的是为了突出图像的红色成分。LRCD是一种用于成熟草莓果实采摘的机器视觉系统颜色模型。我国中国农业大学在2012年的草莓大会上展示的“采摘童一号”草莓采摘机器人对草莓的识别就是利用LRCD识别方法。首先求得其中每个像素色差，在灰度图像上显示以色差值为灰度的图像，选取恰当阈值进行二值化处理，完成图像分割，再提取图像几何特征。LRCD灰度变换的参数模型为：

g(x，y)＝R-Y＝0.701R-0.587G-0.114B

式中g(x，y)为处理后灰度大小，R、G、B分别代表图像某位置的三通道灰度值大小，Y为加权平均值法灰度化所得灰度值。但是在采用此变换之前需要将图像的三个RGB通道分开，分裂成三幅图像再进行运算。具体参见图17，图18、图19、图20，从图20的处理效果可以看出已经能够将红色指针提取出来，但LRCD模型中有乘法运算，对应每个像素点都需要进行三次乘法再做加减，耗费运算时间，因此为了提高运算处理速度，利用红色信道的灰度与蓝色或绿色信道相应灰度做差，如下式所示：

g(x，y)＝R-G或g(x，y)＝R-B

上面两中减法方式所得的结果基本相似，处理后的效果图分别如图21，图22所示。

将LRCD红色指针提取灰度图与改进后的效果图相比较，图像对红色指针的提取效果是十分接近，所以采用改进后的能省去乘法运算时间。

本实施例对于在红色指针提取完整的基础上，为了较为准确的获取读数，首先需要将红色指针提取灰度图转化为二值图像，然后识别指针方向。

在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图。

对于二值化算法，将红色指针灰度图二值化，首先需要对图像的直方图进行分析，普通图像的直方图是对图像的所有灰度值进行统计分析，但由于红色指针灰度图的背景颜色设置为灰度值为0的黑色，且像素数量相对于红色指针像素数量过多，若仍然对所有像素值的灰度进行直方图进行统计就会造成图像的直方图仅在0位置处有很大的值，其他位置基本忽略的现象，无法实现对指针灰度直方图的统计，所以本实施例对红色指针提取的灰度图对应的直方图从1开始统计分析，所得的直方图如图23所示。

从图23中可以看出其为较为理想灰度分布图像，图像灰度统计直方图中存在两个波峰夹一个波谷，这样将波谷所在位置灰度选做阈值就可以取得较好的二值化效果，其结果如图24所示。

在上述实施例中，首先介绍了水表滤波后图像的二值化，然后基于霍夫变换提取子表盘，实现子表盘区域分割，再利用LRCD色差变换模型实现红色指针的提取，最后将提取红色指针图像进行二值化。

本发明的上述实施例，实现了对指针式水表图像特征区域子表盘的提取以及指针提取。在指针子表盘提取过程中并非针对识别子表盘圆内信息提取而是采用扩充提取以保证指针针尖部分的完整性。对于指针的提取过程并非直接对灰度图像进行处理，而是充分利用图像彩色信息，利用LRCD彩色差分原理提取指针。实现了指针圆完整准确检定。利用提取指针信息准确识别指针转轴即指针圆圆心位置。

图25表示本发明实施例所述的基于图像处理的水表检定装置的结构示意图。

本发明提供的基于图像处理的水表检定装置，该基于图像处理的水表检定装置包括：获取模块2501，用于获取预设表盘图像的指针二值化图；提取模块2502，用于从所述指针二值化图中提取图像指针轮廓线，根据所述指针轮廓线计算所有轮廓点到指针圆圆心的第一距离；第一确定模块2503，用于所述第一距离中的最大值，定义所述最大值为longest_diatance；第一计算模块2504，用于计算每个指针轮廓点到所述指针圆圆心的第二距离；搜索模块2505，用于搜索所述第二距离中范围在0.8倍longest_diatance至1倍longest_diatance之间的轮廓点；求取模块2506，用于求取所述轮廓点的平均坐标值，并依据所述平均坐标值确定指针针尖位置；第二确定模块2507，用于根据所述指针针尖位置与所述指针圆圆心确定指针方向；判读模块2508，用于根据所述指针方向对所述预设水表进行示数判读，得到判读结果；第二计算模块2509，用于计算所述判读结果并确定所述预设水表的示值状态。

优选地，所述获取模块2501包括：灰度化处理模块(图中未示)，用于对所述表盘图像进行灰度化处理，得到第一处理结果；增强模块(图中未示)，用于在所述第一处理结果的基础上对所述表盘图像进行表盘图像增强，得到第二处理结果；去噪模块(图中未示)，用于在所述第二处理结果的基础上对所述表盘图像进行滤波去噪，得到第三处理结果；二值化处理模块(图中未示)，用于在所述第三处理结果的基础上对所述表盘图像进行二值化处理，得到第四处理结果；特征区域提取模块(图中未示)，用于在所述第四处理结果的基础上对所述表盘图像进行特征区域提取，得到第五处理结果；指针提取模块(图中未示)，用于在所述第五处理结果的基础上对所述表盘图像进行指针提取，得到第六处理结果；指针二值化处理模块(图中未示)，用于在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图。

本发明提供的基于图像处理的水表检定系统包括上述的水表检定装置。

本发明提高的水表读数识别方法，具有如下技术效果：

第一、对指针式水表图像预处理过程中的各种方法进行了分析，图像预处理尤其对相对清晰度不高干扰较大的表盘图像显得更为重要。

第二、实现了对指针式水表图像特征区域子表盘的提取以及指针提取。在指针子表盘提取过程中并非针对识别子表盘圆内信息提取而是采用扩充提取以保证指针针尖部分的完整性。对于指针的提取过程并非直接对灰度图像进行处理，而是充分利用图像彩色信息，利用LRCD彩色差分原理提取指针。

第三、实现了指针圆完整准确检定。利用提取指针信息准确识别指针转轴即指针圆圆心位置。

以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种基于图像处理的水表检定方法，其特征在于，包括：

获取预设表盘图像的指针二值化图；

从所述指针二值化图中提取图像指针轮廓线，根据所述指针轮廓线计算所有轮廓点到指针圆圆心的第一距离；

确定所述第一距离中的最大值，定义所述最大值为longest_diatance；

计算每个指针轮廓点到所述指针圆圆心的第二距离；

搜索所述第二距离中范围在0.8倍longest_diatance至1倍longest_diatance之间的轮廓点；

求取所述轮廓点的平均坐标值，并依据所述平均坐标值确定指针针尖位置；

根据所述指针针尖位置与所述指针圆圆心确定指针方向；

根据所述指针方向对预设水表进行示数判读，得到判读结果；

计算所述判读结果并确定所述预设水表的示值状态。

2.如权利要求1所述的水表检定方法，其特征在于，所述获取预设表盘图像的指针二值化图包括：

对所述表盘图像进行灰度化处理，得到第一处理结果；

在所述第一处理结果的基础上对所述表盘图像进行表盘图像增强，得到第二处理结果；

在所述第二处理结果的基础上对所述表盘图像进行滤波去噪，得到第三处理结果；

在所述第三处理结果的基础上对所述表盘图像进行二值化处理，得到第四处理结果；

在所述第四处理结果的基础上对所述表盘图像进行特征区域提取，得到第五处理结果；

在所述第五处理结果的基础上对所述表盘图像进行指针提取，得到第六处理结果；

在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图。

3.如权利要求2所述的水表检定方法，其特征在于，所述在所述第三处理结果的基础上对所述表盘图像进行二值化处理包括：

计算经过滤波去噪的所述表盘图像的灰度均值，计算公式如下：

$\stackrel{\&OverBar;}{m}=\frac{1}{height*width}{\mathrm{\&Sigma;}}_{x=0}^{height}{\mathrm{\&Sigma;}}_{y=0}^{width}f(x,y)$

以0为阈值将图像分为大于等于均值的第一部分和小于均值的第二部分；

计算所述第一部分与所述第二部分内像素数目n1、n2，再计算所述第一部分的平均值与所述第二部分的平均值计算公式如下：

$\left\{\begin{array}{cc}\stackrel{\&OverBar;}{m_1}=\frac{{\mathrm{\&Sigma;}}_{n=0}^{n_1}f{(x,y)}_n}{n_1}& f(x,y)<\stackrel{\&OverBar;}{m}\\ \stackrel{\&OverBar;}{m_2}=\frac{{\mathrm{\&Sigma;}}_{n=0}^{n_2}f{(x,y)}_n}{n_2}& f(x,y)\&GreaterEqual;\stackrel{\&OverBar;}{m}\end{array}\right.$

根据所述第一部分的平均值与所述第二部分的平均值计算类间方差，计算公式如下：

$\mathrm{\&sigma;}=\frac{n_1}{height*width}*{(\stackrel{\&OverBar;}{m_1}-\stackrel{\&OverBar;}{m})}^2+\frac{n_2}{height*width}*{(\stackrel{\&OverBar;}{m_2}-\stackrel{\&OverBar;}{m})}^2$

将阈值范围在[0，255]处的类间方差做比较，获取类方差最大时的阈值，并根据所述类方差最大时的阈值处理所述表盘图像；

其中，width为宽度、height为高度、(x，y)位置处灰度大小用f(x，y)表示，表盘为图像的灰度均值σ表示类间方差。

4.如权利要求2所述的水表检定方法，其特征在于，所述在所述第四处理结果的基础上对所述表盘图像进行特征区域提取：

在所述表盘图像二值化的基础上对所述表盘图像进行霍夫变换算法，识别表盘，再识别子表盘；

依据LRCD算法对所述子表盘特征区域提取红色指针，具体方法包括：

采用预设参数模型对所述表盘图像进行灰度变换，所述预设参数模型为：

g(x，y)＝R-Y＝0.701R-0.587G-0.114B

式中g(x，y)为处理后灰度大小，R、G、B分别代表图像某位置的三通道灰度值大小，Y为加权平均值法灰度化所得灰度值；

利用红色信道的灰度与蓝色或绿色信道相应灰度做差，获取灰度差

g(x，y)＝R-G或g(x，y)＝R-B

根据所述灰度差对所述表盘图像进行所述第四处理。

5.如权利要求1所述的水表检定方法，其特征在于，所述根据所述指针方向对所述预设水表进行示数判读包括：

确定所述预设水表的表盘圆心位置以及所述指针圆圆心的位置；

根据预设计算公式以及所述表盘圆心位置、所述指针圆圆心的位置、识别的所述指针方向计算所述预设水表的指针读数；

所述预设计算公式为： $tan\mathrm{\&theta;}=\frac{(q-n)*(x-m_2)-(y-n_2)*(p-m)}{(x-m_2)*(p-m)+(y-n_2)*(q-n)};$

其中，所述表盘圆心位置为(x，y)，所述指针圆圆心的位置按照从左到右的顺序依次为(m1，n1)、(m2，n2)、(m3，n3)、(m4，n4)，对应的针尖的位置依次为(p1，q1)、(p2，q2)、(p3，q3)、(p4，q4)，以向量(x-m2，y-n2)为零向量，计算各个指针方向向量(p-m，q-n)与零向量所成的角度。

6.如权利要求1所述的水表检定方法，其特征在于，在所述根据所述指针方向对所述预设水表进行示数判读之前，所述水表检定方法还包括：

判断所述指针方向向量与所述零向量的数量积是否为零；

在所述指针方向向量与所述零向量的数量积为零时，确定所述角度为90度或270度；

在所述指针方向向量与所述零向量的数量积不为零时，则利用所述预设计算公式求取所述角度的向量夹角θ。

7.如权利要求2所述的水表检定方法，其特征在于，所述在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图包括：

基于霍夫变换算法提取所述子表盘，并对所述子表盘进行区域分割；依据所述LRCD色差变换模型对所述子表盘中的红色指针进行提取，得到所述指针二值化图。

8.一种基于图像处理的水表检定装置，其特征在于，包括：

获取模块，用于获取预设表盘图像的指针二值化图；

提取模块，用于从所述指针二值化图中提取图像指针轮廓线，根据所述指针轮廓线计算所有轮廓点到指针圆圆心的第一距离；

第一确定模块，用于所述第一距离中的最大值，定义所述最大值为longest_diatance；

第一计算模块，用于计算每个指针轮廓点到所述指针圆圆心的第二距离；

搜索模块，用于搜索所述第二距离中范围在0.8倍longest_diatance至1倍longest_diatance之间的轮廓点；

求取模块，用于求取所述轮廓点的平均坐标值，并依据所述平均坐标值确定指针针尖位置；

第二确定模块，用于根据所述指针针尖位置与所述指针圆圆心确定指针方向；

判读模块，用于根据所述指针方向对所述预设水表进行示数判读，得到判读结果；

第二计算模块，用于计算所述判读结果并确定所述预设水表的示值状态。

9.如权利要求8所述的水表检定装置，其特征在于，所述获取模块包括：

灰度化处理模块，用于对所述表盘图像进行灰度化处理，得到第一处理结果；

增强模块，用于在所述第一处理结果的基础上对所述表盘图像进行表盘图像增强，得到第二处理结果；

去噪模块，用于在所述第二处理结果的基础上对所述表盘图像进行滤波去噪，得到第三处理结果；

二值化处理模块，用于在所述第三处理结果的基础上对所述表盘图像进行二值化处理，得到第四处理结果；

特征区域提取模块，用于在所述第四处理结果的基础上对所述表盘图像进行特征区域提取，得到第五处理结果；

指针提取模块，用于在所述第五处理结果的基础上对所述表盘图像进行指针提取，得到第六处理结果；

指针二值化处理模块，用于在所述第六处理结果的基础上对所述表盘图像进行指针二值化处理，得到指针二值化图。

10.一种基于图像处理的水表检定系统，其特征在于，包括权利要求8至9任一项所述的水表检定装置。