In this article we will see how we can reconstruct image from the transformed image of daubechies wavelet in mahotas. In general the Daubechies wavelets are chosen to have the highest number A of vanishing moments, (this does not imply the best smoothness) for given support width 2A. There are two naming schemes in use, DN using the length or number of taps, and dbA referring to the number of vanishing moments. So D4 and db2 are the same wavelet transform.
In this tutorial we will use “luispedro” image, below is the command to load it.
mahotas.demos.load('luispedro')
Below is the luispedro image
In order to do this we will use mahotas.idaubechies method
Syntax : mahotas.idaubechies(img, ‘D8’)
Argument : It takes image object and string i.e one of ‘D2’, ‘D4’, … ‘D20’ as argument
Return : It returns image object
Note : Input image should be filtered or should be loaded as grey
In order to filter the image we will take the image object which is numpy.ndarray and filter it with the help of indexing, below is the command to do this
image = image[:, :, 0]
Example 1:
Python3
# importing various librariesimport numpy as npimport mahotasimport mahotas.demosfrom mahotas.thresholding import soft_thresholdfrom pylab import imshow, showfrom os import path# loading imagef = mahotas.demos.load('luispedro', as_grey = True)# Transform using D8 Wavelet to obtain transformed image tt = mahotas.daubechies(f, 'D8')# Discard low-order bits:t /= 8t = t.astype(np.int8)# showing transformed imageprint("Transformed Image")imshow(t)show()# reconstructed imager = mahotas.idaubechies(t, 'D8')# showing imageprint("Reconstructed Image")imshow(r)show() |
Output :
Example 2:
Python3
# importing required librariesimport mahotasimport numpy as npfrom pylab import imshow, showimport os# loading imageimg = mahotas.imread('dog_image.png')# filtering imageimg = img[:, :, 0]# Transform using D8 Wavelet to obtain transformed image tt = mahotas.daubechies(img, 'D8')# showing transformed imageprint("Transformed Image")imshow(t)show()# reconstructed imager = mahotas.idaubechies(t, 'D8')# showing imageprint("Reconstructed Image")imshow(r)show() |
Output :
