Математическая морфология.
Электронный математический и медико-биологический
журнал. - Т. 8. -
Вып. 1. - 2009. - URL:
http://www.smolensk.ru/user/sgma/MMORPH/TITL.HTМ
http://www.smolensk.ru/user/sgma/MMORPH/N-21-html/TITL-21.htm
http://www.smolensk.ru/user/sgma/MMORPH/N-21-html/cont.htm
УДК 621.391.161
Обнаружение ТЕПЛОВЫХ ОБЪЕКТОВ
на коррелированном атмосферном
фоне
Ó
2009 г. Якименко И. В.,
Жендарев М. В., Гурченков Д. А.
В данной работе изложен способ обнаружения тепловых объектов на
коррелированном атмосферном фоне с помощью инфракрасных теплопеленгаторов с
широким полем зрения. Предлагаемый способ основан на использовании отличий
пространственных спектров излучения точечного теплового объекта и протяженного,
более холодного, атмосферного фона. На основе данного способа обнаружения тепловых объектов разработан алгоритм
двумерного пространственного фильтра, который может быть реализован программно
в видеопроцессоре теплопеленгатора. Алгоритм позволяет автоматизировать процесс
обнаружения воздушных целей (самолетов, вертолетов, крылатых ракет, беспилотных
летательных аппаратов и т. д.) и обеспечить вывод их видеоизображений на монитор
видеоконтрольного устройства в бинарном виде.
Ключевые слова: обнаружение тепловых
объектов, инфракрасные оптико-электронные системы, коррелированный фон
В настоящее время для
теплопеленгаторов применяется амплитудный способ обнаружения точечных тепловых
объектов на атмосферном фоне. Данный способ основан на зрительном разделении
яркости отметок точечного теплового объекта и флюктуаций атмосферного фона.
Получение фоноцелевом изображении
(ФЦИ) на мониторе видеоконтрольного устройства (ВКУ) проводится в соответствии
с видеосигналом, полученным с выхода МПИ. Каждому пикселю ФЦИ в соответствие ставят
ограниченное число дискретных уровней яркости – квантов, пропорциональных
энергетической яркости собственного излучения атмосферы и цели. Чем больше
уровней квантования, тем выше качество изображения.
Особенностью
работы оптико-электронных приборов с широким полем зрения, работающих в
диапазоне 8-13 мкм,
является высокие перепады яркости одновременно наблюдаемых фрагментов
изображения. Например, радиационный контраст кучевого облака на фоне чистого
неба как днём, так и ночью, может достигать 60 К [3]. Аналогичная ситуация возникает
при визировании цели вблизи линии горизонта, когда в угловом поле прицела наблюдается
часть небесной полусферы и наземный ландшафт, который обычно «теплее» небесного.
При чувствительности матрицы <0,1 К и наличии ограничения в выходном
видеосигнале «уровнем белого» осуществляется компрессия сигналов с МПИ. Это, в
свою очередь, приводит к уменьшению сигнала цели на входе ВКУ [4]. ФЦИ, полученное на мониторе
ВКУ теплопеленгатора при наличие цели (отметка от цели находится в разрыве линии),
представлено на рис. 1. Изображение получено в темное время суток при кучевой облачности.
Проблемы видения возникают вследствие потенциально слабой контрастности
ФЦИ, низкого уровня заметности целей и возможного появления пассивных и
активных видов противодействия (маскировка, ложные тепловые цели (ЛТЦ), прожекторы
и т. п.). Шумовые компоненты сигнала, образующиеся вследствие особенностей
применяемых датчиков и неидеальности сигнала передачи данных приводят к
ухудшению изображения точечного теплового объекта на мониторе ВКУ. С целью улучшения
качества изображений были созданы различные способы цифровой фильтрации и обработки двухмерного
изображения. Обработка изображений
проводится в реальном масштабе времени, без прореживания пиксельного потока
видеоданных и без пропуска кадров. Основными способами обработки изображений
являются алгоритмы цифровой коррекции неоднородностей, сглаживание
шумов, повышение контрастов, подчеркивание границ, высокочастотная коррекция,
коррекция дефектных пикселей и т. д., реализованных программно в видеопроцессоре
теплопеленгатора. Отметим, что
перечисленные способы обработки изображений косвенно улучшают
характеристики теплопеленгаторов за счет улучшения визуального восприятия изображения
оператором.
Окончательной оценкой эффективности улучшения качества видеоизображений
должно служить субъективно оцениваемое качество изображения, поскольку конечным
потребителем воспроизводимого изображения является оператор. Зрительные возможности у
разных операторов различны, и на неравномерном изображении облачного атмосферного
фона трудно выделить точечную отметку теплового объекта. Ещё сложнее селектировать
несколько отметок, поэтому при использовании зрительного обнаружения повышается
вероятность ошибок принятия решений.
Рис. 1. Исходное изображение, полученное на мониторе ВКУ
Таким образом, необходимы
новые способы обнаружения воздушных целей, которые способны выделять отметки точечных
тепловых объектов на изображении атмосферного фона без участия оператора. На
конечном этапе обнаружения оператор на мониторе ВКУ должен увидеть изображение,
для которого характерно максимально контрастное представление отметок тепловых
объектов на скомпенсированном равномерном фоне. Для получения такого изображения
достаточно двух уровней квантования. Полученные на экране монитора бинарные
изображения тепловых объектов позволяют определять их угловые координаты и
количественный состав (рис. 2).
Рис. 2. Изображение на мониторе ВКУ, полученное после
фильтрации
Такое выделение тепловых
объектов на изображении атмосферного фона можно отнести к пространственной
фильтрации. Пространственный фильтр предлагается реализовать программно в
видеопроцессоре. Структурная схема теплопеленгатора, реализующая предлагаемую
пространственную фильтрацию, представлена на рис. 3.
Сущность разработанного
авторами способа обнаружения тепловых объектов на атмосферном фоне заключается
в том, что предлагается использовать отличия пространственных спектров
излучения точечных объектов и протяженного, более холодного, атмосферного фона.
Пространственный спектральный анализ атмосферного фона содержит расчет нормированных
пространственных спектральных плотностей мощности флуктуаций излучения. Он
связан с методическими и вычислительными сложностями, которые касаются
оптимизации выбора корреляционного окна.
Согласно теореме Винера - Хинчина, пространственная спектральная
плотность мощности преобразованием Фурье связана с взаимной корреляционной
функцией, поэтому предлагается оценивать спектральную плотность мощности
излучения на основе анализа нормированных коэффициентов взаимной корреляции для
соседних строк (столбцов) цифрового массива ФЦИ 
. Математическое выражение, используемое для расчетов нормированных
коэффициентов взаимной корреляции для соседних строк (столбцов) цифрового
массива изображения, приведено в работе [3, с 132]. Эта формула используется
для обработки экспериментальных данных при изучении пространственных
корреляционных связей флуктуаций излучения атмосферных фонов. Принято считать:
если коэффициент взаимной корреляции второй диагонали корреляционной матрицы
стремится к единице, то корреляционные связи сильны или устойчивы, а если
меньше 0,5, то корреляционные связи слабы или отсутствуют [3].
Рис. 3. Структурная схема
теплопеленгатора, построенного с использованием способа корреляционной компенсации фона
В ходе исследования
корреляционных связей фонового излучения облачной атмосферы установлено, что
нормированные коэффициенты взаимной корреляции соседних строк, столбцов массива
изображения, находящиеся во вторых диагоналях корреляционных матриц, принимают
значения в интервале [0,65-0,99]. При проведении
предварительного математического моделирования по обнаружению теплового объекта
на атмосферном фоне установлено, что присутствие теплового объекта в элементе 
массива изображения 
приведет к ослаблению
корреляционных связей. При этом нормированные коэффициенты взаимной корреляции
соседних строк, столбцов массива изображения, в которых находится изображение
цели, находящиеся во вторых диагоналях корреляционных матриц, снижаются по
сравнению со значениями для строк, столбцов, в которых его нет. Это явление
предлагается считать признаком строк и столбцов, в которых находится отметка теплового
объекта.
Основой алгоритма
пространственной фильтрации является расчет диагональной корреляционной матрицы
по строкам и матрицы 
по столбцам для
массива 
ФЦИ. Первая
(центральная) диагональ каждой корреляционной матрицы содержит элементы, равные
единице, поскольку они получены путем вычисления коэффициентов корреляции
строки (столбца) с самим собой. Вторая диагональ матрицы показывает, каков коэффициент
взаимной корреляции первой строки со второй, второй с третьей и т. д. Из элементов
вторых диагоналей корреляционных диагональных матриц 
и 
получают векторы 
и 
. Вектора 
, 
используются для
формирования массивов 
, в каждый элемент, (имеющий индексы n и m) которого записывался результат умножения в соответствии с формулой
, где 
и 
- величины n-го и m-го
элементов векторов 
и 
соответственно. По
адаптивному порогу 
, в соответствии с решающим правилом в элемент 
записывается 0, если в соответствующем элементе массива
. В элемент 
записывается 1, если в соответствующем элементе массива
. Таким образом, окончательно формируется массив 
. Массив 
используется для
создания бинарного изображения на мониторе ВКУ. По номерам строк и столбцов элементов,
в которых находятся единичные значения, определяется пространственное положение
цели в сегменте полусферы поиска. Изображение на мониторе с отфильтрованной
точечной тепловой целью представлено на рис. 2.
С целью проверки
эффективности предлагаемого алгоритма обнаружения теплового объекта было
проведено математическое моделирование с использованием моделей излучения
воздушных целей и фона. Модели построены на основе измерений собственного излучения
воздушных целей и атмосферных фонов при различных метеорологических ситуациях.
Измерения собственного излучения воздушных целей и атмосферных фон проводились
с помощью измерительно-вычислительного комплекса (ИВК), в состав которого входит
двухспектральный радиометр (3–5 мкм и 8–13 мкм). Характеристики радиометра приведены
в работе [3].
Модели фона представляли
собой цифровые массивы изображений
, в каждом элементе 
которых записаны
энергетические яркости собственного излучения атмосферного фона. Модели
излучения воздушной цели представлены средними контрастами излучения воздушных
целей, полученных на различных по типу облачности, фонах и дальности [5]. Средний
контраст излучения тепловой воздушной цели рассчитывался по формуле
,
(1)
где 
- среднее значение собственного излучения
тепловой цели; 
- среднее значение собственного излучения
фона, на котором наблюдается тепловая цель.
Для проведения
математического моделирования обнаружения теплового объекта на атмосферном фоне
из имеющейся базы моделей фона были отобраны около 400 массивов изображений 
, полученных в спектральных диапазонах 3–5 мкм и 8–13 мкм,
размером 512×1024 пикселей. Массивы изображения 
охватывают верхнюю
полусферу от горизонта до 40° и имеют больший размер по азимуту (80о).
В набор массивов изображений 
включены результаты
дневных и ночных измерений в летний и осенний сезоны для следующих типов
облачности: ясно, кучевая, слоистая, перистая.
С целью определения
вероятности обнаружения, реализуемой при использовании предлагаемого способа,
проведено моделирование процесса обнаружения теплового объектов на исследуемых
фоновых кадрах. Были взяты по 100 фрагментов массивов изображения 
для ясного неба,
кучевой, слоистой и перистой типов облачности, размером 
пикселей.
В качестве полезного сигнала
(изображение теплового объекта) использовалась величина среднего контраста
излучения вертолета, который вычислен по формуле (1). Размер изображения
теплового объекта умещался в одном элемент 
(пикселе) изображения атмосферного фона 
, что соответствует при линейном размере цели 4 м удалению 12
- 15 км. На выбранные массивы изображения
атмосферного фона 
«накладывался» сигнал
от тепловых объектов. «Наложение» т. е. имитация появления изображения
воздушной цели в выбранном пикселе 
массива проводилось в
соответствии с формулой
, (2)
где 
- средний контраст излучения
цели и фона. [5]. В результате моделирования, изображение теплового объекта
наблюдалось, затеняя участок атмосферного фона находящегося за ним. При каждом
«наложении» в случайно выбранный пиксель 
массива изображения
фона 
, создавалось отношение сигнал/шум
,
где 
- дисперсия флуктуаций
излучения атмосферного фона примерно равная дисперсии шума аппаратуры
(радиометра) 
≈ 4 · 10-11 (Вт∙ см-2
∙ ср-1)2 [2].
Обнаружение осуществлялось
по адаптивному порогу p для коэффициента взаимной
корреляции по строкам и столбцам массивов изображений 
. Диапазон изменения адаптивного порога находится в интервале
[0,1 – 0,6]. Значение порога зависит от типа атмосферного фона, которому
присуще свое значение радиуса корреляции, и от величины контраста цели на фоне
атмосферных помех.
Анализ результатов математического
моделирования подтвердил работоспособность предлагаемого способа пространственной
фильтрации в обоих инфракрасных диапазонах. Вероятности обнаружения теплового
объекта на атмосферном фоне при отношении сигнал/шум, равном 2, представлены в
таблице 1.
Таблица
1.
Вероятности обнаружения теплового объекта на коррелированном
атмосферном фоне
|
Диапазон длин волн |
Тип облачности |
|||
|
Ясно |
Кучевая |
Слоистая |
Перистая |
|
|
3 - 5 мкм |
0,82 |
0,67 |
0,73 |
0,82 |
|
8 - 13 мкм |
0,91 |
0,73 |
0,84 |
0,87 |
Как видно из сравнения
видеоизображений, показанных на рисунках 1 и 2, при использовании предлагаемого
способа отметка теплового объекта (целей), ввиду отсутствия фона на бинарном
изображении, видна контрастно. Поэтому оценить их пространственное положение в
сегменте полусферы поиска легко. Полученные на экране монитора бинарные
изображения тепловых объектов позволяют определять их угловые координаты,
количественный состав и другие характеристики без предварительной подготовки
оператора, что приводит к увеличению информационной способности и простоте
эксплуатации теплопеленгатора.
ЛИТЕРАТУРА
1. Якушенков Ю.Г.,
Тарасов В. В. Инфракрасные системы «смотрящего» типа. - М.: «Логос»,
2004. 430 с.
2. Приходько В.Н.,
Хисамов Р.Ш. Обнаружение «точечных» объектов теплопеленгатором на основе
матричного фотоприёмного устройства. //Оборонная техника.// Вып. 1–2, 2007. С.
64-66.
3. Алленов М.И. и
др. Стохастическая структура излучения облачности. - СПб.: Гидрометеоиздат,
2000. 175 с.
4. Чупраков А. М.,
Хитрик А.С. Тепловизионный прицел на основе матричного болометрического
приемника. //Оптико–электронные системы визуализации и обработки оптических
изображений.// Вып. 2. - М.: ЦНИИ «Циклон».
2007. С. 60-71.
5. Левшин В.Л. Пространственная фильтрация в оптических
системах пеленгации. - М.: «Советское
радио», 1971. 199 с.
DETECTION of thermal objects on the
correlated atmospheric background
Jakimenko I.V., Zhendarev M.V.,
Gurthenkov D.A.
In this abstract is stated the way of the thermal objects detection on
the correlated atmospheric background by means of infra-red heat direction
finders with a wide field of vision. The way is based on use of differences of
spatial spectra of dot thermal object radiation and extensional colder atmospheric
background. On the basis of this way is developed the algorithm of the bidimentional
spatial filter which can be realized programmatically in the video processor of
the heat direction finder. The algorithm allows automating process of the air
targets detection (planes, helicopters, cruise missiles, unpiloted aircraft,
etc.) and to provide the display of their video images to the monitor of the
video control device in the binary kind.
Key words: the thermal objects detection, infra-red optiko-electronic systems
Военная академия войсковой
ПВО
Вооруженных Сил Российской Федерации
имени Маршала Советского Союза А.М. Василевского
(ВА ВПВО ВС РФ)
Поступила в редакцию
16.03.2009.