Revista Colombiana de Fsica, vol. 41, o.

2, Abril 2009

Cifrado de Imgenes utilizando La Transformada Wavelet Fraccional

Images Encryption using Fractional Wavelet Transform
Juan M. Vilardy1, Juan Useche1, Cesar O. Torres1, Lorenzo Mattos1
1

Laboratorio de ptica e Informtica, Universidad Popular del Cesar, Valledupar, Colombia. Recibido XXXX; Aceptado XXXX; Publicado en lnea XXXX

Resumen
En este trabajo se desarrolla una tcnica para el cifrado de imgenes digitales utilizando la Transformada Wavelet Fraccional (Fractional Wavelet Transform, FWT) y mascaras de fases aleatorias. La imagen digital a cifrar es transformada fraccionalmente con la FWT, luego los coeficientes resultantes de la FWT (Aproximacin, Detalles: Horizontal, Vertical y Diagonal) son multiplicados cada uno por diferentes mascaras de fases aleatorias (estadsticamente independientes) y a estos ltimos resultados se le aplica una Transformada Wavelet Inversa (Inverse Wavelet Transform, IWT), obteniendo la imagen digital cifrada. La tcnica de descifrado es la misma tcnica de cifrado en el sentido inverso. Esta tcnica provee ventajas inmediatas de seguridad frente a las tcnicas convencionales; pues en dicha tcnica la familia Wavelet madre y los ordenes fraccionales asociados a la FWT, constituyen claves adicionales que dificultaran el acceso a la informacin a una persona no autorizada (adems de las mascaras de fases aleatorias empleadas), con lo cual se eleva extraordinariamente el nivel de seguridad del cifrado. En este trabajo se desarrolla tambin el soporte matemtico para el uso de la FWT en el algoritmo computacional para el cifrado. Palabras claves: Encriptacin, Transformada Wavelet Fraccional, Mascaras de Fases Aleatorias.

Abstract
In this paper a technique for the coding of digital images is developed using Fractional Wavelet Transform (FWT) and random phase masks. The digital image to encrypt is transformed with the FWT, after the coefficients resulting from the FWT (Approximation, Details: Horizontal, vertical and diagonal) are multiplied each one by different random phase masks (statistically independent) and these latest results is applied an Inverse Wavelet Transform (IWT), obtaining the encrypted digital image. The decryption technique is the same encryption technique in reverse sense. This technique provides immediate advantages security compared to conventional techniques, in this technique the mother wavelet family and fractional orders associated with the FWT are additional keys that make access difficult to information to an unauthorized person (besides the random phase masks used), thereby the level of encryption security is extraordinarily increased. In this work the mathematical support for the use of the FWT in the computacional algorithm for the coding is also developed. Keywords: Encryption, Fractional Wavelet Transform, FWT, Random Phase Masks. 2009 Revista Colombiana de Fsica. Todos los derechos reservados.

1. Introduccin En el presente trabajo se expone una aplicacin para codificar imgenes, usando la transformada wavelet pero en su versin fraccional para incluir como llaves adicionales, los rdenes fraccionales de la transformacin fraccional en el proceso de cifrado, conservando las ventajas en cuanto al uso de los parmetros de desplazamiento y dilatacin. El punto de partida de este trabajo de investigacin se encuen1

tra en los trabajos anteriores de cifrado de imgenes realizados al interior del laboratorio de ptica e Informtica de La Universidad Popular del Cesar [1-2]. 2. La Transformada Wavelet De la misma forma como la transformada de Fourier puede representar a una funcin f(t) en una base de funciones armnicas de exponenciales complejas, la Transformada

Juan M. Vilardy et al.: Cifrado de Imgenes utilizando La Transformada Wavelet Fraccional

Wavelet W(a,b), lleva a la representacin de una funcin f(x) en una base de funciones ortonormales denominadas wavelets madre, definidas como [3]:

W ( a, b) =

f ( x)

1 x b *( )dx ; a a

a > 0 (1)

En la ecuacin (1) el parmetro b corresponde al desplazamiento de la funcin wavelet madre, a es el parmetro de escala. Es decir El conjunto de funciones base en la cual se representa la funcin f(x) (wavelets madre), se construye a partir de una funcin generadora o wavelet madre, a la cual se le va cambiando la escala en un factor a y, a la vez se le va trasladando en una cantidad b [3]. 3. La Transformada Wavelet Fraccional Aplicar una transformada Wavelet Fraccional a una funcin f(x), equivale a la aplicacin de la transformada fraccional de Fourier sobre f(x) y luego realizar una descomposicin wavelet a este resultado [4]. La expresin matemtica para la FWT de una seal f(x) es entonces la siguiente:

transformacin inversa wavelet para obtener la imagen cifrada. Las siguientes ecuaciones describen este proceso: {CA, DH , DV , DD} = DWT2 {F x I F y , familia} (7) CAM = CA Masc1; DHM = DH Masc 2 (8) DVM = DV Masc3; DDM = DD Masc 4 (9) Icifrada = IDWT2 {CAM , DHM , DVM , DDM , familia} (10)

Mascn = exp[ 2irn ( x, y )]

(11)

Donde, Icifrada es la imagen cifrada, CA, DH, DV y DD son las matrices de los coeficientes de aproximacin, de detalle horizontal, de detalle vertical y detalle diagonal que resultan al aplicar la transformada wavelet de una determinada familia en dos dimensiones (DWT2). CAM, DHM, DVM y DDM son las matrices CA, DH, DV y DD despus de haberlas multiplicado por las mscaras aleatorias de fase. Fx y Fy son las transformada fraccionales de Fourier unidimensionales de rdenes x, y; finalmente rn(x,y) es una matriz de nmeros aleatorios del tamao adecuado para multiplicarse por las correspondientes CA, DH, DV y DD (n=1, 2, 3, 4). Para la etapa de descifrado, se aplica el procedimiento en sentido inverso, procurando que cada operacin produzca el efecto inverso en la secuencia. Para invertir la transformada fraccional de Fourier se aplica sta con el orden fraccional negativo (-) [4]. De esta manera, tendremos las siguientes ecuaciones:

W ( a , b) =

B ( x, x' ) f ( x' ) ab ( x)dx' dx (2)

e
i (

sgn( sen
4 [ i (

) 2

B ( x, x' ) = sen

x 2 + x '2 xx ' )2 i ] tan sen

(3)

(4)

{CAM ' , DHM ' , DVM ' , DDM ' } = DWT2{Icifrada, familia ' }
(12) CA' = CAM 'inv ( Masc1); DH ' = DHM 'inv ( Masc 2) (13) DV ' = DVM 'inv ( Masc 3); DD ' = DDM 'inv ( Masc 4) (14)
Irecuperada = F x {IDWT2 {CA' , DH ' , DV ' , DD' , familia '}} F y

Aqu es el orden fraccional de la transformada fraccional de Fourier, la definicin de la FWT es peridica en , con periodo 4. Para Implementar digitalmente la transformada fraccional de Fourier se utiliza el conocido algoritmo de Candan [5-6]:

F ( f [ n ]) = F [ m, n] f [ n ]
n =0

(5) (6)

F [ m, n] = u k [ m]e
k =0

i ka 2

u k [ n]

(15) Cuando la familia usada en el descifrado coincide con la familia usada en el cifrado y, se emplea el mismo orden fraccional en las dimensiones x y y, la imagen descifrada coincide con la descifrada, es decir se recupera la imagen original.

Donde los uk [n] son las hermticas gaussianas discretas. 4. Proceso de Cifrado y Descifrado de una Imagen 4.1 Formulacin Matemtica de los Mtodos de Cifrado y Descifrado Si tenemos una imagen real I(x,y) de M x N pixeles (M y N son nmero de filas y columnas de la matriz que representa a la imagen en escala de grises), el proceso de cifrado de esta imagen se realiza transformndola primero con la fraccional de Fourier bidimensional discreta (DFRT), luego aplicando wavelet discreta (DWT), despus se aplican las mscaras aleatorias de fase y finalmente se realiza una 2

4.2 Proceso de Cifrado Digital

Los algoritmos digitales de cifrado y descifrado de este trabajo se implementaron en la plataforma de programacin Matlab v.7.4, debido a sus grandes facilidades y su alto rendimiento en esta clase de procesos. En la fig. No. 1 se muestra el diagrama de bloques del proceso de cifrado.

4.3 Proceso de Descifrado Digital

El proceso de descifrado es el mismo proceso de cifrado, pero en el sentido inverso utilizando los ordenes fraccionales negativos para la DFRT, las mascaras aleatorias de fase inversas e igual familia wavelet madre, para as obtener la imagen que inicialmente fue cifrada.

rev. col. fs.(c), vol. 41, o. 2, (2009)

Si las llaves utilizadas en el proceso de descifrado no son iguales a las llaves utilizadas en el proceso de cifrado, la imagen original no se recuperar. En la fig. No. 4 se muestra la imagen recuperada cuando se introduce un error de una dcima en los rdenes fraccionales, utilizando 0,877 en lugar de 0,777 y dejando la misma familia wavelet (db4).

Fig. 4. Imagen recuperada cometiendo un error de una dcima en los rdenes fraccionales.

Fig. 1. Diagrama de bloques del proceso de cifrado de una imagen.

5. Resultados Experimentales
Al aplicar el algoritmo para cifrar una imagen digital de MxN pxeles, la imagen digital cifrada oculta la totalidad de su informacin contenida, como se aprecia en la fig. No. 2.

En la anterior figura se observa que con cambios pequeos en las llaves correctas, se genera mucha distorsin en la imagen recuperada. El proceso de descifrado es tan sensitivo al orden fraccional como a la familia wavelet utilizada, esto se observa en la fig. No. 5, en la cual se utilizan los rdenes fraccionales correctos 0,777 (el mismo utilizado en el proceso de cifrado) y se cambia la familia wavelet de db4 por db2.

Fig. 2. Imagen original e imagen encriptada.

Fig. 5. Imagen recuperada utilizando una familia wavelet distinta a la usada para cifrarla.

Cuando se realiza el proceso de descifrado con las llaves correctas (ordenes fraccionales, familia wavelet y mascaras de fase), se logr recuperar la imagen digital original sin prdidas visibles, como lo muestra la fig. No. 3.

6. Anlisis de Resultados
El Error Cuadrtico Medio (MSE, Mean Square Error, ecuacin 16) y la Relacin Seal a Ruido (SNR, Signal-toNoise Ratio, ecuacin 17) entre la imagen digital de entrada y nuestra imagen digital cifrada son calculados para validar la confiabilidad del algoritmo implementado:

MSE =
M

1 M

[I ( x, y) I1 ( x, y)] x =1 y =1

(16)
2

S R = I ( x, y )
x =1 y =1

]2

x =1 y =1

[I ( x, y ) I1 ( x, y )] (17)

Fig. 3. Imagen recuperada usando las llaves correctas.

Donde I(x,y) y I1(x,y) son las matrices elemento de la imagen digital de entrada y nuestra imagen digital cifrada en el pxel (x,y), respectivamente, y M es el tamao de la ima3

Juan M. Vilardy et al.: Cifrado de Imgenes utilizando La Transformada Wavelet Fraccional

gen digital. La sensitividad de los ordenes fraccionales de la FWT son examinados introduciendo errores en estos (de manera individual), dejando fijas las mscaras utilizada y la familia wavelet. El MSE y la SNR son empleados para medir el nivel de proteccin sobre las imgenes digitales cifradas. En esta prueba de desviacin de los rdenes fraccionales de la FWT sobre los valores correctos del proceso de descifrado, se introduce un error que vara entre -0.25 y 0.25, posteriormente para cada variacin se calcula el MSE y la SNR, los resultados se observan en la fig. No. 6 (Los resultados para el orden fraccional en la direccin y son muy similares a los resultados para el orden fraccional en la direccin x). De experimentos computacionales se encontr que los ordenes fraccionales son sensible a una variacin de 1 x 10-4.

gen cifrada intervenga sobre toda la matriz de la imagen recuperada.

(a) (b) Fig. 7. Oclusin de una parte de la imagen cifrada y la respectiva imagen descifrada resultante en un cifrado con FWT: (a) Oclusin de 25% sobre la imagen cifrada e imagen descifrada, (b) Oclusin de 50% sobre la imagen cifrada e imagen descifrada.

Para el anlisis del espacio de clave del algoritmo implementado, se considera el universo de posibilidades de introducir distintos valores de las llaves utilizadas para cifrar la imagen. En este caso sera los dos rdenes fraccionales en las dos dimensiones (x y y), cada uno de los cuales vara en un rango de 0 a 4 con pasos sensibles de 1 x 10-4, la familia de la transformada wavelet y la mscara de fase aleatoria.

Conclusiones
Se ha presentado un nuevo sistema de seguridad de imgenes digitales utilizando la transformada Fraccional wavelet para realizar cifrado de imgenes. Para este sistema de cifrado se aument la seguridad contra ataques de fuerza bruta, por el gran tamao de su espacio de clave. El uso de la FWT en el cifrado de imgenes aumenta enormemente los parmetros de seguridad de la imagen cifrada, debido a la sensibilidad en cualquier cambio que se haga sobre los rdenes fraccionales utilizados (clave numrica) y, adicional a esto, el enorme nmero de posibilidades de usar una familia Wavelet y varias mscaras de fase distintas, incrementan considerablemente la dificultad para cualquier persona que intente descifrar la imagen sin ser autorizado.

(a)

Referencias
[1] (b) Fig. 6. Variacin del MSE y SNR cuando vara el orden fraccional de la FWT en direccin x. J. Vilardy, C. Torres, L. Mattos, Encriptacin en Fase Aplicado a Imgenes Digitales a Color. En: Revista Colombiana de Fsica. Vol. 41, (2008) No.1. L. Brito, C. Torres. Encriptamiento de imgenes digitales a color mediante transformada de Fourier. En: Revista Colombiana de Fsica, Vol. 35, (2003) No.1. Mallat, Stphane. A Wavelet Tour of Signal Processing. San Diego California: Academic Press, 1999. 3 p. Mendlovic, David et al. Fraccional Wavelet Transform. En: Applied optics, Vol. 36, (1997), No.20. A. Bultheel, H. Martnez. Computation of the fractional Fourier transform. En: Applied and Computational Harmonic Analysis. Vol. 16, (2004), 182-202. C. Candan, dFRT: The Discrete Fractional Fourier Transform, A Matlab Program, 1998, <http://www.ee.bilkent.edu.tr/~haldun/dFRT.m>

[2]

[3]

En la siguiente prueba, se chequea la tolerancia a la perdida de la informacin de la imagen cifrada cuando esta es descifrada, para esta prueba se realiza la oclusin del 25% y 50% de los pxeles de la imagen cifrada. La fig. No. 7 muestra las imgenes ocluidas y las correspondientes imgenes recuperadas. De est figura se puede decir que no existe una correspondencia de posicin entre los pxeles de la imagen original y la imagen a descifrar, debido a que las multiplicaciones matriciales aplicadas en el proceso de cifrado y descifrado hacen que los valores de cualquier parte de la ima-

[4] [5]

[6]