Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Приборы и устройства оптического и СВЧ диапазонов

..pdf
Скачиваний:
57
Добавлен:
05.02.2023
Размер:
5.41 Mб
Скачать

 

 

 

 

при

x < 0

381

 

 

0

 

lim

I ®

 

¥

при

0 < x < x o

(16.7)

 

x →

0

 

0

при

x > 0

 

 

 

 

 

обладающая следующими двумя свойствами:

 

 

d ( x ) dx

= 1

(16.8,а)

− ∞

 

 

¥

 

 

f ( x ) × d ( x

- x o ) dx = f ( x o )

(16.8,б)

- ¥

 

 

Последнее свойство часто называется фильтрующим свойством и δ - функции. Найдем

¥

 

- j×w x

× x

 

 

S ( w x ) = d

( x ) × e

dx

(16.9)

 

 

- ¥

 

 

 

 

 

формально преобразование Фурье для δ - функции:

 

 

 

Применив формулу Эйлера, получим:

 

 

 

 

 

¥

¥

 

 

 

 

S ( w x ) = d ( x ) × cos( w x × x ) dx -

d ( x ) × sin( w x × x ) dx

(16.10)

- ¥

- ¥

 

 

 

 

Поскольку по определению δ (x) равна нулю на всей оси, кроме точки x = 0 , где она стремится к бесконечности в выражении (16.10) Из первого интеграла можно вынести постоянный множитель

cos(wx × x) = 1

 

x = 0 ,

 

 

 

а из второго интеграла множитель

 

 

 

 

 

 

sin(wx × x) = 0

при

 

x = 0.

 

 

Учитывая свойство (16.8,б), получим

 

 

 

 

S ( ω x ) =

e

j×w x

× x

= 1

x = 0.

(16.11)

 

 

 

Можно δ функцию представить в виде обратного преобразования Фурье от равномерной спектральной плотности

 

 

1

¥

j×w x

× x

 

 

 

1

¥

 

j×w x

× x

 

 

d ( x ) =

 

× S ( w x ) × e

dx

=

 

×

e

d w x

(16.12)

 

× p

 

 

 

× p

 

 

2

- ¥

 

 

 

2

- ¥

 

 

 

 

 

Аналогичным образом можно написать соотношения для двумерной δ - функции;

 

x = x

o

δ(x − x0

 

 

, y − y0 ) =

y = yo

 

0

382

d( x - x o , y - y o )dx dy = 1

∞ − ∞

востальной области

∫ ∫ f ( x , y) × d( x - x o , y - y o )dx dy = f ( x o , y o ) .

− ∞ − ∞

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

Рассмотрим, как осуществляется это преобразование, на примере оптической системы, изображенной на рис. 16.9.

Примем следующие допущения: в системе не происходит поглощения и отражения электромагнитной анергии; система не имеет аберраций.

Предположим, что в передней фокальной плоскости P1 линзы Л расположен транспарант с коэффициентом, прозрачности

S (x1, y1 ) = So (x1, y1 ) ×e j×θ ( x1 , y1 ) .

Обозначим когерентную монохроматическую световую волну, распрост-

раняющуюся

в направлении z

и освещающую оптическую систему

Eo exp(i ×ω × t) .

Используя положение геометрической оптики, мы можем счи-

тать, что заданное пространственное распределение коэффициента прозрачности является модулирующей функцией для световой волны Eo exp(i ×ω × t) . На выходе

такого модулятора, комплексное распределение амплитуды светового поля можно записать в следующем виде

E1 ( x1 , y1 , t ) =

E o × S o ( x1 , y1 ) × e

j×(θ ( x1

, y1 )+ω ×t )

 

(16.13)

383

Найдем распределение амплитуды светового поля на выходе оптической системы - задней фокальной плоскости линзы Л, Эта линза с апертурой D будет осуществлять суммирование дифрагированной энергии в фокальной плоскости. Для оценки результата такого суммирования для некоторого направления V необходимо найти выражение для фаз дифрагированной энергии во всех точках плоскости, перпендикулярной направлению V1 и выразить его через координату плоскости x1, y1 . При этом будем считать, что плоская волна, падающая на линзу Л, образует с её осью очень малый угол. В общем виде, уравнение плоскости, перпендикулярной направлению V имеет вид

a1 × x1 + b1 × y1 + c1 × z1 = 0 ,

(16.14,а)

где а1 , в1 , с1 , - направляющие косинусы плоскости, перпендикулярной направлению V . Тогда расстояние до точки в плоскости Р1 до этой плоскости равно

 

 

D2

= a1 × x1 + b1 × y1

 

(16.14,б)

следовательно, для фаз дифрагированной энергии можно написать

 

q = q(x1

, y1 ) +

2 × p

× D2 = q(x1 , y1 ) +

2 × p

× (a1 × x1 + b1 × y1 ) ,

(15.15)

l

l

 

 

 

 

 

где λ - длина волны света.

Далее определяется положение точки суммирования дифрагированной

энергии, определяемое значениями

Аi,

Вi,

Сi, пользуясь тем, что свет в направ-

лении V, проходящий через центр оптической система, не отклоняется от своего

направления. Если расстояние от центра линзы до этой точки равно

R , то

 

 

a1 =

x2

 

b1

=

y2

 

 

c1 =

F

,

(16.16)

 

 

R

R

 

R

 

 

 

 

 

 

 

 

 

где x2 , y2

-

координаты выходной плоскости.

 

 

 

Обозначим через L оптическую длину пути от рассматриваемой точки до вход-

ной плоскости Рi с координатами

x1

и y1

до точки суммирования в выходной

плоскости

P2

с координатами

 

x2

и

y2 .

Подставив соотношение (16.16) в вы-

ражение (16.15), получим

 

 

 

 

 

 

 

 

 

θ = θ (x , y ) +

2 ×π

× L +

2 ×π × c1

× (x × x

 

+ y × y

 

)

(16.17)

λ

 

2

2

1

1

 

F

1

1

 

 

 

 

 

 

 

 

 

 

Теперь можно написать выражение для распределения освещенности

Е2(x2 , y2 ,

t)в плоскости P2 в виде двойного интеграла

E ( x 2 , y 2 , t) = E o e

j×w×t

∫∫ So

j{q( x1

, y1 )+

2×p

L +

2 pc1

( x1x 2 + y1y 2 )

 

l

 

 

(x1 , y1 )e

 

 

lF

dx 1dy 1

 

 

 

 

 

 

 

 

P2

 

 

 

 

 

 

 

(16.18)

384

Введем вспомогательные обозначения

ω x 2 =

× x2 ,

ω y 2 =

× y

2 ,

 

 

 

 

 

λ × F

 

λ × F

 

μ −1 .

которые носят название пространственных частот и имеют размерность

Они являются пространственными аналогами частоты изменения сигналов во

времени, имеющей размерность

c −1 .

Если

дифрагированные углы

малы, то

c1 ≈ 1 . Можно считать, что для всех значений

x2 и y2 , величина α

будет по-

стоянной. Поэтому член, содержащий

α

в выражении (16.18), можно опустить,

так как он характеризует постоянный сдвиг по всей плоскости P2. Так как обра-

ботки сигналов осуществляются

за

много периодов светового

излучения, то

Eo exp(i ×ω × t) можно заменить его среднеквадратичным значением

Eo

и принять

Еo = 1.

 

 

 

 

 

 

Тогда на основе принятых допущений окончательно получим:

 

 

E ( x2 , y2 ) = ∫∫ So ( x1 , y1 ) × e j ×x 2 ×x1 +ω y 2 × y1 +θ ( x1 , x 2 )) dx1dy1

(16.19)

P2

S(x1, y1) .

Это выражение является двумерным преобразованием Фурье функции

Если осуществлять только пространственную амплитудную модуляцию (напри-

мер, в плоскости Рi задать плоское изображение), то

θ (x1,x2) = const и выра-

жение (16.19) упростится

j ×x ×x +ω y

 

 

 

E x y ) = S x y ) = ∫∫ So ( x, y) × e

× y )

dxdy

(16.20)

 

 

 

 

 

 

Итак, в когерентной идеальной оптической системе комплексное распределение амплитуды освещенности светового поля в передней фокальной плоскости линзы связано с комплексным распределением светового поля в задней фокальной плоскости линзы двумерный преобразованием Фурье. Функция S (ω x y )

называется пространственным спектром объекта S(x,y). Плоскость, в которой дифракционная картина объекта S(x,y) описывается его преобразованием Фурье называется Фурье-плоскостью.

16.2.3. Пространственная фильтрация

Простота осуществления операций прямого и обратного фурьепреобразований над двумерными функциями (в общем случае комплексными) в сочетании с голографическим способом регистрации световых полей открывает широкие возможности для синтеза схем пространственной фильтрации с произвольным видом передаточной функции [47-48].

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

385

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

Ax1 y1 ) , то тогда

когерентный световой поток,

пройдя через транспарант с

функцией пропускания S1(x1,y1)

создаст, благодаря дифракции в фокальной

плоскости линзы Л1,

пространственный спектр объекта S x1 y1 ) . Свет, прой-

дя через линзу Л2 , образует изображение

S(x2,y2)

в плоскости

П2 с координа-

тами x2 и y2 так что изображение

S2(x2 ,y2 ) можно рассматривать как результат

операции интегрирования, происходящей над полем

S x1 y1 )

в зрачке линзы

Л2, т.е. изображение

S2(x2,y2) соответствует обратному преобразованию Фурье.

Теперь введем фильтр с пропусканием

Ax1 y1 ) на который проектируется

спектр сигнала S x1 y1 ) . На выходе плоскости П2

имеем напряженность элек-

трического поля

 

 

 

 

 

S '2 (x1 , y1 ) = Ax1 y1 ) × S x1 y1 ) ,

(16.21)

изменение которого по пространственным частотам изменено в соответствии с характеристикой фильтра.

386

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

имеющего вид сетки. Так в спектральной плоскости П линзы Л (см.рис. 16.10,а) приведено изображение не изменённого спектра сетки. Из-за периодичности предмета появляется набор отдельных спектральных компонент; ширина каждой ив них определяется апертурой падающей световой волны.

При помещении в спектральную плоскость узкой цели, так чтобы через неё проходил только один ряд спектральных компонент, в плоскости изображения П3 получим спектр сетки отфильтрованной этой щелью (рис. 16.12,6). Если в спектральную плоскость поместить ирисовую диафрагму и установить её так, чтобы через нее проходила только осевая Фурье компонента, то при последовательном расширении диафрагмы можно шаг за шагом проследить Фурье-синтез сетки. Если же убрать диафрагму, а вместо неё на оси линзы в фокальной

плоскости поместить маленький экран, который закрывал бы только центральный порядок, или компоненту “ нулевой частоты", то мы подучим изображение сетки с обращенным контрастом (рис. 16.12,б).

387

16.2.4. Согласованная фильтрация

Для решения различных задач обнаружения сигнала в шумах наибольшее распространение получил метод оптической согласованной фильтрации [48]. Достоинством оптической согласованной фильтрации (OCФ) является высокая информационная емкость, банк эталонных сигналов, многоканальность, обработка больших массивов данных, за время, равное времени распространения света в оптической системе от модулятора до фотоприемника. Наряду с достоинствами этот метод имеет и недостатки: необходимости изготовления голографических согласованных фильтров (ГСФ). Юстировка этих фильтров в частной плоскости оптической системы должна проводиться с очень высокой точностью, которая предъявляется также и всей оптической системе.

Создание согласованного голографического фильтра (СГФ)

Рассмотрим вопрос о записи согласованных голографических фильтров. Если полезный известной формы обозначить через S(x) , а помеху через n(x) , то оптический сигнал на входном транспаранте когерентной оптической системы представляется в виде суммы

f ( x) = S ( x) + n( x)

(16.22)

спектральной плоскости спектр f(x) будет состоять из спектров сигнала и шума

F x 2 ) = S x 2 ) + N x 2 )

Если в спектральной плоскости поставить пространственный фильтр с характеристикой пропускания

H x 2 ) = S * x2 ) ,

(16.23)

то суммированный cигнал на выходе фильтра будет равен

 

Sx 2 ) × S* x 2 ) = N2 x2 ) × S* x2 )

(16.24)

Если теперь произвести над этим сигналом обратное преобразование Фурьe, то в фокальной плоскости второй линзы образуется автокорреляционная функция сигнала и корреляционная функция сигнала и шума

 

E3 (x3 , y3 ) = S(x1 ) × S* (x1 - x3 )dx + n(x1 ) × S* (x1 - x3 )dx1

(16.25)

−∞

−∞

 

Поскольку сигнал и шум не коррелированы, то второе слагаемое будет мало по сравнению с первым и сигнал будет выделен из шума.

388

Рассмотрим, как создается комплексный оптимальный фильтр оптическими методами. В частности на принципах голографии (см. рис. 16.13). Эталонный S(x) сигнал устанавливается в виде транспаранта в предметной плоскости P1. Линза 5 осуществляет прямое преобразование Фурье этого сигнала. С помощью призмы 3 и поворотного полупрозрачного зеркала формируются опорный и сигнальный лучи, которые пересекаются в светочувствительной среде под углом θ .Запись голограммы Фурье-образа с помощью лазера (1). Равномерная засветка транспаранта осуществляется и помощью коллиматора 4. Светочувствительная среда зарегистрирует интенсивность интерференционной картины

I (ω) = [S(ω) + e j×ω×xo ]2 = 1 + S 2 + S * × e j×ω×xo + S × e- j×ω×xo

(16.26)

Далее, устанавливая голограмму фильтра в оптическую систему, осуществляют процесс распознания образов.

Фильтрация в оптическом корреляторе. Когерентные оптические кор-

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

389

опустим теперь через согласованный фильтр (6) световой сигнал, имеющий образ входного сигнала f(х ), т.е. F(ω ω ω ) рис. 16.14. В плоскости P1 устанавливается транспарант анализируемого сигнала. В плоскости P2 - СГФ. Линза Л2 осуществляет обратное преобразование Фурье. На выходе фильтра, возникает световой сигнал с амплитудой

H (ω ) = F (ω ) × I (ω ) = F [1 + S 2 ]+ F × S * × e j×ω×xo + F × S × e- j×ω×xo

(16.27)

Этот сигнал имеет три составляющих. Основной пучок, отображаемый первым членом в выражении (16.27), распространяется в направлении оптической оси. Второй пучок, описываемый членом F × S * ×exp( j ×ω × x0 ) идет в направлении

опорного пучка. Третий член F × S × exp( - j ×ω × x0 ) описывает пучок, отклоняю-

щийся относительно оси в противоположном направлении. Сигнал, передаваемый опорным пучком

F × S * = (S + N ) × S* = S 2 + N × S *

(16.28)

состоит из двух слагаемых. В результате преобразования Фурье, функция

S2

превратится в сигнал S × S * , т.е. в сигнал автокорреляции функции

S(x) . Этот

сигнал появится в центре верхнего бокового изображения, в виде яркой точки. Первое слагаемое нижнего бокового пучка

F × S = (S * + N) × S = S × S * + N × S

(16.29)

превратится в сигнал автосвертки S Ä S . Сигналы N × S * и N × S

вследствие

случайной природы шума и определенного рассогласования в фильтре, изготовленном для обнаружения только полезного сигнала S(x) значения свертки nÄ S кросс корреляции N × S * будут малы.

390

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

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

Согласованный пространственный фильтр осуществляет параллельный корреляционный анализ входного изображения. Результат анализа может быть представлен в виде световой картины, т. е. поля функций корреляции входного и эталонного изображений. Схема согласованной пространственной фильтрации не чувствительна к положению опознаваемого объекта в пределах анализируемого изображения, однако допустимые изменения его масштаба и угловой ориентации не должны превышать соответственно ±(5÷10)% и 2÷5° по отношению к эталону в зависимости от структурной сложности объекта. Отношение сигнал/шум на выходе фильтра колеблется от нескольких десятков до нескольких сотен в зависимости от структуры объекта и окружающего его фона. Неинвариантность схемы фильтрации к масштабу и ориентации объекта может быть устранена путем использования многоканального фильтра, либо осуществлением поиска по масштабу и ориентации. В последнем случае увеличивается время анализа. Применения согласованной пространственной фильтрации связаны как с формированием меры близости в виде функции взаимной корреляции при решении задач, связанных с необходимостью выполнения операций обнаружения или идентификации, изображений, так и с необходимостью измерения интенсивности функции взаимной корреляции при решении задач измерения или контроля. С помощью согласованной пространственной фильтрации могут быть решены следующие задачи:

обнаружение объектов на некотором сложном фоне; навигационная привязка летательных аппаратов по наземным ориентирам;

исследование скорости движения облаков по фотоснимкам, получаемым с метеорологических спутников;

автоматический дактилоскопический поиск (опознавание полных отпечатков пальцев); автоматическая идентификация ИК спектров веществ;

ввод буквенно-цифровой информации в ЦВМ; интегральный контроль качества печатных плат и интегральных схем;

интегральный контроль предела усталости лопаток турбин; идентификация и распознавание электрических сигналов;

оптимальная фильтрация электрических сигналов; В оптике могут быть реализованы не только корреляционные методы распозна-

вания. Весьма вероятно, что распознающий автомат будущего будет гибридной