Wavelet

De GirinoWiki

Wavelet é uma função capaz de decompor e descrever outras funções no domínio da freqüência, de forma a podermos analisar estas funções em diferentes escalas de freqüência e de tempo. A decomposição de uma função com o uso de wavelets é conhecida como transformada de wavelet e tem suas variantes contínua e discreta. Graças a capacidade de decompor as funções tanto no domínio da freqüência quanto no domínio do tempo, as funções wavelet são ferramentas poderosas para a análise de sinais e compressão de dados.

[editar] Características

Para ser considerada uma wavelet, uma função tem de atender as seguintes características:

A área total sob a curva da função é 0, ou seja $\int_{-\infty}^{\infty} \psi (t) dt = 0$
A energia da função é finita, ou seja $\int_{-\infty}^{\infty} |\psi (t)|^2 dt < \infty$

Estas condições são equivalentes a dizer que $ψ(t)$ é quadro integrável ou que pertence ao conjunto $\mathcal{L}^2(\R)$ das funções quadro integráveis. As propriedades acima sugerem que $ψ$ tende a oscilar acima e abaixo do eixo $t$ , e que tem sua energia localizada em uma certa região, já que ela é finita.

Essa característica de energia concentrada em uma região finita é que diferencia a análise usando wavelets da análise de Fourier, já que esta última usa as funções de seno e cosseno que são periódicas. Na análise de Fourier podemos extrair apenas informações sobre o domínio da freqüência, enquanto na análise com wavelets podemos extrair também informações da função no domínio do tempo. Entretanto as wavelets não são tão bem localizadas no domínio da freqüência como as funções da base de Fourier (em uma analogia ao Princípio da incerteza de Heisenberg, chamamos esta relação entre os domínios da freqüência e do tempo de relação de incerteza ou simplesmente de princípio da incerteza). As wavelets provêem uma forma de balancear a incerteza entre o domínio do tempo e o domínio da freqüência.

Para ser utilizada na análise de sinais uma função wavelet precisa também de outra característica que chamamos de condição de admissibilidade, e que permite a existência da transformada inversa de wavelet. Esta característica será discutida mais abaixo.

Imagem:Wavelet - Morlet.png

Uma wavelet de Morlet.

Alguns exemplos de funções que atendem estas características são a função wavelet de Morlet (ver figura ao lado):

$\psi(t)=e^{-t^2}cos \left ( \pi t \sqrt{\frac{2}{\ln{2}}} \right ) \approx e^{-t^2} \cos{(2.885 \pi t)}$

e a curva conhecida como chapéu mexicano (do inglês mexican hat), definida por:

$\psi(t) = (1 - 2t^2)e^{-t^2}$

que é a segunda derivada da função Gaussiana $-0.5e^{-t^2}$

[editar] Transformada de wavelet

A transformada de wavelet decompõe uma função definida no domínio do tempo em outra função, definida no domínio do tempo e no domínio da freqüência. Ela é definida como:

$W(a,b) = \int_{-\infty}^{\infty} f(t) \frac{1}{\sqrt{|a|}}\psi^* \left ( \frac{t-b}{a} \right ) dt.$

que é uma função de dois parâmetros reais, a e b. o símbolo * indica o conjugado complexo. Se definirmos $ψ a, b (t)$ como:

$\psi_{a,b}(t) = \frac{1}{\sqrt{|a|}}\psi \left ( \frac{t-b}{a} \right ),$

Podemos reescrever a transformada como o produto interno das funções $f (t)$ e $ψ a, b (t)$ :

$W(a,b) = \left \langle f(t), \psi_{a,b}(t) \right \rangle = \int_{-\infty}^{\infty} f(t) \psi_{a,b}(t) dt.$

A função $ψ(t)$ que equivale a $ψ 1,0 (t)$ é chamada de wavelet mãe (do inglês mother wavelet) enquanto as outras funções $ψ a, b (t)$ são chamadas de wavelets filhas. O parâmetro b indica que a função $ψ(t)$ foi transladada no eixo t de uma distância equivalente a b, sendo então um parâmetro de translação. Já o parâmetro a causa uma mudança de escala, aumentando (se $a > 1$ ) ou diminuindo (se $a < 1$ ) a wavelet formada pela função. Por isto o parâmetro a é conhecido como parâmetro de escala (do inglês scaling parameter). O termo $\frac{1}{\sqrt{|a|}}$ é um fator de normalização que garante que a energia de $ψ a, b (t)$ seja independente de a e de b, tal que:

$\int_{-\infty}^{\infty} |\psi_{a,b} (t)|^2 dt = \int_{-\infty}^{\infty} |\psi (t)|^2 dt$

[editar] Transformada inversa de wavelet

Como usamos wavelets para transformar uma função, precisamos também da transformada inversa, de forma a recompor o sinal no domínio do tempo a partir da sua decomposição. Se chamarmos de $Ψ(ω)$ a transformada de Fourier da função $ψ(t)$ :

$\Psi(\omega) = \int_{-\infty}^{\infty} \psi(t)e^{-i\omega t}dt$

e se $W (a, b)$ for a transformada de wavelet da função $f (t)$ usando a wavelet $ψ(t)$ , então temos que a transformada inversa é dada por:

$f(t) = \frac{1}{C} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{1}{|a|^2}W(a,b)\psi_{a,b}(t)da \ db,$

onde $C = \int_{-\infty}^{\infty} \frac{|\Psi(\omega)|^2}{|\omega|}dw.$

Este parâmetro $C$ necessita ser finito e positivo, o que nos leva a uma nova restrição. Esta restrição sobre o valor de $C$ é a condição de admissibilidade citada anteriormente.

[editar] Análise de wavelet

A análise de wavelet é feita pela aplicação sucessiva da transformada de wavelet com diversos valores para a e b, representando a decomposição do sinal original em diversos componentes localizados no tempo e na freqüência, de acordo com estes parâmetros. Cada wavelet possui melhor ou pior localização nos domínios da freqüência e do tempo, por isso a análise pode ser feita com wavelets diferentes de acordo com o resultado desejado.

A análise wavelet traz consigo uma análise em resoluções múltiplas, onde o nível de resolução é dado pelo índice $a$ . Nesta análise em resoluções múltiplas, geramos uma seqüência de subespaços encaixantes, onde as funções de base numa escala $a 0$ não "enxergam" detalhes de tamanho menor que $2^{-a_0}$ .

[editar] Aplicações

Em geral a transformada contínua de wavelet é usada na análise de sinais, enquanto a sua versão discreta é usada na compressão de dados. Logo, a transformada discreta é usada em engenharia e ciência da computação, enquanto a transformada contínua e usada na pesquisa científica.

As transformadas de wavelet são hoje empregads numa vsta gama de aplicações, substituindo com freqüência a tradicional transformada de Fourier. Diversas áreas da física viram esta mudança de paradigma, incluindo a dinâmica molecular, cálculos ab initio, astrofísica, localização de matriz de densidade, geofísica sísmica, óptica, turbulência e mecânica quântica. Esta mudança também vem ocorrendo no processamento de imagem, análises de pressão sanguínea, ritmo cardíaco e ECG, análise de DNA e proteínas, climatologia, processamento de sinais em geral, reconhecimento de voz, computação gráfica e análise multifractal. Na visão computacional e no processamento de imagens a noção de escala de espaço e operadores de derivadas Gaussianas são vistos como representações multi-escalares canônicas.

Um uso que vem crescendo é na compressão de dados. Como outras transformadas (a Transformada discreta de cosseno sendo a mais comum), as transformações de wavelet podem ser usadas para transformar dados e codificar de forma eficiente os dados transformados, resultando em compressão. O padrão JPEG 2000 para compressão de imagens usa wavelets biortogonais. Isso quer dizer que o memso conjunto de funções pode ser usado tanto na transformada quanto na sua inversa.

[editar] Bibliografia

de Oliveira, Hélio Magalhães. ANALISE DE SINAIS PARA ENGENHEIROS: Uma abordagem via Wavelets. 1.ed. BRASPORT, 2007. ISBN 857452283X
((en)) SALOMON, David. Data Compression: The Complete Reference. 2.ed. Nova Iorque: Springer, 2000.

Retirado de "http://www.girino.org/girinowiki/Wavelet"

Categorias: Computação | Compressão de dados