Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Luciv / МКОИ_Lab

.pdf
Скачиваний:
68
Добавлен:
21.03.2016
Размер:
1.6 Mб
Скачать

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qk = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.17)

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µ = qk µk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.18)

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично (4.10) дисперсия межкластерных отклонений равна:

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ2B = qk (µk − µ)2

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.19)

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя уравнения (4.17) и (4.18), межкластерная дисперсия (4.19)

может быть записана как:

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

2 − µ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ2B = qk µk

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.20)

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку не зависит от от выбора порогов, то для выполнения

критерия (4.16) вместо σ2

 

достаточно анализировать только σ2 :

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

b

 

M

 

 

 

 

2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ2

= q

k

µ

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.21)

b

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Согласно критерию, чтобы найти оптимальные пороги, диапазон поиска

для максимальной

σ2 составляет 1 ≤

f <L-M+1,

f + 1 ≤ f

2

<L-M+2, …, и

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

1

 

 

1

 

 

f M 1 + 1≤

 

f M 1<(L-1).

На рис. 4.2

показаны зоны поиска порогов в

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

 

 

 

 

 

 

 

 

 

 

 

...

fM 1

 

 

 

 

L-1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L-(M-1)+2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f1

 

 

 

 

f2

 

 

 

 

 

 

 

 

 

L-(M-1)+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L-(M-1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

 

2

 

 

3

 

 

 

 

 

 

L-1

 

 

 

 

 

1

 

f1

f2

 

 

 

 

 

 

fM 1

 

L

 

 

Рисунок 4.2 Соотношение порогов в методе полного перебора.

Каждое следующее значение порогового значения должно быть хотя бы на 1 больше предыдущего. Это условие формирует левую границу диапазона. После каждого i-го порога должно быть установлено L-M+i следующих пороговых значений. Это условие определяет правую границу возможного диапазона порога.

Нахождение оптимального выбора порогов выполняется методом полного перебора.

21

4.2 Порядок выполнения работы

1)Разработать алгоритм и код программы квантования на два уровня при определении порога по методу Оца.

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

3)Синтезировать тестовые изображения

4)Выполнить сегментацию тестовых изображений.

5)Оценить эффективность сегментации.

6)Проанализировать полученные результаты, составить отчет о выполненных исследованиях, сделать выводы по работе.

4.3 Библиографический список

1)В.Т. Фисенко, Т.Ю. Фисенко, Компьютерная обработка и распознавание изображений: учеб. пособие.- СПб: СПбГУ ИТМО, 2008.-192 с.

2)Гонсалес, Р. Цифровая обработка изображений / Р.Гонсалес, Р.Вудс - М.:

Техносфера, 2005.-1072 с.

3)A Fast Algorithm for Multilevel Thresholding. Ping-Sung Liao, Tse-Sheng Chen* and Pau-Choo Chung. Department of Electrical Engineering. National Cheng Kung University.Tainan, 701 Taiwan. P. 15

4.4 Вопросы для самопроверки:

1)Какие методы пороговой обработки Вы знаете?

2)Какой критерий положен в основу порогового метода Оцу?

3)Какие допущения сделаны при разработке метода кластеризации Оцу?

4)Нарисуйте схему алгоритма Оцу?

5)Как связаны критерии минимума внутрикластерных отклонений и максимума межкластерных расстояний?

6)Какой параметр используется для оценки внутрикластерных отклонений и как он вычисляется?

7)Опишите схему алгоритма многопороговой обработки.

8)Какой критерий положен в основу метода многопороговой обработки?

9)Как осуществляется поиск оптимального набора порогов при многопороговой обработке?

22

5Лабораторная работа № 5

Исследование методов морфологической сегментации

Цель работы – ознакомиться с методами морфологической сегментации на примере метода морфологической реконструкции и метода выращивания областей.

Задание по работе:

Изучить теоретическую часть работы.

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

Разработать алгоритм выращивания областей и написать код программы.

Выполнить алгоритмы морфологической сегментации.

Сделать выводы по выполненной работе.

5.1 Теоретическая часть

Сегментация – одна из самых важных проблем в задачах классификации изображений [1]. Рассмотрим методы морфологической сегментации на примере двух методов: метода реконструкции и метода выращивания областей.

5.1.1Метод реконструкции

Метод реконструкции [2] применяется в тех случаях, когда в бинарном изображении необходимо выделить связные области, помеченные маркерами. Бинарное изображение, содержащее области, назовем изображением масок g. Маркеры представлены в изображении, назовем его изображением маркеров f, в виде областей, принадлежащих объекту интереса так, что маркер является подмножеством масок:

f g .

(5.1)

Изображение маркеров может быть получено при анализе характеристик объекта, отличных от тех, по которым были выделены маски. Например, в изображении g выделены области одного тона. Известно, что изображения, имеющие один тон, но разные яркости, воспринимаются ЗС человека как разные. В качестве маркера, может быть использовано бинарное изображение 2, полученное по компоненту яркости и последующему применению операции логического «И» с изображением 1.

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

 

1

 

 

1

1

 

 

 

 

 

B = 1

1

1 .

(5.2)

1

1

1

 

 

 

 

 

Алгоритм реконструкции выполняется как итеративная процедура: Инициализация:

23

присвоить h1 изображение маркеров f: h1=f

задать cтруктурирующий элемент B.

Далее выполняются итерации алгоритма:

h k +1

= (hk B)Ig .

(5.3)

Критерий останова:

 

h k +1

= h k .

(5.4)

На рис. 5.1 Представлен пример морфологической реконструкции.

маска g

маркер f

результат

 

1

 

1

1

B = 1

1

1

1

1

1

 

 

 

а) б) в) г)

Рисунок 5.1 Морфологическая реконструкция. а) изображение маски g; б) изображение контуров областей (масок), точками внутри контуров представлены маркеры; в) изображение маркеров f; г) выходное изображение.

5.1.2Метод выращивания областей

Второй метод – метод выращивания областей,- состоит в наращивании областей от некоторых заданных маркеров [3] («центров кристаллизации») путем добавления соседних элементов, имеющих заданные свойства. Исходными данными для алгоритма являются:

входное полутоновое изображение f,

множество «центров кристаллизации» S,

множество пороговых значений T.

Множество «центров кристаллизации» S и множество пороговых значений T в общем случае могут представлять собой изображения, размер которых равен размеру сегментируемого изображения f, и трактоваться как множество пороговых значений для формирования изображений масок и

24

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

Алгоритм выращивания областей следующий. Формирование массива центров кристаллизации SI (маркера):

 

если f = S

1,

SI =

(5.5)

0,

иначе

Формирование бинарного массива TI (изображение масок), в котором единичным элементам соответствуют элементы изображения f, в которых выполняется условие однородности:

 

если abs( f S ) T

1,

TI =

(5.6)

0,

иначе

Выполнение морфологической реконструкции, для которой маркером является SI, а маской TI. Структурирующий элемент B задается в соответствии с заданным критерием связности.

Выполнение селекции связных компонентов и получение изображения связных компонентов g и их количества.

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

Вкачестве примера выполнена сегментация изображения при S=98, T=8, как это показано на рис. 5.2.

На рис. 5.2 б) представлено изображение центров кристаллизации, полученное простым выбором тех отсчетов, в которых яркость изображения равна 98. На рис. 5.2 в) представлено изображение масок, в котором единичные значения присвоены элементам изображения, в которых модуль разности яркости и значения S не превышает порогового значения T. На рис.

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

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

5.2 Порядок выполнения работы

1)Разработать алгоритм и код программы морфологической сегментации областей по методу реконструкции.

2)Разработать алгоритм и код программы сегментации областей по методу выращивания областей.

3)Выбрать тестовые изображения.

4)Разработать алгоритм формирования изображения маркеров.

25

5)Выполнить сегментацию тестовых изображений.

6)Оценить эффективность сегментации.

7)Проанализировать полученные результаты, составить отчет о выполненных исследованиях, сделать выводы по работе.

а) б) в) г)

Рисунок 5.2 Выращивание областей. а) исходное изображение б) изображение центров кристаллизации SI; в) изображение маcок TI; г) выходное изображение.

5.3 Библиографический список

1)В.Т. Фисенко, Т.Ю. Фисенко, Компьютерная обработка и распознавание изображений: учеб. пособие.- СПб: СПбГУ ИТМО, 2008.-192 с.

2)Vincent L. “Morphological Grayscale Reconstruction in Image Analysis: Applications and Efficient Algorithms, “IEEE Trans. on Image Processing”,1973, Vol. 2, № 2, pp.176-201.

3)Гонсалес, Р. Цифровая обработка изображений / Р.Гонсалес, Р.Вудс - М.:

Техносфера, 2005.-1072 с.

5.4 Вопросы для самопроверки:

1)Что такое сегментация изображения?

2)Каким условиям должна удовлетворять сегментация изображения?

3)Какие морфологические методы сегментации областей вы знаете?

4)Для чего предназначена морфологическая реконструкция?

26

5)Как выполняется морфологическая реконструкция?

6)В каких случаях выполняется метод выращивания областей?

7)Как формируется изображение центров кристаллизации?

8)Как формируется изображение областей?

9)Для чего может использоваться алгоритм утончения при выращивании областей?

10)Приведите пример, когда сегментацию областей можно выполнить только с помощью алгоритма выращивания.

11)Какие признаки используются при сегментации областей по исследованным методам?

12)Как оценивается эффективность сегментации?

27

6Лабораторная работа № 6

Исследование морфологических методов изменения бинарных областей

Цель работы – ознакомиться с морфологическими методами изменения бинарных областей. Исследовать алгоритм «утончение» областей.

Задание по работе:

Изучить теоретическую часть работы.

Разработать алгоритм «утончения» областей и написать код программы утончения области.

Выполнить алгоритм морфологического утончения.

Сделать выводы по выполненной работе.

6.1 Теоретическая часть

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

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

Утончение множества A по примитиву B определяется с помощью преобразования «успех/неудача» в соответствии с уравнением:

A B = A

 

(A * B) =A I (A * B)C .

(6.1)

 

Для симметричного утончения используется не один примитив, а множество примитивов{B}:

{B}= {B1,B2 ,..., Bn }.

(6.2)

Изображения примитивов представлены на рис. 6.1. На рис. 6.1 а) представлено подмножество примитивов, полученных при вращении примитиваB1, на рис. 6.1 б) представлено подмножество примитивов, полученных при вращении примитива B2 . Операция утончения по множеству примитивов B производится в соответствии с уравнением:

A {B}= ((...((A B1 ) B2)...) Bn ).

(6.3)

Алгоритм выполняется как итерационный.

 

Каждая итерация состоит из выполнения (6.1) с одним примитивом. К результату этого преобразования применяется (6.1) со следующим примитивом и т.д. в соответствии с (6.3).

28

Критерий останова: результаты двух последовательных итераций совпадают или число итераций превысило заданное число итераций.

центр

 

=-1

 

 

 

= 1

 

 

 

 

 

 

 

 

 

X

X

X

 

 

 

X

 

B1

B3

 

 

X

X

X

 

 

 

X

X

X

B2

X

X

B5

B7

B6

X

= 0

 

X

X

 

B4

 

X

X

 

B8

 

а) б)

Рисунок 6.1 Множество структурирующих элементов, используемых для утончения области.

На рис. 6.2 приведен пример применения алгоритма утончения к области.

Исходное изображение

1 итерация

2 итерация

3 итерация

 

а) б) в) г)

Рисунок 6.2 Преобразование «утончение области» для числа итераций от 1(б) до трех (г).

29

Третья итерация является последней, поскольку изображение на третьей итерации не производит изменений области, сформированной на 2-ой итерации.

6.2 Порядок выполнения работы

1)Изучить теоретическую часть работы.

2)Разработать алгоритм «утончения» областей и написать код программы утончения области.

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

4)Выполнить алгоритм морфологического утончения для разных множеств примитивов

5)Сделать выводы по проделанной работе.

6.3

Библиографический список

1)

В.Т. Фисенко, Т.Ю. Фисенко, Компьютерная обработка и распознавание

 

изображений: учеб. пособие.- СПб: СПбГУ ИТМО, 2008.-192 с.

2)

Гонсалес, Р. Цифровая обработка изображений / Р.Гонсалес, Р.Вудс - М.:

 

Техносфера, 2005.-1072 с.

6.4

Вопросы для самопроверки:

1)Изложите схему алгоритма преобразования «успех/неудача».

2)В чем особенность применения этого алгоритма в методе утончения области?

3)Нарисуйте схему алгоритма утончения области.

4)Какие еще методы утончения области Вы знаете? Опишите их достоинства и недостатки.

5)Какие критерии останова применяются в методе утончения области?

6)Укажите область применения метода утончения

7)Какие методы включают преобразование успех/неудача?

8)Какие преобразования необходимо выполнить после применения метода утончения

9)Какой метод изменения области может быть реализован на использовании метода утончения и как?

30