alexa Image Enhancement by Wavelet with Principal Component Analysis | OMICS International
ISSN: 2277-1891
International Journal of Advance Innovations, Thoughts & Ideas
Make the best use of Scientific Research and information from our 700+ peer reviewed, Open Access Journals that operates with the help of 50,000+ Editorial Board Members and esteemed reviewers and 1000+ Scientific associations in Medical, Clinical, Pharmaceutical, Engineering, Technology and Management Fields.
Meet Inspiring Speakers and Experts at our 3000+ Global Conferenceseries Events with over 600+ Conferences, 1200+ Symposiums and 1200+ Workshops on
Medical, Pharma, Engineering, Science, Technology and Business

Image Enhancement by Wavelet with Principal Component Analysis

Vikas D .Patil*, Sachin D Ruikar

Department of Electronics & Telecommunication Dept., Sinhgad Academy of Engg, Pune, India

*Corresponding Author:
Vikas D .Patil
Department of Electronics & Telecommunication Dept
Sinhgad Academy of Engg, Pune, India
E-mail:[email protected]

Visit for more related articles at International Journal of Advance Innovations, Thoughts & Ideas

Abstract

 This paper demonstrate the dimensionality of image sets with Wavelet using principal component analysis on wavelet coefficients to maximize edge energy in the reduced dimension images. Large image sets, for a better preservation of image local structures, a pixel and its nearest neighbors are modeled as a vector variable, whose training samples are selected from the local window by Local Pixel Grouping (LPG).
The LPG algorithm guarantees that only the sample blocks with similar contents are used in the local statistics calculation for PCA transform estimation, so that the image local features can be well preserved after coefficient shrinkage in the PCA domain to remove the random noise. The LPG-PCA Enhance procedure is used to improve the image quality.
The wavelet thresholding methods used for removing random noise has been researched extensively due to its effectiveness and simplicity. However, not much has been done to make the threshold values adaptive to the spatially changing statistics of images. Such adaptivity can improve the wavelet thresholding performance because it allows additional local information of the image (such as the identification of smooth or edge regions) to be incorporated into the algorithm of a damaged or target region in addition to shape and texture properties

Keywords

Wavelet, Wavelet Transform (WT), Local Pixel Grouping (LPG), Principal Components Analysis (PCA).

Introduction

Principal Components Analysis (PCA) is the way of identifying patterns in data, and expressing the data in such a way as to highlight their similarities and differences. Since patterns in data can be hard to find in data of high dimension, where the luxury of graphical representation is not available PCA is a powerful tool for analyzing data. The other main advantage of PCA is that once you have found these patterns in the data and you compress the data i.e. by reducing the number of dimensions without much loss of information. This technique used in image compression. To improve the image quality image enhancement is an essential step. As a primary low-level image processing procedure, random noise removal has been extensively studied and many enhancement schemes have been proposed, from the earlier smoothing filters and frequency domain denoising methods [1] to the lately developed wavelet [2][3][11], curvelet [12] and ridgelet [13] based methods, sparse representation [14].With the rapid development of modern digital imaging devices and their increasingly wide applications in our daily life, to overcome the problem with conventional WT [15] we use Wavelet PCA. Wavelets are an efficient and practical way to represent edges and image information at multiple spatial scales. Image features at a given scale, such as houses or roads, can be directly enhanced by filtering the wavelet coefficients. Wavelets may be a more useful image representation than pixels. Hence, we consider PCA dimensionality reduction of wavelet coefficients in order to maximize edge information in the reduced dimensionality set of images. The wavelet transform will take place spatially over each image band, while the PCA transform will take place spectrally over the set of images. Thus, the two transforms operate over different domains. Still, PCA over a complete set of wavelet and approximation coefficients will result in exactly the same eigenspectra as PCA over the pixels.

To overcome the problem of WT, in [16] Muresan and Parks proposed a spatially adaptive principal component analysis (PCA) based denoising scheme, which computes the locally fitted basis to transform the image. Elad and Aharon [14] [17] proposed sparse redundant representation and (clustering -singular value decomposition) K-SVD based denoising algorithm by training a highly over-complete dictionary. Foi et al. [18] applied a shape-adaptive discrete cosine transform (DCT) to the neighborhood, which can achieve very sparse representation of the image and hence lead to effective denoising. All these methods show better denoising performance than the conventional WT-based denoising algorithms. The recently developed non-local means (NLM) approaches use a very different philosophy from the above methods in random noise removal. The idea of NLM can be traced back to [19], where the similar image pixels are averaged according to their intensity distance. Similar ideas were used in the bilateral filtering methods [20] [21], where both the spatial and intensity similarities are exploited for pixel averaging. In [22], the NLM denoising frame- work was well established. Each pixel is estimated as the weighted average of all the pixels in the image, and the weights are determined by the similarity between the pixels. This scheme was improved in [23], where the pair-wise hypothesis testing was used in the NLM estimation. Inspired by the success of NLM methods, recently Dabov et al. [24] proposed a collaborative image denoising scheme by patch matching and sparse 3D transform. They searched for similar blocks in the image by using block matching and grouped those blocks into a 3D cube.

In the system we give an input as a blur image then apply algorithm Local Pixel Grouping then performs the principal component analysis (PCA) on the blur image. The output from PCA is the eigenimage and the eigenvectors. Apply soft thresholding on a PCA component, the number of them can be selected by the user, and the reconstruction quality in the inverse PCA (IPCA) depends on that number. The eigenimages and the eigenvectors are entropy coded. Then the wavelet transform (WT) is applied to that residual. Then apply thresholding on wavelet coefficient and the reconstruction quality in the Inverse Wavelet Transform (IWT). And finally get the output enhanced image. The rest of the paper is structured as follows. Section 2 briefly reviews the experimental results. Section 3 presents discussion of procedure of PCA, Wavelets transform , Implementation paper. Section 4 concludes the paper.

Experimental Results

In experimenting result, we try the different experiment to prove the superiority of proposed method. The experimental processes are conducted: it shows Enhancement images. We use the several different characteristic of the images, prove that our Local Pixel Grouping be used provides better results than other existing technique. In order to test the quality of our proposed image Enhancement method, we used various images, including photos, scenery, and artistic compositions. Also we have used different wavelet transforms.

Principle Component Analyses

Denote by x=[x1,x2 ,. . . xm]T an m-component vector variable and denote by

image (1)

.the sample matrix of x, where xji , j=1,2,..,n, are the discrete samples of variable xi , i=1,2,..,m. The ith row of sample matrix X, denoted by

xi=[xi1,xi2 ,. . . xin] (2)

is called the sample vector of xi. The mean value of Xi is calculated as

image (3)

and then the sample vector Xi is centralized as

image (4)

Where image accordingly, the centralized matrix of X is

image (5)

Finally, the co-variance matrix of the centralized dataset is calculated as

image (6)

The goal of PCA is to find an ortho-normal transformation matrix P to de-correlate image, i.e.image so that the covariance matrix of Y is diagonal. Since the covariance matrix X is symmetrical, it can be written as:

image (7)

where ø = [ ø 1 ø 2 ……. ø n] is the m*m orthonormal eigenvector matrix and ^ = diag{λ1λ2..λm}is the diagonal eigenvalue matrix with λ1>= λ2 >=……….>=λm the terms ø1 ø2 ……. øn and the λ1,λ2..λm are the eigenvectors and eigenvalues of X. By setting

P= ø T (8)

image can be decorrelated.

Local Pixel Grouping (LPG)

The LPG, for an under lying pixel to be denoised, we set a K*K window centered on it and denote by x=[x1………xm]T,m=K2,the vector containing all the components within the window. Since the observed image is noise corrupted, we denote by

Xv=X+v (9)

Grouping the training sample similar to the central K*K block in the L*L training window is indeed a classification problem and thus different grouping methods, such as block matching, correlation-based matching , K-mean s clustering, etc, can be employed based on different criteria. Among them the block matching method may be the simplest yet very efficient one. There are totally (L-K+1)2 possible training blocks of xn in the L*L training window. We denote by image the column sample vector containing the pixel sin the central K*K block an denote by image , i=1,2,…, (L-K+1)2-1,the sample vectors corresponding to the other blocks. Letimagerand image be the associated noise less sample vectors of imager and image and, respectively. It can be easily calculated that

image (10)

In above equation we used the fact that noise u is white and uncorrelated with signal. with equation if

image (11)

Where T is a preset threshold, then we select image as a sample vector of Xv. Suppose we select n sample vectors of Xv including the central vector image For the convenience of expression, we denote these sample vectors asimage The noise less counter parts of these vectors are denoted asimageimage accordingly. The training dataset for xv is then formed byIMAGE The noiseless counterpart of Xv is denoted as image

Wavelet

For the wavelet transform, the coefficients at the course level represent a larger time interval but a narrower band of frequencies. This feature of the wavelet transform is very important for image coding. In the active areas, the image data is more localized in the spatial domain, while in the smooth areas, the image data is more localized in the frequency domain. With traditional transform coding, it is very hard to reach a good compromise. The target region (damaged or lost data or object to be removed) information of the image can be divided into two kinds of conditions. The first class, the distribution of the target information of the image is the local and concentration.

Wavelet transform has been used for various image analysis problems due to its nice multi-resolution properties and decoupling characteristics. The proposed algorithm utilizes the advantages of wavelet transforms for image Enhancement. Wavelet transform has been used as a good image representation and analysis tool mainly due to its multi-resolution analysis, data reparability, and compaction and sparsely features in addition to statistical properties [13]. A wavelet function (t) is a small wave, which must be oscillatory in some way to discriminate between different frequencies. The wavelet [5] [6] [10] contains both the analysing shape and the window. In order to observe the degree of influence of image textural on the reconstructed composition, we applied the two-level wavelet transformation to separate an image into three frequency components: high, medium, and low, as shown in Fig.3 (a).

innovations-thoughts-ideas-Block-diagram

Figure 1: Block diagram for enhancement of image

innovations-thoughts-ideas-resulting-images

Figure 2: resulting images

innovations-thoughts-ideas-LPG-Example

Figure 3: LPG Example

innovations-thoughts-ideas-Wavelet-method

Table 1: The PSNR (dB) results of the image enhancement in the Wavelet PCA method.

The original image was processed through a secondary-level wavelet transformation analysis, as illustrated in Fig4 (a), where the highlighted image in the uppermost left hand corner is represented by the section LL2 illustrated in Fig.4 (b). Where analysis is concerned, the components of the overall image composition are all taken into consideration. This procedure can also be utilized as preliminary image analysis. The four components LL2, LH2, LH2, and HH2 are then processed through reversed wavelet transformations to heighten the resolution of the image.

innovations-thoughts-ideas-layers-wavelet

Figure 4: Results of the wavelet transformation analysis derived from various layers of a given image

Discrete Wavelet Transform

Filters are one of the most widely used signal processing functions. Wavelets [7][8][9] can be realized by iteration of filters with rescaling. The resolution of the signal, which is a measure of the amount of detail information in the signal, is determined by the filtering operations, and the scale is determined by up sampling and down sampling (sub sampling) operations [25].

The DWT is computed by successive low pass and high pass filtering of the discrete time-domain signal as shown in figure 5. This is called the Mallat algorithm or Mallat-tree decomposition. Its significance is in the manner it connects the continuous time mutiresolution to discrete-time filters. In the figure, the signal is denoted by the sequence x[n], where n is an integer. The low pass filter is denoted by G0 while the high pass filter is denoted by H0. At each level, the high pass filter produces detail information d[n], while the low pass filter associated with scaling function produces coarse approximations, a[n].

innovations-thoughts-ideas-Three-level-wavelet

Figure 5: Three-level wavelet decomposition tree

innovations-thoughts-ideas-reconstruction-tree

Figure 6: Three-level wavelet reconstruction tree

A highest frequency of ω, which requires a sampling frequency of 2ω radians, then it now, has a highest frequency of ω/2 radians. It can now be sampled at a frequency of ω radians thus discarding half the samples with no loss of information. This decimation by 2 halves the time resolution as the entire signal is now represented by only half the number of samples. Thus, while the half band low pass filtering removes half of the frequencies and thus halves the resolution, the decimation by 2 doubles the scale.

Figure 6 shows the reconstruction of the original signal from the wavelet coefficients. Basically, the reconstruction is the reverse process of decomposition. The approximation and detail coefficients at every level are up sampled by two, passed through the low pass and high pass synthesis filters and then added. This process is continued through the same number of levels as in the decomposition process to obtain and H the original signal. The Mallat algorithm works equally well if the analysis filters, G0 and H0, are exchanged with the synthesis filters, G1 and H1.Then the filters have to satisfy the following two conditions as given in [26]:

IMAGE (12)

image (13)

Implementation

Steps for the image enhancement by wavelet with PCA analysis are carried out as below.

I. Apply algorithm of LPG on input blur image.

II. Do the PCA transform of LPG output image as below.

a. Get some data.

b. Subtract the mean.

c. Calculate the covariance matrix.

d. Calculate the eigenvectors and eigenvalues of the covariance matrix.

e. Choosing components and forming a feature vector.

f. Deriving the new data set.

III. Apply the soft thresholding of PCA output image.

IV. Take Inverse PCA.

V. Do Wavelet transform of IPCA output image.

VI. Apply the soft thresholding of Wavelet output image.

VII. Take Inverse WT.

VIII. Calculate PSNR for the evaluation of the algorithm.

Conclusion

To preserve the local image structures when Enhancement, we modeled a pixel and its nearest neighbors as a vector variable, and the Enhancement of the pixel was converted into the estimation of the variable from its noisy observations. The PCA technique was used for such estimation and the PCA transformation matrix was adaptively trained from the local window of the image. However, in a local window there can have very different structures from the underlying one; therefore, a training sample selection procedure is necessary.

The block matching based local pixel grouping (LPG) was used for such a purpose and it guarantees that only the similar sample blocks to the given one are used in the PCA transform matrix estimation. The PCA transformation coefficients were then shrunk to remove noise.

Refrences

Select your language of interest to view the total content in your interested language
Post your comment

Share This Article

Relevant Topics

Article Usage

  • Total views: 11978
  • [From(publication date):
    May-2012 - Feb 23, 2018]
  • Breakdown by view type
  • HTML page views : 8114
  • PDF downloads : 3864
 

Post your comment

captcha   Reload  Can't read the image? click here to refresh

Peer Reviewed Journals
 
Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals
International Conferences 2018-19
 
Meet Inspiring Speakers and Experts at our 3000+ Global Annual Meetings

Contact Us

Agri & Aquaculture Journals

Dr. Krish

[email protected]

1-702-714-7001Extn: 9040

Biochemistry Journals

Datta A

[email protected]

1-702-714-7001Extn: 9037

Business & Management Journals

Ronald

[email protected]

1-702-714-7001Extn: 9042

Chemistry Journals

Gabriel Shaw

[email protected]

1-702-714-7001Extn: 9040

Clinical Journals

Datta A

[email protected]

1-702-714-7001Extn: 9037

Engineering Journals

James Franklin

[email protected]

1-702-714-7001Extn: 9042

Food & Nutrition Journals

Katie Wilson

[email protected]

1-702-714-7001Extn: 9042

General Science

Andrea Jason

[email protected]

1-702-714-7001Extn: 9043

Genetics & Molecular Biology Journals

Anna Melissa

[email protected]

1-702-714-7001Extn: 9006

Immunology & Microbiology Journals

David Gorantl

[email protected]

1-702-714-7001Extn: 9014

Materials Science Journals

Rachle Green

[email protected]

1-702-714-7001Extn: 9039

Nursing & Health Care Journals

Stephanie Skinner

[email protected]

1-702-714-7001Extn: 9039

Medical Journals

Nimmi Anna

[email protected]

1-702-714-7001Extn: 9038

Neuroscience & Psychology Journals

Nathan T

[email protected]

1-702-714-7001Extn: 9041

Pharmaceutical Sciences Journals

Ann Jose

[email protected]

1-702-714-7001Extn: 9007

Social & Political Science Journals

Steve Harry

[email protected]

1-702-714-7001Extn: 9042

 
© 2008- 2018 OMICS International - Open Access Publisher. Best viewed in Mozilla Firefox | Google Chrome | Above IE 7.0 version