Mahotas: reconstrucción de la imagen a partir de la imagen wavelet transformada de Daubechies

En este artículo veremos cómo podemos reconstruir la imagen a partir de la imagen transformada de daubechies wavelet en mahotas. En general, las ondículas de Daubechies se eligen para tener el mayor número A de momentos de fuga (esto no implica la mejor suavidad) para un ancho de soporte dado 2A. Hay dos esquemas de nombres en uso, DN que usa la longitud o el número de toques y dbA que se refiere al número de momentos de fuga. Entonces D4 y db2 son la misma transformada wavelet. 
En este tutorial usaremos la imagen «luispedro», a continuación se muestra el comando para cargarla. 
 

mahotas.demos.load('luispedro')

A continuación se muestra la imagen de luispedro 
 

Para ello utilizaremos el método mahotas.idaubechies 
 

Sintaxis: mahotas.idaubechies (img, ‘D8’)
Argumento: toma el objeto de imagen y la string, es decir, uno de ‘D2’, ‘D4’, … ‘D20’ como argumento
Retorno: Devuelve el objeto de imagen 
 

Nota: la imagen de entrada debe filtrarse o debe cargarse como gris
. Para filtrar la imagen, tomaremos el objeto de imagen que es numpy.ndarray y lo filtraremos con la ayuda de la indexación, a continuación se muestra el comando para hacer esto.
 

image = image[:, :, 0]

Ejemplo 1: 
 

# importing various libraries
import numpy as np
import mahotas
import mahotas.demos
from mahotas.thresholding import soft_threshold
from pylab import imshow, show
from os import path
 
# loading image
f = mahotas.demos.load('luispedro', as_grey = True)
 
 
 
# Transform using D8 Wavelet to obtain transformed image t
t = mahotas.daubechies(f, 'D8')
 
# Discard low-order bits:
t /= 8
t = t.astype(np.int8)
 
 
# showing transformed image
print("Transformed Image")
imshow(t)
show()
 
 
# reconstructed image
r = mahotas.idaubechies(t, 'D8')
 
# showing image
print("Reconstructed Image")
imshow(r)
show()

Producción : 
 

Ejemplo 2: 
 

# importing required libraries
import mahotas
import numpy as np
from pylab import imshow, show
import os
 
 
# loading image
img = mahotas.imread('dog_image.png')
 
# filtering image
img = img[:, :, 0]
 
 
# Transform using D8 Wavelet to obtain transformed image t
t = mahotas.daubechies(img, 'D8')
 
 
# showing transformed image
print("Transformed Image")
imshow(t)
show()
 
# reconstructed image
r = mahotas.idaubechies(t, 'D8')
 
# showing image
print("Reconstructed Image")
imshow(r)
show()

Producción :