

ORIGINAL ARTICLE 

Year : 2019  Volume
: 9
 Issue : 3  Page : 145157 

The ellipselet transform
Zahra Khodabandeh^{1}, Hossein Rabbani^{2}, Alireza Mehri Dehnavi^{2}, Omid Sarrafzadeh^{3}
^{1} Department of Bioelectrics and Biomedical Engineering; Student Research Committee, School of Advanced Technologies in Medicine, Isfahan University of Medical Sciences, Iran ^{2} Department of Bioelectrics and Biomedical Engineering; Medical Images and Signal Processing Research Center, School of Advanced Technologies in Medicine, Isfahan University of Medical Sciences, Iran ^{3} Department of Biomedical Engineering, Mashhad Branch, Islamic Azad University, Mashhad, Iran
Date of Submission  06Oct2017 
Date of Acceptance  22May2019 
Date of Web Publication  29Aug2019 
Correspondence Address: Dr. Alireza Mehri Dehnavi Department of Bioelectrics and Biomedical Engineering, School of Advanced Technologies in Medicine, Isfahan University of Medical Sciences, Isfahan; Medical Images and Signal Processing Research Center, School of Advanced Technologies in Medicine, Isfahan University of Medical Sciences, Isfahan Iran
Source of Support: None, Conflict of Interest: None
DOI: 10.4103/jmss.JMSS_42_17
Background: A fair amount of important objects in natural images have circular and elliptical shapes. For example, the nucleus of most of the biological cells is circular, and a number of parasites such as Oxyuris have elliptical shapes in microscopic images. Hence, atomic representations by twodimensional (2D) basis functions based on circle and ellipse can be useful for processing these images. The first researches have been done in this domain by introducing circlet transform. Methods: The main goal of this article is expanding the circlet to a new one with elliptical basis functions. Results: In this article, we first introduce a new transform called ellipselet and then compare it with other Xlet transforms including 2Ddiscrete wavelet transform, dualtree complex wavelet, curvelet, contourlet, steerable pyramid, and circlet transform in the application of image denoising. Conclusion: Experimental results show that for noises under 30, the ellipselet is better than other geometrical Xlets in terms of Peak Signal to Noise Ratio, especially for Lena which contains more circular structures. However, for Barbara which has fine structures in its texture, it has worse results than dualtree complex wavelet and steerable pyramid.
Keywords: Basis functions, circlet, ellipselet, image denoising, Xlets
How to cite this article: Khodabandeh Z, Rabbani H, Dehnavi AM, Sarrafzadeh O. The ellipselet transform. J Med Signals Sens 2019;9:14557 
Introduction   
In recent years, one of the most considerable approaches in image modeling is using transforms with different basis functions. For example, the atoms of discrete Fourier transform (DFT) are sine and cosine functions which can analyze the frequency component of stationary signals; however, they are not appropriate for analyzing nonstationary signals such as local visual properties in natural images (e.g., edge detection). Hence, although DFT has a highfrequency resolution, it suffers from lack of time/spatial resolution. In order to solve this weakness, the wavelet transform (WT) has been introduced for joint time–frequency analysis of signals.^{[1]} This multiresolution transform is a powerful tool in signal/image analysis and has many applications in denoising, enhancement, and feature extraction. Indeed, DFT and WT are defined for onedimensional (1D) signals, and they cannot efficiently represent geometries in natural images. In fact, these transforms can efficiently model the pointsingularities, and even though they have many applications in signal processing tasks, they are not efficient enough for representing twodimensional (2D) singularities such as edges and curves in natural images. Note that extending 1Ddiscrete wavelet transform (DWT) to 2DDWT does not solve this problem because 2Dseparable DWT which is the tensor product of 1DDWT reconstructs the 2D singularities by aggregation of point singularities around the edges which cannot model the smoothness along the direction of edges. In order to obviate the weaknesses of wavelets in larger dimensions, nonseparable 2D multiscale transforms called geometrical Xlets have been proposed. Ridgelet transform is one of these Xlets which was first introduced by Candes and Donho in 1999.^{[2]} The ridgelet transform represents line singularity in 2D by using Radon transform and maps these singularities to point singularities. In addition to scale and translation parameters which are considered by DWT, the ridgelet transform also considers angle parameter. In addition, to represent the curve singularities in an image, the ridgelet transform has been extended to a new Xlet called curvelet transform.^{[3]} The curvelet transform is one of the multiscale transforms which provides analysis in some windows with different sizes to segment curves as a set of straights in subimages. The main idea behind this transform is decomposing an image to multiresolution subbands and then analyzing each partitioned subband by ridgelet transform.^{[4]} One of the main drawbacks of this transform is its need to discretization (e.g., rotating the curvelet basis functions can be easily performed in polar domain, but it is a challenging task in the Cartesian domain). To solve this problem, contourlet transform has been developed ^{[5]} which is directly defined in the discrete domain. The contourlet transform, at first, employs a multiscale transform to the images for detecting the edges, and then a local directional filter bank is applied for detecting contour segments.^{[6]} In addition to lines and curves, circular and elliptical shapes are frequently seen in some medical images. For example, the nucleus of most of the biological cells is circular, and a number of parasites such as Leishmania and Oxyuris have elliptical shapes. Hence, introducing 2D basis functions based on circle and ellipse in atomic representations can be useful for processing these images. The first researches have been done in this domain by introducing circlet transform.^{[7]} This transform decomposes an image using a set of circles with different radii and a determined width via a DFT filter bank. The main goal of this article is extending the theory of circlet transform to produce a new Xlet transform with elliptical basis functions, i.e., ellipselet transform. In this article, first, the Xlet transforms are explained briefly in Xlet Transforms section. The Ellipselet Transform section is dedicated for introducing ellipselet transform by changing the circular basis functions to elliptical ones, and the results of ellipse detection using ellipselet transform on a simple image is explained in this section. For comparing the new transform with other Xlets, image denoising application with presented Xlets is compared in XLets for Image Denoising section. Finally, this article is discussed in Discussion and Conclusion section.
XLet Transforms   
Usually, transforms decompose an image to a series of elementary waveforms called basis functions or dictionary atoms. Different directional time–frequency dictionaries provide various geometrical Xlet transforms in two or higher dimensions. In this article, we provide a snapshot of a number of geometrical Xlet transforms including 2DDWT, dualtree complex WT (DTCWT), curvelet transform, contourlet transform and steerable pyramid, and circlet transform and introduce a new one namely ellipselet transform which is an extension of the circlet transform by modifying the circular basis function to the elliptical atoms.
Discrete wavelet transform
DWT is a powerful tool for joint time–frequency analysis of signals which decomposes a signal based on a series of basis functions called wavelets. The wavelets are produced by dilation and translation of a mother wavelet. This transform can show point singularities properly, but it is not an optimal tool for the representation of 2D singularities. For image processing, a simple way for using DWT is using tensor product for the extension of 1DDWT to 2DDWT. 2DDWT provides a multiresolution representation by creating four subbands in each scale [Figure 1]. As shown in [Figure 2], the basis functions of 2DDWT can recognize vertical, horizontal, and diagonal spectral features. Although this transform has been used vastly in different applications of signal processing, it has some limitations such as lack of shift invariance and poor directional selectivity.^{[1]}  Figure 1: (a) Filter structure of twodimensionaldiscrete wavelet transform. (b) Twodimensionaldiscrete wavelet transform decomposition in three levels^{[3]}
Click here to view 
Dualtree complex wavelet transform
DTCWT is one of the extensions of WT, which is first introduced by Kingsbury in 1998.^{[9]} This transform uses two real DWTs in parallel as shown in [Figure 3] to create real and imaginary parts of the transform. By producing six subbands in six directions, it has partly improved the deficiencies of 2DDWT. In addition, the redundancy factor of this transform is 2^{d} for ddimensions.^{[10]} [Figure 4] shows the subbands of DTCWT in the direction of ±15, ±45, and ±75.  Figure 3: Analysis of filter bank for the dualtree complex discrete wavelet transform^{[10]}
Click here to view 
 Figure 4: The subbands of real oriented twodimensional dualtree complex discrete wavelet transform in (a) space domain (b) Fourier spectrum^{[10]}
Click here to view 
Discrete curvelet transform
The ridgelet transform, introduced by Cand'es and Donho,^{[2]} is an efficient tool for representing line singularities. Moreover, for demonstrating curves in images, this transform has been developed to curvelet transform. The curvelet is a multiscale transform which uses scaling, translating, and rotating parameters to create its basis functions. The main idea behind the curvelet transform is decomposing an image into sets of wavelet subbands and then analyzing each subband with the local ridgelet transform.^{[2]} The block diagram of the first version of curvelet transform and its subbands are shown in [Figure 5] and [Figure 6], respectively.
Discrete contourlet transform
The contourlet transform was first introduced by Do and Vetterli ^{[5]} to overcome the limitations of curvelet transform in discretization. This transform can detect 2Dgeometrics in images in two steps: first, it decomposes an image into a set of subbands by Laplacian Pyramid and then applies a series of directional filter banks to the image. The basis functions are oriented in different scales and orientations.^{[6]} The block diagram of the Contourlet transform and its subbands are illustrated in [Figure 7] and [Figure 8], respectively.
Steerable pyramids
The steerable pyramid is a linear multiscale and multiresolution transform based on angular and radial decompositions.^{[13]} According to [Figure 9], first, the image is divided into low and highpass subbands. Second, the lowpass subband is separated into a series of oriented bandpass subbands and a lowpass subband. Then, the lowpass subband is downsampled by a factor of 2, and this process will be repeated. If the transform has k orientation bands, it would be overcomplete by a factor of 4k/3, showing redundancy.^{[14]}
Circlet transform
The circlet transform is a robust tool to detect circular objects in images in which the binary image segmentation is not needed. The transform decomposes an image to a set of circles with different radii and a determined width using a DFT filter bank. This decomposition is defined in Fourier domain using the following definitions which are very close to the one introduced by Chauris et al.^{[7]} The circlet parameters are described by a central position (x_{0}, y_{0}), radius (r_{0}), and central frequency content (f_{0}). All circlet components C_{μ} (x, y) can be obtained by a reference circlet C_{ref} (x, y) which can be shifted or changed in radius and central frequency content of the circlet. The circlet function is defined by the following equation:
Where, and Ω is an oscillating function such as wavelet function to distinguish discontinuities. Practically, C_{μ} is defined in 2D Fourier domain.^{[7]} In circlet decomposition, an image f(x, y) is broken down into a sum of basic functions C_{μ} as shown in Eq. 2:
The circlet transform is a tight frame system, and the amplitudes A_{μ} are created by a scalar product as shown in Eq. 3:
From a practical point of view, the circlet coefficients can be defined in the Fourier domain using Parseval's theorem as shown below:
Where f^ is the 2D Fourier transform of f and f* is the conjugate of f. Hence, the circlet transform is defined in 2D Fourier domain with the definition of; Ĉμ (w_{1}, w_{2}), the Fourier transform of Cμ plays an important role in the performance of circlet transform.
In order to construct circularshaped filters, 1D filters F_{k} and 2D filters G_{k} are used such that for all ω and (w_{1}, w_{2}), the following perfect reconstruction conditions are satisfied:
In these conditions, the filter F_{k} is defined as shown in Eq. 6, where N is the number of filters and
In order to construct circular shape in spatial domain, the 2D filter G_{k} is defined by a phase delay of 1D filter as F_{k} shown in Eq. 7.
In Eq. 7, ω (w_{1}, w_{2}) = ω and is defined as follows:
By the definition of G_{k}, the Fourier transform of a circlet is described in Eq. 9.
Where x_{c}= (x 0, y 0) is the central position, and r_{0} is the radius of the circlet. [Figure 10] shows the magnitude and phase of basis functions in circlet transform.^{[16]}  Figure 10: (a) The magnitude of basis functions of circlet transform. (b) The phase of basis functions of circlet transform
Click here to view 
Ellipselet Transform   
The main goal of this study is expanding the circlet transform to a new one with elliptical basis functions. Therefore, we have to design a new atomic representation system with elliptical basis functions. As mentioned in Eq. 7, in F_{k}ω and for ω = (w_{1}, w_{2}), ω is defined by the Eq. 8 that it constructs the circular shape of basis functions in circlet. Now, we introduce a new norm definition as follows:
In this equation, by considering different values of k_{1} and k_{2} and α parameters, elliptical basis functions in different sizes and directions are built. We assume that we have a new transform with four basis functions in four different directions as shown in [Figure 11].  Figure 11: Basis functions of the ellipselet transform in four directions and four frequency contents. (ad) the basis function at angles 0°, 90°, −45°, and 45°, respectively
Click here to view 
For detecting ellipses in images, we need to know the position, angle, and size of the major and minor axes of ellipse. We assume a simple ellipse as shown in [Figure 12].  Figure 12: A simple ellipse with major axis: 50 pixels, minor axis: 20 pixels, angle: 90, center position: (250, 250), size of image: (512, 512)
Click here to view 
The produced subbands of [Figure 12] using the introduced basis functions in [Figure 11] are represented in [Figure 13]. As shown, the first row which is related to basis functions at 0° is similar to the original image. By converting the subbands to binary images and using morphological operations, the new subimages as shown in [Figure 14] are achieved.  Figure 13: The subbands of ellipselet transform on the ellipse of Figure 12
Click here to view 
 Figure 14: Produced subimages of Figure 13 after binerization and applying morphological operators
Click here to view 
It is obvious from [Figure 14] that only subbands in the first row are similar to the original image which is related to the basis functions with zero direction. In order to find the desired size and location of the ellipse, radon transform can be employed.^{[17]} If we use radon transform for each subband in the angle of perpendicular direction of its basis, we achieve signals as shown in [Figure 15]. We consider a percentage of maximum of peaks in signals for thresholding of signals as shown in [Figure 15]. By choosing the subband in which all the four subimages are above the threshold, the first subband is extracted that is related to the basis of zero direction and shows that the main ellipse is in vertical direction at 90°. Now, for detecting the position of the main ellipse, all row/column signals of the selected subband are plotted and the middle of two maximum points is selected as the x/y position of the center of the desired ellipse [Figure 16].  Figure 15: (a) Radon transform of the first row of Figure 14 in the direction of 90°. (b) Radon transform of the second row of Figure 14 in the direction of 0°. (c) Radon transform of the third row of Figure 14 in the direction of +45°. (d). Radon transform of the fourth row of Figure 14 in the direction of −45°
Click here to view 
After determining the center position and angle of the ellipse, we calculate the dimensions of the ellipse by using two radon transforms in the direction of the ellipse and its perpendicular line. Considering the obtained parameters of ellipse, we can draw the extracted ellipse by ellipselet transform. [Figure 17] shows the contour of the produced ellipse in blue color. This figure shows that the position and angle of the ellipse have been detected correctly.  Figure 17: (a) Detected ellipse with the new transform. (b) The zoomedin detected ellipse of part
Click here to view 
XLets for Image Denoising   
In general, digital images suffer from noise due to acquisition/transmission process and shortcoming of capturing modalities and receivers. In addition to linear methods such as Wiener filtering, nonlinear techniques including applying thresholding/shrinkage functions in transform domains have been reported in recent years.^{[3],[5],[9],[18],[19]} In general, in image transformbased denoising approaches, two important issues should be considered. The first one is choosing proper transform and the other is selecting proper thresholding function. In this article by using a proper thresholding function, the performance of geometrical Xlets such as 2DDWT, DTCWT, curvelet, contourlet, steerable pyramid, circlet transform, and the proposed transform of ellipselet was compared in reducing additive white Gaussian noise from natural images.
Experimental Results   
Different frameworks can be considered for image denoising. Usually, reported studies in Bayesian frameworks outperform others. For example, Rabbani in 2009^{[19]} used local Laplace pdf and maximum a posteriori (MAP) estimator in steerable pyramid domain and showed better performance in comparison with the other stateoftheart denoising methods such as Bayes least squared Gaussian scale mixture technique.^{[20]} The statistical features of images can be simplified in sparse domains because of some properties of transforms such as sparsity. Hence, the main features of transformed image can be represented by a few large coefficients, and the remained coefficients are approximately around zero. Hence, there is a large peak at zero in the histogram of subbands and its tails goes to zero slower. It means that their distribution is close to the Laplace pdf and far from the Gaussian pdf. In this study, we first show the histogram of the Xlet coefficients in a specific subband and then the Laplacian and Gaussian pdfs are fitted to the histogram [Figure 18], and the goodness of fit of each pdf is reported. Because Laplacian pdf is well fitted to the histograms, using MAP estimator and Laplacian prior, the soft thresholding with a threshold of is obtained,^{[19]} where σ_{n} is the standard deviation of noise and σ_{k} is the standard deviation of noisefree image in the k ^{th} subband. Experiments performed on standard gray scale images of Lena, Barbara, and Boat at a resolution of 512 × 512 pixels corrupted by additive Gaussian noise with different levels to compare the performance of geometrical Xlets transforms in image denoising.  Figure 18: The histogram of the XLet coefficients in a specific subband of Lena image and the best Gaussian (green plot) and Laplace pdf (red plot) fitted to this histogram with Chisquare test
Click here to view 
According to [Figure 19], by testing different combinations of ellipselet subbands (for n = 4 in Eq. 6), the best signaltonoise ratio is obtained on the first, third, and fourth subbands of ellipselet ( [1, 3, 4]). Hence, for denoising by circlet and ellipselet transforms, the second subband, which includes the main global information of the image, remains unchanged, and soft thresholding is applied on the first, third, and fourth subbands. However, similar to the usual procedure of transformbased denoising methods for other Xlets, in this article, the lowpass subband is kept unchanged and the thresholding function is applied on other subbands.  Figure 19: Comparison among thresholding on different subbands of ellipselet. (a) Peak signaltonoise ratio of denoized Lena in different subband thresholding. (b) Peak signaltonoise ratio of denoized Barbara in different subband thresholding. (c) Peak signaltonoise ratio of denoized boat in different subband thresholding
Click here to view 
The efficiency of the presented methods is first visually evaluated as shown in [Figure 20] and [Figure 21]. Then, the experimental results are numerically compared in terms of peak signaltonoise ratio (PSNR) [Figure 22] and Structural Similarity Index (SSIM) [Table 1] using the following definitions:  Figure 20: Results of image denoising of Lena with different Xlets for sigma = 15
Click here to view 
 Figure 21: Results of image denoising of boat with different Xlets for sigma = 25
Click here to view 
 Figure 22: The results of peak signaltonoise ratio for Lena and Barbara images denoized by different Xlets
Click here to view 
 Table 1: The results of structural similarity index for different denoising methods
Click here to view 
Where MSE is the mean squared error and MAX_{t} is the maximum possible pixel value of the image.
The SSIM shows structural similarities between two images and defines as follows.^{[21]}
Discussion   
In this article, we introduced a new Xlet transform namely ellipselet transform by using elliptical basis functions. We showed that this transform can correctly detect a simple ellipse at 90°. For ellipses in three other directions of 0°, −45°, and +45°, we also did the same work described above and achieved the desired answers. Hence, the new transform with elliptical basis functions introduced in this article can detect the direction and position of simple ellipses in images. However, this new transform has some limitations. For example, it has limited number of basis functions, and we just tested it on simple images with one ellipse. Considering new criterions and expanding this transform by using more basis functions (e.g., 8 basis functions in angles 0°, ±22.5°, ±45°, ±67.5°, and 90°) is suggested in future studies.
In this article, we also compared the introduced ellipselet transform with other Xlet transforms in reducing noise from natural images. As shown in [Figure 22], DTCWT usually outperforms others in terms of PSNR. The main reason is that the DTCWT has a good directional selectivity and perfect reconstruction property which makes it a proper tool for denoising applications. According to [Table 1], DTCWT has the maximum SSIM for Boat and Barbara, and steerable pyramid has the maximum SSIM for Lena in different levels of noise. The reason of good performance of steerable pyramid is that the steerable pyramid operates based on a polarseparable decomposition in the frequency domain that allows independent representation of scale and orientation, and the representation is translation and rotation invariant which makes it proper to analyze image structures and edge preservation. Moreover, as shown in [Figure 22], for noises under 30, the ellipselet is better than others, especially for Lena which contains more circular structures. However, for Barbara which has fine structures in its texture, it has worse results than DTCWT and steerable pyramid.
Different results in different images by different levels of noise show that Xlets' performance is related to the image content and level of noise. As illustrated in [Figure 23], the basis functions of different Xlets produce different frequency features, and each of them partitions the 2D frequency plane in a different way. Based on the location of components of image in 2D frequency plane and how they match with basis functions, any of these transforms could be more beneficial. About the newly introduced transform, ellipselet, we can see that it has better results than circlet transform and both of them have good results in lower levels of noise, especially for images containing circular patterns.  Figure 23: Frequency partitioning of: (a) Dualtree complex wavelet, (b) curvelet,^{[22]} (c) contourlet, (d) steerable pyramid, (e) circlet, (f) ellipselet in the directions of 0, 90, +45, −45
Click here to view 
Financial support and sponsorship
None.
Conflicts of interest
There are no conflicts of interest.
Zahra Khodabandeh is a PhD candidate at Isfahan University of Medical Sciences. She received her BSc and MSc degree in biomedical engineering from Isfahan University and Isfahan University of Medical Sciences, Isfahan, Iran in 2011 and 2017, respectively. Her current research interests include medical image analysis, machine learning and deep learning.
Email: shimakhodabandeh@yahoo.com
Hossein Rabbani received his BSc degree in electrical engineering from Isfahan University of Technology, Isfahan, Iran, in 2000, and his MSc and PhD degrees in bioelectrical engineering from Amirkabir University of Technology, Tehran, Iran, in 2002 and 2008, respectively. He is now a full professor in the Biomedical Engineering Department, Isfahan University of Medical Sciences, Isfahan. His research interests are medical image analysis and modeling, signal processing, sparse transforms, and image restoration.
Email: h_rabbani@med.mui.ac.ir
Alireza Mehri Dehnavi received his BSc degree in electronic engineering from Isfahan University of Technology in 1988. He received his MSc degree in measurement and instrumentation from Indian Institute of Technology Roorkee, India, in 1992 and his PhD in medical engineering from Liverpool University in 1996. He is a full professor in the Biomedical Engineering Department, Isfahan University of Medical Sciences, Isfahan, Iran. His research interests are medical optics, devices and signal, and image processing.
Email: mehri@med.mui.ac.ir
Omid Sarrafzadeh received his PhD in biomedical engineering from Isfahan University of Medical Sciences, Isfahan, Iran. He is assistant professor in Islamic Azad University Mashhad Branch. His research interests are image and signal processing and pattern recognition.
Email: o.sarrafzade@gmail.com
References   
1.  Daubechies I, Bates BJ. Ten Lectures on Wavelets. J Acoust Soc Am 1993;93:1671. 
2.  Candes EJ, Donoho DL. Ridgelets: A key to higherdimensional intermittency? Philos Trans R Soc A Math Phys Eng Sci 1999;357:2495509. 
3.  Alzubi S, Islam N, Abbod M. Multiresolution analysis using wavelet, ridgelet, and curvelet transforms for medical image segmentation. Int J Biomed Imaging 2011;2011:136034. 
4.  Starck JL, Candès EJ, Donoho DL. The curvelet transform for image denoising. IEEE Trans Image Process 2002;11:67084. 
5.  Do MN, Vetterli M. The contourlet transform: An efficient directional multiresolution image representation. IEEE Trans Image Process 2005;14:2091106. 
6.  Yang L, Guo BL, Ni W. Multimodality medical image fusion based on multiscale geometric analysis of contourlet transform. Neurocomputing 2008;72:20311. 
7.  Chauris H, Karoui I, Garreau P, Wackernagel H, Craneguy P, Bertino L. The circlet transform: A robust tool for detecting features with circular shapes. Comput Geosci 2011;37:33142. 
8.  Amini Z, Rabbani H. Classification of medical image modeling methods: A review. Curr Med Imaging Rev 2016;12:1. [Last cited on 2019 May 20]. 
9.  Kingsbury N. The DualTree Complex Wavelet Transform: A New Efficient Tool for Image Restoration and Enhancement; 1998. 
10.  Selesnick IW, Baraniuk RG, Kingsbury NC. The dualtree complex wavelet transform. IEEE Signal Process Mag 2005;22:12351. 
11.  Lambert P, Pires S, Ballot J, García RA, Starck JL, TurckChièze S. Curvelet analysis of asteroseismic data I: Method description and application to simulated sunlike stars. Astron Astrophys 2006;454:10217. 
12.  Satheesh S, Prasad K. Medical image denoising using adaptive threshold based on contourlet transform. Adv Comput Int J 2011;2:528. 
13.  Simoncelli EP, Freeman WT. The steerable pyramid: A flexible architecture for multiscale derivative computation. In: Proceedings, International Conference on Image Processing. Washington, DC. IEEE Comput Soc Press; 1995. p. 4447. 
14.  Simoncelli EP, Freeman WT, Adelson EH, Heeger DJ. Shiftable multiscale transforms. IEEE Trans Inf Theory 1992;38:587607. 
15.  Wang J. Exposure fusion based on steerable pyramid for displaying high dynamic range scenes. Opt Eng 2009;48:117003. 
16.  Sarrafzadeh O, Dehnavi AM, Rabbani H, Ghane N, Talebi A. Circlet based framework for red blood cells segmentation and counting. In: 2015 IEEE Workshop on Signal Processing Systems (SiPS). IEEE; 2015. p. 16. 
17.  Comaniciu D, Meer P. Mean shift: A robust approach toward feature space analysis. IEEE Trans Pattern Anal Mach Intell 2002;24:60319. 
18.  Hamdi M. A comparative study in wavelets, curvelets and contourlets as denoising biomedical Images. Image Process Commun 2011;16:1320. 
19.  Rabbani H. Image denoising in steerable pyramid domain based on a local Laplace prior. Pattern Recognit 2009;42:218193. 
20.  Portilla J, Strela V, Wainwright MJ, Simoncelli EP. Image denoising using scale mixtures of Gaussians in the wavelet domain. IEEE Trans Image Process 2003;12:133851. 
21.  Wang Z, Bovik AC, Sheikh HR, Simoncelli EP. Image quality assessment: From error visibility to structural similarity. IEEE Trans Image Process 2004;13:60012. 
22.  Candès E, Demanet L, Donoho D, Ying L. Fast discrete curvelet transforms. Multiscale Model Simul 2006;5:86199. 
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8], [Figure 9], [Figure 10], [Figure 11], [Figure 12], [Figure 13], [Figure 14], [Figure 15], [Figure 16], [Figure 17], [Figure 18], [Figure 19], [Figure 20], [Figure 21], [Figure 22], [Figure 23]
[Table 1]
