Войти
Логин:
Пароль:
Забыли пароль?
научная деятельность
структура институтаобразовательные проектыпериодические изданиясотрудники институтапресс-центрконтакты
русский | english
Научная деятельность >> Прикладные разработки >> ПРИМЕНЕНИЕ МЕТОДОВ ТЕОРИИ ГИББСОВСКИХ ПО...

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

     В работе [1] был разработан алгоритм, позволяющий автоматически решать задачу детектирования. В основе этого алгоритма лежит стохастическая равновесная динамика рождения-гибели в пространстве конфигураций в непрерывном пространстве, стационарная мера которой имеет гиббсовский вид, и в которой заложены все необходимые требования на искомую конфигурацию. Этот алгоритм представляет собой эволюцию макро-объектов, находящихся в непрерывном пространстве. Существенным отличием этого алгоритма от существующих в настоящее время стохастических алгоритмов для непрерывных моделей, таких как, например, RJMCMC (Random Jump Monte Carlo Markov Chain), является возможность за один итерационный шаг произвести кардинальное изменение текущей конфигурации, что приводит к значительному ускорению сходимости. Обзор существующих алгоритмов и сравнение с новым подходом можно найти в наших работах [2, 5].

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

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

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

     Дальнейшее направление исследований связано с системами, содержащими макро-объекты с одним или более внутренними параметрами (степенями свободы), см. [3, 4]. Такие системы называются моделями маркированных точечных полей, а внутренние параметры объекта называются марками. Например, в случае динамики эллипсов, марками являются величины полуосей и ориентация эллипсов. При этом новая динамика помимо процесса рождения-гибели содержит диффузионную часть, отвечающую изменению этих марок. Здесь в зависимости от задачи возможно сочетание детерминированных и стохастических диффузионных динамик.

Литература

1. Descombes X., Minlos R.A., Zhizhina E., Object extraction using stochastic birth-and-death dynamics in continuum, Journal of
Mathematical Imaging and Vision, Vol.33 (3), p. 347, 2009.
2. Descombes X., Zhizhina E., The Gibbs fields approach and related dynamics in image processing, Condensed Matter Physics, 2008, Vol. 11, No.2, p. 1-20.
3. M. Kulikova, I. Jermyn, X. Descombes, J. Zerubia, E. Zhizhina, Extraction of arbitrarily shaped objects using stochastic multiple birth-and-death dynamics and active contours, IS&T/SPIE Electronic Imaging 2010 meeting, 17-21 January 2010, San Jose, California, USA.
4. M. Kulikova, I. Jermyn, X. Descombes, J. Zerubia, E. Zhizhina, A marked point process model with strong prior shape information for the extraction of multiple, arbitrarily-shaped objects, SITIS, The 5-th International Conference on Signal-Image Technology & Internet-based Systems, November 29- December 4, 2009, Marrakech.
5. E. Zhizhina, X. Descombes, Gibbs point field models for extraction problems in image analysis, Proceedings of Dobrushin International Conference, Moscow, Russia, July 15-20 2009, ISBN 978-5-901158-10-4.

НОВОСТИ И ОБЪЯВЛЕНИЯ
Открытый семинар Лаборатории зрительных систем №11: 20.01.2017 (пятница), 17:00, ауд. 615 ИППИ РАН. ...
В новом номере газеты "Троицкий вариант - Наука" в рамках проекта "Математические прогулки" опублико...
Защиты кандидатских диссертаций в совете Д.002.077.05: 23.01.2017 (понедельник), ауд.615 ИППИ...
Внимание! 24.01.2017, во вторник, на общем собрании коллектива ИППИ РАН состоятся выборы директора ...
Семинар по теории кодирования: 17.01.2017 (вторник),19:00, ауд.307 ИППИ. Кабатянский-Лебедев Г.-В.ru...
Семинар лаборатории № 8: 19 января в 14:30 в ИПЭЭ РАН. В.М. Ольшанский: Про закон...
Врио директора ИППИ РАН Андрей Соболевский посетил БГПУ им. М. Акмуллы в Уфе. Целью визита стало обс...
Cеминар "Дискретная и вычислительная геометрия": 17.01.2017 (вторник), 13:45, ауд. 307 ИППИ РАН. Вла...
18.01.17 (среда), ФИАН - Научная сессия ОФН РАН, посвященная 100-летию со дня рождения...
Семинар лаборатории № 8: 12 января в 14:30 в ИПЭЭ РАН. Е.М. Максимова. Коннектом глицинергических ам...
Семинар <<Глобус>> Независимого Московского Университета: 13 января в 15.40 в НМУ, конференц-зал. И....
Семинар <<Глобус>> Независимого Московского Университета: 12 января в 15.40 в НМУ, конференц-зал. Ле...
Итоги-2016: Дюжина самых ярких событий ИППИ РАН уходящего...
Сотрудники Центра распределенных вычислений ИППИ РАН об облачной платформе "Эверест" в программе "Ч...
Поздравление врио директора ИППИ РАН Андрея Соболевского с 55-летием Института и наступающими праздн...
В совместной работе сотрудников ИППИ РАН, МФТИ, НИУ ВШЭ, компании Яндекс и Université catholi...
Семинар "Структурные модели и глубинное обучение", внеочередное заседание: 26.12.2016 (понедельник),...
Семинар по теории кодирования: 27.12.2016 (вторник),19:00, ауд.307 ИППИ. Александр Варди (UCSD Jaco...
Семинар Добрушинской математической лаборатории: 27.12.2016 (вторник), 16:00, ауд. 307 ИППИ РАН. Дан...
Сотрудники Сектора анализа данных в нейронауках №10.3 Анвар Курмуков и Юлия Додонова - лауреаты отк...
Все новости   
 

 

  © Федеральное государственное бюджетное учреждение науки
Институт проблем передачи информации им. А.А. Харкевича Российской академии наук, 2017
Об институте  |  Контакты  |  Старая версия сайта