Книги по разным темам
Pages: | 1 | 2 | Электронный журнал ИССЛЕДОВАНО В РОССИИ 1355 Метод оценки пороговой частотно-энергетической характеристики оптико-электронного тракта Сикорский Д. А. (old@tut.by) НПО УЛептонФ В качестве показателя информативности оптико-электронного тракта (ОЭТ) принято линейное разрешение на местности, которое оценивается с использованием стандартной трехшпальной миры. Метод априорной оценки этого показателя сводится к аналитическому описанию энергетических и частотных преобразований входного воздействия элементами тракта [1]. По результатам аналитического описания составляется сквозное частотноэнергетическое уравнение, решение которого определяет предельноразрешаемую пространственную частоту пред. Уравнение имеет вид:
n 4 Sшп SNR ()=Фпор(Рзад), (1) Т =0 i Sкш i=1 где Sшп - площадь изображения шпалы миры;
Sкш - площадь корреляции шума;
SNR=0 - аппаратное отношение сигнала к шуму на нулевой пространственной частоте, определяемое выражением:
Э =0 SNR=0=, (2) Э здесь Э=0 - разность экспозиций от полуплоскостей функции Хевисайда с заданными альбедо её бьефов;
Э - среднее квадратическое отклонение шумовой составляющей сквозного ОЭТ, выраженное в единицах экспозиции;
n T () - сквозная частотно-контрастная характеристика (ЧКХ) ОЭТ, i i=1 являющаяся произведением ЧКХ отдельных его элементов (звеньев). Далее будет показано, что наиболее критичными элементами тракта, с точки зрения описаний их частотных и энергетических свойств, являются звено анализирующей дискретизации и зрительный анализатор.
Фпор(Рзад) - пороговое отношение сигнала к шуму в изображении шпалы миры, обеспечивающее дешифрирование элемента миры по принятому критерию с заданной вероятностью Рзад.
Искомое линейное разрешение на местности LRM=Dн/(2предFCos), (3) где Dн - дальность наблюдения;
F - фокусное расстояние ОС;
- угол визирования.
Электронный журнал ИССЛЕДОВАНО В РОССИИ 1356 Подробно рассматривая элементы сквозного ОЭТ, с акцентированием внимания на зрительном анализаторе (ЗА), можно отметить его явно нелинейные свойства как в пространственной, так и в энергетической областях. Применение теории линейной фильтрации к описанию этих свойств ЗА возможно лишь для узкого круга задач, с множеством ограничений.
Тем не менее, проектировщики ОЭА все же используют понятие ЧКХ глаза. Изучению ЧКХ глаза или его реакции на синусоидальную волну, т.е.
на тест-объект с синусоидальным распределением яркости по пространственной координате, посвящен ряд работ, выполненных в СССР, США, ФРГ, Франции, Японии и других странах. Один из лучших обзоров этих работ сделан М. А. Островской [2].
Из [1] следует, что форма ЧКХ глаза такова, что она имеет малую величину на очень низких пространственных частотах, затем резко возрастает до максимума и далее спадает примерно по закону Гаусса. В угловых величинах на частотах, превышающих 0,2 мрад-1, она приближенно описывается следующим выражением:
Г( )=exp(-22гл 2), (4) где гл =0,2Е0,3 мрад.
Фактически ЧКХ ЗА можно интерпретировать как огибающую множества узкополосных резонансных фильтров пространственных частот.
В каждом конкретном случае ЗА является оптимальным фильтром для решения данной задачи.
В [7] была экспериментально оценена ЧКХ ЗА в плоскости цифровых изображений объектов простой военной техники (танк, самолет и др.) при оптимальных условиях наблюдения. В эксперименте вероятность идентификации объектов по их изображениям составляла не более 0.8, и синтезирующая апертура описывалась двумерной rect-функцией. Полученная ЧКХ описывалась следующим выражением:
Г( )=exp(-6,4422), (5) где - линейный размер синтезирующей апертуры, совпадающий с шагом дискретизации (периодом УрешёткиФ фотоприемника).
С целью изучения частотных свойств зрительного анализатора при наблюдении цифровых изображений стандартной трёхшпальной миры (тестобъект) автором был проведён аналогичный психофизический эксперимент.
Наблюдателя заранее (до проведения эксперимента) ознакомили с тестобъектом. В ходе эксперимента определялась такая дистанция наблюдения, при которой обеспечивалась наибольшая вероятность идентификации теста по заданному критерию. Такая дистанция называлась оптимальной (Lопт).
Рисунок 1 иллюстрирует схему проведения эксперимента.
Электронный журнал ИССЛЕДОВАНО В РОССИИ 1357 Lорт 1,5Т в среднем ЗА=глаз+мозг Цифровое изображение Рис.1. Наблюдатель устанавливает оптимальное расстояние до изображения теста, при котором обеспечивается наибольшая вероятность его идентификации Эксперимент показал, что в этом случае оптимальное расстояние обеспечивало наблюдение элемента разложения изображения в среднем под углом 2.9Т, что резко отличается от результата (1.5Т ) полученного в [7].
Можно предположить, что выявленное различие объясняется наличием у наблюдателя априорных знаний об объекте наблюдения и явно подтверждает нелинейные частотные свойства ЗА.
Результат эксперимента может быть наглядно интерпретирован в частотной области. На рис.2 приводится одномерная иллюстрация, поясняющая различие в процессе фильтрации спектров дискретизированных безшумовых изображений объектов простой военной техники и элемента стандартной трехшпальной миры зрительным анализатором.
Из рисунка видно, что ЗА при наличии априорной информации об объекте допускает значительные искажения, обусловленные выборкой.
Даже при наблюдении изображений одного и того же объекта, но сформированных с разными периодами выборки и одинаковым уровнем шума, ЧКХ ЗА различны. Так, на рис.3 показано изменение ЧКХ ЗА при наблюдении изображений элемента миры, при различных шагах дискретизации и аппаратном отношении сигнала к шуму в сигнале теста, равного двум.
В энергетической области ЗА также является нелинейным элементом.
Так, многочисленные исследования [3] показывают, что ЗА осуществляет пространственную суммацию некоррелированных шумов изображения.
Разработчики ОЭА оперируют так называемым эффективным (видимым) отношением сигнала к шуму в изображении объекта (SNRэфф). При этом, между SNRэфф и аппаратным отношением сигнала к шуму (SNRапп) выявлена следующая зависимость:
SNRэфф= SNRапп(Sоб/ Sкш)0.5, (6) где Sоб - площадь изображения объекта наблюдения;
Электронный журнал ИССЛЕДОВАНО В РОССИИ 1358 Мультиплицированный спектр Г() изображения объекта простой военной D() техники - Qоб(), л/мм 1/ 2/ -4/ -3/ -2/ -1/ 0 3/ 4/ Спектр Qоб() после НЧ фильтрации визуализатором с апертурой типа rectфункции и ЗА с ЧКХ Г()=exp(-6.4422) 1/ 2/ -4/ -3/ -2/ -1/ 0 3/ 4/ Мультиплицированный спектр изображения элемента стандартной Г() трехшпальной миры - Qмира() D() 1/ 2/ -4/ -3/ -2/ -1/ 0 3/ 4/ Спектр Qмиры() после НЧ фильтрации визуализатором с апертурой типа rectфункции и ЗА с ЧКХ Г()=exp(-1.6622) 1/ 2/ -4/ -3/ -2/ -1/ 0 3/ 4/ Рис.2. Одномерная иллюстрация, поясняющая различие в процессе НЧ фильтрации спектров дискретизированных безшумовых изображений объектов простой военной техники и миры Электронный журнал ИССЛЕДОВАНО В РОССИИ 1359 Изображение элемента стандартной трехшпальной миры, искаженное аддитивным гауссовым шумом при Ш/=3.5, где Ш - ширина шпалы.
Мультиплицированный спектр изображения элемента стандартной трехшпальной миры - Qмира(), при Ш/=3.5.
Г() D() 1/ -1/ Изображение элемента стандартной трехшпальной миры, искаженное аддитивным гауссовым шумом при Ш/==1.41.
Мультиплицированный спектр изображения элемента стандартной трехшпальной миры - Qмира(), при Г() Ш/==1.41.
D() 1/ -3/ -2/ -1/ 0 2/ 3/ Рис.3. Одномерная иллюстрация, отражающая различие ЧКХ ЗА при наблюдении изображений элемента стандартной трехшпальной миры, искаженных аддитивным гауссовым шумом и дискретизированных с разными шагами выборки Электронный журнал ИССЛЕДОВАНО В РОССИИ 1360 Sкш - площадь корреляции шума.
При этом эффект суммации прекращается, когда угловой размер анализируемого изображения более 40 угл. мин.
Более подробно следует остановиться на анализе порогового значения отношения сигнала к шуму в изображении шпалы миры, обеспечивающего дешифрирование её элемента по принятому критерию с заданной вероятностью Рзад. В [7] было определено, что для достижения вероятности дешифрирования элемента миры по Укритерию предельно-разрешаемого элементаФ, равной 0.8, видимое (эффективное) отношение сигнала к шуму в изображении её шпалы должно составлять 8.5 - 9.0.
Данный результат нельзя распространять на весь диапазон условий наблюдения, поскольку он получен при определенных ограничениях. Так, тест-объект в эксперименте представлял собой только вертикальную группу элемента миры, и пределы изменения отношения ширины шпалы к шагу выборки составляли 1.1 - 1.3. Не рассматривались пространственные частоты максимально приближенные к частоте Найквиста. Более того, результат получен для наилучшей фазы выборки.
Также следует отметить, что в существующих методиках эффект действия выборки оценивается частотно-контрастной характеристикой фазы в случае её наихудшего соотношения. С точки зрения теории линейной фильтрации это не совсем корректно, хотя для пространственных частот менее 0.9 частоты Найквиста физический смысл ЧКХ фазы может быть достаточно строго определен. На более высоких частотах форма оптимально восстановленного по выборкам изображения миры резко отличается от синусоидального распределения, и использование ЧКХ в этой области не правомерно.
Проведенный анализ показывает, что раздельное описание частотных и энергетических свойств таких звеньев ОЭТ как УвыборкаФ, с присущей ей фазой и ЗА не всегда корректно и удобно.
Автором предложено осуществлять описание частотно-энергетических свойств звена выборки, визуализатора и зрительного анализатора единым функционалом. Предложенный функционал должен отражать зависимость требуемого значения аппаратного отношения сигнала к шуму в изображении центрального штриха группы элемента миры (в общем случае - в изображении структурной единицы теста) от пространственной частоты элемента. Причем, требуемое значение отношения сигнала к шуму выбирается исходя из обеспечения заданной вероятности дешифрирования элемента миры (теста) по заданному критерию. Такой функционал будем называть пороговой частотно-энергетической характеристикой (ПЧЭХ) оптико-электронного тракта.
Далее по тексту приводится описание метода оценки ПЧЭХ. В качестве иллюстрационного материала использованы изображения элемента миры.
Формировалась серия цифровых изображений элемента миры, которая затем обрабатывалась оператором. На рис. 4 показана блок-схема алгоритма формирования одной такой серии.
Электронный журнал ИССЛЕДОВАНО В РОССИИ 1361 Ввод исходных данных:
SNR - аппаратное отношение сигнала к шуму в изображении шпалы миры;
шп/ - отношение ширины шпалы элемента миры к шагу выборки;
- параметр сквозной апертурной функции ОЭТ, заданной гауссоидой вращения;
Nзад - заданное количество синтезируемых изображений.
N=Синтез изображения элемента стандартной трехшпальной миры Поворот элемента миры на случайный угол =rnd*360 относительно центра изображения, где rnd - датчик случайных чисел равномерно распределенных на интервале 0-1.
Низкочастотная фильтрация повернутого изображения элементами ОЭТ, определение величины полезного сигнала Sig Выборка с шагом и фазой Хf=int(rnd*), Yf=int(rnd*) Синтезирующая дискретизация апертурой вида rect(x/)*rect(y/) incr N Зашумление результата выборки аддитивным, с гауссовым распределением амплитуд шумом с СКО=Sig/SNR да нет N Сечение изображения y группы x S-S S-S Sig Рис. 5. Иллюстрация, отражающая эффект низкочастотной фильтрации элемента стандартной трехшпальной миры и поясняющая определение величины сигнала Итак, для конкретных значений шп и заданного SNR была сформирована серия цифровых изображений элемента миры. Изображения серии отличались друг от друга фазой выборки и реализациями шумовой компоненты, формируемыми одним и тем же генератором случайных чисел. На рис. 6 приводится фрагмент серии изображений, сформированных при шп =2.33 и SNR=2.7. Серия предъявлялась эксперту (наблюдателю), предварительно ознакомленному с критерием дешифрирования группы миры. К эксперименту привлекались наблюдатели, обладающие практическим опытом анализа цифровых изображений и не имеющие аномальных отклонений в зрении. При проведении эксперимента наблюдателю разрешалось подбирать наилучшую комфортность наблюдения цифровых изображений: изменять яркость, контрастность изображения, дистанцию до него, ракурс наблюдения и уровень внешнего фона. Электронный журнал ИССЛЕДОВАНО В РОССИИ 1363 Эксперт из предложенной к анализу серии отбирал изображения, удовлетворяющие критерию дешифрирования. По полученным результатам оценивалась вероятность дешифрирования элемента, а именно: Электронный журнал ИССЛЕДОВАНО В РОССИИ 1364 Рис. 6. Фрагмент серии цифровых изображений элемента стандартной трехшпальной миры _n Pдеш = Niдеш / (Nс*M), (6) i= где - Niдеш - количество изображений элемента миры, отобранных из анализируемой серии i-ым экспертом по принятому критерию дешифрирования; Nс - количество изображений в серии; M - количество привлекаемых экспертов. Результаты эксперимента в зависимости от SNR, при конкретном шп сводились в график, приведенный в качестве примера на рис. 7. Книги по разным темам