CC BY
22
РАДИОФИЗИКА, ЭЛЕКТРОНИКА, АКУСТИКА
Компенсация хроматических аберраций в акустооптических системах
спектрального анализа изображений
В. Б. Волошинов1'0, К. Б. Юшков2'6, Т. В. Юхневич1
1 Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра физики колебаний. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.
2Национальный исследовательский технологический университет «МИСиС», Научно-технологический и учебный центр акустооптики. Россия, 119049, Москва. E-mail: а [email protected], b [email protected] Статья поступила 18.11.2011, подписана в печать 29.05.2012.
В работе рассмотрены усовершенствованные оптические схемы акустооптических систем спектрального анализа изображений, характеризующиеся низким уровнем продольных и поперечных хроматических аберраций. Рассмотрены физические механизмы компенсации хроматизма акустооп-тической дифракции за счет оптимального использования дисперсионных свойств пассивных оптических элементов системы. Снижение продольных хроматических аберраций рассмотрено для системы с одним объективом и дополнительной рассеивающей линзой, а также в конфокальной системе, сформированной двумя объективами. Экспериментальное исследование в видимом диапазоне спектра показало снижение продольных смещений в 2.5 раза и поперечных смещений на два порядка и более.
Ключевые слова: акустооптика, акустооптические фильтры, обработка изображений, продольные и поперечные смещения изображений.
УДК: 534.23:535.42:535.243.3:681.7.012.2. PACS: 42.15.Fr; 42.25.Fx; 42.30.Va; 42.79.Jq.
Введение
Известно, что брэгговская дифракция света в анизотропных средах используется для управляемой спектральной фильтрации световых пучков [1-4]. Перестраиваемые спектральные фильтры являются одним из основных классов акустооптических (АО) устройств, находящих применение при создании уникальных приборов и систем. Неколлинеарные АО фильтры, в которых свет распространяется под небольшим углом к акустическому волновому фронту, позволяют обрабатывать неколлимированные световые пучки [1-5]. Таким образом, на основе широкоапертурных АО фильтров становится возможным создание спектральных систем для анализа изображений с электронной перестройкой и возможностью синтеза многополосной функции пропускания, с рабочим диапазоном, превышающим октаву, и характерным быстродействием менее миллисекунды [5-15].
При перестройке длины волны света, пропускаемого АО фильтром, в обрабатываемых изображениях возникают хроматические аберрации, которые можно разделить на продольные и поперечные [4]. Поперечные аберрации вызваны зависимостью угла дифракции от длины волны света [6-9]. Общая схема АО системы обработки изображений и возникающие в ней хроматические аберрации представлены на рис. 1. Влияние поперечных аберраций заключается в боковом смещении изображения в плоскости дифракции на угол 7', как показано на рисунке. Продольные хроматические аберрации обусловлены изменением оптической длины пути в системе из-за дисперсии показателей преломления АО кристалла [7]. В результате плоскость наилучшей резкости перемещается вдоль оптической оси системы
Рис. 1. Принципиальная оптическая схема системы обработки изображений с АО фильтром: / — объект с широкополосным спектром излучения; 2 — объектив; 3 — поляризатор; 4 — перестраиваемый АО фильтр; 5 — скрещенный поляризатор; 6 — изображение в первом дифракционном порядке на длине волны А]; 7 — изображение в первом дифракционном порядке на длине волны Аг (А] <Аг). На вставке схема хода лучей в акустооптическом фильтре
на расстояние 5Ь, что является причиной дефокусировки изображения. Негативное влияние хроматических аберраций на характеристики АО систем обработки изображений особенно возрастает для фильтров со сверхшироким рабочим диапазоном, имеющим два пье-зопреобразователя [12]. В подобных системах диапазон перестройки длин волн света превышает две октавы, что неизбежно приводит к поперечным и продольным смещениям отфильтрованного изображения. Компенсация аберраций необходима и в устройствах, работающих в спектральном диапазоне, близком к границе про-
зрачности кристалла [13, 14]. Например, в кристалле парателлурита, выбранном для исследования в настоящей работе, аберрации наиболее сильны в ультрафиолетовом диапазоне на длинах волн А > 350 нм [15].
В настоящей работе рассмотрены некоторые методы, позволяющие компенсировать продольные и поперечные хроматические аберрации в широком спектральном диапазоне. Общей чертой использованных методов является компенсация хроматических аберраций, вызванных акустооптической дифракцией, за счет дисперсии оптических элементов, приводящих к аберрациям противоположного знака. Рассмотрены особенности конструкции АО ячеек, позволяющие заметно снизить поперечные аберрации, а также оптические схемы, позволяющие уменьшить продольные или поперечные аберрации. Численные оценки в настоящей работе сделаны для АО фильтра на основе парателлурита с рабочим диапазоном длин волн от Аш1п = 400 нм ДО Атах = 1000 нм и с центральной длиной волны Ао = 633 нм. Падающая на кристалл волна является необыкновенной, поляризация которой лежит в плоскости АО взаимодействия в кристалле.
1. Снижение поперечной аберрации
Для снижения поперечной хроматической аберрации в АО фильтрах наиболее часто используется метод, при котором выходная грань фильтра наклоняется по отношению к входной грани на некоторый угол [6-9]. Физический принцип данного метода заключается в том, что дисперсия угла дифракции и угла рефракции света на выходной грани могут иметь противоположные знаки и, следовательно, частично компенсировать друг друга. Входная грань кристалла выполняется ортогонально падающему пучку для снижения оптических потерь на отражение и облегчения юстировки прибора. Известно, что при широкоугольной геометрии дифракции лучевые векторы нулевого и первого порядков параллельны. Поэтому угол отклонения волнового вектора дифрагировавшего луча 7 зависит от угла падения света в¡, отсчитываемого от кристаллографической оси по формуле [10]
1 = в^агс\ё[(пЦп2е)\ёв/],
(1)
где п0 и пе — главные показатели преломления кристалла. Если плоскость выходной оптической грани развернута относительно входной на угол ¡3, то после преломления на выходе из фильтра дифрагировавшая волна распространяется под углом
7' = /3 + агсэт [п0 эт(7 - /3)]
(2)
матической аберрации угол наклона выходной грани ¡3[ определяется выражением
п0 (Агпт М^тш) - «о(А тах Ы-^тах)
А
(3)
«л(Атт) ^о(Атах)
в котором угол падения света в-, выступает в качестве параметра, определяющего геометрию АО взаимодействия. При разработке фильтров также часто используется конфигурация кристалла, при которой дифрагировавший луч распространяется параллельно падающему. Условие 7' = 0 позволяет получить уравнение для величины угла наклона выходной грани кристалла эт/З = п0 эт(/3 — 7), откуда в приближении малых углов величина ¡3 определяется выражением
А> » П,л/(п0 - 1). (4)
Сравнение расчетных спектральных изменений угла распространения света при дифракции приведено на рис. 2 для выражений (3) и (4), а также для АО фильтра без наклона выходной грани. Для анализа был выбран широкоапертурный фильтр на кристалле парателлурита с плоскостью акустооптического взаимодействия (110) и медленной сдвиговой акустической волной, распространяющейся под углом а = 10° к оси [110] [15]. При вычислении величина угла падения света была выбрана равной в-, = 23.4°, что соответствует углу а = 10° между направлением волнового фронта ультразвука и осью 2. При этом на центральной длине волны Ао = 633 нм угол отклонения волнового вектора обыкновенной волны от оси 2 равен 7 = 2.6°. Соответствующие значения углов наклона выходной грани равны ¡3\ = 6.6° и Д> = 4.8°. Из представленных данных видно, что для данной геометрии кристалла выбор оптимального угла наклона выходной грани (кривая 1) приводит к снижению углового отклонения луча на 2 порядка по сравнению с плоскопараллельным исполнением кристалла (кривая 3). Если наклон выходной грани соответствует преломлению дифрагировавшего луча в направлении, параллельном падающему лучу (кривая 2), то снижение угла отклонения происходит приблизительно в 3 раза относительно конфигурации
к направлению распространения падающей волны. Именно этот угол определяет поперечные смещения отфильтрованного изображения.
Оптимальная величина угла ¡3 определяется численными методами на основании уравнений (1) и (2) с учетом дисперсии показателей преломления кристалла п0{А) и пе{А). Критерием выбора значения ¡3 является достижение минимального по абсолютному значению сканирования дифрагировавшего луча в диапазоне рабочих длин волн фильтра от Аш1п до Атак. Для наиболее эффективного снижения поперечной хро-
0.6 0.7 0.8 Длина волны X, мкм
Рис. 2. Поперечное смещение луча при дифракции немонохроматического излучения: / — при максимальной компенсации поперечной аберрации {[5 = [5\); 2 — при сохранении направления распространения луча (¡3 = /%); 3 — для плоскопараллельного кристалла без компенсации (,0 = 0). На вставке — кривая / в увеличенном масштабе
с углом наклона (3 = 0. С другой стороны — угол отклонения лучей при сохранении направления света оказывается больше минимально достижимого, реализуемого при ¡3 = /3[.
Для компенсации поперечных аберраций АО фильтров может быть использована конфокальная оптическая схема [11]. В такой системе используются две линзы, формирующие изображение объекта на фотоприемной матрице. Акустооптическая ячейка размещается в плоскости изображения, формируемого первой линзой, и эта же плоскость является объектной плоскостью второй линзы. Для разделения нулевого и первого порядков дифракции и в такой системе необходимо использовать телецентрические диафрагмы. Более подробно конфокальная система рассматривается в следующем разделе, при этом основное внимание уделяется возникающим в ней продольным хроматическим аберрациям.
2. Оптические схемы с компенсацией продольной аберрации
Продольная хроматическая аберрация существенным образом зависит от оптической схемы, используемой для формирования изображения. В качестве линз в системах формирования изображений, как правило, используются анахроматические объективы. Также на практике вместо линз часто используется зеркальная оптика, полностью лишенная собственных хроматических аберраций. Поэтому в настоящей работе при расчетах параметров системы собственные хроматические аберрации используемой оптики не учитывались.
Задержка оптического пути в кристалле равняется
Д/ = (л0-1)1,, (5)
где Ьс — длина кристалла. При проведении расчетов был выбран показатель преломления обыкновенной волны, так как при малых углах падения света, наиболее часто реализуемых в широкоугольной геометрии дифракции, п,'йл,„ где щ — показатель преломления падающего на кристалл света. Важно, что при сделанных предположениях величина продольной аберрации не зависит от фокусного расстояния объектива Анализ показал, что в простейшей системе, модель которой представлена на рис. 1, продольное смещение изображения при перестройке длины волны определяется выражением (5). Для количественного описания дисперсии показателей преломления АО материала удобно использовать величину относительной дисперсии [4]
п0(А) - по(\0)
d{ Х) =
(6)
п0(А0) - 1
Таким образом, продольное смещение или продольная аберрация, обусловленная дисперсией кристалла, определяется выражением [7]
дп0 1
6ЬС{\) ■■
ЭХ ril
I,(A-Ao).
(7)
Величина продольной аберрация в данной схеме может быть существенно снижена, если в оптической схеме устройства дополнительно использовать рассеивающую линзу, не имеющую компенсации собственной хроматической аберрации. Модель такой системы представлена на рис. 3, а. Анализ доказывает, что
L-TU-
2 3
Рис. 3. Усовершенствованные оптические схемы систем спектрального анализа изображений со снижением продольной хроматической аберрации с дополнительной рассеивающей линзой (а) и двухлинзовая конфокальная система (б): / — объект; 2 — рассеивающая линза; 3 — объектив; 4 — поляризатор; 5 — перестраиваемый АО фильтр; 6 — скрещенный поляризатор; 7 — дополнительная линза; 8 — изображение в первом дифракционном порядке
рассеивающая линза 2 на рис. 3 с показателем преломления щ{А) и фокусным расстоянием /7(А) = 0.5Й/Х х(я/(А) — I)-1 при одинаковых радиусах кривизны поверхностей й/ < 0 вносит отрицательную продольную хроматическую аберрацию
fliF0
^'-"'W -(..:<,-ww «<*■>■
(8)
которая зависит от расстояния между объектом и рассеивающей линзой fli>0, а также от расстояния между рассеивающей линзой и объективом 1\ >0, где di(X) = (rt/(А) — rt/(A0))/(rt/(A0) — 1). Таким образом, выбор двух параметров а^ и Ii определяет величину 5bi{\m\n) — öb/(\max) <0, что позволяет скомпенсировать положительное смещение 3bc(Xmin) - 3bc(Xmax), определяемое уравнением (7). Необходимо отметить, что на величины ci\ и 1\ накладывается следующее ограничение, обусловленное требованием формирования объективом действительного изображения на конечном расстоянии:
Ь-Т^гЖ- (9)
а 1 — -п
Таким образом, полная величина продольной хроматической аберрации в рассматриваемой системе определяется выражением
5Ь = 5ЪС + Sbi.
(10)
При определении оптимальных параметров оптической системы с учетом аберраций, рассчитанных по формулам (7) и (8), эффективное фокусное расстояние основного объектива считалось постоянным, ^о = 8.5 см. Также учитывалось, что фокусное расстояние рассеивающей линзы на центральной длине волны составляло величину /-}(Ао) = ^35 см, равную
фокусному расстоянию линзы, использованной в эксперименте. При анализе в качестве параметров оптимизации были выбраны расстояния между оптическими элементами ct\ и 1\.
В системах спектрального анализа изображений наряду с простой оптической системой, изображенной на рис. 1 и сформированной единственным объективом 2, также зачастую применяются конфокальные системы, основанные на двух объективах [11]. К преимуществам таких систем относятся гибкость характеристик (возможность изменять коэффициент линейного увеличения изображения при постоянных фокусных расстояниях линз) и возможность компенсировать продольную хроматическую аберрацию АО фильтра [8]. Оптическая схема конфокальной системы приведена на рис. 3, б, в которой по сравнению со схемой на рис. 3, а использована фокусирующая линза 8 и исключена рассеивающая линза 2. Кроме того, возможно использование афокальной двухлинзовой системы. В этой системе объект расположен в передней фокальной плоскости первого объектива, а изображение формируется в задней фокальной плоскости выходного объектива. Такая система лишена продольной хроматической аберрации, однако поперечные аберрации проявляются в ней наиболее сильно. Другим недостатком афокальной системы является фиксированный коэффициент линейного увеличения изображения, что ограничивает возможности согласования поля зрения прибора с размерами светочувствительной матрицы.
Ниже подробно рассматривается продольная аберрация в конфокальной системе. Считается, что система сформирована двумя одинаковыми объективами с фокусным расстоянием Fq, расстояние между смежными главными плоскостями которых составляет величину Ц). Объект расположен на расстоянии от передней главной плоскости первого объектива, а изображение при центральной длине волны А0 сфокусировано в плоскости на расстоянии Ь0 от задней главной плоскости второго объектива. Изменение оптического пути между линзами за счет дисперсии показателей преломления кристалла приводит как к изменению положения главных плоскостей системы в целом, так и к изменению ее эффективного фокусного расстояния. В результате смещение положения плоскости изображения может быть найдено как функция длины волны света:
Ви(А) _ aF0dAli3aV - 2/7о) + 4F0 + adAl} (11)
0 [a(l - 2F0) - F|]2(/ + dAl) '
где приняты обозначения / = /о + AI и а = а о -- IF0/(1 - 2F0), а относительная дисперсия d = d{А) определяется выражением (6). Величины AI и F0 конструктивно определены выбором линз и размеров акустооптической ячейки, в то время как расстояния Щ), /о и Ь{) могут меняться при юстировке системы.
Выражение (11) получено в приближении малой дисперсии d{А) <с 1 при Amin < А < Ашах. При определенных параметрах оптической системы зависимость 5Ь0(Х) является немонотонной, что позволяет провести эффективную компенсацию продольной хроматической аберрации в широком спектральном диапазоне. Проведенный анализ доказывает, что наилучшая ком-
пенсация аберрации достигается в том случае, если выполнено условие &o(Amin) = MAmax) ■ Это условие означает, что величине относительной дисперсии при перестройке длины волны света в диапазоне длин волн от Amin до Атах соответствует нулевой сдвиг изображения 5Ьо = 0. Из последнего условия было получено уравнение, связывающее параметры оптической схемы а и I, с величиной относительной дисперсии d в качестве параметра:
4/f
3(1 - 2F0) + dAl'
(12)
Оказалось, что выбор значений а и I в соответствии с уравнением (12) оставляет возможность свободного изменения коэффициента увеличения системы.
3. Экспериментальное исследование продольного смещения изображения на выходе фильтра
Измерение продольной хроматической аберрации проводилось как для простой однолинзовой системы с дополнительной рассеивающей линзой, показанной на рис. 3, а, так и для конфокальной системы, сформированной двумя объективами и представленными на рис. 3, б. Измерения проводились в видимом диапазоне спектра на длинах волн от А1 = 450 нм до Аг = 650 нм.
Для измерения продольной хроматической аберрации был использован АО фильтр на основе парател-лурита с углом наклона волнового вектора ультразвука а = 10° и длиной кристалла 1С = 2.7 см [15]. При данной конфигурации кристалла, имеющего при А0 = 633 нм главные показатели преломления п0 = 2.26 и пе = 2.41, показатель преломления для необыкновенной волны в направлении = 23.6° составлял величину /г; = 2.28. В качестве основных объективов были использованы два одинаковых фотографических объекти-ва-анахромата с эффективным фокусным расстоянием
= 8.5 см, а также рассеивающая линза из плавленого кварца с фокусным расстоянием ^(Ао) = ^35 см. В качестве источника света с непрерывным спектром применялась лампа накаливания, а отфильтрованное изображение регистрировалось с помощью приемной камеры или фиксировалось визуально с помощью микроскопа.
Основные результаты экспериментального исследования представлены на рис. 4. На рис. 4, а представлено сравнение величины продольной хроматической аберрации в системе с одним анахроматичесим объективом, полученные без компенсации смещения (кривая 1) и с компенсацией при помощи дополнительной рассеивающей линзы (кривая 2). Сплошными линиями на рисунке показаны результаты расчетов. Из рисунка следует, что в интервале длин волн от А1 = 450 нм до А2 = 650 нм отфильтрованное изображение перемещается на выходе фильтра без компенсации аберраций на расстояние 5Ь0 около 0.1 см. Применение рассеивающей линзы привело к снижению продольных смещений до величины 5&о~0.03 см. Таким образом, применение рассеивающей линзы заметно уменьшило нежелательный эффект продольного смещения изображения на выходе фильтра. Наблюдаемое на рис. 4, а незначительное расхождение наблюдаемых данных с расчетными объясняется тем, что при расчетах не могли быть
1.0
0.5
| 0
0
g -0.5
g 0.4 0.6 0.8 1.0
1 Длина волны X, мкм
о
-1 -2
0.4 0.6 0.8 1.0
Длина волны А,, мкм
Рис. 4. Измерение продольного смещения изображения в видимой области спектра в системе с дополнительной рассеивающей линзой (а) и в двухлинзовой конфокальной системе (б). На рисунках: / — без компенсации; 2 — с компенсацией
учтены остаточные хроматические аберрации использованных в экспериментах объективов.
Результаты измерения продольной аберрации в конфокальной системе представлены на рис. 4,6. Кривая 1 соответствует оптической системе фильтрации без компенсации продольной аберрации. Приведенные на рис. 4,6 данные доказывают, что после проведения компенсации по методике, описываемой формулой (11), величина продольной аберрации и видимом диапазоне была снижена более чем на 40% (кривая 2). При этом также в 1.5 раза возросла глубина резкости изображения, что удобно с практической точки зрения при создании систем обработки изображений. Сплошные линии на рис. 4,6 соответствуют расчетному значению вторичного спектра аберрации, полученному в соответствии с выражением (11). Как и в случае применения рассеивающей линзы, компенсация смещений в конфокальной схеме значительно улучшает рабочие характеристика системы фильтрации изображений. С другой стороны, исследование также доказывает, что компенсация продольных смещений изображений
становится проблематичной в том случае, если длина волны оптического излучения близка к краю полосы пропускания кристалла парателлурита при А < 450 нм.
Заключение
В работе были теоретически исследованы и экспериментально проверены методы снижения хроматических аберраций в АО системах спектрального анализа изображений. В качестве основного метода снижения поперечной аберрации использован наклон выходной оптической грани кристалла. Для снижения продольной хроматической аберрации рассмотрено использование дополнительной длиннофокусной рассеивающей линзы в системе с одним основным объективом, а также предложен метод компенсации продольной хроматической аберрации в конфокальной оптической системе, сформированной двумя анахроматическими объективами. Измерение продольных аберраций в системах с компенсацией показало существенное снижение влияния аберрации на качество отфильтрованного изображения.
Работа выполнена при финансовой поддержке РФФИ (гранты 10-07-00683 и 10-07-000623), а также Программы создания и развития НИТУ «МИСиС».
Список литературы
1. Балакший В.И., Парыгин В.Н., Чирков Л.Е. Физические основы акустооптики. М., 1985.
2. Хи L, Stroud R. Acousto-optic Devices. N. Y., 1992.
3. Goutzoulis A., Pape D. Design and Fabrication of Acousto-optic Devices. N.Y., 1994.
4. Борн M., Вольф Э. Основы оптики. M., 1970.
5. Changl.C. // Appl. Phys. Lett. 1974. 25, N 7. P. 370.
6. Yano T., WatanabeA. // Appl. Opt. 1976. 15, N 9. P. 2250.
7. Волошинов В.Б., Кулаков Л.А., Миронов О.В. // Радио-техн. и электрон. 1988. 33, № 10. С. 2177.
8. Епихин В.М., Калинников ¡O.K. // ЖТФ. 1989. 59, № 2. С. 160.
9. Voloshinov V.B., Yushkov К.В., Linde В. // J. Opt. A: Pure and Appl. Opt. 2007. 9, N 4. P. 341.
10. Волошинов В.Б., Москера X.C. // Опт. спектр. 2006. 101, № 4. С. 675.
11. Suhre D.R., De nés L.J., Gupta N. // Appl. Opt. 2004. 43, N 6. P. 1255.
12. Gupta N., Voloshinov V.B. // Appl. Opt. 2007. 46, N 7. P. 1081.
13. Gupta N., Voloshinov F. // Opt. Lett. 2005. 30, N 9. P. 985.
14. Gupta N., Voloshinov V.B. // Appl. Opt. 2004. 43, N 13. P. 2752.
15. Волошинов В.Б., Юшков К.Б. // Радиотехн. и электрон. 2007. 52, № 6. С. 727.
Chromatic aberrations compensation in acousto-optic systems of spectral analysis of images V.B. Voloshinov1,0, K.B. Yushkov2ft, T.V. Yukhnevich1
1 Department of Oscillation Physics, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.
2 Acousto-Optical Research Center, National University of Science and Technology «MISIS», Moscow 119049, Russia.
E-mail: a [email protected], b [email protected].
The paper considers novel optical schemes related to acousto-optic systems of spectral analysis of images. The systems are characterized by a low level of transverse and longitudinal chromatic aberrations. Mathematical expressions describing magnitudes of the transverse aberration in various configurations of acousto-optic filters are derived. A decrease of the longitudinal chromatic aberrations is examined in an optical scheme based on a single objective lens and an additional negative lens. A similar analysis is carried out for a confocal system formed by two objective lenses. The experimental investigations carried out in the visible range of spectrum demonstrated a decrease in the longitudinal shifts of images to a factor of 2.5 and in the transverse image shifts to as much as one hundred of magnitude.
Keywords: acousto-optics, acousto-optic tunable filters, processing of images, longitudinal and transverse shifts of images.
PACS: 42.15.Fr; 42.25.Fx; 42.30.Va; 42.79.Jq.
Received 28 November 2011.
English version: Moscow University Physics Bulletin 5(2012).
Сведения об авторах
1. Волошинов Виталий Борисович — канд. физ.-мат. наук, доцент; тел.: (495) 939-44-04, e-mail: [email protected].
2. Юшков Константин Борисович — канд. физ.-мат. наук, ст. науч. сотрудник; тел.: (495) 959-99-84, e-mail: [email protected].
3. Юхневич Татьяна Викторовна — аспирантка; тел.: (495) 939-44-04, e-mail: [email protected].
