скачать На правах рукописи БЕСПАЛОВ Дмитрий Анатольевич РАЗРАБОТКА И ИССЛЕДОВАНИЕ МЕТОДА ОБНАРУЖЕНИЯ ОБЪЕКТОВ ИСКУССТВЕННОГО ПРОИСХОЖДЕНИЯ СРЕДСТВАМИ ВЕЙВЛЕТ-АНАЛИЗА Специальность: 05.13.17 – Теоретические основы информатики АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук Таганрог – 2007 Работа выполнена в Технологическом институте Южного федерального университета в г. Таганроге. Научный руководитель: доктор технических наук, профессор Гузик Вячеслав Филиппович. Официальные оппоненты: доктор технических наук, профессор Финаев Валерий Иванович. кандидат технических наук, старший научный сотрудник Заковоротнов Евгений Анатольевич. Ведущая организация: Федеральный научно-производственный центр ОАО «Научно-производственное объединение «МАРС» г. Ульяновск. Защита диссертации состоится 28 июня 2007 г. в 1420 на заседании диссертационного совета Д 212.208.21 при Южном федеральном университете по адресу г. Таганрог. 347928, Таганрог, пер. Некрасовский, 44, ауд. Д-406. С диссертацией можно ознакомиться в Зональной научной библиотеке ЮФУ по адресу Пушкинская, 148, г. Ростов-на-Дону. Автореферат разослан « 25 » мая 2007 г. Ученый секретарь диссертационного совета Д 212.208.21, доктор технических наук, профессор Бабенко Л.К. ^ Актуальность темы. Задача совершенствования методов цифровой обработки сигналов, в общем, и методов обнаружения объектов искусственного происхождения, в частности, является одной из важнейших во многих отраслях науки и техники: гидролокации, сейсморазведке, обработки изображений, видео и звука, медицине и т.п. Актуальность работы в таком случае, обуславливается расширением сферы применения методов обнаружения, внедрением новой теоретической и алгоритмической базы, использованием новых методов оптимизации на уровне аппаратных и программных средств поддержки вычислительного процесса. До последнего времени, одним из основных способов анализа первичной или вторичной информации, отражающей состояние интересующей исследователя среды, являлось визуальное наблюдение. Однако с ростом производительности вычислительных систем сбора и обработки данных возрастает роль автоматизированных методов выделения и классификации разного рода объектов. При применении таких систем делается попытка заменить человека системой с учетом требований к скорости и точности решения поставленной задачи, а также, если такая замена диктуется невозможностью выполнения заданных требований человеком. Важный научный интерес и практическое значение в таком случае представляет собой проблема построения вычислителей, способных выделять и классифицировать изображения объектов. Одной из малоисследованных задач является задача построения системы, способной производить автоматическое выделение двумерных и трехмерных объектов из окружающей среды с последующей их классификацией при наличии неопределенности в исходных данных. Нечеткость данных, ровно как неполнота их, требуют применения нетривиальных методов анализа. При этом необходимо создание систем, способных синтезировать определенный интуитивный и логический подход к анализу ситуации и принимать решения относительно поставленной задачи. Таким образом, задача построения систем, оперирующих методами выделения и первичной классификации искусственных объектов с учетом наличия неопределенности исходных данных, является актуальной и востребованной. Данная работа является развитием результатов исследований, проводимых на кафедре вычислительной техники ТТИ ЮФУ под руководством д.т.н., профессора Гузика В.Ф. и д.т.н., профессора Золотовского В.Е. в рамках научного направления: «Многопроцессорные потоковые вычислительные системы, компьютерные сети и сетевые технологии региональной информатизации». ^ В диссертации решаются теоретические проблемы решения задач обнаружения и первичной классификации объектов искусственного происхождения при неопределенности исходных данных средствами вейвлет-анализа. В ней проводится разработка и исследование комплекса алгоритмов обработки данных различной физической природы в задаче обнаружения, а также рассматриваются варианты их параллельной реализации на базе многопроцессорных кластеров процессоров цифровой обработки сигналов и системах на кристаллах ПЛИС. Для достижения этой цели необходимо решение следующих задач:
^ Для теоретических исследований использовались методы теории вейвлет-анализа, цифровой обработки сигналов, численные методы, методы математического моделирования, методы параллельных вычислений, распознавания образов и решения задач классификации. Кроме того, проводились экспериментальные исследования на реальных и модельных объектах и сигналах и апробация разработанных алгоритмов на базе ЭВМ. ^ заключается в создании нового метода обнаружения объектов искусственного происхождения средствами вейвлет-анализа, синтезе комплекса алгоритмов обработки данных и их реализации на базе аппаратных и программных средств поддержки вычислений. К числу наиболее важных результатов диссертации относятся:
^ полученных в работе результатов подтверждается полнотой и корректностью исходных посылок, проведенным теоретическим обоснованием и практической реализацией. Исследования и эксперименты были осуществлены при решении широкого класса задач обнаружения объектов (обнаружение эталонных объектов, выделение характерных особенностей сигналов в отсутствии эталонов, анализ сильно зашумленных сигналов) с использованием как универсальных ЭВМ, так и действующих прототипов и моделей многопроцессорных комплексов ПЦОС и потоковых вычислителей на базе ПЛИС. ^ . В диссертационной работе решена важная научно-техническая задача, заключающаяся в развитии методов обработки сигналов различной физической природы в задаче обнаружения и первичной классификации объектов искусственного происхождения. ^ заключается в создании нового метода обнаружения объектов в условиях неопределенности исходных данных, позволяющего проводить комплексный анализ сигналов различной физической природы, его очистку от шумов адаптивными порогами, выделение компактного кода, разделение сигнала на области вероятного нахождения объектов, как в масштабе, так и в пространстве, синтез векторов характеристических признаков и первичную классификацию локализованных объектов. Использование этого метода снижает трудоемкость этапов подготовки данных и сокращает время выполнения анализа сигналов, переводя расчеты в реальный масштаб времени за счет использования быстрых алгоритмов, реализованных на базе современных аппаратных средств методами параллельных вычислений. ^ заключается в синтезе методов, алгоритмов и комплекса программ обработки сигналов различной физической природы; а также в разработке ряда конфигураций аппаратных средств для их реализации методами параллельных вычислений. Их использование позволяет сократить время разработки новых комплексов цифровой обработки сигналов различной физической природы, а также расширить возможности существующих, посредством внедрения предложенных в диссертационной работе результатов. ^ . Основные теоретические и практические результаты диссертационной работы использованы при проведении научно-исследовательских работ кафедрой Вычислительной техники Технологического института Южного федерального университета (ВТ ТТИ ЮФУ) совместно с ФНПЦ ФГУП «НПО «Марс», ОАО НИИ «Системотехники», а так же в ряде госбюджетных и хоздоговорных работ. Материалы диссертации использованы в учебном процессе на кафедре ВТ ТТИ ЮФУ при проведении практических занятий в цикле лабораторных работ по курсу: «Проблемно-ориентированные высокопроизводительные системы». ^ . Основные результаты диссертации докладывались и обсуждались на всероссийских и международных научно-технических конференциях и выставках. В частности на научно-техническом семинаре с международным участием «Микропроцессорные системы мониторинга, диагностики и управления сложными техническими объектами, организационными техническими системами и комплексами». - Таганрог, 2004 г., VIII Всероссийской научной конференции студентов и аспирантов. «Техническая кибернетика, радиоэлектроника и системы управления». - Таганрог, 2004, L Научно-технической конференции профессорско-преподавательского состава, аспирантов и сотрудников ТРТУ. – Таганрог, 2004, III Всероссийской научной конференции молодых ученых, аспирантов и студентов «Информационные технологии, системный анализ и управление». - Таганрог, 2005, Всероссийской научно-технической конференции с международным участием «Компьютерные технологии в науке, инженерии и управлении». - Таганрог, 2005, Всероссийской научной конференции студентов и аспирантов «Техническая кибернетика, радиоэлектроника и системы управления». – Таганрог, 2006, Всероссийской научно-технической конференции с международным участием «Компьютерные технологии в науке, инженерии и управлении» - Таганрог, 2006. ^ , выносимые на защиту:
^ Все научные и практические результаты, а именно: анализ методов обнаружения, разработка требований к синтезируемому методу обнаружения, синтез комплекса алгоритмов первичной обработки данных, шумоподавления адаптивными порогами, многомасштабной сегментации, выделения векторов характеристических признаков и первичной классификации локализованных объектов, разработка комплексного метода обнаружения, создание ряда конфигурационных проектов для реализации на базе аппаратных средств поддержки вычислительного процесса методами параллельных вычислений, получены автором лично. Публикации. По теме диссертационной работы опубликовано 14 печатных работ, из них 9 статьей и научных сообщений, 4 тезиса докладов на российских и международных научно-технических конференциях, а также 1 методическое указание по выполнению лабораторных работ. ^ . Материал основной части диссертационной работы изложен на 222 страницах машинописного текста. Диссертация состоит из введения, четырех разделов, списка литературы из 135 наименований, содержит 5 таблиц, 48 рисунков и 2 приложения на 26 листах. ^ Во введении приводиться обоснование актуальности темы диссертационной работы, описывается цель работы, определяется круг задач, формулируется объект и предмет исследования, указываются методы исследования, сообщается научная новизна и практическая значимость, а так же подтверждается достоверность полученных результатов. ^ работы посвящен аналитическому обзору существующих подходов к решению задачи обнаружения и носит постановочный характер. В нем также проводится исследование большинства стандартных методов обнаружения с целью выявления их основных возможностей, особенностей реализации, а также основных достоинств и недостатков. В разделе проводится постановка задачи обнаружения объектов искусственного происхождения, а также вводятся основные понятия, используемые при решении задачи обнаружения. В частности, вводятся понятия объекта искусственного происхождения и фона на множестве двумерных сигналов и поверхностей, понятие характеристических точек объектов обнаружения и неопределенности исходных данных. В качестве основных характеристик исследуемых методов обнаружения объектов искусственного происхождения принимаются: инвариантность относительно пространственного и масштабного сдвига, локализация геометрических размеров и координат объектов, устойчивость к шуму, локализация и выделение свойств сигналов, наличие быстрых вычислительных алгоритмов, адаптивность к особенностям анализируемых данных. Специфика применения методов обнаружения в составах сложных вычислительных комплексов обработки сигналов различной физической природы, а также наличие высоких требований к их адаптивности, помехоустойчивости, скорости работы и возможности многопрофильного анализа данных требует разработки нового метода, удовлетворяющего перечисленным выше требованиям. В соответствии с этим, был проведен анализ современных методов обнаружения объектов искусственного происхождения, выполнена их классификация, выявлены их основные достоинства и недостатки. Результаты сравнительного анализа сведены в единую таблицу, представленную на рисунке 1. ![]() Рисунок 1 – Таблица сравнительного анализа представленного решения с существующими методами обнаружения. Предлагаемый метод базируется на средствах вейвлет-анализа и ориентирован на комплексную обработку сигналов, включающую в себя: определение наилучшего базиса разложения, вейвлет-фильтрацию адаптивными порогами, разделение на области вероятного нахождения искомых объектов, выделение характерных пространственных и масштабных свойств и векторов признаков, а также их первичную классификацию. Наряду с решением прямой задачи обнаружения, формируется компактный код данных, позволяющий в несколько раз сократить объемы передаваемой, хранимой и обрабатываемой информации, проводится дополнительный анализ данных с выделением множеств значений локальной гладкости, локализацией нарушений регулярности и закономерностей текстур объектов и фона, а также достигается параллельная организация вычислительного процесса на базе современных аппаратных средств. Таким образом, по результатам анализа можно заключить, что предлагаемый метод обладает практически всеми ключевыми характерными признаками, реализующими основные этапы обработки данных в задаче обнаружения. При этом, в полной мере данными свойствами не обладает ни один из рассмотренных методов. ^ является ядром теоретического обоснования предложенного метода обнаружения. Он содержит основы самой теории вейвлет-анализа, теоретические предпосылки создания отдельных элементов метода, а также его алгоритмический и методологический базис в целом. Проводится исследование возможностей методов и алгоритмов теории вейвлет-анализа. На основании этого делается вывод о целесообразности применения методов цифровой вейвлет-обработки сигналов для решения задачи обнаружения объектов искусственного происхождения. Большое внимание в разделе уделено вопросам объединения преимуществ отдельных алгоритмов вейвлет-анализа с целью поиска новых путей их использования. В частности делается вывод, что основным источником информации о структуре и свойствах сигналов является ограниченное множество максимумов ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Совокупность последовательно повторяемых преобразований компонент сигнала создают многоуровневую древовидную структуру его ![]() ![]() ![]() Вследствие особых свойств базисных функций ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Рисунок 2 – Алгоритм Малла для 3-х уровней разложения. Известно, что характерные черты сигналов сосредоточены в областях с резкими изменениями его амплитуды и частоты, то есть в местах разрывов, неоднородностей и перепадов. Показано, что именно на такие участки сигналов лучше всего «реагирует» вейвлет-преобразование. Полученное множество вейвлет-коэффициентов большой амплитуды чрезвычайно важно, так как определяет основной состав и свойства сигнала, другая часть коэффициентов (малой амплитуды) может быть опущена, так как является результатом преобразования его высокочастотных (мелкомасштабных) компонент, наиболее часто являющихся шумовыми. ![]() ![]() Рисунок 3 – Линии максимумов вейвлет-коэффициентов в точках локальных неоднородностей сигналов. В классической теории вейвлет-анализа удаление коэффициентов малой амплитуды играет роль фильтра, носит название трешолдинга и выполняется пороговым отсечением деталей на каждом уровне декомпозиции. Основной проблемой является выделение такого значения отсекающегося порога, чтобы потери полезной информации при обработке были минимальны. Таким образом, для решения задачи адаптивного шумоподавления необходимо провести адаптивную генерацию локальных порогов, позволяющих уменьшить влияние аддитивного шума на чистую форму сигнала, и сохранить значимые коэффициенты большой амплитуды, характеризующие локальные особенности анализируемых данных Вследствие этого была предложено улучшить классическую схему трешолдинга путем введения процедуры адаптации к особенностям анализируемых сигналов. Нахождение величины данного порога ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() Результатом работы предложенной процедуры является ограниченное множество вейвлет-максимумов, практически полностью определяющих структуру и свойства сигнала. Использование именно этого ограниченного множества данных позволяет выделить так называемый компактный код данных и сократить объемы передаваемой, хранимой и обрабатываемой информации на 70-90% для каждого уровня декомпозиции, то есть фактически обеспечить коэффициент сжатия до 10-ти раз без применения дополнительных мер сверхкомпрессии, уменьшая тем самым суммарную трудоемкость предложенного метода и загруженность вычислительных ресурсов. Кроме того, выделенные коэффициенты позволяют локализовать места неоднородностей и перепадов анализируемых сигналов, то есть провести его пространственное разделение на области с характерными особенностями. Это дает возможность проводить обнаружение объектов искусственного происхождения не во всем поступающем сигнале, а только в отдельных его сегментах, свойства которых отличаются от свойств фона. С другой стороны, вейвлет-преобразование обладает способностью к масштабному и частотному разделению сигналов. Частотно-временное разрешение преобразования ![]() ![]() Каждое разложение анализируемого сигнала ![]() ![]() ![]() ![]() ![]() ![]() ![]() Переход вниз или вверх по дереву декомпозиции сопровождается потерей или приобретением некоторого количества информации. Добавляемые или исключаемые детали составляют элементы некоторого дополнительного множества функциональных подпространств ![]() ![]() ![]() ![]() ![]() ![]() ![]() Добавление этой компоненты в состав сигнала ![]() ![]() ![]() ![]() ![]() В результате исследования был сделан вывод, что каждый объект оставляет свой «след» на множестве уровней декомпозиции. Таким образом, для каждого объекта можно выделить верхнюю границу ![]() На основании этого была предложена процедура многомасштабной сегментации, позволяющей на основании анализа линий максимумов, выделенных на этапе адаптивной фильтрации, производить разделение сигналов на области «наибольшего интереса» с характерными свойствами, как во временном, так и в масштабном (частотном) пространствах с локализацией перепадов, выделением многомасштабной геометрии и контуров объектов ![]() ![]() Полученные данные достаточно хорошо описывают особенности сигнала вследствие инвариантности вейвлет-преобразования относительно масштабного и пространственного сдвига. Они могут быть использованы для идентификации типопринадлежности локализованного объекта в процедуре первичной классификации. На основании проведенного исследования было выделено множество векторов характеристических признаков локализованных объектов, определяемое на основании результатов обработки данных предложенными процедурами выделения вейвлет-максимумов, фильтрации и сегментации, а также проведена их классификация и определены границы их изменений для локализуемых объектов. Раздел завершается выделением этапов предложенного метода первичной обработки анализируемых сигналов в задаче обнаружения объектов искусственного происхождения средствами вейвлет-анализа:
^ диссертационной работы посвящен алгоритмической составляющей разработанного метода обнаружения объектов искусственного происхождения, анализу предложенных алгоритмов с точки зрения их вычислительной трудоемкости и способам параллельной реализации вычислительного процесса на программном и аппаратном уровне. Введением в данный раздел является анализ требований к алгоритмам обнаружения, определяющих общие черты вычислительного процесса и процесс обнаружения в целом. Далее, в разделе последовательно рассмотрены следующие этапы предложенного метода обнаружения с точки зрения их алгоритмической реализации: быстрые алгоритмы адаптивной декомпозиции/реконструкции сигнала в базисе вейвлет-функций; алгоритм адаптивной вейвлет-фильтрации с оптимальными порогами; алгоритм многомасштабной сегментации анализируемого сигнала; выделения векторов масштабных локальных и глобальных признаков сегментов; алгоритм первичной классификации обнаруженных объектов, присутствующих в базе эталонов. Рассмотрены структурные схемы данных алгоритмов и особенности их реализации. В результате работы над данным разделом была разработана структурная схема и описание комплексного алгоритма обнаружения объектов искусственного происхождения соответствующие предложенному методу. Также выделены его особенности и варианты реализации. Далее проводится анализ трудоемкости разработанных алгоритмов согласно классической теории и выполняется оценка каждого этапа в отдельности с выделением аналитических зависимостей для наихудшего, наилучшего и среднего случаев с графической демонстрацией полученных результатов. На основании проведенного исследования был сделан вывод, что трудоемкость предложенного комплексного алгоритма является линейной. Квадратичная зависимость наблюдается только на этапе генерации порога фильтрации и анализа геометрии объектов, однако доля данных этапов вследствие использованного ограниченного множества вейвлет-коэффициентов в общем алгоритме невелика. На основании исследования особенностей структуры и операционного базиса алгоритмов вейвлет-анализа был сделан вывод о возможности и целесообразности реализации предложенной процедуры обнаружения методами параллельных вычислений. Исследование показало наличие естественного параллелизма всех предложенных этапов обработки как на уровне операционного базиса, алгоритмической организации потоков данных, так и на уровне использования параллельных структур аппаратных средств поддержки вычислительного процесса. Проведенное экспериментальное исследование показало эффективность предложенного подхода. ^ поддержки вычислительного процесса с точки зрения целесообразности и эффективности реализации на их базе разработанных алгоритмов. После чего был разработан и реализован ряд конфигурационных проектов и экспериментов, подтверждающих высказанные предположения. В качестве альтернативных решений были выбраны направления процессоров цифровой обработки сигналов типа ADSP 21xxx и систем на базе программируемых интегральных схем типа Cyclone. Для первого варианта были реализованы и оценены варианты параллельного исполнения для 2-х и 4-х процессов, распределенных на вычислительных мощностях четырехпроцессорного кластера ADP160QPCI (рисунок 4). ![]() Рисунок 4 – Структурная схема программного обеспечения многопроцессорного кластера. Для второго варианта был также реализован конфигурационный проект 2-х и 4-х потокового вычислителя на базе Cyclone (рисунок 5). Ограничениями второго варианта являлись: ограниченная разрядность данных, целочисленная форма представления сигналов и самого вейвлет-преобразования, ограничение ресурсов памяти, возможности расширения и реконфигурации. ![]() Рисунок 5 – Структурная схема конфигурации ПЛИС. В результате проведенного исследования, эффективность параллельного решения, вычисленная как время выполнения вычислительного процесса относительно однопроцессорного варианта для двух- и четырехпроцессорных комплексов ПЦОС, составила 50,2% и 25,4% соответственно. Для вариантов целочисленного вейвлет-преобразования с соответствующими ограничениями для систем на кристаллах ПЛИС - 3,76% и 2,35% соответственно. В четвертом разделе приводятся результаты экспериментального исследования и практического применения разработанного метода обнаружения. Введением в раздел является анализ требований к программной реализации предложенного метода, а также постановка задачи обнаружения с позиции его экспериментального исследования. ^ , включающая использование разработанных алгоритмов обнаружения объектов искусственного происхождения как модели обнаружителя в составе комплекса анализа гидроакустической обстановки (рисунок 6). Для этого был определен состав блоков системы и разработана методика исследования, состоящая из трех этапов. Для каждого этапа было выделено множество оцениваемых параметров, оказывающих наибольшее влияние на результат работы алгоритмов, определены зависимости и связи между ними (рисунок 6). В результате реализации предложенной схемы экспериментального исследования был проведен анализ предложенных алгоритмов обработки данных для выделенного множества тестовых и реальных сигналов с заданными характеристиками в рамках поставленной задачи обнаружения объектов искусственного происхождения. ![]() Рисунок 6 – Структурная схема модели системы. В результате проведения экспериментов, было проведено выделение наиболее значимых параметров алгоритмов, определены их оптимальные значения и ограничены диапазоны их изменения. ![]() Рисунок 7 – Структура связей между оцениваемыми параметрами внутри этапов экспериментального исследования. ^ , что: наилучшим базисом для анализируемого класса сигналов является базисе вейвлетов Добеши порядков 1-4, погрешность восстановления при мягкой и жесткой схемах составили не более 5%-7%, количество ненулевых коэффициентов в среднем 10%-15%, коэффициент компрессии 7,2-9,3 раза и вероятность обнаружения не менее 95% для среднего уровня шумов. После этого, была решена задача поиска навигационных опасностей по данным гидроакустического сканирования акватории с затонувшим теплоходом «Адмирал Нахимов», в ходе которой проверена работоспособность предложенного метода и получены желаемые результаты. Были четко выделены места неоднородностей рельефа донной поверхности, проведена очистка сигналов, выделение контуров корпуса и основных элементов затонувшего судна, а также выделен ряд техногенных объектов, расположенных вдоль границ основного контура. Результаты работы были использованы при выполнении х/д работ кафедры ВТ, что подтверждается соответствующими актами внедрения. ^ сформулированы основные выводы и результаты диссертационной работы. В приложении приводятся примеры программной реализации синтезированных алгоритмов, функциональные схемы, данные конфигурационных проектов и описание блоков для аппаратных средств, а также акты об использовании результатов диссертационной работы. ^
^
В работах, опубликованных в соавторстве, лично автору принадлежат следующие результаты: в работе [3] – метод обработки гидроакустической информации на основе множеств коэффициентов вейвлет-преобразования, в работе [4] – комплексный алгоритм обработки сигналов различной физической природы, а также практическая реализация предложенного алгоритма на базе ПЦОС ADSP-21xxx, в работе [5] – метод масштабной и пространственной сегментации сигналов на основе анализа линий максимумов коэффициентов вейвлет-преобразования, в работе [6] – метод и основной алгоритм компрессии сигналов на основе ограниченной выборки коэффициентов вейвлет-преобразования, в работе [7] – модель распределения вычислительных потоков комплекса первичной и вторичной обработки гидроакустических данных. Аспирант Беспалов Д.А.
|