高 麗

（山西省計量科學(xué)研究院，太原 030032）

0 引言

指針式儀表由于結(jié)構(gòu)簡單、使用方便、安全可靠、易于維護(hù)、價格低廉等優(yōu)勢，被廣泛應(yīng)用于生產(chǎn)、生活等各領(lǐng)域[1]。在指針式儀表的校準(zhǔn)過程中，采用圖像自動識別技術(shù)可以大幅度降低校準(zhǔn)人員的工作強(qiáng)度，提高勞動效率，降低人眼識別誤差。但對識別精度要求較高，如果識別失敗或者識別準(zhǔn)確度不高，需要人工進(jìn)行復(fù)核，反而會加重校準(zhǔn)人員的工作負(fù)擔(dān)。

由于表盤樣式、背景顏色、指針數(shù)量等因素的不同，指針式儀表多種多樣，但其識別流程大體是一致的。主要識別步驟有：儀表區(qū)域定位、刻度區(qū)域識別、指針檢測、極坐標(biāo)轉(zhuǎn)換和計算識別結(jié)果等5 步，但識別過程中的細(xì)節(jié)算法對識別精度有較大的影響。為提高指針式儀表的識別準(zhǔn)確度，本文主要針對儀表盤圓心精確定位算法、指針檢測算法、指針刻度精確識別算法等細(xì)節(jié)，分別進(jìn)行探討。

圖1 指針儀表識別模塊流程圖Fig.1 Flow chart of pointer meter identification module

1 儀表盤圓心精確定位算法

圓形表盤中，儀表刻度一般為圓弧形狀，儀表指針穿過圓心。因此，刻度圓弧圓心的精確定位是決定儀表示值能否正確識別的關(guān)鍵，圓心定位失誤會影響儀表識別的準(zhǔn)確度，甚至造成識別失敗。

對于圓形表盤圓心定位常用的方法是霍夫圓變換，霍夫圓變換是將二維圖像空間中1 個圓轉(zhuǎn)換為該圓半徑、圓心橫縱坐標(biāo)所確定的三維參數(shù)空間中1 個點的過程。因此，圓周上任意3 點所確定的圓，經(jīng)霍夫圓變換后在三維參數(shù)空間應(yīng)對應(yīng)一點。該過程類似于選舉投票過程，圓周上任意3 個點為一選舉人，而這3 個點所確定的圓則為一侯選人（以下稱為候選圓）。遍歷圓周上所有點，任意3 個點所確定的候選圓進(jìn)行投票。

遍歷結(jié)束后，得票數(shù)最高點（理論上圓周上任意3 點確定的圓在霍夫（Hough）變換后均對應(yīng)三維參數(shù)空間中的同一點）所確定的圓即為該圓周上絕大多數(shù)點所確定的圓（以下稱為當(dāng)選圓）。即絕大多數(shù)點均在該當(dāng)選圓的圓周上，以此確定該圓。

在圖1 中可以看到，大數(shù)量的綠色的圓就是通過霍夫圓變換獲得的結(jié)果?；舴驁A檢測的參數(shù)調(diào)整對于結(jié)果的影響很大，檢測階段圓心的累加器閾值越大，判定條件越嚴(yán)苛，對于品質(zhì)不高的儀表照片就可能識別失??；累加器閾值越小，檢測到的圓就越多。而圖1 中紅色的圓弧才是刻度圓弧，圓弧的圓心與指針重合，是需要精確定位的目標(biāo)，需要從眾多的綠色圓中選出目標(biāo)。常見方法是采用置信度和投票的方法實現(xiàn)這一選擇，這種方法對表盤布局和圖片質(zhì)量要求較高，圖2 中的儀表采用此方法就造成了表盤圓心定位誤差。

圖2 圓心定位失敗的儀表Fig.2 Instrument with center positioning failure

圖3 儀表區(qū)域識別Fig.3 Identification of instrument area

本文設(shè)計的方法將儀表區(qū)域定位和刻度區(qū)域識別兩個步驟合并，利用刻度區(qū)域為霍夫變換增加了圓心和半徑的范圍限制，實現(xiàn)表盤圓心的精確定位。這是一種形態(tài)學(xué)處理方法，根據(jù)表盤刻度區(qū)域的定位，估算對應(yīng)的圓心和半徑的約束范圍以擬合弧形刻度線，從而完成對表盤圓心的定位。

首先，需要對儀表區(qū)域圖像進(jìn)行灰度化、二值化處理，然后再進(jìn)行腐蝕，主要是為了刪除文字和其他區(qū)域，保留刻度區(qū)域。然后，再進(jìn)行膨脹操作，目的是刻度線能夠連接在一起，形成最大連通域。

圖3 中可以看出，最大的兩個連通域就是儀表的兩個表盤圓弧。使用ROI 裁剪出這兩個最大連通域如圖4 所示。在這一段刻度線所在的圓弧上，可以知道圓弧左右兩端點A 和B 的坐標(biāo)，根據(jù)圓的定理，圓心應(yīng)該是在這兩端點連成的線段的垂直平分線L2 上，然后通過在一定范圍內(nèi)枚舉圓的半徑，計算枚舉的圓與這一區(qū)域的重疊面積，取最大重疊面積，即可估計圓弧所對應(yīng)的圓心和半徑。

圖4 估計圓心Fig.4 Estimating the center of a circle

由于這兩個數(shù)值是對圖像進(jìn)行膨脹處理之后的結(jié)果，可能不夠精確，將其作為約束條件與霍夫圓變換結(jié)合，可實現(xiàn)儀表盤圓心的精確定位。

2 指針檢測算法

霍夫變換還可以識別圖像中的直線。歐式空間中一條直線上的點在霍夫參數(shù)空間中為一條正弦曲線；歐式空間中同一條直線上的多個點在霍夫參數(shù)空間中為1 個正弦曲線簇且曲線簇相交于一點，稱此點為峰值點。而霍夫參數(shù)空間下的峰值點，則對應(yīng)了歐式空間下的一條直線。 霍夫變換提取的指針直線存在著指針形狀不規(guī)則，出現(xiàn)偏離指針線的情況。

借助精準(zhǔn)的圓心定位，克服了指針擬合偏差問題。本文的指針檢測方法是采用枚舉的方式，如圖6 所示。通過圓心向外延伸構(gòu)造一條線段，將線段順時針在一定夾角范圍內(nèi)進(jìn)行旋轉(zhuǎn)，旋轉(zhuǎn)過程中計算線段與刻度線的相交面積比，即重疊區(qū)域像素點個數(shù)與旋轉(zhuǎn)線段的像素點個數(shù)的比值，取最大值對應(yīng)的線段為最終刻度線所在的線段。

圖6 的儀表中有兩根指針，這是一種常見情況。因此，需要先檢測大指針，然后將大指針刪除，再檢測小指針，以防止大指針干擾小指針的識別。為避免大小指針重合的極端情況，本文設(shè)定一個針對這個指針的相交面積比值的最低閾值，在計算相交面積比值的時候，如果所有角度完成之后發(fā)現(xiàn)沒有超過這個閾值的角度，則認(rèn)為這個指針與大指針重疊了，這時候?qū)⒋笾羔樳€原回去重新尋找正確的指針旋轉(zhuǎn)角即可。

3 指針刻度精確識別算法

圖5 箭頭所指為指針擬合直線產(chǎn)生偏差Fig.5 The arrow indicates the deviation of the pointer fitting line

圖6 指針檢測過程圖Fig.6 Pointer detection process diagram

在圓心和指針的精確定位之后，就可以識別指針刻度了。為了精確識別指針刻度，需要通過極坐標(biāo)系轉(zhuǎn)換，將儀表表盤區(qū)域進(jìn)行清晰展現(xiàn)。將直角坐標(biāo)系的圓心與極坐標(biāo)系的圓心一一對應(yīng)，圓弧儀表盤圓弧可以通過極坐標(biāo)變換到極坐標(biāo)系的一條直線上，實現(xiàn)由圓形到直線的轉(zhuǎn)換。為提高儀表盤清晰度，本文進(jìn)行了相應(yīng)的處理：把刻度區(qū)域圖長寬放大8 倍，由于是canny 邊緣檢測后的圖片，所以放大對圖片沒有很大的影響，但是極坐標(biāo)轉(zhuǎn)換后的圖片刻度變得很清晰。

在獲得刻度區(qū)域之后，采用橫向投影的方式將數(shù)字和刻度線區(qū)域分離。

實現(xiàn)指針式儀表圖像的自動識別讀數(shù)算法，按基本原理分可以歸結(jié)為兩大類：角度法和距離法。

施健等提出了一種指針式化工儀表的遠(yuǎn)程抄表設(shè)計方法，首先利用減影法獲得指針位置，其次利用霍夫變換獲得指針擬合直線的夾角，該方法受光照的影響較大。

董保同提出了一種距離法，即通過計算指針擬合直線和相鄰刻度線的距離得到指針的數(shù)值，此類方法對于表盤內(nèi)部無干擾因素的儀表有較好的識別效果，但對于表盤內(nèi)部存在字符或其他干擾因素的情況下識別效率較低。

圖7 極坐標(biāo)轉(zhuǎn)換效果圖Fig.7 Polar coordinate conversion effect

有兩種方式進(jìn)行指針位置的確定，一種是采用縱向投影法之后尋找最大值坐標(biāo)和兩端點的坐標(biāo)，另一種是縱向投影之后尋找極大值和最大值，計算最大值是在第幾個極大值位置處，即可得到最后的結(jié)果。

縱向投影法，維護(hù)一個x 坐標(biāo)的數(shù)組，掃描圖像每一列，累積這一列中所有非0 像素點的個數(shù)，計入數(shù)組對應(yīng)列的橫坐標(biāo)處。方法一是在統(tǒng)計完投影結(jié)果之后，可以記錄首尾兩邊的橫坐標(biāo)x1，x2，以及上一步驟首尾的數(shù)字num1，num2，然后再記錄投影數(shù)組中的最大值橫坐標(biāo)x0（也就是指針處的橫坐標(biāo)），然后根據(jù)公式（1）得到最后結(jié)果。

圖8 兩種結(jié)果計算方法Fig.8 Two calculation methods of results

另一個方法是在統(tǒng)計完投影結(jié)果之后，尋找圖像中的極大值和極小值點，將極大值認(rèn)為是每一根刻度線，記錄下最大值處為第幾根刻度線?？潭染€總數(shù)N，指針在第n個刻度線處，兩端點數(shù)字num1，num2，則最后結(jié)果通過公式（2）可以得到。

兩種結(jié)果計算方法對比如圖8 所示，第一種方法最后結(jié)果會精確一些。這是由于第二種方法極大值極小值的檢測很容易丟失一兩條刻度線，而且指針在兩根刻度線中間時誤差比較大。

4 結(jié)論

經(jīng)過測試實驗驗證，以上識別算法在筆者承擔(dān)的壓力式溫度計自動校準(zhǔn)系統(tǒng)科研項目中（進(jìn)行了為期半年的測試，取得數(shù)據(jù)200 組，山西省回國留學(xué)人員科研資助項目，項目編號：2016-134），識別率達(dá)到95%，對光線及角度的適應(yīng)大大加強(qiáng)，能作為正式的算法寫入識別程序中。