EP4172655A1 - Procede de detection de cible pour radar a penetration de sol et radar associe - Google Patents

Procede de detection de cible pour radar a penetration de sol et radar associe

Info

Publication number
EP4172655A1
EP4172655A1 EP21733471.3A EP21733471A EP4172655A1 EP 4172655 A1 EP4172655 A1 EP 4172655A1 EP 21733471 A EP21733471 A EP 21733471A EP 4172655 A1 EP4172655 A1 EP 4172655A1
Authority
EP
European Patent Office
Prior art keywords
target
hypothesis
ground
area
radar
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
EP21733471.3A
Other languages
German (de)
English (en)
Inventor
Raffaele D'ERRICO
Jean-Baptiste Dore
Gloria MAKHOUL
Luc Maret
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Commissariat a lEnergie Atomique CEA
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Commissariat a lEnergie Atomique CEA, Commissariat a lEnergie Atomique et aux Energies Alternatives CEA filed Critical Commissariat a lEnergie Atomique CEA
Publication of EP4172655A1 publication Critical patent/EP4172655A1/fr
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/411—Identification of targets based on measurements of radar reflectivity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414—Discriminating targets with respect to background clutter

Definitions

  • the invention relates to the field of ground penetrating radars or ground penetrating radars which cover all the techniques for detecting, locating or identifying underground targets by means of a radiofrequency system.
  • Underground targets are, for example, pipes of different diameter and nature (steel, PVC, cement, concrete, etc.) which can be buried at different depths.
  • ground penetrating radars One objective of ground penetrating radars is to locate such objects with precision in order to be able to correctly map a subsoil, for example for safety requirements during work.
  • a general problem targeted by the invention relates to the location of underground targets with precision regardless of the dielectric media traversed by the radar signals.
  • a first type of compatible localization method for ground-penetrating radars relates to so-called migration techniques derived from imagery. These methods are based on solving the wave equation by interpolation in the frequency domain.
  • the document [1] describes a radar imaging method for systems with one transmitting antenna and one receiving antenna (of the SISO type for "Single Input Single Output” in English).
  • Documents [4] and [6] give examples of radar imaging methods applied to systems having several transmit and receive antennas, that is to say antenna arrays (of the MIMO type for " Multiple Input Multiple Output ”in English).
  • a drawback of imaging or migration methods is that they generally only use the phase information of the signals picked up.
  • the amplitude information is not used neither to determine the position of a target nor to estimate the transmission losses of the signals in the medium traversed.
  • these methods aim to construct an image of an area and require post-processing to determine the position of a target.
  • Document [5] describes another method based on an antenna array to virtually create a plane wave. This solution also has the drawback of requiring transmission preprocessing.
  • Document [2] describes a MIMO-type radar technique which exploits the amplitude of the signals picked up but does not consider the influence of transmission losses during the propagation of signals in the ground. This simplification has the effect of limiting the location precision, in particular in the case of heterogeneous media and areas comprising several nearby targets to be detected.
  • the present invention aims to remedy the drawbacks of the aforementioned techniques of the state of the art by proposing a location method which takes into account the dielectric losses in the ground in order to improve the precision of detection and location of targets. .
  • the subject of the invention is a method for detecting at least one target buried in an area of the ground, using a ground penetrating radar, the radar comprising at least one transmitting antenna and at least one receiving antenna, the method comprising the steps of:
  • - define a first hypothesis H0 corresponding to an absence of target in the ground area, the channel matrix under this first hypothesis being modeled by noise, - define a second hypothesis H1 corresponding to the presence of at least one target in the ground area, the channel matrix under this second hypothesis being modeled by a propagation channel model depending at least on the characteristics of the transmitting antennas and receiving antennas and signal propagation losses in the ground,
  • the step of performing a test of the likelihood ratio between the first hypothesis H0 and the second hypothesis H1 comprises the sub-steps of:
  • the step of determining at least one estimator of the maximum likelihood of the position and of the radar equivalent area of a target, under the second hypothesis H1 is carried out by searching for at least an extremum of a function of three variables depending on the position of a target in the ground area, on the radar equivalent area of a target and on a signal transmission loss coefficient in the ground area.
  • the search for at least one extremum of said function of three variables is carried out by:
  • the likelihood of the channel matrix under the second hypothesis H1 is determined by an extremum of said function.
  • the likelihood of the channel matrix under the first hypothesis H0 is determined at least by calculating the energy of the channel matrix.
  • a target position is sought in a two or three-dimensional grid having a defined initial pitch, the method being iterated over several iterations by reducing the area and the initial pitch at each iteration.
  • the method according to the invention comprises the steps of:
  • the temporal windowing is defined so as to filter the contributions of the signals corresponding to a distance interval between the radar and the ground area.
  • the process is repeated for several different filters so as to cover the entire area of the ground.
  • the method according to the invention further comprises a step of canceling the contribution of a target detected in the measurements of signals acquired, for example by means of a cancellation algorithm. echo, the method being iterated over several successive iterations, canceling on each new iteration the echo of the target detected in the previous iteration.
  • the invention also relates to a ground penetrating radar comprising at least one transmitting antenna and at least one receiving antenna and a device for detecting at least one target buried in an area of the ground configured to execute the steps of the detection method according to the invention.
  • the invention also relates to a computer program comprising instructions for the execution of the method according to the invention, when the program is executed by a processor as well as a recording medium readable by a processor on which is recorded a program comprising instructions for the execution of the method according to the invention, when the program is executed by a processor.
  • Figure 1 is an illustrative diagram of the operation of a ground penetrating radar
  • FIG. 2 represents a general flowchart describing the steps of a target detection method according to one embodiment of the invention
  • FIG. 3 represents a flowchart describing an embodiment of the target detection step of the method according to the invention
  • FIG. 4 represents an illustrative diagram of an example of iterative operation of the invention
  • FIG. 5 represents a flowchart describing the steps of an alternative embodiment of the method according to the invention
  • FIG. 6 represents a diagram of a temporal response of a propagation channel to illustrate an alternative embodiment of the invention
  • FIG. 7 represents a flowchart describing the steps of another variant embodiment of the method according to the invention.
  • FIG. 8 is a diagram of an example of a ground penetrating radar capable of carrying out the invention.
  • FIG. 1 shows a simplified diagram of the operating principle of an RDR ground penetrating radar.
  • a radar comprises one or more transmitting antennas and one or more receiving antennas.
  • the radar has five transmitting antennas and four receiving antennas. Without departing from the scope of the invention, the number of transmitting and receiving antennas can take any non-zero integer value.
  • An area of the ground broken down into a grid G is analyzed by ground penetrating radar RDR by performing a sequence of transmissions and signal acquisitions for each pair consisting of a transmitting antenna and a receiving antenna.
  • an acquisition sequence consists in transmitting a radar signal from the first antenna Tx1 in transmission, this signal is reflected on a possible target P located in the zone G then back-propagated towards each reception antenna Rx1, Rx2, Rx3 , Rx4. This sequence is repeated for each transmitting antenna Tx2, Tx3, Tx4, Tx5.
  • the target search area can be a two or three dimensional area.
  • the radar can be mobile, for example on board a vehicle, to cover different areas of the ground.
  • the invention is implemented by a radar of the type of Figure 1.
  • Figure 2 illustrates the main steps of the target detection method according to the invention.
  • the first step 201 consists in acquiring the signals during a transmission-reception sequence as described above.
  • the acquisition sequence performs a frequency sweep according to a set of predefined frequencies.
  • the transmission-reception sequences on all the antennas are repeated for several frequencies defining a frequency band to be analyzed.
  • the acquisition sequence can be sequential on transmission as described above but can also be carried out by simultaneously transmitting several orthogonal signals on all the transmission antennas so as to directly carry out a simultaneous acquisition of the various echoes of each transmitting antenna for a given frequency.
  • the scanning in successive frequencies can be replaced by the use of a waveform of the OFDM (Orthogonal Frequency Division Multiplexing) type which makes it possible to code in transmission the signals by distribution on orthogonal frequency channels in the form multiple sub-carriers.
  • OFDM Orthogonal Frequency Division Multiplexing
  • the width of the channels or sub-bands and their spacing can be configured. These parameters can vary for each sub-band according to the signal-to-noise ratio specific to this sub-band.
  • a step 202 of estimating channel matrices is performed.
  • Signal acquisition is triggered synchronously with transmission.
  • the signals are transposed to baseband and then demodulated via OFDM demodulation.
  • a channel estimation is carried out on each subcarrier, for example by correlating the symbols received with the known symbols transmitted on this subcarrier.
  • the channel estimation is performed on several OFDM symbols by consistently averaging the result to reduce the impact of thermal noise or interference.
  • a channel estimation is carried out for each propagation channel defined by a link between a transmission antenna Txi and a reception antenna Rxj but also for each frequency sub-band.
  • the complete frequency response is obtained by concatenation of the responses of each sub-band.
  • a set of channel matrices H r (f) is obtained, for several frequencies, each matrix having a dimension equal to the number of reception antennas by the number d 'transmitting antennas.
  • K is the number of targets for which an echo is received.
  • n (f) denotes the noise matrix for the frequency f.
  • the invention is in particular based on a particular modeling of the channel matrix as detailed below, which takes into account the dielectric losses in the medium of which the soil is made.
  • the noise-free channel matrix H (f, p k ) can be expressed using the following relation, where is the term-by-term multiplication operator between two matrices.
  • U (f, p k , L) is the propagation delay matrix whose elements (i, j) are expressed using the following relation:
  • v is the speed of propagation of the wave in the medium crossed by the supposedly homogeneous wave.
  • the medium traversed is partly made up of air and soil layers or only soil layers if the transmitting antennas are positioned directly against the ground surface.
  • L is a loss coefficient which depends on the medium crossed by the wave. Depending on the medium (earth, sand, rock) of which the soil is made, the emitted wave undergoes more or less significant attenuations that should be taken into account to correctly model the channel matrix.
  • the matrix A (f, p k , L) represents the matrix of the amplitudes of the echoes reflected on the target. This matrix integrates the losses in the medium, the antenna gains and the radar equivalent area of the target.
  • Y (/, p fc ) models the radar equivalent surface of the target at the position p k .
  • G t (f, p k ) models the transmission antenna pattern in the direction of the position of the target.
  • G r (f, p k ) models the reception antenna pattern in the direction of the position of the target.
  • the two matrices G t (f, p k ) and G r (f, p k ) are known and are determined from the characteristics of the antennas.
  • B 0 (f, p k ) and B 1 (f, p k , L) are two matrices which together model the propagation losses of spherical waves in the ground.
  • the invention differs in particular from the methods of the state of the art by taking into account the loss coefficient L in the middle of the soil.
  • An objective of the detection method by radar processing, according to the invention is to determine the position of a target having knowledge of the frequency response of the monitored area.
  • the third step 203 of the method according to the invention consists in carrying out a target detection test from the modeling of the propagation channel.
  • H 0 corresponds to a hypothesis of the absence of a target
  • Hi corresponds to an hypothesis of the presence of a target.
  • the channel matrix can be expressed under the two assumptions using the following relationships:
  • H 0 ) denotes the probability that the received channel matrix corresponds to the hypothesis H 0 .
  • H 1 ) denotes the probability that the received channel matrix corresponds to the hypothesis H 1 .
  • the sign oc means "proportional to”.
  • the Neyman-Pearson criterion is used which was developed to maximize the probability of detection while ensuring a given probability of false alarm. This criterion consists in comparing the ratio of the likelihoods of the channel matrices, under each of the hypotheses, with a predetermined threshold l. This is a so-called LRT test for "Likelihood Ratio Test".
  • a target is detected when the following relationship is verified:
  • Relation (11) is also written in the form:
  • FIG. 3 details the steps necessary for carrying out the detection test 203.
  • the first step 301 consists in estimating the unknown parameters, that is to say the position of the target and its radar equivalent surface, as well as the loss coefficient L, in the sense of maximum likelihood.
  • the observation zone or space is discretized in the form of a grid G having a predefined pitch (see Figure 1).
  • a detection is declared for any cell of the grid validating the test under the H-i hypothesis.
  • the maximum likelihood estimator of the parameters y and p is obtained by finding the minimum of the following function:
  • step 301 consists in finding the minimum of a function g (y, p k , L, /), the function g being defined from relation (14).
  • the function g can have several local minima in the space defined by the three parameters (Y, p k , L).
  • a first variant of the invention consists in searching for several local minima (for example the N preponderant minima, with N a predefined integer) which will each correspond to a potential target.
  • Another variant of the invention consists in being limited to the search for a single local minimum (the most important), the method described in Figures 2 and 3 being applied for the search for a single target.
  • the search for one or more local minima can be carried out by an optimization function available to those skilled in the art, for example one of the digital resolution functions described in one of the references [7], [8], [9].
  • the estimate of y (f) is calculated using relation (18) for a set of assumptions of position p k and of loss coefficient L.
  • positions P k vary over all the positions of the grid G superimposed on the zone to be analyzed. These positions can be defined in two or three dimensions.
  • the values of the coefficient L can be chosen from among several hypotheses, taking into account the a priori knowledge of the composition of the soil.
  • Step 301 can be refined by successive iterations by redefining, at each iteration, a new grid G2 of finer pitch around the estimate of the position p k obtained in the previous iteration as illustrated in figure 4.
  • step 301 one (or more) estimate is obtained which is then used to carry out a detection test via equation (13).
  • a step 302 the likelihood of the channel matrix received under the Hi hypothesis is calculated. This calculation is in fact directly given by the value of the estimate of
  • a step 303 the likelihood of the channel matrix received is calculated under the hypothesis H 0 . This calculation is carried out by calculating the term which corresponds to the sum or the average, over all the frequencies, of the squared modulus of the channel matrix. It corresponds to the energy in frequency of the channel matrix.
  • the norm operator of a complex matrix is defined by the trace of the product of the matrix and the conjugated transpose matrix:
  • a step 304 the likelihood test of equation (13) is carried out by comparing the difference between the terms calculated respectively in steps 302 and 303, with a detection threshold l.
  • the speed v of propagation of the wave in the media crossed can be estimated by an external method or else can be considered as an additional variable of the minimization problem.
  • Figure 5 shows schematically an alternative embodiment of the invention suitable for the detection of multiple targets.
  • the target detection method described in Figures 2 and 3 is applied to search for a single target (step 501), that is to say the target having the preponderant echo.
  • an echo cancellation algorithm is applied to the signals received (or directly to the channel matrix H r (f)) in order to cancel the contribution of the reflections of the signal on the target detected in step 501.
  • One possible algorithm is the algorithm known by the acronym CLEAN, for example described in reference [3].
  • the iterations are stopped via a stop criterion 503 which can consist of a fixed number of iterations or of a criterion depending on the radar equivalent area or the energy of the echoes of the detected targets.
  • the criterion can also be defined from a difference in energy or radar equivalent area between two targets detected in two successive iterations.
  • the image of the analyzed area is supplied as the input of an optimization algorithm 504 to improve the accuracy of the image.
  • the optimization algorithm is, for example, chosen from those described in references [7], [8] or [9].
  • the algorithm can be based on the stochastic gradient method.
  • the basic assumption for the detection of a target via a test under two assumptions as described previously, is that the noise statistic corresponds to an additive Gaussian white noise. In the case of several targets, this assumption is no longer true, which leads to bias in the detection algorithm.
  • This variant consists in applying a filter in the form of temporal windowing to the channel matrices in order to separately treat different sub-zones of the zone of the ground to be analyzed and thus to overcome the problems of dynamics.
  • h a vector consisting of the estimate of the frequency propagation channel for a pair of transmit and receive antennas
  • this vector is converted into the time domain via an inverse Fourier transform IFFT (step 701).
  • a filter is applied in the time domain (step 702) to the time response of the channel.
  • any type of window can be used, for example a Flanning window, or even a Blackman-Flarris window.
  • the result is converted (step 703) into the frequency domain via a direct Fourier transform then the new frequency vector obtained is used to constitute the channel matrix H r (f) at the input of the detection method according to one of the variants described above.
  • the process described in FIG. 7 is applied for each propagation channel connecting one of the transmitting antennas to one of the receiving antennas.
  • the detection method in its entirety, is iterated for several different time windows, for example by varying the windows over the entire area to be analyzed by circularly permuting the filter w of p / 2 positions at each iteration. where p is equal to the number of 1.
  • FIG. 6 shows schematically, on a diagram, an example of the temporal response 601 of a propagation channel between a transmitting antenna and a receiving antenna. On the same diagram, several time windows have been shown staggered and corresponding to several iterations of the method described above.
  • FIG. 8 is a diagram of an example of a ground penetrating radar configured to implement the invention.
  • the radar consists of an EM transmitter and a REC receiver, each provided with at least one antenna A1, A2.
  • the EM transmitter comprises a generator of a sequence of GEN bits or symbols, an OFDM MOD, a frequency converter RF1 to transpose the signal to a carrier frequency and an amplifier AMP.
  • the receiver REC comprises a frequency converter RF2 for transposing the received signal into baseband, an OFDM DEMOD demodulator and a processing unit consisting of a channel estimation module EST which determines the channel matrices from of the demodulated signals and of the sequence emitted by the generator GEN, an integrator INT which carries out an integration (average) and a concatenation of the frequency sub-bands and a target detection module DET which implements the processing steps described above for perform any one of the embodiments of the invention.
  • a channel estimation module EST which determines the channel matrices from of the demodulated signals and of the sequence emitted by the generator GEN
  • an integrator INT which carries out an integration (average) and a concatenation of the frequency sub-bands
  • a target detection module DET which implements the processing steps described above for perform any one of the embodiments of the invention.
  • the EST channel estimation, INT integration and DET detection modules can be produced in software and / or hardware form, in particular by using one or more processor (s) and one or more memory (s).
  • the processor can be a generic processor, a specific processor, an integrated circuit specific to an application (also known under the English name of ASIC for “Application-Specific Integrated Circuit”) or an array of programmable gates in situ (also known as the English name of FPGA for “Field-Programmable Gâte Array”).

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

Procédé de détection d'au moins une cible enterrée dans une zone du sol, à l'aide d'un radar à pénétration de sol, le procédé comprenant les étapes de : - Acquérir (201) une mesure d'un signal émis par chaque antenne d'émission et réfléchi dans la zone du sol, - Déterminer (202) une estimée, dans le domaine fréquentiel, du canal de propagation pour chaque couple constitué d'une antenne d'émission et d'une antenne de réception, - définir une première hypothèse H0 correspondant à une absence de cible dans la zone du sol, - définir une seconde hypothèse H1 correspondant à la présence d'au moins une cible dans la zone du sol, - réaliser (203) un test du rapport de vraisemblance entre la vraisemblance de la matrice de canal sous la seconde hypothèse H1 et la vraisemblance de la matrice de canal sous la première hypothèse H0, pour conclure à la présence d'une cible en une position donnée dans la zone du sol.

Description

DESCRIPTION
Titre de l’invention : Procédé de détection de cible pour radar à pénétration de sol et radar associé
[0001] L’invention concerne le domaine des radars à pénétration de sol ou géoradars qui couvrent l’ensemble des techniques permettant de détecter, localiser ou identifier des cibles souterraines au moyen d’un système radiofréquence.
[0002] Les cibles souterraines sont, par exemple, des canalisations de diamètre et nature différentes (acier, PVC, ciment, béton...) qui peuvent être enterrées à différentes profondeurs.
[0003] Un objectif des radars à pénétration de sol est de localiser de tels objets avec précision afin de pouvoir cartographier correctement un sous-sol, par exemple pour des besoins de sécurité lors de travaux.
[0004] Un problème général visé par l’invention concerne la localisation de cibles souterraines avec précision quels que soient les milieux diélectriques traversés par les signaux radars.
[0005] Un premier type de méthode de localisation compatible des radars à pénétration de sol concerne des techniques dites de migration issues de l’imagerie. Ces méthodes sont basées sur la résolution de l’équation d’onde par une interpolation dans le domaine fréquentiel.
[0006] Le document [1] décrit une méthode d’imagerie radar pour des systèmes à une antenne d’émission et une antenne de réception (de type SISO pour « Single Input Single Output » en anglais). Les documents [4] et [6] donnent des exemples de méthodes d’imagerie radar appliquées à des systèmes ayant plusieurs antennes d’émission et de réception, c’est-à-dire des réseaux d’antennes (de type MIMO pour « Multiple Input Multiple Output » en anglais).
[0007] Un inconvénient des méthodes d’imagerie ou de migration est qu’elles n’exploitent généralement que l’information de phase des signaux captés. L’information d’amplitude n’est pas utilisée ni pour déterminer la position d’une cible ni pour estimer les pertes de transmission des signaux dans le milieu traversé. Par ailleurs, ces méthodes ont pour objectif de construire une image d’une zone et nécessitent un post traitement pour déterminer la position d’une cible.
[0008] Le document [5] décrit une autre méthode basée sur un réseau d’antennes pour créer virtuellement une onde plane. Cette solution présente également l’inconvénient de nécessiter un prétraitement en émission.
[0009] Le document [2] décrit une technique radar de type MIMO qui exploite l’amplitude des signaux captés mais ne considère pas l’influence des pertes de transmission lors de la propagation des signaux dans le sol. Cette simplification a pour effet de limiter la précision de localisation notamment dans le cas de milieux hétérogènes et de zones comportant plusieurs cibles proches à détecter.
[0010] La présente invention vise à remédier aux inconvénients des techniques de l’état de l’art précitées en proposant une méthode de localisation qui prend en compte les pertes diélectriques dans le sol afin d’améliorer la précision de détection et localisation des cibles.
[0011] L’invention a pour objet un procédé de détection d’au moins une cible enterrée dans une zone du sol, à l’aide d’un radar à pénétration de sol, le radar comprenant au moins une antenne d’émission et au moins une antenne de réception, le procédé comprenant les étapes de :
- Acquérir, sur chaque antenne de réception, une mesure d’un signal émis par chaque antenne d’émission et réfléchi dans la zone du sol,
- Déterminer, à partir desdites mesures, une estimée, dans le domaine fréquentiel, du canal de propagation pour chaque couple constitué d’une antenne d’émission et d’une antenne de réception, lesdites estimées étant regroupées dans une matrice de canal ayant une première dimension égale au nombre d’antennes de réception, une seconde dimension égale au nombre d’antennes d’émission et une troisième dimension égale au nombre de fréquences,
- définir une première hypothèse H0 correspondant à une absence de cible dans la zone du sol, la matrice de canal sous cette première hypothèse étant modélisée par du bruit, - définir une seconde hypothèse H1 correspondant à la présence d’au moins une cible dans la zone du sol, la matrice de canal sous cette seconde hypothèse étant modélisée par un modèle de canal de propagation dépendant au moins des caractéristiques des antennes d’émission et des antennes de réception et des pertes de propagation du signal dans le sol,
- réaliser un test du rapport de vraisemblance entre la vraisemblance de la matrice de canal sous la seconde hypothèse H1 et la vraisemblance de la matrice de canal sous la première hypothèse H0, pour conclure à la présence d’une cible en une position donnée dans la zone du sol.
[0012] Selon un aspect particulier de l’invention, l’étape de réaliser un test du rapport de vraisemblance entre la première hypothèse H0 et la seconde hypothèse H1 comprend les sous-étapes de :
- Déterminer au moins un estimateur du maximum de vraisemblance de la position et de la surface équivalente radar d’une cible, sous la seconde hypothèse H1 ,
- Déterminer la vraisemblance de la matrice de canal sous la seconde hypothèse H1 pour chaque estimateur,
- Déterminer la vraisemblance de la matrice de canal sous la première hypothèse H0,
- comparer le ratio entre les deux vraisemblances à un seuil de détection prédéterminé pour en déduire la présence ou l’absence d’une cible pour chaque position estimée.
[0013] Selon un aspect particulier de l’invention, l’étape de déterminer au moins un estimateur du maximum de vraisemblance de la position et de la surface équivalente radar d’une cible, sous la seconde hypothèse H1 est réalisée en recherchant au moins un extremum d’une fonction de trois variables dépendant de la position d’une cible dans la zone du sol, de la surface équivalente radar d’une cible et d’un coefficient de pertes de transmission du signal dans la zone du sol. [0014] Selon un aspect particulier de l’invention, la recherche d’au moins un extremum de ladite fonction de trois variables est réalisée en :
- déterminant un estimateur du maximum de vraisemblance de la surface équivalente radar en calculant une corrélation normalisée entre la matrice de canal et une matrice de pertes incluant les fonctions de transfert des antennes d’émission et de réception, les pertes de propagation du signal dans l’espace libre et dans la zone du sol,
- recherchant au moins un extremum selon la position et le coefficient de pertes de ladite fonction évaluée pour la surface équivalente radar estimée.
[0015] Selon un aspect particulier de l’invention, la vraisemblance de la matrice de canal sous la seconde hypothèse H1 est déterminée par un extremum de ladite fonction.
[0016] Selon un aspect particulier de l’invention, la vraisemblance de la matrice de canal sous la première hypothèse H0 est déterminée au moins en calculant l’énergie de la matrice de canal.
[0017] Selon un aspect particulier de l’invention, une position de cible est recherchée dans une grille en deux ou trois dimensions ayant un pas initial défini, le procédé étant itéré sur plusieurs itérations en réduisant la zone et le pas initial à chaque itération.
[0018] Dans une variante de réalisation, le procédé selon l’invention comprend les étapes de :
- Convertir les estimées de canal de propagation pour chaque couple constitué d’une antenne d’émission et d’une antenne de réception, dans le domaine temporel,
- Appliquer un filtre aux estimées sous la forme d’un fenêtrage temporel,
- Convertir les estimées filtrées dans le domaine fréquentiel,
- Regrouper les estimées filtrées dans une matrice de canal [0019] Selon un aspect particulier de l’invention, le fenêtrage temporel est défini de manière à filtrer les contributions des signaux correspondant à un intervalle de distance entre le radar et la zone du sol.
[0020] Selon un aspect particulier de l’invention, le procédé est répété pour plusieurs filtres différents de manière à couvrir toute la zone du sol.
[0021] Dans une variante de réalisation, le procédé selon l’invention comprend en outre une étape d’annulation de la contribution d’une cible détectée dans les mesures de signaux acquis, par exemple au moyen d’un algorithme d’annulation d’écho, le procédé étant itéré sur plusieurs itérations successives, en annulant à chaque nouvelle itération l’écho de la cible détectée à l’itération précédente.
[0022] L’invention a aussi pour objet un radar à pénétration de sol comprenant au moins une antenne d’émission et au moins une antenne de réception et un dispositif de détection d’au moins une cible enterrée dans une zone du sol configuré pour exécuter les étapes du procédé de détection selon l’invention.
[0023] L’invention a aussi pour objet un programme d'ordinateur comportant des instructions pour l'exécution du procédé selon l’invention, lorsque le programme est exécuté par un processeur ainsi qu’un support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution du procédé selon l’invention, lorsque le programme est exécuté par un processeur.
[0024] D’autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit en relation aux dessins annexés suivants.
[0025] [Fig. 1] la figure 1 représente un schéma illustratif du fonctionnement d’un radar à pénétration de sol,
[0026] [Fig. 2] la figure 2 représente un organigramme général décrivant les étapes d’un procédé de détection de cible selon un mode de réalisation de l’invention,
[0027] [Fig. 3] la figure 3 représente un organigramme décrivant un mode de réalisation de l’étape de détection de cible du procédé selon l’invention,
[0028] [Fig. 4] la figure 4 représente un schéma illustratif d’un exemple de fonctionnement itératif de l’invention, [0029] [Fig. 5] la figure 5 représente un organigramme décrivant les étapes d’une variante de réalisation du procédé selon l’invention,
[0030] [Fig. 6] la figure 6 représente un diagramme d’une réponse temporelle d’un canal de propagation pour illustrer une variante de réalisation de l’invention,
[0031] [Fig. 7] la figure 7 représente un organigramme décrivant les étapes d’une autre variante de réalisation du procédé selon l’invention,
[0032] [Fig. 8] la figure 8 représente un schéma d’un exemple de radar à pénétration de sol apte à exécuter l’invention.
[0033] La figure 1 représente un schéma simplifié du principe de fonctionnement d’un radar à pénétration de sol RDR. Un tel radar comporte une ou plusieurs antennes d’émission et une ou plusieurs antennes de réception. Sur l’exemple de la figure 1 , le radar comporte cinq antennes en émission et quatre antennes en réception. Sans sortir du cadre de l’invention, le nombre d’antennes en émission et en réception peut prendre n’importe quelle valeur entière non nulle.
[0034] Une zone du sol décomposée en une grille G est analysée par le géoradar RDR en effectuant une séquence d’émissions et d’acquisitions de signaux pour chaque couple constitué d’une antenne d’émission et d’une antenne de réception. Par exemple, une séquence d’acquisition consiste à émettre un signal radar depuis la première antenne Tx1 en émission, ce signal est réfléchi sur une cible éventuelle P située dans la zone G puis rétro-propagé vers chaque antenne de réception Rx1 ,Rx2,Rx3,Rx4. Cette séquence est répétée pour chaque antenne d’émission Tx2,Tx3,Tx4,Tx5.
[0035] La zone de recherche de cible peut être une zone en deux ou trois dimensions.
[0036] Le radar peut être mobile, par exemple embarqué sur un véhicule, pour couvrir différentes zones du sol.
[0037] L’invention est mise en oeuvre par un radar du type de la figure 1.
[0038] La figure 2 illustre les principales étapes du procédé de détection de cibles selon l’invention.
[0039] La première étape 201 consiste à acquérir les signaux lors d’une séquence d’émission-réception telle que décrite ci-dessus. [0040] La séquence d’acquisition réalise un balayage en fréquence selon un ensemble de fréquences prédéfini. Autrement dit les séquences d’émission-réception sur l’ensemble des antennes sont répétées pour plusieurs fréquences définissant une bande de fréquences à analyser.
[0041] La séquence d’acquisition peut être séquentielle à l’émission comme décrit ci- dessus mais peut aussi être réalisée en émettant simultanément plusieurs signaux orthogonaux sur toutes les antennes d’émission de sorte à réaliser directement une acquisition simultanée des différents échos de chaque antenne émettrice pour une fréquence donnée.
[0042] Par ailleurs le balayage en fréquences successives peut être remplacé par l’utilisation d’une forme d’onde de type OFDM (Orthogonal Frequency Division Multiplexing) qui permet de coder en émission les signaux par répartition sur des canaux fréquentiels orthogonaux sous forme de multiples sous-porteuses. La largeur des canaux ou sous-bandes et leur espacement peut être paramétré. Ces paramètres peuvent varier pour chaque sous-bande selon le rapport signal à bruit propre à cette sous-bande.
[0043] A partir des signaux reçus, une étape 202 d’estimation de matrices de canal est réalisée. L’acquisition des signaux est déclenchée de manière synchrone avec émission. Après acquisition, les signaux sont transposés en bande de base puis démodulés via une démodulation OFDM.
[0044] Après démodulation OFDM, une estimation de canal est réalisée sur chaque sous-porteuse, par exemple en corrélant les symboles reçus avec les symboles connus émis sur cette sous-porteuse.
[0045] L’estimation de canal est réalisée sur plusieurs symboles OFDM en moyennant de façon cohérente le résultat pour réduire l’impact du bruit thermique ou des interférences.
[0046] Une estimation de canal est réalisée pour chaque canal de propagation défini par une liaison entre une antenne d’émission Txi et une antenne de réception Rxj mais aussi pour chaque sous-bande de fréquence.
[0047] La réponse en fréquence complète est obtenue par concaténation des réponses de chaque sous-bande. [0048] Ainsi, à l’issue de l’étape 202, on obtient un ensemble de matrices de canal Hr(f), pour plusieurs fréquences, chaque matrice ayant une dimension égale au nombre d’antennes de réception par le nombre d’antennes d’émission.
[0049] On peut écrire l’expression de la matrice de canal estimée en réception en fonction de la fréquence au moyen de la relation suivante :
[0050]
[0051] H(f,pk) est la matrice de canal, hors bruit, pour la fréquence f et une cible située en position pk = { xk, yk, zk}. K est le nombre de cibles pour lesquelles un écho est reçu. n(f) désigne la matrice de bruit pour la fréquence f.
[0052] L’invention est notamment basée sur une modélisation particulière de la matrice de canal telle que détaillée ci-dessous, qui prend en compte les pertes diélectriques dans le milieu dont est constitué le sol.
[0053] La matrice de canal hors bruit H(f,pk) peut s’exprimer à l’aide de la relation suivante, où est l’opérateur de multiplication termes à termes entre deux matrices.
[0055] U(f, pk, L) est la matrice des retards de propagation dont les éléments (i,j) s’expriment à l’aide de la relation suivante :
[0056] désigne l’unité imaginaire.
[0058] v est la vitesse de propagation de l’onde dans le milieu traversé par l’onde supposé homogène. Selon le positionnement du radar, le milieu traversé est en partie constitué par l’air et par les couches du sol ou uniquement par les couches du sol si les antennes d’émission sont positionnées directement contre la surface du sol.
[0059] est la somme des distances respectives entre l’antenne d’émission d’indice i et la cible en position pk d’une part et cette cible et l’antenne de réception d’indice j d’autre part.
[0060] L est un coefficient de pertes qui dépend du milieu traversé par l’onde. Selon le milieu (terre, sable, roche) dont le sol est constitué, l’onde émise subit des atténuations plus ou moins importantes qu’il convient de prendre en compte pour modéliser correctement la matrice de canal.
[0061] La matrice A(f,pk,L ) représente la matrice des amplitudes des échos réfléchis sur la cible. Cette matrice intègre les pertes dans le milieu, les gains d’antennes et la surface équivalente radar de la cible.
[0062] Cette matrice est exprimée à l’aide de la relation suivante :
[0063]
[0064] y(/,pfc) modélise la surface équivalente radar de la cible à la position pk.
[0065] Gt(f,pk) modélise le diagramme d’antenne d’émission dans la direction de la position de la cible.
[0066] Gr(f,pk) modélise le diagramme d’antenne de réception dans la direction de la position de la cible.
[0067] Les deux matrices Gt(f,pk ) et Gr(f,pk ) sont connues et sont déterminées à partir des caractéristiques des antennes.
[0068] B0(f, pk) et B1(f,pk,L) sont deux matrices qui modélisent ensemble les pertes de propagation des ondes sphériques dans le sol.
[0070] est le produit des distances respectives entre l’antenne d’émission d’indice i et la cible en position pk d’une part et cette cible et l’antenne de réception d’indice j d’autre part.
[0072] L’invention se différencie notamment des méthodes de l’état de l’art par la prise en compte du coefficient de pertes L dans le milieu du sol.
[0073] Un objectif de la méthode de détection par traitement radar, selon l’invention, est de déterminer la position d’une cible ayant connaissance de la réponse fréquentielle de la zone surveillée.
[0074] Ainsi, la troisième étape 203 du procédé selon l’invention consiste à réaliser un test de détection de cible à partir de la modélisation du canal de propagation introduite ci-dessus, des matrices de canal déterminées à l’étape 202 et en considérant un problème de détection à deux hypothèses H0 et H1. H0 correspond à une hypothèse d’absence de cible et H-i correspond à une hypothèse de présence de cible.
[0075] La matrice de canal peut être exprimée sous les deux hypothèses à l’aide des relations suivantes :
[0076] H0: H(f,pk) = n(f) (7)
[0077]
[0078] Sous l’hypothèse d’un bruit blanc non corrélé, on peut exprimer les log- rapports de vraisemblance de la matrice de canal sous chacune des deux hypothèses via les relations suivantes :
[0079]
[0080]
[0081]Pr( Hr| H0) désigne la probabilité que la matrice de canal reçue corresponde à l’hypothèse H0.
[0082] Pr(Hr| H1) désigne la probabilité que la matrice de canal reçue corresponde à l’hypothèse H1.
[0083] Le signe oc signifie « proportionnel à ».
[0084] On utilise le critère de Neyman-Pearson qui a été développé pour maximiser la probabilité de détection tout en assurant une probabilité de fausse alarme donnée. Ce critère consiste à comparer le rapport des vraisemblances des matrices de canal, sous chacune des hypothèses, à un seuil l prédéterminé. Il s’agit d’un test dit LRT pour « Likelihood Ratio Test ».
[0085] Une cible est détectée lorsque la relation suivante est vérifiée :
[0086]
[0087] La relation (11 ) s’écrit aussi sous la forme :
[0088]
[0089] En pratique, les matrices A et U sont inconnues et le log rapport de vraisemblance pour l’hypothèse H1 ne peut pas être calculé directement. [0090] Ainsi, on propose d’utiliser un test dit GLRT pour « Generalized Likelihood Ratio Test » consistant à injecter dans le test LRT défini par l’équation (12) une estimation des paramètres inconnus du rapport de vraisemblance pour l’hypothèse H-i au sens du maximum de vraisemblance.
[0091] Le test défini à l’équation (11) devient alors :
[0092]
[0093] La figure 3 détaille les étapes nécessaires à la réalisation du test de détection 203.
[0094] La première étape 301 consiste à estimer les paramètres inconnus, c’est-à- dire la position de la cible et sa surface équivalente radar, ainsi que le coefficient de pertes L, au sens du maximum de vraisemblance.
[0095] En l’absence de connaissance de la position exacte de la cible, la zone ou l’espace d’observation est discrétisé sous forme d’une grille G ayant un pas prédéfini (cf figure 1). Une détection est déclarée pour toute case de la grille validant le test sous l’hypothèse H-i.
[0096] L’estimateur du maximum de vraisemblance des paramètres y et p est obtenu en recherchant le minimum de la fonction suivante :
[0097]
[0098] En effet la minimisation de la fonction donnée par la relation (14) est équivalente à la maximisation du log rapport de vraisemblance logPr(Hr\H1 ) (tel que défini à l’équation (13)).
[0099] En développant la relation (14) à partir du modèle des matrices A et U (relations (3) et (4)), l’étape 301 consiste à rechercher le minimum d’une fonction g(y, pk, L, /), la fonction g étant définie à partir de la relation (14).
[0100 \ g(y,pk,L,f) =
[0101 ]
[0102] La fonction g peut présenter plusieurs minima locaux dans l’espace défini par les trois paramètres ( Y,pk,L ). [0103] Une première variante de l’invention consiste à rechercher plusieurs minima locaux (par exemple les N minima prépondérants, avec N un entier prédéfini) qui correspondront chacun à une cible potentielle.
[0104] Une autre variante de l’invention consiste à se limiter à la recherche d’un seul minimum local (le plus important), le procédé décrit aux figures 2 et 3 étant appliqué pour la recherche d’une seule cible.
[0105] La recherche d’un ou plusieurs minima locaux peut être réalisée par une fonction d’optimisation à la disposition de l’Homme du métier, par exemple une des fonctions de résolution numérique décrites dans l’une des références [7], [8], [9].
[0106] Une autre solution moins coûteuse en opérations, consiste à effectuer cette recherche en plusieurs étapes.
[0107] Tout d’abord, on différencie la fonction g par rapport à la variable y(/), en recherchant la valeur de y(/) qui est solution de l’équation suivante :
[0109] Une solution de l’équation (17) est donnée par :
[0111] Ainsi une estimée de la surface équivalente radar peut être calculée à partir de la relation (18) en réalisant une corrélation normalisée entre la matrice de canal reçue et la matrice M définie par l’équation (16).
[0112] L’estimée de y(f) est calculée à l’aide de la relation (18) pour un ensemble d’hypothèses de position pk et de coefficient de pertes L.
[0113] Par exemple, les positions Pk varient sur l’ensemble des positions de la grille G superposée à la zone à analyser. Ces positions peuvent être définies en deux ou en trois dimensions.
[0114] De même les valeurs du coefficient L peuvent être choisies parmi plusieurs hypothèses en tenant compte de la connaissance a priori de la composition du sol.
[0115] Ensuite, pour chaque estimée de y(/) calculée, on détermine la valeur de [0116] Ensuite on retient la valeur minimale de parmi celles calculées avec les estimées de pk et L correspondantes.
[0117] L’étape 301 peut être affinée par itérations successives en redéfinissant, à chaque itération, une nouvelle grille G2 de pas plus fin autour de l’estimée de la position pk obtenue à l’itération précédente tel qu’illustré à la figure 4.
[0118] A l’issue de l’étape 301 on obtient une (ou plusieurs) estimée de qui est ensuite utilisée pour réaliser un test de détection via l’équation (13).
[0119] Pour cela, dans une étape 302, on calcule la vraisemblance de la matrice de canal reçue sous l’hypothèse H-i. Ce calcul est en fait directement donné par la valeur de l’estimée de
[0120] Dans une étape 303 on calcule la vraisemblance de la matrice de canal reçue sous l’hypothèse H0. Ce calcul est réalisé en calculant le terme qui correspond à la somme ou à la moyenne, sur l’ensemble des fréquences, du module au carré de la matrice de canal. Il correspond à l’énergie en fréquence de la matrice de canal. L’opérateur de norme d’une matrice complexe est défini par la trace du produit de la matrice et de la matrice transposée conjuguée :
[0121] où h est l’opérateur Hermitien et trace() désigne la trace d’une matrice.
[0122] Enfin, dans une étape 304, on réalise le test de vraisemblance de l’équation (13) en comparant la différence entre les termes respectivement calculés aux étapes 302 et 303, à un seuil de détection l.
[0123] En variante, la vitesse v de propagation de l’onde dans les milieux traversés, peut être estimée par une méthode externe ou bien peut être considérée comme une variable supplémentaire du problème de minimisation.
[0124] La figure 5 schématise une variante de réalisation de l’invention adaptée à la détection de multiples cibles.
[0125] Dans cette variante, le procédé de détection de cible décrit aux figures 2 et 3 est appliqué pour rechercher une seule cible (étape 501), c’est-à-dire la cible ayant l’écho prépondérant.
[0126] Ensuite, dans une étape 502, on applique un algorithme d’annulation d’écho aux signaux reçus (ou directement à la matrice de canal Hr(f)) pour annuler la contribution des réflexions du signal sur la cible détectée à l’étape 501. Un algorithme possible est l’algorithme connu sous l’acronyme CLEAN par exemple décrit dans la référence [3].
[0127] Les itérations sont arrêtées via un critère d’arrêt 503 qui peut consister en un nombre d’itérations fixe ou en un critère dépendant de la surface équivalente radar ou de l’énergie des échos des cibles détectées. Le critère peut aussi être défini à partir d’une différence d’énergie ou de surface équivalente radar entre deux cibles détectées à deux itérations successives.
[0128] Dans une variante de réalisation, l’image de la zone analysée est fournie en entrée d’un algorithme d’optimisation 504 pour améliorer la précision de l’image. L’algorithme d’optimisation est, par exemple, choisi parmi ceux décrits dans les références [7], [8] ou [9]. L’algorithme peut être basé sur la méthode du gradient stochastique.
[0129] La méthode décrite à la figure 5 basée sur l’annulation d’échos par itérations successives présente cependant certains inconvénients.
[0130] Tout d’abord, l’annulation de la contribution d’une cible n’est jamais parfaite et introduit donc des effets parasites. C’est particulièrement vrai lorsque plusieurs échos sont reçus avec des dynamiques très différentes, par exemple parce que les différentes cibles sont éloignées les unes des autres.
[0131] Ensuite, l’hypothèse de base pour la détection d’une cible via un test sous deux hypothèses tel que décrit précédemment, est que la statistique du bruit correspond à un bruit blanc additif Gaussien. Dans le cas de plusieurs cibles, cette hypothèse n’est plus vraie ce qui amène à des biais dans l’algorithme de détection.
[0132] Pour pallier ces inconvénients, on propose une variante de réalisation de l’invention décrite aux figures 6 et 7.
[0133] Cette variante consiste à appliquer un filtre sous forme de fenêtrage temporel aux matrices de canal afin de traiter séparément différentes sous-zones de la zone du sol à analyser et ainsi de s’affranchir des problèmes de dynamique.
[0134] Les étapes supplémentaires mises en oeuvre dans cette variante sont illustrées à la figure 7. Si on note h un vecteur constitué de l’estimée du canal de propagation en fréquence pour un couple d’antennes d’émission et de réception, ce vecteur est converti dans le domaine temporel via une transformée de Fourier inverse IFFT (étape 701). Ensuite, un filtre est appliqué dans le domaine temporel (étape 702) à la réponse temporelle du canal. Le filtre w modélise une fenêtre temporelle, il est par exemple de la forme w=[0,0,0, ...0,1 ,1 ,1 ,0,0,0] , le nombre de 1 consécutifs et leur position dans le filtre permettant de définir une tranche de la zone du sol à analyser en profondeur (dans la direction principale de de propagation de l’onde). Sans sortir du cadre de l’invention, tout type de fenêtre peut être utilisée, par exemple une fenêtre de Flanning, ou encore une fenêtre de Blackman-Flarris.
[0135] Après filtrage, le résultat est converti (étape 703) dans le domaine fréquentiel via une transformée de Fourier directe puis le nouveau vecteur fréquentiel obtenu est utilisé pour constituer la matrice de canal Hr(f) en entrée du procédé de détection selon l’une des variantes décrites précédemment. Le processus décrit à la figure 7 est appliqué pour chaque canal de propagation reliant l’une des antennes d’émission à l’une des antennes de réception.
[0136] Le procédé de détection, dans son intégralité, est itéré pour plusieurs fenêtres temporelles différentes, par exemple en faisant varier les fenêtres sur l’intégralité de la zone à analyser en permutant circulairement à chaque itération le filtre w de p/2 positions où p est égal au nombre de 1.
[0137] La figure 6 schématise, sur un diagramme, un exemple de réponse temporelle 601 d’un canal de propagation entre une antenne d’émission et une antenne de réception. Sur le même diagramme, on a représenté plusieurs fenêtres temporelles disposées en quinconce et correspondant à plusieurs itérations du procédé décrit ci- dessus.
[0138] L’application de fenêtres temporelles permet de s’affranchir des problèmes de dynamique et d’analyser le sol par couches successives. Cette méthode permet de détecter des cibles dont l’écho est faible (par exemple du fait d’une cible profonde) dans une scène avec la présence d’échos forts espacés en profondeur.
[0139] La méthode proposée permet, en plus de détecter la position d’une cible, d’estimer les pertes de propagation dans le milieu du sol et la surface équivalente radar de la cible. Ces informations peuvent servir à caractériser les sols ou la nature de la cible. [0140] La figure 8 schématise un exemple de radar à pénétration de sol configuré pour mettre en oeuvre l’invention.
[0141] Le radar est constitué d’un émetteur EM et d’un récepteur REC, chacun pourvu d’au moins une antenne A1 ,A2.
[0142] L’émetteur EM comporte un générateur d’une séquence de bits ou symboles GEN, un modulateur OFDM MOD, un convertisseur de fréquence RF1 pour transposer le signal vers une fréquence porteuse et un amplificateur AMP.
[0143] Le récepteur REC comporte un convertisseur de fréquence RF2 pour transposer le signal reçu en bande de base, un démodulateur OFDM DEMOD et une unité de traitement constituée d’une module d’estimation de canal EST qui détermine les matrices de canal à partir des signaux démodulés et de la séquence émise par le générateur GEN, un intégrateur INT qui réalise une intégration (moyenne) et une concaténation des sous-bandes fréquentielles et un module de détection de cible DET qui met en oeuvre les étapes de traitement décrites précédemment pour réaliser l’un quelconque des modes de réalisation de l’invention.
[0144] Sans sortir du cadre de l’invention, le schéma de la figure 8 est donné à titre purement illustratif et le radar peut comporter d’autres équipements non identifiés à la figure 8.
[0145] Les modules d’estimation de canal EST, d’intégration INT et de détection DET peuvent être réalisés sous forme logicielle et/ou matérielle notamment en utilisant un ou plusieurs processeur(s) et une ou plusieurs mémoire(s). Le processeur peut être un processeur générique, un processeur spécifique, un circuit intégré propre à une application (connu aussi sous le nom anglais d’ASIC pour « Application-Specific Integrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour « Field-Programmable Gâte Array »).
[01461 Références
[0147] [1 ] M. Ambrosanio, M. T. Bevacqua, T. Isernia, and V. Pascazio. The tomographie approach to ground-penetrating radar for underground exploration and monitoring : A more user-friendly and unconventional method for subsurface investigation. IEEE Signal Processing Magazine. [0148] [2] Chin Yuan Chong. Signal Processing for MIMO radars : détection under gaussian and non-gaussian environments and application to STAP. In PhD thesis, 2011.
[0149] [3] J. A. Hogbom. Aperture Synthesis with a Non-Regular Distribution of Interferometer Baselines. AAPS, 15 :417, Jun 1974.
[0150] [4] D. H. N. Marpaung and Yilong Lu. A comparative study of migration algorithms for UWB GPR images in SISO-SAR and MIMO-array configurations. In 2014 15th International Radar Symposium (1RS), pages1-4, 2014.
[0151] [5] Z. Zeng, J. Li, L. Huang, and F. Liu. Improvement of target imaging quality by multi-polarization MIMO GPR. In 2012 14th International Conférence on Ground Penetrating Radar (GPR), pages 119-124, 2012.
[0152] [6] X. Zhuge and A. G. Yarovoy. Three-dimensional near-field mimo array imaging using range migration techniques. IEEE Transactions on Image Processing, 21 (6) :3026-3033, 2012.
[0153] [7] Byrd, R. H., J. C. Gilbert, and J. Nocedal. “A Trust Région Method Based on Interior Point Techniques for Nonlinear Programming.” Mathematical Programming, Vol 89, No. 1 , 2000, pp. 149-185.
[0154] [8] Byrd, R. H., Mary E. Hribar, and Jorge Nocedal. “An Interior Point Algorithm for Large-Scale Nonlinear Programming.” SIAM Journal on Optimization, Vol 9, No. 4, 1999, pp. 877-900.
[0155] [9] Coleman, T. F. and Y. Li. “An Interior, Trust Région Approach for Nonlinear Minimization Subject to Bounds SIAM Journal on Optimization, Vol. 6, 1996, pp. 418-445.

Claims

REVENDICATIONS
1. Procédé de détection d’au moins une cible enterrée dans une zone du sol, à l’aide d’un radar à pénétration de sol, le radar comprenant au moins une antenne d’émission et au moins une antenne de réception, le procédé comprenant les étapes de :
- Acquérir (201), sur chaque antenne de réception, une mesure d’un signal émis par chaque antenne d’émission et réfléchi dans la zone du sol,
- Déterminer (202), à partir desdites mesures, une estimée, dans le domaine fréquentiel, du canal de propagation pour chaque couple constitué d’une antenne d’émission et d’une antenne de réception, lesdites estimées étant regroupées dans une matrice de canal ayant une première dimension égale au nombre d’antennes de réception, une seconde dimension égale au nombre d’antennes d’émission et une troisième dimension égale au nombre de fréquences,
- définir une première hypothèse H0 correspondant à une absence de cible dans la zone du sol, la matrice de canal sous cette première hypothèse étant modélisée par du bruit,
- définir une seconde hypothèse H1 correspondant à la présence d’au moins une cible dans la zone du sol, la matrice de canal sous cette seconde hypothèse étant modélisée par un modèle de canal de propagation dépendant au moins des caractéristiques des antennes d’émission et des antennes de réception et des pertes de propagation du signal dans le sol,
- réaliser (203) un test du rapport de vraisemblance entre la vraisemblance de la matrice de canal sous la seconde hypothèse H1 et la vraisemblance de la matrice de canal sous la première hypothèse H0, pour conclure à la présence d’une cible en une position donnée dans la zone du sol.
2. Procédé de détection selon la revendication 1 dans lequel l’étape de réaliser un test du rapport de vraisemblance entre la première hypothèse H0 et la seconde hypothèse H1 comprend les sous-étapes de :
- Déterminer (301) au moins un estimateur du maximum de vraisemblance de la position et de la surface équivalente radar d’une cible, sous la seconde hypothèse H1 ,
- Déterminer (302) la vraisemblance de la matrice de canal sous la seconde hypothèse H1 pour chaque estimateur,
- Déterminer (303) la vraisemblance de la matrice de canal sous la première hypothèse H0,
- comparer (304) le ratio entre les deux vraisemblances à un seuil de détection prédéterminé pour en déduire la présence ou l’absence d’une cible pour chaque position estimée.
3. Procédé de détection selon la revendication 2 dans lequel l’étape de déterminer (301) au moins un estimateur du maximum de vraisemblance de la position et de la surface équivalente radar d’une cible, sous la seconde hypothèse H1 est réalisée en recherchant au moins un extremum d’une fonction de trois variables dépendant de la position d’une cible dans la zone du sol, de la surface équivalente radar d’une cible et d’un coefficient de pertes de transmission du signal dans la zone du sol.
4. Procédé de détection selon la revendication 3 dans lequel la recherche d’au moins un extremum de ladite fonction de trois variables est réalisée en :
- déterminant un estimateur du maximum de vraisemblance de la surface équivalente radar en calculant une corrélation normalisée entre la matrice de canal et une matrice de pertes incluant les fonctions de transfert des antennes d’émission et de réception, les pertes de propagation du signal dans l’espace libre et dans la zone du sol,
- recherchant au moins un extremum selon la position et le coefficient de pertes de ladite fonction évaluée pour la surface équivalente radar estimée.
5. Procédé de détection selon la revendication 4 dans lequel la vraisemblance de la matrice de canal sous la seconde hypothèse H1 est déterminée par un extremum de ladite fonction.
6. Procédé de détection selon l’une des revendications 2 à 5 dans lequel la vraisemblance de la matrice de canal sous la première hypothèse H0 est déterminée au moins en calculant l’énergie de la matrice de canal.
7. Procédé de détection selon l’une quelconque des revendications précédentes dans lequel une position de cible est recherchée dans une grille (G1 ,G2) en deux ou trois dimensions ayant un pas initial défini, le procédé étant itéré sur plusieurs itérations en réduisant la zone et le pas initial à chaque itération.
8. Procédé de détection selon l’une quelconque des revendications précédentes comprenant les étapes de :
- Convertir (701) les estimées de canal de propagation pour chaque couple constitué d’une antenne d’émission et d’une antenne de réception, dans le domaine temporel,
- Appliquer (702) un filtre aux estimées sous la forme d’un fenêtrage temporel,
- Convertir (703) les estimées filtrées dans le domaine fréquentiel,
- Regrouper les estimées filtrées dans une matrice de canal
9. Procédé de détection selon la revendication 8 dans lequel le fenêtrage temporel est défini de manière à filtrer les contributions des signaux correspondant à un intervalle de distance entre le radar et la zone du sol.
10. Procédé de détection selon la revendication 9, le procédé étant répété pour plusieurs filtres différents de manière à couvrir toute la zone du sol.
11. Procédé de détection selon l’une quelconque des revendications précédentes comprenant en outre une étape d’annulation (501) de la contribution d’une cible détectée dans les mesures de signaux acquis, par exemple au moyen d’un algorithme d’annulation d’écho, le procédé étant itéré sur plusieurs itérations successives, en annulant à chaque nouvelle itération l’écho de la cible détectée à l’itération précédente.
12. Radar à pénétration de sol comprenant au moins une antenne d’émission et au moins une antenne de réception et un dispositif de détection d’au moins une cible enterrée dans une zone du sol configuré pour exécuter les étapes du procédé de détection selon l’une quelconque des revendications précédentes.
13. Programme d'ordinateur comportant des instructions pour l'exécution du procédé selon l’une quelconque des revendications 1 à 11 , lorsque le programme est exécuté par un processeur.
14. Support d'enregistrement lisible par un processeur sur lequel est enregistré un programme comportant des instructions pour l'exécution du procédé selon l’une quelconque des revendications 1 à 11 , lorsque le programme est exécuté par un processeur.
EP21733471.3A 2020-06-30 2021-06-22 Procede de detection de cible pour radar a penetration de sol et radar associe Pending EP4172655A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR2006870A FR3111994B1 (fr) 2020-06-30 2020-06-30 Procédé de détection de cible pour radar à pénétration de sol et radar associé
PCT/EP2021/067038 WO2022002700A1 (fr) 2020-06-30 2021-06-22 Procede de detection de cible pour radar a penetration de sol et radar associe

Publications (1)

Publication Number Publication Date
EP4172655A1 true EP4172655A1 (fr) 2023-05-03

Family

ID=73013559

Family Applications (1)

Application Number Title Priority Date Filing Date
EP21733471.3A Pending EP4172655A1 (fr) 2020-06-30 2021-06-22 Procede de detection de cible pour radar a penetration de sol et radar associe

Country Status (4)

Country Link
US (1) US12392889B2 (fr)
EP (1) EP4172655A1 (fr)
FR (1) FR3111994B1 (fr)
WO (1) WO2022002700A1 (fr)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114488062B (zh) * 2022-01-18 2025-03-04 中国铁道科学研究院集团有限公司 一种道床脏污状态的确定方法、装置、设备和存储介质
FR3134190B1 (fr) 2022-03-29 2024-05-03 Commissariat Energie Atomique Procédé de détection de structures longitudinales enterrées au moyen d’un radar à pénétration de sol
FR3134191B1 (fr) 2022-03-29 2024-09-06 Commissariat Energie Atomique Procédé de détection de structures longitudinales enterrées au moyen d’un radar à pénétration de sol
CN115517638B (zh) * 2022-08-01 2025-04-01 电子科技大学长三角研究院(衢州) 基于分布式mimo雷达的异常心率检测方法
CN118282440A (zh) * 2022-12-30 2024-07-02 北京三星通信技术研究有限公司 一种通信系统中的节点设备及其方法
CN115877385B (zh) * 2023-03-02 2023-05-09 中国电子科技集团公司信息科学研究院 基于无人飞艇平台的分布式雷达探测系统目标检测方法
CN116660889A (zh) * 2023-04-19 2023-08-29 杭州腓腓科技有限公司 多基地雷达系统及其探测方法、装置和存储介质
CN116595862B (zh) * 2023-04-20 2024-02-13 西安电子科技大学 基于高斯过程回归的自适应建模方法
CN117233750B (zh) * 2023-08-28 2026-04-24 中煤浙江勘测设计有限公司 一种探地雷达监测地质的方法、系统及存储介质
CN118837846A (zh) * 2024-07-08 2024-10-25 中国科学院声学研究所 一种目标位置确定方法、装置、介质、设备和程序产品
CN118604776B (zh) * 2024-08-09 2024-10-29 北京中建建筑科学研究院有限公司 一种基于分层遍历能量最大的探地雷达C-scan切片深度自动确定方法
FR3166979A1 (fr) * 2024-10-02 2026-04-03 Commissariat à l'Energie Atomique et aux Energies Alternatives Méthode de détermination de propriétés radioélectriques d’un milieu stratifié au moyen d’un système de détection radiofréquences
FR3166980A1 (fr) * 2024-10-02 2026-04-03 Commissariat à l'Energie Atomique et aux Energies Alternatives Méthode de détermination de propriétés radioélectriques d’un milieu non-homogène au moyen d’un système de détection radiofréquence
FR3166978A1 (fr) * 2024-10-02 2026-04-03 Commissariat à l'Energie Atomique et aux Energies Alternatives Méthode de détermination de propriétés radioélectriques d’un milieu non-homogène au moyen d’un système de détection radiofréquence
CN119717060B (zh) * 2024-12-05 2025-12-09 中国人民解放军海军工程大学 一种水下电磁探测系统最佳接收机推导方法及系统
CN121805991A (zh) * 2026-03-11 2026-04-07 自然资源部珠海海洋中心(自然资源部珠海海洋预报台) 基于无人机ris辅助的地波雷达定位方法、系统及装置

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4783753A (en) * 1986-05-07 1988-11-08 Environmental Research Inst. Of Michigan System for reducing speckle noise
US4839654A (en) * 1987-03-11 1989-06-13 The Tokyo Electric Power Co., Inc. System for detecting underground objects
US5327139A (en) * 1992-09-11 1994-07-05 The Boeing Company ID microwave holographic sensor
FR2751087B1 (fr) * 1996-07-09 1998-11-06 Thomson Csf Procede et dispositif de detection de cibles pour radar doppler a impulsions non ambigu a large bande
US6700526B2 (en) * 2000-09-08 2004-03-02 Witten Technologies Inc. Method and apparatus for identifying buried objects using ground penetrating radar
US6466155B2 (en) * 2001-03-30 2002-10-15 Ensco, Inc. Method and apparatus for detecting a moving object through a barrier
US7173560B2 (en) * 2003-08-28 2007-02-06 University Of Florida Research Foundation, Inc. Land mine detector
US7834801B2 (en) * 2003-11-25 2010-11-16 Metrotech Corporation, Inc. Sensor fusion for model-based detection in pipe and cable locator systems
US7382312B2 (en) * 2004-10-12 2008-06-03 George Mason Intellectual Properties, Inc. Syntactic target classifier
US7511654B1 (en) * 2006-01-12 2009-03-31 The United States Of America As Represented By The Secretary Of The Army Systems and methods for mine detection
US9075129B2 (en) * 2008-12-10 2015-07-07 The United States Of America As Represented By The Secretary Of The Army Method and system for forming images by comparing subsets of image data
US8730085B2 (en) * 2010-08-26 2014-05-20 Lawrence Livermore National Security, Llc Spot restoration for GPR image post-processing
US8970429B2 (en) * 2012-06-14 2015-03-03 Raytheon Company Systems and methods for tracking targets by a through-the-wall radar using multiple hypothesis tracking
US9291710B2 (en) * 2012-10-31 2016-03-22 Board Of Regents, The University Of Texas System Method and apparatus for detecting subsurface targets using data inversion and a temporal transmission line model
EP2972252B1 (fr) * 2013-03-14 2023-10-18 Robert Ernest Troxler Systèmes et procédés permettant de mesurer les propriétés d'ensemble d'une chaussée et l'humidité d'un sol au moyen d'un géoradar
US10175350B1 (en) * 2018-06-28 2019-01-08 University Of South Florida Systems and methods for detecting buried objects

Also Published As

Publication number Publication date
US20230266461A1 (en) 2023-08-24
FR3111994B1 (fr) 2023-02-24
US12392889B2 (en) 2025-08-19
FR3111994A1 (fr) 2021-12-31
WO2022002700A1 (fr) 2022-01-06

Similar Documents

Publication Publication Date Title
WO2022002700A1 (fr) Procede de detection de cible pour radar a penetration de sol et radar associe
EP2100161B1 (fr) Procede de traitement radar passif multivoies d'un signal d'opportunite en fm
EP0665665B1 (fr) Procédé et dispositif permettant à un modem de se synchroniser sur un transmetteur de données numériques par voie hertzienne en présence de brouilleurs
EP2579063B1 (fr) Procédé et système de localisation d'interferences par sous-bande fréquentielle
EP2523020A1 (fr) Procédé et système de localisation d'interférences affectant un signal de radionavigation par satellite
EP1582888B1 (fr) Procédé de localisation aveugle large bande d'un ou plusieurs émetteurs à partir d'un porteur défilant
WO2015177172A1 (fr) Contraintes conjointes de transitivite de differences temporelles et effet dopler multibandes pour la separation, caracterisation, et localisation de sources sonores par acoustique passive
EP4375663A1 (fr) Méthode et dispositif pour déterminer la permittivité relative d'un matériau à l'aide d'un radar à pénétration de sol
EP4359818A1 (fr) Méthode d'imagerie géoradar et géoradar associé
EP3356840B1 (fr) Procédé de détermination de caractéristiques d'une onde électromagnétique
EP2851703A1 (fr) Procédé réalisant conjointement la synchronisation, l'identification, la mesure, l'estimation du filtre de propagation et la localisation d'émetteurs utiles et interferants
EP3036840B1 (fr) Procédé de détection d'un signal électromagnétique par un réseau antennaire et dispositif mettant en oeuvre le procédé
EP3236280B1 (fr) Calcul automatique d'une dimension d'une plateforme mobile
EP1229696B1 (fr) Estimation des paramètres de canal au moyen d'estimation de maximum de vraisemblance
EP4082122B1 (fr) Procédé de détection de la présence d'un signal de brouillage et dispositifs associés
EP3400457B1 (fr) Procédé d'élimination d'un signal en provenance d'un radar de bord
Shin et al. Distributed blind deconvolution of seismic signals under sparsity constraints in sensor networks
EP4722763A1 (fr) Methode de determination de proprietes radioelectriques d'un milieu non-homogene au moyen d'un systeme de detection radiofrequence
EP4254008A1 (fr) Procede de detection de structures longitudinales enterrees au moyen d'un radar a penetration de sol
EP4722765A1 (fr) Methode de determination de proprietes radioelectriques d'un milieu non-homogene au moyen d'un systeme de detection radiofrequence
EP4722764A1 (fr) Methode de determination de proprietes radioelectriques d'un milieu stratifie au moyen d'un systeme de detection radiofrequences
CA2586008A1 (fr) Procede pour caracteriser les emetteurs par association de parametres lies a un meme emetteur radio-electrique

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20221222

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
RAP3 Party data changed (applicant data changed or rights of an application transferred)

Owner name: COMMISSARIAT A L'ENERGIE ATOMIQUE ET AUX ENERGIESALTERNATIVES

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20250307