Например, Бобцов

ИССЛЕДОВАНИЕ МУАР-ИНТЕРФЕРЕНЦИОННОЙ КАРТИНЫ, ВОЗНИКАЮЩЕЙ ПРИ ДИФРАКЦИИ СВЕТА НА ЗАЗОРЕ МЕЖДУ «ТОЛСТЫМ» КРАЕМ И ПЛОСКОЙ ОТРАЖАЮЩЕЙ ПОВЕРХНОСТЬЮ

ИССЛЕДОВАНИЕ МУАР-ИНТЕРФЕРЕНЦИОННОЙ КАРТИНЫ …
УДК 531.7.082.5:535.42/44
ИССЛЕДОВАНИЕ МУАР-ИНТЕРФЕРЕНЦИОННОЙ КАРТИНЫ, ВОЗНИКАЮЩЕЙ ПРИ ДИФРАКЦИИ СВЕТА
НА ЗАЗОРЕ МЕЖДУ «ТОЛСТЫМ» КРАЕМ И ПЛОСКОЙ ОТРАЖАЮЩЕЙ ПОВЕРХНОСТЬЮ
А.Н. Иванов
Рассмотрена проблема извлечения информации о геометрических параметрах контролируемого объекта из муар-интерференционной картины, которая сформирована наложением частотных спектров этого объекта. Предложена конкретная методика для ее решения. Ключевые слова: дифракция, муар, измерения.
Введение
При интерференционных измерениях информация о геометрических параметрах контролируемых объектов содержится в интерференционной картине – интерферограмме. Некоторую качественную информацию об этих параметрах можно получить по виду этой интерферограммы. Однако в большинстве случаев наибольший интерес представляет количественная информация. При этом для автоматизации контроля, увеличения точности и скорости получения результатов интерферограмма регистрируется ПЗС-матрицей и затем обрабатывается на компьютере. Для оценки измеряемой характеристики необходимо на стадии численного моделирования определить зависимость, связывающую изменение этой характеристики с изменением параметров интерференционной картины. Другим условием получения точных результатов является предварительная обработка интерферограмм, направленная на снижение шумов в ней без потери полезной информации об исследуемом объекте.
Формирование муар-интерференционных полос
d

А а
θ
А’

Б Б’

Рис. 1. Схема формирования дифракционной картины от «толстого» края
В работе [1] было предложено использовать для контроля геометрических параметров объекта дифракционную картину от зазора, образованного краем этого объекта и эталонной плоской отражающей поверхностью. Проведенные исследования показали, что при дифракции на таком зазоре в дифракционной картине формируется дополни-
16 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)

А.Н. Иванов

тельная система муар-интерференционных полос. Эти полосы образуются наложением

частотных спектров, сформированных при дифракции света на краях объекта А и Б

(рис. 1). Так как прямая и отраженная волны падают на края объекта с разными по зна-

ку углами падения, то при смещении этих краев между их частотными спектрами появ-

ляется фазовый сдвиг, пропорциональный величине смещения

ψ = k θ Δa / 2 + k θ Δb / 2 ,

(1)

где k – волновое число, θ – угол падения волны на край, Δа и Δb – смещения краев А и

Б. Если объект наклонить относительно поверхности зеркала, то разность фаз будет из-

меняться вдоль зазора, и в дифракционной картине возникнет дополнительная система

полос равной ширины.

В случае, когда толщина края d ≈ 0 (тонкий край), распределение амплитуды от

зазора в дальней области будет иметь вид

U (ωx , y ') = Asinc(k ωx a( y) / 2) cos((k ωx a( y) + k Δa( y) θ) / 2) ,

(2)

где ωx = x '/ z – пространственная частота, z – расстояние от зазора до плоскости реги-

страции. Рассматривая множители выражения (2) как пространственные решетки, уда-

лось получить параметрическое уравнение полос p = (2 Δ a( y) θ − λ) / 2 λ и определить,

что они ориентированы строго перпендикулярно поверхности зеркала [1]. Когда d ≠ 0 , мы имеем «толстый» край. Полагая Δa = Δb и пренебрегая срезани-

ем части волнового фронта, проходящего через зазор, получим выражение для дальней области:

U (ωx , y ') ≈ Asinc(k ωx a( y) / 2) cos((k ωx a( y) + k Δa( y) θ) / 2 − k ω2x d / 4) .

(3)

Наличие «толстого» края приводит к частичному срезанию отраженного от зерка-

ла волнового фронта в соответствии с выражением a ' =a − d θ . Мы не учитываем это,

так как анализ выражений (2) и (3) показал, что равномерный сдвиг края не приводит к изменению параметров муаровых полос.

Рис. 2. Муар-интерференционная картина от зазора, образованного краем толщиной 0,9 мм и плоской отражающей поверхностью

Продольное смещение краев А и Б приводит к фазовому сдвигу между частотны-

ми спектрами ϕ = k d ω2x / 2 [2] и искривлению муаровых полос (рис. 2). Параметрическое уравнение полос в этом случае имеет вид

p = (2 Δa( y) θ − ω2x d − λ) / 2 λ .

(4)

Наклоняя объект относительно поверхности зеркала, получаем линейный закон

изменения ширины зазора Δa = α y , где α − угол наклона. Тогда уравнение (4) примет

вид

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)

17

ИССЛЕДОВАНИЕ МУАР-ИНТЕРФЕРЕНЦИОННОЙ КАРТИНЫ …

y = (ω2x d + λ (2 p +1)) / 2 α θ .

(5)

Из выражения (5) видно, что, определив координаты центра полосы и порядок по-

лосы p, можно определить толщину объекта.

Определение толщины объекта по муар-интерференционной картине

Анализ полученных муаровых картин показывает, что определить непосредственным измерением координаты полос невозможно. Невозможно определить и порядок регистрируемых полос. Кроме того, сама муаровая картина требует предварительной обработки для удаления из нее фоновых дифракционных полос. Поэтому была разработана методика, позволяющая определить нужный параметр косвенным методом.
x’

x’1 1 x’2 2 x’3 3

y’1 y’2 y’3

y’

Рис. 3. Вид муар-интерференционной картины после ее обработки низкочастотным фильтром с функцией пропускания типа Гаусса

Рис. 4. Определение минимума полосы как центра тяжести фигуры
Сначала из муаровой картины удалялся фон. Для этого выполнялась свертка полученного изображения с одномерной функцией Гаусса, ширина купола которой определялась из выражения σx ≥ 3ν , где ν − частота фоновых полос [3]. На следующем эта-
18 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)

А.Н. Иванов

пе из муаровой картины удалялся высокочастотный аддитивный шум. Так как, согласно (5), ширина муаровых полос равна

S = y( p +1) − y( p) = λ / α θ ,

то из Фурье-спектра изображения можно удалить все гармоники, частоты которых бу-

дут больше граничной частоты νгр = 2π / S . Пример обработанной таким образом муа-

ровой картины приведен на рис. 3.

Затем было выполнено три сечения зарегистрированной муаровой картины с ко-

ординатами x '1 , x '2 и x '3 таким образом, что Δ x ' = x '3 − x '2 = x '2 − x '1 (рис. 3). По этим сечениям определялось смещение минимума муаровой полосы порядка p:

Δ y'2 = y'3 − y'2 и Δ y'1 = y'2 − y'1 . Положение минимума полосы определялось как центр

тяжести заштрихованной фигуры, изображенной на рис. 4. Для этого сначала прибли-

женно оценивалось положение минимума по наименьшей интенсивности сигнала и по-

ложение границ фигуры в пикселях изображения: n1 ≈ m − (S / 4t) и n2 ≈ m + (S / 4t) , где t – размер одного пикселя. Затем приближенно оценивалось интенсивность сигнала в

максимуме полосы порядка p − Imax. Центр тяжести фигуры определялся по формуле

n2 n2
min = ∑ ((Imax / 2 − I ( j)) j) / ∑ (Imax / 2 − I ( j)) ,

j =n1

j =n1

где I ( j) – интенсивность сигнала в j-м пикселе сечения.

Так как определить абсолютные координаты точек минимумов полос невозможно,

попытаемся определить толщину объекта через разности этих координат. Для этого со-

ставим две системы уравнений, связывающих координаты точек минимумов, восполь-

зовавшись для этого выражением (5):

⎪⎧ ⎨ ⎪⎩

y y

'3 '2

= =

x '32 x '22

d d

/ /

2 2

z2α θ z2α θ

+ +

λ (2 λ (2

p p

+ 1) + 1)

/ /

2α θ, 2α θ;

⎪⎧ ⎨ ⎪⎩

y y

'2 '1

= =

x '12 x '12

d d

/ /

2 2

z2α θ + z2α θ +

λ (2 λ (2

p p

+ 1) + 1)

/ /

2α θ, 2α θ.

.

Вычитая уравнения, входящие в системы, друг из друга, получаем:

Δ y '2 = y '3 − y '2 = Δx '(x '3 + x '2 ) d / 2 z2α θ ,

(6)

Δ y '1 = y '2 − y '1 = Δ x '(x '2 + x '1) d / 2 z2α θ , и, вычитая из уравнения (6) уравнение (7), получаем:

(7)

Δ y '2 − Δ y '1 = Δ x '(x '3 − x '1) d / 2 z2α θ .

С учетом того, что x'3 −x'1 = 2Δx , упрощаем выражение:

Δ y '2 − Δ y '1 = Δ x '2 d / z2α θ .

(8)

Таким образом, удалось связать разность координат центра муаровой полосы с

толщиной края.

Экспериментальная часть

Для экспериментальной оценки предложенного метода был собран макет измерительной установки, изображенный на рис. 5. В роли контролируемого объекта был использован калиброванный цилиндр диаметром 13 мм с неотражающей поверхностью, изготовленный с допуском h9. Согласно методу эквивалентных диафрагм [4], его воздействие на падающую волну аналогично действию двух тонких краев, смещенных от-

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)

19

ИССЛЕДОВАНИЕ МУАР-ИНТЕРФЕРЕНЦИОННОЙ КАРТИНЫ …

носительно друг друга на расстояние d = D θ , где D – диаметр цилиндра. Поэтому,
оценив по муаровым полосам величину d, можно перейти от нее к диаметру цилиндра. Регистрация муаровых полос производилась цифровой камерой со встроенным
объективом. Поэтому, чтобы оценить размер пикселя изображения в плоскости матового экрана, на последний наносился эталонный объект. Расстояние от объекта до мато-
вого экрана выбиралось из условия наблюдения дифракции Фраунгофера z ≥ a2 / λ .

z d

56

14

a
2
3
Рис. 5. Схема экспериментальной установки для измерения диаметров цилиндров: 1 – лазер, 2 – формирователь пучка, 3 – плоская отражающая поверхность, 4 – контролируемый цилиндр, 5 – матовый экран, 6 – цифровая камера
На данной установке были проведены измерения диаметра цилиндра по приведенной выше методике. Погрешность измерения составила порядка 3 %. Такая достаточно большая погрешность может быть объяснена ошибкой установки дистанции δz
и достаточно большой ошибкой измерения координат минимумов полосы. Хотя метод, основанный на определении минимума как центра тяжести фигуры, теоретически дает ошибку измерения сдвига полосы порядка 0,8 %, на практике эта ошибка сильно зависит от уровня шума. Ошибка установки дистанции обусловлена трудностью определения положения эквивалентных плоскостей. Эта ошибка может быть уменьшена использованием цилиндрического объектива, выполняющего преобразование Фурье. В этом случае положение эквивалентных плоскостей относительно линзы не оказывает влияния на параметры муаровой картины.
Заключение
Предложен дифракционный метод контроля толщин объектов, основанный на использовании муаровых полос, которые возникают при наложении частотных спектров краев объекта. Их использование позволяет значительно увеличить чувствительность измерений.
Получено параметрическое уравнение муаровых полос, позволившее связать координаты полосы с толщиной края. Так как абсолютные координаты полос не могут быть определены при измерениях, была разработана методика, позволяющая использовать для измерения толщины объекта разность координат нескольких точек центра полосы.
Для проверки правильности сделанных выводов был собран макет установки для измерения диаметров цилиндров предложенным методом. Получено хорошее совпадение между результатами численного моделирования и эксперимента. Доказана возможность применения предложенной методики для оценки толщины объекта. Увеличение точности предложенного метода может быть достигнуто путем применения более точного метода определения координат минимумов муаровых полос.

20 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)

С.А. Миронов

Работа выполнена при финансовой поддержке Правительства Санкт-Петербурга в виде гранта № 28-04/18.

Литература

1. Иванов А.Н., Назаров В.Н. Использование явления муара для увеличения точности дифракционных методов контроля геометрических параметров и пространственного положения объектов // Оптический журнал. – 2009. – Т. 76. – № 1. – С. 46–50.
2. Гудмен Дж. Введение в Фурье-оптику. – М.: Мир, 1970. 3. Грейсух Г.И., Ежов Г.И., Земцов А.Ю., Степанов С.А. Разработка методов и про-
гаммных средств подавления шумов в интерферограммах на этапе их предварительной обработки // Компьютерная оптика. – 2005. – Т. 28. 4. Лемешко Ю.А, Чугуй Ю.В. Размерный контроль круговых отражающих цилиндров интерференционным методом // Автометрия. – 2003. – Т. 39. – № 5. – С. 42–52.

Иванов Александр Николаевич

– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кандидат технических наук, доцент, i_off@mail/ru

Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)

21