58
Sinal unidimensional 1 X=10 f(x) |F(u)|

Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Embed Size (px)

Citation preview

Page 1: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Sinal unidimensional 1

X=10

f(x) |F(u)|

Page 2: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

espectro centrado em u = 0 espectro centrado em u = 60

Page 3: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Sinal unidimensional 2

f(x) |F(u)|

Page 4: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

espectro centrado em u = 0 espectro centrado em u = 60

Page 5: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Sinal unidimensional 2

f(x) |F(u)|

Page 6: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

espectro centrado em u = 0 espectro centrado em u = 60

Page 7: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

original Fourier

Fourier com shiftFourier com log(1+|F(u)|)

200000

47

Page 8: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtro ideal

F

H

Page 9: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtro ideal ].[1 HF

Page 10: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

F

H

Filtro de Butterworth

Page 11: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtro de Butterworth].[1 HF

Page 12: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtro passa-altas ideal

F

H

Page 13: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtro passa-altas ideal ].[1 HF

Page 14: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtro passa-altas Butterworth

F

H

Page 15: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtro passa-altas Butterworth ].[1 HF

Page 16: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Convolução espacial

0 200 400 m

f(m)

0 200 400 m

0 200 400 m

h(-m)

h(m)

0 200 400

m

h(x-m)

x

0 200 400 x

f(x)*h(x)

800600

Page 17: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Convolução levando-se em conta a periodicidade da DFT

0 200 400 m

h(m)

0 200 400 m

h(-m)

0 200 400 m

h(x-m)

0 200 400 m

f(m)

m

x

100 400

f(x)*h(x)

x0

Page 18: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Convoluçao com as funções estendidas: T = 800 (400+400)

0 200 400 600 800m

f(m)

0 200 400 600 800m

h(m)

0 200 400 600 800m

h(-m)

0 200 400 600 800m

h(x-m)

x

x

f(x)*h(x)

500 800

Page 19: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Etapas da implementação da filtragem por DFT

• Estender o período MxN de f considerando o período PxQ do filtro h Novo período RxS

MxN:f PxQ :h

1)- QN S e 1-PM (R • Calcular a TF de f estendida

• Gerar a função H do filtro de tamanho RxS

• Multiplicar ponto a ponto a transformada de f estendida pelo filtro G = F.H

• Obter a parte real da transformada inversa discreta de Fourier

• Recortar o canto superior esquerdo de g de dimensão MxN g = crop(g)

(G)g 1

(h) H

)(f F RxS

Page 20: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Exemplo: filtragem sem extensão (padding)

f F

H

Page 21: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

].[1 HF

.

.

.

.

.

.

……

sem padding

Page 22: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Exemplo: filtragem com extensão (padding)

f F

H

Page 23: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

.

.

.

.

.

.

Com padding

].[1 HF

Page 24: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtragem gaussiana

original imagem com ruído

Page 25: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

imagem com ruído espectro de Fourier

Page 26: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Máscara h de convolução gaussiana 15x15 (sigma = 3)

Page 27: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

imagem com ruído f

Filtragem no domínio espacial

imagem filtrada: f*h

Page 28: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

espectro de Fourier H da máscara de convolução

Page 29: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

H estendido e com origem no centro

Page 30: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

H estendido e com origem não transladada

Page 31: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

imagem com ruído f

Filtragem no domínio da frequência

).(1 HFg

Page 32: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

g = f*h

Comparação dos resultados

).(1 HFg

Page 33: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

imagem de diferença entre as filtragens

13-

-13

3.41x10- valor mínimo

1.13x10 valor máximo

Page 34: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

H = Espectro máscara de Sobel

1 0 -1

2 0 -2

1 0 -1

h = Sobel

Filtro passa-altas

Page 35: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

H = Sobel

Page 36: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

imagem f )( fF

Page 37: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

g = f * h ).(1 HFg

Page 38: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Contornos após limiarização (10% do maior valor na imagem)

por convolução por TF

Page 39: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Original f

Exemplo: Filtro passa-baixas ideal

Filtros definidos diretamente no domínio da frequência

)( fF

Page 40: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

H ideal

Imagem filtrada

F

).(1 HFg

Page 41: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

H ideal F

Imagem filtrada

).(1 HFg

Page 42: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

H ideal F

Imagem filtrada

).(1 HFg

Page 43: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

H = Butterworth

Imagem filtrada).(1 HFg

Page 44: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtragem gaussiana passa-baixas

1 máximoseu valor do 0.607 atinge filtro o D D Quando

),(

0

2/),( 20

2

DvuDevuH

0.3M D com H 0

M

M

Page 45: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

original f ).(1 HFg

Page 46: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

0.05M D com H 0

20

2 2/),(),( DvuDevuH

M

M

Page 47: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

original f ).(1 HFg

Page 48: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

0.02M D com H 0

M

M

Page 49: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

original f ).(1 HFg

Page 50: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtragem passa-altas

• Definida também por:

baixas-passa H ; altas-passa H

v)(u,H - 1 v)u,(H

pbpa

pbpa

Passa-altas ideal (ideal) 1 pbH

Page 51: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

f

H ideal

).(1 HFg

Page 52: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtro passa-altas gaussiano

)(gaussiano 1 pbH

Page 53: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

f).(1 HFg

H gaussiano

Page 54: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Filtragem com ênfase nas altas frequências

• restitui parte da componente DC perdida na filtragem passa-altas e reforça as componentes de alta frequência do filtro Hpa.

),(),( vubHavuH paEAF

Page 55: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Exemplo:

Original Equalização

Page 56: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Original Butterworth

Page 57: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

Butterworth Equalização

Page 58: Sinal unidimensional 1 X=10 f(x) |F(u)|. espectro centrado em u = 0 espectro centrado em u = 60

0.5 + 2*Hpa Equalização