автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Разработка и исследование алгоритмов обнаружения протяженных аномалий на многозональных изображениях

На правах рукописи

005015567

Лучков Николай Владимирович

РАЗРАБОТКА И ИССЛЕДОВАНИЕ АЛГОРИТМОВ ОБНАРУЖЕНИЯ ПРОТЯЖЕННЫХ АНОМАЛИЙ НА МНОГОЗОНАЛЬНЫХ ИЗОБРАЖЕНИЯХ

Специальность: 05.13.18 - «Математическое моделирование, численные методы и комплексы программ»

1 МАР 2П12

Автореферат

диссертации на соискание ученой степени кандидата технических наук

Ульяновск-2012

005015567

Работа выполнена на кафедре «Телекоммуникации» Ульяновского государственного технического университета.

Научный руководитель - доктор технических наук, профессор

Васильев Константин Константинович

Официальные оппоненты - доктор физико-математических наук, профессор,

зав. кафедрой «Прикладная математика и информатика» УлГТУ Валеев Султан Галимзянович

- кандидат технических наук, нач. отдела

информационного и технического обеспечения филиала ОАО «БИНБАНК» в Ульяновске Герчес Владислав Геннадьевич

Ведущая организация - Ульяновский государственный университет

Защита диссертации состоится « 28 » марта 2012 г. в 12 часов 00 минут на заседании диссертационного совета Д 212.277.02 при Ульяновском государственном техническом университете по адресу: 432027, г. Ульяновск, ул. Северный Венец, 32 (ауд. 211, Главный корпус).

С диссертацией можно ознакомиться в библиотеке Ульяновского государственного технического университета.

Автореферат разослан « » февраля 2012 г.

Ученый секретарь диссертационного совета, доктор технических наук, профессор

В.Р. Крашенинников

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы

В связи с широким распространением методов мультиспектральной и гиперспектральной регистрации земной поверхности, все большее значение приобретают задачи вероятностного представления и статистического анализа многомерных данных. Совместная обработка данных в нескольких спектральных диапазонах позволяет эффективно решать задачи радио- и гидролокационных комплексов, систем глобального мониторинга Земли, систем технического зрения и др. Для ряда приложений особый интерес представляют задачи обнаружения и оценивания параметров протяженных аномалий на многозональных изображениях.

В ряде работ предлагаются локальные алгоритмы обнаружения на многомерных изображениях, основанные на оценках градиента поля изображения. Рассматривается применение адаптивных алгоритмов компенсации коррелированных помех и алгоритмов обнаружения с «выбеливанием» многомерного изображения. Несмотря на большое число публикаций по проблемам синтеза алгоритмов обнаружения сигналов, в настоящее время мало исследована их эффективность. Существующие решения не дают удовлетворительного результата в задачах анализа эффективности и сравнения соответствующих алгоритмов обнаружения, особенно алгоритмов обнаружения протяженных аномалий на многомерных изображениях.

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

Цель и задачи исследования

Целью работы является повышение эффективности систем мониторинга земной поверхности на основе разработки, моделирования и анализа вероятностных характеристик алгоритмов обнаружения протяженных неоднородностей на многозональных изображениях.

Для достижения этой цели необходимо решение следующих задач:

- анализ и разработка моделей формирования многозональных изображений, позволяющих при относительно небольших вычислительных затратах формировать изображения и их последовательности, имитирующие реальные многозональные изображения;

- синтез и анализ квазиоптимальных алгоритмов обнаружения протяженных аномалий на фоне мешающих многозональных изображений в условиях априорной неопределенности относительно уровней яркости сигнала;

- проведение сравнительного анализа эффективности известных и синтезированных квазиоптимальных алгоритмов обнаружения протяженных аномалий на фоне мешающих многозональных изображений;

-разработка пакета прикладных программ для анализа эффективности синтезированных обнаружителей на основе статистического моделирования, а также при использовании кадров реальных спутниковых изображений.

Методы исследований. При решении поставленных задач были использованы методы теории вероятностей, математической статистики, математического моделирования, численной аппроксимации функций, средства современной информатики и вычислительной техники, включая язык программирования Delphi 7 и программные средства работы с графическими изображениями Paint.NET и MultiSpec.

Научная новизна результатов

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

2. Синтезированы квазиоптимальные алгоритмы обнаружения протяженных аномалий на фоне многозональных изображений при наличии случайных помех в условиях неизвестных параметров полезного сигнала.

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

Практическая ценность диссертационной работы заключается в следующем:

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

2. Разработанный пакет программ допускает его непосредственное использование при проектировании современных и перспективных систем моделирования и обработки изображений. Предложенные модели многомерных СП могут быть использованы разработчиками для представления реальных данных при решении различных прикладных задач.

Апробация работы.

Основные результаты диссертационной работы докладывались и обсуждались на следующих научно-технических конференциях: 63 Научная сессия, посвященная Дню радио (Москва, 2008); Седьмая Международная конференция «Математическое моделирование физических, экономических, технических, социальных систем и процессов» (Ульяновск, 2009); Шестая Всероссийская научно-практическая конференция «Современные проблемы создания и эксплуатации радиотехнических систем» (Ульяновск, 2009); 13 Международная конференция по распознаванию образов и анализу

изображений (Санкт-Петербург, PRIA, 2010); Научно-техническая конференция «Интегрированные автоматизированные системы управления» (Ульяновск, ФНПЦ ОАО «НПО «МАРС», 2011); ежегодные научно-технические конференции профессорско-преподавательского состава Ульяновского государственного технического университета (2009-2012).

Реализация результатов работы.

Результаты работы использованы ФНПЦ ОАО «НПО «Марс». Реализация результатов работы подтверждена соответствующим актом.

Публикации. Основные научные результаты диссертационной работы отражены в 14 публикациях, среди которых 5 статей, из них 2 статьи в журналах из перечня ВАК, 5 работ в трудах и материалах международных и Всероссийских конференций. По результатам работы подана одна заявка на изобретение.

Структура диссертационной работы: Диссертационная работа написана на русском языке и состоит из введения, четырех глав, заключения, списка использованной литературы и приложения. Работа изложена на 121 страницах, содержит 67 рисунков и список литературы из 151 наименования.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обоснована актуальность работы, сформулированы цели и основные задачи, показана научная новизна и практическая ценность работы.

В первой главе рассмотрены математические модели многозональных изображений, а также предложена методика моделирования многомерных изображений. Дается краткий обзор известных алгоритмов предварительной обработки спутникового материала и существующих методов обнаружения аномалий на изображениях. Формулируются основные нерешенные задачи, определившие тему диссертационной работы.

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

(базовое СП). Далее значения СП в полученном кадре р преобразуются в

совокупность корреляционных параметров {/?,, 7е(.Л >—>Л/)}> где М -

размерность формируемого изображения. Затем формируется изображение как модель СП с изменяющимися корреляционными параметрами р-. Для

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

Рис. 1,а Рис. 1,6

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

размеры и формы объектов на моделируемом изображении, увеличение этих параметров приведет к увеличению размеров объектов на моделируемом изображении.

Для решения задачи идентификации параметров дважды стохастических авторегрессионных моделей на основе реальных данных проведены аналитические исследования вероятностных характеристик этих моделей.

Допустим необходимо подобрать параметры дважды стохастической модели СП так, чтобы она имела характеристики реального изображения с

полем значений яркости / = 1,...,ЛГ; /=1,...,Л". Опишем его с помощью авторегрессионной модели Хабиби:

2,1 = + - РщРяА,-„о-,) + '

= Vя*' / = и.,М,; 7=1 ,...,М2.

Определим значение статистических параметров изображения |г?| в

каждой точке. Для этого выберем размер скользящего окна Л'хЛ' так, чтобы оно не превышало размеры однородных областей изображения. Далее в скользящем окне оценим среднестатистическое значение коэффициента корреляции для каждого пикселя по строке рХ!] и по столбцу р ., среднестатистические значения

математического ожидания т:1/ и дисперсии о~1;.

Рис. 2, а

Рис. 2, б

Рис. 2. в

Рис. 2, г

Таким образом, формируем СП корреляционных параметров и

||, СП значений и СП значений {с^}. на основе которых можно

имитировать изображения с изменяющимися корреляционными параметрами.

На рис. 2 представлен пример использования предложенной методики имитации реального спутникового снимка: а - спутниковый снимок размеров 300*300 пикселей; б - дважды стохастическая модель изображения; в - СП значений дисперсии; г - поле корреляционных параметров. Статистическая оценка параметров спутникового снимка проводилась в скользящем окне размером 10x10.

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

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

Рассмотрение проблем обнаружения аномалий привело к необходимости определения ковариационных матриц ошибок фильтрации. Чтобы найти необходимые для решения задачи обнаружения оптимальные оценки СП и ковариационную матрицу ошибок фильтрации, необходимо решить систему линейных уравнений Винера-Хопфа относительно вектора а , составленного из весовых коэффициентов оптимального оценивания.

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

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

Вычислительная сложность предложенного метода расчета ковариационной матрицы ошибок фильтрации СП двухмерного изображения по сравнению с прямым методом, приводит к уменьшению объема вычислений в М1 раз. Для СП многомерного изображения - в N раз, где N количество кадров в изображении.

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

Как показывает анализ, реальные возможности получения компактных соотношений для спектральной плотности связаны с полями, КФ которых

п

могут быть представлены в виде произведения i?(/C>) = (р, j или линейной

¡=1

комбинации таких произведений. При этом относительная дисперсия ошибки фильтрации многомерного СП найдена в виде Д-мерного интеграла.

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

Для оценки сходимости аппроксимирующего соотношения к точному результату были произведены численные расчеты дисперсии ошибки фильтрации СП непрерывного аргумента и дискретного для СП различной размерности. Вычисление интегралов производилось адаптивным методом Симпсона в среде программирования Delphi 7.

В результате анализа полученных результатов можно сделать вывод, что с увеличением размерности СП ошибка аппроксимации интегральными соотношениями для непрерывного аргумента увеличивается. На рис. 3 представлены

1-Р """" ~~

зависимости граничных

значений д и (1-р) для СП

размерности п=2 при

р{ = р2 = р, обеспечивающие

заданную относительную

погрешность аппроксимации нг

0 = \а2 -<х; |/<т2 .

| нспр ОНСКр \1 лепр

На основе результатов

10 |

расчетов представляется

Рис. 3

возможным определить

область значений ц и Р) обеспечивающую заданную относительную погрешность аппроксимации. Например, чтобы обеспечить ошибку аппроксимации не более 20 % для СП с коэффициентом корреляции р = 0.999, уровень сигнала не должен превышать 7, а для ограничения ошибки аппроксимации на уровне 5 % значение сигнал/шум не должен превышать 0.6.

Проведена оценка вычислительной сложности предложенного алгоритма вычисления дисперсии ошибки в сравнении с проведением статистического моделирования. При вычислении интеграла с точностью 10'6 и /V = М~3, количество операций умножения уменьшиться примерно в 10 раз для изображения размером 100*100, и в 10*раз для изображения 1000x1000.

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

построить характеристики оптимальных алгоритмов обнаружения аномалий на многозональных изображениях.

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

Существенным недостатком алгоритмов обнаружения сигналов с неизвестными уровнями является достаточно сложный анализ их эффективности. В процессе поиска путей упрощения такого анализа удалось получить еще одну процедуру обнаружения протяженных сигналов с неизвестными уровнями :

1 = Ы ~) [> А) ~сшиал есть'

а2в [< ¿0 - сигнала нет. где Хук - оптимальные оценки СП в области (;'./) = О'; Д = 1 .../V, сделанная на основе всех наблюдений |г|у,,| (оценка в «точку»). При этом оценки уровней

к = \....И, находятся из системы ¿V линейных уравнений.

Для определения эффективности предложенного алгоритма обнаружения протяженных аномалий с неизвестными уровнями яркости проведен расчет зависимости вероятности правильного обнаружения Рс1 от уровня полезного сигнала Б, обеспечивающего заданные значения вероятности ложной тревоги Р/=0.0()1. Рассмотрены случаи, когда полезный сигнал имеет одинаковые уровни значений на всех кадрах изображения (рис. 4, а) и различные уровни (рис. 4, б). При этом графики с номером 1 построены для оптимального алгоритма обнаружения, а графики с номером 2 - для алгоритма с предварительным оцениванием полезного сигнала. Использовались следующие параметры моделируемого СП: = 1;<тв2 = \;р = 0.9;рк = 0.9; размер сигнала 7x7.

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

0.3 1 Рис. 4, а

Реї

О 3 0.6 04 02 0

: ? .......

.........! / ! 1. ........;/-..... "УТ^Г........

¿'І--"" і

0.5

2.6 5

1 1.5

Рис. 4, б

правильного обнаружения увеличивается примерно на 60 % при уровне сигнала равном 1.5. При этом применение предложенного алгоритма приводит к увеличению операций умножения примерно на Л,-,2/У3, Л'5 - размер сигнала; Ы- количество кадров СП.

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

Анализ графиков, представленных на рис. 5, позволяет сделать вывод, что при увеличении изображения на один кадр и г = 4 вероятность правильного обнаружения увеличивается примерно на 10 %. Таким образом, представляется возможным выбрать необходимое количество кадров наблюдаемой сцены, обеспечивающее заданные характеристики эффективности обнаружения.

гнс. :>

1.2

0.4 0.2

0.8 0.6

10

Рис. 6

100

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

порогового полезного сигнала от внутрикадрового интервала корреляции СП и от размера полезного сигнала. Расчеты проводились при вероятности ложной тревоги ^=0.001 и вероятности правильного обнаружения Рв=0.5. На основании результатов проведенных расчетов можно сделать вывод, что при небольшой протяженности полезного сигнала значения порогового уровня обнаружения принимает максимальное значение при г « аМ5, а ~ 0,5. При увеличении значения № коэффициент а уменьшается и стремится к значению = 0,2.

Для уменьшения вычислительных затрат были произведены исследования зависимости вероятности правильного обнаружения Р0 от размера выбранного окна М, в котором производится обнаружение полезного сигнала с фиксированной протяженностью Ж Анализ полученных результатов позволяет сделать вывод, что при размерах окна, больших значения

дальнейшее увеличение размера практически перестает влиять на величину вероятности правильного обнаружения.

В четвертой главе представлены результаты экспериментальных исследований, иллюстрирующие эффективность синтезированных методов и алгоритмов при работе с реальными многозональными спутниковыми снимками. Изложены основные принципы организации разработанного программного обеспечения.

N=N5 +

Разработанный программный пакет предназначен для решения следующих задач: математическое моделирование отдельных изображений и их последовательностей; варьирование параметров моделей в широких пределах; тестирование алгоритмов фильтрации и обнаружения протяженных аномалий, как на реализациях моделей, так и на реальных изображениях; изучение сравнительных статистических характеристик алгоритмов фильтрации и обнаружения. Программный пакет состоит из следующих основных частей: программа моделирования СП; программа фильтрации СП; программа обнаружения аномалий на СП; пакет программ для исследования эффективности предложенных алгоритмов. Пакет программ реализован на объектно-ориентированном языке Delphi 7 и работает на платформе Windows.

На рис. 7 представлен пример эффективного использования спутниковых данных в организации дистанционного мониторинга лесных пожаров.

Рис. 7

Результаты работы алгоритма обнаружения, представленные на рис. 7, позволили выявить появление на спутниковых снимках четырех новых очагов пожаров. Причем анализ эффективности показал, что при нормированных значениях полезного сигнала s>0.03 вероятность правильного обнаружения близка к единице.

Полученные аналитические результаты были использованы в корабельных АСУ при решении задач автоматического обнаружения с помощью радио и гидролокационных комплексов.

В заключении сформулированы основные научные и практические результаты диссертационной работы.

Основными результатами являются следующие.

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

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

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

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

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

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

Основные результаты диссертации изложены в следующих публикациях

В изданиях из списка ВАК РФ:

1. Васильев, К. К. Эффективность цифрового обнаружения протяженных аномалий с неизвестными уровнями на фоне мешающих помех / К. К. Васильев, Н. В. Лучков// Цифровая обработка сигналов. -2011. -№3. - С. 24-27.

2. Васильев, К. К. Обнаружение протяженных аномалий с неизвестными уровнями на фоне мешающих помех / К. К. Васильев, Н. В. Лучков // Автоматизация процессов управления. -2011.- №4. - С. 3-7.

В других изданиях:

3. Васильев, К. К. Эффективность фильтрации случайных полей на многомерных сетках / К. К. Васильев, Н. В. Лучков // Вестник УлГТУ. - 2008. -№2.-С. 28-31.

4. Васильев, К. К. Эффективность фильтрации случайных полей на многомерных сетках / К. К. Васильев, Н. В. Лучков // Труды 63 научной сессии, посвященной Дню Радио. - М., 2008. - С. 378-380.

5. Лучков, Н. В. Моделирование обнаружения протяженных сигналов на фоне коррелированных помех / Н. В. Лучков // Математическое моделирование физических, экономических, технических, социальных систем и процессов: труды седьмой Международной конф. - Ульяновск, 2009. - С. 172-173.

6. Лучков, Н. В. Исследование эффективности алгоритмов обнаружения протяженных аномалий на фоне коррелированных помех / Н. В. Лучков // Современные проблемы создания и эксплуатации радиотехнических систем: труды шестой Всероссийской науч.-практ. конф. - Ульяновск, 2009. - С. 224-227.

7. Васильев, К. К. Анализ эффективности обнаружения протяженных сигналов на фоне коррелированных помех / К. К. Васильев, Н. В. Лучков // Автоматизация процессов управления. -2009. -№2. - С. 53-57.

8. Васильев, К. К. Моделирование алгоритма обнаружения протяженных сигналов на многозональных изображениях / К. К. Васильев, Н. В. Лучков // Автоматизация процессов управления. - 2009. — №4. - С. 85-88.

9. Лучков, Н. В. Анализ эффективности обнаружения протяженных сигналов на многомерных сетках / Н. В. Лучков // Современные проблемы проектирования, производства и эксплуатации радиотехнических систем: сборник науч. трудов. - Ульяновск : УлГТУ, 2010. Вып. 7. - С. 97-101.

10. Лучков, Н. В. Эффективность обнаружения протяженных сигналов различной формы / Н. В. Лучков // Современные проблемы проектирования, производства и эксплуатации радиотехнических систем : сборник науч. трудов. -Ульяновск: УлГТУ, 2010. Вып. 7. - С. 251-254.

11. Лучков, Н. В. Эффективность обнаружения протяженных сигналов на фоне коррелированных помех / Н. В. Лучков // Итоги диссертационных исследований: труды II Всероссийского конкурса молодых ученых. - М.: РАН, 2010.-С. 131-143.

12. Васильев, К. К. Анализ эффективности обнаружения протяженных сигналов на изображении морского дна / К. К. Васильев, Н. В. Лучков // Интегрированные автоматизированные системы управления : сборник докладов науч.-техн. конф. - Ульяновск: ФНПЦ ОАО «НПО «МАРС», 2011 - С. 100-103.

13. Васильев, К. К. Анализ эффективности фильтрации случайных полей на многомерных сетках / К. К. Васильев, В. Е. Дементьев, Н. В. Лучков // Физика и технические приложения волновых процессов : материалы Международной науч.-техн. конф. - Самара: 2011. - С. 346-348.

14. Vasil'ev, К. Detection effectiveness analysis for elongated signals on multidimensional grids / K. Vasil'ev, N. Luchkov, V. Dementjev // lOth International Conference on Pattern Recognition and Image analysys. - St. Petersburg, 2010. T.l. -P. 371-373.

ЛУЧКОВ Николай Владимирович

РАЗРАБОТКА И ИССЛЕДОВАНИЕ АЛГОРИТМОВ ОБНАРУЖЕНИЯ ПРОТЯЖЕННЫХ АНОМАЛИЙ НА МНОГОЗОНАЛЬНЫХ ИЗОБРАЖЕНИЯХ

Автореферат

Подписано в печать 13.02.2012. Формат 60x84/16. Усл. печ. л. 2.16. Тираж 100 экз. Заказ 163.

Типография УлГТУ, 432027, г. Ульяновск, Северный Венец, 32.

61 12-5/2275

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ федеральное государственное бюджетное образовательное учреждение высшего профессионального образования

«УЛЬЯНОВСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»

На правах рукописи

Лучков Николай Владимирович

РАЗРАБОТКА И ИССЛЕДОВАНИЕ АЛГОРИТМОВ ОБНАРУЖЕНИЯ ПРОТЯЖЕННЫХ АНОМАЛИЙ НА МНОГОЗОНАЛЬНЫХ

ИЗОБРАЖЕНИЯХ

Специальность: 05.13.18 - «Математическое моделирование, численные

методы и комплексы программ»

Диссертация на соискание ученой степени кандидата технических наук

Научный руководитель д. т. н., профессор К.К. Васильев

Ульяновск - 2012

ОГЛАВЛЕНИЕ

ВВЕДЕНИЕ..................................................................................................................4

ГЛАВА 1. МЕТОДЫ МОДЕЛИРОВАНИЯ И АНАЛИЗА ЭФФЕКТИВНОСТИ

АЛГОРИТМОВ ОБНАРУЖЕНИЯ АНОМАЛИЙ.................................................10

1.1 Постановка задачи...............................................................................................10

1.2. Обзор задач дистанционного зондирования Земли........................................11

1.3 Методы имитации изображений........................................................................18

1.4. Дважды стохастическая модель СП.................................................................25

1.5 Анализ характеристик дважды стохастической модели СП..........................30

1.6. Методика формирования изображений...........................................................37

1.7. Обзор алгоритмов обнаружение протяженных аномалий на

многозональных изображениях...............................................................................40

1.8 Выводы.................................................................................................................55

ГЛАВА 2. МОДЕЛИРОВАНИЕ И АНАЛИЗ АЛГОРИТМОВ ФИЛЬТРАЦИИ СЛУЧАЙНЫХ ПОЛЕЙ............................................................................................57

2.1. Постановка задачи..............................................................................................57

2.2. Эффективность фильтрации СП на многомерных сетках.............................58

2.3. Аппроксимация дисперсий ошибок фильтрации СП дискретного аргумента с помощью СП непрерывного аргумента...............................................................63

2.4. Определение условий приемлемой аппроксимации дисперсий ошибок фильтрации СП дискретного аргумента.................................................................68

2.5. Выводы................................................................................................................72

ГЛАВА 3. ОБНАРУЖЕНИЕ ПРОТЯЖЕННЫХ АНОМАЛИЙ НА ФОНЕ МНОГОМЕРНЫХ ИЗОБРАЖЕНИЙ......................................................................73

3.1. Постановка задачи..............................................................................................73

3.2. Исследование зависимостей вероятности правильного обнаружения от характеристик СП......................................................................................................74

3.3. Анализ эффективности обнаружения протяженных аномалий на многомерных СП.......................................................................................................78

3.4. Обнаружение протяженных аномалий с неизвестными уровнями...............82

3.5. Обнаружение появления аномалий на последовательности кадров многозональных изображений.................................................................................88

3.6. Эффективность обнаружения протяженных сигналов различной формы... 91

3.7. Выводы................................................................................................................94

ГЛАВА 4. ОСОБЕННОСТИ ПРОГРАММНОЙ РЕАЛИЗАЦИИ АЛГОРИТМА ОБНАРУЖЕНИЯ............................................................................95

4.1 Постановка задачи....................................... .........................................................95

4.2 Программный пакет для исследования алгоритмов обнаружения протяженных аномалий на фоне коррелированных помех..................................95

4.3 Пакет программ для анализа эффективности алгоритмов обнаружения протяженных аномалий на многомерных изображениях.....................................98

4.4 Применение предложенных алгоритмов........................................................101

4.5 Выводы...............................................................................................................109

ЗАКЛЮЧЕНИЕ.......................................................................................................110

СПИСОК ЛИТЕРАТУРЫ.......................................................................................111

ПРИЛОЖЕНИЕ.......................................................................................................123

ВВЕДЕНИЕ

Актуальность работы

В связи с широким распространением методов мультиспектральной (до 10 спектральных диапазонов) и гиперспектральной (до 300 диапазонов) регистрации земной поверхности, все большее значение приобретают задачи вероятностного представления и статистического анализа многомерных данных. Совместная обработка данных в нескольких спектральных диапазонах позволяет эффективно решать задачи радио- и гидролокационных комплексов, систем глобального мониторинга Земли, систем технического зрения и др. Для ряда приложений особый интерес представляют задачи обнаружения и оценивания параметров протяженных аномалий на многозональных изображениях. Вместе с тем, для выбора наиболее подходящего правила обнаружения актуальным является исследование и сопоставление эффективности соответствующих алгоритмов, и потенциальные возможности реальных систем обработки многозональных изображений.

В ряде работ предлагаются локальные алгоритмы обнаружения на многомерных изображениях, основанные на оценках градиента поля изображения. Рассматривается применение адаптивных алгоритмов компенсации коррелированных помех и алгоритмов обнаружения с «выбеливанием» многомерного изображения. Несмотря на большое число публикаций по проблемам синтеза алгоритмов обнаружения сигналов, в настоящее время мало исследована их эффективность. Существующие решения не дают удовлетворительного результата в задачах анализа эффективности и сравнения соответствующих алгоритмов обнаружения, особенно алгоритмов обнаружения протяженных аномалий на многомерных изображениях. Это обусловливает необходимость дополнительного исследования эффективности алгоритмов обнаружения аномалий.

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

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

Цель и задачи исследования

Целью работы является повышение эффективности систем мониторинга земной поверхности на основе разработки, моделирования и анализа вероятностных характеристик алгоритмов обнаружения протяженных неоднородностей на многозональных изображениях.

Для достижения этой цели необходимо решение следующих задач:

- анализ и разработка моделей формирования многозональных изображений, позволяющих при относительно небольших вычислительных затратах формировать изображения и их последовательности, имитирующие реальные многозональные изображения;

- синтез и анализ квазиоптимальных алгоритмов обнаружения протяженных аномалий на фоне мешающих многозональных изображений в условиях априорной неопределенности относительно уровней яркости сигнала;

- проведение сравнительного анализа эффективности известных и синтезированных квазиоптимальных алгоритмов обнаружения протяженных аномалий на фоне мешающих многозональных изображений;

- разработка пакета прикладных программ для анализа эффективности синтезированных обнаружителей на основе статистического моделирования, а также при использовании кадров реальных спутниковых изображений.

Методы исследований. При решении поставленных задач были использованы методы теории вероятностей, математической статистики, математического моделирования, численной аппроксимации функций, средства современной информатики и вычислительной техники, включая язык программирования Delphi 7 и программные средства работы с графическими изображениями Paint.NET и MultiSpec.

Научная новизна результатов

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

2. Синтезированы квазиоптимальные алгоритмы обнаружения протяженных аномалий на фоне многозональных изображений при наличии случайных помех в условиях неизвестных параметров полезного сигнала.

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

Практическая ценность диссертационной работы заключается в следующем:

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

2. Разработанный пакет программ допускает его непосредственное использование при проектировании современных и перспективных систем моделирования и обработки изображений. Предложенные модели многомерных СП могут быть использованы разработчиками для представления реальных данных при решении различных прикладных задач.

Апробация работы.

Основные результаты диссертационной работы докладывались и обсуждались на следующих научно-технических конференциях: 63 Научная сессия, посвященная Дню радио (Москва, 2008); Седьмая Международная конференция «Математическое моделирование физических, экономических, технических, социальных систем и процессов» (Ульяновск, 2009); Шестая Всероссийская научно-практическая конференция «Современные проблемы создания и эксплуатации радиотехнических систем» (Ульяновск, 2009); 13 Международная конференция по распознаванию образов и анализу изображений (Санкт-Петербург, РША, 2010); Научно-техническая конференция «Интегрированные автоматизированные системы управления» (Ульяновск, ФНПЦ ОАО «НПО «МАРС», 2011); ежегодные научно-технические

конференции профессорско-преподавательского состава Ульяновского государственного технического университета (2009-2012).

Реализация результатов работы.

Результаты работы использованы ФНПЦ ОАО «НПО «Марс». Реализация результатов работы подтверждена соответствующими актами.

Публикации. Основные научные результаты диссертационной работы отражены в 14 публикациях, среди которых 5 статей, из них 2 статьи в журналах из перечня ВАК, 5 работ в трудах и материалах международных и Всероссийских конференций. По результатам работы подана одна заявка на изобретение.

Структура диссертационной работы: Диссертационная работа написана на русском языке и состоит из введения, четырех глав, заключения, списка использованной литературы и приложения. Работа изложена на 123 страницах, содержит 67 рисунков и список литературы из 151 наименования.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обоснована актуальность работы, сформулированы цели и основные задачи, показана научная новизна и практическая ценность работы.

В первой главе рассмотрены математические модели, которые наиболее часто применяются для описания многозональных изображений, а также предложена методика моделирования многомерных изображений. Дается краткий обзор известных алгоритмов предварительной обработки спутникового материала и существующих методов обнаружения аномалий на изображениях. Формулируются основные нерешенные задачи, определившие тему диссертационной работы.

Во второй главе предложены алгоритмы анализа эффективности предварительной фильтрации многозональных изображений, направленные на

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

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

В четвертой главе представлены результаты экспериментальных исследований, иллюстрирующие эффективность синтезированных методов и алгоритмов при работе с реальными многозональными спутниковыми снимками. Изложены основные принципы организации разработанного программного обеспечения.

В заключении сформулированы основные научные и практические результаты диссертационной работы.

ГЛАВА 1. МЕТОДЫ МОДЕЛИРОВАНИЯ И АНАЛИЗА ЭФФЕКТИВНОСТИ АЛГОРИТМОВ ОБНАРУЖЕНИЯ АНОМАЛИЙ

1.1 Постановка задачи

Задачи обнаружения и оценивания параметров аномалий на многомерных изображениях представляют интерес для целого ряда приложений. Среди них можно выделить радио- и гидролокационные комплексы с пространственными антенными решетками, аэрокосмические системы глобального мониторинга Земли, системы технического зрения и др. [59Д01Д05,]. В последнее время задачи обработки многомерных данных приобретают особую актуальность в связи с широким распространением методов мультиспектральной и гиперспектральной регистрации участков земной поверхности. Это обусловливает необходимость разработки новых методов качественного и количественного анализа аэрокосмических наблюдений как единой многомерной совокупности. Для ряда приложений особый интерес представляют задачи обнаружения и оценивания параметров аномалий на многозональных изображениях.

Получение достаточного объема информации в натурных экспериментах довольно сложно и требует значительных вычислительных затрат. Одним из подходов к обеспечению исходными данными задач проектирования аэрокосмических систем является моделирование стохастических полей яркости [109,110,113,11]. Весьма желательно, чтобы математические модели позволяли решать задачи анализа и синтеза. Т.е., должна быть возможность по известным параметрам модели определить, какие вероятностные свойства имеет описываемое этой моделью многомерное изображение, и наоборот, можно было бы так выбрать параметры модели, чтобы она описывала многомерное изображение с заданными свойствами.

Особенно остро проблемы стоят в задачах, требующих не только правильно извлекать информацию из регулярно поступающего материала, но и

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

В настоящей главе рассмотрен ряд важных задач дистанционного зондирования Земли (п. 1.2), а также проанализированы основные известные методы имитации изображений, заданных в дискретном пространстве-времени (п. 1.3). В п. 1.4 представлен новый способ моделирования неоднородных изображений. В п. 1.5 приведен анализ характеристик дважды стохастической модели СП. В п. 1.6. предложена методика формирования изображений. В п. 1.7. дан аналитический обзор существующих алгоритмов оценки многозональных изображений и обнаружения аномалий на них, а также определен вид исследуемых кадров многозонального изображения.

1.2. Обзор задач дистанционного зондирования Земли

Дистанционное зондирование природной среды охватывает область человеческой деятельности, задачей которой является исследование природы неинвазивными методами как на расстоянии от исследуемого объекта, так и в непосредственной близости от него (вплоть до контакта) [59]. Примером могут служить аэрогеофизические �