On the introduction of canny operator in an advanced imaging algorithm for real-time detection of hyperbolas in ground-penetrating radar data

This paper focuses on the use of the Canny edge detector as the first step of an advanced imaging algorithm for automated detection of hyperbolic reflections in ground-penetrating radar (GPR) data. Since the imaging algorithm aims to work in real time; particular attention is paid to its computational efficiency. Various alternative criteria are designed and examined, to fasten the procedure by eliminating unnecessary edge pixels from Canny-processed data, before such data go through the subsequent steps of the detection algorithm. The effectiveness and reliability of the proposed methodology are tested on a wide set of synthetic and experimental radargrams with promising results. The finite-difference time-domain simulator gprMax is used to generate synthetic radargrams for the tests, while the real radargrams come from GPR surveys carried out by the authors in urban areas. The imaging algorithm is implemented in MATLAB

Measurement of snow water equivalent using drone-mounted ultra-wide-band radar

The use of unmanned aerial vehicle (UAV)-mounted radar for obtaining snowpack parameters has seen considerable advances over recent years. However, a robust method of snow density estimation still needs further development. The objective of this work is to develop a method to reliably and remotely estimate snow water equivalent (SWE) using UAV-mounted radar and to perform initial field experiments. In this paper, we present an improved scheme for measuring SWE using ultra-wide-band (UWB) (0.7 to 4.5 GHz) pseudo-noise radar on a moving UAV, which is based on airborne snow depth and density measurements from the same platform. The scheme involves autofocusing procedures with the frequency–wavenumber (F–K) migration algorithm combined with the Dix equation for layered media in addition to altitude correction of the flying platform. Initial results from field experiments show high repeatability (R > 0.92) for depth measurements up to 5.5 m, and good agreement with Monte Carlo simulations for the statistical spread of snow density estimates with standard deviation of 0.108 g/cm3. This paper also outlines needed system improvements to increase the accuracy of a snow density estimator based on an F–K migration technique

Time-Lapse Monitoring of Two-Dimensional Non-Uniform Unsaturated Flow Processes Using Ground Penetrating Radar

Unsaturated flow in the vadose zone often manifests as preferential flow resulting in transport of water and solutes through the soil much faster than would occur for uniform matrix flow. Time-lapse ground penetrating radar (TLGPR) shows promise as a non-invasive means to monitor unsaturated flow and here is used to monitor lab-scale forced infiltration events for capturing evidence of non-uniform and preferential flow phenomena directly from arrivals in the GPR images while simultaneously characterizing parameters of the flow system, such as bulk water content and rates of wetting front movement. This was accomplished by 1) directly interpreting transient arrivals in GPR profiles for evidence of ono-uniform flow and 2) with the aid of migration processing techniques to improve the quality of GPR images for identification and tracking of transient arrivals related to wetting in the soil. A novel method is described and evaluated to characterize the 2D velocity structure of a soil and used to migrate the GPR images. This method incorporates multi-offset measurements to characterize the depth to a potentially unknown static reflector and root mean square (RMS) velocity above the reflector with incremental changes in travel time to the static reflector and a transient reflector (i.e. the wetting front) determined from single-offset constant offset profiles to determine incremental changes in velocity above and below the transient arrival. The method is applied to TLGPR data during infiltration experiments in a 60 cm deep sand-filled tank and monitored with water content probes. To verify the approach the methodology is applied to GPR data simulated using transient water contents generated by the unsaturated flow simulator HYDRUS 2D given lab-measured hydraulic properties of the soil. For both the empirical and simulated data, we found that the 2D velocity analysis was effective in monitoring changes in the wetting front and that migration of the reflection profiles was able to improve the interpretation of non-uniform flow

SPOT-GPR: a freeware tool for target detection and localizationin GPR data developed within the COST action TU1208

SPOT-GPR (release 1.0) is a new freeware tool implementing an innovative Sub-Array Processing method, for the analysis of Ground-Penetrating Radar (GPR) data with the main purposes of detecting and localizing targets. The software is implemented in Matlab, it has a graphical user interface and a short manual. This work is the outcome of a series of three Short-Term Scientific Missions (STSMs) funded by European COoperation in Science and Technology (COST) and carried out in the framework of the COST Action TU1208 “Civil Engineering Applications of Ground Penetrating Radar” (www.GPRadar.eu). The input of the software is a GPR radargram (B-scan). The radargram is partitioned in subradargrams, composed of a few traces (A-scans) each. The multi-frequency information enclosed in each trace is exploited and a set of dominant Directions of Arrival (DoA) of the electromagnetic field is calculated for each sub-radargram. The estimated angles are triangulated, obtaining a pattern of crossings that are condensed around target locations. Such pattern is filtered, in order to remove a noisy background of unwanted crossings, and is then processed by applying a statistical procedure. Finally, the targets are detected and their positions are predicted. For DoA estimation, the MUltiple SIgnal Classification (MUSIC) algorithm is employed, in combination with the matched filter technique. To the best of our knowledge, this is the first time the matched filter technique is used for the processing of GPR data. The software has been tested on GPR synthetic radargrams, calculated by using the finite-difference timedomain simulator gprMax, with very good results

High-resolution coherency functionals for improving the velocity analysis of ground-penetrating radar data

We aim at verifying whether the use of high-resolution coherency functionals could improve the signal-to-noise ratio (S/N) of Ground-Penetrating Radar data by introducing a variable and precisely picked velocity field in the migration process. After carrying out tests on synthetic data to schematically simulate the problem, assessing the types of functionals most suitable for GPR data analysis, we estimated a varying velocity field relative to a real dataset. This dataset was acquired in an archaeological area where an excavation after a GPR survey made it possible to define the position, type, and composition of the detected targets. Two functionals, the Complex Matched Coherency Measure and the Complex Matched Analysis, turned out to be effective in computing coherency maps characterized by high-resolution and strong noise rejection, where velocity picking can be done with high precision. By using the 2D velocity field thus obtained, migration algorithms performed better than in the case of constant or 1D velocity field, with satisfactory collapsing of the diffracted events and moving of the reflected energy in the correct position. The varying velocity field was estimated on different lines and used to migrate all the GPR profiles composing the survey covering the entire archaeological area. The time slices built with the migrated profiles resulted in a higher S/N than those obtained from non-migrated or migrated at constant velocity GPR profiles. The improvements are inherent to the resolution, continuity, and energy content of linear reflective areas. On the basis of our experience, we can state that the use of high-resolution coherency functionals leads to migrated GPR profiles with a high-grade of hyperbolas focusing. These profiles favor better imaging of the targets of interest, thereby allowing for a more reliable interpretation

Stochastic hyperbola fitting, probabilistic inversion, reverse-time migration and clustering: A novel interpretation toolbox for in-situ planetary radar

Ground-penetrating radar (GPR) is becoming a mainstream tool in planetary exploration, and one of the few in-situ planetary geophysical methods. There are currently three missions (Perseverance, Tianwen-1, Chang'E-4) with GPR-equipped rovers, and two future missions (Chang'E-7, ExoMars) that will include GPR in their scientific payload. The large number of GPR data, combined with the novel setup of the measurements, creates the need for new data processing and interpretation techniques to address the unique challenges of in-situ planetary radar. The current paper proposes an interpretation pipeline that starts with a novel stochastic hyperbola fitting that estimates the probability kernel density of the bulk permittivity at different depths. Subsequently, the bulk permittivity distribution is transformed via a novel probabilistic inversion to a 1-dimensional (1D) permittivity profile. The inverted 1D permittivity profile is then used as an input to a bespoke reverse-time migration (RTM) using the finite-difference time-domain (FDTD) method. RTM using FDTD does not assume a clinical homogeneous half-space; instead, it accounts for the expected layered structure of the investigated medium. Lastly, the migrated radargram is clustered in order to identify subsurface targets and distinguish them from the background medium. Each of the processing steps has never been reported in planetary radar; and together act as a complete processing toolbox tuned for planetary science. The suggested interpretation pipeline is validated numerically in a 1D case study with a complex layered structure and multiple subsurface targets. The proposed processing scheme is then applied to the GPR data from the Chang'E-4 mission at the Von Karman crater, revealing a previously unseen layered structure and a complex distribution of rocks/boulders

Stochastic hyperbola fitting, probabilistic inversion, reverse-time migration and clustering : A novel interpretation toolbox for in-situ planetary radar

Peer reviewedPublisher PD

Characterization of components of water supply systems from GPR images and tools of intelligent data analysis

[EN] Over time, due to multiple operational and maintenance activities, the networks of water supply systems (WSSs) undergo interventions, modifications or even are closed. In many cases, these activities are not properly registered. Knowledge of the paths and characteristics (status and age, etc.) of the WSS pipes is obviously necessary for efficient and dynamic management of such systems. This problem is greatly augmented by considering the detection and control of leaks. Access to reliable leakage information is a complex task. In many cases, leaks are detected when the damage is already considerable, which brings high social and economic costs. In this sense, non-destructive methods (e.g., ground penetrating radar - GPR) may be a constructive response to these problems, since they allow, as evidenced in this thesis, to ascertain paths of pipes, identify component characteristics, and detect primordial water leaks. Selection of GPR in this work is justified by its characteristics as non-destructive technique that allows studying both metallic and non-metallic objects. Although the capture of information with GPR is usually successful, such aspects as the capture settings, the large volume of generated information, and the use and interpretation of such information require high level of skill and experience. This dissertation may be seen as a step forward towards the development of tools able to tackle the problem of lack of knowledge on the WSS buried assets. The main objective of this doctoral work is thus to generate tools and assess their feasibility of application to the characterization of components of WSSs from GPR images. In this work we have carried out laboratory tests specifically designed to propose, develop and evaluate methods for the characterization of the WSS buried components. Additionally, we have conducted field tests, which have enabled us to determine the feasibility of implementing such methodologies under uncontrolled conditions. The methodologies developed are based on techniques of intelligent data analysis. The basic principle of this work has involved the processing of data obtained through the GPR to look for useful information about WSS components, with special emphasis on the pipes. After performing numerous activities, one can conclude that, using GPR images, it is feasible to obtain more information than the typical identification of hyperbolae currently performed. In addition, this information can be observed directly, e.g. more simply, using the methodologies proposed in this doctoral work. These methodologies also prove that it is feasible to identify patterns (especially with the preprocessing algorithm termed Agent race) that provide fairly good approximation of the location of leaks in WSSs. Also, in the case of pipes, one can obtain such other characteristics as diameter and material. The main outcomes of this thesis consist in a series of tools we have developed to locate, identify and visualize WSS components from GPR images. Most interestingly, the data are synthesized and reduced so that the characteristics of the different components of the images recorded in GPR are preserved. The ultimate goal is that the developed tools facilitate decision-making in the technical management of WSSs, and that such tools can even be operated by personnel with limited experience in handling non-destructive methodologies, specifically GPR.[ES] Con el paso del tiempo, y debido a múltiples actividades operacionales y de mantenimiento, las redes de los sistemas de abastecimiento de agua (SAAs) sufren intervenciones, modificaciones o incluso, son clausuradas, sin que, en muchos casos, estas actividades sean correctamente registradas. El conocimiento de los trazados y características (estado y edad, entre otros) de las tuberías en los SAAs es obviamente necesario para una gestión eficiente y dinámica de tales sistemas. A esta problemática se suma la detección y el control de las fugas de agua. El acceso a información fiable sobre las fugas es una tarea compleja. En muchos casos, las fugas son detectadas cuando los daños en la red son ya considerables, lo que trae consigo altos costes sociales y económicos. En este sentido, los métodos no destructivos (por ejemplo, ground penetrating radar - GPR), pueden ser una respuesta a estas problemáticas, ya que permiten, como se pone de manifiesto en esta tesis, localizar los trazados de las tuberías, identificar características de los componentes y detectar las fugas de agua cuando aún no son significativas. La selección del GPR, en este trabajo se justifica por sus características como técnica no destructiva, que permite estudiar tanto objetos metálicos como no metálicos. Aunque la captura de información con GPR suele ser exitosa, la configuración de la captura, el gran volumen de información, y el uso y la interpretación de la información requieren de alto nivel de habilidad y experiencia por parte del personal. Esta tesis doctoral se plantea como un avance hacia el desarrollo de herramientas que permitan responder a la problemática del desconocimiento de los activos enterrados de los SAAs. El objetivo principal de este trabajo doctoral es, pues, generar herramientas y evaluar la viabilidad de su aplicación en la caracterización de componentes de un SAA, a partir de imágenes GPR. En este trabajo hemos realizado ensayos de laboratorio específicamente diseñados para plantear, elaborar y evaluar metodologías para la caracterización de los componentes enterrados de los SAAs. Adicionalmente, hemos realizado ensayos de campo, que han permitido determinar la viabilidad de aplicación de tales metodologías bajo condiciones no controladas. Las metodologías elaboradas están basadas en técnicas de análisis inteligentes de datos. El principio básico de este trabajo ha consistido en el tratamiento adecuado de los datos obtenidos mediante el GPR, a fin de buscar información de utilidad para los SAAs respecto a sus componentes, con especial énfasis en las tuberías. Tras la realización de múltiples actividades, se puede concluir que es viable obtener más información de las imágenes de GPR que la que actualmente se obtiene con la típica identificación de hipérbolas. Esta información, además, puede ser observada directamente, de manera más sencilla, mediante las metodologías planteadas en este trabajo doctoral. Con estas metodologías se ha probado que también es viable la identificación de patrones (especialmente el pre-procesado con el algoritmo Agent race) que proporcionan aproximación bastante acertada de la localización de las fugas de agua en los SAAs. También, en el caso de las tuberías, se puede obtener otro tipo de características tales como el diámetro y el material. Como resultado de esta tesis se han desarrollado una serie de herramientas que permiten visualizar, identificar y localizar componentes de los SAAs a partir de imágenes de GPR. El resultado más interesante es que los resultados obtenidos son sintetizados y reducidos de manera que preservan las características de los diferentes componentes registrados en las imágenes de GPR. El objetivo último es que las herramientas desarrolladas faciliten la toma de decisiones en la gestión técnica de los SAAs y que tales herramientas puedan ser operadas incluso por personal con una experiencia limitada en el manejo[CA] Amb el temps, a causa de les múltiples activitats d'operació i manteniment, les xarxes de sistemes d'abastament d'aigua (SAAs) se sotmeten a intervencions, modificacions o fins i tot estan tancades. En molts casos, aquestes activitats no estan degudament registrats. El coneixement dels camins i característiques (estat i edat, etc.) de les canonades d'aigua i sanejament fa evident la necessitat d'una gestió eficient i dinàmica d'aquests sistemes. Aquest problema es veu augmentat en gran mesura tenint en compte la detecció i control de fuites. L'accés a informació fiable sobre les fuites és una tasca complexa. En molts casos, les fugues es detecten quan el dany ja és considerable, el que porta costos socials i econòmics. En aquest sentit, els mètodes no destructius (per exemple, ground penetrating radar - GPR) poden ser una resposta constructiva a aquests problemes, ja que permeten, com s'evidencia en aquesta tesi, per determinar rutes de canonades, identificar les característiques dels components, i detectar les fuites d'aigua quan encara no són significatives. La selecció del GPR en aquest treball es justifica per les seves característiques com a tècnica no destructiva que permet estudiar tant objectes metàl·lics i no metàl·lics. Tot i que la captura d'informació amb GPR sol ser reeixida, aspectes com ara la configuració de captura, el gran volum d'informació que es genera, i l'ús i la interpretació d'aquesta informació requereix alt nivell d'habilitat i experiència. Aquesta tesi pot ser vista com un pas endavant cap al desenvolupament d'eines capaces d'abordar el problema de la manca de coneixement sobre els actius d'aigua i sanejament enterrat. L'objectiu principal d'aquest treball doctoral és, doncs, generar eines i avaluar la seva factibilitat d'aplicació a la caracterització dels components de los SAAs, a partir d'imatges GPR. En aquest treball s'han dut a terme proves de laboratori específicament dissenyats per proposar, desenvolupar i avaluar mètodes per a la caracterització dels components d'aigua i sanejament soterrat. A més, hem dut a terme proves de camp, que ens han permès determinar la viabilitat de la implementació d'aquestes metodologies en condicions no controlades. Les metodologies desenvolupades es basen en tècniques d'anàlisi intel·ligent de dades. El principi bàsic d'aquest treball ha consistit en el tractament de dades obtingudes a través del GPR per buscar informació útil sobre els components d'SAA, amb especial èmfasi en la canonades. Després de realitzar nombroses activitats, es pot concloure que, amb l'ús d'imatges de GPR, és factible obtenir més informació que la identificació típica d'hipèrboles realitzat actualment. A més, aquesta informació pot ser observada directament, per exemple, més simplement, utilitzant les metodologies proposades en aquest treball doctoral. Aquestes metodologies també demostren que és factible per identificar patrons (especialment el pre-processat amb l'algoritme Agent race) que proporcionen bastant bona aproximació de la localització de fuites en SAAs. També, en el cas de tubs, es pot obtenir altres característiques com ara el diàmetre i el material. Els principals resultats d'aquesta tesi consisteixen en una sèrie d'eines que hem desenvolupat per localitzar, identificar i visualitzar els components dels SAAS a partir d'imatges GPR. El resultat més interessant és que els resultats obtinguts són sintetitzats i reduïts de manera que preserven les característiques dels diferents components registrats en les imatges de GPR. L'objectiu final és que les eines desenvolupades faciliten la presa de decisions en la gestió tècnica de SAA, i que tals eines poden fins i tot ser operades per personal amb poca experiència en el maneig de metodologies no destructives, específicament GPR.Ayala Cabrera, D. (2015). Characterization of components of water supply systems from GPR images and tools of intelligent data analysis [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/59235TESISPremios Extraordinarios de tesis doctorale

Development of Data Processing Tools for the Analysis of Radargrams in Utility Detection Using Ground Penetrating Radar, Journal of Telecommunications and Information Technology, 2017, nr 3

The extraction of quantitative information from Ground Penetrating Radar (GPR) data sets (radargrams) to detect and map underground utility pipelines is a challenging task. This study proposes several algorithms included in the main stages of a data processing chain associated with radargrams. It comprises preprocessing, hyperbola enhancing, hyperbola detection and localization, and parameter extraction. Additional parameters related to the GPR system such as the frequency band and the polarization bring data sets additional information that need to be exploited. Presently, the algorithms have been applied step by step on synthetic and experimental data. The results help to guide future developments in signal processing for quantitative parameter estimation