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

книги / Методы и средства цифровой обработки пространственно-временных сигналов

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
5.9 Mб
Скачать

Министерртво высшего и- оредпего специального ойравования Й&ВР Уральский ордена Трудового Краоного Знамени Нолитехнический институт им.О*М.Кирова

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

Межвузовский' оборник научая* трудов

Свердловск 198В

Включает отатьи, поовященные развитию и использованию ма­ тематического аппарата для решения вадач обработки простран­ ственно-временных сигналов.

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

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

Помет быть рекомендован студентам радиотехнических, электротех­ нических и приборостроительных специальностей.

Рецензент: кафедра вычислительной техники Уральского лесо­

технического института

 

Утвержден редакционно-изда эльским ооветом института (вып.1)

Р е д а к ц и о н н а я

к о л л е г и я :

проф., д -р техн.наук Ю.И.Алимов (ДОМ УНЦ АН С00Р)з чл .-кор. АН

Кав.СОР В.М.АмербаеВ (МИЭТ); проф.,

д -р техн.наук И.Я.Билин-

ский

(ИЭВТ Латв.ССР)* от .н а у ч .оотр .,

канд. техн.наук И .Б.Гуре-

вич

(ИПК АН СССР)з чл.-кор. АН СССР Ю.И.Куравлев (ВЦ АН СССР);

акад. АН УССР А.Л.Кухтенко (МК АН УССР)* доц. , канд.техн.наук

И.И.Канатов (Л?ТИ)| от.неуч. о о т р .,

канд.техн.наук Г .0 .Колмого­

ров

(УЛИ): Д оц.,

вавд.техн.наук В.Г'.Лвбунец - отв.редактор

(УЛИ); мл.науЧ.

ротр. В.Н.Круглов -

отв.секретарь (УШЮ* доц .,

канд.техн.наук В.А.Лексаченко

Ш АЛ ): доц .,

каНД.техН.наук

А.П.Мальцев (УЛИ); проф., д -р

техн.наук П.МЛеголин (НТК АН

БССР); йроф ., д -р физ.-мат.наук А.Г.Ченцов

(ИДО УНЦ АН СССР)

© Уральский политехнический институт

им. С. М. Кирова,

1908

И.И.Канатов, В.Л.Литвинов, В.П.Ситов (Ленинградский электро­ технический институт)

АЛГОРИТМ ВЫДЕЛЕНИЯ ЛОКАЛЬНЫХ ОБЪЕКТОВ НА ОДНОРОДНОМ ФОНЕ

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

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

ление его координат на изображении. Для этой

цели предлагается

использовать харакл 'тише

особенности спектров

таких объектов

в базисе Уолша. С од. *>й

стороны, применение спектральных мето­

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

антность ПУ к

 

сдвигу,

например, октавная фильтраций!

 

 

 

 

 

 

У ” '*/(К+1№ т*Ч ...

 

 

 

 

 

 

в т - 1

[ А . . Щ * ч ]

;

 

 

 

где

6(т). ,

 

А'=о

7 =н-2*~п'

 

 

Ш

\

ас -

-

.

.

с Номером

 

 

 

спектр в октавной полосе

 

 

значение

сигнала на

С -м отсчете;^

М?77* -/ (^)

-

 

С -е

значение

функции Уолша с номером 2 т-1

г ;

' /

КОЛИ49СТ—

во

октав.

 

 

 

 

 

 

 

 

 

Чем медное ШДОЛЯ6АЩЙ объект, тем меньше его энергетичес­ кий вклад в. октавную составляющую спектра. Ооновой модульных преобравоввдай является одрегэление на кеадой отдельной итерации абсолютного эначания спектра Уолша, вычисляемого при помо­

щи быотрнх алгоритмов. <3 ном щью этих методов

невозможно опре­

делить местоположение объектов. Военикает явное

противоречие

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

В основе предлагаемого алгоритма лежит возможность выявления локальных объектов по коррелироввнности .их двумерного спектра в баэиое Уолша и маосива1 преобразованного по Уолшу относитель­ но первой строки и первого столбца вышеупомянутого спектра*

Покажем. Что для любого локального объекта спектр будет

разряды двоичного кода номера оточета. иоратное лреоорааование Уолша

 

« >

Используя гомсморфием переменных 5

я И , поменяем их фи

эический сМыол, т .е . будем считать,

Что .исходный сигнал задан

спектров Уол-а, составляющие которого имеют одинаковые нонуле*

вне значения на

6

ооседних позициях. Это позволяет

наг­

лядно представить

вид преобравованного по <2 ) сигнала как

ре­

зультат суммирования соседних по номеру функций Уолша» Чтобы

упростить дальнейшие выкладки,

обозначим

черев

Р:

при всех

/

* (рТГП .

1Ъопмотр1ш сумму дцух функций

Уолша о

номерами

I

и 1+ 1

при четном

I :

так как

I

к ( С + 1 )-я

кодовые комбинации двоичного

ко­

да ра влнчаютоя только нулевш

раэшдом. Для кода Грея в инверс­

ном порядка

Рл * О для вое:. 2 в

 

и

Р0 = I

-

2 =

^ ^

7

.

Тогда <4> примет следующий вид:

 

 

 

 

 

 

 

 

, 8=0,ап~’-1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

<5

)

 

Иополъвуя

это

свойство, можно показать, что ^ля

суш и $у

,

 

У *Й 7 ^ /

 

функций Уолша о

номерами от ^

^

до

2

 

>

( -

0,2п~у-1

л.^врно олвдг'ющев:

 

 

 

 

 

 

 

 

 

 

,

$ ~ 0 у2 ™ ~ -1 ,

 

 

 

 

 

 

д

 

,

- . » - г

- я

 

<в)

 

 

 

 

 

 

8 - 2 п'г,2 'г~1

 

 

 

Из ( 6) видно|

что

такая суш а

дает спектр, представляющий оо -

бой функцию Уолша размврнооти

 

и располагающийся в

низ­

кочастотной части спектральной области* Необходимо также заме­

тить,

что амплитуда спектральных поставляющих будет равна

2 ^ *

а не

I ,

Что позволяет проводить

идентификацию объекта. №

прак­

тике же

Объекты* спектр которых будет определяться выражением

( 6) ,

встречаются редко. .Как правило»

спектр малопротяженных

объектов

состоит из нвокоЛькйх таких сумм* Дня объекта, занима­

ющей

 

С Позиций»

начиная о

р

+•й и иопбльауя выражение

( 6 ) 4 можно показать*

что

 

 

 

Р < 2 п-’ -Ч < р Н ~ 1 ;

(7 а )

р < 2 л ч ~ '( Ш ) ч < р + с - и

 

{2п-к'Ч Н [ 2 У( , й га + 1 )Ч ];

С7 « )

(2 п~л'*а ч Н ) $ [ & , 2 уа ч ) - ! ] .

 

Замена |/ на П~У~1 осуществлена потому, что суммирова­ ние необходимо начинать с объединения наибольшего возможного количества рассматриваемых функций Уолша в соответствии с ( 6)* При этом каждая функция может входить только в одну единствен­ ную сумму <7 а ). Условия (7 б) характерны для локальных объек­ то в . В этом случаэ основную информацию об объекте несет суш а наибольшего количества 2 V функций Уолша. Она будет "зашум­ лена" искажениями, вносимыми суммами меньшего порядка, которые могут изменить амплитуду спектральных составляющих, но не из­ менят их знак (для ПУ это, по существу, фаза преобразования

[ г ] ) .

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

°^7

Видно, что этот массив сохранит знаковую (фазовую) структуру опектра в соответствии о первым столбцом и первой строкой. Пусть ( $1 , йд ) . - взаимно корреляционная функция сформированного массива и исходного спектра. Тогда, если для низкочастотной части спектра

т ^ Ш У > т ,

где

/

- постоянная,

зависящая от серого фона изображения

оС0 д

, то можно говорить о существовании локального объекта.

 

Из формул (6 ), (7)

можно вывести правило определения мес­

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

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

размером

2 т* 2 т .

Тогда координаты объекта

можно указать о

точностью

до

квадрата

величиной 2 п~т *2 п~т

,

где 2п*2п-

размер всего

изображения. При этом но количеству

указанных

отрезков можно определить координаты правой нижней точки квадра­

та. Ширина зоны нечувствительности

выбирается порядка 2аС0012т=

= оС00 12т’ 1

На второй итерации можно увеличить анализи­

руемую область*

уменьшив соответственно ширину зоны нечувстви­

тельности.

"Квадрат'1 второй итерации должен содержать т в

"квадрате" первой итерация. Однако из-за присутствия шумов на

некоторой

/

-Й итерации

/

-й квадрат выйдет за пре­

делы

~ 1)

- г о , что явл' этоя

критерием окончания проце­

дуры обнаружения. Таким образом,

чем меньше

отношение сигнад/шум,

тем точнее можно определить местоположение

объекта.

Заметим, что при расчете на ЭВМ предложенного алгоритма необходимо хранение в памяти лишь первой строки и первого стол­ бца спектра. Формирование массива по правилу Уолша и вычисление значения ВКФ между !шм и спектром изображения могут быть выпол­ нены построчно-столбцевым методом и не требуют выделения допол­ нительной памяти. Таким обризом, предложенный алгоритм при по­ мощи малых вычислительных затрат позволяет определить местопо­ ложение и размеры мвлопротлженных объектов.

1.ЯРОСЛАВСКИЙ Л.П., Введение в цифровую обработку изобра­ жений. М. , 1979. 312 о.

2.А.В.ОПИЕНХАЙМ, ДНЕ.С. ЛИМ. Важность фазы при обработке сигналов / ТИИЭР. 1961. Т.69, # 6 . С.5 - 17,

УДК 681.327

В.ГЛабунец (Уральокий политехничес­

 

кий институт)

СВЕРХБЫСТРОЕ МНОГОМЕРНОЕ ПРЕОБРАЗОВАНИЕ РАДСЧА-ХАРТЛИ. ПРИМЕНЕНИЕ В СИСТЕМАХ СЕЛЕКЦИИ ДНИЖУ1ДОСЯ ОБЪЕКТОВ

За последние годы получены многочисленные быстрые алгорит­ мы (БА) различных одномерных ортогональных преобразований (ОП). Мнотомерным ОП уделялось мало вк..мвнп, так как считалось, что

- 7 -

одномерный БЛ автоматически давт "многомерный алгоритм", приме­ нением одномерного по всем независимым координатным направлени­ ям (быстрая "строчно-столбцовая" обработка многомерных сигналов). Ив исследований [ 1 - 4 ] , видно, что в многомерном случае* имекггся дополнительные возможности, позволяющие значительно снизите

общие вычислительные

затраты* ЕА. Это

было* доказано

на примере

нового^ БА

\1 вер н ого дискретного преобразования Фурье (Д О ),

ари^метячеокая оложнооть которого

в

^

рае меньше

оложности

классических

"строчноотолбцовых" БА.

 

 

 

 

 

 

 

 

Покажем,

что аналогичные результаты верны1и для любых ОД

с ядром вида

Н(<(1)1Х>)

,

где

<и)/Х> в^Щ Х^

-

окалйр-

ное произведение в

$ -мерном евклидовом пространстве*

Ху

Пусть

 

Я) - ироо^айотво квадратично интегрируемых,

функций

/ (X ) ' /?у

— Я

♦ Введем вг

 

С г,

 

Я)

пару

Н -

преобразований:

 

 

 

 

 

 

 

 

 

 

 

 

№ ) ' ! № № № > №

 

 

 

 

,

 

 

 

т

 

/(X) * $ т )Щ < ю / х >)Ш й'л [Р (Ю )I ,

 

 

 

где

Я*

-

оопрчженное *к

 

/?у

пространство,

 

-рав­

но по обозначению. Функцию

 

Г(IV)

назовем многомерный

И -

спектром оигнала

/ (X ) .

 

 

 

 

 

 

 

 

 

 

 

Заметим,

что

К(<Ю!Х>)

 

представляет

собой

 

плоскую

я

-волну,

распространяющуюся в

Яу

вдоль

вектора

Ш

.Б ол и

К(<Ю/Х>) = ^Хр(1<Ю1Х>)

то

Н -волна

превращается в рбычную гармоническую волну. Разложение

(2>

 

функции

у (х )

есть

ее представление

в

виде* оуперпоэиции раз­

личных ПЛОСКИХ'

К

-ВОЛН;

 

 

 

/ (X )

 

 

 

 

 

 

Определение

I . Интеграл Функции

 

по .всем

гипер-

* плоскостям < ^

/X > р

 

 

 

 

 

 

 

 

 

 

/ о »

,Р) = /

/(х )# (р -< у х ) с/х йЯ [/(х ) }

 

(э)

называется преобразованием Радона функций

?(Х )

■ где

и *

; , - , ,

 

 

-

единичная сфера

в

Я^

 

 

48010

К!ш ,Р)

называют проекцией функция

{(X) вдоль

гиперплоскости <

> .

л

 

 

Теорема_1 . Одномерный /г -спектр проекции х(^и)Х)

^ -м ер­

ной функции

х (Х )

есть центральное

сечение

V

-мерного

Н-спектра этой же функции.

 

Доказательство. Интеграл

(I)

можно вычислить

в

следующей

последовательности. Так 1г к

<М/Х > ^ < $Ш/Х>

,

где

 

~

у

 

,

то интеКирова юте можно сна­

чала провести по

гиперплоскости

р к<^ш1Х> ,

а

затем по

переменной

/7

•*сю

 

 

 

 

 

 

 

 

 

 

 

^

- ^ (/>

!ш) -

/ ( / / Г - г у ф

=

 

 

7 (/от) (Г(р - ^

 

Н ( р р ) ^ (

,рщ/>Р)Ф

Следовательно,

 

 

 

 

 

 

 

 

 

 

А

 

 

 

 

 

^ ( / ^

= 1

Ш ш>р) К Ш Ф

,

 

(4)

что и требовалось доказать.

 

 

 

 

Аналогичное

утверждение возможно и для дискретных Н -пре­

образований.

Пусть

= ( П*,

, . , /2 и,) е 7^ *2.#'

••* 2^ ,

 

 

 

 

 

 

скалярн#е^произве­

дение

целочисленных векторов'2

й / 2 , где

 

,/2^2^.,

щ а * т

 

-

произвольные

целые числа,

2^ .

* -группа*

сложений целых чисел по модулю

/К* • Введем обозначения

2*& 7

*2

*

*2

( 7 * г = 7 * *7* *

* 7 * .

 

где

/К - (/4^

УК,

 

• Тогда ш есто

(I) и

(2 ) имеем

Р М = 2 ,/ (п )К Ы т > ) ,

(5)

Пусть

 

[об0]

-

такая минимальная совокупность век­

торов, что

лучи в и ^

{ а ^ }

заполняют вое дискретное

пространство

( 2 * У .

 

 

 

 

Определение 2. Функция

/

(об0 , Р)

, определенная на

множестве

евклидовых

/)

-мерных плоскостей

 

<оС°!Л >

 

оС?Л; * Р

равенством

 

 

М

Л

' , у . ^

- 1 г ,Г О !)1(р -< *°т > ),

т

называется дискретным преобразованием Радона (Д1Р) сигнала $(п) Для дискретного -преобразования (5) верна теорема,

аналогичная первой.

Теорема 2. Дискретное Н -преобразование сигнала У(Л) связано с дискретным преобразованием Радона равенством

 

 

 

Г(аоС°) = 1 р ( А Р ) Н ( №

 

 

(в)

 

Из (8 ) следует, что для вычисления

\)

-мерного дискрет­

ного

К -преобразования (5) необходимо одно ДПР и такое ко­

личество

одномерных

И -преобразований,

какова

мощность

множества

К

} , а

она, как оказывается,

в существенной

мере зависит от канонических разложений чисел

} / 1^

 

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

 

 

9) ,

 

 

^

 

Теорема 3. Общее число лучей

 

накрывающих

О Д .

а следовательно, и общее число одномерных

И -преобразова­

ний,

необходимых для вычисления

у

-мерного

К -преобра­

зования,

равно

 

 

 

 

 

 

 

М 9, у) =(^0/(^~0 )

воли

/\1-с

9 VI

УТ;

 

 

 

/77^ - /;

 

 

 

 

 

 

 

Ш

9т) Яд*-Г_/ 1Г ~ Г '

еоли

 

 

V I я!,').

Следствие

I . При

9 У

 

, (I = ГУ)

 

имеем

 

 

 

-о т<пт* птгЛ [ Ч У П'1 У ~ ---------

/

Соседние файлы в папке книги