Usually, in ground-penetrating radar (GPR) datasets the user defines the limits between the useful signal and the noise through standard filtering to isolate the effective signal as much as possible. However, there are true reflections that mask the coherent reflectors that can be considered noise. In archaeological sites these clutter reflections are caused by scattering with origin in subsurface elements (e.g., isolated masonry, ceramic objects and archaeological collapses). Its elimination is difficult because the wavelet parameters similar to coherent reflections and there is a risk of creating artifacts. In this study a procedure to filtering the clutter reflection noise (CRN) from GPR datasets is presented. The CRN filter is a singular value decomposition-based method (SVD), applied in the 2D spectral domain. This CRN filtering was tested in a dataset obtained from a controlled laboratory environment, to establish a mathematical control of this algorithm. Also, it has been applied in a 3D-GPR dataset acquired in the Roman villa of Horta da Torre (Fronteira, Portugal), which is an uncontrolled environment. The results show an increase in the quality of archaeological-GPR planimetry that were verified via archaeological excavation.

INT-FFT algorithm presented in this work uses an interpolation methodology to densify 3D-GPR datasets to sharpen images obtained in GPR surveys obtained in an archaeological context. It allows the reconstruction of missing data from the combined use of mathematical transforms (e.g., the Fourier and Curvelet transform) and predictive filters. This technique makes it possible to calculate the missing signal simply by meeting two requirements: the data in the frequency domain must be limited in a range of values and must be able to be represented by a distribution of Fourier coefficients (verified conditions). The INT-FFT algorithm uses an open-access routine (Suinterp, Seismic Unix) to interpolate the GPR profiles based on seismic trace interpolation. This process uses automatic event identification routines by calculating spatial derivatives to identify discontinuities in space by detecting very subtle changes in the signal, thus allowing for more efficient interpolation without artifacts or signal deterioration. We successfully tested the approach using GPR datasets from the Roman Villa of Horta da Torre (Fronteira, Portugal). The results show an increase in the geometric sharpness of the GPR reflectors and have not produced any numerical artifacts. The tests performed to apply the methodology to GPR-3D data allowed for assessing the interpolation efficiency, the level of recovery of missing data, and the level of information lost when one chooses to increase the distance between profiles in the acquisition stage of the data.