
Приборы и устройства оптического и СВЧ диапазонов
..pdf
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 = |
2π |
× x2 , |
ω y 2 = |
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].
Оптические системы способны выполнять общие линейные преобразования поступающих на вход данных. Объём данных бывает настолько большим, что существенно превышает возможности наблюдателя. Тогда для уменьшения количества данных решающую роль будет играть линейное преобразование, выделяющие нужные группа, таким образом, чтобы они могли восприниматься наблюдателем.


386
В качестве иллюстрации общих положений рассмотрим несколько примеров оптических пространственных фильтров. Если в плоскости спектра сигнала П1 стоит диафрагма о отверстием, как показано на рис. 16.11. то такая оптическая система соответствует полосовому фильтру в области верхних частот. Ширина полосы фильтра определяется шириной отверстая диафрагмы. Это явление особенно хорошо можно продемонстрировать на фильтрации пространственного спектра объекта,
имеющего вид сетки. Так в спектральной плоскости П линзы Л (см.рис. 16.10,а) приведено изображение не изменённого спектра сетки. Из-за периодичности предмета появляется набор отдельных спектральных компонент; ширина каждой ив них определяется апертурой падающей световой волны.
При помещении в спектральную плоскость узкой цели, так чтобы через неё проходил только один ряд спектральных компонент, в плоскости изображения П3 получим спектр сетки отфильтрованной этой щелью (рис. 16.12,6). Если в спектральную плоскость поместить ирисовую диафрагму и установить её так, чтобы через нее проходила только осевая Фурье компонента, то при последовательном расширении диафрагмы можно шаг за шагом проследить Фурье-синтез сетки. Если же убрать диафрагму, а вместо неё на оси линзы в фокальной
плоскости поместить маленький экран, который закрывал бы только центральный порядок, или компоненту “ нулевой частоты", то мы подучим изображение сетки с обращенным контрастом (рис. 16.12,б).


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° по отношению к эталону в зависимости от структурной сложности объекта. Отношение сигнал/шум на выходе фильтра колеблется от нескольких десятков до нескольких сотен в зависимости от структуры объекта и окружающего его фона. Неинвариантность схемы фильтрации к масштабу и ориентации объекта может быть устранена путем использования многоканального фильтра, либо осуществлением поиска по масштабу и ориентации. В последнем случае увеличивается время анализа. Применения согласованной пространственной фильтрации связаны как с формированием меры близости в виде функции взаимной корреляции при решении задач, связанных с необходимостью выполнения операций обнаружения или идентификации, изображений, так и с необходимостью измерения интенсивности функции взаимной корреляции при решении задач измерения или контроля. С помощью согласованной пространственной фильтрации могут быть решены следующие задачи:
обнаружение объектов на некотором сложном фоне; навигационная привязка летательных аппаратов по наземным ориентирам;
исследование скорости движения облаков по фотоснимкам, получаемым с метеорологических спутников;
автоматический дактилоскопический поиск (опознавание полных отпечатков пальцев); автоматическая идентификация ИК спектров веществ;
ввод буквенно-цифровой информации в ЦВМ; интегральный контроль качества печатных плат и интегральных схем;
интегральный контроль предела усталости лопаток турбин; идентификация и распознавание электрических сигналов;
оптимальная фильтрация электрических сигналов; В оптике могут быть реализованы не только корреляционные методы распозна-
вания. Весьма вероятно, что распознающий автомат будущего будет гибридной