Распознавание непрерывных комплекснозначных контуров изображений Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ И СИСТЕМЫ

УДК 621.391

Р. Г. Хафизов, С. А. Охотников

РАСПОЗНАВАНИЕ НЕПРЕРЫВНЫХ КОМПЛЕКСНОЗНАЧНЫХ КОНТУРОВ ИЗОБРАЖЕНИЙ

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

Ключевые слова: непрерывный контур, элементарный контур, дискретизация контура, вероятность правильного распознавания.

Введение. В системах технического зрения обработка изображений и распознавание зрительных образов связаны с решением задач контурного анализа. Изображение, подлежащее обработке, вводится в запоминающее устройство, имеющее матричную организацию. Размер клетки (пиксела) сетчатки обычно выбирается в соответствии с возможностью восстановления исходного изображения *(х,у) по оцифрованному ш(п1,^)• Возникающие

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

размеров пиксела, но и от критерия принятия решения о принадлежности каждого из пикселов фону или изображению [1].

Для обработки контура аналитическим или компьютерным способом производится его кодирование, т.е. в соответствие каждому контурному элементу ставится определенное число. Контур Г , состоящий из элементов у(п), п = 0,1,..., * -1, описывается выражением

Г = Мп)}о,,_1 =М0), Н1)^ у(* -1)}, где у(п) — элементарный вектор, соединяющий центры или узлы соседних контурных ячеек сетчатки.

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

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

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

Математическая модель непрерывного контура. Представим контур X = {х(l)}0 L

как непрерывную замкнутую кривую, заданную на комплексной плоскости (рис. 1):

X(l) = Xi (l) + iX2 (l) = |x(l)| exp |/Ф(/)}, l e [0,2n],

где Xi(l) = Re x(l) и X2(l) = Im X(l) — действительная и мнимая компоненты функции x(l);

|x(l) = ^xa(1 ) + X2(l) и ф(1) = argx(l) — амплитудное и фазовое представления функции X(l). Примеры непрерывных контуров изображений приведены на рис. 1.

ImX(/)

4 2 0 -2

-4

X

Im x(l)

10

-5

X

-10

-5

5 Re x(l)

-20 -10

-L*

10 Re x(l)

Рис. 1

При многократном обходе вдоль замкнутой кривой контур X можно представить в виде функции с периодом 2п, т.е. [2]

Х(/) = х(/ + 2па), а = 0,1,2,...

Таким образом, если функция х(/) изменяется, подчинялась тем или иным условиям, например X = {х(/)}о ^ задает контур изображения, то можно говорить о комплексной переменной, заданной на плоскости С комплексной переменной х . Представление контуров изображений в виде функции комплексной переменной позволяет использовать для их анализа и обработки теорию функции комплексной переменной. Так, например, дробно-линейное отображение вида

ау + Ь

v = ■

bc -ad Ф 0,

cy + d'

устанавливает взаимно-однозначное соответствие между точками у -плоскости и точками v -плоскости. Частный случай дробно-линейного отображения вида v = a у + b, где a = ai + ia2 и b = bi + b — произвольные комплексные числа, соответствующий вращению на угол arg a, растяжению в |a| раз и параллельному сдвигу на вектор b , служит для моделирования процесса формирования наблюдаемого контура X(c) путем растяжения исходного контура X в раз, поворота на угол Дф и сдвига начальной точки на d элементарных векторов, т.е.

5

0

0

0

Х(с) ={И- |х (( + й)| ехр {/• [ф (( + й) + Аф]}}0 2п .

Шумовой контур Z = {£(/)}о 2П = {^1 (I) + ^2 (I)}о 2П представляет собой комплексно-значный аналог флуктуационного шума, где (^(0 и (^(0 — действительные случайные функции действительного аргумента /. Закон распределения функций <^(/) и £ 2(1) предполагается нормальным с математическим ожиданием М^ = М^ = М^ = 0 и дисперсией

Зашумленный непрерывный контур К = {к(/)]о 2П образуется как аддитивная совокупность исходного X = {х(/)}0 2п и шумового Z = {£(/)}0 2п контуров, т.е. к(/) = х(/) + С(/),

/ = [0,2п]. Элемент к(/) представляет собой сумму постоянной составляющей детерминированной функции х(/) и случайной функции <^(/), имеющей нулевое математическое ожидание, дисперсию и подчиняющейся нормальному закону распределения. Поэтому функция к(/) будет распределена по нормальному закону с дисперсией Дт и математическим ожиданием М£ =х(/) .

Дискретизация непрерывных контуров изображений. Если X = {х(/2П — произвольный контур, спектр Р = {р(ш)) которого отличен от нуля лишь в полосе частот

-шв <ш <шв , то его можно разложить в обобщенный ряд Фурье по базису Котельникова 8п . С учетом комплекснозначного характера функции X произведем дискретизацию раздельно для действительной и мнимой частей, в результате получим

ю

х(/)= Е с1,п8п(/;®в)+сп8п(/;®в),

п=-ю

где с1;п =(х1 (/),Sn (/;шв)) и с2,п =(х2 (/),8п (/;шв)) — коэффициенты ряда, причем величины с п и е1п — не что иное, как У1 (п) = Х1 (/п ) и У2 (п) = х2 (/п ), т.е. мгновенные значения непрерывного контура х(/) в п-й отсчетной точке /п = птс/шв .

На рис. 2 представлен пример, иллюстрирующий процесс дискретизации непрерывного элементарного контура Хт = {хт (/)}0 2п = { |х|ехр{гда/))02п порядка т = 1. Совокупность

Хт при т = -ю, ю задает семейство замкнутых непрерывных элементарных контуров. При |х| = 1 любые два контура из полного семейства непрерывных элементарных контуров являются ортогональными, т.е. скалярное произведение контуров равно

, \ Ь ( л [2п при т = п,

(Хт, Хп ) = I ехр {(т - п)/)й/ = \

0 [0 при т ^ п.

В результате дискретизации непрерывного элементарного контура Х1 порядка т = 1 получаем дискретный элементарный контур Г порядка т = 1.

На рис. 3 представлены контуры Г, полученные путем дискретизации непрерывных контуров X (см. рис. 1), при этом шв = 5 (рис. 3, а, в) и шв = 10 (рис. 3, б, г).

Xl(l) 1

o -1

X2(/) 1

o -1

1 2 3

1,1

o,55

5 б l

1 2 3

o

-o,55

-1,1

/ — \

( Xl л

V ;

\ /

-1,1 -0,55 0 0,55 1

Рис. 2

а)

Im Г,Х б

4

2

o

-2 -4

б)

Im Г, X

i

б 4 2 o -2

-10 -5

5 Re

Г, X

-4

5 б l

в)

Im Г, X

10

-5

г)

Im Г, X

10

_L

-5

-lo -5

-20 -10

10

Re Г, X Рис. 3

-20 -10

Г1

5 Re

Г, X

■А*

1o Re Г, X

o

4

o

4

o

o

5

5

o

o

o

o

Распознавание непрерывных контуров изображений. Задача распознавания дискретных контуров изображений достаточно подробно рассмотрена в работах [1, 3]. Задача исследования потенциальной эффективности распознавания дискретных комплекснозначных векторных сигналов решалась в работе [4]. В качестве меры различия двух сигналов было выбрано расстояние между задающими их векторами. Помехоустойчивость распознающей системы характеризовалась вероятностью правильного распознавания /^р в предположении, что

величина /р одинакова для всех классов распознаваемых сигналов.

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

Рассмотрим задачу распознавания зашумленных непрерывных контуров на примере

двух классов. Пусть Х^ = {х(1) (/ь и х(2) = {х(2) (/ь — непрерывные контуры соответственно классов и ^2 . На устройство распознавания подается либо зашумленный кон-

тур первого класса К(1) = Х(С)) + Z, либо зашумленный контур второго класса К(2) = Х(2 + Z .

(1) (1)

_ лт(с)

(2) (2)

Задача распознавания заключается в обоснованном отнесении зашумленного контура К к классу Д или А . Для принятия решения необходимо сформировать безусловное отношение правдоподобия и сравнить его с пороговым значением сигнала. Функции правдоподобия контуров К (1) и К (2) имеют соответственно следующий вид:

Л (к, ||да) |, Афа> й(1)) = е^ | Л (, ||(2) |, АФ(2), й(2)) = ®Ф<

1

2а2 1

2а2

ПК2 +

||К||2 +

Х( с) Х(1)

Х(с) Х(2)

- 2Re (, Х^

- 2Re (, Х(2)^

тогда логарифм безусловного отношения правдоподобия определяется как

1П

Л (, ||(1) |, ||(2) |, АФ(1), АФ(2), d(1), й(2))

2а2

Х(с) Х(1)

Х( 1) Х(2)

-2

^ (, Х(1)) )-Re (, Х(2)'

где а — дисперсия шума.

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

На рис. 4 приведены графики зависимости вероятности правильного распознавания /^р

зашумленных непрерывных Хт = {хт(/)} = {exp{iml}}о (сплошные кривые) и дискретных

1

Гт = {хт (и))о 5-1 = {ехр/тп}} 5-1 (штриховые кривые) элементарных контуров от отношения

сигнал/шум q при объемах алфавита Q = 5, 10 и 15.

Величина q равна отношению средней энергии элементарного вектора элементарного контура к дисперсии шумового контура Z :

— для непрерывного контура Хт q = ■

Хт

— для дискретного контура Гт q = •

Г

2па

|2

2

2

Рп

Q = 5 Q = 10 Q = 15

Q = 5 Q = 10 Q = 15

0,01 0,02 0,03 0,04 0,05 0,06 0,07 0,08 0,09 1 q

Рис. 4

Скалярное произведение двух элементарных контуров Хт и Хп или Гт и Гп, являющихся ортогональными сигналами, равно нулю. Поэтому расстояние между нормированными элементарными контурами постоянно. Уменьшение вероятности при распознавании за-

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

Заключение. В статье рассмотрены вопросы представления непрерывных контуров изображений, заданных в комплекснозначном виде. Исследованы вопросы дискретизации непрерывных контуров изображений.

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

Получены сравнительные характеристики вероятности правильного распознавания непрерывных и дискретных элементарных контуров.

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований, проект №11-07-00585а.

2

т

Коммутация структуры пространства линейных устройств

9

СПИСОК ЛИТЕРАТУРЫ

1. Введение в контурный анализ и его приложение к обработке изображений и сигналов / Под ред. Я. А. Фурмана. М.: ФИЗМАТЛИТ, 2002. 592 с.

2. Хафизов Р. Г. Анализ непрерывных комплекснозначных сигналов, задающих контуры изображений плоских объектов // Вестн. КГТУ им. А. Н. Туполева. 2006. № 4. С. 24—27.

3. Егошина И. Л., Хафизов Р. Г. Распознавание плоских зашумленных изображений по их форме // Изв. вузов. Приборостроение. 2006. Т. 49, № 4. С. 46—51.

4. Комплекснозначные и гиперкомплексные системы в задачах обработки многомерных сигналов / Под ред. Я. А. Фурмана. М.: ФИЗМАТЛИТ, 2004. 456 с.

5. Хафизов Р. Г., Охотников С. А. Линейная фильтрация непрерывных контуров изображений, заданных в комплекснозначном виде // Компьютерная оптика. 2010. Т. 34, № 3. С. 408—416.

Сведения об авторах

Ринат Гафиятуллович Хафизов — д-р техн. наук, профессор; Марийский государственный технический

университет, кафедра радиотехнических и медико-биологических систем, Йошкар-Ола; E-mail: [email protected] Сергей Аркадьевич Охотников — аспирант; Марийский государственный технический университет, ка-

А. В. Ушаков, Е. С. Яицкая

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

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

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

Введение. Постановка задачи. Устройства рекуррентного преобразования двоичных кодов являются двоичными динамическими системами. В теории двоичных динамических систем (ДДС) сформировались два направления [1—3] их синтеза, которые развиваются по законам системной диалектики практически независимо друг от друга.

Первое направление, основанное на теории конечных автоматов, связано с решением задач вычислительной техники при синтезе микропрограммных автоматов, а также частично задач общетехнической дискретной автоматики. С точки зрения общесистемного подхода „автоматная" ветвь ДДС представляет собой класс нелинейных двоичных динамических систем (НДДС), аналитически описываемых выражениями

федра радиотехнических и медико-биологических систем, Йошкар-Ола; E-mail: [email protected]

Рекомендована кафедрой радиотехнических и медико-биологических систем

Поступила в редакцию 13.12.11 г.

УДК [517.938 + 519.713 / .718]: 621.398

(1) (2)