Improvement of a Low Cost MEMS Inertial-GPS Integrated System Using Wavelet Denoising Techniques
Improvement of a Low Cost MEMS Inertial-GPS Integrated System Using Wavelet Denoising Techniques
International Journal of Aeronautical and Space Sciences. 2011. Dec, 12(4): 371-378
Copyright ©2011, The Korean Society for Aeronautical Space Science
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( 3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : November 03, 2011
  • Accepted : December 07, 2011
  • Published : December 30, 2011
Export by style
Cited by
About the Authors
Chang Ho Kang
Sun Young Kim
Chan Gook Park

In this paper, the wavelet denoising techniques using thresholding method are applied to the low cost micro electromechanical system (MEMS)-global positioning system(GPS) integrated system. This was done to improve the navigation performance. The low cost MEMS signals can be distorted with conventional pre-filtering method such as low-pass filtering method. However, wavelet denoising techniques using thresholding method do not distort the rapidly-changing signals. They can reduce the signal noise. This paper verified the improvement of the navigation performance compared to the conventional pre-filtering by simulation and experiment.
1. Introduction
By the recent development of the micro electromechanical system (MEMS) technology, many applications of low cost MEMS-global positioning system (GPS) integrated navigation system are being popularly researched. These MEMS-based inertial sensors have been integrated with the GPS to provide reliable positioning solutions in case of GPS outages that commonly occur in the urban areas. Especially, low cost MEMS-GPS integrated navigation system is used for mobile robots, unmanned aerial vehicle or micro-aerial vehicle and pedestrian navigation system (Hasan et al., 2010; Noureldin et al., 2004; Yoon and Vaidyanathan, 2004).
In the MEMS-GPS integrated system, inertial sensor data includes large signal noise which makes significant error on the position result. Sensor noise is well compensated on a good observable trajectory in general integrated navigation system. But it cannot be exactly compensated on low cost systems. If the noise component could be removed, the overall inertial navigation accuracy is expected to improve considerably. The resulting position errors are proportional to the existing sensor noise.
In this paper, the wavelet denoising technique is implemented to eliminate the sensor noise for improving the accuracy of the navigation results (Hasan et al., 2010; Kang and Park, 2009; Nassar and El-Sheimy, 2005; Noureldin et al., 2004). Wavelet transforms have been successfully applied for denoising, classification, recognition, compression, and other applications. It can decompose the signal to a frequency component in local time. By using this characteristic, the wavelet denoising method shrinks the signal noise by eliminating the frequency component which contains only noise. Furthermore, Denoising quality is also improved by wavelet thresholding techniques. These are the schemes which are used to remove the noise in wavelet transform domain by using thresholding operator. Weaver and later DeVore proposed a technique known as the soft threshold. Donoho and several others coupled this result, known as wavelet shrinkage, with threshold selection techniques. This resulted in the methods such as, WaveShrink, SureShrink, and VisuShrink (Donoho 1995, Donoho and Johnstone, 1994, 1995). In this paper, the wavelet thresholding method is verified using the collected real data from the field test. It is implemented to the MEMS/GPS integrated system. This will be shown in the conclusion.
This paper is organized as follows: Section 2 describes the basic information of the wavelet denoising technique. Section 3 presents the wavelet denoising method implemented to the MEMS/GPS integrated system. It contains the structure of the MEMS/GPS integrated system and Kalman filter which are used in the integration of the system. The experimental result is analyzed in Section 4. Finally, conclusion is given in Section 5.
2. Wavelet Denoising Technique
- 2.1 Discrete wavelet transform
Wavelet transform is a signal transform technique which is popularly used in several areas such as image processing and audio signal processing. The comparison of the wavelet transform and Fourier transform is shown in Fig. 1 . As specified in the left of Fig. 1 , the Fourier transform decomposes the signal into each frequency component over the entire time interval. It means that the time domain information is lost in transforming to the frequency domain. When looking at a Fourier transform of a signal, it is impossible to recognize at what point a particular event takes place. To overcome this drawback, the transform is adapted to analyze a window of the signal at a particular instance of time. Furthermore, it is necessary to have multiple resolutions in time and frequency domain in order not to tradeoff corresponding to the choice of the window function’s width. The wavelet analysis is based
Lager Image
Various types of time-frequency domain sampling. Fourier, Wavelets.
on a windowing technique with variable-sized windows as shown in Fig. 1 . The wavelet transform applies the wide window (long time intervals) to low frequency and the narrow window (short time intervals) to high frequency.
Discrete wavelet transform is executed as in Eq. (1).
Lager Image
Where, d is an n ×1 vector comprising of both the discrete scaling coefficients, and discrete wavelet coefficients. W is an orthogonal n × n matrix associated with the orthonormal wavelet basis chosen. x is a vector of function values, x = [ x ( t 1 ), x ( t 2 ), …, x ( t n )] T at equally spaced points ti .
The scaling coefficient can be expressed as the convolution of the scaling coefficient and the wavelet coefficient in the lower level. It means that the low-frequency area can be decomposed to the high-frequency area and low-frequency area. The relationship is shown as in Fig. 2 .
Wavelet transform is progresses stage by stage. When the signal is divided into low-frequency waves, it requires twice the amount of data. In addition, the lowest possible decomposable frequency area matches the DC value of the Fourier transform. This is calculated by using the entire data (Goswami and Chan, 1999).
- 2.2 Wavelet thresholding technique
Wavelet denoising technique was developed to estimate a signal which is corrupted by additive noise. Wavelet thresholding technique is a signal estimation technique that exploits the capabilities of a wavelet transform for signal denoising. It removes the noise, by eliminating coefficients that are insignificant compared to some threshold. Wavelet thresholding technique assumes that the magnitude of the actual signal is greater than the noise level, and the noise is white noise.
General low-pass filter (LPF) has the characteristics of
Lager Image
Wavelet decomposition.
removing all the frequencies over a certain threshold. If sensor data exists in the high-frequency component, low pass filter results in the loss of the sensor signal. For example, the accelerometer signal has sudden change when the vehicle accelerates. The LPF distorts this change of signal. It turns it into a gradually changing signal. It is clearly shown in Fig. 3 . In Fig. 3 , the actual signal is expressed with a solid line. The result of the low-pass signal is expressed with the 1-dot chain line, which is found to smoothly follow the signal change. For this reason, the LPF cannot be used as the pre-processing filter in the inertial navigation system.
However, the wavelet thresholding technique reduces the noise level with almost no distortion for the sudden change in signal, as shown in Fig. 3 . The result of the thresholding technique has almost no distortion and is accurate. So that it sits right on the actual signal almost indistinguishably. Therefore, it can be used by the pre-processing filter for the inertial sensor signal and overcomes the shortages of the existing LPF.
There are some kinds of thresholding methods namely, hard thresholding scheme, soft thresholding scheme (Antoniadis, 2007), nonnegative garrote function scheme (Gao, 1998) and polynomial thresholding scheme (Smith et al., 2008), which are illustrated in Figure 4 . The figure shows the results of the thresholding schemes when the input signal u is, x ∈ [?10 10]and threshold λ is 5. Thresholding allows the signal itself to decide which wavelet coefficients are significant.
The soft thresholding scheme generally results in systematically biased estimates. The hard thresholding estimates are less biased but denoising performance is worse. Equation (2) shows the hard thresholding function and Eq. (3) shows the soft thresholding function.
Lager Image
Comparison of the signal using wavelet and low-pass filter.
Lager Image
Lager Image
Where, u is the wavelet coefficient, and λ is threshold
To overcome the drawbacks and to achieve better performance, compromise thresholding schemes have been suggested in previous studies. As an example, nonnegative garrote function is mentioned. The resulting wavelet thresholding scheme off ers generally smaller mean squared error (MSE), less sensitive perturbations in the data, smaller bias and overall MSE. Equation (4) shows the nonnegative garrote function.
Lager Image
These schemes are functions dependent on the statistical models which give zero for small values of u and close values to u for larger values of the argument.
Polynomial thresholding scheme is also one of the solutions to remedy the drawbacks of both hard and soft thresholding scheme. The threshold operator used in this scheme is in the form of polynomials of u , which can be written in Eq. (5).
Lager Image
Where a = [ a 0 , a 1 , … , an ] T is an array of the coefficients for the polynomial. The flexibility of this scheme allows many possibilities for the removal of some coefficients while
Lager Image
Wavelet thresholding methods.
preserving others.
Among these methods, polynomial thresholding scheme has better performance of reconstruction of the signal than other schemes. So, the polynomial thresholding scheme is used for this experiment. The performance of reconstruction of the signal using thresholding schemes is shown in Fig. 5 . The MSE of reconstruction is presented in Table 1 . Doppler signal varying in its center frequency is used in the performance test. Figure 5 and Table 1 indicate the simulation results which demonstrate the superiority of the polynomial thresholding approach over the other schemes.
One of the important elements influencing the performance in thresholding technique is how the standard value of λ is set (Chan and Peng, 2003; Hasan et al., 2010). Generally, it is determined by Eq. (6). Here, σ is the standard deviation of the signal, and n is the number of signal samples.
Lager Image
a) The results of the signal reconstruction. LPF: low-pass filter.
Lager Image
b) Enlargement of the reconstruction results. LPF: low-pass filter.
Lager Image
The value determined by Eq. (6), however, does not lead to the optimal result. So, the appropriate value must be determined through experimentation. The thresholding algorithm is executed as follows. First, the wavelet transform is executed for the signal to acquire the wavelet coefficients. Next, the thresholding operation is executed for each wavelet coefficient. Then, the original coefficients are replaced to the coefficients from the result of the thresholding operation. Finally, inverse transform is carried out (Chan and Peng, 2003; Yoon and Vaidyanathan, 2004).
3. MEMS Inertial-GPS Integrated Navigation System using Wavelet technique
The MEMS-GPS integrated navigation system is used in the experiment. It is the general 15th order loosely coupled model (Titterton et al., 1997). The construction of the system is explained in Fig. 6 .
The state variables of the Kalman filter consist of the position error, speed error, attitude error, accelerometer error, and gyro error. The measurement of the Kalman filter is the position information of the GPS which updates every second. The inertial sensor signal using wavelet thresholding technique is updated for every 0.01 second. It is used to calculate the integrated navigation result. If the LPF is applied
MSE of the signal reconstructionMSE: mean square error, LPF: low-pass filter.
Lager Image
MSE of the signal reconstruction MSE: mean square error, LPF: low-pass filter.
Lager Image
Construction of micro electromechanical system inertial-global positioning system (GPS) integrated system. INS: inertial navigation system.
to the GPS signal, it will yield better performance. However, the GPS signal was not filtered separately, as the objective of this research is to find out how much performance is improved by the denoising of the inertial sensor signal.
In this paper, the Daubechies wavelet (Daubechies, 1992) is used for wavelet transform as shown in Fig. 7 . The setting of the wavelet does not have a big difference no matter which one is used. MEMS inertial sensor signal includes the actual vehicle motion dynamics and the sensor noise as well as some other undesirable noise such as the vehicle engine vibration. Therefore, the criterion for the selection of the appropriate wavelet level of decomposition (LOD) will be different from the stable motion case. Furthermore, the choice of threshold is crucial to maintain the quality of the denoising process. This should be made carefully (Hasan et al., 2010).
The wavelet transform was carried out to the 3th order. The threshold value is selected from Eq. (6) for the accelerometer and the gyro. These values are determined experimentally through multiple simulations.
The thresholding algorithm is implemented to the inertial
Lager Image
Daubechies scaling function and wavelet function.
Lager Image
Flowchart of the wavelet thresholding algorithm. MEMS: micro electromechanical system, GPS: global positioning system.
sensor as follows. The wavelet transform is executed for the gyro. The accelerometer signals of each axis and the wavelet coefficients are acquired. Next, the wavelet coefficients are replaced by the coefficients resulting from the thresholding operation. The inverse transform is carried out. The inverse transformed signal is the noise-mitigated signal, which is used to execute the MEMS inertial-GPS integrated navigation. The process is described briefly in the flowchart shown in Fig. 8 .
4. Experimental Result
The trajectory of the experiment is the belt way of the Seoul National University campus shown in the Fig. 9 . It also shows the position data of the MEMS-GPS integrated system. The true trajectory is measured by differential GPS for twenty minutes of travelling. In this test, a low grade MEMS inertial
Lager Image
The trajectory of the experiment.
Typical performance of the MTi-GRoll/Pitch and heading accuracy means dynamic accuracy.GPS: global positioning system, INS: inertial navigation system, RMS: root mean square.
Lager Image
Typical performance of the MTi-G Roll/Pitch and heading accuracy means dynamic accuracy. GPS: global positioning system, INS: inertial navigation system, RMS: root mean square.
Lager Image
a) Accelerometer signal. LPF: low-pass filter.
Lager Image
b) Gyroscope signal. LPF: low-pass filter.
Lager Image
Enlargement of the inertial navigation system data. MEMS: micro electromechanical system, LPF: low-pass filter.
Lager Image
a) The position data.
Lager Image
b) Micro electromechanical system inertial-global positioning system navigation error.
RMSE of the MEMS-GPS integrated systemRMSE: root mean square error, MEMS: micro electromechanical system, GPS: global positioning system, LPF: low-pass filter.
Lager Image
RMSE of the MEMS-GPS integrated system RMSE: root mean square error, MEMS: micro electromechanical system, GPS: global positioning system, LPF: low-pass filter.
navigation system (INS) (MTi-G) and GPS integrated system are used. The sensor specifications used in the test are shown in Table 2 . The minimum number of available satellites was 8. The average vehicle speed was 40 km/h. All the analysis results in this paper are implemented using the MATLAB computer-aided design software including the wavelet analysis.
Figure 10 illustrates the experimental result of the MEMS output signal arranged x, y, and z axis in order. It shows original MEMS sensor signal. It also shows the signals applied to the filtering methods. Comparison was made between the wavelet threshold method and the noise reduction method used in the research of Sameh and Naser (2005). In the research of Sameh and Naser (2005), the method of controlling LOD eliminates the noise of the INS signal. Its results are almost identical to that of the LPF, as the frequency response of the scale function is similar to the LPF. Thus, this result is marked as LPF in Figs. 10 - 12 .
Figure 10 shows that result of the comparison between the wavelet thresholding technique and LPF. The original signal shown in cyan line means the raw data without any de noising technique. The signal applied to wavelet thresholding technique and LPF is shown as a red line and blue line, respectively. LPF means the signal applied to wavelet transform without thresholding technique.
In the stable state, both Wavelet thresholding technique and LPF do not distort the sensor data. On the other hand, LPF is more efficient than the Wavelet thresholding technique in terms of signal denoising.
In the case of INS signals, passing through the low pass filter often cause the distortion of sensor signal as shown in Fig. 11 . This shows the data influenced by the dynamic motion of the vehicle for a period of one second. On the other hand, the signal applied to Wavelet thresholding technique maintains good performance during this period.
Consequently, using the pre-filtering method such as the low pass filter to remove the noise of the MEMS signal, it results in the decreased navigation performance.
Figure 12 shows the MEMS-GPS Navigation position output data. The navigation error is arranged latitude, longitude, and altitude in order. The data in the altitude position component applied to LPF is verified so that the LPF distorts the sudden change signal.
Table 3 shows the root mean square error (RMSE) of the MEMS-GPS integrated system’s position output data based on DGPS position data. It compares the RMSE between the wavelet thresholding method and the LPF in the direction of north, east and down. The RMSE of wavelet thresholding technique is smaller than LPF. This verifies whether the wavelet thresholding method has better performance than LPF or not.
Therefore, the chosen wavelet filter (Daubechies wavelet) and thresholding technique contributed to eliminate the undesired noise in the system and maintain better performance than the result used by LPF.
5. Conclusions
The LPF used in the MEMS-GPS integrated navigation system has limitations such as signal distortion and delay. Such issues cause the deterioration of the navigation performance. However, the wavelet denoising techniques which use the thresholding method can reduce the signal distortion. This technique is useful for denoising MEMS INS signal which rapidly changes according to the motion of the vehicle.
By applying the wavelet denoising techniques using polynomial thresholding scheme in this paper, it was proved that the overall navigation performance was enhanced by 19.72% more than the performance using conventional prefiltering technique.
This work has been supported by National GNSS Research Center program of Defense Acquisition Program Administration and Agency for Defense Development.
Antoniadis A. (2007) Wavelet methods in statistics: some recent developments and their applications Statistic Surveys 1 16 - 55    DOI : 10.1214/07-SS014
Chan A. K. , Peng C. (2003) Wavelets for Sensing Technologies Artech House Boston
Daubechies I. (1992) Ten Lectures on Wavelets Society for Industrial and Applied Mathematice Philadelphia
Donoho D. L. (1995) De-noising by soft-thresholding IEEE Transactions on Information Theory 41 613 - 627    DOI : 10.1109/18.382009
Donoho D. L. , Johnstone J. M. (1994) Ideal spatial adaptation by wavelet shrinkage Biometrika 81 425 - 455    DOI : 10.1093/biomet/81.3.425
Donoho D. L. , Johnstone I. M. (1995) Adapting to unknown smoothness via wavelet shrinkage Journal of the American Statistical Association 90 1200 - 1224    DOI : 10.2307/2291512
Gao H. (1998) Wavelet shrinkage denoising using the nonnegative garrote Journal of Computational and Graphical Statistics 7 469 - 488    DOI : 10.2307/1390677
Goswami J. C. , Chan A. K. (1999) Fundamentals of Wavelets: Theory, Algorithms, and Applications Wiley New York
Hasan A. M. , Samsudin K. , Ramli A. R. , Azmir R. S. (2010) Comparative study on wavelet filter and thresholding selection for GPS/INS data fusion International Journal of Wavelets, Multiresolution and Information Processing 8 457 - 473    DOI : 10.1142/S0219691310003572
Kang C. W. , Park C. G. (2009) Improvement of INSGPS integrated navigation system using wavelet thresholding Journal of the Korean Society for Aeronautical and Space Sciences 37 767 - 773    DOI : 10.5139/JKSAS.2009.37.8.767
Nassar S. , El-Sheimy N. (2005) Wavelet analysis for improving INS and INS/DGPS navigation accuracy Journal of Navigation 58 119 - 134    DOI : 10.1017/S0373463304003005
Noureldin A. , Osman A. , El-Sheimy N. (2004) A neuro-wavelet method for multi-sensor system integration for vehicular navigation Measurement Science and Technology 15 404 - 412    DOI : 10.1088/0957-0233/15/2/013
Titterton D. H. , Weston J. L. Institution of Electrical Engineers (1997) Strapdown Inertial Navigation Technology Peter Peregrinis Ltd. on behalf of the Institution of Electrical Engineers London, UK
Yoon B. J. , Vaidyanathan P. P. (2004) Wavelet-based denoising by customized thresholding Montreal, Canada Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing II925 - II928