МЕТОД РАСПОЗНАВАНИЯ ПОКАЗАНИЙ СТРЕЛОЧНЫХ ИНДИКАТОРОВ
МЕТОДЫ И СРЕДСТВА РЕАЛИЗАЦИИ ОПТИКО-ЭЛЕКТРОННЫХ УСТРОЙСТВ И СИСТЕМ
УДК 004.932
Е. В. ГОЛОВЕНКОВ, С. В. ДЕГТЯРЕВ
МЕТОД РАСПОЗНАВАНИЯ ПОКАЗАНИЙ СТРЕЛОЧНЫХ ИНДИКАТОРОВ
Описана математическая модель распознавания показаний стрелочного индикатора; представлены алгоритмы вычисления показаний стрелочного индикатора и угла отклонения стрелки. Проведены расчет вычислительной сложности описанных алгоритмов и оценка точности разработанного метода.
Ключевые слова: распознавание изображений, контрольно-измерительный прибор, указатель, стрелочный индикатор.
Введение. Результаты исследований Р. Саблатнига [1], Д. Б. Волегова, Д. В. Юрина [2],
Ю. Г. Свинолупова [3], М. С. Ройтмана, В. П. Войтко, Д. Л. Удута, Э. Н. Седова, С. К. Кисе-
лева [4], посвященные автоматизации контроля производства стрелочных индикаторов с по-
мощью оптико-электронных систем, распознающих показания на передней панели прибора,
выявляют высокий уровень развития технических и алгоритмических средств, которые, одна-
ко, не удовлетворяют требованиям быстродействия.
Так, в работе [1] рассмотрен метод распознавания показаний стрелочного индикатора,
основанный на применении модифицированного преобразования Хафа над бинаризованным
изображением указателя. В [2] применяются алгоритмы вычисления контурного изображения
указателя, преобразования в полярную систему координат, свертки полученного изображения
при помощи преобразования Хартли. Наиболее быстродействующим является метод распо-
знавания показаний стрелочного индикатора [3], основанный на вычислении угла наклона
стрелочного указателя с помощью метода наименьших квадратов. Однако перечисленные ме-
тоды обладают высокой вычислительной сложностью, таким образом, актуальной является
задача сокращения времени распознавания показаний стрелочного индикатора.
Математическая модель распознавания показаний стрелочного индикатора со-
держит:
— бинарное изображение шкалы стрелочного индикатора
I (i, j) = {0,1} ; i = 0,h −1, j = 0,w −1 ,
(1)
где h и w — соответственно высота и ширина изображения в пикселах. Массив пикселов
черного цвета I (i, j ) = 1 определяет область на изображении, принадлежащую объекту указа-
теля; массив пикселов белого цвета I (i, j ) = 0 определяет область на изображении, принад-
лежащую иным объектам, например, фону, штриховым и числовым отметкам. Алгоритм получения бинарного изображения аналогичен описанному в работе [3], он характеризуется
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2
Метод распознавания показаний стрелочных индикаторов
37
тем, что в процессе вычитания формируются два изображения I0 (i, j ) , I f (i, j ) , первое из
которых соответствует начальному показанию стрелочного индикатора, второе — конечному;
— функцию преобразования S ( I , ϕ) массива пикселов I (i, j) = 1 в угол отклонения
указателя ϕ ;
— функцию преобразования P (ϕ,V ) угла ϕ в показания индикатора V .
Для определения показаний необходимо распознавать изображения шкалы индикатора. При этом числовое значение показания индикатора определяется на основе двух изображе-
ний I0 (i, j) , I f (i, j) .
Определение угла отклонения указателя. Следующие ниже операции описывают
функцию преобразования S ( I , ϕ) .
Сначала изображение шкалы стрелочного индикатора I (i, j ) сканируется, определяют-
ся координаты O ( x0 , y0 ) центра стрелочного указателя по следующим формулам:
⎧ x0 = ∑ i
⎪⎪ ∀Iij =1; i=0,h−1; j=0,w−1
⎨ ⎪
y0
=
∑
j
⎪⎩ ∀Iij =1; i=0,h−1; j=0,w−1
где n — число черных пикселов на изображении I (i, j) .
n; n,
(2)
Далее изображение шкалы стрелочного индикатора I (i, j) сканируется, определяются
координаты xA , yA , xB , yB центров образованных отрезков по следующим формулам:
⎧ xA = ∑ i
⎪⎪ ∀Iij =1; i=0, y0; j=0,x0
⎨ ⎪
yA
=
∑
j
⎪⎩ ∀Iij =1; i=0, y0; j=0,x0
k; k,
(3)
⎧ xB = ∑ i
⎪⎪ ∀Iij =1; i= y0 ,h−1; j=x0 ,w−1
⎨ ⎪
yB
=
∑
j
⎩⎪ ∀Iij =1; i= y0 ,h−1; j=x0 ,w−1
s; s,
(4)
где k и s — число черных пикселов в первом и во втором отрезках изображения указателя,
причем
k+s= n.
(5)
Угол наклона указателя вычисляется по координатам двух точек A( xA , yA ) , B ( xB , yB )
по формуле:
tg (ϕ) =
yB xB
− −
yA xA
.
(6)
( )Полученные значения tg (ϕ0 ) и tg ϕ f , соответствующие начальному ϕ0 и конечному
ϕ f углу отклонения указателя, запоминаются.
Определение численного значения показаний. По характеру зависимости линейных или угловых расстояний между соседними отметками шкалы от измеряемой величины различают равномерные и неравномерные шкалы. Равномерная шкала характеризуется одинаковым
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2
38 Е. В. Головенков, С. В. Дегтярев
расстоянием между отметками, в то время как на неравномерной шкале расстояние между отметками изменяется по определенному закону.
Необходимо найти функции преобразования P (ϕ,V ) угла ϕ отклонения указателя в
показания индикатора V , т.е. определить вид функции Z (ϕ) = V .
Для случая равномерной шкалы имеем:
Z (ϕ) = κϕ;
( )κ = Z (ϕ1 ) = Z (ϕ2 ) = " = Z ϕq
ϕ1 ϕ2
ϕq
; ϕ1 < ϕ2 < … < ϕq ; q ∈ N ,
(7)
где κ — постоянный коэффициент. Для случая неравномерной шкалы имеем:
Z (ϕi ) = {Z (ϕ1 ), Z (ϕ2 ), …, Z (ϕd )}, ϕi = {ϕ1, ϕ2 ,…, ϕd }, i = 1, d,
Z
(ϕ1 )
ϕ1
≠
Z
(ϕ2
ϕ2
)
≠"
≠
Z
(ϕd
ϕd
);
ϕ1
<
ϕ2
<
…<
ϕd
;
d
∈
N,
(8)
где Z (ϕi ) — дискретная функция, заданная таблично.
Каждому значению ϕi , i = 1, d соответствует показание стрелочного указателя Z (ϕi ) .
Значение d зависит от количества делений l на шкале стрелочного указателя; очевидно
d ≤ l . Угол отклонения указателя относительно начального положения рассчитывается по формуле:
ϕ = ϕ f − ϕ0 ,
(9)
где ϕ0 и ϕ f — угол начального и конечного отклонения указателя.
Таким образом, в зависимости от типа шкалы стрелочного указателя выбирается функ-
ция Z (ϕ) , найденное значение tg (ϕ0 ) сравнивается с заранее определенными табличными
значениями, запоминается начальный угол ϕ0 отклонения указателя, найденное значение
( )tg ϕ f сравнивается с заранее определенными табличными значениями, запоминается ко-
нечный угол ϕ f отклонения указателя, вычисляется угол отклонения указателя по формуле
(9), выполняется преобразование P (ϕ,V ) , значение V запоминается.
Расчет вычислительной сложности. Формула для оценки вычислительной сложности будет иметь следующий вид:
{ }D = d(a) , d(m) , d(d) ,
(10)
где d(a) , d(m) , d(d) — соответственно число операций сложения, умножения, деления.
Вычислительная сложность D(1) разработанного метода рассчитывается по следующим формулам:
D(1) = d (0) + d (1) ,
{ }d (0) = d((a0)) , d((m0)) , d((d0)) ,
(11) (12)
{ }d (1) = d((a1)) , d((m1)) , d((d1)) ,
(13)
где d (0) — вычислительная сложность алгоритма бинаризации; d (1) — вычислительная сложность алгоритмов, описанных выше. По сравнению с существующим методом [3] вы-
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2
Метод распознавания показаний стрелочных индикаторов
39
числительная сложность разработанного метода уменьшена за счет исключения из формул (2)—(6) операций умножения.
Эксперимент. По результатам исследования была написана программа для ЭВМ, моделирующая процесс распознавания показаний стрелочного индикатора предлагаемым и существующим методом [3]. На рисунке проиллюстрировано вычисление угла отклонения указателя; показано окно с результатами распознавания, на котором отображены начальное и конечное положения стрелочного указателя, точки с координатами, вычисленными по формулам (2)—(4), прямые, проходящие через данные точки.
A(xA, yA)
O(xO, yO)
B(xB, yB)
B2(xB2, yB2) o2(xO2, yO2)
A2(xA2, yA2)
ϕ
Экспериментально рассчитанная абсолютная погрешность распознавания показаний стрелочного индикатора разработанным методом по сравнению с методом, приведенным в [3], составляет ∆ ≤ 1° .
СПИСОК ЛИТЕРАТУРЫ
1. Sablatnig R., Kropatsch W.G. Automatic Reading of Analog Display Instruments // Proc. of the 12th IAPR International Conference on Pattern Recognition. Jerusalem, 1994. Vol. 1. P. 794—797.
2. Волегов Д. Б., Юрин Д. В. Считывание показаний стрелочного прибора при помощи web-камеры // ГрафиКон2008. М., 2008. С. 242—247.
3. А. с. 1383242 СССР, МКИ 4 G 01 R 35/00. Способ автоматической поверки стрелочных измерительных приборов и устройство для его осуществления / Ю. Г. Свинолупов, В. П. Войтко, Н. М. Степаненко, Д. Л. Удут. Заявл. 21.10.86; опубл. 23.03.88. Бюл. № 11. 5 с.
4. Киселев С. К. Система автоматизации поверки электроизмерительных приборов // Датчики и системы. 2003. № 6. С. 33—37.
Сведения об авторах
Евгений Владимирович Головенков — аспирант; Юго-Западный государственный университет, кафедра
информационных систем и технологий, Курск;
E-mail: theaswert@yandex.ru
Сергей Викторович Дегтярев
— д-р техн. наук, профессор; Юго-Западный государственный универ-
ситет, кафедра информационных систем и технологий, Курск;
E-mail: sergeyd12@gmail.com
Рекомендована Юго-Западным государственным университетом
Поступила в редакцию 24.10.11 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2
УДК 004.932
Е. В. ГОЛОВЕНКОВ, С. В. ДЕГТЯРЕВ
МЕТОД РАСПОЗНАВАНИЯ ПОКАЗАНИЙ СТРЕЛОЧНЫХ ИНДИКАТОРОВ
Описана математическая модель распознавания показаний стрелочного индикатора; представлены алгоритмы вычисления показаний стрелочного индикатора и угла отклонения стрелки. Проведены расчет вычислительной сложности описанных алгоритмов и оценка точности разработанного метода.
Ключевые слова: распознавание изображений, контрольно-измерительный прибор, указатель, стрелочный индикатор.
Введение. Результаты исследований Р. Саблатнига [1], Д. Б. Волегова, Д. В. Юрина [2],
Ю. Г. Свинолупова [3], М. С. Ройтмана, В. П. Войтко, Д. Л. Удута, Э. Н. Седова, С. К. Кисе-
лева [4], посвященные автоматизации контроля производства стрелочных индикаторов с по-
мощью оптико-электронных систем, распознающих показания на передней панели прибора,
выявляют высокий уровень развития технических и алгоритмических средств, которые, одна-
ко, не удовлетворяют требованиям быстродействия.
Так, в работе [1] рассмотрен метод распознавания показаний стрелочного индикатора,
основанный на применении модифицированного преобразования Хафа над бинаризованным
изображением указателя. В [2] применяются алгоритмы вычисления контурного изображения
указателя, преобразования в полярную систему координат, свертки полученного изображения
при помощи преобразования Хартли. Наиболее быстродействующим является метод распо-
знавания показаний стрелочного индикатора [3], основанный на вычислении угла наклона
стрелочного указателя с помощью метода наименьших квадратов. Однако перечисленные ме-
тоды обладают высокой вычислительной сложностью, таким образом, актуальной является
задача сокращения времени распознавания показаний стрелочного индикатора.
Математическая модель распознавания показаний стрелочного индикатора со-
держит:
— бинарное изображение шкалы стрелочного индикатора
I (i, j) = {0,1} ; i = 0,h −1, j = 0,w −1 ,
(1)
где h и w — соответственно высота и ширина изображения в пикселах. Массив пикселов
черного цвета I (i, j ) = 1 определяет область на изображении, принадлежащую объекту указа-
теля; массив пикселов белого цвета I (i, j ) = 0 определяет область на изображении, принад-
лежащую иным объектам, например, фону, штриховым и числовым отметкам. Алгоритм получения бинарного изображения аналогичен описанному в работе [3], он характеризуется
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2
Метод распознавания показаний стрелочных индикаторов
37
тем, что в процессе вычитания формируются два изображения I0 (i, j ) , I f (i, j ) , первое из
которых соответствует начальному показанию стрелочного индикатора, второе — конечному;
— функцию преобразования S ( I , ϕ) массива пикселов I (i, j) = 1 в угол отклонения
указателя ϕ ;
— функцию преобразования P (ϕ,V ) угла ϕ в показания индикатора V .
Для определения показаний необходимо распознавать изображения шкалы индикатора. При этом числовое значение показания индикатора определяется на основе двух изображе-
ний I0 (i, j) , I f (i, j) .
Определение угла отклонения указателя. Следующие ниже операции описывают
функцию преобразования S ( I , ϕ) .
Сначала изображение шкалы стрелочного индикатора I (i, j ) сканируется, определяют-
ся координаты O ( x0 , y0 ) центра стрелочного указателя по следующим формулам:
⎧ x0 = ∑ i
⎪⎪ ∀Iij =1; i=0,h−1; j=0,w−1
⎨ ⎪
y0
=
∑
j
⎪⎩ ∀Iij =1; i=0,h−1; j=0,w−1
где n — число черных пикселов на изображении I (i, j) .
n; n,
(2)
Далее изображение шкалы стрелочного индикатора I (i, j) сканируется, определяются
координаты xA , yA , xB , yB центров образованных отрезков по следующим формулам:
⎧ xA = ∑ i
⎪⎪ ∀Iij =1; i=0, y0; j=0,x0
⎨ ⎪
yA
=
∑
j
⎪⎩ ∀Iij =1; i=0, y0; j=0,x0
k; k,
(3)
⎧ xB = ∑ i
⎪⎪ ∀Iij =1; i= y0 ,h−1; j=x0 ,w−1
⎨ ⎪
yB
=
∑
j
⎩⎪ ∀Iij =1; i= y0 ,h−1; j=x0 ,w−1
s; s,
(4)
где k и s — число черных пикселов в первом и во втором отрезках изображения указателя,
причем
k+s= n.
(5)
Угол наклона указателя вычисляется по координатам двух точек A( xA , yA ) , B ( xB , yB )
по формуле:
tg (ϕ) =
yB xB
− −
yA xA
.
(6)
( )Полученные значения tg (ϕ0 ) и tg ϕ f , соответствующие начальному ϕ0 и конечному
ϕ f углу отклонения указателя, запоминаются.
Определение численного значения показаний. По характеру зависимости линейных или угловых расстояний между соседними отметками шкалы от измеряемой величины различают равномерные и неравномерные шкалы. Равномерная шкала характеризуется одинаковым
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2
38 Е. В. Головенков, С. В. Дегтярев
расстоянием между отметками, в то время как на неравномерной шкале расстояние между отметками изменяется по определенному закону.
Необходимо найти функции преобразования P (ϕ,V ) угла ϕ отклонения указателя в
показания индикатора V , т.е. определить вид функции Z (ϕ) = V .
Для случая равномерной шкалы имеем:
Z (ϕ) = κϕ;
( )κ = Z (ϕ1 ) = Z (ϕ2 ) = " = Z ϕq
ϕ1 ϕ2
ϕq
; ϕ1 < ϕ2 < … < ϕq ; q ∈ N ,
(7)
где κ — постоянный коэффициент. Для случая неравномерной шкалы имеем:
Z (ϕi ) = {Z (ϕ1 ), Z (ϕ2 ), …, Z (ϕd )}, ϕi = {ϕ1, ϕ2 ,…, ϕd }, i = 1, d,
Z
(ϕ1 )
ϕ1
≠
Z
(ϕ2
ϕ2
)
≠"
≠
Z
(ϕd
ϕd
);
ϕ1
<
ϕ2
<
…<
ϕd
;
d
∈
N,
(8)
где Z (ϕi ) — дискретная функция, заданная таблично.
Каждому значению ϕi , i = 1, d соответствует показание стрелочного указателя Z (ϕi ) .
Значение d зависит от количества делений l на шкале стрелочного указателя; очевидно
d ≤ l . Угол отклонения указателя относительно начального положения рассчитывается по формуле:
ϕ = ϕ f − ϕ0 ,
(9)
где ϕ0 и ϕ f — угол начального и конечного отклонения указателя.
Таким образом, в зависимости от типа шкалы стрелочного указателя выбирается функ-
ция Z (ϕ) , найденное значение tg (ϕ0 ) сравнивается с заранее определенными табличными
значениями, запоминается начальный угол ϕ0 отклонения указателя, найденное значение
( )tg ϕ f сравнивается с заранее определенными табличными значениями, запоминается ко-
нечный угол ϕ f отклонения указателя, вычисляется угол отклонения указателя по формуле
(9), выполняется преобразование P (ϕ,V ) , значение V запоминается.
Расчет вычислительной сложности. Формула для оценки вычислительной сложности будет иметь следующий вид:
{ }D = d(a) , d(m) , d(d) ,
(10)
где d(a) , d(m) , d(d) — соответственно число операций сложения, умножения, деления.
Вычислительная сложность D(1) разработанного метода рассчитывается по следующим формулам:
D(1) = d (0) + d (1) ,
{ }d (0) = d((a0)) , d((m0)) , d((d0)) ,
(11) (12)
{ }d (1) = d((a1)) , d((m1)) , d((d1)) ,
(13)
где d (0) — вычислительная сложность алгоритма бинаризации; d (1) — вычислительная сложность алгоритмов, описанных выше. По сравнению с существующим методом [3] вы-
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2
Метод распознавания показаний стрелочных индикаторов
39
числительная сложность разработанного метода уменьшена за счет исключения из формул (2)—(6) операций умножения.
Эксперимент. По результатам исследования была написана программа для ЭВМ, моделирующая процесс распознавания показаний стрелочного индикатора предлагаемым и существующим методом [3]. На рисунке проиллюстрировано вычисление угла отклонения указателя; показано окно с результатами распознавания, на котором отображены начальное и конечное положения стрелочного указателя, точки с координатами, вычисленными по формулам (2)—(4), прямые, проходящие через данные точки.
A(xA, yA)
O(xO, yO)
B(xB, yB)
B2(xB2, yB2) o2(xO2, yO2)
A2(xA2, yA2)
ϕ
Экспериментально рассчитанная абсолютная погрешность распознавания показаний стрелочного индикатора разработанным методом по сравнению с методом, приведенным в [3], составляет ∆ ≤ 1° .
СПИСОК ЛИТЕРАТУРЫ
1. Sablatnig R., Kropatsch W.G. Automatic Reading of Analog Display Instruments // Proc. of the 12th IAPR International Conference on Pattern Recognition. Jerusalem, 1994. Vol. 1. P. 794—797.
2. Волегов Д. Б., Юрин Д. В. Считывание показаний стрелочного прибора при помощи web-камеры // ГрафиКон2008. М., 2008. С. 242—247.
3. А. с. 1383242 СССР, МКИ 4 G 01 R 35/00. Способ автоматической поверки стрелочных измерительных приборов и устройство для его осуществления / Ю. Г. Свинолупов, В. П. Войтко, Н. М. Степаненко, Д. Л. Удут. Заявл. 21.10.86; опубл. 23.03.88. Бюл. № 11. 5 с.
4. Киселев С. К. Система автоматизации поверки электроизмерительных приборов // Датчики и системы. 2003. № 6. С. 33—37.
Сведения об авторах
Евгений Владимирович Головенков — аспирант; Юго-Западный государственный университет, кафедра
информационных систем и технологий, Курск;
E-mail: theaswert@yandex.ru
Сергей Викторович Дегтярев
— д-р техн. наук, профессор; Юго-Западный государственный универ-
ситет, кафедра информационных систем и технологий, Курск;
E-mail: sergeyd12@gmail.com
Рекомендована Юго-Западным государственным университетом
Поступила в редакцию 24.10.11 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2