МОРФОЛОГИЧЕСКИЕ СИСТЕМЫ и АНАЛИЗ ИЗОБРАЖЕНИЙ

Содержание

Слайд 2

Что такое морфология? СПОСОБЫ - эвристики, эксперименты НАУКА МЕТОДЫ - математические

Что такое морфология?

СПОСОБЫ - эвристики, эксперименты
НАУКА
МЕТОДЫ - математические

модели
- формализованные критерии
- решения, обладающие доказанными свойствами
- оптимальность
МАШИННОЕ ЗРЕНИЕ - весь комплекс проблем, связанных с
получением пространственой инфорамции,
включая сенсоры, вычислители и алгоритмы
КОМПЬЮТЕРНОЕ ЗРЕНИЕ - математические и
алгоритмические аспекты машинного зрения
МАТЕМАТИЧЕСКОЕ ЗРЕНИЕ
ФОТОГРАММЕТРИЯ МОРФОЛОГИЯ
(геометрия простр.распред. данных) (модели данных и процедур)
Слайд 3

Что такое морфология? Термин: Морфология – (1) «наука о форме»; (2)

Что такое морфология?

Термин: Морфология – (1) «наука о форме»;
(2) методы анализа

изображений, основанные на содержательных яркостно-геометрических моделях и критериях.
Источники: Морфология Серра, Морфология Пытьева.
Обобщение 1: Морфологический анализ – схема анализа данных, которая в качестве обязательного этапа предполагает обоснованное (в некотором смысле оптимальное) построение модельного описания гипотетического (скрытого) прообраза наблюдаемых данных (сегментация + реконструкция).
Обобщение 2: "Морфология" или "морфологическая система"- это такой формализм анализа данных (изображений), в котором любые образы (изображения) рассматриваются как элементы некоторого пространства (алгебры), любые задачи формулируются в терминах этого пространства, и операции осуществляются над элементами этого пространства (целыми изображениями), а не над отдельными пикселями.
Слайд 4

Формальная морфология ϑ Λ M ⊆ ϑ ε δ ϕεδ =

Формальная морфология

ϑ

Λ

M ⊆ ϑ

ε

δ

ϕεδ = δε

множество образов

множество описаний

модельное
множество

Морфологическая сегментация ε:

ϑ → Λ
Морфологическая реконструкция δ: Λ → ϑ
Морфологический фильтр ϕεδ(E)=δ(ε(E)): ϑ→Λ→ϑ
Слайд 5

Сегментация + Реконструкция множество образов множество описаний Искусственный изоморфизм Естественный гомоморфизм

Сегментация + Реконструкция

множество образов

множество описаний

<5,6,8>
<1,0,1,0,0>
<2,0,0,1,0,…>

Искусственный
изоморфизм

Естественный
гомоморфизм

568=5×102+6×101+8×10
40=1×25+1×23
28=22 × 71
f(x)


a4 x4+a3 x3+a2 x2+a1 x+a0

Естественный
изоморфизм

Позиционные системы счисления

Разложение на простые множители

Аппроксимация полиномами

Слайд 6

Пример 1. Морфологический анализ Пытьева Сравнение по форме, выделение отличий ─

Пример 1. Морфологический анализ Пытьева Сравнение по форме, выделение отличий


=


=

Алгоритм сравнения

изображений по форме:
Выделить связные области на изображении A.
Вычислить среднюю яркость по областям A на B.
Сформировать C по форме A с яркостями из B.
Найти разность С и B.

A

B

B

C

форма A и С

форма

яркость

(реконструкция B)

Слайд 7

Морфологический анализ Пытьева A1 A2 A3 A4 Индикаторная функция множества Ai

Морфологический анализ Пытьева

A1

A2

A3

A4

Индикаторная функция множества Ai

Значение яркости(интенсивности) пикселя с координатами

(x, y) (функция изображения)

с1,…, с4 – яркости областей A1,…, A4 фона и граней куба

X

Модель изображения - кусочно-постоянная функция

X – поле зрения

Слайд 8

Морфологический анализ Пытьева Проекция изображения на форму Изображение , определенное на

Морфологический анализ Пытьева

Проекция изображения на форму

Изображение , определенное на поле

Х

Форма изображения V(f) в виде множества

- интеграл яркостей по области Ai

Проекция вектора g на плоскость V(f)

Определение коэффициентов ci*

Дифференцируя по ci, получим решение задачи в виде

- площадь области Ai

Слайд 9

Морфологический анализ Пытьева α = arccos ku, β = arccos km

Морфологический анализ Пытьева

α = arccos ku, β = arccos km

Нормированный
коэфициент корреляции
Морфологический
коэффициент корреляции:
0

≤ km ≤ 1
Km не зависит от Морфологический проектор преобразования яркости F(f(x,y)).
Сравнение форм:

Описание формы

>

Слайд 10

Морфологический анализ Пытьева

Морфологический анализ Пытьева

Слайд 11

Пример 2. Математическая морфология Серра Обработка с учетом формы, выделение деталей

Пример 2. Математическая морфология Серра Обработка с учетом формы, выделение деталей


MIN
эрозия

MAX
дилатация

MIN
эрозия

MAX
дилатация

реконструкция A

A

Слайд 12

Математическая морфология Серра B BT T Структурирующий элемент Исходный образ Трансляция

Математическая морфология Серра

B

BT

T

Структурирующий элемент Исходный образ

Трансляция

Opening (открытие)

Closing (закрытие)

XoB =

∪{Bz | Bz ⊆ X}
Слайд 13

Математическая морфология Серра Трансляция A по z: Az = {y| a∈A,

Математическая морфология Серра

Трансляция A по z:
Az = {y| a∈A, y=a=z}.
Сложение

Минковского (дилатация):
A⊕B = {a=b| a∈A, b∈B} =
= U{Ba} = U{Ab}
Вычитание Минковского (эрозия):
AB = {z| Bz ⊆ A} = U{Az}

Оператор, сохраняющий включение: X⊆Y ⇒ Ψ(X) ⊆ Ψ(Y)
Экстенсивный оператор: Ψ(X)⊇X Антиэкстенсивный оператор: Ψ(X)⊆X
Усиливающий оператор (Ψ(Ψ(X))⊇Ψ(X))
Ослабляющий оператор (Ψ(Ψ(X))⊆Ψ(X))

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

Открытие X по B:
XoB = (XB)⊕B = U{Bz| Bz ⊆ X}.
Закрытие X по B:
X∙B = (X⊕B)B

Слайд 14

Математическая морфология Серра ММ-операторы: ММ-проекторы: Эрозия (сжатие) Дилатация (расширение) ММ-открытие ММ-закрытие

Математическая морфология Серра

ММ-операторы:

ММ-проекторы:

Эрозия (сжатие)

Дилатация (расширение)

ММ-открытие

ММ-закрытие

ММ-фильтрация

Учет формы путем выбора структурирующих элементов:

Слайд 15

Математическая морфология Серра

Математическая морфология Серра

Слайд 16

Пример 3. Бинарная морфология на базе скелетов A Фигурой называется связная

Пример 3. Бинарная морфология на базе скелетов

A

Фигурой называется связная замкнутая

область плоскости, ограниченная конечным числом неперсекающихся жордановых кривых.
Пусть P – евклидова плоскость, D(p,r) – открытый круг радиуса r с центром в точке p.
Пустым или вписанным кругом фигуры A называется круг D(p,r)⊆A. Максимальным пустым кругом называется пустой круг, который не содержится целиком ни в одном другом пустом круге данной фигуры. Скелетом S(A) фигуры A называется множество центров всех ее максимальных пустых кругов.
Радиальной или дистанционной функцией rA(p) точки p∈P для фигуры A называется максимальная величина радиуса пустого круга с центром в данной точке.

S(A)

Слайд 17

Бинарная морфология на базе скелетов

Бинарная морфология на базе скелетов

Слайд 18

ПРОЕКТИВНОСТЬ

ПРОЕКТИВНОСТЬ

Слайд 19

Проекторы как распознающие операторы (М. Павель) Структурный фильтр: процедура преобразования образа

Проекторы как распознающие операторы (М. Павель)

Структурный фильтр: процедура преобразования образа к виду,

соответствующему заданному классу структур.
Алгебраический проектор:
F(X)=F(F(X))

Геометрическая интерпретация:

Два способа описания класса:
Проектор – оператор, ставящий в соответствие любому образу образ из модельного множества.
Модель (модельное множество) – множество стабильных элементов проектора.

Слайд 20

Сравнение форм по сложности (Пытьев) M1 M2 M3 M3 ⊆ M2

Сравнение форм по сложности (Пытьев)

M1

M2

M3

M3 ⊆ M2 ⊆ M1

Морфологическая сложность:


Если одно модельное множество целиком принадлежит другому,
то соответствующая форма изображения не сложнее (проще).

>

>

Морфология Пытьева

Морфология Серра


>

>



Структурная сложность:
Чем больше элементов в модели, тем сложнее описание.

Слайд 21

Морфологический спектр (Maragos) r ∂ || A o D(r) || / ∂r монотонное убывание сложности

Морфологический спектр (Maragos)

r

∂ || A o D(r) || / ∂r

монотонное

убывание сложности
Слайд 22

АЛГЕБРА ИЗОБРАЖЕНИЙ

АЛГЕБРА ИЗОБРАЖЕНИЙ

Слайд 23

Морфологические алгебры Общие особенности рассмотренных методов Проектор в алгебраическом смысле Проекция

Морфологические алгебры

Общие особенности рассмотренных методов

Проектор в алгебраическом смысле
Проекция образа

на образ имеет линейный вид: Pr(A,B)=r(A,B)*B
Проекция образа на набор образующих представляет собой комбинацию проекций на каждую образующую

Оператор проекции - фильтр со следующими свойствами:

1, если B⊆A
r(A,B) =
0, если B⊄A

Pr(A,B)=∅

Пример - ММ Серра

B

A

Pr(A,B)=B

B

A

Слайд 24

Морфологические алгебры

Морфологические алгебры

Слайд 25

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

Пример. Проективные морфологии на базе методов интерполяции

Интерполяционные многочлены
Пусть на отрезке

[a,b] заданы n+1 опорных (узловых) точек: a≤x0In(x) = ∑i f(xi)∙L(xi,x),
L(xi,x) = ((x-x0)…(x-xi-1)(x-xi+1)…(x-xn)) /
((xi-x0)…(xi-xi-1)(xi-xi+1)…(xi-xn)).
Слайд 26

Пример. Проективные морфологии на базе методов интерполяции Проективная морфология на базе

Пример. Проективные морфологии на базе методов интерполяции

Проективная морфология на базе

интерполяционных многочленов
Проективная морфология:
A=Vk=1..n(akEk),
Pr(A,E) = Vk=1..n(Pr(A,Ek)) = Vk=1..n(r(A,Ek)∙Ek).
Морфология на базе интерполяции:
A=f(x); E={Ek}={L(xk,x)}; Pr(A,E) = In(x); V=∑;
rk(A,Ek) = f(xk); Prk(A,Ek) = rk(A,Ek)∙Ek = f(xk)∙L(xk,x).
Ортогональный базис:
∀Ei,Ek∈E, k≠i: Prk(Ei,Ek)=0, Pri(Ek,Ei)=0.
Слайд 27

Пример. Проективные морфологии на базе методов интерполяции Кусочно-линейная интерполяция функции на отрезке I(x)=f(xi)(x-xi)+f(xi+1)(x-xi+1).

Пример. Проективные морфологии на базе методов интерполяции

Кусочно-линейная интерполяция функции на

отрезке
I(x)=f(xi)(x-xi)+f(xi+1)(x-xi+1).
Слайд 28

Пример. Проективные морфологии на базе методов интерполяции Морфология на базе кусочно-линейной

Пример. Проективные морфологии на базе методов интерполяции

Морфология на базе кусочно-линейной

проекции
I(x) = ∑i f(xi)∙L(xi,x),
L(xi,x) = {(x-xi-1)/(xi-xi-1): x∈[xi-1,xi]; (x-xi+1)/(xi-xi+1): x∈[xi,xi+1];
0: x∉[xi-1,xi+1]}.
Вид опорной функции:
Слайд 29

Пример. Проективные морфологии на базе методов интерполяции Вид опорной функции L(xi,yi,x,y) = {0: (x,y)∉P(xi,yi); ((x-xj)(y-yk)-(y-yj)(x-xk))/((xi-xj)(yi-yk)-(yi-yj)(xi-xk)):(x,y)∈Tijk∈P(xi,yi)}.

Пример. Проективные морфологии на базе методов интерполяции

Вид опорной функции L(xi,yi,x,y) =

{0: (x,y)∉P(xi,yi); ((x-xj)(y-yk)-(y-yj)(x-xk))/((xi-xj)(yi-yk)-(yi-yj)(xi-xk)):(x,y)∈Tijk∈P(xi,yi)}.
Слайд 30

Пример. Проективные морфологии на базе методов интерполяции Проективная морфология на базе

Пример. Проективные морфологии на базе методов интерполяции

Проективная морфология на базе

двумерной кусочно-линейной интерполяции:
I(x,y) = ∑i f(xi,yi)∙L(xi,yi,x,y),
A=f(x,y); E={Ek}={L(xk,yk,x,y)};Pr(A,E) = I(x,y); V=∑;
rk(A,Ek) = f(xk,yk); Prk(A,Ek) = rk(A,Ek)∙Ek = f(xk,yk)∙L(xk,yk,x,y).
Слайд 31

Пример. Проективные морфологии на базе методов интерполяции Морфологический коэффициент корреляции: K(g,f)

Пример. Проективные морфологии на базе методов интерполяции

Морфологический коэффициент корреляции:
K(g,f) =

min(||Pr(g,f)||,||g||) / max(||Pr(g,f)||,||g||),
Pr(g,f) = ∑i g(xi,yi)∙L(xi,yi,x,y)
«Форма» двумерной функции:
Z={∑i ai∙L(xi,yi,x,y): ai∈R},
{(xi,yi)}=lextr2(f(x,y)),
lextr2(f(x,y)) – множество локальных экстремумов кусочно-линейной функции f(x,y).
Слайд 32

Пример. Проективные морфологии на базе методов интерполяции Сравнение по форме: Проекция на сходную форму

Пример. Проективные морфологии на базе методов интерполяции

Сравнение по форме: Проекция

на сходную форму
Слайд 33

Пример. Проективные морфологии на базе методов интерполяции Сравнение по форме: Проекция на отличающуюся форму:

Пример. Проективные морфологии на базе методов интерполяции

Сравнение по форме: Проекция

на отличающуюся форму:
Слайд 34

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

Пример. Проективные морфологии на базе методов интерполяции

Форма двумерной функции

Вычисление интерполяционной проекции

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

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

Пример. Проективные морфологии на базе методов интерполяции

Сравнение по форме двумерных функций

Примеры

вычисления интерполяционной проекции на кусочно-линейную форму:
форма, определяемая триангуляцией по значимым экстремумам эталона;
эталонное изображение объекта и проекция эталона на форму;
тестовое изображение объекта и проекция тестового изображения на форму;
изображение другого объекта и проекция другого объекта на форму эталона.
Слайд 36

ПРОСТРАНСТВО РАЗЛОЖЕНИЙ

ПРОСТРАНСТВО РАЗЛОЖЕНИЙ

Слайд 37

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

Проективные морфологические разложения

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

описаний этих образов является обоснованным.
Отсюда и все полезные практические свойства таких разложений.
Слайд 38

Проективные морфологические разложения Типы морфологических разложений Условие разложимости: ∃ E⊆Ω: Pr(A,E) = Vk=1..n(Pr(A,Ek)) = Vk=1..n(r(A,Ek)∙Ek)

Проективные морфологические разложения

Типы морфологических разложений
Условие разложимости:
∃ E⊆Ω: Pr(A,E) = Vk=1..n(Pr(A,Ek)) =

Vk=1..n(r(A,Ek)∙Ek)
Слайд 39

Проективные морфологические разложения Ортогональные морфологии на базе БПФ и т.п. БПФ,

Проективные морфологические разложения

Ортогональные морфологии на базе БПФ и т.п.

БПФ, ДКП, фильтры

НЧ, ВЧ

Вейвлет-преобразование и фильтры

Слайд 40

Проективные морфологические разложения

Проективные морфологические разложения

Слайд 41

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

Проективные морфологические разложения Морфологический анализ изображений

Переход от образов к изображениям (двумерным

функциям):
Введем пространства параметров изображения P и разложения Q:
A→A(p); Ek→Ek(p)→ϕ(p,q); E→E(p,q).
Морфологические разложения изображений:
Морфо-геометрическая проекция:
Pr(A(p),E(p,q))=Vq∈Q(A(q)∙ϕ(p,q)).
Морфо-геометрическое разложение:
dec(A(p))=A(q): Ω(P)→Θ(Q).
Проекция разложения на разложение:
Pr(A(q),B(q))=r(A(p),B(p))∙B(q).
Нормированный коэффициент линейной корреляции разложений:
K(A(q),B(q)) = ||Pr(A(q),B(q))||/||A(q)||.
Слайд 42

Проективные морфологические разложения Морфологический анализ изображений Фильтрация изображений с использованием разложений:

Проективные морфологические разложения Морфологический анализ изображений

Фильтрация изображений с использованием разложений:

Слайд 43

Проективные морфологические разложения Морфологический анализ изображений Структурное сравнение изображений (обобщение методики

Проективные морфологические разложения Морфологический анализ изображений

Структурное сравнение изображений (обобщение методики Ю.П.

Пытьева)
Структурный проектор = Морфологический фильтр, применяемый к образу A, область пропускания которого согласована с образом B.
Характеристический базис образа B:
Eχ(B)={χ(bk)∙Ek, Ek∈E},
χ(x)={0, если x=0; 1 – в противном случае}, где E –исходный базис, χ(x) - индикатор структурной связи.
Морфологическая проекция образа A на модель образа [B]:
Pr(A,[B]) = Vk=1..n(ak∙χ(bk)∙Ek) = Pr(A,Eχ(B)).
Морфологическая проекция разложений:
Pr(a,[b]) = Pr({ak},[{bk}]) = {ak∙χ(bk)}.
Структурный морфологический коэффициент корреляции:
Kстр(A,B)= ||Pr(a,[b])|| / ||a||, где A,B∈Ω; a=dec(A),b=dec(B)∈Θ, со стандартными свойствами:
(a) 0 ≤ Kстр(A,B) ≤ 1; (b) Kстр(A,A) = 1; (c) Kстр(A,B) = 0 ⇔ Pr(A,[B]) = ∅.
Класс морфологически эквивалентных структур:
В={X∈Ω: Kстр(X,B)=1}.
Отношение «более простой/более сложный по структуре»:
(Kстр(A,B) = 1, Kстр(B,A) < 1) ⇔ («A сложнее B», «B проще A»).
Слайд 44

Проективные морфологические разложения Конструирование алгоритмов обнаружения объектов Обнаружение объектов с использованием разложений:

Проективные морфологические разложения Конструирование алгоритмов обнаружения объектов

Обнаружение объектов с использованием разложений:


Слайд 45

Проективные морфологические разложения Морфологические операторы сегментации и сжатия данных Морфологическая сегментация на базе проективных разложений

Проективные морфологические разложения Морфологические операторы сегментации и сжатия данных

Морфологическая сегментация на

базе проективных разложений
Слайд 46

Проективные морфологические разложения Проективная сегментация без потерь Морфология Пытьева «форма» Пытьева

Проективные морфологические разложения Проективная сегментация без потерь

Морфология Пытьева «форма» Пытьева

Морфология Серра

морфологический скелет

Минимальное число областей

Минимальное число дисков

Полное писксельное разбиение

Полное дисковое представление

Слайд 47

МОДУЛЬНОСТЬ (в поисках нетривиального описания)

МОДУЛЬНОСТЬ (в поисках нетривиального описания)

Слайд 48

Формальная морфология ϑ Λ M ⊆ ϑ ε δ ϕεδ =

Формальная морфология

ϑ

Λ

M ⊆ ϑ

ε

δ

ϕεδ = δε

множество образов

множество описаний

модельное
множество

Морфологическая сегментация ε:

ϑ → Λ
Морфологическая реконструкция δ: Λ → ϑ
Морфологический фильтр ϕεδ(E)=δ(ε(E)): ϑ→Λ→ϑ
Слайд 49

Способ описания: Преобразование Хафа и GHT Голосование точек в аккумулятор Анализ

Способ описания: Преобразование Хафа и GHT

Голосование точек в аккумулятор

Анализ аккумулятора

Параметризация:

Обобщение:

НТ

GНТ

Пространство параметров

Слайд 50

Модульные проективные морфологии

Модульные проективные морфологии

Слайд 51

Монотонные проективные морфологии на базе преобразования Хафа и GHT H-открытие -

Монотонные проективные морфологии на базе преобразования Хафа и GHT

H-открытие - объединение

проекций изображения A(p) на отдельные прямые линии: Pr(A(p),t) = MAXq∈Q(A(q,t)∙Pr(A(p),ϕ(p,q))) = MAXq∈Q(A(q,t)∙A(p)∙ϕ(p,q)),
где p=(x,y); q=(ρ,θ) – параметры нормальной параметризации прямой; Q – пространство параметров; ϕ(p,q)∈{0,1} – характеристическая функция прямой с параметрами q; A(q,t)∈{0,1} – аккумулятор преобразования Хафа, бинаризованный по порогу t.
(а) (b) (с)
Пример морфологического H-открытия: a – исходное бинарное изображение;
b – аккумулятор пространства Хафа c – результат H-открытия.
На исходном контурном препарате выделены глобальные прямолинейные структуры.
Аналогичным образом строится монотонная проективная морфология на базе
обобщенного преобразования Хафа (GHT).
Слайд 52

Морфологическая фильтрация на базе рекуррентного преобразования Хафа в скользящем окне Вычислительно-эффективная

Морфологическая фильтрация на базе рекуррентного преобразования Хафа в скользящем окне
Вычислительно-эффективная реализация алгоритма.
(a) (b)

(c)
Пример морфологического RHT-открытия. Выделены локальные прямолинейные структуры.
a) (b) (c)
Пример морфологической RHT-фильтрации с различными параметрами размера окна.
Выдеелены линеаменты различных размеров.
Слайд 53

МОДУЛЬНОСТЬ (комбинирование процедур)

МОДУЛЬНОСТЬ (комбинирование процедур)

Слайд 54

Альтернативные модульные морфологии Идея: Построение различных модульных морфологических операторов путем комбинирования

Альтернативные модульные морфологии

Идея: Построение различных модульных морфологических операторов путем комбинирования

разных операторов сегментации с разными операторами реконструкции.
Пример: селективные морфологии на базе операторов ММ Серра
Слайд 55

Селективные морфологии (a) (b) (c) (a) (b) (a) (a) Открытие MM

Селективные морфологии

(a)

(b)

(c)

(a)

(b)

(a)

(a)

Открытие MM

(a)

(a)

O(Object)=

∅, if E(Object) = ∅ (a)

Object’⊆Object, if E(Object)=∅

(b)

Object, if ∃Im:O(Im)=Object (c)

∅, if E(Object) = ∅ (a)

Object, if E(Object) ≠ ∅ (b)

SO(Object)=

Стандартная ММ

Селективная ММ

Открытие SM

Эрозия MM

ε

δ

δ′

Слайд 56

Селективные морфологии Im E(Im) O(Im) Im-O(Im) Полутоновое MM-открытие

Селективные морфологии

Im

E(Im)

O(Im)

Im-O(Im)

Полутоновое
MM-открытие

Слайд 57

Селективные морфологии Полутоновое SM-открытие Im E(Im) SO(Im) Im-SO(Im)

Селективные морфологии

Полутоновое
SM-открытие

Im

E(Im)

SO(Im)

Im-SO(Im)

Слайд 58

Селективные морфологии Im D(Im) C(Im) Im-C(Im) Полутоновое MM-закрытие

Селективные морфологии

Im

D(Im)

C(Im)

Im-C(Im)

Полутоновое
MM-закрытие

Слайд 59

Селективные морфологии Im D(Im) SC(Im) Im-SC(Im) Полутоновое SM-закрытие

Селективные морфологии

Im

D(Im)

SC(Im)

Im-SC(Im)

Полутоновое
SM-закрытие

Слайд 60

Селективные морфологии Полутоновое MM-открытие Im E(Im) O(Im) Im-O(Im)

Селективные морфологии

Полутоновое
MM-открытие

Im

E(Im)

O(Im)

Im-O(Im)

Слайд 61

Селективные морфологии Полутоновое SM-открытие Im E(Im) SO(Im) Im-SO(Im)

Селективные морфологии

Полутоновое
SM-открытие

Im

E(Im)

SO(Im)

Im-SO(Im)

Слайд 62

Селективные морфологии Полутоновое MM-закрытие Im D(Im) C(Im) Im-C(Im)

Селективные морфологии

Полутоновое
MM-закрытие

Im

D(Im)

C(Im)

Im-C(Im)

Слайд 63

Селективные морфологии Полутоновое SM-закрытие Im D(Im) SC(Im) Im-SC(Im)

Селективные морфологии

Полутоновое
SM-закрытие

Im

D(Im)

SC(Im)

Im-SC(Im)

Слайд 64

Селективные морфологии Im E1D(Im) SO1D(Im) Im-SO1D(Im) Контурная селективная морфология (на базе

Селективные морфологии

Im

E1D(Im)

SO1D(Im)

Im-SO1D(Im)

Контурная селективная морфология
(на базе оператора удаления заданного числа концевых точек)

Слайд 65

МОРФОЛОГИЧЕСКИЙ АНАЛИЗ НА БАЗЕ КРИТЕРИЕВ (сегментация с регуляризацией)

МОРФОЛОГИЧЕСКИЙ АНАЛИЗ НА БАЗЕ КРИТЕРИЕВ (сегментация с регуляризацией)

Слайд 66

Критерии в классических морфологиях Морфология Пытьева: Pr(A,M(B)) = argminC∈M(B) || A

Критерии в классических морфологиях
Морфология Пытьева:
Pr(A,M(B)) = argminC∈M(B) || A -

C ||
Морфология Серра:
AoB = argminC∈M(B) {|| A - С ||: C⊆A}
или
AoB = argmaxC∈M(B) {||С ||: C⊆A}

M

A

C∈M

|| A - C ||

B

A

AoB

B∈M

Слайд 67

Критериальная морфология: Модель: M(λ): Λ→[0,1] ⇔ M(L):ϑ→[0,1] Критерий соответствия: K(E,λ): ϑ×Λ→[0,1]

Критериальная морфология:
Модель: M(λ): Λ→[0,1] ⇔ M(L):ϑ→[0,1]
Критерий соответствия:
K(E,λ): ϑ×Λ→[0,1] ⇔ K(E,L): ϑ×ϑ→[0,1]
Критериальный

морфологический фильтр ϕФ на базе (ε,δ):
εФ(E)=λ, ϕФ(E)=δ(λ): Ф(E,λ)=K(E,δ(λ))×M(λ)→max(λ∈Λ)
ℑ(ℜ)={ϑ,Λ,δ,K,M}⇔ℑ′(ℜ)={ϑ,Λ,εФ,δ} – ℑ-морфология.
Проективные критериальные морфологии: ϕФ(E)=ϕФ(ϕФ(E)).

Морфологический подход к анализу данных

Слайд 68

Морфологическое решение задач анализа данных: 1. Фильтрация ϕФ(E)=δ(λ): 2. Сегментация εФ(E)=λ:

Морфологическое решение задач анализа данных:
1. Фильтрация ϕФ(E)=δ(λ):
2. Сегментация εФ(E)=λ:
Ф(E,λ)=K(E,δ(λ))×M(λ)→max(λ∈Λ)
3. Распознавание cФ(E)=H:
Ф(E,λ,H)=K(E,δ(λ))×M(λ,H)×M(H)→max(λ∈Λ,H∈Θ)
4.

Обнаружение/локализация επФ(E)=λ:
Параметрическая выборка π(E,λ): ϑ×Λ→ϑ,
Фπ(E,λ,H)=K(π(E,λ),δ(λ))×M(λ,H)×M(H)→max(λ∈Λ,H∈Θ)
селективный морфологический фильтр
ϕπ(E)=π(E,επФ(E)): ϑ×Λ→ϑ.
Вывод: морфологический подход позволяет единым унифицированным способом решать все основные задачи обработки и анализа данных.

Морфологический подход к анализу данных

Слайд 69

Нечеткие модели: [0,1] Максимум достоверности: Ф(A,L)=K(A,L)×M(L)→max(L∈Ω) Вероятностные модели: [0,1] Максимум апостериорной

Нечеткие модели: [0,1]
Максимум достоверности:
Ф(A,L)=K(A,L)×M(L)→max(L∈Ω)
Вероятностные модели: [0,1]
Максимум апостериорной вероятности
ψ(A)=L: P(A,L)=P(A/L)×P(L)→max(L∈Ω).
Четкие или логические

модели: [0,1] → {0,1}.
Морфологическая проекция на модельное множество:
ψ(A,M): K(A,L)→max(L∈M), M = {B∈Ω: M(B)=1}.
Теоретико-информационные критерии: [0,1] → [0,+∞)
Максимум энтропии (минимум информации):
Ф(A,L)=J(A,L)+α×Q(L)→min(L∈Ω) ⇔ K(A,L)×M(L)α→max(L∈Ω).
J(A,L) = – log(P(A/L)); Q(L) = – log(P(L)); α - модельный параметр
Интерпретация: регуляризация задачи сегментации по Тихонову

Форма и семантический смысл критериев

Слайд 70

Регуляризация Ф(A,L) = J(A,L) → min(L∈F(x)) x f(x)

Регуляризация

Ф(A,L) = J(A,L) → min(L∈F(x))

x

f(x)

Слайд 71

Регуляризация Ф(A,L) = J(A,L) → min(L∈F(x)) x f(x)

Регуляризация

Ф(A,L) = J(A,L) → min(L∈F(x))

x

f(x)

Слайд 72

Регуляризация Ф(A,L)=J(A,L)+α×Q(L)→min(L∈F(x)) x f(x)

Регуляризация

Ф(A,L)=J(A,L)+α×Q(L)→min(L∈F(x))

x

f(x)

Слайд 73

Регуляризация ⇒ сегментация с потерями Ф(A,L)=J(A,L)+α×Q(L)→min(L∈F(x)) x f(x)

Регуляризация ⇒ сегментация с потерями

Ф(A,L)=J(A,L)+α×Q(L)→min(L∈F(x))

x

f(x)

Слайд 74

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

Критериальные проективные морфологии

Пусть имеется множество образов Ω, на котором определена операция

сложения (‘+’), задающая на Ω группу с «нулевым образом» ∅. Кроме этого, на множестве образов определена Ω – норма μ(A)=||A||: Ω→R, ||∅||=0, причем норма разности обладает свойствами расстояния. На множестве пар образов задана
Функция-критерий (критерий штрафа)
Ф(A,B): Ω×Ω→R
Морфологический проектор на базе критерия:
Pr(A,Ф)=B: Ф(A,B)→min(B∈Ω), Pr(A)=Pr(Pr(A)), Pr(∅)=∅.
Критериальная морфологическая модель -
М = {A∈Ω: Pr(A,Ф)=A}
множество собственных (стабильных) элементов проектора. Модель M2 по отношению к M1 является более сложной, если M2 ⊆ M1.
Морфологический коэффициент корреляции:
KМ(A,Pr)=KM(A,M)=exp( - ||A-Pr(A,M)||/||Pr(A,M)||),
0≤KM(A,M)≤1; KM(A,M)=1 ⇔ A∈M; Pr(A,M)=∅ ⇔ KM(A,M)=0.
Конкретные морфологии определяются конкретным видом критерия.
Слайд 75

Стандартный критерий штрафа Ф(A,B)= J(A,B) + χ(A,B) + α×Q(B) где J(A,B)

Стандартный критерий штрафа
Ф(A,B)= J(A,B) + χ(A,B) + α×Q(B)
где J(A,B) – критерий

соответствия проекции и образа, причем
∀A∈Ω, B∈V(A,Ф): J(A,A)≤J(A,B),
χ(A,B) – критерий (предикат) допустимости решения, определяющий ОДЗ
χ(A,B) = {0: B∈V(A,Ф); +∞: B∉V(A,Ф)},
Q(B) – критерий качества проекции, характеризующий ее принадлежность модели M;
α≥0 – структурирующий параметр, обеспечивающий компромисс между требованиями соответствия и качества.
Утверждение. С увеличением значения структурирующего параметра α сложность модели, которую определяет проектор, монотонно убывает.
⇒ α - параметр морфологической сложности модели.
Морфологический спектр:
Sp(A,α) = ∂ || Pr(A,J,χ,α,Q) || / ∂α
Коэффициент максимальной морфологической сложности:
αmax(A)=max{α≥0: A=Pr(A,J,χ,α,Q)}.

Критериальные проективные морфологии

Слайд 76

Морфологический спектр по параметру r ∂ || A o D(r) || / ∂r монотонное убывание сложности

Морфологический спектр по параметру

r

∂ || A o D(r) || / ∂r

монотонное

убывание сложности
Слайд 77

Морфологический спектр по сложности α Sp(A,α) = ∂ || Pr(A,J,Q) ||

Морфологический спектр по сложности

α

Sp(A,α) = ∂ || Pr(A,J,Q) || / ∂α

монотонное

убывание сложности

Ф(A,B) = J(A,B) + α×Q(B)

Слайд 78

Достаточные условия построения проективных операторов Ф(A,B)= J(A,B) + χ(A,B) + α×Q(B) Критериальные проективные морфологии

Достаточные условия построения проективных операторов
Ф(A,B)= J(A,B) + χ(A,B) + α×Q(B)

Критериальные проективные

морфологии
Слайд 79

Алгоритмические аспекты морфологии (от математического зрения к компьютерному)

Алгоритмические аспекты морфологии (от математического зрения к компьютерному)

Слайд 80

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

Среднеквадратичная проективная сегментация одномерных функций

Пример СКО-фильтрации. Исходная функция, результаты сегментации для α=500 и

α=2000.

Критериальные проективные морфологии

Слайд 81

Монотонная проективная фильтрация одномерных функций Примеры применения операторов DP-Open (α=200) и

Монотонная проективная фильтрация одномерных функций

Примеры применения операторов DP-Open (α=200) и DP-Close(α=200).

Примеры

применения операторов DP-Open (α=1000) и DP-Close(α=1000).

Критериальные проективные морфологии

Слайд 82

Монотонная проективная сегментация одномерных функций Примеры применения операторов DP-Open: Исходная функция,

Монотонная проективная сегментация одномерных функций

Примеры применения операторов DP-Open: Исходная функция, DP-Open

(α=1000), DP-Open (α=50000)

Примеры применения операторов DP-Close: Исходная функция, DP-Close (α=1000), DP-Close (α=50000).

Критериальные проективные морфологии

Слайд 83

α=700 α=600 α=500 α=400 α=300 α=200 α=100 α=0 Монотонная фильтрация и


α=700
α=600
α=500
α=400
α=300
α=200
α=100
α=0

Монотонная фильтрация и сегментация двумерных кривых (контуров бинарных изображений)

Пример кусочно-линейной

сегментации типа «закрытие».

Критериальные проективные морфологии

Слайд 84

α=600 α=500 α=400 α=0 α=100 α=200 α=300 α=700 Монотонная фильтрация и


α=600
α=500
α=400
α=0
α=100
α=200
α=300
α=700

Монотонная фильтрация и сегментация двумерных кривых (контуров бинарных изображений)

Пример кусочно-линейной

сегментации типа «открытие».

Критериальные проективные морфологии

Слайд 85

Морфологии на базе сопоставления функций 1D-Matching: g(x) = f(x+v(x)) + ξ

Морфологии на базе сопоставления функций
1D-Matching:
g(x) = f(x+v(x)) + ξ
J(f(x),g(x),L(x))=J(f(x),g(x+L(x))).
2D-Matching:
g(x,y) = f(x+vx(x,y),y+vy(x,y))

+ ξ,
L(x,y)={Lx(x,y),Ly(x,y)}
Ф(f(x,у),g(x,y),L(x,у))= J(f(x,у),g(x,y),L(x,у)) +αQ(L(x,у)) →min(L),
J(f(x,y),g(x,y),L(x,у))= ∑x∑y(f(xi,yi)-g(xi+Lx(x,y),yi+Ly(x,y)))2
L(x) - функция-решение = функция диспаратности
Pr(f,g)=g(xi+Lx(x,y),yi+Ly(x,y)) – морфологический проектор

Критериальные проективные морфологии

Слайд 86

Морфология на базе оптимальной кусочно-линейной интерполяции Ф(f,L) = -J(L) + χ(f,L)

Морфология на базе оптимальной кусочно-линейной интерполяции
Ф(f,L) = -J(L) + χ(f,L)

+ α×Q(L) → min(L),
где J(L) – длина графика ломаной; Q(L) – число узловых точек.
Множество опорных точек – переменная оптимизации.
Квазимонотонная ОДЗ: V(f)⊄V(L)

Критериальные проективные морфологии

Слайд 87

Морфология на базе оптимальной кусочно-линейной интерполяции a) b) c) a –

Морфология на базе оптимальной кусочно-линейной интерполяции
a) b) c)
a – исходная

функция; b – результат сегментации α=5; с – результат сегментации α=15.

Критериальные проективные морфологии

Слайд 88

α=1000 α=900 α=800 α=600 α=400 α=200 α=0 Морфология на базе оптимальной


α=1000
α=900
α=800
α=600
α=400
α=200
α=0

Морфология на базе оптимальной кусочно-линейной интерполяции двумерных кривых (контуров бинарных изображений)


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

Критериальные проективные морфологии

Слайд 89

д) α=1000 г) α=100 в) α=10 б) α=5 а) α=0 Морфология

д) α=1000

г) α=100

в) α=10

б) α=5

а) α=0

Морфология на базе оптимальной кусочно-постоянной сегментации

двумерных полутоновых изображений

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

Динамическое программирование: 2-D как 1-D

Слайд 90

д) α=1000 г) α=100 в) α=10 б) α=5 а) α=0 Морфология

д) α=1000

г) α=100

в) α=10

б) α=5

а) α=0

Морфология на базе оптимальной кусочно-постоянной сегментации

двумерных полутоновых изображений

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

Динамическое программирование: 2-D как 1-D

Слайд 91

Срезовый стек: Оператор реконструкции срезового стека δ: Морфологическая фильтрация с использованием стековых представлений

Срезовый стек:

Оператор реконструкции срезового стека δ:

Морфологическая фильтрация с использованием стековых

представлений
Слайд 92

г) в) б) а) з) ж) е) Морфологическая фильтрация с использованием

г)

в)

б)

а)

з)

ж)

е)

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

д)

Стековое свойство элементов срезовой цепи изображения:
а)

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

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

Оператор уменьшающей срезовой фильтрации:

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

Слайд 94

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

Дерево срезовых сегментов:

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

Слайд 95

1. Морфологический фильтр, сохраняющий структуру стекового дерева 2. Морфологический фильтр с

1. Морфологический фильтр, сохраняющий структуру стекового дерева

2. Морфологический фильтр с

преобразованиями стекового дерева:

Удаление ветвей

Слияние узлов

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

Слайд 96

з) α=10000000, n=1. ж) α=1000000, n=2; е) α=100000, n=4; д) α=10000,

з) α=10000000, n=1.

ж) α=1000000, n=2;

е) α=100000, n=4;

д) α=10000, n=12;

Морфологическая фильтрация с

использованием стековых представлений

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

г) α=1000, n=33;

в) α=100, n=45;

б) α=10, n=53;

а) α=0, n=56;

||A - B|| + α×n(B) → min(B: B ≤ A)

Слайд 97

з) α=10000000, n=1. ж) α=1000000, n=1; е) α=100000, n=4; д) α=10000,

з) α=10000000, n=1.

ж) α=1000000, n=1;

е) α=100000, n=4;

д) α=10000, n=12;

г) α=1000, n=32;

в)

α=100, n=45;

б) α=10, n=53;

а) α=0, n=56;

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

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

||A - B|| + α×n(B) → min(B: A ≤ B)

Слайд 98

з) α=10000000, n=1. ж) α=1000000, n=3; е) α=100000, n=8; д) α=10000,

з) α=10000000, n=1.

ж) α=1000000, n=3;

е) α=100000, n=8;

д) α=10000, n=24;

г) α=1000, n=120;

в)

α=100, n=371;

б) α=10, n=686;

а) α=0, n=1584;

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

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

||A - B|| + α×n(B) → min(B: B ≤ A)

Слайд 99

з) α=10000000, n=1. ж) α=1000000, n=2; е) α=100000, n=12; д) α=10000,

з) α=10000000, n=1.

ж) α=1000000, n=2;

е) α=100000, n=12;

д) α=10000, n=38;

г) α=1000, n=156;

в)

α=100, n=480;

б) α=10, n=813;

а) α=0, n=1395;

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

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

||A - B|| + α×n(B) → min(B: A ≤ B)

Слайд 100

НЕПРОЕКТИВНЫЕ МОРФОЛОГИИ МОРФОЛОГИЧЕСКИЙ АНАЛИЗ СВИДЕТЕЛЬСТВ (от математического зрения к компьютерному)

НЕПРОЕКТИВНЫЕ МОРФОЛОГИИ МОРФОЛОГИЧЕСКИЙ АНАЛИЗ СВИДЕТЕЛЬСТВ (от математического зрения к компьютерному)

Слайд 101

Основной современный подход к задаче обнаружения и распознавания – обнаружение и

Основной современный подход к задаче обнаружения и распознавания – обнаружение и

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

Составление яркостно-геометрической конструкции (модели объекта) – неформальный элемент алгоритмизации
Поиск и локализация объекта распознавания – формальный элемент алгоритмизации

Схема модельного подхода к обнаружению объектов

робастные алгоритмы
голосования или
сопоставления (matching)

сцена

сенсоры

изображения

сегментация (выделение)
характерных черт

Неформальная
структурированная
яркостно-геометрическая
модель объекта

МОДЕЛЬНЫЙ ПОДХОД к построению эффективных процедур обнаружения объектов

Слайд 102

ХЧ Точки Углы, соединения линий, точки высокой кривизны градиента яркости, центр

ХЧ

Точки

Углы, соединения линий, точки высокой кривизны градиента яркости, центр тяжести области,

концы линий, точки экстремальных значений признаков

Линии

Прямые или криволинейные структуры, границы областей

Области

Сегментированные области, специфические формы (эллипсы, прямоугольники и т.д.)

Структуры

Комбинации ХЧ

Характерные элементы (черты), используемые
в иерархических алгоритмах обнаружения

Слайд 103

Характерные черты на изображении имеют следующие виды атрибутов: Положение: Концы отрезка,

Характерные черты на изображении имеют следующие виды атрибутов:
Положение: Концы отрезка, центр

отрезка, центр тяжести области, вершины многоугольников;
Геометрические атрибуты: Ориентация, длина, кривизна, площадь, периметр, ширина линии, минимальный и максимальный диаметр области, оси симметрии, число и положение особых точек, показатель компактности;
Радиометрические атрибуты: Контраст, статистика распределения яркости, знак и величина края, автокорреляция;
Текстурные атрибуты: Матрица смежности, показатель однородности, энергия, энтропия, статистика градиентов текстуры, результаты применения текстурных фильтров, моменты;
Топологические атрибуты: Связность, соседство, общие точки, пересечение, параллельность, перекрытие, включение;
Цветовые/многозональные атрибуты: вектор атрибутов для каждого канала;
Динамические атрибуты: атрибуты статических и движущихся объектов;
Временные атрибуты: функции изменения атрибутов со временем.

Свойства (атрибуты) характерных черт

Слайд 104

Примеры атрибутов характерных черт ‘.’ - близость; ‘!’ - коллинеарность; ‘=‘

Примеры атрибутов характерных черт

‘.’ - близость;
‘!’ - коллинеарность;
‘=‘ - параллельность;
‘L’- перпендикулярность;

‘<‘ – угол.

Выбранный отрезок
Другие линии

Исходная информация

Слайд 105

Присутствие/плотность: наличие данных ХЧ на всех используемых изображениях, достаточная плотность ХЧ

Присутствие/плотность: наличие данных ХЧ на всех используемых изображениях, достаточная плотность ХЧ

для покрытия интересующего района;
Редкость/Уникальность: редкость конкретной ХЧ на изображении, уникальность ХЧ в окрестности;
Инвариантность/Устойчивость: робастность по отношению к геометрическим и радиометрическим искажениям, нечувствительность к шуму;
Локализация: возможность точной локализации;
Интерпретация: возможность быстрого распознавания и интерпретации;
Скорость: время выделения данного класса ХЧ из исходного изображения.

Критерии выбора ХЧ и их атрибутов для построения алгоритмов обнаружения

Слайд 106

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

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

Слайд 107

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

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

Контурный препарат

Слайд 108

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

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

Выделение характерных черт - областей

поддержки линеаментов
Слайд 109

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

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

Выделение первичных линеаментов

Слайд 110

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

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

Выделение крупных и/или параллельных линеаментов

Слайд 111

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

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

Фильтрация по размеру

Слайд 112

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

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

Поиск П-образности

Слайд 113

Пример автоматического обнаружения зданий с летательного аппарата Окончательная 2D-модель

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

Окончательная 2D-модель

Слайд 114

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

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

Межкадровое соответствие

Слайд 115

Окончательная модель

Окончательная модель

Слайд 116

Морфологический анализ свидетельств Вероятностная интерпретация методов морфологического анализа изображений Вероятностная модель

Морфологический анализ свидетельств

Вероятностная интерпретация методов морфологического анализа изображений
Вероятностная модель формирования образа
P(M): Ω→[0,1],
Вероятностная

модель регистрации изображения
P(L/M): M→[0,1],
Вероятностная модель искажений
P(A/L): M×Ω→[0,1]
Критерий максимальной вероятности
P(A,L)=P(A/L)×P(L/M)×P(M)→max(L).
Оператор максимально вероятной реконструкции образа
ψ: Ω→M, ψ(A)=L: P(A,L)→max(L).
Слайд 117

Морфологический анализ свидетельств Вероятностное обобщение понятия проективности При вторичном применении оператора

Морфологический анализ свидетельств

Вероятностное обобщение понятия проективности
При вторичном применении оператора максимально вероятной

реконструкции необходимо учитывать изменения в вероятностной модели, происходящие при переходе от реального образа к модельному:
P(L,L)=P(L/L)×P(L/M)×P(M)→max(L),
причем переход L→L не вносит искажений.
Пусть σ - численный параметр мощности искажений. Тогда:
∀A≠L,σ>0: P(A/L,σ)>P(A/L,0);
∀σ≠0 ∃A≠L: P(A/L)>P(L/L) ⇒ проективность может быть нарушена.
Параметрическая модель:
P(A,L,σ)=P(A/L,σ)×P(L/M)×P(M),
ψ(σ): Ω→M, ψ(A,σ)=L: P(A,L,σ)→max(L).
Требование детерминированной проективности
ψ(A)=ψ(ψ(A))
Требование вероятностной проективности
ψ(A,σ)=ψ(ψ(A,σ),0)
Слайд 118

Морфологический анализ свидетельств Признаки как достаточные статистики. Независимые признаки Набор признаков:

Морфологический анализ свидетельств

Признаки как достаточные статистики. Независимые признаки
 Набор признаков: f(A)=∈Ψn,
Критерий

максимальной вероятности:
P(f,L)=P(f/A)×P(A/L)×P(L/M)×P(M)→max(L). Вероятностная модель измерений:
P(f/L)=P(f/A)×P(A/L),
P(f,L)=P(f/L)×P(L/M)×P(M)→max(L). Оператор максимально вероятной реконструкции образа по признакам
ψf: Ω→M, ψf(A)=L: P(f,L)→max(L). Вероятностная обоснованность системы признаков в задаче морфологического анализа:
P(f,L)→max(L) ⇒ P(A,L)→max(L),
∀A∈Ω: ψf(A)=ψ(A), P(L/A)=P(L/f(A)), то есть f(A) является достаточной статистикой для A.
Независимость признаков: существует факторизация
P(f,L)=Πi=1,…,n P(fi,L).
Слайд 119

Морфологический анализ свидетельств Анализ морфологических свидетельств Морфологическое событие: e(p)={f(A,p)=e∈X}. Морфологическая гипотеза:

Морфологический анализ свидетельств

Анализ морфологических свидетельств
Морфологическое событие:
e(p)={f(A,p)=e∈X}.
Морфологическая гипотеза:
h(q)={ψ(A) = L(q)}.
Модель голосования:
P(E(A),h(q))=Πp∈X, P(e(p),h(q))→max(h(q)),
где

E(A) – совокупность морфологических событий или точнее совокупное морфологическое событие, связанное с образом A; h(q)⊆H(Θ), H(Θ) – пространство морфологических гипотез; P(e(p),h(q)) - вероятностная модель морфологического голосования.
Носитель гипотезы (множество влияющих событий):
S(h(q))={e(p): ∃h′(q)≠h(q): P(e(p),h(q))≠P(e(p),h′(q))}.
Носитель события (множество влияющих гипотез):
S(e(p))={h(q): ∃e′(p)≠e(p): P(e(p),h(q))≠P(e′(p),h(q))}.
Полную группу событий, относящихся к одному признаку f(p), будем называть доменом событий, полную группу гипотез, соответствующую различным значениям L(q) – доменом гипотез.
Слайд 120

Морфологический анализ свидетельств Анализ морфологических свидетельств Под анализом морфологических свидетельств понимается

Морфологический анализ свидетельств

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

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

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

Метод анализа морфологических свидетельств

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

Слайд 121

РАЗРАБОТКА ВЫЧИСЛИТЕЛЬНО ЭФФЕКТИВНЫХ АЛГОРИТМОВ АНАЛИЗА МОРФОЛОГИЧЕСКИХ СВИДЕТЕЛЬСТВ (от математического зрения к компьютерному)

РАЗРАБОТКА ВЫЧИСЛИТЕЛЬНО ЭФФЕКТИВНЫХ АЛГОРИТМОВ АНАЛИЗА МОРФОЛОГИЧЕСКИХ СВИДЕТЕЛЬСТВ (от математического зрения к

компьютерному)
Слайд 122

Способы повышения вычислительной эффективности: • независимое аккумулирование свидетельств • декомпозиция вектора

Способы повышения вычислительной эффективности:
•  независимое аккумулирование свидетельств
•  декомпозиция вектора параметров S(θ)=S'(θ')∙S"(θ")


•  редукция вектора параметров S(θ)→S'(θ')
•  загрубление модели объекта M → M'⊇M
•  иерархический анализ свидетельств
Модульная схема алгоритма обнаружения
•  обработка изображения по схеме голосования с целью выделения объектов или их составляющих
•  анализ аккумулятора с целью определения положения и/или ориентации объектов
•  повторный анализ изображения с целью проверки природы обнаруженных объектов и уточнения их параметров

Морфологический анализ свидетельств

Слайд 123

Иерархический анализ свидетельств. Теорема разделения Морфологический анализ свидетельств

Иерархический анализ свидетельств. Теорема разделения

Морфологический анализ свидетельств

Слайд 124

Последовательность шагов разработки алгоритма обнаружения и идентификации объектов 1. описать модели

Последовательность шагов разработки алгоритма
обнаружения и идентификации объектов
1.  описать модели объекта, регистрации и

искажений
2.  определить степень загрубления модели объекта
3.  осуществить необходимую редукцию параметров
4.  определить типы «событий»
5.  составить качественную вероятностную модель
6.  определить процедуру голосования
7.  определить соответствующую процедуру анализа аккумулятора
8. разработать процедуру постпроверки достоверности детектирования

Морфологический анализ свидетельств

Слайд 125

Пример. Метод обнаружения штриховых кодов и текстовых областей на изображениях

Пример. Метод обнаружения штриховых кодов и текстовых областей на изображениях

Слайд 126

Пример. Метод обнаружения штриховых кодов и текстовых областей на изображениях Обнаружение

Пример. Метод обнаружения штриховых кодов и текстовых областей на изображениях

Обнаружение штриховых

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

Модель объекта. Прямоугольная область плоскости, заполненная белыми и черными полосами. Ширина полос может быть различной.
Модель регистрации. Геометрическая - проективные преобразования с малыми углами (≈аффинные). Радиометрическая - отображение черно-белой палитры на шкалу серого цвета (0-15).
Модель искажений. По яркости - блики, перечеркивание, замещение фрагментов. По геометрии - изгибы и "коробление" несущей поверхности.

Обнаружение штриховых кодов. Разработка алгоритма

1. Загрубление модели. Прямоугольная область плоскости (x1,у1, x2,у2, x3,у3, x4,у4) с высокой колинеарностью градиентов.
2. Загрубление и декомпозиция вектора параметров. Модель "четыре угла" преобразуется к модели "четыре прямые" (ρ1,θ,ρ2,θ+90°,ρ3,θ,ρ4,θ+90°). Отсюда декомпозиция модели на кодосодержащую полосу (ρ1,ρ2,θ) и положение в полосе (ρ3, ρ4).
3. Аккумулирование свидетельств. Голоса пикселей с высоким градиентом аккумулируются в двухпараметрическом пространстве Hough (ρ,θ) с учетом направления градиента.

Слайд 127

ρ Параметризация: Голосование точек в аккумулятор Параметризация x cosθ + y

ρ

Параметризация:

Голосование точек в аккумулятор

Параметризация

x cosθ + y sinθ = ρ

ПРЕОБРАЗОВАНИЕ ХАФА

x

y

x

y

ϴ

ρ

ϴ

ρ

ϴ0

ρ0

.

Достоинства:

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

Идея модифицированного преобразования Хафа Метод выделения штриховых и текстовых строк в

Идея модифицированного преобразования Хафа

Метод выделения штриховых и текстовых строк в пространстве

Хафа

Исходное изображение и голосование с учетом направления градиента

Область с выделенными машиносчитываемыми строками

Соответствующий аккумулятор пространства Хафа

MHT для обнаружения штриховых кодов и текстовых строк

1

2

3

4

1

2

3

4

Обрабатывается полутоновое изображение
2. Объекты поиска – полосы (группы прямых)
3. Голосуют точки с высоким модулем градиента яркости
4. Голосование с учетом направления градиента в точке
5. Анализ аккумулятора – поиск сигнала заданной формы

Особенности метода

x

y

ϴ

ρ

x

y

ϴ

ρ

Слайд 129

АВТОМАТИЗИРОВАННОЕ КОНСТРУИРОВАНИЕ МОРФОЛОГИЧЕСКИХ АЛГОРИТМОВ (от математического зрения к компьютерному)

АВТОМАТИЗИРОВАННОЕ КОНСТРУИРОВАНИЕ МОРФОЛОГИЧЕСКИХ АЛГОРИТМОВ (от математического зрения к компьютерному)

Слайд 130

Автоматизированное конструирование алгоритмов обнаружения объектов Постановка задачи P1(R,A) → min(R,A) |

Автоматизированное конструирование алгоритмов обнаружения объектов Постановка задачи

P1(R,A) → min(R,A) | P2(R,A)

≤ P2max, T(R,A) ≤ Tmax
P1 – вероятность необнаружения объекта;
P2 – вероятность ложной тревоги;
T – вычислительная стоимость алгоритма (время, ресурсы);
R – используемая морфологическая система;
A – алгоритм анализа данных.
Слайд 131

Автоматизированное конструирование алгоритмов обнаружения объектов 1. Метод автоматизированного конструирования алгоритмов обнаружения

Автоматизированное конструирование алгоритмов обнаружения объектов

1. Метод автоматизированного конструирования алгоритмов обнаружения объектов,

основанный на преобразованиях модельных описаний
Функции системы
Метод преобразования модельных описаний
Рекурсивные модели и алгоритмы
Нерекурсивные модели и алгоритмы
Проективные морфологии на базе неоднородных структурных моделей, описываемых логическими предикатами
2. Метод автоматизированного конструирования модульных процедур обнаружения объектов, основанный на «генетическом отборе» элементов модельного описания
Общий подход к построению процедур идентификации
Учет информативности опорных элементов
Построение процедур идентификации объектов нескольких классов
Генетический отбор морфологических процедур
Слайд 132

Автоматизированное конструирование алгоритмов обнаружения объектов Конструирование детекторов по модельным описаниям Формальное

Автоматизированное конструирование алгоритмов обнаружения объектов Конструирование детекторов по модельным описаниям

Формальное описание

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

Метод преобразования модельных описаний Модель объекта: Преобразования моделей: перестановка порядка предикатов;

Метод преобразования модельных описаний
Модель объекта:
Преобразования моделей:
перестановка порядка предикатов; ,
декомпозиция (разбиение)

модели на две части и редукция (отсечение) одной из них.
Obj = { | }

Автоматизированное конструирование алгоритмов обнаружения объектов Конструирование детекторов по модельным описаниям

,

Слайд 134

Пример преобразования модельных описаний Модели: Штриховая линия = набор штрихов, лежащих

Пример преобразования модельных описаний
Модели: Штриховая линия = набор штрихов, лежащих на одной

прямой. (М1) Штриховая линия = прямая, состоящая из отдельных штрихов. (М2)
Процедуры: Найти все штрихи, выбрать те, что лежат на одной прямой. (П1) Последовательно находить штрихи, лежащие на одной прямой. (П2)

Автоматизированное конструирование алгоритмов обнаружения объектов Конструирование детекторов по модельным описаниям

Слайд 135

Пример преобразования модельных описаний Рекурсивная модель (M2): Штрих_линия(X)=лежать_на_прямой(х) рядом(х,у)Штрих_линия(Х\х); Штрих_линия({х})=лежать_на_прямой(х). После

Пример преобразования модельных описаний Рекурсивная модель (M2):
Штрих_линия(X)=лежать_на_прямой(х) рядом(х,у)Штрих_линия(Х\х);
Штрих_линия({х})=лежать_на_прямой(х). После перестановки предикатов:
Штрих_линия(X)=лежать_на_прямой(х) Штрих_линия(Х\х) рядом(х,у);
Штрих_линия({х})=лежать_на_прямой(х). После

декомпозиции:
Нерекурсивная модель (M2′):
Штрих_линия1(X)=лежать_на_прямой(х) Штрих_линия1(Х\х);
Штрих_линия1({х})=лежать_на_прямой(х) ( )
Рекурсивная модель (M2′′):
Штрих_линия2(X)=рядом(х,у) Штрих_линия2(Х\х);
Штрих_линия2({х})=ИСТИНА.
Смешанная модель:
M1’=(М2′,M2′)

Автоматизированное конструирование алгоритмов обнаружения объектов Конструирование детекторов по модельным описаниям

Слайд 136

Алгоритм применения построенной модели процедуры голосования: 1. Осуществить все возможные успешные

Алгоритм применения построенной модели процедуры голосования:
1.  Осуществить все возможные успешные индексации

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

Автоматизированное конструирование алгоритмов обнаружения объектов Конструирование детекторов по модельным описаниям

Слайд 137

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

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

(достаточное условие построения проективного морфологического фильтра на базе логической модели): Морфологическое преобразование на базе модели M является морфологическим проектором, если выполняется условие
∀q: A(q)=0 ⇒ ∀A′(q)≠0: M(A(Q)VA(q))=M(A(Q)VA′(q)
(то есть в пользу M(A(Q)) голосуют только ненулевые элементы A(q)).
Морфологический проектор:
Pr(A(p),M)=Φ(A(p),M)=Φ(Φ(A(p),M),M)
Морфологический коэффициент корреляции изображения с моделью:
KM(A(p),M)=min(||Pr(A(p),M)||,||A(p)||) / max(||Pr(A(p),M)||,||A(p)||)

Автоматизированное конструирование алгоритмов обнаружения объектов Конструирование детекторов по модельным описаниям

Слайд 138

Принцип конструирования процедуры идентификации Множество процедур обнаружения: , где - процедура

Принцип конструирования процедуры идентификации
Множество процедур обнаружения:
, где - процедура обнаружения фрагмента,

реализующая один из заданных базовых алгоритмов (j=1,2,…,mk; k=1,2,3…).

Автоматизированное конструирование алгоритмов обнаружения объектов Метод генетического отбора структурных моделей

Слайд 139

Задача условной оптимизации: , где: - время работы процедуры на изображениях

Задача условной оптимизации:
,
где: - время работы процедуры на изображениях из

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

Автоматизированное конструирование алгоритмов обнаружения объектов Метод генетического отбора структурных моделей

Слайд 140

Схема применения генетического алгоритма: 1. Ген = одна из элементарных процедур.

Схема применения генетического алгоритма:
1.     Ген = одна из элементарных процедур.
2.     Хромосома

= последовательность генов ограниченной длины.
3.     Функция качества хромосомы:
,
,
где р - процедура обнаружения заданного объекта; Si - изображение из обучающей выборки; - время работы процедуры р на изображении Si; А – настроечный коэффициент; - штрафная функция.
4.     Операция скрещивания – перегруппировка и обмен составных частей существующих решений (цепочек процедур обнаружения).
5.     Операция мутации позволяет изменить параметры (xj,yj,wj,hj) для выбранной элементарной процедуры.
6.     Генетический отбор осуществляется путем итеративного «размножения», тестирования и селекции в каждом поколении хромосом с наилучшим значением функции качества. При этом на каждом этапе случайным образом осуществляются мутации параметров и скрещивание моделей.

Автоматизированное конструирование алгоритмов обнаружения объектов Метод генетического отбора структурных моделей

Слайд 141

Схема применения генетического алгоритма для формирования морфологических детекторов: 1. Ген =

Схема применения генетического алгоритма для формирования морфологических детекторов:
1.     Ген = один

из возможных структурных примитивов, характеризуемый набором {Mk(u,qk),tk,qk}.
2.     Хромосома = последовательность генов = морфо-геометрическая модель объекта M(p,u).
3.     Функция качества хромосомы - аналогично.
4.     Операция скрещивания – аналогично.
5.     Операция мутации позволяет изменить параметры локализации {Mk(u,qk),qk} для выбранного элемента модели.
6.     Генетический отбор - аналогично.

Автоматизированное конструирование алгоритмов обнаружения объектов Метод генетического отбора структурных моделей

Слайд 142

Интерпретация результата: 1. Процедурная интерпретация = Близкая к оптимальной процедура обнаружения

Интерпретация результата:
1.     Процедурная интерпретация = Близкая к оптимальной процедура обнаружения заданного

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

Автоматизированное конструирование алгоритмов обнаружения объектов Метод генетического отбора структурных моделей

Слайд 143

РАЗНООБРАЗИЕ МОРФОЛОГИЧЕСКИХ СИСТЕМ

РАЗНООБРАЗИЕ МОРФОЛОГИЧЕСКИХ СИСТЕМ

Слайд 144

- Проективные / непроективные - Модульные / унитарные - Алгебраические /

- Проективные / непроективные
- Модульные / унитарные
- Алгебраические / критериальные
- Монотонные

/ ортогональные
- Вложения / заполнения
- Аппроксимирующие / интерполирующие
- Дискретные / непрерывные
- На базе однородных / неоднородных структур
- На базе независимого проецирования
- На базе независимого голосования
- На базе логического программирования
- На базе динамического программирования
- На базе линейного программирования…

Какие бывают морфологии:

Слайд 145

Типы описаний, используемых в анализе изображений

Типы описаний, используемых в анализе изображений

Слайд 146

Ф(A,λ)=K(e(A),δ(λ))×M(λ)→max(λ∈Λ). Выбор ℑ(ℜ)={ϑ,Λ,δ,K,M} определяет конкретную морфологическую систему. Возникает нечто вроде систематической

Ф(A,λ)=K(e(A),δ(λ))×M(λ)→max(λ∈Λ).

Выбор ℑ(ℜ)={ϑ,Λ,δ,K,M} определяет конкретную морфологическую систему.
Возникает нечто вроде систематической таблицы методов

анализа изображений. Часть клеток заполнена, часть – еще нет.
Вывод: на наш век морфологий хватит:)

Типы описаний. Систематическая таблица морфологий

Слайд 147

ПРИМЕР ИСПОЛЬЗОВАНИЯ ТАБЛИЦЫ Морфологическое сравнение бинарных изображений на базе их скелетных

ПРИМЕР ИСПОЛЬЗОВАНИЯ ТАБЛИЦЫ
Морфологическое сравнение бинарных изображений на базе их скелетных представлений


Основная

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

Морфологическое сравнение бинарных изображений на базе их скелетных представлений Пусть дан

Морфологическое сравнение бинарных изображений на базе их скелетных представлений

Пусть дан образ B

и требуется сравнить его форму с формой образа A.
Морфологическая проекция бинарного образа B на форму образа A:
Pr(B,A) = ∪p∈S(A) { G(p, rB(p)) }.
Нормированный морфологический коэффициент корреляции (МКК) Пытьева:
KM(B,A) = || Pr(B,A) || / || B ||.
МКК оценивает степень сходства формы фигуры B с формой фигуры A.
Слайд 149

Морфологическое сравнение бинарных изображений на базе их скелетных представлений Примеры проекции

Морфологическое сравнение бинарных изображений на базе их скелетных представлений


Примеры проекции фигуры

на форму

Pr(B,A) = ∪p∈S(A) { G(p, rB(p)) }

Слайд 150

Морфологическое сравнение бинарных изображений на базе их скелетных представлений Проекция на

Морфологическое сравнение бинарных изображений на базе их скелетных представлений


Проекция на класс

форм «расширение/сжатие»

Pr(B, M(A,b)) = A ⊕ bB,
bB = arg max b∈R {||A ⊕ b||: (A ⊕ b) ⊆ B}.

Слайд 151

Морфологическое сравнение бинарных изображений на базе их скелетных представлений Проекция на

Морфологическое сравнение бинарных изображений на базе их скелетных представлений


Проекция на класс

форм «утончение/утолщение»

Pr(B, M(A,a)) = A ⊗ aB,
aB = arg max a∈R {||A ⊗ a||: (A ⊗ a) ⊆ B}.

Слайд 152

Морфологическое сравнение бинарных изображений на базе их скелетных представлений Инвариантные проекции

Морфологическое сравнение бинарных изображений на базе их скелетных представлений


Инвариантные проекции на

разные классы форм

Максимальный по норме вписанный образ в заданном модельном классе:
PR(B, M(A,Σ)) = argmax p∈P ||Pr(B, M(A(p),Σ))||

Слайд 153

Спасибо за внимание!

Спасибо за внимание!