DCT
DCT
DCT
I NTRODUCCIN
II .
G(r) =
k=0
r = 0,1,...,N-1
Donde:
1
c(r) = para r =0
2
y c(r) = 1 para r=1,2,...,N-1
Esta transformada tiene aplicaciones en procesamiento de
imagenes. El mtodo convensional para aplicar la DCT usa un
algoritmo de la DFT de 2N muestras.
Aplicaremos un mtodo mas eficiente de evaluar la DCT.
Concidere el clculo de:
N 1
h
2c(r) X
r i
x(k) cos (2k + 1)
N
2N
N
1
X
F (r) =
x(k) cos[(2k + 1)
k=0
r
]
2N
Para: r = 0,1,...,N-1
Si asumimos que N es parm definimos una nueva secuencia y(k)
por:
y(k) = x(2k)
y(N 1 k) = x(2k + 1)
con k = 0,1,...,(N/2)-1.
Con esto F(r)se reescribe y nos queda de la siguiente manera:
(N/2)1
F (r) =
y(k) cos[(4k + 3)
k=0
r
]+
2N
(N/2)1
y(N 1 k) cos[(4k + 3)
k=0
r
]
2N
F (r) =
Figura 1. Concentracin de pontencia de una DCT-II bidimensional comparada con la concentracin de potencia de una DFT tambin bidemensional.
N
1
X
y(k) cos[(4k + 3)
k=0
r
]
2N
Donde:
H(r) = ejr/2N
N
1
X
y(k)ej2rk/2N
k=0
H(N r) = j[H(r)]
Donde * denota una conjugacion compleja. Y separando las
partes real e imaginaria tenemos:
F (r) = Re[H(r)]
F (N r) = Im[H(r)] Para:
r=0,1,...,(N/2)
El clculo de la IDTC se puede realizar de la siguiente manera:
x(k) =
N
1
X
r=0
h
r i
c(r)G(r) cos (2k + 1)
2N
,
para k=0,1,2,...,N-1
Las muestras pares son evaluadas como:
x(2k) =
N
1
X
[c(r)G(r)ejr/2N ]ej2rk/N
r=0
,
para k=0,1,2,...,N-1
Para obtener los puntos impares note que:
x(2k + 1) = x[2(N 1 k)]
A LGORITMOS DCT
A.
Cuantificacin
A.
E JEMPLO DCT-1
, n=0 y N-1
2
[n] =
1 ,1 n N 2
Para encontrar
x
1 [n] = x [((n))2N 2 ] + x [((n))2N 2 ]
Se necesita
1
x [n] = [( )(4) 1(3)
2
Con 2N 2 se tiene 8 2 = 6
1
0 0]
2
1
2 3]
x [n6 ] = [2 0 0
2
x
1 [n] = [4 3 2 1 2 3]
x [n6 ] = [2
N
1
X
[n]x[n]cos(
n=0
xc1 [k] = 2
3
X
kn
)
N 1
[n]x[n]cos(
n=0
kn
)
3
x
1 [n] = x [((n))2N 2 ] + x [((n))2N 2 ]
0nN 1
0n3
1
1
xc1 [0] = 2( (4))cos(0)+1(3)cos(0)+1(2)cos(0)+ (1)cos(0) = 15
2
2
2
1
1
c1
x [1] = 2( (4))cos(0)+1(3)cos( )+1(2)cos( )+ (1)cos() = 4
2
3
3
2
1
2
4
1
c1
x [2] = 2( (4))cos(0)+1(3)cos( )+1(2)cos( )+ (1)cos(2) = 0
2
3
3
2
1
1
c1
x [3] = 2( (4))cos(0)+1(3)cos()+1(2)cos(2)+ (1)cos(3) = 1
2
2
xc1 [k] = [15
1
( )(1)]
2
1(2)
1]
A.
N
1
X
xN [n]ej
2kn
N
0k N 1
n=0
x[k] =
5
X
x6 [n]ej
2kn
6
0k5
n=0
x[k] = [15
Confirmando que:
1]
A PLICACIONES
VII .
C ONCLUSIONES
B IBLIOGRAFA