CC BY
64
2. Гуров И.П. Компьютерная фотоника: принципы, проблемы и перспективы II Научно-технический вестник СПбГУ ИТМО. - 2005. - Вып. 21. - С. 5-20.
3. Marquet P., Rappaz B., Magistretti P.J., Cuche E., Emery Y., Colomb T., Depeursinge C. Digital holographic microscopy: a noninvasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy II Optics Letters. - 2005. - V. 30. - № 5. - P. 468-470.
4. Bernhardt I., Ivanova L., Langehanenberg P., Kemper B., von Bally G. Application of digital holographic microscopy to investigate the sedimentation of intact red blood cells and their interaction with artificial surfaces II Bioelectrochemistry. - 2008. - № 73. - P. 92-96.
5. Langehanenberg P., Ivanova L., Bernhardt I., Ketelhut S., Vollmer A., Dirksen D., Georgiev G., von Bally G., Kemper B. Automated three-dimensional tracking of living cells by digital holographic microscopy II Journal of biomedical optics. - 2009. - V. 14. - № 1. - P. 014018-1-014018-7.
6. Schnars U., Jueptner W. Digital holography. - Berlin, Heidelberg: Springer, 2005. - 164 p.
7. Балтийский С.А., Гуров И.П., Де Никола С., Коппола Д., Ферраро П. Современные методы цифровой голографии II Проблемы когерентной и нелинейной оптики I Под ред. И.П. Гурова и С. А. Козлова. -СПб: СПбГУ ИТМО, 2004. - С. 91-117.
8. Горбатенко Б.Б., Гребенюк А.А., Максимова Л.А., Перепелицына О.А., Рябухо В.П. Спекл-фотография и голографическая интерферометрия с цифровой записью дифракционного поля в Фурье-плоскости II Компьютерная оптика. - 2010. - Вып. 34. - № 1. - С. 69-81.
Гребенюк Константин Александрович - Саратовский государственный университет им. Н.Г.Чернышевского,
кандидат физ.-мат. наук, доцент, KonstantinAG@yandex.ru Гребенюк Антон Александрович - Саратовский государственный университет им. Н.Г.Чернышевского,
аспирант, GrebenyukAA@yandex.ru Рябухо Владимир Петрович - Саратовский государственный университет им. Н.Г.Чернышевского,
доктор физ.-мат. наук, профессор, rvp-optics@yandex.ru
УДК 535.4
ОПРЕДЕЛЕНИЕ МАЛЫХ СМЕЩЕНИЙ ПОВЕРХНОСТИ ОБЪЕКТОВ МЕТОДОМ ЦИФРОВОЙ ГОЛОГРАФИИ В.Г. Гендин
Представлены результаты исследования точности определения малых смещений поверхности объектов методом цифровой голографии. Рассматривается процесс регистрации цифровых голограмм и восстановления фазовой информации из голограмм. Сравнивается точность определения нормальных смещений поверхности объекта при использовании двух методов устранения нулевого порядка дифракции с последующим вычислением фазы предметной волны.
Ключевые слова: цифровая голография, восстановление изображений, устранение нулевого порядка дифракции.
Введение
Формирование и анализ изображений трехмерных объектов - актуальная задача в различных областях науки и техники. При этом особенно важно получение количественных оценок геометрических характеристик исследуемых объектов. Известно, что наиболее высокую точность определения геометрических характеристик обеспечивают интерферометрические методы. Однако методы классической интерферометрии применимы только к объектам с оптически гладкими поверхностями и относительно малыми отклонениями от правильной геометрической формы. Исследования объектов произвольной формы с интерферометрической точностью возможны при использовании голографических методов, в том числе методов цифровой голографии [1, 2]. Отличительная особенность голографических методов состоит в том, что голограмма содержит полную информацию об оптическом поле предметной волны, при этом из голограммы можно извлечь всю необходимую фазовую и амплитудную информацию о предметной волне, характеризующей объект.
В настоящей работе представлены результаты исследования точности определения малых смещений поверхности объектов методами цифровой голографии при записи цифровых голограмм с использованием интерферометра Майкельсона. При этом задача устранения нежелательного нулевого порядка дифракции решается двумя методами: с помощью фильтрации высоких частот в голограмме и с использованием метода управляемого фазового сдвига, когда для определения смещения поверхности записываются две или более голограмм в различных положениях объекта и далее проводится восстановление информации об изменениях фазы предметной волны.
Запись цифровых голограмм с использованием интерферометра Майкельсона
Запись голограмм в цифровой голографии осуществляется посредством регистрации с помощью видеокамеры суммы предметной волны, отраженной от объекта, и опорной волны. Для реконструкции
изображения из голограммы используются вычислительные методы, позволяющие определить значения амплитуды и фазы предметной волны [1, 2].
Запись цифровых голограмм можно осуществить с использованием схемы, представленной на рис. 1.
^^ 4
1
1=1-
2
а
Рис. 1. Схема интерферометра Майкельсона для записи цифровых голограмм: 1 - лазерный источник;
2 - светоделитель; 3 - исследуемый объект; 4 - опорное зеркало; 5 - ПЗС-видеокамера;
6 - вычислительная система; 7 - пьезоэлектрический преобразователь
Излучение лазерного источника 1 при помощи светоделителя 2 разделяется на предметную волну, освещающую исследуемый объект 3, и опорную волну, которая отражается от зеркала 4. Голограмма регистрируется при помощи видеокамеры 5. Реконструкция трехмерного изображения из голограммы выполняется в вычислительной системе 6.
Опорное зеркало 4 закреплено на площадке с микрометрическими винтами, которые позволяют задавать угол между направлениями распространения опорной и предметной волн для реализации вне-осевой записи голограмм. Это необходимое условие для точного восстановления амплитудной и фазовой информации, так как в обычной осевой голограмме нулевой порядок дифракции, действительное и мнимое изображения наложены друг на друга, что затрудняет выделение мнимого изображения, характеризующего объект. Объект 3 закреплен на пьезоэлектрическом преобразователе 7, который позволяет задавать малые смещения объекта. При проведении экспериментов использовалась ПЗС-видеокамера 5 с выделением полезной области размером 320*240 пикселей. В качестве источника излучения использовали стабилизированный Ие-№ лазер 1 с длиной волны 632 нм.
Пример голограммы объекта, полученной с использованием описанной схемы, представлен на рис. 2. Как отмечено выше, голограмма содержит полезную информацию о мнимом изображении объекта и две неинформативные составляющие, обусловленные наличием нулевого порядка дифракции и изображения-двойника (действительного изображения). Для выделения полезной информации из цифровой голограммы используются различные методики и вычислительные алгоритмы, рассматриваемые далее.
Рис. 2. Голограмма нерассеивающего объекта (5*5 мм)
Восстановление изображений из голограммы
Восстановленное дифрагированное поле в плоскости записи голограммы на расстоянии ё от плоскости голограммы описывается интегралом Рэлея-Зоммерфельда [1, 2]
2
3
7
5
6
Q(%,л) = exp|'dy-j J JR(x,y)I(x,y)exp i-d[-x)2 + (л-y)
dxdy, (1)
где ^ - координаты в плоскости изображения; х, у - координаты в плоскости голограммы; X - длина волны; d - расстояние от плоскости голограммы до плоскости изображения; Я(х,у) - комплексная амплитуда опорной волны; I (х, у) - распределение интенсивности в плоскости голограммы. Интеграл (1) в цифровой голографии интерпретируется как двумерное преобразование Фурье опорной волны, объектной волны и квадратичной фазовой функции пространственного распространения волны:
л)=-Ъ exp(' -f J exp(' ^(^2+л2) |x
ад ад / чх 2 \
J J R(x,y)I(X,y)expl i (X2 + y2)J expl--^(§x + ^y)Jdxdy.
— ГГ, —Ю ^ ' ^ '
(2)
Выражение (2) можно представить в форме
®(У%) = ТьА 6ХР('Тd)еХР(V + 0]х3+1 [[,У)I(х,уМх,у)], (3)
где 3+1[[(х,у)1 (х,у)^(х,у)] - прямое преобразование Фурье; V.,- пространственные частоты, соответствующие координатам в плоскости изображения (£, п),
V. = , (4)
5 Xd
. (5)
Xd
Из формул (4) и (5) видно, что пространственные частоты зависят от длины волны и расстояния между плоскостями.
В цифровой голографии требуется выполнение критерия Найквиста, а именно, на каждый период наивысшей пространственной частоты в голограмме должны приходиться не менее двух отсчетов, т.е. светочувствительных ячеек (пикселей) оцифрованного изображения. Кроме этого, размер пикселя Ax ограничивает допустимое значение максимального угла amax при интерференции сферических вторичных волн от каждой точки объекта и опорной волны согласно условию _ X
а max _ ~ . .
2Ax
Точное восстановление распределения амплитуды и фазы объектной волны достигается только при условии соблюдения ограничения на расстояние d согласно соотношению
, , NAx
d > dc _-,
c X
где N - количество пикселей в видеокамере (для координаты x). При d _ dc возможно точное восстановление не только амплитуды, но и фазы, при этом размер пикселя восстановленного изображения равен размеру пикселя голограммы.
В настоящей работе при записи голограмм соблюдалось ограничение на значение максимального угла при внеосевой записи голограммы, а восстановление проводилось для расстояния d = 27 см, рассчитанного в соответствии с приведенными выше формулами и обеспечивающего точное восстановление фазы предметной волны.
Определение малых смещений поверхности объекта
Для определения малого нормального смещения поверхности объекта достаточно записать две голограммы (в исходном и конечном состоянии). При этом восстанавливаются фазы обеих голограмм с использованием соотношения
Ф& л) _ arctgl ImQ(^ Л)
^(0(5, л), после чего вычисляются изменения фазы Дф(5, л) =Ф2(5, л) "Ф1(5, л).
Величина смещения в единицах длины определяется по очевидной формуле
X
Дd = —Дф. 4л
Для точного определения фазы предметной волны нужно вначале выделить информацию о мнимом изображении, а для этого, прежде всего, необходимо устранить нулевой порядок дифракции.
Устранение нулевого порядка дифракции
В работе исследовались два метода устранения нулевого порядка дифракции - с помощью фильтрации высоких частот и с использованием метода управляемого фазового сдвига.
Поскольку нулевой порядок дифракции - это непродифрагировавшая на низкочастотных участках голограммы волна, то устранить нулевой порядок дифракции можно путем подавления низких частот в голограмме. Изображение восстанавливается из голограммы без нулевого порядка дифракции в результате фильтрации в исходной голограмме высоких частот с малой частотой среза.
В методе фазового сдвига [1, 2] регистрируют три или более голограмм, смещенных относительно друг друга по фазе, для каждого состояния объекта. В настоящей работе использовался известный квадратурный алгоритм с записью четырех голограмм, сдвинутых относительно друг друга по фазе на л / 2 . Разность фаз задавалась пьезоэлектрическим преобразователем 7 (рис. 1). При этом набор голограмм можно представить в виде
1(х, у) = (1/4 Я*){1 (х, у; а = 0) -I (х, у; а = л) +1[1 (х, у; а = л/2) -I (х, у; а = 3л/2)]}, (6)
где значок * обозначает комплексно-сопряженную величину, и полученные в (6) значения подставляются в правую часть формулы (3) вместо величины I(х, у).
а б
Рис. 3. Результат восстановления амплитуды предметной волны (а) и фазы предметной волны (б).
Размер изображений 10*10 мм
На рис. 3 представлен пример восстановления амплитуды и фазы предметной волны без нулевого порядка дифракции из голограммы (рис. 2). Объект представляет собой металлическую пластинку с выдавленным на ней символом (в виде римской цифры II).
Экспериментальное исследование точности определения малых смещений объекта
В схеме (рис. 1) с помощью пьезоэлектрического преобразователя 7 задавали малые смещения объекта и определяли эти смещения по изменениям фаз - с использованием рассмотренных выше методик при устранении нулевого порядка дифракции путем фильтрации высоких частот и с помощью метода фазового сдвига. Смещение задавалось в диапазоне 5-150 нм с шагом 10 нм. Погрешность позиционирования пьезоэлектрического преобразователя составляет 0,3 нм, т.е. пренебрежимо мала по сравнению с задаваемыми смещениями.
На рис. 4, 5 представлены графики зависимости отклонения измеренных значений смещения от заданных значений для двух методов устранения нулевого порядка дифракции. При этом рассчитанное смещение представлено как среднее значение смещения по всей плоскости объекта.
Из графиков следует, что среднее квадратическое отклонение (СКО) определения смещения для различных методов различно: для метода с использованием фильтрации СКО полученных смещений от заданных оно составляет 3,3 нм, а для метода фазового сдвига - 2,3 нм. Таким образом, метод цифровой голографии позволяет определять смещения поверхности объекта с погрешностью в единицы нанометров, что соответствует интерферометрической точности. Различные значения погрешностей при разных смещениях объясняются воздействием внешних факторов, таких как микровибрации. Точность измерений можно дополнительно повысить при обеспечении требуемой стабильности оптической схемы и устранении влияния внешних факторов.
10
н
О -6 -1-1-1-1-1-1-1-1-1-1-1-1-1-л
5 15 25 35 45 55 65 75 85 95 105 115 125 135 145 155 Заданный сдвиг, нм
Рис. 4. График зависимости отклонения величины смещения от заданного значения при устранении нулевого порядка дифракции методом фильтрации высоких частот
о
и
о и и
S н о S н о
(D
s и
(D
И
si
W
10 8 6 4 2 0 -2 -4
15 25 35 45 55 65 75 85 95 105 115 125 135 145 155 Заданный сдвиг, нм
Рис. 5. График зависимости отклонения величины смещения от заданного значения при устранении нулевого порядка дифракции методом фазового сдвига
Заключение
Проведенные исследования показали, что метод цифровой голографии позволяет определять геометрические характеристики объектов с погрешностью в единицы нанометров. Применение метода фильтрации высоких частот и метода фазового сдвига для устранения нулевого порядка дифракции дают сопоставимые по точности результаты. Метод фильтрации высоких частот более прост в реализации. Однако метод фазового сдвига обеспечивает несколько более высокую точность, несмотря на необходимость записи и обработки нескольких голограмм, т.е. проведения более длительного эксперимента.
Работа выполнена при финансовой поддержке Министерства образования и науки Российской Федерации.
Литература
1. Schnars U., Jueptner W. Digital Holography: Digital hologram recording, numerical reconstruction, and related techniques. - Berlin: Springer-Verlag, 2005. - 164 c.
2. Балтийский С.А., Гуров И.П., Де Никола С., Коппола Д., Ферраро П. Современные методы цифровой голографии // Проблемы когерентной и нелинейной оптики. - СПб: СПбГУ ИТМО, 2004. - C. 91-117.
Гендин Владислав Геннадьевич
- Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики, аспирант, vlad.gendin@gmail.com
