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

Переработанные лекции (1)

.pdf
Скачиваний:
10
Добавлен:
31.03.2015
Размер:
1.26 Mб
Скачать

Введение

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

Области применения цифровой обработки информации:

Космос — передача цифровых изображений со спутника, например фото Луны, Марса, Венеры, предсказание погоды, «виды» на будущий урожай, пожары, наводнения.

Связь — видеотелефонная связь по каналам, передача графической информации.

Геология — исследование из космоса, авиации, например автоанализ характера местности и исследование природных ресурсов.

Биология — микроскопия, например улучшение качества биологических изображений.

Медицина — рентгенография, томография (рентгеновская, ультразвуковая, ЯМР).

Военное дело — обнаружение целей, распознание противника, например в танковом бою.

Роботизация — техническое зрение.

Криминалистика — идентификация личности.

Области применения цифровой обработки информации в НК:

Оптический метод

Рентгеновский метод

Ультразвуковой метод

Термография

Обязательно требует обработки сигнала:

Томография

Магнитный метод

Вихретоковый метод

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

Основные задачи цифровой обработки информации:

представление изображение (модель восприятия)

улучшение визуального качества

восстановление изображения

анализ изображения

реконструкция изображения

сжатие данных

Восприятие — контраст, пространственные частоты, цвет.

Физические величины — освещенность (фото), затухание (томография), отражение (эхолокация), температурный профиль (термография).

Локальные модели — дискретизация, квантование, детерминированные модели, преобразования — разложение в ряд, статистические модели.

Глобальные модели — анализ сцен, кластеризация, сегментация, распознавание. Улучшение визуального качества — преобразование изображения без увеличения

содержащейся в нем информации, задача — выделение определенных деталей для дальнейшего анализа.

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

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

усиление границ

псевдоцвета

уменьшение шума

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

гистограммные эффекты

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

Распознавание образов — обнаружение образов (или другой информации) и их выделение из сигналов.

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

Специальные системы

Одномерные непрерывные функции (сигналы) от координат или времени обозначаются так: f(x), U*(x), s(t).

Соответствующие им дискретные (цифровые) функции (сигналы): f(n), U*(n),

s(k).

Непрерывное изображение представляет собой функцию двух независимых переменных:

f(x,y), U*(x,y), s(x,t).

Соответствующие им дискретные (цифровые) изображения (двумерные сигналы): f(m,n), U*(m,n), s(m,k).

Графическое представление двумерной последовательности:

Таблица некоторых функций:

Функция

 

Математическое описание

Изображение

Дельта-функция

x 0,

 

x 0

 

 

 

 

 

Дирака

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

x dx 1

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Фильтрующее

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f x

 

свойство

 

f x

x x dx

 

 

 

 

x

 

 

 

 

 

Масштабируемость

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дельта-функция

 

 

0, n 0

 

 

 

 

 

Кронекера

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1, n 0

 

 

 

 

 

Фильтрующее

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

свойство

f m n m f n

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Прямоугольная

 

 

 

 

 

1,

 

 

x

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

функция

rect x

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

0,

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Функция знака

1, x 0

 

 

 

sgn x 0, x 0

 

 

 

1, x 0

Sinc-функция

sinc x

sin x

 

 

x

 

 

Comb-функция

 

(ступенька)

comb x x n

 

n 0

Треугольная функция

1

 

x

 

,

 

x

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tri x

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

0,

 

x

 

 

 

 

 

 

 

 

Двумерные версии функций формируются по правилу: f x, y f1 x f2 y

Двумерная дельта-функция будет иметь вид: Дельта-функция Дирака: x, y x y

Дельта-функция Кронекера: m, n m n

Которые будут удовлетворять свойствам: Для дельта-функции Дирака:

 

 

 

 

 

x

 

 

 

 

 

f x, y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f x , y

x , y

y

dx dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x, y dxdy 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

Для дельта-функции Кронекера:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x m, n

x m , n

m m , n n

 

 

 

 

m n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m, n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m n

 

 

 

 

 

 

 

 

 

 

 

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

1, x y 0x, y

0, x 0, y 0

Большим кружком обозначен отсчет со значением 1, маленькими кружками — отсчеты со значением 0.

Двумерный линейный импульс — последовательность, имеющая постоянное значение в одном направлении и являющаяся импульсной в другом:

 

 

 

1, y 0

1

x, y

 

0, y 0

 

 

 

1, x 0

1

x, y

 

0, x 0

Двумерная ступенька:

1, x 0, y 0 u x, y

0, x 0, y 0

Из свойств дельта-функции Кронекера следует, что любой сигнал можно разложить на сумму взвешенных и сдвинутых единичных импульсов:

xm, n x m , n m m , n n

m n

Линейные системы

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

Примем, что x(m,n) — входной сигнал, y(m,n) — выходной сигнал, тогда: x m, n y m, n

Система является линейной только в том случае, если для любых двух входных сигналов x1(m,n) и x2(m,n), и соответствующих им выходных сигналов линейной системы y1(m,n) и y2(m,n), выполняется выражение:

H x

m, n x

m, n

H x

m, n

 

H x

m, n

1 1

2 2

 

1 1

 

2

2

 

1 y1 m, n 2 y2 m, n

Это называется линейной суперпозицией.

Если на вход линейной системы подается двумерная дельта-функция Кронекера, то реакция системы на это возбуждение описывается передаточной функцией в виде

функции импульсного отклика:

 

 

 

 

 

 

 

m m , n n

h m, n, m , n

 

 

 

 

 

 

 

Импульсная функция является исчерпывающей характеристикой линейной системы. Зная импульсную функцию можно узнать реакцию системы на произвольный сигнал:

 

 

 

 

 

 

 

 

 

 

 

 

 

y m, n

 

 

x m , n m m , n n

 

x m, n

m n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

m , n

 

m m , n n

 

 

 

 

 

m n

 

 

 

 

 

 

 

 

 

 

 

y m, n x m , n h m, n; m , n

 

 

 

 

 

 

 

m

n

 

 

 

 

 

 

 

 

Импульсный отклик системы иначе называется функцией рассеяния точки. Бывают КОНЕЧНЫЕ импульсные характеристики (по аналогии с КИХ-фильтрами) и

— БЕСКОНЕЧНЫЕ (БИХ), в зависимости от того, конечна или бесконечна область

реакции системы на импульсный всплеск.

Система является пространственно инвариантной, или инвариантной к сдвигу,

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

m, n h m, n, 0, 0

Для системы инвариантной к сдвигу должно выполняться выражение:

 

 

 

 

 

 

 

m m , n n

h m m , n n , 0, 0

 

 

 

 

 

 

 

Можно записать, что:

h m, n, m , n h m m , n n

Т. е. импульсная функция является функцией только двух переменных. Это значит, что форма импульсной функции не меняется, если сигнал смещается по n, m осям.

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

вид:

y m, n x m , n h m m , n n

m n

Эта формула представляет собой свертку двух двумерных функций: входного сигнала с импульсным откликом системы.

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

Графическая интерпретация данной процедуры представлена на рисунке:

Массив значений импульсной характеристики был повернут на 180°, сдвинут на (n,m) и наложен на массив x(n′,m′). Сумма значений массивов {x(·,·)} и {h(·,·)} в перекрывающихся областях дает результат.

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

Непрерывный сигнал:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

f

 

 

 

 

x , y

 

g

 

x, y

 

 

x, y

 

 

x, y

 

 

f

 

 

h

 

x x , y y dx dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дискретный сигнал:

y m, n h m, n x m, n x m , n h m m , n n

m n

Преобразование непрерывного сигнала Фурье

Прямое Фурье-преобразование непрерывного одномерного сигнала:

F f x exp j2 x dx

Обратное Фурье-преобразование непрерывного одномерного сигнала:

f x F exp 2 x d

Прямое Фурье-преобразование непрерывного двумерного изображения:

 

 

f x, y exp j2 1x 2 y dxdy

F 1, 2

 

Обратное Фурье-преобразование непрерывного двумерного изображения:

 

 

1, 2 exp j2 1x 2 y d 1d 2

 

 

 

 

 

 

 

 

 

 

f x, y

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Исходная функция f x, y

 

 

 

 

 

Преобразование Фурье F 1, 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x, y

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x0 , y y0

 

 

 

 

 

 

exp j2 x0 1 exp j2 y0 2

 

 

exp j2 x 1 exp j2 y 2

 

 

 

1

1, 2

2

 

 

 

 

 

 

exp x2

y

2

 

 

 

 

 

 

exp

12 22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rect x, y

 

 

 

 

 

 

 

 

sinc 1, 2

 

 

 

 

 

 

 

 

tri x, y

 

 

 

 

 

 

 

 

sinc2 ,

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

comb x, y

 

 

 

 

 

 

 

 

comb 1, 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Свойства непрерывного Фурье-преобразования:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Свойство

 

 

Исходная функция

 

 

 

 

Преобразование Фурье

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отображение

 

f x, y

 

 

 

 

 

 

 

 

 

F 1, 2

 

 

 

 

Линейность

 

 

1 f1 x, y 2 f2 x, y

 

 

 

 

1F1 1, 2 2 F2 1, 2

 

 

Комплексная

 

f * x, y

 

 

 

 

 

 

 

 

 

F 1, 2

 

 

 

 

сопряженность

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сепарабельность

f1 x

f2 y

 

 

 

 

 

 

 

 

 

F1 1 F2 2

 

 

 

 

(разделимость)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F 1 / a, 2 / b

 

 

 

 

Масштабирование

f ax,by

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ab

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Смещение

 

 

f x x0 , y y0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp j2 x0 1 y0 2 F 1, 2

 

Модуляция

 

 

exp

 

j2 x

y

 

 

f

 

x, y

 

 

F

 

,

2

 

 

2

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Свертка

g x, y h x, y f x, y

G 1, 2 H 1, 2 F 1, 2

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

g x, y h x, y f x, y

G 1, 2 H 1, 2 F 1, 2

Пространственная

g x, y h x, y f x, y

G 1, 2 H 1, 2 F 1, 2

корреляция

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I H 1, 2 F 1, 2 d 1d 2

Интегрирование

I

 

 

h

 

x, y

f

x, y dxdy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Некоторые выводы, следующие из свойств:

При применении преобразования Фурье не происходит потеря информации.

По определению преобразование Фурье разделимо, это значит, что преобразование по любой координате можно провести отдельно.

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

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

экспоната

exp j2

 

x

 

y

.

 

 

 

 

 

 

 

 

 

 

 

 

g x, y

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h x x , y y

exp j2 x

 

y

dx dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y y y ):

Если применить свойство коммутативности свертки ( x x x ,

g x, y H ,

2

exp j2

x

2

y

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

Функция

 

H 1, 2

является

Фурье

 

преобразованием функции импульсного

отклика и называется частотный отклик или частотной характеристикой системы. Она показывает комплексную амплитуду ответа системы на определенную пространственную частоту.

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

корреляция между двумя реальными функциями h x, y и

f x, y :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

f

 

 

 

 

 

 

 

c

 

x, y

 

 

x, y

 

 

x, y

 

 

 

h

 

x , y

 

f

 

x x , y y dx dy

Смена переменных показывает, что функция c x, y так же представляет собой свертку h x, y f x, y . В результате получаем формулу:

Ñ1, 2 H 1, 2 F 1, 2

Еще одним важным свойством Фурье преобразования является то, что интегрирование двух функций равно их интегралу их Фурье образов, т.е.:

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

1

2

1

2

 

h

 

 

 

H

 

 

 

 

x, y

 

f

 

x, y dxdy

 

 

,

 

F

 

,

 

d d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если h f , то данная формула превращается в теорему Парсеваля для энергии:

 

 

 

f x, y

 

 

 

 

F 1, 2

 

 

 

 

 

 

2dxdy

 

 

 

2 d 1d 2

 

 

 

 

 

 

 

 

 

 

 

 

 

Т.е. полная энергия функции такая же, как у её преобразования Фурье.

Преобразование дискретного сигнала Фурье

Преобразование Фурье (прямое и обратное) одномерного дискретного сигнала:

X x n exp jn

 

n

 

1

 

x n

X exp jn d

2

 

 

 

 

Преобразование Фурье (прямое и обратное) двумерного дискретного сигнала функции x(m,n) определяется как:

 

 

 

 

 

 

 

X 1,2 x m, n exp j m 1 n 2 ,

1,2

 

m n

 

 

 

 

1

 

 

exp j m 1 n 2

d 1d 2

x m, n

X 1, 2

2

 

4

 

 

 

 

 

 

 

 

 

Здесь функция X 1, 2 является периодичной с периодом каждого аргумента:

X 1 2 ,2 2 X 1 2 ,2 X 1,2 2 X 1,2

Аналогично с преобразованием Фурье непрерывных функций, H 1, 2 ,

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

Свойства дискретного Фурье-преобразования:

Свойство

Исходная функция

Преобразование Фурье

 

 

 

 

 

 

x(m, n), y(m, n), h(m, n),

X ( 1, 2 ),Y ( 1, 2 ), H ( 1, 2 ),

 

Линейность

1x1 m, n 2 x2 m, n

1 X1 1, 2 2 X2 1, 2

 

Комплексная

x* m, n

F 1, 2

 

 

сопряженность

 

 

 

 

 

Сепарабельность

x1 m x2 n

X 1, 2

 

 

(разделимость)

 

 

 

 

 

Смещение

x m m0 , n n0

 

 

 

exp j2 m0 1 n0 2

X 1, 2

Модуляция

exp j

01m 02n x m, n

X 1 01, 2 02

 

 

 

 

 

 

 

 

Свертка

y m, n

h m, n x m, n

Y 1, 2 H 1, 2 X 1, 2