УЧЕТ ХАРАКТЕРИСТИК ФЛУКТУАЦИЙ ФОНОВОГО ИЗЛУЧЕНИЯ ПРИГОРИЗОНТНОЙ ОБЛАСТИ МОРЯ В АЛГОРИТМАХ ОБРАБОТКИ СИГНАЛОВ ИНФРАКРАСНЫХ ПРИБОРОВ
МОДЕЛИРОВАНИЕ ФОНОЦЕЛЕВОЙ ОБСТАНОВКИ В ОПТИКО-ЭЛЕКТРОННЫХ СИСТЕМАХ РАЗЛИЧНОГО НАЗНАЧЕНИЯ
УДК 551.463.5:535.31
Ю. И. БЕЛОУСОВ, Д. В. ИВАНОВ
УЧЕТ ХАРАКТЕРИСТИК ФЛУКТУАЦИЙ ФОНОВОГО ИЗЛУЧЕНИЯ ПРИГОРИЗОНТНОЙ ОБЛАСТИ МОРЯ В АЛГОРИТМАХ ОБРАБОТКИ
СИГНАЛОВ ИНФРАКРАСНЫХ ПРИБОРОВ
Исследована зависимость закона распределения флуктуаций инфракрасной пригоризонтной области моря от условий наблюдения и параметров оптической системы применительно к разработке алгоритмов автоматической обработки сигналов в оптико-электронных системах. Даны рекомендации относительно выбора методов адаптации алгоритмов обработки к текущим условиям наблюдения.
Ключевые слова: обработка изображений, оптико-электронные системы, фоноцелевые модели.
В настоящее время повышение эффективности оптико-электронных систем (ОЭС) напрямую связывается с увеличением их чувствительности и разрешающей способности. Однако наращивание только потенциальной пороговой чувствительности оптических каналов уже не дает пропорционального прироста их эффективности. Одновременно с увеличением чувствительности каналов ОЭС снижается их помехозащищенность, поскольку возрастает вклад фоновых и организованных помех, на фоне которых необходимо выделять пороговые сигналы от целей.
При разработке алгоритмов обнаружения оптических сигналов необходимо в максимальной степени ориентироваться на реальные статистические характеристики фоновых помех, поскольку использование простейших моделей фоновых помех может привести к результату, который будет оптимален в абстрактном случае, но в реальных условиях приведет к снижению эффективности ОЭС в целом.
В настоящей статье рассматривается вопрос использования известной „модели ИКизлучения взволнованной поверхности моря“ [1] при оценке зависимости статистических характеристик флуктуаций фонового излучения от условий наблюдения пригоризонтной области моря с помощью ИК-приборов применительно к разработке алгоритмов автоматической обработки сигналов.
Общая сложность задачи состоит в том, что целый ряд параметров ИК-канала ОЭС зависит от того фона, на котором наблюдается цель. Эта зависимость неоднозначна и определяется как характеристиками условий наблюдения, так и характеристиками самого прибора. Например, в зависимости от фоновой засветки изменяется дисперсия шумов при регистрации сигнала, т.е. один и тот же по амплитуде оптический сигнал, наблюдаемый на различном
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
44 Ю. И. Белоусов, Д. В. Иванов
фоне, обусловливает разное отношение сигнала к внутреннему шуму на выходе фотоприемного устройства, причем статистические данные флуктуаций фонового излучения и внутренних шумов могут либо быть близкими, либо существенно различаться.
При этом преобладающими могут оказаться либо фоновые шумы, либо внутренние. Их совокупное воздействие заметным образом различается. Эта проблема становится все более существенной по мере повышения чувствительности ИК-прибора. Наиболее характерно ее проявление при наблюдении пространства вблизи линии морского горизонта.
Непосредственно у линии горизонта флуктуации фонового излучения относительно малы по амплитуде и соизмеримыми с флуктуациями внутренних шумов прибора. Статистические характеристики суммарного шума соответствуют характеристикам суммы двух независимых случайных процессов с законами распределения, близкими к нормальному.
Периферийные зоны поля обзора в большинстве ситуаций отличаются от области линии горизонта как по среднему значению, так и по характеру и размаху флуктуаций яркости. Изза этого изменяется соотношение вкладов внешних и внутренних шумов. Доминирующими становятся фоновые шумы, а их распределение отличается от нормального закона, и для них алгоритмы, оптимальные при выделении сигналов из шумов с нормальным законом распределения, могут давать существенный проигрыш по сравнению с другими. При изменении соотношения внутренних и внешних шумов вероятность правильного обнаружения может изменяться незначительно, в то время как вероятность ложных тревог изменяется на несколько порядков. Это положение в принципе характерно для всех ситуаций наблюдения в пригоризонтной области моря, но степень различия характеристик фоноцелевой обстановки в центре и на периферии поля зрения ИК-канала зависит от его величины и будет разной, например, при разной кратности увеличения оптической системы одного и того же прибора. В режиме наблюдения за надводной обстановкой в широком поле зрения амплитуда фоновых помех велика, а в режиме сопровождения и выработки координат цели в узком поле зрения флуктуации существенно сглаживаются. Кроме того, в ИК-приборах, как правило, сохраняется количество элементов в изображении. И тогда при изменении поля обзора происходит разное усреднение в пределах одного пиксела (элементарного поля зрения). Одновременно изменяются и амплитуда полезного сигнала, и характеристики воспринимаемых фоновых помех. Если в приборе производится автоматическая обработка сигналов, то она должна быть адаптивной, иначе приводит в каком-то случае к проигрышу.
Предметом исследования настоящей статьи являются вид и параметры закона распределения флуктуаций ИК-излучения взволнованной поверхности моря (ВПМ). В общем виде полученные результаты применимы как для сканирующих, так и для несканирующих ИКприборов. Три особенности существенным образом влияют на применение разработанных для них в [1] аналитических выражений к оценкам для конкретных параметров ОЭС.
Первая особенность заключается в резкой зависимости вида и параметров закона распределения флуктуаций фонового ИК-излучения морской поверхности от угла места вблизи линии горизонта.
Вторая особенность связана с тем, что в реальных приборах происходит усреднение пространственных статистических характеристик в пределах конечного мгновенного поля зрения прибора, причем степень усреднения также зависит от угла места линии визирования. В общем случае для каждой строки изображения флуктуации фонового излучения моря имеют разные статистические характеристики, так как каждый элемент фотоприемника усредняет флуктуации ИК-излучения с разных по площади участков ВПМ, находящихся на различных расстояниях от прибора.
Третья особенность состоит в том, что при изменении углового поля обзора существенным образом изменяется гистограмма распределения яркости фонового сигнала в кадре как по количеству и положению максимумов, так и по ее ширине (динамическому диапазону).
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
Учет характеристик флуктуаций фонового излучения
45
В настоящей статье приводятся результаты исследования зависимости вида и параметров закона распределения флуктуаций ИК-излучения ВПМ от размеров мгновенного поля зрения (пиксела) прибора при различных углах места линии визирования пригоризонтной области моря.
Расчеты базируются на приведенных в [1—3] выражениях для плотности распределения вероятности w(L) и параметров закона распределения флуктуаций спектральной энергетической яркости (СЭЯ) L [Вт/см2⋅ср⋅мкм] ИК-излучения участка ВПМ. Они получены в рамках стохастической фацетной модели в [1—3] и относятся к „точечным“ характеристикам излучения моря, т.е. описывают параметры оптического поля столь малого участка поверхности, что единственное требование к его размерам состоит в том, чтобы он содержал достаточное для применения статистических методов количество фацетов. В системе координат, связанной с ИК-прибором, этот участок наблюдается под углами Ф и α. Угол Ф отсчитывается в вертикальной плоскости относительно угла визирования линии горизонта Фг. Угол α в горизонтальной плоскости отсчитывается относительно главного направления распространения волнения.
Применительно к рассматриваемой задаче в пределах поля зрения прибора (δв — в вертикальной и δг — в горизонтальной плоскости) происходит оптическое усреднение статистических характеристик флуктуаций излучения ВПМ. Формально это можно выразить интегри-
рованием значения СЭЯ в пределах углового мгновенного поля зрения ω= δв ×δг , но тогда
полученная величина имеет размерность спектральной энергетической светимости М [Вт/см2⋅мкм] в направлении (Ф, α) :
M (Φ, α) = ∫∫ L (Φ, α) dδвdδг .
ω
При малых значениях ω процедуру интегрирования с определенной погрешностью можно заменить умножением подынтегральной функции на значение области интегрирова-
ния δв ×δг . Для значений углов ω, в пределах которых подынтегральная функция изменяется
заметно, необходимо осуществлять интегрирование без упрощений. Из исследований [1] известно, что в горизонтальной плоскости характеристики флук-
туаций излучения моря изменяются достаточно монотонно по характеру и весьма мало по значению при изменениях угла визирования α в довольно больших пределах. Поэтому для всех характеристик допустимо интегрирование по α заменить умножением на δг.
Однако в вертикальной плоскости все статистические параметры имеют сложный характер в зависимости от угла наблюдения Ф, а воспринимаемое наблюдателем ИК-излучение
подстилающей поверхности — поверхности моря — Mвосп (Φ, α) ослабляется в τ(λ, Φ) раз
атмосферной трассой, длина которой также зависит от угла визирования Ф и высоты наблюдателя. Кроме того, собственное ИК-излучение этой трассы добавляется к ослабленному ею же излучению моря, причем коэффициент излучения атмосферы ε связан с коэффициентом ее пропускания
ε(λ, Φ) =1−τ(λ,Φ) ,
поэтому при больших значениях δв необходимо осуществлять интегрирование без упрощений
{ }Mвосп (Φ, α) = ∫∫ L (Φ, α) τ(λ,Φ)+ L0в ⎡⎣1−τ(λ,Φ)⎦⎤ dδвdδг , ω
где L0в — яркость излучения абсолютно черного тела (АЧТ) с температурой воздуха.
Чем ниже расположена ОЭС над уровнем моря, тем сильней зависят статистические характеристики наблюдаемого фона от угла места линии визирования. Например, для высоты 25 м над уровнем моря при изменении угла визирования непосредственно ниже линии гори-
зонта на 0,2° расстояние от поверхности моря до наблюдателя изменяется от 17 до 4 км.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
46 Ю. И. Белоусов, Д. В. Иванов
Поэтому вариации характеристик фонового излучения этой зоны во многом обусловлены из-
менениями оптических свойств атмосферной трассы.
В зависимости от состояния небосвода, соотношения температуры моря (Тs) и воздуха (Та), радиационная температура излучения морского фона может быть больше, равна или
меньше температуры излучения АЧТ с температурой воды ( L0s ). Поэтому, например, слабо
нагретая цель может изменять знак своего оптического контраста по мере сближения из-за
того, что наблюдается на фоне разных по яркости фрагментов подстилающей поверхности —
поверхности моря.
Анализ результатов расчетов статистических характеристик флуктуаций воспринимае-
мого ИК-излучения фона проводится для условий сплошной облачности (СО) и ясного небо-
свода (ЯН) при соотношениях Тs < Та и Тs > Та. Из совокупности статистических характеристик ИК-излучения ВПМ как фона наиболее
просто выражается среднее значение флуктуаций ИК-излучения моря (первого момента) [3]:
{( ) }M
СО восп
(Φ,α)= ω
L0s −[L0s − L0a ]A(Φ, α) τ(λ, Φ)+ L0а ⎣⎡1−τ(λ, Φ)⎦⎤ ,
{( ) }M
ЯН восп
(Φ,α) =ω
L0s −[L0s −L0a ]A(Φ,α)−[L0a −Lз ]Z(Φ,α) τ(λ,Φ)+L0a ⎡⎣1−τ(λ,Φ)⎤⎦ .
Здесь L0a — СЭЯ АЧТ при температуре воздуха; Lз — СЭЯ излучения небосвода в зените;
А(Ф,α) и Z(Ф,α) — специальные функции, описывающие эффективное спектральное отраже-
ние ИК-излучения от ВПМ с учетом трехмерного характера волнения и эффекта затенения и
в общем случае — поляризацию ИК-излучения.
Выражение для второго момента — дисперсии флуктуаций фонового излучения моря
D(L) — полученное в статье [1], преобразуется в соответствии с рекомендациями [4]: „по-
стоянный множитель может быть вынесен за знак дисперсии, если этот множитель возвести
при этом в квадрат “. Если не рассматривать (в первом приближении) флуктуации собствен-
ного излучения атмосферной трассы, то „постоянными множителями“ являются значения
мгновенного поля зрения прибора и коэффициент спектрального пропускания атмосферы.
Формально при преобразованиях они играют роль масштабных коэффициентов и не влияют
на вид пространственно-временных зависимостей дисперсии флуктуаций воспринимаемого
излучения D[Mвосп] при наблюдении объектов на фоне морской поверхности:
D[
M
СО восп
(Φ, α)]=[ωτ(λ, Φ)]2 (L0s − L0a )2 D1 ,
D[M
ЯН восп
(Φ, α)]=[ωτ(λ, Φ)]2[(L0s − L0a )2 D1 + 2(L0s − L0a )(L0a − Lз )D2 +(L0a − Lз )2 D3 ] .
D1—D3 — специальные функции [1], имеющие смысл дисперсий — флуктуаций коэффициентов отражения отдельных фацетов относительно „средних“ („эффективных“) коэффици-
ентов отражения единичных участков ВПМ (D1 и D3) и второго смешанного момента (D2) флуктуаций коэффициентов отражения фацетов относительно соответствующих средних
значений с учетом трехмерного характера волнения и эффекта затенения скользящих лучей.
Следует отметить связанную с поляризацией ИК-излучения особенность оценки стати-
стических характеристик оптического поля ВПМ. В расчетные выражения среднего значения
СЭЯ входят формулы коэффициентов отражения для каждой из двух поляризационных со-
ставляющих, а выражение для „неполяризованного излучения“ записывается как полусумма
модулей комплексных коэффициентов отражения воды. Для дисперсии флуктуаций яркости
коэффициент отражения входит во второй степени, поэтому методически правильно рассчи-
тывать по отдельности характеристики флуктуаций для каждой из составляющих поляриза-
ции, а дисперсию флуктуаций суммарного излучения определять как среднее геометрическое
значение дисперсий флуктуаций поляризационных составляющих ИК-излучения ВПМ.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
Учет характеристик флуктуаций фонового излучения
47
Наиболее сложные трансформации претерпевает выражение для закона распределения
флуктуаций воспринимаемого ИК-излучения морского пригоризонтного фона w(Мвосп) по отношению к найденному в [1] выражению для плотности распределения вероятности флуктуаций
энергетической яркости ВПМ w(Lвосп). Помимо введения масштабных множителей ω и τ(λ, Φ)
дополнительно нужно учесть собственное излучение атмосферной трассы. Полагая его независи-
мым от физических факторов, обусловливающих флуктуации ИК-излучения ВПМ, оно формаль-
но трактуется как „сдвиг“ плотности распределения вероятности по оси абсцисс (Мвосп). Причем в зависимости от соотношения значений температуры воздуха и воды этот „сдвиг“ может быть на-
правлен как в сторону уменьшения, так и в сторону увеличения абсолютных значений
Mвосп (Φ, α) :
( )w
⎡ ⎣
M
восп
(Φ , α
)⎤⎦
=
w{[l τ−1
−
L0a
(1−
τ(λ,
Φ
) )]ω −1}
τ (λ,Φ )ω
−1 .
Здесь l — нормированная СЭЯ для подансамбля облачности СО или ЯН [1, 2].
Приведенные выше аналитические выражения позволяют расчетным путем выявить
особенности фоноцелевой обстановки применительно к выработке наиболее рациональных
алгоритмов обработки ИК-изображений. Результаты расчетов, послужившие основой для
анализа и выводов, представлены на рис. 1—4.
d, м h, м
16 000
3000
12 000
2000
8 000 4 000
1000
00
89,6 89,65 89,7 89,75 89,8 89,85 α, …° Рис. 1
На рис. 1 приведена зависимость расстояния (d, сплошная кривая) до наблюдаемого участка поверхности и продольного размера пиксела от угла места линии визирования (h, пунктир). Высота размещения прибора 25 м, мгновенное поле зрения 0,2 мрад.
DЯН, 103(мВт2⋅с)/ (см2⋅мкм⋅ср)
DСО, 103(мВт2⋅с)/ (см2⋅мкм⋅ср)
12
0,8
8
4 0,4
00
70
75 80
85 α, …°
Рис. 2
На рис. 2 представлена зависимость дисперсии флуктуаций яркости от угла места линии визирования (сплошная кривая — DЯН, пунктир — DСО).
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
48 Ю. И. Белоусов, Д. В. Иванов
На рис. 3 приведена зависимость плотности вероятности СЭЯ ИК-излучения ВПМ от угла места линии визирования при скорости ветра 5 м⋅с–1, ненаправленном волнении и сплошной облачности, Та = 285, Тs = 290 K (Φ — 1 — 89, 2 — 85, 3 — 80, 4 — 70°).
w
120
4 80
3 40 2
1
0 0,76 0,78 0,8 0,82 0,84 L, мВт/(мкм⋅ср⋅см2)
Рис. 3
На рис. 4 представлено распределение воспринимаемой яркости излучения ВПМ для различной длины атмосферной трассы (1 — 0, 2 — 1, 3 — 2, 4 — 5 км). Сплошная облачность, Ф=85°.
w
120 4
80 23
1 40
0
0,76 0,78
0,8 0,82 0,84 L, мВт/(мкм⋅ср⋅см2)
Рис. 4
Выводы можно сформулировать в следующем виде.
1. При наблюдении надводной поверхности с помощью дискретного ФПУ (матричного
или сканирующего) для каждой строки изображения пригоризонтной области моря статисти-
ческие характеристики флуктуаций излучения фона различны как по виду, так и по парамет-
рам закона распределения флуктуаций.
2. В ОЭС с изменением углового поля зрения (дискретным или плавным) необходима
адаптация алгоритмов обработки изображений к изменяющимся статистическим характери-
стикам флуктуаций излучения фона.
3. Поскольку большинство алгоритмов обработки сигналов ориентируется на гауссово
распределение шумов, то важно количественное определение параметров условий, при кото-
рых они становятся негауссовыми и встроенные алгоритмы обработки сигналов становятся
неоптимальными.
4. Для задач автоматического выделения сигналов на фоне внутренних и внешних шу-
мов изменение закона распределения помех при выборе процедуры оптимальной фильтрации
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
Система отображения и наблюдения для тренажерного комплекса
49
и обработки сигналов или при выработке решающих правил обусловливает необходимость использования адаптивных алгоритмов.
СПИСОК ЛИТЕРАТУРЫ
1. Белоусов Ю. И., Иванов Д. В., Утенков А. Б. Модель ИК излучения взволнованной поверхности моря // Оптич. журн. 1994. № 10. С. 22—26.
2. Белоусов Ю. И., Иванов Д. В., Утенков А. Б. Метод расчета плотности вероятности яркости ИК излучения взволнованной поверхности моря в диапазоне 3—14 мкм // Оптич. журн. 1999. № 8. С. 36—42.
3. Белоусов Ю. И., Демидов Е. Ф. Средняя составляющая ИК излучения моря // Изв. АН СССР. Сер. Физика атм. и океана. 1982. Т. 18, № 2. С. 206—211.
4. Левин Б. Р. Теоретические основы статистической радиотехники. М.: Сов. радио, 1969.
Сведения об авторах
Юрий Иванович Белоусов
— д-р техн. наук, профессор; Центральный научно-исследовательский
институт „Комета“, Научно-проектный центр оптоэлектронных ком-
плексов наблюдения, Санкт-Петербург; зам. директора;
E-mail: lab6@eoss.ru
Дмитрий Владимирович Иванов — канд. техн. наук; Центральный научно-исследовательский институт
„Комета“, Научно-проектный центр оптоэлектронных комплексов
наблюдения, Санкт-Петербург; ст. науч. сотр.;
E-mail: lab6@eoss.ru
Рекомендована кафедрой оптико-электронных приборов
Поступила в редакцию 29.05.08 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
УДК 551.463.5:535.31
Ю. И. БЕЛОУСОВ, Д. В. ИВАНОВ
УЧЕТ ХАРАКТЕРИСТИК ФЛУКТУАЦИЙ ФОНОВОГО ИЗЛУЧЕНИЯ ПРИГОРИЗОНТНОЙ ОБЛАСТИ МОРЯ В АЛГОРИТМАХ ОБРАБОТКИ
СИГНАЛОВ ИНФРАКРАСНЫХ ПРИБОРОВ
Исследована зависимость закона распределения флуктуаций инфракрасной пригоризонтной области моря от условий наблюдения и параметров оптической системы применительно к разработке алгоритмов автоматической обработки сигналов в оптико-электронных системах. Даны рекомендации относительно выбора методов адаптации алгоритмов обработки к текущим условиям наблюдения.
Ключевые слова: обработка изображений, оптико-электронные системы, фоноцелевые модели.
В настоящее время повышение эффективности оптико-электронных систем (ОЭС) напрямую связывается с увеличением их чувствительности и разрешающей способности. Однако наращивание только потенциальной пороговой чувствительности оптических каналов уже не дает пропорционального прироста их эффективности. Одновременно с увеличением чувствительности каналов ОЭС снижается их помехозащищенность, поскольку возрастает вклад фоновых и организованных помех, на фоне которых необходимо выделять пороговые сигналы от целей.
При разработке алгоритмов обнаружения оптических сигналов необходимо в максимальной степени ориентироваться на реальные статистические характеристики фоновых помех, поскольку использование простейших моделей фоновых помех может привести к результату, который будет оптимален в абстрактном случае, но в реальных условиях приведет к снижению эффективности ОЭС в целом.
В настоящей статье рассматривается вопрос использования известной „модели ИКизлучения взволнованной поверхности моря“ [1] при оценке зависимости статистических характеристик флуктуаций фонового излучения от условий наблюдения пригоризонтной области моря с помощью ИК-приборов применительно к разработке алгоритмов автоматической обработки сигналов.
Общая сложность задачи состоит в том, что целый ряд параметров ИК-канала ОЭС зависит от того фона, на котором наблюдается цель. Эта зависимость неоднозначна и определяется как характеристиками условий наблюдения, так и характеристиками самого прибора. Например, в зависимости от фоновой засветки изменяется дисперсия шумов при регистрации сигнала, т.е. один и тот же по амплитуде оптический сигнал, наблюдаемый на различном
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
44 Ю. И. Белоусов, Д. В. Иванов
фоне, обусловливает разное отношение сигнала к внутреннему шуму на выходе фотоприемного устройства, причем статистические данные флуктуаций фонового излучения и внутренних шумов могут либо быть близкими, либо существенно различаться.
При этом преобладающими могут оказаться либо фоновые шумы, либо внутренние. Их совокупное воздействие заметным образом различается. Эта проблема становится все более существенной по мере повышения чувствительности ИК-прибора. Наиболее характерно ее проявление при наблюдении пространства вблизи линии морского горизонта.
Непосредственно у линии горизонта флуктуации фонового излучения относительно малы по амплитуде и соизмеримыми с флуктуациями внутренних шумов прибора. Статистические характеристики суммарного шума соответствуют характеристикам суммы двух независимых случайных процессов с законами распределения, близкими к нормальному.
Периферийные зоны поля обзора в большинстве ситуаций отличаются от области линии горизонта как по среднему значению, так и по характеру и размаху флуктуаций яркости. Изза этого изменяется соотношение вкладов внешних и внутренних шумов. Доминирующими становятся фоновые шумы, а их распределение отличается от нормального закона, и для них алгоритмы, оптимальные при выделении сигналов из шумов с нормальным законом распределения, могут давать существенный проигрыш по сравнению с другими. При изменении соотношения внутренних и внешних шумов вероятность правильного обнаружения может изменяться незначительно, в то время как вероятность ложных тревог изменяется на несколько порядков. Это положение в принципе характерно для всех ситуаций наблюдения в пригоризонтной области моря, но степень различия характеристик фоноцелевой обстановки в центре и на периферии поля зрения ИК-канала зависит от его величины и будет разной, например, при разной кратности увеличения оптической системы одного и того же прибора. В режиме наблюдения за надводной обстановкой в широком поле зрения амплитуда фоновых помех велика, а в режиме сопровождения и выработки координат цели в узком поле зрения флуктуации существенно сглаживаются. Кроме того, в ИК-приборах, как правило, сохраняется количество элементов в изображении. И тогда при изменении поля обзора происходит разное усреднение в пределах одного пиксела (элементарного поля зрения). Одновременно изменяются и амплитуда полезного сигнала, и характеристики воспринимаемых фоновых помех. Если в приборе производится автоматическая обработка сигналов, то она должна быть адаптивной, иначе приводит в каком-то случае к проигрышу.
Предметом исследования настоящей статьи являются вид и параметры закона распределения флуктуаций ИК-излучения взволнованной поверхности моря (ВПМ). В общем виде полученные результаты применимы как для сканирующих, так и для несканирующих ИКприборов. Три особенности существенным образом влияют на применение разработанных для них в [1] аналитических выражений к оценкам для конкретных параметров ОЭС.
Первая особенность заключается в резкой зависимости вида и параметров закона распределения флуктуаций фонового ИК-излучения морской поверхности от угла места вблизи линии горизонта.
Вторая особенность связана с тем, что в реальных приборах происходит усреднение пространственных статистических характеристик в пределах конечного мгновенного поля зрения прибора, причем степень усреднения также зависит от угла места линии визирования. В общем случае для каждой строки изображения флуктуации фонового излучения моря имеют разные статистические характеристики, так как каждый элемент фотоприемника усредняет флуктуации ИК-излучения с разных по площади участков ВПМ, находящихся на различных расстояниях от прибора.
Третья особенность состоит в том, что при изменении углового поля обзора существенным образом изменяется гистограмма распределения яркости фонового сигнала в кадре как по количеству и положению максимумов, так и по ее ширине (динамическому диапазону).
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
Учет характеристик флуктуаций фонового излучения
45
В настоящей статье приводятся результаты исследования зависимости вида и параметров закона распределения флуктуаций ИК-излучения ВПМ от размеров мгновенного поля зрения (пиксела) прибора при различных углах места линии визирования пригоризонтной области моря.
Расчеты базируются на приведенных в [1—3] выражениях для плотности распределения вероятности w(L) и параметров закона распределения флуктуаций спектральной энергетической яркости (СЭЯ) L [Вт/см2⋅ср⋅мкм] ИК-излучения участка ВПМ. Они получены в рамках стохастической фацетной модели в [1—3] и относятся к „точечным“ характеристикам излучения моря, т.е. описывают параметры оптического поля столь малого участка поверхности, что единственное требование к его размерам состоит в том, чтобы он содержал достаточное для применения статистических методов количество фацетов. В системе координат, связанной с ИК-прибором, этот участок наблюдается под углами Ф и α. Угол Ф отсчитывается в вертикальной плоскости относительно угла визирования линии горизонта Фг. Угол α в горизонтальной плоскости отсчитывается относительно главного направления распространения волнения.
Применительно к рассматриваемой задаче в пределах поля зрения прибора (δв — в вертикальной и δг — в горизонтальной плоскости) происходит оптическое усреднение статистических характеристик флуктуаций излучения ВПМ. Формально это можно выразить интегри-
рованием значения СЭЯ в пределах углового мгновенного поля зрения ω= δв ×δг , но тогда
полученная величина имеет размерность спектральной энергетической светимости М [Вт/см2⋅мкм] в направлении (Ф, α) :
M (Φ, α) = ∫∫ L (Φ, α) dδвdδг .
ω
При малых значениях ω процедуру интегрирования с определенной погрешностью можно заменить умножением подынтегральной функции на значение области интегрирова-
ния δв ×δг . Для значений углов ω, в пределах которых подынтегральная функция изменяется
заметно, необходимо осуществлять интегрирование без упрощений. Из исследований [1] известно, что в горизонтальной плоскости характеристики флук-
туаций излучения моря изменяются достаточно монотонно по характеру и весьма мало по значению при изменениях угла визирования α в довольно больших пределах. Поэтому для всех характеристик допустимо интегрирование по α заменить умножением на δг.
Однако в вертикальной плоскости все статистические параметры имеют сложный характер в зависимости от угла наблюдения Ф, а воспринимаемое наблюдателем ИК-излучение
подстилающей поверхности — поверхности моря — Mвосп (Φ, α) ослабляется в τ(λ, Φ) раз
атмосферной трассой, длина которой также зависит от угла визирования Ф и высоты наблюдателя. Кроме того, собственное ИК-излучение этой трассы добавляется к ослабленному ею же излучению моря, причем коэффициент излучения атмосферы ε связан с коэффициентом ее пропускания
ε(λ, Φ) =1−τ(λ,Φ) ,
поэтому при больших значениях δв необходимо осуществлять интегрирование без упрощений
{ }Mвосп (Φ, α) = ∫∫ L (Φ, α) τ(λ,Φ)+ L0в ⎡⎣1−τ(λ,Φ)⎦⎤ dδвdδг , ω
где L0в — яркость излучения абсолютно черного тела (АЧТ) с температурой воздуха.
Чем ниже расположена ОЭС над уровнем моря, тем сильней зависят статистические характеристики наблюдаемого фона от угла места линии визирования. Например, для высоты 25 м над уровнем моря при изменении угла визирования непосредственно ниже линии гори-
зонта на 0,2° расстояние от поверхности моря до наблюдателя изменяется от 17 до 4 км.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
46 Ю. И. Белоусов, Д. В. Иванов
Поэтому вариации характеристик фонового излучения этой зоны во многом обусловлены из-
менениями оптических свойств атмосферной трассы.
В зависимости от состояния небосвода, соотношения температуры моря (Тs) и воздуха (Та), радиационная температура излучения морского фона может быть больше, равна или
меньше температуры излучения АЧТ с температурой воды ( L0s ). Поэтому, например, слабо
нагретая цель может изменять знак своего оптического контраста по мере сближения из-за
того, что наблюдается на фоне разных по яркости фрагментов подстилающей поверхности —
поверхности моря.
Анализ результатов расчетов статистических характеристик флуктуаций воспринимае-
мого ИК-излучения фона проводится для условий сплошной облачности (СО) и ясного небо-
свода (ЯН) при соотношениях Тs < Та и Тs > Та. Из совокупности статистических характеристик ИК-излучения ВПМ как фона наиболее
просто выражается среднее значение флуктуаций ИК-излучения моря (первого момента) [3]:
{( ) }M
СО восп
(Φ,α)= ω
L0s −[L0s − L0a ]A(Φ, α) τ(λ, Φ)+ L0а ⎣⎡1−τ(λ, Φ)⎦⎤ ,
{( ) }M
ЯН восп
(Φ,α) =ω
L0s −[L0s −L0a ]A(Φ,α)−[L0a −Lз ]Z(Φ,α) τ(λ,Φ)+L0a ⎡⎣1−τ(λ,Φ)⎤⎦ .
Здесь L0a — СЭЯ АЧТ при температуре воздуха; Lз — СЭЯ излучения небосвода в зените;
А(Ф,α) и Z(Ф,α) — специальные функции, описывающие эффективное спектральное отраже-
ние ИК-излучения от ВПМ с учетом трехмерного характера волнения и эффекта затенения и
в общем случае — поляризацию ИК-излучения.
Выражение для второго момента — дисперсии флуктуаций фонового излучения моря
D(L) — полученное в статье [1], преобразуется в соответствии с рекомендациями [4]: „по-
стоянный множитель может быть вынесен за знак дисперсии, если этот множитель возвести
при этом в квадрат “. Если не рассматривать (в первом приближении) флуктуации собствен-
ного излучения атмосферной трассы, то „постоянными множителями“ являются значения
мгновенного поля зрения прибора и коэффициент спектрального пропускания атмосферы.
Формально при преобразованиях они играют роль масштабных коэффициентов и не влияют
на вид пространственно-временных зависимостей дисперсии флуктуаций воспринимаемого
излучения D[Mвосп] при наблюдении объектов на фоне морской поверхности:
D[
M
СО восп
(Φ, α)]=[ωτ(λ, Φ)]2 (L0s − L0a )2 D1 ,
D[M
ЯН восп
(Φ, α)]=[ωτ(λ, Φ)]2[(L0s − L0a )2 D1 + 2(L0s − L0a )(L0a − Lз )D2 +(L0a − Lз )2 D3 ] .
D1—D3 — специальные функции [1], имеющие смысл дисперсий — флуктуаций коэффициентов отражения отдельных фацетов относительно „средних“ („эффективных“) коэффици-
ентов отражения единичных участков ВПМ (D1 и D3) и второго смешанного момента (D2) флуктуаций коэффициентов отражения фацетов относительно соответствующих средних
значений с учетом трехмерного характера волнения и эффекта затенения скользящих лучей.
Следует отметить связанную с поляризацией ИК-излучения особенность оценки стати-
стических характеристик оптического поля ВПМ. В расчетные выражения среднего значения
СЭЯ входят формулы коэффициентов отражения для каждой из двух поляризационных со-
ставляющих, а выражение для „неполяризованного излучения“ записывается как полусумма
модулей комплексных коэффициентов отражения воды. Для дисперсии флуктуаций яркости
коэффициент отражения входит во второй степени, поэтому методически правильно рассчи-
тывать по отдельности характеристики флуктуаций для каждой из составляющих поляриза-
ции, а дисперсию флуктуаций суммарного излучения определять как среднее геометрическое
значение дисперсий флуктуаций поляризационных составляющих ИК-излучения ВПМ.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
Учет характеристик флуктуаций фонового излучения
47
Наиболее сложные трансформации претерпевает выражение для закона распределения
флуктуаций воспринимаемого ИК-излучения морского пригоризонтного фона w(Мвосп) по отношению к найденному в [1] выражению для плотности распределения вероятности флуктуаций
энергетической яркости ВПМ w(Lвосп). Помимо введения масштабных множителей ω и τ(λ, Φ)
дополнительно нужно учесть собственное излучение атмосферной трассы. Полагая его независи-
мым от физических факторов, обусловливающих флуктуации ИК-излучения ВПМ, оно формаль-
но трактуется как „сдвиг“ плотности распределения вероятности по оси абсцисс (Мвосп). Причем в зависимости от соотношения значений температуры воздуха и воды этот „сдвиг“ может быть на-
правлен как в сторону уменьшения, так и в сторону увеличения абсолютных значений
Mвосп (Φ, α) :
( )w
⎡ ⎣
M
восп
(Φ , α
)⎤⎦
=
w{[l τ−1
−
L0a
(1−
τ(λ,
Φ
) )]ω −1}
τ (λ,Φ )ω
−1 .
Здесь l — нормированная СЭЯ для подансамбля облачности СО или ЯН [1, 2].
Приведенные выше аналитические выражения позволяют расчетным путем выявить
особенности фоноцелевой обстановки применительно к выработке наиболее рациональных
алгоритмов обработки ИК-изображений. Результаты расчетов, послужившие основой для
анализа и выводов, представлены на рис. 1—4.
d, м h, м
16 000
3000
12 000
2000
8 000 4 000
1000
00
89,6 89,65 89,7 89,75 89,8 89,85 α, …° Рис. 1
На рис. 1 приведена зависимость расстояния (d, сплошная кривая) до наблюдаемого участка поверхности и продольного размера пиксела от угла места линии визирования (h, пунктир). Высота размещения прибора 25 м, мгновенное поле зрения 0,2 мрад.
DЯН, 103(мВт2⋅с)/ (см2⋅мкм⋅ср)
DСО, 103(мВт2⋅с)/ (см2⋅мкм⋅ср)
12
0,8
8
4 0,4
00
70
75 80
85 α, …°
Рис. 2
На рис. 2 представлена зависимость дисперсии флуктуаций яркости от угла места линии визирования (сплошная кривая — DЯН, пунктир — DСО).
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
48 Ю. И. Белоусов, Д. В. Иванов
На рис. 3 приведена зависимость плотности вероятности СЭЯ ИК-излучения ВПМ от угла места линии визирования при скорости ветра 5 м⋅с–1, ненаправленном волнении и сплошной облачности, Та = 285, Тs = 290 K (Φ — 1 — 89, 2 — 85, 3 — 80, 4 — 70°).
w
120
4 80
3 40 2
1
0 0,76 0,78 0,8 0,82 0,84 L, мВт/(мкм⋅ср⋅см2)
Рис. 3
На рис. 4 представлено распределение воспринимаемой яркости излучения ВПМ для различной длины атмосферной трассы (1 — 0, 2 — 1, 3 — 2, 4 — 5 км). Сплошная облачность, Ф=85°.
w
120 4
80 23
1 40
0
0,76 0,78
0,8 0,82 0,84 L, мВт/(мкм⋅ср⋅см2)
Рис. 4
Выводы можно сформулировать в следующем виде.
1. При наблюдении надводной поверхности с помощью дискретного ФПУ (матричного
или сканирующего) для каждой строки изображения пригоризонтной области моря статисти-
ческие характеристики флуктуаций излучения фона различны как по виду, так и по парамет-
рам закона распределения флуктуаций.
2. В ОЭС с изменением углового поля зрения (дискретным или плавным) необходима
адаптация алгоритмов обработки изображений к изменяющимся статистическим характери-
стикам флуктуаций излучения фона.
3. Поскольку большинство алгоритмов обработки сигналов ориентируется на гауссово
распределение шумов, то важно количественное определение параметров условий, при кото-
рых они становятся негауссовыми и встроенные алгоритмы обработки сигналов становятся
неоптимальными.
4. Для задач автоматического выделения сигналов на фоне внутренних и внешних шу-
мов изменение закона распределения помех при выборе процедуры оптимальной фильтрации
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8
Система отображения и наблюдения для тренажерного комплекса
49
и обработки сигналов или при выработке решающих правил обусловливает необходимость использования адаптивных алгоритмов.
СПИСОК ЛИТЕРАТУРЫ
1. Белоусов Ю. И., Иванов Д. В., Утенков А. Б. Модель ИК излучения взволнованной поверхности моря // Оптич. журн. 1994. № 10. С. 22—26.
2. Белоусов Ю. И., Иванов Д. В., Утенков А. Б. Метод расчета плотности вероятности яркости ИК излучения взволнованной поверхности моря в диапазоне 3—14 мкм // Оптич. журн. 1999. № 8. С. 36—42.
3. Белоусов Ю. И., Демидов Е. Ф. Средняя составляющая ИК излучения моря // Изв. АН СССР. Сер. Физика атм. и океана. 1982. Т. 18, № 2. С. 206—211.
4. Левин Б. Р. Теоретические основы статистической радиотехники. М.: Сов. радио, 1969.
Сведения об авторах
Юрий Иванович Белоусов
— д-р техн. наук, профессор; Центральный научно-исследовательский
институт „Комета“, Научно-проектный центр оптоэлектронных ком-
плексов наблюдения, Санкт-Петербург; зам. директора;
E-mail: lab6@eoss.ru
Дмитрий Владимирович Иванов — канд. техн. наук; Центральный научно-исследовательский институт
„Комета“, Научно-проектный центр оптоэлектронных комплексов
наблюдения, Санкт-Петербург; ст. науч. сотр.;
E-mail: lab6@eoss.ru
Рекомендована кафедрой оптико-электронных приборов
Поступила в редакцию 29.05.08 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 8