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

Методичка_Сигналы

.pdf
Скачиваний:
22
Добавлен:
02.05.2015
Размер:
690.02 Кб
Скачать

71

Обозначив GT Ck и повторив процедуру предельного перехода при T для

формулы (2.143), можно получить каноническое разложение стационарной случайной функции X t :

 

 

X t G ei t d ,

(2.152)

где функция Sx есть дисперсия случайной функции G d .

Отметим основные свойства функции спектральной плотности. Для этого воспользуемся формулой Эйлера и перейдем от формулы (2.150) к одностороннему спектру в тригонометрической форме представления:

 

1

 

 

1

 

 

i

 

 

Sx

Kx e i

 

Kx cos d

Kx sin d .

(2.153)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как функция

Kx четная,

то второе слагаемое равно нулю, а

первое

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

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Sx

Kx cos d ,

 

(2.154)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда видно, что функция Sx является действительной и четной функцией, то есть

Sx Sx .

Поэтому в выражении (2.149) тоже можно ограничиться только положительными частотами:

 

 

Kx Sx ei d .

(2.155)

Выражения (2.149) и (2.150), а также (2.154) и (2.155) являются парами интегрального преобразования Фурье (прямые и обратные). Исходя из свойств преобразования Фурье корреляционная функция Kx и спектральная плотность

Sx подчиняются закономерности: чем уже одна из них, тем протяженнее вторая, и наоборот.

Интересно отметить, что площадь, ограниченная непрерывной кривой Sx на спектральной диаграмме, равняется дисперсии Dx случайной функции X t .

Действительно, положив в формуле (2.155) 0 , получим:

72

 

 

Kx 0 Dx Sx d .

(2.156)

0

 

Если под случайной функцией X t подразумевать напряжение, то Dx

можно

рассматривать как среднюю мощность, выделяемую этим напряжением на активном сопротивлении 1 Ом:

Dx М X 2 t Px .

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

dPx Sx d

представляет собой долю средней мощности, выделяемой составляющими спектра,

относящимися к

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

, d . Поэтому

спектральную плотность

Sx называют

еще

спектральной плотностью мощности или энергетическим

спектром стационарной

случайной

функции, поскольку

Sx имеет размерность

энергии.

 

 

 

 

§2.18. Случайные поля при исследовании природных образований

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

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

 

 

 

 

Sx 1, 2

Kx 1 , 2

e i 1 1 2 2 d 1d 2 .

(2.157)

 

 

 

 

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

плотности

Sx 1, 2 найти корреляционную функцию Kx 1 , 2 .

73

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

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

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

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

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

Исследование профилей изображений лесных массивов показывает, что соответствующие спектральные плотности обычно имеют максимум при частотах около 0,1 м-1 .

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

74

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

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

§2.19. Математическое описание непрерывных двумерных сигналов на примере изображений

Распределение энергии источника светового излучения по пространственным координатам x, y времени t и длинам волн можно описать функцией

C x, y, t, .

(2.158)

Энергия излучения пропорциональна

квадрату амплитуды и

поэтому

представляет собой действительную положительную величину.

 

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

в ограниченной области

 

x

 

Lx ,

 

y

 

Ly и

наблюдаются в течение

конечного

 

 

 

 

промежутка времени T t T .

 

 

 

 

 

 

75

Таким образом, функция C x, y, t, является ограниченной функцией четырех ограниченных переменных.

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

 

 

 

Y x, y, t

C x, y, t, Vs d ,

(2.159)

 

0

 

где Vs - спектральная чувствительность человеческого зрения.

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

 

 

 

Fi x, y, t

C x, y, t, Si d ,

(2.160)

 

0

 

где Si - спектральная чувствительность i-того датчика.

Функция F x, y, t представляет собой распределение яркости или какой-либо другой физической величины, связанной с яркостью. Очень часто, особенно при работе с изображениями, аргумент t меняется дискретно, поэтому в дальнейшем аргумент t можно опустить.

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

 

Q a1F1 x, y a2 F2 x, y a1Q F1 x, y a2Q F2 x, y ,

(2.161)

где a1 , a2 - некоторые постоянные;

 

 

Q - оператор, описывающий преобразование входного сигнала;

 

F1 x, y ,

F2 x, y - двумерные функции.

 

 

Используя

свойство дельта-функции, сигнал на

входе системы F x, y

можно

представить как взвешенную сумму дельта-функций

 

 

 

 

 

 

 

F x, y F , x , y d d ,

(2.162)

 

 

 

 

где F ,

- весовой множитель дельта-функции,

имеющего координаты

, на

плоскости x, y .

Если на входе линейной системы действует сигнал F x, y , то на выходе имеет место сигнал

76

 

 

G x, y Q F x, y

(2.163)

или

 

 

 

 

 

 

 

 

G x, y Q

F , x , y d d F , Q x , y d d .

(2.164)

 

 

 

 

Можно поменять местами операции линейного преобразования и интегрирования, так как линейный оператор действует только на тот множитель подинтегрального выражения, который зависит от пространственных переменных x, y .

Запишем второй множитель подинтегрального выражения в виде

H x, y, , Q x , y

(2.165)

и будем называть эту функцию импульсным откликом двумерной системы. Подстановка H x, y, , в соотношение (2.164) дает интеграл суперпозиции

 

 

 

G x, y

F , H x, y, , d d .

(2.166)

Линейная двумерная система называется пространственно-инвариантной, если ее

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

y .

Для такой

системы

 

 

 

 

H x, y, , H x , y

 

(2.167)

и интеграл суперпозиции имеет особую форму, называемую интегралом свертки

 

 

 

 

 

G x, y F , H x , y d d .

 

(2.168)

 

 

 

 

Операция свертки символически записывается так

 

 

 

G x, y F x, y H x, y .

 

(2.169)

 

Свойства интеграла свертки.

 

 

1. Интеграл свертки симметричен, то есть

 

 

 

 

 

 

G x, y F

, H x , y d d F x , y H , d d . (2.170)

 

 

 

 

2. Спектр функции, полученной в результате свертки двух функций

F x, y и

H x, y равен произведению спектров исходных функций, то есть

 

 

 

QF F x, y H x, y F x , y H x , y ,

 

(2.171)

где F x , y QF F x, y - спектр функции F x, y ;

 

 

H x , y QF H x, y - спектр функции H x, y .

 

 

77

3. Спектр функции, полученной в результате перемножения двух функций F x, y и H x, y , пропорционален свертке спектров исходных функций, то есть

Q

 

F x, y H x, y

1

F

,

 

H

,

 

,

(2.172)

 

4 2

 

 

 

F

 

x

 

y

x

 

y

 

 

где F x , y QF F x, y - спектр функции; H x , y QF H x, y - спектр функции.

Наряду с одномерными сигналами, которые рассмотрены ранее, при моделировании различных полей, например, изображениями, приходится иметь дело с двумерными сигналами. К таким сигналам также может быть применено преобразование Фурье, которое в этом случае носит название двумерного. В результате двумерного преобразования Фурье сигнал F x, y , описывающий, скажем, изображение, преобразуется в двумерный спектр (двумерную спектральную плотность) G x , y этого изображения, который определяется из соотношения

 

 

 

 

 

 

 

 

 

G x , y F x, y e i x x y y dxdy ,

(2.173)

 

 

 

 

 

 

 

где x , y - пространственные частоты вдоль оси x и y.

 

В общем

случае спектр G x , y есть комплексная величина.

Ее можно

разложить на действительную и мнимую части:

 

 

 

G x , y Re x , y i Im x , y ,

(2.174)

и представить с помощью амплитуды и фазы

 

 

 

 

G x , y A x , y ei x , y ,

(2.175)

где A x , y

Re2 x , y Im2 x , y ;

 

x , y arctg

Im x , y

.

 

 

 

 

(2.176)

 

 

 

Re x , y

 

Достаточным условием существования Фурье-спектра функции F x, y является

абсолютная интегрируемость этой функции:

 

 

 

 

 

 

 

 

 

 

 

F x, y

 

dxdy .

(2.177)

 

 

 

 

 

 

 

 

 

 

 

По известному Фурье-спектру может быть восстановлено исходное изображение

F x, y :

78

1

 

G x , y e

i x x y y

 

F x, y

 

 

d x d y .

(2.178)

4 2

Это преобразование называется обратным двумерным преобразованием Фурье.

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

 

1. Если функция F x, y

разделима по пространственным переменным, так что

F x, y fx x f y y , то

 

 

 

 

G x , y gx x g y y ,

(2.179)

где gx x и g y y - одномерные Фурье-спектры функций fx x и f y y .

 

2. Оператор двумерного преобразования Фурье линеен

 

Q a F1 x, y b F2 x, y a G1 x , y b G2 x , y .

(2.180)

3. Изменение масштаба пространственных переменных приводит к обратному изменению масштаба пространственных частот и пропорциональному изменению значений спектра

 

F ax,by

1

 

 

 

 

 

 

 

 

 

 

Q

 

 

G

x

a

,

 

y

.

(2.181)

a b

 

 

F

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4. Сдвиг (изменение координат) на исходной плоскости приводит к фазовым

изменениям на частотной плоскости

 

 

 

 

 

 

 

 

 

 

 

 

Q F x a, y b G

x

,

y

e i a x b y .

(2.182)

F

 

 

 

 

 

 

 

 

 

 

 

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

Q 1 G

x

a

,

y

b

F x, y ei ax x by y .

(2.183)

F

x

 

y

 

 

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

Пусть функция FH x, y описывает непрерывное изображение, представляющее собой распределение яркости, оптической плотности или какого-либо другого

79

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

 

 

S x, y x l x, y k y ,

(2.184)

l k

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

x, y (рис.2.17).

y

x

Рис. 2.17

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

Fg x, y FH x, y S x, y FH l x, k y x l x, y k y , (2.185)

lk

вкотором учитывается, что функцию FH x, y можно внести под знак суммирования и

задать ее значения только в точках отсчета.

Для анализа процесса дискретизации удобно использовать спектр Fg x , y ,

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

Fg x , y

 

 

Fg x, y e j x x, y y dxdy .

(2.186)

80

В силу свойства спектра свертки, спектр дискретизированного изображения

можно представить в виде свертки

спектра

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

FH x , y и

спектра дискретизирующей функции

J x , y , то есть

 

 

 

 

 

F

 

 

,

 

 

 

1

F

 

 

,

 

J

 

,

 

.

(2.187)

 

 

4 2

 

 

 

 

g

 

x

 

y

 

H

 

x

 

y

 

x

 

y

 

 

Двумерное преобразование Фурье дискретизирующей функции дает бесконечный набор дельта-функций в плоскости пространственных частот с шагомxs 2 x и ys 2 y . Таким образом,

J x , y

4

2

x l xs , y

k ys .

(2.188)

 

 

 

 

 

 

x y

l k

 

 

Так как спектр реального

изображения ограничен, то

есть FH x , y 0 при

 

x

 

x.max

 

 

Fg x , y

и

y

 

y.max , то вычисляя с помощью (2.187) свертку, найдем

 

 

1

 

 

 

 

 

FH x , y l xs , k ys d d . (2.189)

x

y

 

l k

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

Fg x , y 1 FH x l xs , y k ys . (2.190)

x y l k

На рис.2.18а изображен типичный спектр исходного изображения, а на рис.2.18б - спектр дискретизированного исходного изображения.