Из условия (5) следует, что для оптимальной фильтрации полезного сигнала, осуществляемой в плоскости изображений, необходимо такое пропускание по полю пространственного фильтра-растра, чтобы оно соответствовало закону изменения освещенности в изображении объекта. Например, для обнаружения точечного излучателя необходимо в плоскости анализа изображения установить полевую диафрагму очень малых размеров с пропусканием h(a), соответствующим распределению освещенности s(a) в изображении точки. Поскольку в большинстве случаев вид функции s(a) либо трудно определить, либо он меняется для различных условий работы прибора, обычно применяют диафрагму малых размеров с резким переходом от прозрачной части, по форме повторяющей контур изображения, к непрозрачной.
Такая конструкция пространственного фильтра практически очень неудобна. Во-первых, для просмотра большого поля обзора малым мгновенным угловым полем (диафрагмой анализатора) затрачивается сравнительно много времени. Во-вторых, как уже отмечалось выше, для исключения потери информации при переходе от двумерного (пространственного) представления сигнала к одномерному (временному) необходимо преобразовать пространственный сигнал во временной, что осуществляется при относительном перемещении объекта или его изображения и фильтра. При таком перемещении поток от объекта, т.е. сигнал, модулируется. Для уменьшения до минимума полосы временных частот, занимаемой сигналом, и, следовательно, для уменьшения влияния шумов целесообразно получить гармоническую модуляцию сигнала. При использовании указанной выше конструкции узкопольного пространственного фильтра, состоящего из одной прозрачной ячейки, при просмотре всего поля обзора возникает импульсная модуляция потока с широкой полосой спектра сигнала. Достаточно хорошее приближение к непрерывной гармонической модуляции достигается с помощью периодической структуры растра, ячейки (полупериоды) которого близки по форме и размерам изображению объекта.
Если рассмотреть выражение (12) для частотной характеристики оптимального пространственного фильтра, то его можно представить в виде двух основных последовательно включенных звеньев: согласованного фильтра с характеристикой S*(jwx, jwy), который обеспечивает преимущественное пропускание спектра полезного сигнала, и помехоподавляющего звена с характеристикой 1/Фш(wx, wy).
Последнее должно выполнять функции пространственного дифференцирующего звена, причем порядок дифференцирования определяется видом шума.
При изотропном фоне, спектр которого описывается кубической гиперболой вида Фш(wa)~, помехоподавляющее звено должно быть дифференциатором третьего порядка. При фоне со спектром типа квадратической гиперболы (Фш(wa)~) требуется дифференцирование второго порядка. При изменении модели спектра фона, в том числе при учете анизотропии фона, частотная характеристика оптимального пространственного фильтра и максимально достижимое значение отношения сигнал-шум ? могут весьма заметно меняться.
Наибольшие трудности при оптимальной пространственной фильтрации или приближении к ней (квазиоптимальная фильтрация) вызывает необходимость иметь пространственное дифференцирование по полю.
При дифференцировании электрического сигнала, снимаемого с выхода приемника, можно получить непрерывное дифференцирование. Однако при этом перевод сигнала в электрическую форму приводит к добавлению шумов приемника и предварительного усилителя к сигналам от объекта и помех. Кроме того, для идеального непрерывного дифференцирования (выборки d-функции) необходимо иметь бесконечно большую полосу пропускания, что невозможно и нецелесообразно с точки зрения помехозащищенности.
При использовании в качестве пространственного фильтра растра или мозаичного приемника непрерывное дифференцирование (т.е. дифференцирование в каждой точке поля) осуществить невозможно. Поэтому здесь используется образование конечных разностей сигналов в дискретных точках поля, находящихся на малых, но конечных расстояниях друг от друга. Конечные разности образуются путем членения поля на отдельные элементы, размер которых обычно согласовывается с размером изображения источника полезного сигнала, и придания этим элементам различных весов (коэффициента пропускания ячеек растра, чувствительности отдельных площадок мозаичного приемника).
В качестве примера на рис. 8 приведены распределения весов:
- для подавления фона с ярко выраженной анизотропией (рис. 8, а);
- для двух вторых разностей и ортогонально симметричного фона (рис. 8, б);
- для разностей четвертого порядка по трем осям симметрии (рис. 8, в).
Рис. 8. Распределение весов отдельных элементов пространственных фильтров
В качестве простейшего примера пространственной фильтрации с помощью подобных фильтров рассмотрим задачу выделения малоразмерного излучателя на фоне крупноразмерных помех с помощью трехэлементного фильтра с распределением весов -1/2, 1, - 1/2 (рис. 9). Стрелкой на рис. 9, а показано направление сканирования, в процессе которого объект Об и помехи последовательно перекрываются элементами фильтра Ф. Если сигнал uп1 от близкой к изотропной помехи П1 подавляется таким фильтром достаточно эффективно, то этого нельзя сказать про сигнал uп2 от анизотропной помехи П2. Для подавления последнего требуется либо ориентировать такой фильтр по направлению, в котором вытянута помеха, либо применять более сложные фильтры, например, с распределением весов, представленным на рис. 8, б, в. Как это следует из рис. 9, помеха не подавляется и в том случае, если амплитуда сигнала от помехи uп1/2 больше порога срабатывания uпор, т.е. сравнима с амплитудой uоб или больше ее. Это имеет место, когда контраст между помехой и фоном больше чем в 2 раза контраста между объектом и фоном.
Рис. 9. Простейший пример пространственной фильтрации:
а - схема сканирования фильтром поля обзора; б - сигналы на выходе фильтра
Во многих случаях приходится увеличивать порядок дифференцирования, т.е. использовать большее число элементов фильтра с меньшими значениями разностей весов периферийных элементов, а также применять другие методы селекции, например спектральную и временную.
Отметим, что поскольку пропускание оптических растров не может быть отрицательным, получение отрицательных весов для них невозможно. Поэтому реализовать структуры, представленные на рис. 8 и им подобные, удается только с помощью многоэлементных приемников, используя для отдельных элементов приемника усилители с изменяющимися в соответствии со знаками и величинами весовых коэффициентов (весов) коэффициентами усиления.
Оптимальные значения весов отдельных элементов пространственного фильтра, определяемые функцией взаимной двумерной ковариации смеси сигнала и шума и среднего значения сигнала, подбирают в соответствии с приведенным выше правилом так, чтобы минимизировать среднее квадратическое значение шума.
Так как источник сигнала может занимать произвольное положение в угловом поле прибора, фаза сигнала, описываемая экспоненциальным членом в выражении для частотной характеристики оптимального фильтра (6), оказывается неизвестной. Поэтому максимум отношения сигнал-шум обеспечивается путем установки большого числа параллельно работающих оптимальных пространственных фильтров, различающихся только фазовыми множителями, что при работе в угловом поле, превышающем размер изображения источника полезного сигнала, достигается также с помощью многоэлементных структур (растров, мозаичных приемников и др.). Часто вместо сплошной двумерной структуры используется один оптимальный фильтр (один элемент), совершающий последовательный просмотр требуемого поля. Однако такому фильтру свойственны существенные недостатки, отмеченные выше.
Наряду с использованием пространственных фильтров-дифференциаторов в виде многоэлементных приемников, достаточно часто на практике применяют пространственные фильтры в виде растров с чередующимися прозрачными и непрозрачными участками (см. § 9.5), т.е. с распределением весов +1, 0, +1, 0 и т.д., или многоэлементные приемники, пикселы которых включены в чередующемся порядке в противофазе, т.е. распределение весов имеет вид+1, -1, +1, -1 и т.д. Даже такие упрощенные фильтры часто достаточно эффективно селектируют малоразмерные объекты на фоне крупноразмерных помех, хотя они критичны к большим перепадам яркости помех («внутри» помехи или на границах ее).
Фильтры на мозаичных приемниках более эффективны, чем фильтры-растры. В этом легко убедиться, если провести графическое построение сигналов аналогично тому, как это сделано на рис. 9. Кроме того, фронт и спад сигналов, снимаемых с мозаичных приемников, при «входе» фильтра на изображение помехи и «сходе» с него значительно меньше, а амплитуда входного сигнала при использовании многоэлементных приемников с противофазным включением вдвое больше амплитуды сигнала от того же малоразмерного объекта при применении аналогичного по рисунку растра. В то же время мозаичным приемникам свойственны уже отмеченные выше недостатки (дороговизна, сложность изготовления малоразмерных элементов с одинаковыми параметрами, внесение дополнительных шумов схемы считывания и обработки сигналов и ряд других).
При увеличении размеров поля обзора, которое должен исследовать ОЭП, шумы от протяженного фона увеличиваются, а энергия малоразмерного сигнала остается постоянной. Эффективность фильтрации в этом случае понижается. Действительно, рассмотренные выше типовые пространственные фильтры-растры (см. рис. 9.6) обладают хорошим пропусканием в области низких пространственных частот, т.е. там, где велико влияние протяженного фона. Отличительной особенностью пространственной фильтрации в некогерентных системах по сравнению, например, с фильтрацией, применяемой в радиотехнике и радиолокации, является то, что она осуществляется в пространственной области, а не в частотной. При этом техническая реализация фильтра с пространственной структурой со сложным, как правило, законом изменения прозрачности или чувствительности по его площади (в соответствии со сложной двумерной структурой сигнала) - задача трудная, а порой и невыполнимая.
Важно отметить, что сложность и возможные изменения форм и размеров изображений в процессе работы ОЭП, например, вследствие изменения ракурсов наблюдения, дальности и других причин приводят к нестационарности пространственно-частотных спектров таких излучателей. В этом заключается существенное отличие ОЭП от радиоэлектронных приборов, для которых случайные изменения спектров сигналов и помех, как правило, менее значительны. По этой же причине в ОЭП часто невозможно с помощью пространственных фильтров полностью «отсечь» спектр помех. Учитывая возможные изменения размеров изображения, на практике размер ячеек растров пространственных фильтров выбирают больше размера изображения (порой не менее чем в 3 раза).
Другая специфическая особенность пространственной фильтрации - нестационарность передаточных функций оптических систем по их угловому полю или полю обзора. Полевые аберрации оптических систем приводят к тому, что размер кружка рассеяния меняется по полю, как правило, увеличиваясь к краю. Поэтому в растрах, используемых в качестве пространственных фильтров, часто размер ячеек возрастает от центра к краю.
Пространственно-частотные спектры многих фонов и внешних излучающих помех, например, неба, облачности, наземных ландшафтов, имеют четко выраженный низкочастотный характер [13, 24, 26]. Так, основная энергия в спектре облаков содержится на гармониках не выше восьмой, если за первую гармонику принимается частота просмотра всего углового поля ОЭП, например, частота вращения полудискового анализатора-модулятора. В то же время спектр малоразмерных объектов типа самолетов имеет гармоники выше двенадцатой [1, 30]. Это позволяет эффективно использовать полосовые фильтры в электронном тракте ОЭП (после растра-модулятора и приемника).
В приборах или системах, где требуется улучшить качество изображения наблюдаемого поля, пространственная фильтрация помогает увеличить контраст в изображении информативных структур этого поля. С помощью пространственной фильтрации разделяются низкочастотная и высокочастотная составляющие сигнала, а затем одна из них используется для управления другой.
Если динамический диапазон яркостей отдельных участков наблюдаемого поля велик и превышает динамический диапазон всего ОЭП или отдельных его элементов, например системы отображения, то целесообразно ослабить малоинформативную часть спектра пространственных частот, как правило, низкочастотную, и усилить высокочастотную. Для этого из текущего значения сигнала U (x, y), получаемого на выходе приемника излучения при сканировании поля обзора мгновенным угловым полем ОЭП в сумматоре С1 (рис. 10) вычитается низкочастотная составляющая, например его среднее значение Ucp (x, y), которое соответствует среднему значению яркости по полю обзора или по его части, заметно превышающей мгновенное поле, «образующее» сигнал U (x, y).
Рис. 10. Схема коррекции динамического диапазона яркости изображения
Среднее значение Ucp(x, y) можно получить, применяя фильтр низких частот ФНЧ, например, полевую диафрагму с расположенным за ней приемником, работающим в линейной зоне своей энергетической (или фоновой) характеристики и создающим сигнал, пропорциональный общему, осредненному по полю обзора (или его части) потоку излучения. Осреднение может вестись и в электронном тракте, например на выходе мозаичного приемника.