CC BY
77
МЕТРОЛОГИЧЕСКИЙ КОНТРОЛЬ И ИЗМЕРЕНИЯ МЕХАНИЧЕСКИХ ВЕЛИЧИН ОПТИЧЕСКИМИ МЕТОДАМИ
УДК 531.7.082: 004.932
А. Н. Иванов, Ю. А. Каракулев, В. М. Михайлов
АЛГОРИТМ ИЗМЕРЕНИЯ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ ОБЪЕКТА ПО ЕГО МУАР-ИНТЕРФЕРЕНЦИОННОЙ КАРТИНЕ
Рассмотрен метод обработки муаровых полос, основанный на использовании параметрического уравнения, связывающего изменение формы муаровых полос с геометрической формой контролируемого объекта. Предложена методика предварительной обработки изображения для подавления шумов.
Ключевые слова: муаровая картина, обработка изображений.
Введение. Современные методы контроля геометрических параметров объектов должны обладать высокой точностью, быстродействием и возможностью регистрации контролируемых параметров в режиме реального времени. Во многом эти характеристики зависят от выбора метода регистрации получаемого изображения.
Для формирования изображения в современных приборах широко используются ПЗС-приемники, преобразующие изображение в цифровой сигнал, который необходимо расшифровать. Процесс расшифровки включает:
— определение зависимости между изменением контролируемого параметра объекта и изменением регистрируемого параметра изображения;
— предварительную обработку изображения в целях устранения искажений;
— регистрацию изменения выбранного параметра полученного изображения.
Модель обработки муар-интерференционной картины. Для измерения, например,
диаметра металлических цилиндров в работе [1] предложено использовать метод, основанный на анализе муар-интерференционной картины на щели между эталонной отражающей поверхностью и образующей цилиндра. Пример такой картины, полученной на цилиндре диаметром 13 мм, приведен на рис. 1. Распределение амплитуды на щели описывается выражением
и(юх, у') « А бшс (к юх а(у) / 2)соБ((к юх а(у) + к Аа(у) 0) / 2 - к&2х Б 0 /4), (1)
где = х' / г — пространственная частота; а(у) = а0 + А а(у) — функция, описывающая изменение ширины щели (а) между зеркалом и краем объекта; к — волновое число; Б — диаметр цилиндра; 0 — угол падения волны на щель; А = а(у)/-^Тк , у = у'.
Рассматривая функции втс(-) и соб(-) выражения (1) как пространственные амплитудно-частотные решетки, можно получить параметрическое уравнение муаровых полос:
р = (2 Аа( у) 0-А,)/2 ^, (2)
где р — порядок полосы.
При наклоне цилиндра относительно зеркала ширина щели будет меняться по линейному закону А а = а у, где а — угол наклона. Раскрывая уравнение (2), получаем
у' = ((х')2 В0/г2 + Л,(2р +1))/2а0 . (3)
Из выражения (3) видно, что измерив координаты х' и у' центра полосы, можно определить диаметр цилиндра. К сожалению, вычислить координаты полос непосредственным » Л V ф 1Ш измерением невозможно. Невозможно определить и поря-
А. 1 V 1 ♦ Жл док регистрируемых полос. Кроме того, сама муаровая I | Т ' картина требует предварительной обработки для удаления
А^Д ^^ щр , из нее фоновых дифракционных полос. Поэтому авторами ^ I ж ^^ " ^^ ш» была разработана методика, позволяющая вычислить
диаметр цилиндра относительным методом.
Сначала из муаровой картины удалялся фон. Для в __ , этого выполнялась свертка полученного изображения с
одномерной функцией Гаусса, ширина купола которой определялась из выражения ох = 3 V, где V — частота фоновых полос [2]. На следующем этапе из муаровой картины удалялся высокочастотный аддитивный шум. Согласно выражению (2) ширина муаровых полос определяется как
5 = у(р +1)-у(р) = Х/а0 . (4)
Выражение (4) демонстрирует возможность удаления из фурье-спектра изображения всех гармоник, частоты которых превышают граничную частоту vГр = 2п / 5 . Пример обработанной таким образом муаровой картины приведен на рис. 2.
Рис. 1
х
х 2
х' 3
у 1 у 2 у 3 Рис. 2
Затем было получено три сечения зарегистрированной муаровой картины с координатами х{, х2 и х3 таким образом, что Ах ' = х3 - х2 = х2 - х{ (см. рис. 2). По этим сечениям определялось смещение точки минимума муаровой полосы порядка р: Ау2 = у3 - у2 и Ау1 = у2 - у1. Положение минимума полосы определялось как центр тяжести заштрихованной фигуры (рис. 3). Для этого сначала приближенно оценивались положение точки минимума (т) по наименьшей интенсивности сигнала и положение границ фигуры в пикселах изо-
1
Алгоритм измерения параметров объекта по его муар-интерференционной картине
35
бражения: n1 « m - (S /41) и n2 « m + (S /41), где t — размер одного пиксела. Затем приближенно оценивалась интенсивность сигнала в максимуме полосы порядка p — /max. Тогда центр тяжести фигуры определялся по формуле
m = ] ((/max / 2 - /(j)) j) ] (/max / 2 - /( j)), (5)
j=ni / j=ni где / (j) — интенсивность сигнала в j-м пикселе сечения.
т
0 п1 т п2 ]
Рис. 3
Так как вычислить абсолютные координаты точек минимумов полос невозможно, определим диаметр объекта через разности этих координат. Для этого составим две системы уравнений, связывающих координаты точек минимумов, воспользовавшись для этого выражением (5):
у3 = (х3)2 Б/2г2а + Л,(2р + 1)/2а0, у2 = (х2)2 Б / 2 г2 а + Л (2 р + 1)/2а0; у2 = (х2)2 Б/2г 2а + Л (2 р + 1)/2а0, у[ = (х{)2 Б/2г 2а + Л (2 р + 1)/2а0. Вычитая уравнения, входящие в системы, друг из друга, получаем
Ау2 = у3 - у2 = Ах'(х3 + х2)Б/2г2а ; (6)
Ау[ = у2 -у[ = Ах'(х2 + х1)Б/2г2а, (7)
а вычитая из уравнения (6) уравнение (7) —
Ау2 - Ау{ = Ах' (х3 - х{) Б /2 г2а.
С учетом того, что х3 - х{ = 2Ах', упрощаем последнее выражение:
А у2 -А у{ = (А х')2 Б / г2 а ,
Б = (Ау2 -А у!) г2 а /(Ах ')2. (8)
Таким образом, определив по полученному изображению разности координат трех точек, можно относительным способом определить диаметр цилиндра.
Эксперимент. Для экспериментальной оценки предложенного метода был собран макет измерительной установки (рис. 4). В качестве контролируемого объекта был использован калиброванный цилиндр диаметром 13 мм с неотражающей поверхностью, изготовленный с допуском Ь9.
Регистрация муаровых полос производилась цифровой камерой со встроенным объективом. Поэтому для того чтобы оценить размер пиксела изображения, на матовый экран наносился тест-объект. Расстояние от объекта до матового экрана выбиралось исходя из усло-
вия наблюдения дифракции Фраунгофера: z > a / X .
Г
Формирователь лазерного пучка
А
Л
Цифровая камера
а )
///'//////
Рис. 4
На данной установке были проведены измерения диаметра цилиндра по приведенной выше методике. Погрешность измерения составила порядка 3 %. Такая достаточно большая погрешность может быть объяснена погрешностью измерения расстояния г и существенной ошибкой измерения координат точек минимумов полосы. Метод, основанный на определении минимума как центра тяжести фигуры, теоретически дает погрешность измерения сдвига полосы порядка 0,8 %, однако на практике эта погрешность в значительной степени зависит от уровня шума. Погрешность измерения расстояния г обусловлена трудностью определения положения эквивалентных плоскостей. Эта ошибка может быть уменьшена путем использования цилиндрического фурье-объектива. В этом случае положение эквивалентных плоскостей относительно линзы не оказывает влияния на параметры муаровой картины.
Заключение. Предложенный метод обработки муаровых полос основан на использовании параметрического уравнения, связывающего изменение формы муаровых полос с геометрической формой контролируемого объекта. Экспериментальным путем доказана возможность применения данного метода для измерения диаметров цилиндров.
Работа выполнена при финансовой поддержке Правительства Санкт-Петербурга, грант № 28-04/18.
z
список литературы
1. Назаров В. Н., Иванов А. Н. Использование муарового эффекта для создания высокоточных дифракционных систем контроля геометрических параметров объектов // Тр. VIII Междунар. конф. „Прикладная оптика — 2008". СПб: ГОИ им. С. И. Вавилова, 2008. Т. 1. С. 96—100.
2. Грейсух Г. И., Ежов Г. И., Земцов А. Ю., Степанов С. А. Разработка методов и программных средств подавления шумов в интерферограммах на этапе их предварительной обработки // Компьютерная оптика. 2005. Т. 28.
Сведения об авторах
Александр Николаевич Иванов — канд. техн. наук, доцент; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра компьютеризации и проектирования оптических приборов; E-mail: [email protected]
Юрий Александрович Каракулев — канд. техн. наук, доцент; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра компьютеризации и проектирования оптических приборов
Применение дисперсионного элемента при настройке интерферометров
37
Виталий Михайлович Михайлов
студент; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра компьютеризации и проектирования оптических приборов; E-mail: [email protected]
Рекомендована кафедрой компьютеризации и проектирования оптических приборов
Поступила в редакцию 26.04.11 г.
УДК 681.4.07
А. М. Бурбаев, В. В. Зюзин, М. Г. Лукин, С. А. Марусов, В. И. Мищенко
ПРИМЕНЕНИЕ ДИСПЕРСИОННОГО ЭЛЕМЕНТА ПРИ НАСТРОЙКЕ ИНТЕРФЕРОМЕТРОВ ТИПА ЦЕНДЕРА — МАХА
Изложены принципы спектральной интерферометрии, рассмотрен механизм формирования муаровых полос при введении дисперсионного элемента на выходе интерферометра. Показано, как использование предлагаемого метода позволяет расширить динамический диапазон регистрации интерференционных полос при поиске нулевой разности хода, а также упростить и ускорить процесс настройки интерферометра на нулевую разность хода интерферирующих пучков.
Ключевые слова: принципы спектральной интерферометрии, интерферометр Цендера — Маха, муаровые полосы, волновой цуг, длина когерентности.
Введение. Интерференционные методы измерений и исследований находят широкое применение в науке и технике благодаря высокой чувствительности и точности. Интерферометры используются для точных измерений длин и углов, угловых размеров звезд, характеристик оптических поверхностей, показателей преломления сред, спектрального состава исследуемого излучения [1].
Наибольшее распространение получили интерферометры, работающие по принципу двухлучевой интерференции. К классическим интерферометрам этой группы относятся интерферометры Майкельсона, Тваймана — Грина, Жамена, Цендера — Маха и Рождественского.
Предложенный почти одновременно (в 1891 и 1892 гг.) и независимо друг от друга Цен-дером и Махом интерферометр, принципиальная схема которого приведена на рис.1, обладает рядом преимуществ по сравнению с интерферо- ;
метрами Майкельсона и Жамена. Прежде всего, а ч это возможность получения контрастной и хоро- т шо освещенной интерференционной картины (с полосами заданной ширины) в плоскости объекта, установленного на произвольном расстоянии от зеркал интерферометра [2]. Данное обстоятельство обусловило возможность использования интерферометра Цендера — Маха для измерения
показателей преломления газов, для исследования а'2
неоднородностей воздушного потока в аэродина- Рис.1
мической трубе, для изучения распределения плотностей газовых потоков вокруг моделей, исследуемых в аэродинамических установках, для изучения температурного поля объекта, для измерения напряжений в моделях при наличии статических и динамических нагрузок [3].
2
