Выделение малоразмерных объектов алгоритмами сегментации на основе модели системы со случайной скачкообразной структурой

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

А.Н. Малое, Б.М.Миронов, В.А.Кузнецов Иркутское высшее военное авиационное инженерное училище

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

Введение

Одной из важных задач обработки изображений является их сегментация. Задача сегментации изображений практически неизбежно возникает при построении различных систем обработки двумерных и трехмерных изображений, систем обнаружения и систем распознавания образов [1]. В большинстве случаев под сегментацией понимают выделение на изображении однородных участков (ОУ). Понятие "однородность участка" определяется условиями задачи, например, однородными могут считаться области с одинаковыми параметрами распределения яркости, одинаковыми параметрами текстуры, типом поверхности и т.д. В данной работе под сегментацией будем понимать разделение локационного изображения подстилающей поверхности на ОУ, соответствующие определенным типам (классам) подстилающей поверхности (например- лес, поле и т.д.) со средним значением на изображении величины т(/), где /-номер класса подстилающей поверхности.

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

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

Цель данной работы - построение алгоритмов сегментации изображений со сглаживанием достаточной статистики и исследование свойств алгоритмов на основе модели системы со случайной скачкообразной структурой по выделению малоразмерных объектов.

1. Алгоритмы сегментации со сглаживанием достаточной статистики

Аналогично [2] процесс формирования яркост-ных отсчетов Ак строки изображения опишем моделью одномерной авторегрессии. С учетом обозначений, принятых в [2], процесс смены классов подстилающей поверхности описывается дискретной марковской последовательностью {{к,/ = 0,1,...}, состояния которой / являются номерами классов (/ = 1,2,..., М).

Для случая, когда возможно запоминание всей строки изображения А = {Я0,Я1,...,ЯЫ}, задача состоит в оценивании А0 А,...,А на основе наблюдений аА1 . Оценку номера класса будем определять по критерию максимума апостериорной вероятности. Сглаженная достаточная статистика удовлетворяет рекуррентному уравнению в обратном направлении отсчетов (к = N -1,N - 2,...0) [3]:

WI,

А = в) = wk (вк =в)£-

р

Х^к +1,N

=1 X рАА А+1 = ЗА = /)п

с начальным условием:

^ ^ (А) = ^ (А):

где Жк N Ак = в) = Р(А = в^А ) - сглаженная апостериорная вероятность принадлежности к-го ЭИ к в-му классу (в = 1,2,...,М); Wк(вк = в) =

Р(вк = вАА) - апостериорная вероятность принадлежности к-го ЭИ к в -му классу, полученная при обработке строки изображения в прямом направлении отсчетов [2]; р(Ак+1\Ак,вк+ъвк), п(вк+х\вк) -условные плотности вероятности и переходные вероятности соответственно. Алгоритм сегментации, определяемый выражением (1), назовем одностроч-

1=1

ным алгоритмом сегментации со сглаживанием достаточной статистики - сокращенно ОА2 в отличие от однострочного алгоритма, представленного в работе [2], который сокращенно назовем ОА1.

При обработке изображения по строке и столбцу в соответствии с выражением (1) путем усреднения соответствующих апостериорных вероятностей получим комбинированный алгоритм сегментации со сглаживанием достаточной статистики, назовем его КОА2. Комбинированный алгоритм сегментации, полученный в работе [2], назовем КОА1.

2. Выделение алгоритмами сегментации малоразмерных объектов

Под малоразмерными объектами будем понимать ОУ, сравнимые по площади с размерами элемента разрешения локатора.

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

уровнями яркостей т(1) и т(2). В качестве модели изображения выбиралось тестовое поле размером 150х150 ЭИ, представленное совокупностью ОУ фона и объектов, сформированных в соответствии с выражением (4) в [2]. Параметры моделей определялись следующим образом: т(1) = т0 = 150 ; т(2) = т0 +Ат ; рХ'2 = 0,1; ст(1,2) = 2 . Ошибка в

определении номера класса, аналогично [2], оценивалась значением РОШ, называемой ошибкой распознавания состояния:

Р =

1 ОШ

1

N L

l ,k

КЬ к=1 г=1

где Ь, N - размеры изображения по вертикали и горизонтали, величина, принимающая значение

ноль в случае, когда оценка номера класса в1к и

номер класса в^к 1,к-го ЭИ совпадают, и равная единице в противном случае.

Свойства подстилающей поверхности определяются элементами матрицы переходных вероятностей (МПВ):

п(0к) =

п.

Для описания результатов исследований введем коэффициенты К12 и К22, определяемые МПВ (3):

K12

K 22 =-

22

'21

Значение К12 увеличивается при росте числа объектов на изображении и уменьшении их размеров. В

процессе исследований линейные размеры объектов изменялись от 49х49 ЭИ до 3х3 ЭИ. Коэффициент К22 характеризует соотношение площади объектов и фона на изображении. С целью детального исследования свойств алгоритмов ошибка распознавания состояния рассчитывалась как по всему изображению (ALL IMAGE), так и по площади выделяемых объектов (object).

Сначала исследовалась зависимость Рош от К12 при фиксированных значениях К22. Графики зависимости представлены на рис. 1, 2 для К22=0,1 и К22=0,01, при которых площадь объектов составляет 10% и 1% от площади фона на изображении соответственно. Из графиков видно, что с увеличением К12, то есть с увеличением числа объектов и уменьшением их размеров, ошибка распознавания растет. При этом ошибка по всему изображению больше в том случае, когда больше площадь объектов на изображении. Ошибка по площади объектов больше при меньшем значении К22. Такое поведение кривых на графиках объясняется влиянием МПВ (3) на работу алгоритмов. Наиболее эффективны алгоритмы сегментации со сглаживанием достаточной статистики ОА2 и КОА2. Как видно из графиков, при малых размерах объектов (К12=0,5) эффективность алгоритмов сегментации становится низкой - РОШ достигает значений более 0,5. Поэтому размеры сегментируемых участков объектов не могут быть менее 5x5 ЭИ.

ALL IMAGE

-*— OA1JC22=0,1

*..... KOAl

-*— ОА2

Рис.1. Зависимость ошибки распознавания состояния по всему изображению от коэффициента К12

Далее проводилось исследование зависимости РОШ от К22. Графики зависимости представлены на рис. 3, 4 для К12=0,02 и К12=0,5. Из графика на рис. 3 видно, что при увеличении соотношения площади объектов и фона на изображении РОШ по всему изображению увеличивается, а по площади объектов уменьшается (рис. 4). При этом при достаточно больших значениях К12 (К12=0,5) ошибка имеет значения более 0,5. Таким образом, удовлетворительная работа алгоритмов сегментации возможна при ограничении размеров сегментируемых участков объектов.

п

л

л

21

22

Ж

Л

Л

21

22

12

X

ж

OBJECT

-х-— ОА1.К22=0,1   *...... KOÄ1

0,25 0,50 К12

Рис.2. Зависимость ошибки распознавания состояния по площади объектов от коэффициента К12

ALLIMAGE

Рис.3. Зависимость ошибки распознавания состояния по всему изображению от коэффициента К22

разности средних уровней яркостей классов объектов и фона. Из рис. 5,6,7,8 видно, что с увеличением Am ошибка уменьшается, причем при Am = 10 -15 уменьшается до нуля в зависимости от применяемого алгоритма. При увеличении соотношения площади объектов и фона на изображении при фиксированных значениях К12 (рис. 5) ошибка распознавания по всему изображению увеличивается, а по площади объектов - уменьшается. Как отмечалось выше, это объясняется влиянием МПВ (3) на работу алгоритмов. При уменьшении размеров объектов и увеличении их числа (с увеличением К12) при фиксированных значениях К22 ошибка распознавания состояния увеличивается, достигая больших значений по площади объектов (рис. 7, 8). Из сравнения графиков, приведенных на рис. 5, 6 и 7, 8, следует, что эффективность сегментации изображений в большей степени определяется размерами и числом объектов на изображении, чем соотношением площади объектов и фона.

ALL IMAGE

ОА1.К12=0,02,К22=0,1

KOÀ1

ОА2

КОА2

ОА1.К12=0,02,К22=0,01

KOÀ1

ОА2

КОА2

15

20

2 0,7

S,

OBJECT

-х— ОА1.К22=0,02

х............KOÀ1

-к— ОА2

- КОА2

— ОА1,К22=0,5 ..........KOÀ1

— ОА2

— КОА2

Рис.4. Зависимость ошибки распознавания состояния по площади объектов от коэффициента К22

Рис.5. Зависимость ошибки распознавания состояния по всему изображению от Ат при К12=Сош1

3 0, £

Ki s

1°' л

U

о s

Iй

о,

to,

07

Об

05

04.

03

OBJECT

V

л-

X \

Ц \ \

; 0,02-i--—------

o,oi{

-х-—

X......

-X--

ОА1,К12=0,02; К22=0,1

KOÂ1

ОА2

КОА2

ОА1.К12=0,02; К22=0,01

KOÂ1

ОА2

КОА2

10

15

20

Следующие графики отражают зависимость РОШ (рис. 5, 6, 7, 8) от параметра Ат, равного

Рис.6. Зависимость ошибки распознавания состояния по площади объектов от Ат при К12=Сош1

ALL IMAGE

i 0,6

40,1 ¡о,с I 0,1

10,1

-К— ОА1,К12=0,02; К22=0,1 к..... KOÀ1 -к— ОА2 -к— КОА2 -#— ОА1,К12=0,5; К22=0,1 ..... KOÀ1 ОА2 -ф— КОА2

Тч\

К\______

tvA 1

г \\ \ \\\

Г

г- ^

P^Tirii

10

15

20

Рис.7. Зависимость ошибки распознавания состояния по всему изображению от Ат при К22=Сот1

OBJECT

-к— ОА1.К12=0,02; К22=0,1

х...... KOÄ1

-к— ОА2

— КОА2

— ОА1.К12=0,5; К22=0,1 .....KOÄ1

— ОА2

— КОА2

20

Рис.8. Зависимость ошибки распознавания состояния по площади объектов от Ат при К22=Сош1

Для иллюстрации работоспособности рассмотренных алгоритмов ниже приведен пример обработки реального изображения городской застройки, полученного когерентным локатором и представленного на рис. 9, с помощью алгоритма сегментации КОА2 (рис. 10) при следующих условиях: М=3; т(1) = 54,

а(1) = 14,

да(2) = 99, m(3) = 197;

= 0,1;

сг(2) = 30 , сг(3) = 51. Из рисунков видно, что на обработанном изображении малоразмерные объекты, соответствующие типу подстилающей поверхности с

т(3) = 197 , выделены с сохранением их границ.

Заключение

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

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

Рис.9. Реальное изображение когерентного локатора

Рис.10. Обработанное алгоритмом сегментации КОА2 изображение

Благодарность Работа выполнена при поддержке Российского фонда фундаментальных исследований, грант РФФИ № 06-08-00596-а.

Литература

1. Верденская Н.В. Сегментация изображений- статистические модели и методы // Успехи современной радиоэлектроники, 2002.№ 12. С.33- 47.

2. Скрыпник О.Н., Формирование классификационной карты подстилающей поверхности по изображениям от когерентного локатора // Скрыпник О.Н., Лежанкин Б.В., Малов А.Н., Миронов Б.М., Галиев С.Ф.- Компьютерная оптика, 2006. В. 29. С. 151-159.

3. Клекис Э.А. Оптимальное обнаружение скачкообразных изменений структуры дискретных динамических систем по незашумленным наблюдениям // Статистические проблемы управления. Вильнюс: Ин-т матем. и киберн. АН ЛитССР, 1986. Вып. 73. С. 89-99.