Математическая морфология.

Электронный математический и медико-биологический журнал. - Т. 8. -

Вып. 2. - 2009. - URL:

http://www.smolensk.ru/user/sgma/MMORPH/TITL.HTМ

http://www.smolensk.ru/user/sgma/MMORPH/N-22-html/TITL-22.htm

http://www.smolensk.ru/user/sgma/MMORPH/N-22-html/cont.htm

 

 

УДК 621.391.161

 

СПОСОБ  ЧЕРЕСПЕРИОДНОЙ  КОМПЕНСАЦИИ  АТМОСФЕРНОГО ФОНА ПРИ  ПРОСТРАНСТВЕННОЙ  ФИЛЬТРАЦИИ  ТЕПЛОВЫХ  ОБЪЕКТОВ

 

Ó 2009 г. Якименко И. В., Гурченков Д. А.

 

(jakimenko-1.doc)

 

В данной работе изложен способ обнаружения тепловых объектов на атмосферном фоне с помощью инфракрасных теплопеленгаторов с широким полем зрения. Предлагаемый способ основан на использовании отличий в амплитудах яркости излучения точечного теплового объекта и холодного, атмосферного фона. На основе данного  способа обнаружения тепловых объектов разработан алгоритм двумерного пространственного фильтра, который может быть реализован программно в видеопроцессоре теплопеленгатора. Алгоритм позволяет  обнаруживать воздушные цели (самолеты, вертолеты, крылатые ракеты, беспилотные летательные аппараты и т. д.) на максимальной дальности  и обеспечивает вывод их видеоизображений на монитор видеоконтрольного устройства в бинарном виде.

Ключевые слова: обнаружение тепловых объектов.

 

Рассматривается способ обнаружение точечных тепловых объектов (ТТО). на естественных атмосферных фонах (АФ), что соответствует условиям обнаружения малоразмерного объекта на больших дальностях. Под ТТО понимается малоразмерный объект, изображение которого вписывается в элементарное поле зрения (пиксель) пеленгатора.

В настоящее время для теплопеленгаторов применяется амплитудный способ обнаружения ТТО на фоноцелевом изображении (ФЦИ). Данный способ основан на зрительном разделении яркости отметок ТТО и флюктуаций атмосферного фона (АФ). Получение ФЦИ на мониторе видеоконтрольного устройства (ВКУ) проводится в соответствии с видеосигналом, полученным с выхода МПИ. Каждому пикселю ФЦИ в соответствие ставят ограниченное число дискретных уровней яркости – квантов, пропорциональных энергетической яркости собственного излучения атмосферы и цели. Чем больше уровней квантования, тем выше качество изображения.

Особенностью работы ОЭС с широким полем зрения, работающих в диапазоне 8-13 мкм, является высокие перепады яркости одновременно наблюдаемых фрагментов изображения. Например, радиационный контраст кучевого облака на фоне чистого неба как днём, так и ночью, может достигать 60 К [3]. Аналогичная ситуация возникает при визировании цели вблизи линии горизонта, когда в угловом поле прицела наблюдается часть небесной полусферы и наземный ландшафт, который обычно «теплее» небесного. При чувствительности матрицы  <0,1 К и наличии ограничения в выходном видеосигнале «уровнем белого» осуществляется компрессия сигналов с МПИ. Это, в свою очередь, приводит к уменьшению сигнала цели на входе ВКУ [4]. ФЦИ, полученное на мониторе ВКУ теплопеленгатора при наличие цели (отметка от цели находится в разрыве линии), представлено на рис. 1. Изображение получено в темное время суток при кучевой облачности.

Проблемы видения возникают вследствие потенциально слабой контрастности ФЦИ, низкого уровня заметности целей и возможного появления пассивных и активных видов противодействия (маскировка, ложные тепловые цели (ЛТЦ), прожекторы и т. п.). Шумовые компоненты сигнала, образующиеся вследствие особенностей применяемых датчиков и искажении сигнала передачи данных приводят к ухудшению изображения ТТО на мониторе ВКУ. С целью улучшения качества изображений были созданы различные способы цифровой фильтрации и обработки двухмерного изображения. Обработка изображений проводится в реальном масштабе времени, без прореживания пиксельного потока видеоданных и без пропуска кадров. Основными способами обработки изображений являются алгоритмы цифровой коррекции неоднородностей, сглаживание шумов, повышение контрастов, подчеркивание границ, высокочастотная коррекция, коррекция дефектных пикселей и т. д., реализованных программно в видеопроцессоре теплопеленгатора. Отметим, что перечисленные способы обработки изображений косвенно улучшают характеристики теплопеленгаторов за счет улучшения визуального восприятия изображения оператором.

 

 

Рис. 1. Исходное изображение, полученное на мониторе ВКУ

 

Окончательной оценкой эффективности улучшения качества видеоизображений должно служить субъективно оцениваемое качество изображения, поскольку конечным потребителем воспроизводимого изображения является оператор. Зрительные возможности у разных операторов различны, и на неравномерном изображении облачного АФ трудно выделить ТТО. Ещё сложнее селектировать несколько отметок, поэтому при использовании зрительного обнаружения повышается вероятность ошибок принятия решений.

Таким образом, необходимы новые способы обнаружения воздушных целей, которые способны выделять отметки ТТО на изображении АФ без участия оператора. На конечном этапе обнаружения оператор на мониторе ВКУ должен увидеть изображение, для которого характерно максимально контрастное представление отметок ТТО на скомпенсированном равномерном АФ. Для получения такого изображения достаточно двух уровней квантования. Полученные на экране монитора бинарные изображения ТТО позволяют определять их угловые координаты и количественный состав (рис. 2).

Такое выделение ТТО на изображении АФ можно отнести к пространственной фильтрации. Для компенсации мешающих излучений предлагается использовать свойство излучения АФ, который имеет примерно одинаковый уровень флуктуаций яркости в пределах сегмента полусферы с угловыми размерами 1÷3о [3]. Размер поля зрения теплопеленгатора составляет 15о по углу места и 20о  по  азимуту. Размер  массива  череспериодной  компенсации  (МЧПК)  составляет  240 320 элементов, т. е. при переводе его в визуальное изображение оно будет содержать 240 320 = 76800 пикселей [4]. Это позволяет выделить 48 сегментов углового поля зрения прибора размером 2,5 2,5о, что соответствует делению МЧПК на фрагменты размером 40 × 40 элементов. Таким образом, МЧПК может быть разделён на 48 или менее фрагментов, в пределах которых существует одинаковое число облучённых элементов положительной и отрицательной полярности. Хаотическое размещение в сегменте обнаружения областей с различной интенсивностью излучения фона, а также их различная пространственная протяжённость не оказывают существенного влияния на результат усреднения излучения в том или ином столбце (строке). Неполная компенсация излучения фона может быть обусловлена нечетным числом элементов, принадлежащих частной области излучения фона, отличающейся по уровню излучения от других частных областей. Тем не менее, при попадании в сегмент углового поля зрения теплопеленгатора протяженной помехи в любой момент времени, число облучённых элементов по строкам и столбцам МЧПК, дающих сигнал положительной полярности, будет примерно равно числу элементов, дающих сигнал отрицательной полярности. При сложении амплитуд разнополярных сигналов, снятых по строкам или столбцам матрицы МЧПК, будет происходить их взаимосокращение. Величина среднего арифметического по строкам или столбцам, содержащим сигналы от АФ, будет близка к нулю. При усреднении разнополярных значений элементарных сигналов по строке (столбцу) МЧПК, содержащей амплитуду сигнала от ТТО, взаимной компенсации происходить не будет, т.е. среднее арифметическое будет значительно отличаться от нуля.

 

Рис. 2. Изображение на мониторе ВКУ, полученное после фильтрации

 

Таким образом, существует возможность, установив порог принятия решения, близкий к среднему арифметическому, рассчитанному по строкам и столбцам МЧПК, различать строки и столбцы, в которых находится изображение фона, от строк и столбцов, в которых записана амплитуда сигнала от ТТО, а также определять их номера. В этом состоит сущность способа череспериодной компенсации фона.

Порог устанавливается каждый раз перед началом обработки массива изображения. Для этого используется результат расчета среднего арифметического по всему массиву изображения. Таким образом, меняющийся порог на каждом этапе обработки для каждого массива, адаптивно учитывает средний уровень энергетической яркости холодного протяжённого АФ. Это возможно, поскольку элементов с отражениями фона на несколько порядков больше элементов с отражением от ТТО. При равномерном фоне в сегменте обнаружения порог по модулю уменьшается. В случае наличия участков с разной яркостью уровень порога p растёт. То есть величина порога адаптивно реагирует на условия обнаружения ТТО. Большое число строк и столбцов МЧПК и малое число неоднородных областей фона (более или менее тёмных) приводят к тому, что при усреднении по числу элементов МЧПК величина порога p изменяется несущественно, однако адаптивно реагирует на наличие неоднородностей в анализируемом сегменте небесной полусферы. Чтобы величина порога p была менее чувствительна к наличию неоднородных (разнояркостных) областей, существующих в пределах сегмента обнаружения, необходимо увеличивать число элементов МЧПК за счёт применения более чувствительных приёмников с большей дискретностью опроса.

Пространственный фильтр предлагается реализовать программно в видеопроцессоре. Структурная схема теплопеленгатора, реализующая предлагаемую пространственную фильтрацию, представлена на рис. 3. 

 

 

Рис. 3. Структурная схема теплопеленгатора, построенного с использованием 
способа череспериодной компенсации фона

 

С целью проверки эффективности предлагаемого алгоритма двумерной пространственной фильтрации ТТО было проведено математическое моделирование с использованием моделей излучения воздушных целей и фона. Модели построены на основе оцифрованных, при помощи специальной программы [6], кадров видеосигнала, полученных с видеовыхода  теплопеленгатора.  Примерный вид кадра видеосигнала, из которого получили массив представлен на рисунке 1. Каждый элемент  массива  содержит информацию о дискретных уровнях квантования, пропорциональных яркости излучения АФ, снятого с ячейки многоэлементного приемника в n-й строке на m-м шаге.  Сигнал в каждом пикселе квантован в 16777162 уровней яркости и занимает 24 бита. Характеристики теплопеленгатора приведены в работе [7].

Модели фона представляли собой цифровые массивы изображений , в каждом элементе  которых записаны энергетические яркости собственного излучения атмосферного фона без ТТО. Модели излучения ТТО представлены средними контрастами излучения воздушных целей, полученных на различных по типу облачности, фонах и дальности [1]. Средний контраст излучения ТТО рассчитывался по формуле

,                                                                (4)

где  - среднее значение собственного излучения тепловой цели;  - среднее значение собственного излучения фона, на котором наблюдается тепловая цель.

Для проведения математического моделирования обнаружения теплового объекта на атмосферном фоне из имеющейся базы кадров фона были отобраны около 400 массивов изображений , полученных в спектральных диапазонах 8–13 мкм, размером 320×240 пикселей. В набор массивов кадров  включены результаты ночных измерений в летний и осенний сезоны для следующих типов облачности: ясно, кучевая, слоистая, перистая.

С целью определения вероятности обнаружения, реализуемой при использовании предлагаемого способа, проведено моделирование процесса обнаружения теплового объекта на исследуемых фоновых кадрах. Были взяты по 100 кадров  для названных типов облачности.

В качестве полезного сигнала (изображение теплового объекта) использовалась величина среднего контраста излучения вертолета, который вычислен по формуле (4). Размер изображения ТТО умещался в одном элемент (пиксель) изображения атмосферного фона , что соответствует при линейном размере цели 4 м удалению 12 - 15 км. На выбранные массивы изображения атмосферного фона  «накладывался» сигнал от тепловых объектов. «Наложение», т. е. имитация появления изображения воздушной цели в выбранном пикселе  массива, проводилось в соответствии с формулой

,

где - средний контраст излучения цели и фона. [2]. В результате моделирования, изображение теплового объекта наблюдалось, затеняя участок атмосферного фона находящегося за ним. При каждом «наложении» в случайно выбранный пиксель  массива изображения фона  создавалось отношение сигнал/шум

,

где  - пиковое значение квантованной яркости в  - дисперсия флуктуаций излучения атмосферного фона [5].

Математическое моделирование включало предварительные вычисления по расчёту среднего арифметического значения  яркости фонового шума массива изображения  по формуле

,

где -  элемент массива , содержащий информацию о сигнале, пропорциональном яркости излучения фона, снятого с ячейки многоэлементного приёмника в n-й строке на m-м шаге; N - чётное число строк; M - чётное число столбцов массива.

Модуль адаптивного порогового значения по строкам  рассчитывается по формуле

,

где k1 - весовой коэффициент, позволяющий адаптировать порог принятия решения по строкам.

Модуль адаптивного порогового значения по столбцам рассчитывается по формуле

,

где k2 - весовой коэффициент, позволяющий адаптировать порог принятия решения по столбцам.

Основная часть моделирования включала формирование МЧПК , в каждый элемент (имеющий индексы n и m) которого записывали результат умножения в соответствии с формулой

.

Элементы массива изображения , обязательно должны иметь четное количество N строк и четное число M столбцов.

Для массивов  были рассчитаны средние значения  по столбцам по формуле

.

Затем значения  заносились в соответствующие элементы вектора .

В результате пороговой обработки, при которой сравнивались величины элементов  вектора  с порогом принятия решения , был получен выходной вектора  в соответствии с правилом:

при попадании величины элемента  вектора  в интервал принятия решения, величину элемента обнуляют: ;

если значение элемента  вектора  не попало в интервал принятия решения, ему присваивают единичное значение .

Далее для массивов  были рассчитаны средние значения  по столбцам по формуле

.

Затем значения  заносились в соответствующие элементы вектора .

В результате пороговой обработки, при которой сравнивались величины элементов  вектора  с порогом принятия решения , был получен выходной вектор  в соответствии с правилом:

при попадании величины элемента  вектора  в интервал принятия решения, величину элемента обнуляют: ;

если значение элемента  вектора  не попало в интервал принятия решения , ему присваивают единичное значение .

Полученные векторы ,  были использованы для получения массивов , в каждый элемент (имеющий индексы n и m) которого записывался результат умножения в соответствии с формулой

,

где  и  - величины n-го и m-го элементов векторов  и  после пороговой обработки соответственно.

Массив  был использован для создания фоноцелевого изображения на мониторе видеоприемного устройства. На мониторе формируется бинарное изображение фоноцелевой обстановки, где местоположение цели в сегменте полусферы обозначается на черном фоне отметкой белого цвета. Бинарное изображение фоноцелевой обстановки представлено на рис. 2.

При попадании четного количества точечных отметок нескольких ТТО в одну строку или столбец возможен пропуск цели. Однако после завершения первого этапа обнаружения процесс селекции цели на фоне неба не прекращается, а повторяется многократно. Неизменное положение ТТО, в роли которого выступает воздушная цель, в одном и том же пикселе изображения является маловероятным событием. Поэтому пропуск цели на одном из этапов, за счет случайного попадания отметок ТТО в одну строку или столбец, не приводит к полной потере цели.

Анализ результатов математического моделирования подтвердил работоспособность предлагаемого способа пространственной фильтрации в инфракрасном диапазоне 8-13 мкм. Вероятности обнаружения теплового объекта на атмосферном фоне при отношении сигнал/шум, равном 2, представлены в таблице 1.

 

Таблица 1.

Вероятности обнаружения теплового объекта на атмосферном фоне

 

Диапазон

длин волн

Тип облачности

Ясно

Кучевая

Слоистая

Перистая

8 - 13 мкм

0,88

0,82

0,85

0,89

 

Как видно из сравнения видеоизображений, приведенных на рисунках 1 и 2, при использовании предлагаемого способа отметка ТТО (целей), ввиду отсутствия АФ на бинарном изображении, видна контрастно. Следовательно оценить их пространственное положение в сегменте полусферы поиска легко. Полученные на экране монитора бинарные изображения тепловых объектов позволяют определять их угловые координаты, количественный состав и другие характеристики без предварительной подготовки оператора, что приводит к увеличению информационной способности и простоте эксплуатации теплопеленгатора.

 

Литература

 

1. Якушенков Ю.Г., Тарасов В. В. Инфракрасные системы «смотрящего» типа. - М.: «Логос», 2004. 430 с.

2. Приходько В.Н., Хисамов Р.Ш. Обнаружение «точечных» объектов теплопеленгатором на основе матричного фотоприёмного устройства. //Оборонная техника.// Вып. 1–2, 2007. С. 64-66.

3. Алленов М.И. и др. Стохастическая структура излучения облачности. - СПб.: Гидрометеоиздат, 2000. 175 с.

4. Чупраков А. М., Хитрик А.С. Тепловизионный прицел на основе матричного болометрического приемника. //Оптико–электронные системы визуализации и обработки оптических изображений.// Вып. 2. - М.: ЦНИИ «Циклон». 2007. С. 60-71.

5. Левшин В.Л.  Пространственная  фильтрация  в  оптических  системах пеленгации. - М.: «Советское радио», 1971. 199 с.

6. Якименко И.В., Коваль С.Н., и др. Цифровая обработка сигналов тепловизионных устройств перспективных образцов вооружений // Программа ЭВМ. Зарегистрирована ФГУП «Всероссийский научно-технический центр» Инв. Номер 50200900390

7. Чупраков А. М., Хитрик А.С. Тепловизионный прицел на основе матричного болометрического приемника. //Оптико–электронные системы визуализации и обработки оптических изображений.// Вып. 2. - М.: ЦНИИ «Циклон». 2007. С. 60-71.

 

Way of a spatial filtration of the THERMAL

OBJECTS on the atmospheric background

 

Jakimenko I. V., Gurchenkov D. A.

 

In this abstract is stated the way of the thermal objects detection on the 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.

Кеу words: the thermal objects detection

 

 

 

Академия войсковой ПВО Вооруженных Сил РФ

им. Маршала Советского Союза А. М. Василевского, г. Смоленск

Поступила в редакцию 11.04.2009.