Реконструкция изображения по пространственному распределению интенсивности дифракционного спекл-модулированного поля

РЕКОНСТРУКЦИЯ ИЗОБРАЖЕНИЯ ПО ПРОСТРАНСТВЕННОМУ РАСПРЕДЕЛЕНИЮ ИНТЕНСИВНОСТИ ДИФРАКЦИОННОГО СПЕКЛ-МОДУЛИРОВАННОГО ПОЛЯ

Б. Б. Горбатенко, Л. А. Максимова, В. П. Рябухо, Ю. В. Норов Саратовский государственный технический университет

Саратовский государственный университет Институт проблем точной механики и управления РАН

Аннотация

Рассматриваются физические механизмы записи и восстановления информации об амплитудно-фазовой структуре оптического поля в голографическом процессе. Обсуждаются способ создания голограммноподобной структуры и способ восстановления изображения по записи интенсивности спекл-модулированного дифракционного поля рассеивающего объекта без опорного пучка. Анализируются результаты натурных, численных экспериментов и численного моделирования процессов синтеза голограммноподобной структуры и восстановления изображения плоских рассеивающих объектов симметричной и несимметричной форм.

Введение

Восстановление изображения объекта по зарегистрированному в дифракционной зоне распределению интенсивности рассеянного когерентного поля представляет интерес в голографии и дифракционной оптике, в оптических измерениях и диагностике, в оптических методах и устройствах записи и обработки информации, в когерентной микроскопии и т. д. [1-3]. Основная проблема решения этой задачи состоит в том, что при регистрации пространственного распределения интенсивности без использования когерентного опорного пучка теряется информация о пространственном распределении фазы объектного поля. Поэтому задача восстановления изображения объекта по записи интенсивности дифракционного поля с математической точки зрения является некорректной [1-2]. Решение этой задачи сводится к нахождению способа реконструкции фазы дифракционного поля при записи интенсивности. Известны разные подходы к решению этой проблемы [4-13].

В [14-18] показано, что для плоских рассеивающих объектов распределение интенсивности, по поверхности которых описывается четной функцией координат, или, другими словами, для объектов, обладающих вращательной симметрией четного порядка, требуемое восстановление пространственного распределения фазы дифракционного поля может быть получено достаточно простым способом. Если поверхность такого объекта служит источником 5-коррелированного диффузно-рассеянного когерентного излучения, то в дальней зоне дифракции формируется спекл-модулированное поле с бинарным пространственным распределением фазы - в пределах каждого спекла фаза постоянна, а при переходе к соседнему спеклу фаза меняется на п радиан [14, 15].

Данное обстоятельство определяет простой закон пространственного распределения фазы дифракционного поля и, следовательно, простую процедуру реконструкции этого фазового распределения. Следующим шагом будет простой алгоритм восстановления изображения объекта по записи

распределения интенсивности дифракционного поля [16-18].

В настоящей работе представлены результаты детального анализа процессов восстановления изображений по записи спекл-структур дифракционного поля рассеивающих объектов различной формы. Рассматриваются натурные (аналоговые) и компьютерные эксперименты восстановления изображений по записи цифровых Фурье-спеклограмм.

Голограммная структура

Достоинством голографических методов записи и восстановления информации является полное восстановление информации об амплитудно-фазовом распределении объектного поля. Голограммную структуру, - зарегистрированную картину интерференции объектной и опорной волн, обычно представляют в виде суперпозиции синусоидальных дифракционных решеток или зонных структур Френеля [19, 20].

Такие представления с формальной точки зрения позволяют вполне адекватно описать процессы записи и восстановления волнового поля в голографии. Однако реальная голограммная структура, фактически, представляет собой то, что можно увидеть при ее микроскопическом исследовании (рис. 1): набор элементарных дифракционных ячеек, которыми являются спеклы объектного поля, промоду-лированные несущими интерференционными полосами. При переходе от спекла к спеклу интерференционные полосы смещаются, поскольку фаза поля в спеклах разная. Контраст интерференционных полос определяется соотношением амплитуд объектной и опорной волн в каждом спекле. Следовательно, амплитуда объектного поля кодируется значением контраста интерференционных полос, а фаза - их пространственным положением. Поперечные размеры спеклов определяются шириной углового (пространственного) спектра объектного поля.

Восстанавливаемое с голограммы объектное волновое поле представляет собой суперпозицию элементарных волн, дифрагированных на элементарных дифракционных ячейках голограммы. Такие элементарные дифракционные решетки отчетливо

1РН

ММК

прослеживаются на рис. 1, на котором представлен сильно увеличенный фрагмент реальной голо-граммной структуры рассеивающего объекта.

ЛМи/И

I " * 1 ¡11111 1\

№1, 1ЫН 11 РЛ

11 НИ;,А

х .11.!! Ни

Рис. 1. Фрагмент голограммной структуры рассеивающего объекта

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

Вышеописанный процесс показывает, как с помощью реальной голограммы восстанавливается амплитудно-фазовая структура объектного поля в плоскости голограммы, а значит, и в дифракционном поле. На рис. 2 схематически показан процесс формирования действительного голографического изображения в дифракционном поле в результате суперпозиции элементарных волн, дифрагированных на элементарных дифракционных ячейках голограммы.

Таким образом, реальная голограмма представляет собой сложную дифракционную структуру и состоит из множества элементарных дифракционных решеток, служащих наименьшими информационными ячейками голограммы. При этом пространственное положение этих решеток на голограмме полностью определяется спекл-структурой объектного поля. Полное волновое поле, восстановленное с голограммы, формируется в результате суперпозиции всех элементарных волн, дифрагировавших на ячейках.

Вышеуказанные особенности дифракционной структры голограммы предопределяют простой способ синтеза голограммы по записи интенсивности объектного дифракционного поля без использования опорного пучка света. В каждом спекле записанной структуры необходимо тем или иным способом создать несущие полосы с относительным сдвигом, соответствующим фазе поля в этом спекле.

Такая процедура восстановления изображения обсуждается в [16-18]. Она основана на записи спек-лограммы рассеивающего объекта в дальнем поле дифракции - Фурье-спеклограммы и на бинарном распределении фазы в спеклах этого поля для объектов симметричной формы.

Рис. 2. Схематическое представление процесса восстановления с голограммы действительного изображения объекта при использовании пучка света, сопряженного опорному: 1 - восстанавливающая волна; 2 - голограмма; 3 - элементарные дифракционные ячейки;

4 - элементарные волны, дифрагирующие на элементарных ячейках; 5 - действительное изображение объекта

Схемы записи цифровых Фурье-спеклограмм

Для практического решения задачи восстановления изображения по распределению интенсивности дифракционного поля необходимо зарегистрировать цифровую спеклограмму в дальней зоне дифракции. Для записи цифровых спеклограмм нужно, чтобы регистрируемая спекл-структура уверенно разрешалась ПЗС-матрицей - размеры спеклов должны в несколько раз превышать размеры пикселов ПЗС-матрицы. При использовании Фурье-преобразующей линзы для записи Фурье-спеклограммы это достигается путем уменьшения размера объекта - источника спекл-модулированного поля с использованием линзы с достаточно большим фокусным расстоянием Наименьший поперечный размер спеклов в этом случае

определяется соотношением [21] „ /

где < - наибольший размер апертуры объекта - источника спекл-поля.

В экспериментах по записи безлинзовой Фурье-спеклограммы это достигается путем увеличения расстояния от рассеивающего объекта до плоскости регистрации и уменьшения размеров объекта. В этом случае размер спеклов определяется расстоянием г от объекта до плоскости регистрации

X г

<Лп

Оптические схемы записи линзовой и безлинзовой Фурье-спеклограмм представлены на рис. 3. В безлинзовой схеме записи для более точного выполнения условия дальнего поля дифракции и выполнения Фурье-преобразования объектного поля необходимо использовать квазисферическую освещающую волну, сходящуюся в плоскости регистрации (рис. 3б). В этих схемах записи цифровых Фурье-спеклограмм предполагается использование цифровых фотокамер с открытой ПЗС-матрицей.

а)

1 3 4

-»--

->--

] \ 5

Рис. 3. Схемы записи линзовой (а) и безлинзовой (б) Фурье-спеклограмм рассеивающего объекта:

1 - освещающий лазерный пучек; 2 - рассеиватель;

3 - апертурная диафрагма (транспарант);

4 - ПЗС-матрица; 5 - линза

При использовании фотокамер с несъемным объективом, этот объектив можно применять в качестве Фурье-преобразующей оптической системы (рис. 3 а). Однако, в силу малого, как правило, фокусного расстояния фотографических объективов, регистрируемый объект должен иметь достаточно малые размеры для формирования на ПЗС-матрице разрешаемых спеклов. Например, для объектива с /и100 мм, размер объекта й0 должен быть ё0 < X / / е± и 4 мм при использоваии ПЗС-матрицы с 8 мкм размерами пикселов.

Для регистрации Фурье-спеклограмм с несъемным объективом возможно использование специальных оптических схем с дополнительным рассеивающим экраном или с линзой. На рис. 4 и 5 представлены возможные схемы записи таких Фурье-спеклограмм. Лазерный пучок 1 освещает рассеиватель 2 и транспарант 3. Согласно схеме на рис. 4 цифровой фотоаппарат 6 сфокусирован на рассеивающий экран 5, изображение которого формируется на сопряженной плоскости ПЗС-матрицы. Эта схема записи имеет два недостатка. Интенсивность картины изображения падает с расстоянием от центра ПЗС-матрицы из-за ограниченности индикатрисы рассеяния излучения на экране 5. Поскольку в

этой схеме присутствует второй рассеиватель - экран 5, то в плоскости ПЗС-матрицы образуется двойная спекл-структура 7 (рис. 4). Каждый крупный спекл в этой структуре оказывается промоду-лированным более мелкими спеклами.

Рис. 4. Схема регистрации спеклограммы

с помощью цифрового фотоаппарата с несъемным объективом с использованием рассеивающего экрана: 1 - лазерный пучок; 2 - рассеиватель; 3 - транспарант, 4 - продольный разрез спекл-структуры; 5 - рассеивающий экран; 6 - цифровой фотоаппарат с несъемным объективом; 7 - схематический вид спекл-структуры в плоскости ПЗС-матрицы фотоаппарата

1 з __и £

2 т 2!

Рис. 5. Схема регистрации Фурье-спеклограммы при использовании коллективной линзы с помощью цифрового фотоаппарата с несъемным объективом: 1 - лазерный пучок; 2 - рассеиватель; 3 -транспарант;

4 - продольный разрез спекл-структуры; 5 -линза; 6 - цифровой фотоаппарат с несъемным объективом

Размер крупных спеклов определяется угловыми размерами апертуры транспаранта 3 и коэффициентом увеличения объектива фотоаппарата р = г,/ г0.

X 2Х

где 6 и г - угол, под которым наблюдается источник излучения в плоскости экрана; г 0 - расстояние от рассеивающей поверхности рассеивателя 5 до объектива (передней главной плоскости) фотоаппарата 6, - расстояние от объектива до ПЗС-матрицы.

Размер мелких спеклов, которыми промодулиро-ваны крупные спеклы, определяется формулой

X г,

(1„

где - диаметр апертуры объектива фотоаппарата.

Чтобы исключить модулирование спеклов спек-лами и увеличить равномерность распределения средней интенсивности спекл-структуры в плоско-

±

0

сти ПЗС-матрицы в схеме на рис. 4 можно заменить рассеиватель 5 на коллективную линзу 5. Такая схема представлена на рис. 5. Фотоаппарат 6 сфокусирован на плоскость апертуры линзы, которая в этом случае оптически сопряжена с плоскостью ПЗС-матрицы фотоаппарата. Линза 5 собирает рассеянное лазерное излучение и перенаправляет его в объектив (изображение источника спекл-поля формируется в плоскости апертуры диафрагмы объектива). Изображение поперечного сечения спекл-структуры в плоскости линзы отображается на ПЗС-матрице фотоаппарата.

Цифровое изображение спекл-структуры, создаваемой источником в форме кольца, зарегистрированной с помощью схемы на рис. 5, представлено на рис. 6.

Рис. 6. Цифровое изображение спекл-структуры,

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

Кольцевая область с повышенной интенсивностью на рис. 6а обусловлена проявлением нерассеянной компоненты лазерного излучения, прошедшей через рассеиватель 2 (рис. 5). Размер спеклов в плоскости ПЗС-матрицы в этом случае определяется по формуле (3).

Дифракция дальнего поля на амплитудной Фурье-спеклограмме

Пространственное распределение интенсивности объектного дифракционного поля несет частичную информацию об изображении объекта. Для рассеи-

вающего объекта отдельные спеклы дифракционного поля при голографической регистрации служат элементарными дифракционными ячейками, несущими информацию о фазе и амплитуде объектной волны. В отсутствии опорной волны эти информационные ячейки утрачивают фазовую информацию и, как следствие, утрачивается возможность восстановления полной информации о волновом поле и изображении объекта.

На рис. 7 представлены примеры Фурье-образов Фурье-спеклограмм рассеивающих объектов различной формы. Нетрудно видеть соответствие дифракционных картин автокорреляционным распределениям средней интенсивности поля в плоскости объекта.

+

д

Ш IV

Рис. 7. Фурье-образы Фурье-спеклограмм рассеивающих объектов различной формы: кольцевого квадрата

Формальный анализ, численные и экспериментальные результаты показывают, что при записи амплитудной Фурье-спеклограммы картина дифракции дальнего поля на такой спеклограмме отображает автокорреляцию распределения интенсивности поля в плоскости объекта. Аналогично, при цифровой записи Фурье-спеклограммы ее Фурье-образ отображает автокорреляцию распределения интенсивности поля в плоскости объекта.

Восстановление фазовой информации по спеклограмме

В [14,15] теоретически и экспериментально показано, что для 5 -коррелированных рассеивающих объектов с осесимметричным распределением интенсивности I (х0, у0) разность фаз поля в дальней области дифракции в соседних спеклах с наибольшей вероятностью равна п радиан. Следовательно, пространственный сдвиг интерференционных полос в соседних спеклах интерференционной структуры Фурье-голограммы на рис. 1 с наибольшей вероятностью равен половине периода полос Л . Таким образом, имеем следующий алгоритм для восста-

а

в

г

е

новления фазы объектного поля и, следовательно, для восстановления изображения объекта по зарегистрированной интенсивности дальнего дифракционного поля.

Пусть в эксперименте записывается цифровая Фу-рье-спеклограмма Iп) без использования опорного пучка света. Затем, с использованием графических цифровых технологий в этом распределении интенсивности в пределах спеклов наносятся искусственные интерференционные полосы с некоторым периодом Л. При этом, тем или иным способом, обеспечивается сдвиг полос на полпериода Л при переходе от одного спекла к соседнему. Этим сдвигом полос восстанавливается информация о фазовом пространственном распределении в объектном поле. Фурье-преобразование искусственной голограммы, полученной вышеописанным способом, должно привести к формированию изображения объекта.

На рис. 8 приведен результат численного Фурье-преобразования цифровых Фурье-спеклограмм объектов, имеющих осевую симметрию в распределении средней интенсивности I(х0,у0) (крест (рис. 86),

квадрат (рис. 8в), эллипс (рис. 8г)). Для сравнения на рис. 8а приведен результат численного Фурье-преобразования цифровой Фурье-голограммы, полученной с использованием опорного пучка света.

наносить картины несущих полос в спеклах с требуемым сдвигом. Пользователь этой программы вручную применял бинарную цветовую гамму (красный и зеленый цвета) для определения фазы объектного поля в спеклах: красный цвет - 0 рад, зеленый - п рад. Затем программа автоматически определяла и наносила систему несущих интерференционных полос заданного периода Л во всех отмеченных спеклах. После этого выполнялось Фурье-преобразование от числового файла, соответствующего полученному графическому файлу.

Изображения, приведенные на рис. 8, убедительно показывают эффективность вышеописанного цифрового алгоритма записи и восстановления изображений. В [16,17] нами показана возможность аналогового процесса восстановления изображения по цифровой Фурье-спеклограмме. В этом случае графический файл, получаемый после компьютерной обработки цифровой Фу-рье-спеклограммы, с требуемым увеличением распечатывался, а затем, с необходимым уменьшением фотографически переносился на фотопластинку. Таким образом получали голографически подобный дифракционный оптический элемент (ДОЭ). При освещении такого ДОЭ лазерным пучком в дальней зоне дифракции наблюдались изображения записанного объекта. Пример такого восстановления изображения приведен на рис. 9.

Рис. 8. Картины дифракции дальнего поля, получаемые при численном Фурье-преобразовании цифровой Фурье-голограммы (а) и безопорных цифровых Фурье-спеклограмм с искусственно нанесенными системами несущих

интерференционных полос и соответствующие объекты-транспаранты

В экспериментах использовалась ПЗС-камера с размером матрицы 1360*1024 пикселей (6,3*4,8 мм). Графической обработке подвергались приблизительно 2000 спеклов с помощью программы обработки изображений, позволяющей в ручном режиме

Рис. 9. Картина дифракции лазерного пучка на искусственном голографически подобном дифракционном оптическом элементе

Для сравнения на рис. 10 приведен пример численного Фурье-преобразования цифровой Фурье-спеклограммы эллипса с системой несущих полос без сдвига на половину периода при переходе между соседними спеклами. Как видно из рисунка, распределения интенсивности в первых порядках дифракции не соответствуют изображениям объекта, а представляют собой автокорреляцию распределения поля по объекту, как и в нулевом порядке дифракции.

Рис. 10. Картина дифракции дальнего поля, полученная при Фурье-преобразовании спеклограммы эллипса с нанесенной системой несущих полос без сдвига на половину периода между соседними спеклами

в

г

Именно сдвиг несущих полос на половину периода при переходе от спекла к спеклу обеспечивает формирование изображений объекта в первых порядках дифракции при Фурье-преобразовании спекл-картины с системой несущих полос.

Вышеописанный способ записи и восстановления изображений применим для объектов с осесим-метричным распределением интенсивности I(х0,у0). Только в этом случае фаза дифракционного поля в различных спеклах с наибольшей вероятностью принимает значения 0 и п рад.

На рис. 11 показан результат численного восстановления изображений с использованием вышеописанного алгоритма для несимметричных объектов в форме букв Б и треугольника. В дифракционных порядках распределения интенсивностей не соответствуют изображениям объектов. Этим экспериментом подтверждается теоретическое положение о том, что для 5-коррелированного источника диф-фузно-когерентного излучения разность фаз поля в дальней области дифракции в соседних спеклах с наибольшей вероятностью равна п радиан только для центрально-симметричного объекта [15].

Рис. 11. Распределения интенсивностей, полученные при численном Фурье-преобразовании искусственных Фурье-спеклограмм несимметричных объектов: буква Ш

Особый интерес представляет случай записи спеклограммы объекта с незначительным (частичным) нарушением симметрии. На рис. 11г представлен результат восстановления такого объекта. Видно, что в целом предложенный алгоритм позволяет приближенно восстановить изображение объекта. Однако в восстановленном изображении наблюдается дополнение до симметричного вида - появилась деталь, симметрично расположенная детали, нарушающей симметрию в исходном объекте. Этот эффект пока не нашел у нас объяснения и нуждается в дальнейшем исследовании.

Численное моделирование

Процесс регистрации спекл-полей требует специального оборудования, поэтому его моделирование на компьютере существенно расширяет возможности исследователя. Спекл-структура, смоделированная на компьютере, представлена на рис. 12.

Для обработки спеклограмм с целью создания голограммноподобных структур с несущими полосами были использованы компьютерные программы обработки спеклов, позволяющие наносить на спек-лограмму систему несущих полос со сдвигом на половину периода при переходе от одного спекла к соседнему. Алгоритм восстановления фазового распределения в спекл-поле состоит в следующем:

1. 2.

3.

4.

5.

перевод изображения в массив чисел; определение среднего уровня яркости и проведение бинаризации по этому уровню (рис. 13); разделение слипшихся спеклов; светлые пятна, соединенные тонким перешейком, скорее всего, представляют разные спеклы (рис. 14); разделение спекл-картины на соседние области (рис. 15);

нанесение на спеклограмму системы несущих полос, сдвинутых на половину периода между соседними спеклами (рис. 16).

Рис. 12. Спекл-структура, смоделированная на компьютере

Рис. 13. Бинаризация изображения

Рис. 14. Разделение слипшихся спеклов

Рис. 15. Области разделения спекл-картины

Восстановленные изображения, полученные с помощью Фурье-преобразования спеклограммы с такой системой несущих полос, представлены на рис. 17. Из рисунка видно, что восстановленные изображения в некоторой степени повторяют форму источника.

а

б

в

г

Рис. 16. Спекл-картина с системой несущих полос

Рис. 17. Восстановленные изображения объектов с четным порядком центральной симметрии: креста

На рис. 17а, в представлены восстановленные изображения для случая разделения соседних областей в процессе обработки в автоматическом режиме, на рис. 17б, г - для случая разделения соседних областей в ручном режиме. Сравнение результатов применения этих режимов обработки спеклограмм показывают несколько более высокое качество восстановления изображений при использовании ручного режима определения соседних спеклов с необходимым сдвигом полос. Для повышения качества восстановления изображений автоматическая программа требует дальнейшей доработки в плане более точного определения областей, представляющих соседние спеклы. Как видно из рис. 17, в целом, предложенный алгоритм позволяет приближенно восстанавливать изображения объекта. Графическая обработка спеклограммы с применением специальных программных средств является менее трудоемким процессом по сравнению с ручной обработкой, но она сопровождается существенными ошибками при определении границ спеклов и поэтому требует дальнейшей доработки.

Заключение

Обсуждаемый в работе метод записи и восстановления изображения базируется на записи интенсивности дифракционного поля рассеивающего объекта при его когерентном освещении. Последнее обстоятельство обеспечивает спекл-модуляцию дифракционного поля, которая несет часть информации об источнике этого поля, иными словами, об изображении объекта. Теряемая при записи интен-

сивности без когерентного опорного пучка фазовая информация в рассматриваемом методе восстанавливается на основе представления о почти бинарном (0 и п рад) пространственном распределении фазы в спекл-модулированном поле в дальней области дифракции для объектов симметричной формы (апертуры). Такое восстановление фазового распределения обеспечивается созданием в записанной спекл-структуре высокочастотной квазипериодической модуляции, аналогичной модуляции несущими интерференционными полосами в голографической записи волнового поля. В итоге, получается квазипериодическая пространственная несущая - система полос, модулированная по амплитуде и фазе спекл-структурой записанного дифракционного поля; спекл-структура выступает в качестве огибающей пространственной несущей. При этом требуемая фазовая модуляция обеспечивается сдвигом несущих полос на половину их периода при переходе от одного спекла к соседнему.

Вышеописанная процедура восстановления амплитудно-фазовой структуры дифракционного поля, как показали натурные и численные эксперименты, обеспечивает восстановление изображения удовлетворительного качества, позволяющее сделать вывод о правильности применяемых подходов. Относительно высокие искажения восстановленного изображения (относительно низкое отношение сигнал/шум) по отношению к голографическому, записанному и восстановленному в тех же условиях, обусловлено, на наш взгляд, следующими недостатками обсуждаемой процедуры восстановления амплитудно-фазового распределения. В первую очередь это то, что фаза 0 и п у соседних спеклов понимается в вероятностном смысле, а не точно. Во вторую очередь следует указать на имеющую место неточность определения соседних спеклов дифракционного поля, в которых должен быть фазовый сдвиг в п радиан и, соответственно, сдвиг несущих полос на половину их периода.

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

К четвертой причине следует отнести отсутствие модуляции контраста создаваемой несущей системы полос, отвечающего за запись и восстановление амплитуды объектного поля в отдельных его спеклах. В используемой процедуре мы создавали несущую систему полос с единичным контрастом (с единичным коэффициентом модуляции), тогда как в голо-графической записи контраст несущих полос определяется соотношением интенсивностей опорной и

объектной волн в каждом спекле. Таким образом, в нашей процедуре информация о распределении амплитуды объектного поля восстанавливалась не в полной мере, что приводило, очевидно, к искажению восстановленного изображения. Учет вышеописанных факторов в процедуре восстановления амплитудно-фазовых информации дифракционного поля позволит существенным образом повысить качество (четкость) восстанавливаемых изображений.

Таким образом, в работе показана принципиальная возможность реконструкции фазы спекл-модулированного дифракционного поля симметричного плоского рассеивающего объекта. Показана также возможность восстановления изображений плоских объектов по записи безопорной Фурье-спеклограммы. Экспериментально подтверждена работоспособность предложенного метода восстановления изображения. Реализован вариант цифровой безопорной Фурье-голограммы на основе разработанной методики восстановления изображения объекта по зарегистрированной интенсивности спекл-модулированного дифракционного поля.

Результаты работы позволяют расширить представления о свойствах диффузно-когерентного поля и могут быть использованы для создания новых технологий и подходов в области оптической обработки информации. На основе проведенных исследований возможна разработка новых методов и устройств голографической интерферометрии, спекл-интерферометрии, вычитания изображений для центрально-симметричных объектов, создание специальных дифракционных элементов.

Благодарности Авторы благодарят В.В. Лычагова, к.ф.-м.н.

0.А. Перепелыцыну, к.ф.-м.н. Д.В. Лякина за помощь в экспериментальных исследованиях и Д.М. Масленникова за создание программы обработки цифровых изображений. Исследования были проведены при поддержке грантов РФФИ №06-08-00987а и CRDF REC-006.

Литература

1. Обратные задачи в оптике // Под ред. Болтса Г.П. -М.: Машиностроение, 1984. - 200 с.

2. Гончарский А.В., Попов В.В., Степанов В.В. Введение в компьютерную оптику. - М.: Издательство МГУ, 1991. - 312 с.

3. Компьютеры в оптических исследованиях // Под ред. Б. Фридена - М.: Мир, 1983. - 488 с.

4. Fernández-Guasti M., Jiménez J.L., Granados-Augustín F., Cornejo-Rodríguez A. Amplitude and phase representation of monochromatic fields in physical optics

// Journal of the Optical Society of America A, 2003. -V. 20. - P. 1629-1634.

5. Методы компьютерной оптики // Под ред. В.А. Сой-фера. - М.: Физматлит, 2000. - 688 с.

6. Kotlyar V.V., Serafimovich P.G.,Soifer V.A. Regu-larisated iterative algorithm for the phase retrieval // Optik, 1993. - V. 94. - P. 96-99.

7. Nitta K., Shogenji R., Miyatake S., Tanida J. Image reconstruction for thin observation module by bound optics by using the iterative back projection method // Applied Optics, 2006. - V. 45. - I. 13. - P. 2893-2901.

8. Wu J. S., Weierstall U., Spence J. C. H., Koch C. T. Iterative phase retrieval without support // Opt. Lett., 2004. - V. 29. - I. 23. - P. 2737-2745.

9. Bates R.H.T., Tan D.G.H. Toward reconstructing phases of inverse-scattering signals // Journal of the Optical Society of America А, 1985. - V. 2. - P. 2013-2019.

10. Вахрушева М.В., Власов Н.Г. Сведение фазовой проблемы к расчету интерферограмм сдвига // Прикладная математика и техническая физика, 2003. - Т. 2. -№ 4. - С. 3-4.

11. Vlasov N.G., SazhinA.V., Kalenkov S.G. Solution of phase problem // Laser Physics, 1996. V. 6(2). P. 401-404.

12. Bastiaans M. J., Wolf K. B. Phase reconstruction from intensity measurements in linear systems. // Journal of the Optical Society of America A, 2003. - V. 20. - P. 1046-1049.

13. Kolenovic E. Correlation between intensity and phase in monochromatic light // Journal of the Optical Society of America A, 2005. - V. 22. - I. 5. - P. 899-906.

14. Горбатенко Б.Б., Клименко И.С., Максимова Л.А., Ря-бухо В. П. Статистические свойства пространственного распределения фазы развитого спекл-поля // Письма в ЖТФ, 1992. - Т. 18. - В. 2. - С. 26-28.

15. Горбатенко Б.Б., Клименко И.С., Максимова Л.А., Рябу-хо В.П. О некоторых статистических свойствах разности фаз в развитом спекл-модулированном поле // Оптика и спектроскопия, 1995. - Т. 78. - В. 2. - С. 316-319.

16. Горбатенко Б.Б., Рябухо В.П., Максимова Л.А. Метод восстановления изображения предмета по спекл-структуре его дифракционного поля // Письма в ЖТФ, 2004. - Т. 30. - В. 17. - С. 68-75.

17. Горбатенко Б.Б., Рябухо В.П., Максимова Л.А. Статистические свойства разности фаз в спекл-модулированном поле и метод восстановления изображения предмета по спекл-структуре его дифракционного поля // Компьютерная оптика, 2004. - В. 26. - С. 48-52.

18. Горбатенко Б.Б., Рябухо В.П., Максимова Л.А. Реконструкция пространственного фазового распределения в дифракционном спекл-поле и восстановление изображения объекта по записи интенсивности // Оптика и спектроскопия, 2006. - Т. 101. - № 5. - С. 861-865.

19. Кольер Р., Беркхарт К., Лиин Л. Оптическая голография // - М.: Мир, 1973. - 688 с.

20. Шепелевич В.В. Введение в когерентную оптику и голографию // - Минск: Выш.Шк, 1985. - 144 с.

21. Франсон М. Оптика спеклов // - М.: Мир, 1980. - 171 с.