Aula 5 Transformada de Fourier

Fazer download em pdf ou txt
Fazer download em pdf ou txt
Você está na página 1de 36

Curso de Fsica Mdica

Disciplina: Processamento de Imagens Digitais


Prof. Valnir de Paula
[email protected]

Bases Matemticas do
Processamento de Imagens:

Sries de Fourier e
Transformadas de Fourier

Sries de Fourier
Jean Baptiste Joseph Fourier (1768-1830)
viveu na poca de Napoleo. Mas, seu nome
foi imortalizado pelas sries trigonomtricas
que introduziu em 1807.
Fourier foi levado a desenvolver suas
sries ao estudar a propagao de calor
em corpos slidos.

Admitindo que essa propagao deveria ocorrer por ondas


de calor e levando em conta que a forma mais simples de
uma onda uma funo senoidal, Fourier mostrou que
qualquer funo, por mais complicada que seja, pode ser
decomposta como uma soma de senos e cossenos.
2

F1

Sries de Fourier
A figura mostra o grfico da funo sen(x), onde x um
ngulo medido em radianos.
Essa funo PERIDICA, isto , sua forma se repete a cada
PERODO. Neste caso, a funo seno se repete a cada perodo de 2.
O valor mximo da funo, ou a sua AMPLITUDE, 1.

F1

Sries de Fourier
A funo cosseno tambm peridica, com o mesmo
perodo e amplitude que o seno, mas deslocada de /2 em
relao ao seno.
Portanto as funes seno e cosseno diferem na FASE e a diferena
de fase entre elas de /2.

F1

Sries de Fourier

A soma das funes sen(x) e cos(x) obtida traando-se, em cada


ponto x, a soma dos valores de sen(x) e cos(x).

O ponto da curva na
regio x=5,5 zero, pois
o valor de sen(x) igual
e de sinal oposto ao
valor de cos(x) nesse
ponto.

F1

Sries de Fourier

Uma funo peridica pode ser bem mais complicada que


uma senide. No exemplo da funo f(x) mostrada abaixo, a
curva tambm peridica, mas no apenas um seno ou
um cosseno. Como achar uma funo matemtica que
descreva uma curva como essa?
Segundo Fourier, qualquer
funo peridica, por mais
complicada que seja, pode ser
representada como a soma de
vrias funes seno e cosseno
com amplitudes, fases e
perodos escolhidos
convenientemente.
6

F1

Sries de Fourier
Na prtica, no possvel o trabalho com infinitas parcelas e um nmero
limitado deve ser empregado. Na figura, o sinal retangular dado pela funo f(x)
resultante de:

f(x) = 5 + (4)sen x + (4/3)sen 3x + (4/5)sen 5x + (4/7)sen 7x + ...


Na figura, f(x) a soma das 5 parcelas explcitas nessa relao, que j
produzem uma certa aproximao. Se fossem infinitas, o resultado seria uma
forma geomtrica perfeita ( linha tracejada).

F 20

Sries de Fourier

Em resumo, qualquer funo f(x) pode, segundo Fourier, ser


escrita na forma da soma de uma srie de funes seno e
cosseno da seguinte forma:

f(x) = a0+ a1sen(x)+ a2sen(2x)+a3sen(3x)+....


.....+b1cos(x)+b2cos(2x)+b3cos(3x)+...
Os termos tipo seno e cosseno podem se estender
indefinidamente, se necessrio, para melhor representao da
funo original f(x).
Os coeficientes a0, a1, a2... b1, b2, etc. so as amplitudes de cada
onda componente do desenvolvimento em srie.
8

F1

Sries de Fourier

O primeiro termo (a0) constante e pode-se dizer que ele


o componente de corrente contnua do sinal.
O segundo termo ( sen x) tem o mesmo perodo ou mesma
frequncia do sinal original. Por essa igualdade,
denominada oscilao fundamental (ou primeiro harmnico)
do sinal.
Os termos seguintes tm frequncias mltiplas (sen 2x, sen
3x, ...) da fundamental. So chamados harmnicos do sinal.
Portanto, pode-se dizer que todo sinal peridico formado
por um componente contnuo (que pode ser nulo), uma
oscilao fundamental e oscilaes harmnicas.

F 20

Sries de Fourier

Frequncia fundamental e correspondentes harmnicos


10

Sries de Fourier
A serie de Fourier de uma funo f definida no intervalo (-l; l)
e dada por

a0
f ( x) an cos(n0 x) bn sen(n0 x)
2 n1
onde 0 a frequncia (0

2
)
T

11

F 21

Sries de Fourier
Os coeficientes a0, an e bn so dados pelas frmulas de EulerFourier:

2
a0 f ( x)dx
TT
2
an f ( x) cos(n0 x)dx, n 1, 2, ......
TT
2
bn f ( x) sen(n0 x)dx, n 1, 2, ......
TT
12

F 21

Espectro de Fourier de uma Imagem


A imagem espacial de um pulso quadrado e a sua representao no espectro
de frequncias de Fourier. As frequncias so dispostas em um padro de
distribuio.

13

F 13

Espectro de Fourier de uma Imagem


A princpio parece difcil entender a
visualizao da imagem, pois um
ponto de uma imagem representada
no domnio Fourier (ou da frequncia)
pode conter informaes sobre toda a
imagem no domnio espacial,
indicando quanto desta frequncia h
na imagem.
No espectro de Fourier possvel
ocultar algumas frequncias e alterar
a imagem original.

14

F6

Transformada de Fourier

A transformada de Fourier um limite formal da srie de


Fourier de uma funo f peridica no caso em que f .
Um processo fsico pode ser descrito em seu domnio
espacial, tomando valores de uma determinada
quantidade h em funo do espao, neste caso h(x);
Ou no domnio da frequncia, cujo processo especificado
pela amplitude H como funo da frequncia [H(f)].

15

F 4, p 12

Transformada de Fourier

Deve-se relacionar h e H como duas representaes da mesma


funo, utilizando como operador a TF. Neste caso, as equaes
podem ser escritas como:

H ( f ) h( x)e 2ifx dx

Transformada de Fourier

h( x )

2ifx
H
(
f
)
e
df

Transformada Inversa de Fourier

16

F 4, p 12

Transformada de Fourier
A Transformada de Fourier muda um sinal do domnio do
tempo (ou do espao) para o domnio da frequncia.

Mudamos um sinal de um domnio o qual no temos


controle para outro domnio que pode ser facilmente
modificado (a frequncia).

17

F 11, p 158

Transformada de Fourier

Interpretao Fsica

18

F2

Transformada de Fourier de uma imagem


A chamada "onda quadrada", ou "funo degrau", obtida
pela soma dos termos de srie de Fourier e ilustra uma
sucesso de "bits" com valores 1 e 0.

Quanto mais termos somados na srie de Fourier, melhor ser a aproximao


da onda quadrada.
19

F 13

Transformada de Fourier de uma imagem


y

Seja uma imagem bidimensional


qualquer, monocromtica

Podemos representar qualquer


ponto na imagem por uma
intensidade luminosa I(x,y)

x
20

F 13

Transformada de Fourier de uma imagem


A ideia bsica para o processamento neste domnio pode ser
observada no esquema apresentado na figura adiante. Em
um primeiro passo, a imagem I(x,y) transformada
para o domnio de Fourier, atravs da sua transformada
discreta.
A imagem no domnio de Fourier representada por F(u,v)
e esta convoluida com o filtro H(u,v).
Ao produto F (u,v) H(u,v) aplicada a inversa da
transformada de Fourier para retornar ao domnio
espacial, onde se tem a imagem processada.

21

F6

Transformada de Fourier de uma imagem

No caso bidimensional, basta


decompor em duas frequncias, uma
para cada dimenso da imagem

y (1 , 2 )

I ( x, y)e j 1x 2 y dxdy

Neste caso, ao invs de fazer um


grfico unidimensional, a
transformada de Fourier corresponde
a um grfico bidimensional cujo valor
no 3 eixo corresponde a y.
22

F 13

Transformada inversa de Fourier

Se y(w1,w2) conhecida,
podemos recuperar a informao
de intensidade espacial.

Isto chamado transformada


inversa de Fourier e nada mais
que a transformada da
transformada de Fourier (mas note
o sinal trocado na exponencial).

23

F 13

Algumas transformadas de Fourier


Imagem radiolgica

Domnio espacial

Domnio da frequncia
24

F 13

Algumas transformadas de Fourier


Imagem de TC

Domnio espacial

Domnio da frequncia
25

F 13

Algumas transformadas de Fourier


Imagem de RM

Domnio espacial

Domnio da frequncia
26

F 13

Algumas transformadas de Fourier


Imagem Mamogrfica

Domnio espacial
F 13

Domnio da frequncia

27

O uso de transformadas de Fourier como


mtodo de edio de imagens
Em algumas circunstncias, o uso da TF pode ser
bastante til na edio de imagens:
Remoo de rudos e artefatos
Quando estes possuem frequncia muito bem definida, sendo
bem localizada na TF.

Remoo de padres
Por exemplo, uma cerca pode ter um padro de frequncias
bem definidas.

Filtros de efeitos especiais


A remoo de algumas frequncias pode criar efeitos
interessantes
28

F 13

O uso de transformadas de Fourier como


mtodo de edio de imagens
Os picos no domnio das frequncias (estrelas, linhas ou pontos fora
de padro peridico) indicam artefatos ou ruido da imagem.
Eles podem ser ocultados do
espectro para eliminar estes
defeitos.

29

Alguns exemplos:

Remoo de sombras
Neste caso, a sombra possui
estrutura muito caracterstica em
frequncia

30

F 13

Alguns exemplos

Outros mtodos
Por exemplo, remoo de uma
estrutura espria

31

Alguns exemplos
Imagem original

Transformada de Fourier

O padro de linhas da imagem original representado por frequncias fora do


centro do espectro de Fourier.
32

Alguns exemplos
Supresso de frequncias

Imagem transformada

Retirada das frequncias do espectro e imagem corrigida atravs da inversa da


Transformada inversa de Fourier.
33

Alguns exemplos
(a) Imagem obtida por adio de onda senoidal e o seu espectro de
Fourier equivalente ( b ) .

34

Alguns exemplos
A interferncia indesejada causada pelo padro de brilho sinusoidal pode ser
removida editando-se a informao da frequncia espacial, como mostrado pelas
reas escurecidas em (c). A TF inversa, recupera a imagem original (d).

35

Referncias
BOYCE, William E.; DIPRIMA, Richard C. Equaes Diferenciais
Elementares e Problemas de Valores de Contorno. 3 Edio. Guanabara
Koogan. 1990.

GONZALEZ, R. C., WOODS, R. E. Processamento de Imagens Digitais,


Editora Edgard Blucher Ltda, 2000.
http://dfn.if.usp.br/~suaide/ November 12, 2007.
As Sries de Fourier. Disponvel em
http://www.seara.ufc.br/tintim/matematica/fourier/fourier2.htm
GONALVES, Louis Augusto. Um estudo sobre a Transformada Rpida de
Fourier e seu uso em processamento de imagens. Dissertao de
mestrado. Porto Alegre. 2007

36

Você também pode gostar