can see in Figure 4(c), after applying the filter, noise is removed and the layers are more visible, however, the many of the boundaries have become murky. Structures with low amplitude are most affected by denoising, as opposed to the high-amplitude structures which are more visible, and therefore, reconstruction of the layers has not been done efficiently.

6.3. Applying the AR Filter in the Wavelet Domain

Here we have used the undecimated discrete wavelet transform in applying the autoregressive filter to the noise section in the wavelet domain. We note that here, the applied filter is linear which is a great advantage in noise reduction due to the linearity of the wavelet space.

The denoised section by using this method is shown in Figure 4(d) which looks promising, since the layer boundaries are distinguishable and there is logical smooth trend throughout the section. Both the low and high amplitude structures are more well-defined compared to the autoregressive filter in the f-x domain. We also note that the layer trends are efficiently reconstructed.

6.4. Applying the Method on Real Data

Here, we apply the linear regressive model to real data first in the f-x and then in the wavelet domain in order to remove random noise and eventually compare the results. We show traces #450 onward for better comparison in Figure 6.

We first, suppress the noise in the f-x domain using the autoregressive filter on the noisy real data shown in Figure 5(a) where the boundaries and traces (especially after trace #450) are difficult to distinguish and the bulges are faded in the noise. As shown in Figure 6(c), the autoregressive filter in the f-x domain has properly reduced the noise and the layer trends are visible.

As mentioned before, the noise is random and does not have a specific source. As we can see in Figure 5(a), the right (and specially lower right) portion of the section is very noisy as noise has covered all the trends and layers. Figure 5(c) shows the denoised section after applying the autoregressed filter in the wavelet domain. As we can see in Figure 5(c), the trends are well-defined, especially in the lower right corner of the section and the previously vague parts can now be distinguished to a much greater extent. Therefore the method has successfully reduced the noise.

6.5. Comparing the Application of AR in the Wavelet and f-x Domains

Considering the previous sections on wavelet and f-x domains, as well as the better performance of the filter in the wavelet domain, here we compare denoising on synthetic and real data in both domains. By comparing Figure 5(b) and Figure 5(c), we notice the better performance of the filter in denoising the synthetic data in the wavelet domain. Especially, in the later arrival times which correspond to greater depths, noise reduction is more evident. We note that the boundaries are better represented and are more distinguishable. Also, by comparing

Figure 5. The GPR data section; (a) raw data; (b) after applying autoregression in the f-x domain; (c) after applying autoregression in the wavelet domain.

Figure 6. The GPR data section from trace 450 to the end; (a) raw data; (b) after applying autoregression in the f-x domain; (c) after applying autoregression in the wavelet domain.

Figure 6(b) and Figure 6(c), we can see that again the wavelet domain has had a much better performance in reducing the noise and retrieving the real signal (the layer trends are more visible). Overall, the autoregressive filter in the wavelet domain had done a better job in retrieving the signal due to the linearity of the wavelet domain.

7. Conclusions

As discussed above, various factors such as phone networks, power posts, utility poles, etc. cause contamination in the GPR data. Since the goal of this study was to increase the signal to noise ratio, a method was chosen to damage the signal as least as possible while reducing the noise. In this study, the noise suppression procedure was applied to both synthetic and real GPR data in f-x and wavelet domains through using autoregressive filter. As we see, noise reduction improves interpretation of data and Autoregressive filter bears good results in both f-x and wavelet domains. Which means that Linear regression in the wavelet domain leads to better results, compared to those of the f-x domain, due to the local nature of the wavelet transform and the imposed linearity on the events on different scales.

We should also note that in contrast with the f-x domain, the autoregressive filter is linear just as the wavelet domain which is why the filter does not work as well in the f-x domain.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] Annan, A.P. (2001) Ground Penetrating Radar. Workshop Notes, No. September, 192.
[2] Daniels, D.J. (2004) Ground Penetrating Radar.
[3] Burrus, S.C. and Gopinath, R.A. (2013) Wavelets and Wavelet Transforms. Rice Univ. Houston, Texas, 1-281.
[4] Leonard, M. and Kennatt, B.L. (2008) Multicomponent Autoregressive Techniques for the Analysis of Seismograms. Physics of the Earth and Planetary Interiors, 113.
[5] Oskooi, B., Julayusefi, M. and Goudarzi, A. (2015) GPR Noise Reduction Based on Wavelet Thresholdings. Arabian Journal of Geosciences, 2937-2951.
[6] Kim, J.H. and Cho, S.J. (2007) Removal of Rining Noise in GPR Data by Signal Processing. Geosciences Journal, 11, 75-81.
[7] Canales, L.L. (1984) Random Noise Reduction. 54th Ann. Mtg. Society of Exploration Geophysicists, 3.
[8] Yilmaz, O. (2000) Seismic Data Analysis: Processing, Inversion, and Interpretation of Seismic Data. Soc. Explor. Geophys, 2000, 1, 7. Society of Exploration Geophysicists, 1-7.
[9] Hashemi Amroabadi, M. (2010) Data Denoising in Analog & Digital Domain. A Thesis Present to Ryerson Univ., 11.
[10] Harrison, M.P. (n.d.) f-x Linear Prediction Filtering of Seismic Images. 147-165.
[11] Luetkepohl, H. (2011) Vector Autoregressive Models, Economics Departments, In-stitutes and Research Centers in the World. Retrieved from IDEAS.
[12] Fuss, R. (2007) Vector Autoregressive Models. Albert-Ludwigs-University of Freiburg.
[13] Naghizadeh, M. and Sacchi, M. (2012) Multicomponent f-x Seismic Random Noise Attenuation via Vector Autoregressive Operators. Geophysics, 99, 15.
[14] Wang, Y. (2015) Frequencies of the Ricker Wavelet. Geophysics, 80, A31-A37.

comments powered by Disqus
JSIP Subscription
E-Mail Alert
JSIP Most popular papers
Publication Ethics & OA Statement
Frequently Asked Questions
Recommend to Peers
Recommend to Library
Contact Us

Copyright © 2020 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.