Wavelets Based Simulation and Visualization Approach for Unmixing of Hyperspectral Data

This paper discusses the use of wavelet based simulation and visualization methods to unmix multi-component, intimate


Introduction
Unmixing hyperspectral signals and prole is still a challenge faces environmental, remote sensing, and geoscience communities.Hyperspectral data provide wealth of information about the physicochemical conditions of earth and extraterrestrial targets.However, the purity level of the spectra is a major challenge.Very often the hyperspectral proles contain several end members or components.In theory if a mixed spectra can be simulated using a mathematical approach then an inversion technique would reveal the end members.This approach was tested to simulate mixed spectra of different ratios of this binary system of Na 2 SO 4 -MgSO 4 .On earth these chemicals salts of high enough concentrations can be considered economic minerals and salts that can be used in several industrial applications.For example, Na 2 SO 4 is used in manufacturing wood pulp, glass industry, thermal storage, drying agent.Whereas magnisum sulfate is used commonly as epsom salt for medical purposes, and used in agriculture to increase sulfer or magnesium concentrations in soil, and it is used as brewing salt in beer among other uses.On Mars these salts have been investigated as an analogue to explorer the presence of microbes that possibly inhabit cold MgSO 4 rich brines in Mars.This binary system is also an adequate analogue of the martian salts and its salty regolith.For more details about the endmember extraction and unmixing, we refere to [8][9][10][11][12][13][14][18][19][20][21].
An example from spectral prole extracted from AVRIS image of White Sand in New Mexico, USA is shown in Figure 1.However, the data used in simulation here were obtained using the HR-1024i from the Spectra Vista Corp. (SVC) which is the latest model from their next generation of high performance single-beam eld spectroradiometer measuring over the visible to short-wave infrared wavelength range (350-2500nm).
Reflectance spectroscopy and hyperspectral imaging analysis has picked up because it is reliable, fast, and less expensive and not intrusive [7].However, spectral patterns for mixed materials can't be visually understood.Spectral pretreatment techniques must be applied to smooth spectral graph such as data normalization, continuum removal among other methods.We hypothesis that if the mixed spectra can be simulated mathematically, then we can obtain the end members that make up the mixed spectra.Toward this end, mathematical approaches or models that describe the reectance  The simulation is done by performing three methods, Lagrange interpolation, and Daubechies and Haar wavelet threshold denoising methods, to simulate the given data.The original data can be recovered rst from Lagrange coffiecients and secondly by the (Daubechies and Haar) wavelets coffiecients through so called discrete threshold wavelet lter banks.In order to facilitate the simulation analysis and processing, and as a sample, we considered the reflectances of the wave lenght range 250-2500 nm of the binary system of Na 2 SO 4 -MgSO 4 , then its Lagrange and wavelet threshold transform/interpolation should be discretized.The advantages of the interpolating curves using the methods above is that the entire data can be inferred as if we observe a dense enough set of points and only small amount of coffiecients are needed to achieve accurate approximation.

Metheodolgy
In this study, laboratory experiments under controlled conditions have been carried out to prepare pure Na 2 SO 4 and MgSO 4 crusts and their mixtures.Analytical grade compounds of Na 2 SO 4 and MgSO 4 were considered in this study.The salts were dissolved in water in 1000-ml volumetric asks.Total quantity of 200 ml of each solution were removed and placed in glass Petri dishes.The water was removed from the Petri dishes by evaporation at 40° C in an electric oven for 12 hours, then the samples were removed from the oven.Immediately following the removal from the oven, the samples were placed in a desicator, after which the reectance readings were made.A HR-1024i spectroradiometer was used to obtain the reflectance readings.Most of salts encountered in soils are a mixture of two or more type of salts.Sets of chemically mixed salt samples were prepared from which the pure salts were made using the same steps as described above.The mixing proportions are in 1 : 0, 0:75 : 0:25, 0:5 : 0:5, 0:25 : 0:75, and 0 : 1 weight mixing ratios.The spectra of the mixed samples were compared with the spectra of the pure samples solution fraction.All the pure and mixed samples were examined under petrographic and binocular microscopes for crystal size and morphology observations.Another runs were done by placing soil in the Petri dishes allowing the salt crus to grow on the top.In the ternary system consisting of sodium sulfate and magnesium sulfate, the following solid phases appear in the temperature range from -10° to 110° C: The temperatures and concentration ranges at which these solids appear are shown in the Figure 2. The phase diagrams shown on these pages are calculated with the Extended UNIQUAC thermodynamic model (http://www.phasediagram.dk/default.htm).The equilibrium lines and the experimental data in the diagram represent compositions and temperatures at which two solid phases are in equilibrium with the same solution.The water content of the solutions are not shown.
It can be thought of as a third dimension in the diagram.The green lines in the gure above are tie lines indicating phases and compositions in equilibrium with each other.
The study will apply Lagrange interpolation and wavelet coffiecient analysis.The Lagrange interpolating polynomial is the polynomial P(x) of degree ≤ n-1 that passes through the n points (x k , f(x k )), for k=1,...,n and is dened by (1) where (2) When constructing interpolating polynomials, there is a tradeoff between having a better fit and having a smooth well-behaved fitting function.The more data points that are used in the interpolation, the higher the degree of the resulting polynomial, and therefore, a highdegree interpolation may be a poor predictor of the function between points, although the accuracy at the data points will be "perfect".
Daubechies was first to construct compactly supported orthogonal wavelets with a preassigned degree of smoothness [2].She intended to construct a wavelet with N vanishing moments and supported in [0; 2N-1]; where a wavelet function ψ is said to have N (≥ 2) vanishing moments if, When N = 1, it is Haar wavelet and has only one vanishing moment.In fact, all systems built by using the unitary extension principle (UEP) of Ron and Shen [17] have only one vanishing moment.

By choosing
, a scaling function ϕ has support [0, 2] : Daubechies wavelet does not have a closed form, but instead, can be obtained recursively by (6) (7) It is known that the smoothness of the wavelets increases with N. For an application in numerical analysis, Coifman asked Daubechies to construct a family of wavelets ψ having N vanishing moments, minimum size support and (8) in such a way we have a smooth orthogonal system [2].Daubechies wavelet has vanishing moment for wavelet function ψ only.However, Daubechies designed, at that request, a wavelet (Coiets) that has vanishing moment for both wavelet and scaling functions ψ and ϕ.The wavelet is near symmetric such that wavelet function ψ has 2N vanishing moments and the scaling function has 2N-1 vanishing moments.The two functions have a support width of length 6N-1.Also there is no closed form for Coiets and it can be obtained recursively.

Denition [3]
A compactly supported function is said to be renable if (9) for some nite supported sequence .The sequence h 0 is called the low pass filter of ϕ.

Let
be of the form (10) where is a nitely supported sequence and is called a high pass filter of the system.Dene the wavelet system (11) where . The general setup is to construct a wavelet for and of the form of X(Ψ) [4].
The Haar function is dened by Note that the Haar wavelet can be written as, ( Already in 1910 it was proved by Haar that the functions constitute an orthonormal basis for .The Daubechies wavelets functions for different order are given in Figure 3.

Lagrange and Wavelet Simulations
To illustrate the simulation, assume that X(Ψ) forms a wavelet system for .For any signal/data , f can be expanded as follows, (17) Our simulation is based on signal classication approach using discrete wavelet transform (DWT).A DWT-based linear (un) mixing system is designed specially for estimation.Using the wavelet transform, the original signal is represented by a set of wavelet transform coffiecients, and features are extracted from these coffiecients.Wavelet theory provides a good technique for signal approximation using scaled basis functions.The inverted data of the hyperspectral imaging illustrated in Figure 4.Note that the XY-axes ranges presented in Figures 4-11 are not refer to the actual wavelength and reflectance ranges of this study.The ranges presented here came from the ListLinePlot command using Mathematica software that plots a line through the points {1, y 1 }, {2, y 2 },..., {647, y 647 } where the set {y 1 , y 2 ,..., y 647 } presenting the spectral image data of the binary system of Na 2 SO 4 -MgSO 4 .
The first approach utilized the Lagrange Interpolation Technique.The simulation using Lagrange interpolation is illustrated in Figure 5.A closer view of the simulation using this approach is depicted in Figure 6.
The second approach is by using wavelets interpolation, the first will be done using Haar wavelets and the second by Daubechies wavelets.The simulations using these functions are depicted in Figure 7 and 8.

Splitting the spectra using vector projection
For this section, we represent the hyperspectral data as a sum of two functions that makes it.One way doing that is to use the vector projection method.The vector projection method of a vector a on (or onto) a nonzero vector b is the orthogonal projection of a onto a straight line parallel to .It is a vector parallel to , defined as , , ,      To illustrate the idea, we choose a specic vector such that where and .
Figure 10 shows the graphs of these splitted simulations.

Conclusion
The study used two approaches 1) semantic one by tracing the diagnostic spectral features such as the location, shape and depth of the band for certain endmembers and through study the association of the diagnostic feature between the mixed spectra and their endmembers (semantic or diagnostic approach), 2) mathematically through wavelet approach that simulate and unmix the reflectance spectra.To this end, a one-dimensional threshold method for hyperspectral images has been proposed based on Lagrange interpolation and discreet wavelet threshold denoising method.The performance of Lagrange interpolation method and the wavelet approach based on Haar and Daubechies threshold methods to simulate Na 2 SO 4 -MgSO 4 binary system yielded satisfactory results.Vector projection method was used successfully to split the hyperspectral data as a sum of two data functions that simulate it.

Figure 6 :
Figure 6: A closer view of the simulation using Lagrange interpolation.

Figure 8 :Page 6 of 8 Figure 9 :
Figure 8: A closer view for the simulations using Haar and Daubechies wavelets thresholds, respectively.

Figure 10 :
Figure 10: Splitting the the binary system of Na 2 SO 4 -MgSO 4 using the vector projection method.