УДК 531.7/004.932
ЦИФРОВАЯ ОБРАБОТКА ДИФРАКЦИОННЫХ КАРТИН ФРАУНГОФЕРА ПРИ КОНТРОЛЕ ГЕОМЕТРИЧЕСКИХ ПАРАМЕТРОВ ОБЪЕКТОВ
А.Н. Иванов, В.М. Михайлов
Разработан новый метод цифровой обработки дифракционных картин, который позволит увеличить точность методов дифракционного контроля. Предложено регистрировать координаты точек перегиба главного максимума дифракционной картины, обладающих высокой чувствительностью к изменению размера контролируемого объекта. Ключевые слова: дифракция Фраунгофера, обработка изображения.
Введение
Определение размера объекта по дифракционной картине Фраунгофера заключается в том, что по функции пространственного распределения интенсивности I(х) определяются координаты ее характерных точек - минимумов, максимумов, точек перегиба. Затем определяется коэффициент, связывающий эти координаты с размером объекта. Обычно в качестве характерных точек используют минимумы дифракционной картины, но такой метод сильно ограничен тем, что точно определить координаты этих точек очень трудно - чувствительность сигнала й1(х)/ йх в их окрестности близка к нулю [1]. Именно эта ошибка в определении координат является доминирующей в методах дифракционного контроля и сильно ограничивает их точность. Особенно сильно она сказывается при малых (менее 50 мкм) размерах объекта, когда происходят падение освещенности в дифракционной картине и «размазывание» ее минимумов.
Устранить указанный недостаток можно двумя способами. Первый из них основан на регистрации точек инверсии фазы частотного спектра объекта, координаты которых совпадают с координатами минимумов дифракционной картины. В работе [1] для этого предлагается метод измерения, основанный на последовательном преобразовании Фурье функции амплитудного пропускания объекта и преобразовании Френеля его частотного спектра. Это позволяет уменьшить погрешность измерения на порядок (с 1% до 0,1%). Недостатком данного метода является сложность схем его практической реализации.
Другой способ основан на регистрации характерных точек изображения дифракционной картины после ее преобразования по какому-либо закону [2, 3]. Особый интерес представляет работа [3], в которой предлагается преобразовать главный максимум исходной дифракционной картины в соответствии с выражением I(х)и = 21(х)2п1 -1, где I(х) - распределение интенсивности в дифракционной картине, п -количество преобразований дифракционной картины, х - линейная координата в плоскости регистрации. Определив координаты нулей функции I(х)п, можно вычислить размер контролируемого объекта. Эти нули определяются с высокой точностью, кроме того, использование только центрального максимума позволяет расширить диапазон контролируемых размеров. К сожалению, практическое применение данного метода затруднено тем, что он весьма чувствителен к шуму. Наличие мультипликативного шума приводит к сильному смещению нулей преобразованной функции (ошибка измерения составила 20% при наличии фона в 3 % от величины сигнала).
Анализ главного максимума дифракционной картины от объекта в форме щели
Для решения указанной проблемы предлагается использовать другие характерные точки дифракционной картины - точки перегиба. Исследование дифракционной картины от щели, которая описывается функцией I(и) = Бт(и)2 / и2, где и = к ах/ 2/, а - ширина щели, к = 2 л / X - волновое число, /-фокусное расстояние объектива, показало, что ее главный максимум имеет две такие точки. Чтобы определить их положение, необходимо взять вторую производную от исходной функции I(и):
I" (и) = й21 (и)/йи2. (1)
Приравняв выражение (1) к нулю, найдем значения и, соответствующие положениям точек перегиба в области главного максимума: и2 = 1,303 и и! = -1,303 (рис. 1). Чувствительность функции I (и) к изменению аргумента и в этих точках имеет максимальное значение й!(и) / йи « 0,54.
Если зарегистрировать главный максимум дифракционной картины и определить через его вторую производную й21(х)/ йх2 разность линейных координат точек перегиба Ах = х2 - х1, то она будет связана с разностью А и = и 2 - и1 = 2,606 простой зависимостью Аи = к а Ах /2 / , воспользовавшись которой, можно определить размер объекта а:
а = 2 Аи//к Ах = 5,212 / / к Ах . (2)
А.Н. Иванов, В.М. Михайлов
*к, 1(1/), отн. ед.
Аи
Рис.1. Результат дифференцирования численной модели главного максимума дифракционной картины: 1 - исходное распределение интенсивности, 2 - первая производная, 3 - вторая производная
Для проверки предложенного алгоритма была собрана экспериментальная установка (рис. 2), состоящая из лазера ЛГ-105, спектральной щели с шагом раскрытия 1 мкм, объектива, выполняющего преобразование Фурье, с фокусным расстоянием / = 25 мм и ПЗС-камеры, имеющей диагональ и разрешение 792x576 пикселей.
В ходе эксперимента регистрировались дифракционные картины от щели, ширина которой изменялась в диапазоне от 30 до 55 мкм с шагом 5 мкм. Зарегистрированные максимумы дифракционных картин подвергались предварительной обработке с целью подавления мультипликативного и высокочастотного аддитивного шумов в соответствии с методами, описанными в работе [4].
Рис. 2. Экспериментальная установка для измерения ширины щели дифракционным методом: 1 - лазер, 2 - щелевая диафрагма, 3 - Фурье-объектив, 4 - ПЗС-камера, 5 - компьютер с платой видеозахвата
Затем изображения численно дифференцировались два раза, что позволяло надежно определить координаты их точек перегиба х. Для численного дифференцирования использовалась интерполяционная формула Ньютона [5]
где / - номер элемента ПЗС-приемника, I - интенсивность сигнала в /-м элементе, И - шаг изменения значений (И = 1 пиксель). Результаты численного дифференцирования дифракционной картины приведены на рис. 3, а результаты измерений ширины щели - в таблице. Чтобы показать преимущества предложенного метода, ширина щели была также определена через разность координат минимумов дифракционной картины.
Относительная погрешность измерения предложенным методом оказалась в пределах 0,5-5%. С уменьшением размера объекта погрешность увеличивалась, что обусловлено уменьшением высоты главного максимума и, как следствие, уменьшением отношения сигнал-шум. Искажение формы сигнала приводило к трудности при определении координат точек перегиба. Также источниками погрешности служили ошибка установки щели на нуль, ошибка установки отсчета шкалы щели, разворот дифракционной картины в плоскости приемника.
Экспериментальное исследование
V
- ((Л+1 - I ) - ((+1 - I )2 /2 + ((+1 - Л )3 /3) / И ,
Рис. 3. Результат дифференцирования дифракционной картины от щели шириной 50 мкм: 1 - распределение интенсивности в главном максимуме; 2 - первая производная; 3 - вторая производная
Задаваемое значение ширины щели, мкм Ширина щели, определенная из разности координат точек перегиба, мкм Ширина щели, определенная из разности координат минимумов, мкм
55 55,7 57,8
50 51,4 53,7
45 46,1 46,7
40 41,8 36,6
35 38,7 28,2
30 35,2 20,5
Таблица. Результаты измерения ширины щели
Заключение
Предложен метод цифровой обработки дифракционных картин. Его преимуществом по сравнению с существующими методами является использование точек перегиба главного максимума дифракционной картины, обладающих высокой чувствительностью к изменению размера объекта. Несложные преобразования исходного измерительного сигнала позволяют надежно и с хорошей точностью определить их координаты. Кроме того, центральный максимум менее подвержен влиянию высокочастотного шума, возникающего при рассеянии излучения на микродефектах оптических элементов измерительной схемы. Данный метод удобен для контроля объектов малого размера, когда на площадку ПЗС- приемника укладывается только главный максимум дифракционной картины.
Литература
1. Назаров В.Н., Линьков А.Е. Дифракционные методы контроля геометрических параметров и пространственного положения объектов // Оптический журнал. - 2002. - Т. 69. - № 2. - С. 76-81.
2. Фефилов Г.Д. Оценка погрешности определения моментов экстремума в дифракционной картине Эйри при воздействии аддитивной помехи // Научно-технический вестник СПБГУИТМО. - 2006. -№ 31. - С. 55-58.
3. Власов Н.Г., Кулиш С.М. Нелинейная цифровая обработка дифракционных картин // Мир техники кино. - 2008. - № 8. - С. 41-42.
4. Грейсух Г.И., Ежов Е.Г., Земцев А.Ю., Степанов С.А. Разработка методов и программных средств подавления шумов в интерферограммах на этапе их предварительной обработки // Компьютерная оптика. - 2005. - Т. 28. - С. 140-144.
5. Исаков В.Н. Элементы численных методов. - М.: ЛСЛБЕМЛ, 2003. - 192 с.
Иванов Александр Николаевич - Санкт-Петербургский государственный университет информационных
технологий, механики и оптики, кандидат технических наук, доцент, i_off@mail.ru
Михайлов Виталий Михайлович - Санкт-Петербургский государственный университет информационных
технологий, механики и оптики, студент, mynameisvitalik@gmail.ru