ИСПОЛЬЗОВАНИЕ АЛГОРИТМА НЕПРЕРЫВНОГО ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ В СИСТЕМАХ ТЕХНИЧЕСКОГО ЗРЕНИЯ
10 В. Е. Махов, А. И. Потапов
УДК 621.762
В. Е. МАХОВ, А. И. ПОТАПОВ
ИСПОЛЬЗОВАНИЕ АЛГОРИТМА НЕПРЕРЫВНОГО ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ
В СИСТЕМАХ ТЕХНИЧЕСКОГО ЗРЕНИЯ
Рассматривается вопрос использования алгоритма непрерывного вейвлетпреобразования в измерительной системе технического зрения. Исследуются возможности применения различных алгоритмов измерения в условиях неравномерности освещения и расфокусировки изображения. Представлен способ повышения точности определения линейных размеров изображения, базирующийся на использовании метода непрерывного вейвлет-преобразования для детектирования границ изображения.
Ключевые слова: системы технического зрения, измерение линейных размеров, вейвлет-преобразование, виртуальный прибор.
Важным вопросом при создании измерительных систем технического зрения остается выбор алгоритма измерения [1]. В реальных условиях трудно получить идеальное изображение контролируемого изделия. На результаты измерений влияют неравномерность освещения, неточность позиционирования объекта контроля, дифракция, дефекты оптики и приемника изображения. В задачах с изменяемыми параметрами качества получаемого изображения подбор и настройки функций алгоритма не всегда обеспечивают предсказуемый результат. Кроме того, использование цифровых камер сравнительно невысокого разрешения (640×480, 1280×1024) совместно с известными алгоритмами измерения не позволяет обеспечить требуемую точность определения размеров изображения контролируемого изделия. Поэтому наибольший интерес представляет разработка помехоустойчивого алгоритма измерений, обладающего высокой точностью, с возможностью определения метрологических характеристик системы технического зрения.
В настоящей статье исследуются возможности системы технического зрения для определения линейных размеров изделий на базе промышленной смарт-камеры NI 1722 фирмы “National Instruments” (NI — США) с использованием алгоритма детектирования границ изображения методом непрерывного вейвлет-преобразования.
Методы вейвлет-анализа [2], в отличие от методов фурье-анализа (частотно-временного анализа), обладают более широкими возможностями исследования непериодичных функций и детектирования их границ (с выделением областей скачкообразных изменений функции) путем динамического сравнения их с короткими масштабируемыми функциями (вейвлетами). Это дает основания полагать о новых возможностях использования вейвлетов для детектирования границ изображения при построении измерительных систем технического зрения.
Исследования проводились на экспериментальной установке, разработанной на кафедре приборов контроля и систем экологической безопасности Северо-Западного государственного заочного технического университета (СЗТУ — Санкт-Петербург). Внешний вид установки представлен на рис. 1.
Для освещения объекта контроля использовалась осветительная система на базе либо лазерного источника света 1, либо лампы накаливания. Дополнительно для ослабления светового пучка применялся светофильтр 2, а для формирования освещающего пучка лучей — конденсор 3 и линза 4. В качестве контролируемого объекта использовалась регулируемая щель монохроматора 5 ( d = 0…0, 4 мм, δd = 0, 001 мм). В состав оптической проецирующей системы входит микрообъектив 6, закрепленный в тубусе 7. В качестве тубуса микроскопа использо-
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Использование алгоритма непрерывного вейвлет-преобразования
11
вался осветитель ОИ-21 отраженного света, который применяется при наблюдении непрозрачных и полупрозрачных объектов в отраженном свете в светлом и темном поле. В схеме использовался эпиообъектив 9××0,20 (длина тубуса 190 мм) с рабочим расстоянием 5,25 мм, что позволяло реализовать подсветку объекта контроля через объектив с помощью дополнительного тубуса 8. В качестве фотоприемника использовалась цифровая смарт-камера 9 (NI 1722). Система освещения и контролируемый объект (щель монохроматора) жестко закреплялись на двухкоординатном столике 10 в рейтере 11 для обеспечения фокусировки микроскопа (наводки на резкость). Оптическая система микроскопа и цифровая камера 9 крепились на отдельном (юстировочном) рейтере 12, что позволяло осуществлять продольную и поперечную подвижку (регулировку по высоте) смарт-камеры, изменяя поле зрения микроскопа.
67
45 23 1
10 11
89 12
Рис. 1
В ходе эксперимента использовалась операционная система Windows 7, установленная на ноутбук (MacBook Pro 15″, 2,4 Гц, 4 Гбайт ОЗУ), среда разработки виртуального прибора — NI LabVIEW 2009, модуль IMAQ Vision, NI Vision Development Module 2009.
Для определения размеров теневого изображения в СЗТУ был разработан исследовательский виртуальный прибор (ВП), структурная схема которого представлена на рис. 2. В виртуальном приборе реализованы следующие измерительные алгоритмы:
— дифференцирование интегрального распределения освещенности (рис. 2, б, блоки 2—10);
— непрерывное вейвлет-преобразование (НВП — CWT — Continuous Wavеlet Transform) интегрального распределения освещенности (рис. 2, в, блоки 11—17);
— определение координат прямой края изображения (рис. 2, г, блоки 18—21) — функция модуля IMAQ Vision (Adv. Straight Edge).
Реализация всех алгоритмов предполагает предварительное усредняющее интегрирова-
{ }ние (блок 4) распределения освещенности в изображении I (xi , y j по высоте в заданной об-
ласти 3:
c+h
∫gi = g(xi ) = I (xi , y j )dy , c
где c — положение области (по координате Y), h — высота области.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
12 В. Е. Махов, А. И. Потапов
ВП осуществляет измерения одновременно по всем указанным алгоритмам (рис. 2, б—г), кроме того, в приборе предусмотрена дополнительная возможность ввода данных с ранее сохраненных файлов изображений (рис. 2, а), что обеспечивает идентичность условий для последующего анализа и определения размеров теневого изображения различными методами при изменении настроек параметров алгоритмов.
а) 1
Ввод изображения
9
d=x(2)–x(1)
б) 2
4
56
78
3 Метод сh
Метод
в) 11 12
Вейвлет
Порог
10
13 14 15 Порог
d=x(2)–x(1)
16 d=x(2)–x(1) 17
d=x(2)–x(1)
20
г) 18 19
21 d=x(2)–x(1)
Параметры
Рис. 2
Для
определения
координат
границ
изображения
функция
{gi
}N i=1
разбивается
на
две
ку-
сочные функции (блоки 8, 15), которые соответствуют правому и левому краям теневого изо-
{ } { } { } { }{ }бражения:
gi(1)
N1 ,
i=1
gi(2)
N2 i=1
( N = N1 + N2 ):
gi
N i=1
=
gi(1)
N1 ∪
i=1
gi(2)
N2 . Для каждой функ-
i=1
ции gik определяются координаты центра тяжести функций (блоки 10, 17) и координаты мак-
симума µk лучшей аппроксимирующей гауссианы [3] (блоки 9, 16):
( ) ( )∫xk = xgk (x)dx
∫ gk (x)dx ;
g
k
[i]
=
ak
exp
⎛ ⎝⎜⎜
( xk
[i] − µk 2σk2
)2
⎞ ⎟⎟⎠
+
ck
,
здесь ak , σk , ck — параметры гауссианы. Линейный размер теневого изображения d (qi ) при заданных настройках qi алгоритма
(см. рис 2, б) определяется как для центра тяжести функций: daverage = x(2) + x(N1) − x(1) , так и
для координат гауссианы: dgauss = µ(2) + x(N1) − µ(1) .
На результат определения координат влияют такие факторы, как методы интегрирования и дифференцирование табличной функции, что вызывает скачки детектирующей функции (шум). Удаление шума достигается компарированием и экстраполяцией (блоки 7, 14).
На рис. 3 представлен пример аппроксимации гауссианой табличных детектирующих функций, полученных методом дифференцирования интегрального распределения освещенности (рис. 3, а, б) и методом НВП (рис. 3, в, г).
Измерения проводились для каждой установленной ширины щели монохроматора при выборе порога отсечки шума (блоки 7, 14 на рис. 2), а также исходя из условия временной устойчивости результата измерения относительно заданного уровня освещения (экспозиции).
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Использование алгоритма непрерывного вейвлет-преобразования
13
Для каждого измерения задавалось положение и высота зоны измерения изображения (область усреднения): c + ∆с и h соответственно ( ∆с = 20 пкс, h = 20 пкс). Характерные зависимости результатов измерения ширины щели d (c) представлены на графиках (рис. 4): кри-
вые 1—3 соответствуют измерениям, проведенным с использованием стандартного алгоритма модуля технического зрения IMAQ Vision — Adv. Straight Edge для сфокусированного и
двух расфокусированных изображений (∆z = +0,5 мм, ∆z = 0, ∆z = –0,5 мм соответственно); кривые 4, 5 — измерениям с использованием алгоритма дифференцирования интегрального распределения освещенности (математическое ожидание и максимум аппроксимирующей гауссианы соответственно); кривые 6, 7 — измерениям с использованием алгоритма НВП (математическое ожидание и максимум аппроксимирующей гауссианы соответственно).
а) Дифференцирование: левая граница Y Исходные данные
б) Дифференцирование: правая граница Y
40
10
Гауссиана
Гауссиана
20
Исходные данные
0 2 4 6 8 10 12 14 16 18 X
в) Y
1500
Вейвлет: левая граница Исходные данные
1000
Гауссиана
500
0 2 4 6 8 10 12 14 16 18 20 22
г) Y
3000
Вейвлет: правая граница Исходные
данные
X
2000 1000
Гауссиана
0 5 10 15 20 X 0 2 4 6 8 10 12 14 16 18 20 22 X Рис. 3
d, пкс
335 6
334 7 2
333
332
13 5
331 4
330
0 40 80 120 160 200 240 280 320 360 400 с, пкс Рис. 4
Как видно из графиков, характер изменения всех кривых достоверно отражает дисторсию оптической системы микроскопа. Кроме того, результат, полученный путем определения
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
14 В. Е. Махов, А. И. Потапов
максимума аппроксимирующей гауссианы коэффициентов НВП, почти совпадает со стандартным алгоритмом Adv. Straight Edge.
Исследования точности определения размеров изображения проводились также посредством изменения его освещенности: изменением экспозиции (E) и коэффициента усиления (G). При использовании методов дифференцирования и НВП результат был стабилен во всем диапазоне изменения параметров, тогда как при использовании стандартного алгоритма наблюдалась существенная разница (0,05 пкс) результатов в зависимости от освещенности изо-
бражения. На трехмерном графике (рис. 5, а) представлена зависимость результата измерения (σ) от координат области усреднения (OX) и освещенности (OY); двухмерный график (рис. 5, б) позволяет оценить точность измерения в диапазоне изменения освещенности изображения объекта контроля.
а) б)
σ, пкс 0,67 0,50 0,34
σ, пкс 0,67
0,50 0,34
1,00 1,96
2,92 G, о.е.
368
0,17 0,17 184 Х, пкс 0
92 184 276 368 Х, пкс
Рис. 5
Таким образом, применение НВП позволяет осуществлять более качественное детекти-
рование краевой границы изображения. Функция CWT-коэффициентов более гладкая, чем
функция, полученная дифференцированием. Ширина пиков (всплесков) детектирующей
функции, полученной посредством НВП, может быть в несколько раз шире зоны перехода
светлого поля в темное. Это обеспечивает более высокую точность определения координат ее
максимумов. Дополнительное использование методов регрессионного анализа позволяет в
несколько раз повысить точность определения максимумов (границ изображения). Координа-
ты максимумов менее зависимы от значения порога компарирования фонового шума, явления
дифракции на границе, изменения уровня подсветки, а также от однородности изображения.
Величина максимумов функции CWT-коэффициентов связана с расфокусировкой оптической
системы. В широком диапазоне изменения освещенности объекта контроля (при наличии
различных дефектов изображения) при различных значениях порога компарирования шума
гарантированная точность измерения составляет 0,2 пкс. Метод НВП менее чувствителен к
настройкам параметров алгоритма, что позволяет сократить их количество и сделать процесс
настройки более прозрачным: это, как следствие, упрощает методику калибровки измери-
тельного прибора.
Метод дифференцирования границ совместно с регрессионным анализом аппроксими-
рующей функции также позволяет повысить точность определения линейного размера (0,5 пкс),
что в два раза выше, чем при использовании стандартных алгоритмов измерений.
Комплексный алгоритм, реализованный в виртуальном приборе контроля и измерения,
иллюстрирует работу стандартных алгоритмов детектирования границ изображения, а также
возможность применения различных математических методов для обработки изображения.
Виртуальный прибор можно использовать с различными цифровыми камерами (USB), на-
пример с цифровыми лупами.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Особенности преобразования дискретных сигналов
15
Результаты проведенных исследований могут быть использованы при построении контрольно-измерительных систем технического зрения для длительных технологических процессов, например формирования силикатных покрытий [4]. Кроме того, представляется целесообразным использовать результаты исследований для создания лабораторных практикумов, например, по дисциплине „Компьютерные технологии в приборостроении“.
СПИСОК ЛИТЕРАТУРЫ
1. Визильтер Ю. В., Желтов С. Ю., Князь В. А. и др. Обработка и анализ цифровых изображений с примерами на LabVIEW IMAQ Vision. М.: ДМК Пресс, 2007. 464 с.
2. Дьяконов В. П. Вайвлеты. От теории к практике. М.: СОЛОН-Р, 2002. 448 с.
3. Химмельблау Д. Анализ процессов статистическими методами: Пер. с англ. М.: Мир, 1973. 957 с.
4. Махов В. Е. Использование алгоритмов вейвлет-анализа в исследовании кинетики формирования порошковообжиговых покрытий // Конструкции из композиционных материалов. 2010. № 3. С. 28—36.
Сведения об авторах Владимир Евгеньевич Махов — канд. техн. наук, доцент; Северо-Западный государственный заочный
технический университет, кафедра приборов контроля и систем экологической безопасности, Санкт-Петербург; E-mail: wmahov@rambler.ru Анатолий Иванович Потапов — д-р техн. наук, профессор; Северо-Западный государственный заочный технический университет, кафедра приборов контроля и систем экологической безопасности, Санкт-Петербург; E-mail: potapov@nwpi.ru
Рекомендована кафедрой приборов контроля и систем экологической безопасности
Поступила в редакцию 03.11.10 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
УДК 621.762
В. Е. МАХОВ, А. И. ПОТАПОВ
ИСПОЛЬЗОВАНИЕ АЛГОРИТМА НЕПРЕРЫВНОГО ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ
В СИСТЕМАХ ТЕХНИЧЕСКОГО ЗРЕНИЯ
Рассматривается вопрос использования алгоритма непрерывного вейвлетпреобразования в измерительной системе технического зрения. Исследуются возможности применения различных алгоритмов измерения в условиях неравномерности освещения и расфокусировки изображения. Представлен способ повышения точности определения линейных размеров изображения, базирующийся на использовании метода непрерывного вейвлет-преобразования для детектирования границ изображения.
Ключевые слова: системы технического зрения, измерение линейных размеров, вейвлет-преобразование, виртуальный прибор.
Важным вопросом при создании измерительных систем технического зрения остается выбор алгоритма измерения [1]. В реальных условиях трудно получить идеальное изображение контролируемого изделия. На результаты измерений влияют неравномерность освещения, неточность позиционирования объекта контроля, дифракция, дефекты оптики и приемника изображения. В задачах с изменяемыми параметрами качества получаемого изображения подбор и настройки функций алгоритма не всегда обеспечивают предсказуемый результат. Кроме того, использование цифровых камер сравнительно невысокого разрешения (640×480, 1280×1024) совместно с известными алгоритмами измерения не позволяет обеспечить требуемую точность определения размеров изображения контролируемого изделия. Поэтому наибольший интерес представляет разработка помехоустойчивого алгоритма измерений, обладающего высокой точностью, с возможностью определения метрологических характеристик системы технического зрения.
В настоящей статье исследуются возможности системы технического зрения для определения линейных размеров изделий на базе промышленной смарт-камеры NI 1722 фирмы “National Instruments” (NI — США) с использованием алгоритма детектирования границ изображения методом непрерывного вейвлет-преобразования.
Методы вейвлет-анализа [2], в отличие от методов фурье-анализа (частотно-временного анализа), обладают более широкими возможностями исследования непериодичных функций и детектирования их границ (с выделением областей скачкообразных изменений функции) путем динамического сравнения их с короткими масштабируемыми функциями (вейвлетами). Это дает основания полагать о новых возможностях использования вейвлетов для детектирования границ изображения при построении измерительных систем технического зрения.
Исследования проводились на экспериментальной установке, разработанной на кафедре приборов контроля и систем экологической безопасности Северо-Западного государственного заочного технического университета (СЗТУ — Санкт-Петербург). Внешний вид установки представлен на рис. 1.
Для освещения объекта контроля использовалась осветительная система на базе либо лазерного источника света 1, либо лампы накаливания. Дополнительно для ослабления светового пучка применялся светофильтр 2, а для формирования освещающего пучка лучей — конденсор 3 и линза 4. В качестве контролируемого объекта использовалась регулируемая щель монохроматора 5 ( d = 0…0, 4 мм, δd = 0, 001 мм). В состав оптической проецирующей системы входит микрообъектив 6, закрепленный в тубусе 7. В качестве тубуса микроскопа использо-
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Использование алгоритма непрерывного вейвлет-преобразования
11
вался осветитель ОИ-21 отраженного света, который применяется при наблюдении непрозрачных и полупрозрачных объектов в отраженном свете в светлом и темном поле. В схеме использовался эпиообъектив 9××0,20 (длина тубуса 190 мм) с рабочим расстоянием 5,25 мм, что позволяло реализовать подсветку объекта контроля через объектив с помощью дополнительного тубуса 8. В качестве фотоприемника использовалась цифровая смарт-камера 9 (NI 1722). Система освещения и контролируемый объект (щель монохроматора) жестко закреплялись на двухкоординатном столике 10 в рейтере 11 для обеспечения фокусировки микроскопа (наводки на резкость). Оптическая система микроскопа и цифровая камера 9 крепились на отдельном (юстировочном) рейтере 12, что позволяло осуществлять продольную и поперечную подвижку (регулировку по высоте) смарт-камеры, изменяя поле зрения микроскопа.
67
45 23 1
10 11
89 12
Рис. 1
В ходе эксперимента использовалась операционная система Windows 7, установленная на ноутбук (MacBook Pro 15″, 2,4 Гц, 4 Гбайт ОЗУ), среда разработки виртуального прибора — NI LabVIEW 2009, модуль IMAQ Vision, NI Vision Development Module 2009.
Для определения размеров теневого изображения в СЗТУ был разработан исследовательский виртуальный прибор (ВП), структурная схема которого представлена на рис. 2. В виртуальном приборе реализованы следующие измерительные алгоритмы:
— дифференцирование интегрального распределения освещенности (рис. 2, б, блоки 2—10);
— непрерывное вейвлет-преобразование (НВП — CWT — Continuous Wavеlet Transform) интегрального распределения освещенности (рис. 2, в, блоки 11—17);
— определение координат прямой края изображения (рис. 2, г, блоки 18—21) — функция модуля IMAQ Vision (Adv. Straight Edge).
Реализация всех алгоритмов предполагает предварительное усредняющее интегрирова-
{ }ние (блок 4) распределения освещенности в изображении I (xi , y j по высоте в заданной об-
ласти 3:
c+h
∫gi = g(xi ) = I (xi , y j )dy , c
где c — положение области (по координате Y), h — высота области.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
12 В. Е. Махов, А. И. Потапов
ВП осуществляет измерения одновременно по всем указанным алгоритмам (рис. 2, б—г), кроме того, в приборе предусмотрена дополнительная возможность ввода данных с ранее сохраненных файлов изображений (рис. 2, а), что обеспечивает идентичность условий для последующего анализа и определения размеров теневого изображения различными методами при изменении настроек параметров алгоритмов.
а) 1
Ввод изображения
9
d=x(2)–x(1)
б) 2
4
56
78
3 Метод сh
Метод
в) 11 12
Вейвлет
Порог
10
13 14 15 Порог
d=x(2)–x(1)
16 d=x(2)–x(1) 17
d=x(2)–x(1)
20
г) 18 19
21 d=x(2)–x(1)
Параметры
Рис. 2
Для
определения
координат
границ
изображения
функция
{gi
}N i=1
разбивается
на
две
ку-
сочные функции (блоки 8, 15), которые соответствуют правому и левому краям теневого изо-
{ } { } { } { }{ }бражения:
gi(1)
N1 ,
i=1
gi(2)
N2 i=1
( N = N1 + N2 ):
gi
N i=1
=
gi(1)
N1 ∪
i=1
gi(2)
N2 . Для каждой функ-
i=1
ции gik определяются координаты центра тяжести функций (блоки 10, 17) и координаты мак-
симума µk лучшей аппроксимирующей гауссианы [3] (блоки 9, 16):
( ) ( )∫xk = xgk (x)dx
∫ gk (x)dx ;
g
k
[i]
=
ak
exp
⎛ ⎝⎜⎜
( xk
[i] − µk 2σk2
)2
⎞ ⎟⎟⎠
+
ck
,
здесь ak , σk , ck — параметры гауссианы. Линейный размер теневого изображения d (qi ) при заданных настройках qi алгоритма
(см. рис 2, б) определяется как для центра тяжести функций: daverage = x(2) + x(N1) − x(1) , так и
для координат гауссианы: dgauss = µ(2) + x(N1) − µ(1) .
На результат определения координат влияют такие факторы, как методы интегрирования и дифференцирование табличной функции, что вызывает скачки детектирующей функции (шум). Удаление шума достигается компарированием и экстраполяцией (блоки 7, 14).
На рис. 3 представлен пример аппроксимации гауссианой табличных детектирующих функций, полученных методом дифференцирования интегрального распределения освещенности (рис. 3, а, б) и методом НВП (рис. 3, в, г).
Измерения проводились для каждой установленной ширины щели монохроматора при выборе порога отсечки шума (блоки 7, 14 на рис. 2), а также исходя из условия временной устойчивости результата измерения относительно заданного уровня освещения (экспозиции).
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Использование алгоритма непрерывного вейвлет-преобразования
13
Для каждого измерения задавалось положение и высота зоны измерения изображения (область усреднения): c + ∆с и h соответственно ( ∆с = 20 пкс, h = 20 пкс). Характерные зависимости результатов измерения ширины щели d (c) представлены на графиках (рис. 4): кри-
вые 1—3 соответствуют измерениям, проведенным с использованием стандартного алгоритма модуля технического зрения IMAQ Vision — Adv. Straight Edge для сфокусированного и
двух расфокусированных изображений (∆z = +0,5 мм, ∆z = 0, ∆z = –0,5 мм соответственно); кривые 4, 5 — измерениям с использованием алгоритма дифференцирования интегрального распределения освещенности (математическое ожидание и максимум аппроксимирующей гауссианы соответственно); кривые 6, 7 — измерениям с использованием алгоритма НВП (математическое ожидание и максимум аппроксимирующей гауссианы соответственно).
а) Дифференцирование: левая граница Y Исходные данные
б) Дифференцирование: правая граница Y
40
10
Гауссиана
Гауссиана
20
Исходные данные
0 2 4 6 8 10 12 14 16 18 X
в) Y
1500
Вейвлет: левая граница Исходные данные
1000
Гауссиана
500
0 2 4 6 8 10 12 14 16 18 20 22
г) Y
3000
Вейвлет: правая граница Исходные
данные
X
2000 1000
Гауссиана
0 5 10 15 20 X 0 2 4 6 8 10 12 14 16 18 20 22 X Рис. 3
d, пкс
335 6
334 7 2
333
332
13 5
331 4
330
0 40 80 120 160 200 240 280 320 360 400 с, пкс Рис. 4
Как видно из графиков, характер изменения всех кривых достоверно отражает дисторсию оптической системы микроскопа. Кроме того, результат, полученный путем определения
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
14 В. Е. Махов, А. И. Потапов
максимума аппроксимирующей гауссианы коэффициентов НВП, почти совпадает со стандартным алгоритмом Adv. Straight Edge.
Исследования точности определения размеров изображения проводились также посредством изменения его освещенности: изменением экспозиции (E) и коэффициента усиления (G). При использовании методов дифференцирования и НВП результат был стабилен во всем диапазоне изменения параметров, тогда как при использовании стандартного алгоритма наблюдалась существенная разница (0,05 пкс) результатов в зависимости от освещенности изо-
бражения. На трехмерном графике (рис. 5, а) представлена зависимость результата измерения (σ) от координат области усреднения (OX) и освещенности (OY); двухмерный график (рис. 5, б) позволяет оценить точность измерения в диапазоне изменения освещенности изображения объекта контроля.
а) б)
σ, пкс 0,67 0,50 0,34
σ, пкс 0,67
0,50 0,34
1,00 1,96
2,92 G, о.е.
368
0,17 0,17 184 Х, пкс 0
92 184 276 368 Х, пкс
Рис. 5
Таким образом, применение НВП позволяет осуществлять более качественное детекти-
рование краевой границы изображения. Функция CWT-коэффициентов более гладкая, чем
функция, полученная дифференцированием. Ширина пиков (всплесков) детектирующей
функции, полученной посредством НВП, может быть в несколько раз шире зоны перехода
светлого поля в темное. Это обеспечивает более высокую точность определения координат ее
максимумов. Дополнительное использование методов регрессионного анализа позволяет в
несколько раз повысить точность определения максимумов (границ изображения). Координа-
ты максимумов менее зависимы от значения порога компарирования фонового шума, явления
дифракции на границе, изменения уровня подсветки, а также от однородности изображения.
Величина максимумов функции CWT-коэффициентов связана с расфокусировкой оптической
системы. В широком диапазоне изменения освещенности объекта контроля (при наличии
различных дефектов изображения) при различных значениях порога компарирования шума
гарантированная точность измерения составляет 0,2 пкс. Метод НВП менее чувствителен к
настройкам параметров алгоритма, что позволяет сократить их количество и сделать процесс
настройки более прозрачным: это, как следствие, упрощает методику калибровки измери-
тельного прибора.
Метод дифференцирования границ совместно с регрессионным анализом аппроксими-
рующей функции также позволяет повысить точность определения линейного размера (0,5 пкс),
что в два раза выше, чем при использовании стандартных алгоритмов измерений.
Комплексный алгоритм, реализованный в виртуальном приборе контроля и измерения,
иллюстрирует работу стандартных алгоритмов детектирования границ изображения, а также
возможность применения различных математических методов для обработки изображения.
Виртуальный прибор можно использовать с различными цифровыми камерами (USB), на-
пример с цифровыми лупами.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9
Особенности преобразования дискретных сигналов
15
Результаты проведенных исследований могут быть использованы при построении контрольно-измерительных систем технического зрения для длительных технологических процессов, например формирования силикатных покрытий [4]. Кроме того, представляется целесообразным использовать результаты исследований для создания лабораторных практикумов, например, по дисциплине „Компьютерные технологии в приборостроении“.
СПИСОК ЛИТЕРАТУРЫ
1. Визильтер Ю. В., Желтов С. Ю., Князь В. А. и др. Обработка и анализ цифровых изображений с примерами на LabVIEW IMAQ Vision. М.: ДМК Пресс, 2007. 464 с.
2. Дьяконов В. П. Вайвлеты. От теории к практике. М.: СОЛОН-Р, 2002. 448 с.
3. Химмельблау Д. Анализ процессов статистическими методами: Пер. с англ. М.: Мир, 1973. 957 с.
4. Махов В. Е. Использование алгоритмов вейвлет-анализа в исследовании кинетики формирования порошковообжиговых покрытий // Конструкции из композиционных материалов. 2010. № 3. С. 28—36.
Сведения об авторах Владимир Евгеньевич Махов — канд. техн. наук, доцент; Северо-Западный государственный заочный
технический университет, кафедра приборов контроля и систем экологической безопасности, Санкт-Петербург; E-mail: wmahov@rambler.ru Анатолий Иванович Потапов — д-р техн. наук, профессор; Северо-Западный государственный заочный технический университет, кафедра приборов контроля и систем экологической безопасности, Санкт-Петербург; E-mail: potapov@nwpi.ru
Рекомендована кафедрой приборов контроля и систем экологической безопасности
Поступила в редакцию 03.11.10 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2011. Т. 54, № 9