FR 3061309 A1 20180629 FR 1601875 A 20161227 1 Méthode de traitement d'un signal issu d'un lidar cohérent pour réduire le bruit et système lidar associé DOMAINE DE L'INVENTION L'invention concerne le domaine des lidars cohérents. ETAT DE LA TECHNIQUE Le principe d'un lidar cohérent est bien connu de l'état de la technique et illustré figure 1. Un lidar cohérent comprend une source cohérente L, typiquement un laser, qui émet une onde lumineuse cohérente (domaine IR, visible ou UV proche), un dispositif d'émission DE qui permet d'illuminer un volume de l'espace, et un dispositif de réception DR, qui collecte une fraction de l'onde lumineuse rétrodiffusée par une cible T. Le décalage de fréquence Doppler VIDop de l'onde rétrodiffusée est fonction de la vitesse radiale v de la cible T : VDop = 2v/À À longueur d'onde du laser A la réception on effectue un mélange entre l'onde lumineuse rétrodiffusée reçue Sig de fréquence signal fs et une partie de l'onde émise SOL dénommée onde OL pour « oscillateur local » qui présente une fréquence d'oscillateur local foL. L'interférence de ces deux ondes est détectée par un photodétecteur D, et le signal électrique en sortie du détecteur présente un terme oscillant dénommé signal de battement Sb, en plus des termes proportionnels à la puissance reçue et à la puissance d'oscillateur local. Une unité de traitement UT numérise ce signal et en extrait une information de localisation de la cible T. Cette information de localisation de la cible est préférentiellement une information de vitesse y, et/ou une information de position (avec un signal émis particulier décrit plus loin), voire une information de présence et/ou de vibration. Préférentiellement l'unité de traitement UT filtre électroniquement le signal de 30 battement Sb dans une bande étroite centrée sur la fréquence nulle. Pour les lidars cohérents, les dispositifs d'émission et de réception utilisent préférentiellement une même optique (lidar monostatique), tel qu'illustré 3061309 2 figure 2. Cette caractéristique permet d'obtenir une bonne stabilité mécanique et de réduire l'influence des turbulences atmosphériques à longue distance, les trajets de propagation des ondes incidente et rétrodiffusée étant confondus. 5 Le signal optique à l'émission est amplifié par un amplificateur EDFA, puis transporté pour l'émission par une fibre optique monomode FM. Les voies d'émission et de réception utilisent une même optique 0 et sont séparées à l'aide d'un circulateur C. Ce signal optique peut éventuellement être décalé en fréquence, par exemple à l'aide d'un modulateur acousto-optique qui est 10 préférentiellement positionné avant l'amplificateur EDFA mais peut également être positionné sur le trajet de l'oscillateur local. Une ligne à retard LR permet d'égaliser les chemins optiques de l'oscillateur local et du signal d'émission de manière à filtrer, dans le domaine RF, les défauts des composants optiques placés après l'amplificateur EDFA (défaut de cross-talk 15 du circulateur C, imperfections des traitements anti-reflet de l'optique d'émission/réception 0,...). Le signal de battement présente une composante d'intérêt S et un bruit B: Sb = S+B. 20 Une fois détecté, il est numérisé selon une fréquence d'échantillonnage. Afin de remonter à l'information de vitesse, on calcule, sur un intervalle de temps donné, une transformée dans le domaine fréquentiel du signal de battement numérisé, typiquement une transformée de Fourier, puis on détermine la densité spectrale de puissance DSP, correspondant à la norme 25 (module au carré) de cette transformée. DSP = TF[S + B]12 Le signal d'intérêt est constitué d'un pic P (illustré figure 3). La densité DSP présente également un « plancher » de bruit qui correspond, pour un lidar dimensionné correctement, au bruit de photon (bruit relié à la statistique d'arrivée des photons détectés, également dénommé bruit de grenaille lié à la puissance de l'oscillateur local. Le signal issu de la cible rétrodiffusée étant très faible devant le signal d'oscillateur local, seul le bruit de photon de ce dernier intervient Le bruit de photon Bph s'exprime comme : 3061309 3 Bph=2.e.q.Popt.BA Avec e charge de l'électron, ri rendement du détecteur, Popt puissance optique incidente sur le détecteur et BA bande d'analyse de la densité spectrale de 5 puissance, typiquement comprise entre 0.1 Hz et 100 MHz en fonction du type de détecteur. Ce bruit lorsqu'exprimé dans le domaine spectral, présente une variation aléatoire et : 10 -une ligne de base correspondant à une valeur moyenne ou espérance BO, -une variance V qui limite la sensibilité du lidar (liée à la capacité du système à détecter le ou les pics de fréquence). Pour une meilleure lisibilité compte tenu de la dynamique accessible au lidar, 15 on calcule une densité spectrale DSPN normalisée par la moyenne de la densité spectrale de puissance de bruit, et on la représente selon une échelle logarithmique : DSPN =1TF[S + ..=1TF[S + BO 20 On définit classiquement deux grandeurs caractéristiques du bruit : -Le SNR ou Signal Noise Ratio défini comme la puissance du pic P sur la variance du bruit V: SNR = P/V 25 -Le contraste C du pic défini comme la valeur du pic sur la moyenne du bruit C = P/BO Afin de réduire le bruit, on réalise n calculs de DSP (ou de DSPN) sur n intervalles de temps ti et on effectue une moyenne. Du fait du caractère 30 aléatoire du bruit, il est connu que le fait de calculer une moyenne sur n valeurs de DSP permet de réduire la variance du bruit V, mais pas son espérance (valeur moyenne). L'évolution de la densité spectrale normalisée DSPN en fonction du nombre de valeurs n utilisées pour calculer la moyennes est illustré figure 3. Les 3061309 4 figures 3a, 3b, 3c 3d et 3e correspondent respectivement à n=1, n=2, n=5, n=10 et n=20. Le pic P sort progressivement du bruit. La variance V du bruit diminue également, mais pas sa valeur moyenne. Donc le SNR est proportionnel à n alors que le contraste reste constant. 5 Avec le traitement selon l'état de la technique décrit précédemment, lorsque n augmente, la variance V diminue, et donc le SNR augmente, ce qui améliore la performance de détection, tandis que l'espérance, indépendante de n, ne diminue pas, et le contraste n'est pas amélioré par le traitement. 10 En effet les échantillons d'une transformée de Fourier sont décrits par une variable aléatoire complexe centrée. Un périodogramme (densité spectrale de puissance élémentaire) présente donc une statistique exponentielle décroissante (une loi du x2 d'ordre 2) et l'accumulation de n périodogrammes indépendants est donc décrite par une distribution x2 d'ordre 2n (de fonction 15 de répartition la fonction Gamma incomplète d'ordre 2n). L'espérance de la distribution x2 d'ordre 2n est 2n et sa variance est de 4n. En moyennant en puissance les n périodogrammes (on divise l'accumulation par n et donc chaque densité spectrale de puissance par n), on divise la variance du bruit par n2. La variance de la moyenne en puissance de n 20 densités spectrales de puissance DSPi est donc proportionnelle à 1/n (ou l'écart type est proportionnel à 1/sqrt(n)) mais son espérance est indépendante de n. Une solution pour limiter le bruit de photon est de limiter la puissance de 25 l'oscillateur local mais cela conduit également à réduire la sensibilité du lidar, car réduire la puissance réduit également l'intensité des pics. Et le rapport signal sur bruit est indépendant de la puissance de l'oscillateur local puisque la puissance du signal de battement est elle-même proportionnelle à la puissance de l'oscillateur local. 30 A partir de la densité spectrale de puissance finale, on détermine la ou les fréquences correspondant au(x) pic(s), et on en déduit de manière classique une information de vitesse axiale v de la cible, et éventuellement une information de distance d (voir plus loin modulation de fréquence). 35 3061309 5 On connait de l'état de la technique l'utilisation d'une détection dite équilibrée ou « balancée » («balanced» en anglais) pour supprimer le bruit d'intensité de la source cohérente BOL. Ce bruit provient de l'émission spontanée et des perturbations extérieures de la cavité laser. La variance de ce bruit est 5 proportionnelle au carré de la puissance de l'oscillateur local POL : Il OBLVLfi '°PLdf Pour cela, on utilise un détecteur D équilibré ou « balancé », constitué de deux détecteurs PD1 et PD2 tel qu'illustré figure 4. Les signaux Soi_ et Sig 10 sont répartis entre les deux détecteurs de la manière décrite sur la figure 4, et on effectue la différence des intensités reçues par les deux détecteurs : ppi= 71-4(Sig + S B oL)2 114(Sig 2 ± SoL2 BoL2 2Sig.S + 2Sig.B + 2.S 6,L,.BOL) 15 I pp 2 = (Sig - S - B 01)2 (Sig 2 ± S 2 ± BoL - 2Sig.S 01, - 2Sig.B 01, + 2.S 0L.B 0L) 4 4 Avec ri rendement du détecteur. On calcule l'intensité I : 20 I = Ipp1-Ipp2= ri (Sig. SOL + Sig. BOL) On réduit par ce montage l'influence du bruit d'intensité du laser. Celui-ci restant inférieur au signal SOL associé à l'oscillateur local, le produit Sig.SoL peut être négligé devant Sig.SoL. On retrouve donc le battement Sig.SoL. Ce type de détecteur balancé nécessite un équilibrage précis entre les 25 sensibilités de deux photodiodes et le taux de couplage du coupleur permettant de diviser l'intensité en deux voies. Il peut être utilisé pour tout type de lidar, et plus particulièrement pour tout type de lidar cohérent, mais il ne réduit que le bruit d'intensité de la source et pas le bruit de photon issu de la détection de l'oscillateur local. 30 On connait également de l'état de l'art une solution pour la télémétrie/vélocimétrie par lidar qui consiste à réaliser un système lidar à modulation de fréquence. Cette technique, classique en radar, trouve un intérêt particulier actuellement compte tenu des progrès des sources lasers 3061309 6 fibrées. Grâce à la modulation de fréquence, une analyse temps/fréquence permet de remonter à la distance d (le retard -ro du signal rétrodiffusé par la cible par rapport à l'oscillateur local est fonction de la distance d) et à sa vitesse v. Ce type de lidar permet également de réaliser une fonction 5 d'anémométrie laser. Un exemple d'architecture optique d'un lidar à modulation de fréquence utilisant un détecteur balancé D (comprenant un détecteur PD1 et un détecteur PD2), est décrit figure 5. La source cohérente est modulée en fréquence de manière à ce que la fréquence de l'oscillateur local foL(t) soit 10 modulée selon une fonction prédéterminée f,,,a(t) dénommée forme d'onde dont le contrôle est assuré par le module WFC, synchronisé avec l'unité de traitement UT. La fréquence de l'oscillateur local s'exprime alors : foL(t) =f0 + filma(t) Avec f0 fréquence moyenne du laser L 15 Un exemple d'un lidar cohérent à modulation de fréquence est décrit dans le document « Lidar systems for precision navigation and safe landing on planetary bodies » Farzin Amzajerdian et al, Proc. SPIE 8192, International Symposium on Photoelectronic Detection and Imaging 2011: Laser Sensing 20 and Imaging; and Biological and Medical Applications of Photonics Sensing and Imaging, 819202 (August 19, 2011). La fréquence de l'oscillateur local foL est modulée linéairement selon deux pentes de fréquence ao et ai de manière périodique de période TFo. La figure 6 illustre la variation au cours du temps des fréquences d'oscillateur local foL(t) et de signal fs(t). Comme 25 illustré sur la figure 6a, le signal rétrodiffusé de fréquence f(t) est décalé temporellement d'un temps t du fait de la propagation jusqu'à la zone de mesure (cible T) et donc relié à la distance de la cible d, et est décalé en fréquence d'une valeur vpap du fait de l'effet Doppler par rapport à la fréquence oscillateur local foL(t). 30 Le signal de battement détecté Sb présente une composante de fréquence positive fs-foL. La figure 6b illustre l'évolution au cours du temps de fs -fou On constate que cette différence de fréquence comprend en fonction du temps deux séries de plateaux aux fréquences caractéristiques vao et Val, directement reliées à la distance D de la cible et à sa vitesse y radiale par les 35 équations : 3061309 7 2v 2a0D 2v 2a1D = et v ao En mesurant ces deux fréquences caractéristiques vao et vai du signal de battement Sb, par exemple en réalisant une transformée de Fourier de celui-5 ci, tel que décrit ci-dessus, on remonte à d et v. L'homme de l'art connait également une autre architecture de lidar cohérent illustrée figure 7 dans laquelle le signal d'oscillateur SOL est prélevé grâce à une réflexion en bout de fibre. 10 Un but de la présente invention est de pallier les inconvénients précités en proposant une méthode de traitement du signal de battement et une architecture de détection spécifique associée permettant de réduire le bruit de photons provenant de la détection du signal d'oscillateur local compatible 15 avec l'utilisation d'un détecteur balancé et pouvant s'implémenter sur tout lidar cohérent. DESCRIPTION DE L'INVENTION 20 La présente invention a pour objet une méthode de traitement d'un signal issu d'un lidar cohérent comprenant une source cohérente, la méthode comprenant les étapes consistant à : -générer un premier signal de battement et un deuxième signal de battement 25 respectivement à partir d'un premier ensemble de détection et d'un deuxième ensemble de détection, chaque signal de battement étant issu d'une interférence entre un signal d'oscillateur local issu de la source cohérente et un signal rétrodiffusé par une cible éclairée par le lidar, puis numériser ces signaux de battement, 30 -pour une pluralité de n intervalles temporels, déterminer n valeurs respectives de densité spectrale à partir d'une transformée dans le domaine fréquentiel de la corrélation croisée entre le premier et le deuxième signal de battement, -déterminer une valeur moyenne de la densité spectrale à partir desdites n 35 valeurs de densité spectrale, 3061309 8 -déterminer une information de localisation de la cible à partir de la valeur moyenne de ladite densité spectrale. Selon un mode de réalisation l'étape de détermination de n valeurs de la 5 densité spectrale comprend les sous étapes consistant à: -déterminer une première valeur d'une transformée dans le domaine fréquentiel du premier signal de battement, -déterminer une deuxième valeur du conjugué d'une transformée dans le domaine fréquentiel du deuxième signal de battement, 10 la valeur de la densité spectrale étant déterminée à partir du produit de la première et de la deuxième valeur. Selon une variante la source cohérente est modulée périodiquement en fréquence de sorte que le signal d'oscillateur local présente une fréquence 15 d'oscillateur local constituée de la somme d'une valeur moyenne et d'une fréquence de modulation issue de la modulation de la source, la fréquence de modulation étant périodique selon une période de modulation, et dans laquelle les intervalles temporels sont inférieurs ou égaux à la période de modulation, la méthode de traitement comprenant en outre une étape 20 consistant à déterminer une information de distance de la cible à partir de la valeur moyenne de la densité spectrale. Préférentiellement chaque période de modulation de la fréquence de modulation comprend n parties linéaires présentant respectivement n pentes de fréquence, n étant supérieur ou égal à 2. 25 L'invention concerne également un système lidar cohérent comprenant : -une source cohérente, -un dispositif d'émission d'un signal optique issu de la source cohérente et un dispositif de réception d'un signal rétrodiffusé par une cible illuminée par le 30 lidar, -un premier et un deuxième ensemble de détection configurés pour générer respectivement un premier et un deuxième signal de battement, chaque signal de battement étant issu d'une interférence entre un signal d'oscillateur local issu de la source cohérente et le signal rétrodiffusé par la cible, 3061309 9 -une unité de traitement configurée pour numériser le premier et le deuxième signal de battement, et configurée pour déterminer pour une pluralité de n intervalles temporels, n valeurs de densité spectrale correspondant à une transformée dans le domaine fréquentiel de la corrélation croisée entre le 5 premier et le deuxième signal de battement, -l'unité de traitement étant en outre configurée pour : *déterminer une valeur moyenne de la densité spectrale à partir des n valeurs de densité spectrale calculées, et *déterminer une information de localisation de la cible à partir de la 10 valeur moyenne de la densité spectrale. Selon un mode de réalisation le système lidar cohérent selon l'invention comprend en outre un dispositif de modulation synchronisé avec l'unité de traitement et configuré pour moduler périodiquement en fréquence la source 15 cohérente de sorte que le signal d'oscillateur local présente une fréquence d'oscillateur local constituée de la somme d'une valeur moyenne et d'une fréquence de modulation issue de la modulation de la source, la fréquence de modulation étant périodique selon une période de modulation, chaque période comprenant n parties linéaires présentant respectivement n pentes 20 de fréquence, n étant supérieur ou égal à 2. L'unité de traitement est également configurée pour que les intervalles temporels ti soient inférieurs ou égaux à la période de modulation et pour déterminer une information de distance de la cible à partir de la valeur moyenne de la densité spectrale. 25 Selon une variante préférée, le premier et/ou le deuxième ensemble de détection sont des détecteurs balancés comprenant chacun un premier et un second détecteur, les premiers détecteurs recevant une différence entre le signal d'oscillateur local et le signal rétrodiffusé, les deuxièmes détecteurs recevant une somme du signal d'oscillateur local et du signal rétrodiffusé. 30 Les premier et deuxième signaux de battement sont issus de la différence des intensités reçues par respectivement le premier et le deuxième détecteur du premier ensemble de détection et le premier et le deuxième détecteur du deuxième ensemble de détection. 3061309 10 Préférentiellement le premier et/ou le deuxième ensemble de détection sont placés de sorte que la longueur des trajets empruntés par chacun des signaux jusqu'auxdits ensembles de détection soit sensiblement égale. 5 D'autres caractéristiques, buts et avantages de la présente invention apparaîtront à la lecture de la description détaillée qui va suivre et en regard des dessins annexés donnés à titre d'exemples non limitatifs et sur lesquels : La figure 1 déjà citée décrit le principe d'un lidar cohérent selon l'état de la 10 technique. La figure 2 déjà citée illustre l'architecture d'un lidar monostatique. La figure 3 déjà citée illustre la densité de puissance du signal de battement normalisé. Les figures 3a à 3e illustrent l'évolution de celle-ci en fonction du nombre de valeurs n utilisées pour calculer la valeur moyenne. 15 La figure 4 déjà citée illustre le principe d'une détection équilibrée. La figure 5 déjà citée illustre l'architecture optique d'un lidar cohérent à modulation de fréquence. La figure 6a déjà citée illustre la variation au cours du temps des fréquences oscillateur local et signal, et la figure 6b déjà citée illustre la variation au 20 cours du temps de la composante positive du signal de battement. La figure 7 illustre une autre architecture de lidar cohérent connue de l'homme de l'art. La figure 8 illustre la méthode de traitement d'un signal issu d'un lidar cohérent selon l'invention. 25 La figure 9 représente un mode de réalisation préféré de la méthode de traitement selon l'invention. La figure 10 illustre un exemple d'évolution de la densité spectrale normalisée en fonction de n, les figures 10a, lob, 10c, 10d et 10e correspondant respectivement à n=1, n=2, n=5, n=10 et n=20. 30 La figure 11a illustre la variation au cours du temps des fréquences fodt) et fs(t). La figure llb illustre l'évolution au cours du temps de fs -foi. pour le cas à 4 pentes. La figure 12 représente un exemple de densité spectrale de puissance normalisée exprimée en dB du signal de battement détecté par un détecteur 35 unique selon l'état de la technique. 3061309 11 La figure 13 représente la densité spectrale de puissance normalisée exprimée en dB du signal de battement pour le même cas que la figure 12, le signal de battement étant détecté par deux ensembles de détection et la densité spectrale de puissance normalisée calculée avec la méthode de 5 traitement selon l'invention. La figure 14 schématise un système lidar cohérent selon l'invention. La figure 15 schématise une variante préférée du système lidar cohérent selon l'invention. La figure 16 schématise une variante préférée du système lidar cohérent 10 selon l'invention compatible d'une détection équilibrée. DESCRIPTION DETAILLEE DE L'INVENTION 15 Avant de décrire l'invention nous allons rappeler certains concepts mathématiques connus de l'homme de l'art et nécessaires à la bonne compréhension de l'invention. On définit la transformée de Fourier TF d'une fonction temporelle S(t) , 20 exemple de transformée dans le domaine fréquentiel : +00 TF[S(t)] = S(t). exp(-- 2ievt)dt On définit le produit de convolution entre deux fonctions temporelles 25 complexes S1(t) et S2(t) : S1* S =1° S i(z).S; (t - h- TF [S i(t)* S2 (t)] = TF [S l(t)]. TF [S 2(t)] S étant un nombre complexe : 30 S = Re) -F i IM(S) et s* = Re(s)- im(s) On définit également la fonction de corrélation croisée de Si et S2: 3061309 12 1-co +00 cSi S2 (t) S2 (t) = f (t ).S; )dt = f (t 1").S; )dt (1) Et sa transformée de Fourier : 5 TF[Csis2 = TF[Si (t)* (-- t)] = TF[Si (t)] x TF[S;(- t)] = TF[si (t)] x TF[S2 (t)] r (2) La fonction de corrélation est maximale quand SI et 52 sont identiques à un décalage temporel près. cs1s2 est alors maximale pour une valeur de -co correspondant à ce décalage. -Io La méthode 80 de traitement d'un signal issu d'un lidar cohérent comprenant une source cohérente L selon l'invention est schématisée figure 8. La méthode comprend une étape 100 de génération d'un premier signal de battement Sb1 à partir d'un premier ensemble de détection D1, le signal de 15 battement Sb1 étant issu d'une interférence entre un signal d'oscillateur local SOL issu de la source cohérente et un signal rétrodiffusé par une cible Sig éclairée par le lidar. Puis également dans l'étape 100 le signal Sb1 est numérisé à une fréquence d'échantillonnage. De même une étape 110 consiste à générer un deuxième signal de 20 battement Sb2 à partir d'un deuxième ensemble de détection D2, le signal de battement Sb2 étant issu d'une interférence entre un signal d'oscillateur local Soi_ issu de la source cohérente et un signal rétrodiffusé par une cible Sig éclairée par le lidar. Puis le signal Sb2 est numérisé. Pour générer Sb1 et 5b2 une fraction de chaque signal Sig et SOL est 25 dirigée vers les ensembles de détection D1 et D2 comme illustré plus loin. Puis, pour une pluralité de n intervalles temporels ti indicés i, i variant de 1 à n, une étape 120 détermine n valeurs de densité spectrale correspondantes DSPi, à partir d'une transformée dans le domaine fréquentiel de la corrélation 30 croisée entre le premier signal de battement Sb1(t) et le deuxième signal de battement Sb2(t). Pour un intervalle de temps donné ti, cette opération consiste à réaliser une transformée dans le domaine fréquentiel de la fonction de corrélation Csb1Sb2 définie par la formule (1) sur l'intervalle de temps ti considéré. 3061309 13 Typiquement on réalise une transformée de Fourier de CSblSb2 tel qu'explicité par la formule (2) sur l'intervalle de temps ti considéré, soit numériquement une FFT pour « Fast Fourier Transfomi » en anglais. 5 Selon une variante la densité spectrale calculée DSPi correspond directement à la transformée fréquentielle des signaux temporels de battement: DSPi=TF[Csis,(1-)] 10 Chaque signal de battement Sb1, Sb2 se décompose en un signal d'intérêt S1, S2, et une composante de bruit de photon B1, B2 issus de la détection correspondante D1, D2: Sb1 =S1 +B1 et Sb2 = S2 +B2 15 Selon une autre variante la densité spectrale calculée DSPi correspond à la densité spectrale normalisée par le bruit SPDNi: DSPNi =DSPi BO Puis dans une étape 130 on détermine une valeur moyenne de la 20 densité spectrale à partir des n valeurs DSPi de densité spectrale déterminées pour les intervalles de temps ti. Préférentiellement, on réalise une moyenne linéaire : (Dsp). 'DSPi n 25 Enfin dans une étape 140 on détermine une information de localisation de la cible à partir de la valeur moyenne de ladite densité spectrale. Cette information de localisation de la cible est préférentiellement une information de vitesse y, et/ou une information de position (avec un signal émis 30 particulier décrit plus loin), voir une information de simple présence et/ou de vibration (détection des modes propres de vibration d'une cible, qui ne se déplace pas nécessairement). 3061309 14 Avec la méthode 80 selon l'invention, plutôt que de calculer une norme d'une FFT à partir d'un seul signal détecté selon l'état de la technique, on réalise une transformée fréquentielle à partir de la corrélation de deux signaux de battement issus de deux détecteurs différents, ce pour une pluralité 5 d'intervalles de temps ti. Cette méthode de traitement de signal, et la double détection associée permettent de réduire le plancher de bruit par moyennage, cette réduction s'opérant à la fois sur la variance et sur la valeur moyenne. En effet, les 10 bruits sur chaque détecteur étant principalement liés au bruit de grenaille associé à la puissance de l'oscillateur local, ceux-ci sont essentiellement décorrélés d'un détecteur à l'autre puisqu'ils sont liés à la statistique d'arrivée des photons sur chacun des ensembles de détection Dl et D2. Au contraire, les signaux d'intérêt Si et S2 sont identiques. En calculant la 15 corrélation croisée des signaux issus des deux détecteurs, la contribution du bruit en valeur moyenne diminue en fonction du nombre n de densités spectrales moyennées. Plus précisément, on calcule typiquement pour sur chaque intervalle de 20 temps ti : TF[Si +Bjx{TF[S2 +B2]1* =TF[Si]x {TF[S2]}*+TF[Si]x ITFKlls +TF[Bi]x{TF[S2]1*±TF[Bi]x{TF[B2]}* TF[Si ± B1 ]>( {TF[S2 ± B2 = TF[Si 0 S2 TF[Si 0 B2 + B1 0 S2 ± 25 Si et S2 sont très fortement corrélés, car issus des même sources et égaux (à un décalage temporel près), tandis que les corrélations croisées de Si ou S2 avec 131 ou B2 ne comprennent aucun maximum significatif. Et donc après calcul de la moyenne sur n réalisations, la moyenne : 30 (s, 0 B2 ± B1 0 S2 ± B1 0 B2) va tendre vers 0 (en espérance et en variance). En effet, en moyennant sur n réalisations (n intervalles de temps ti), les produits croisés ci-dessus vont voir leur amplitude réduite. On ne récupère à la fin que la composante d'intérêt correspondant à la transformée de Fourier de la corrélation croisée entre SI 35 et S2. 3061309 15 Plus précisément, seule la partie réelle positive du produit TF(S1).TF(S2)* est d'intérêt. Les parties négatives et imaginaires font aussi partie des termes parasites qui décroissent avec le nombre de moyennes. 5 La réduction de l'espérance du bruit (en plus de sa variance) facilite la détection des pics d'intérêt car le contraste C par rapport au niveau de bruit est fortement augmenté. Une augmentation significative de la portée du lidar est ainsi obtenue. 10 Selon une variante préférée, pratiquement les opérations s'effectuent numériquement par FFT sur les intervalles de temps ti, et on a: DSPi = FFTt1(S1+B1). FFTti(52+B2)* DSPi = FFT(S1).FFT(S2)* + FFT(S1).FFT(B2)* + FFT(B1).FFT(S2)* + 15 FFT(B1).FFT(B2)* De même après une moyenne sur n réalisations seul le terme d'intérêt FFT(S1).FFT(S2)* conserve une amplitude similaire, les autres termes voient leur amplitude se réduire. 20 En effet, si l'on calcule l'espérance de la norme au carré de la moyenne des DSPi. 2 - 1 - E n2 E DSPi E DSPi Dans cette somme, les termes impliquant les bruits B1 et B2 peuvent s'écrire sous la forme suivante compte tenu du fait que les bruits B1 et B2 suivent 25 une distribution complexe centrée : , n zai exp(ivi Où ai et (1); désignent l'amplitude et la phase des termes impliquant les produits du type FFT(S1).FFT(B2)* ou FFT(B1).FFT(S2)*. Or, on peut montrer que : 2 - 3061309 16 2 4121 za i exp(iqyi )\ j=1 Ce majorant décroit clairement vers 0 en fonction de n. Selon un mode de réalisation préféré illustré figure 9, l'étape de 5 détermination 120 de n valeurs de la densité spectrale DSPi comprend, pour chaque intervalle de temps ti, une sous étape 121 consistant à déterminer une première valeur d'une transformée dans le domaine fréquentiel du premier signal de battement TFi(Sb1) et une sous étape 122 consistant à déterminer une deuxième valeur du conjugué d'une transformée dans le 10 domaine fréquentiel du deuxième signal de battement TFi(Sb2)*, la valeur de la densité spectrale étant calculée à partir du produit des première et deuxième valeurs. En effet, pour calculer DSPi, transformée dans le domaine fréquentiel de la 15 corrélation croisée entre le premier et le deuxième signal de battement, on peut soit calculer temporellement la corrélation croisée, puis réaliser la transformée fréquentielle de la fonction calculée, soit calculer la transformée fréquentielle de chaque signal de battement détecté (ou de son conjugué) et en faire le produit. Pour optimiser la vitesse de calcul, on choisit 20 préférentiellement cette dernière méthode de calcul. Selon une variante, chaque DSPi est directement égal audit produit : DSPi = TFi(Sb1).TFi(Sb2)*. 25 Selon une autre variante préférée, chaque DSPi est égal au produit des transformées de Fourier normalisées par le bruit correspondant : DSPNi = TFi(Sb1)/B1 .TFi(Sb2)*/B2. avec norme par rapport à B1 et B2 30 Cela permet de ne pas avoir à adapter le traitement de signal pour l'exploitation des DSPNi par rapport à celui utilisé en l'absence de cette invention. En particulier le calcul des distances et vitesse à partir des DSPNi est alors rigoureusement identique à celui utilisé à partir du SNR pour une architecture mono-détecteur. 3061309 17 La figure 10 illustre un exemple d'évolution de la densité spectrale normalisée DSPN en fonction de n pour un cas similaire à celui de la figure 3, mais en utilisant la méthode 80 de traitement de signal selon l'invention. 5 Les figures 3a, 3b, 3c 3d et 3e correspondent respectivement à n=1, n=2, n=5, n=10 et n=20. Le pic P se détache du bruit avec les comportements suivants : - La variance du bruit diminue toujours en 1/n (donc le SNR diminue en 1/n) 10 - La moyenne du bruit diminue en 1/n. Donc le contraste augmente en n. La méthode 80 de traitement du signal est compatible d'un lidar à modulation de fréquence. 15 Ainsi selon un mode de réalisation, la source cohérente L du lidar cohérent auquel s'applique la méthode 80 selon l'invention est modulée périodiquement en fréquence de sorte que le signal d'oscillateur local SOL présente une fréquence d'oscillateur local foL(t) constituée de la somme d'une valeur moyenne f0 et d'une fréquence de modulation fniod(t) issue de la 20 modulation de la source, la fréquence de modulation étant périodique selon une période de modulation TF0 (voir plus loin figure 11). Les intervalles temporels ti sont dans ce cas inférieurs ou égaux à la période de modulation pour conserver une même pente de fréquence pendant ti. La méthode de traitement 80 comprend en outre une étape consistant à 25 déterminer une information de distance d de la cible à partir de la valeur moyenne de la densité spectrale, selon les méthodes classiques de l'état de la technique. La méthode 80 permet d'accroitre la sensibilité de détection des pics de fréquence et donc d'accroitre la portée de l'instrument. 30 Préférentiellement, chaque période comprend n parties linéaires présentant respectivement n pentes de fréquence ai, n étant supérieur ou égal à 2. Il est connu de l'état de la technique l'utilisation d'un signal à deux pentes. Pour lever les ambiguïtés associées aux éventuels repliements, on utilise préférentiellement une forme d'onde avec 4 rampes de fréquences a0, al, a2, 3061309 18 a3. En effet, la détermination de 4 fréquences caractéristiques conduit à un système de 4 équations, avec 2 inconnues, v et d. Va ' 2v 2a,D Ceci permet d'obtenir une redondance et donc d'utiliser une des équations 5 pour lever les ambiguïtés associées aux éventuels repliements de spectres et une autre comme paramètre de confiance. Ce paramètre de confiance peut par exemple être le résidu de l'inversion entre les fréquences va, et la vitesse radiale et distance. Cette inversion peut être obtenue par une technique de type moindres carrés, éventuellement pondérée de manière 10 itérative (Iteratively reweighted least squares - IRLS). La figure 11a illustre la variation au cours du temps des fréquences foL(t) et fs(t), la fréquence optique moyenne f0 du laser L ayant été soustraite pour plus de clarté. La figure llb illustre l'évolution au cours du temps de fs -foL pour le cas à 4 pentes. On voit sur la figure 10b que cette variation de 15 fréquence au cours du temps présente 4 plateaux correspondant aux 4 fréquences caractéristiques. La figure 12 illustre la densité spectrale de puissance DSPN normalisée exprimée en dB et moyennée avec n = 300 et ti = 16.3 ps, du signal de 20 battement détecté par un détecteur unique D selon l'état de la technique, correspondant au cas d'une cible située à 2507 m et se déplaçant à 200 m/s. La fréquence frima présentant les valeurs de pentes suivantes (laser de fréquence optique f0 = 1.55 pm): a0 = 2 MHz/ps 25 G1 = -2 MHz/ps a2 = 3 MHz/ps a3 = -3 MHz/ps On détecte 4 pics PO, P1, P2 et P3 correspondant respectivement aux 4 fréquences caractéristiques vao, Val, v02 et Va3 30 Ces pics sont symétriques du fait du caractère réel du signal de battement détecté, qui engendre deux composantes lors du calcul des transformées de Fourier. 3061309 19 La figure 13 illustre la densité de puissance DSPN normalisée exprimée en dB et moyennée avec n = 300 et ti = 16.3 ps, pour le même cas que précédemment, mais calculée avec la méthode 80 selon l'invention, c'est-à-dire à partir de la transformée de Fourier de la corrélation croisée des 5 signaux de battement Si et S2 détectés par deux ensembles de détection Dl et D2. On note sur cet exemple, par comparaison de la figure 13 avec la figure 12: -une diminution des pics de signal. Cette réduction est liée à la division du 10 signal et de l'oscillateur local par un facteur 2, du fait de la nécessité d'envoyer une fraction des signaux SOL et Sig (ici la moitié) sur chaque ensemble de détection. -une réduction de la variance V, et donc de l'écart type du plancher de bruit, qui passe sur cet exemple de 0.024 à 0.005, soit également une réduction de 15 6 dB en puissance. Le rapport signal sur bruit SNR reste identique au cas de la figure 12. -une forte augmentation du contraste C de signal par rapport au plancher de bruit : dans cet exemple celui-ci passe de 5.6 dB à 28.8 dB, soit un gain de 23.2 dB. 20 Par rapport au calcul effectué selon l'état de la technique il convient, pour mettre en oeuvre la méthode selon l'invention 80, de réaliser deux FFT au lieu d'une, ce qui allonge légèrement le temps de calcul. De manière pratique le signal de battement est numérisé en temps réel au fil de l'eau et les 25 calculs de traitement de signal sont également réalisés en temps réel au fur et à mesure de l'écoulement du temps, pour une pluralité d'intervalles de temps ti. Cette nouvelle méthode de traitement 80 (associée à une nouvelle architecture de détection) permet un gain de sensibilité du lidar en particulier 30 si le plancher de bruit, et plus particulièrement sa ligne de base est mal connue, ou (ce qui est souvent le cas) si le plancher de bruit évolue temporellement (même à l'échelle d'une période de modulation TFo). La méthode selon l'invention s'applique à l'ensemble des systèmes lidars cohérents tels que ceux visant à détecter de faibles signaux, et notamment 35 pour les applications suivantes : 3061309 20 -télémétrie/vélocimétrie à longue portée, -altimétrie -anémométrie laser. 5 L'invention concerne également un système lidar cohérent 13 illustré figure 14 comprenant : -une source cohérente L, -un dispositif d'émission DE d'un signal optique issu de la source cohérente et un dispositif de réception DR d'un signal rétrodiffusé par une cible T illuminée 10 par le lidar, -un premier ensemble de détection D1 configuré pour générer un premier signal de battement Sb1 issu d'une interférence entre un signal d'oscillateur local SOL issu de la source cohérente et le signal rétrodiffusé Sig par la cible T, 15 -un deuxième ensemble de détection D2 configuré pour générer un deuxième signal de battement Sb2 issu d'une interférence entre un signal d'oscillateur local SOL issu de la source cohérente et le signal rétrodiffusé Sig par la cible T, -une unité de traitement UT configurée pour numériser le premier et le 20 deuxième signal de battement, et configurée pour mettre en oeuvre la méthode selon l'invention décrite précédemment, ainsi que préférentiellement ses variantes et modes de réalisation. Ainsi l'unité de traitement UT est configurée pour : -déterminer, pour une pluralité de n intervalles temporels ti, n valeurs de 25 densité spectrale DSPi correspondant à une transformée dans le domaine fréquentiel de la corrélation croisée entre le premier et le deuxième signal de battement, Et pour : -déterminer une valeur moyenne de la densité spectrale à partir des 30 n valeurs de densité spectrale calculées DSPi, et -déterminer une information de localisation de la cible à partir de la valeur moyenne de la densité spectrale. Il est entendu que le lidar 13 selon l'invention comprend des composants 35 optiques ou de l'optique intégrée permettant de répartir le signal d'oscillateur 3061309 21 local et le signal rétrodiffusé sur deux voies, afin de diriger une fraction de ces signaux sur une première voie vers le premier ensemble de détection et une autre fraction sur une deuxième voie vers le deuxième ensemble de détection. 5 Il est préférable de séparer les signaux en vue de les envoyer sur les ensembles de détection D1 et D2 après le mélange avec d'oscillateur local de manière à garantir un déphasage identique sur les deux voies en particulier dans le cas d'une fréquence variable au cours du temps. 10 Si aucune précaution particulière n'est prise sur l'architecture optique de distribution des signaux vers les ensembles de détection D1 et D2, les trajets effectués par les signaux SOL et Sig, à partir d'un point de référence PREF à partir duquel les signaux SOL et Sig sont recombinés et prêts à être détectés, pour arriver respectivement sur le premier et le deuxième ensemble de 15 détection, sont différents. Les signaux de battement résultants Sb1 et Sb2 sont alors décalés temporellement. Dans le cas d'une modulation de fréquence, ce décalage conduit à réduire l'amplitude des pics obtenus lors du calcul de la transformée de Fourier du produit de corrélation. Pour pallier cet inconvénient, préférentiellement le premier ensemble de 20 détection D1 et le deuxième ensemble de détection D2 sont placés de sorte que la longueur des trajets empruntés par chacun des signaux détectés SOL et Sig, à partir d'un point ou d'un plan de référence PREF, soit sensiblement égale, tel qu'illustré figure 15. Dans ce cas, les signaux de battement Sb1 et Sb2 sont temporellement égaux. 25 Sb1(t) = Sb2(t) Cette condition est remplie exactement si la différence entre les trajets parcourus par SOL et Sig jusqu'à D1 et D2 est très inférieure à la distance que parcourerait la lumière pendant un intervalle de temps égal à l'inverse de la fréquence du signal de battement. 30 Par exemple pour une fréquence de battement de 100 MHz, correspondant à une période de 10 ns, les trajets sont préférentiellement très inférieurs à 30 CM. Selon un mode de réalisation le système lidar 13 selon l'invention comprend 35 en outre un dispositif de modulation WFC synchronisé avec l'unité de 3061309 22 traitement (UT) configuré pour moduler périodiquement en fréquence la source cohérente L de sorte que le signal d'oscillateur local SOL présente une fréquence d'oscillateur local foL(t) constituée de la somme d'une valeur moyenne f0 et d'une fréquence de modulation fn,od(t) issue de la modulation 5 de la source, la fréquence de modulation étant périodique selon une période de modulation TRD, chaque période de la fréquence de modulation comprenant n parties linéaires présentant respectivement n pentes de fréquence ai, n étant supérieur ou égal à 2. Dans ce mode de réalisation l'unité de traitement est en outre configurée 10 pour que les intervalles temporels ti soient inférieurs ou égaux à la période de modulation TRD et pour déterminer une information de distance d de la cible à partir de la valeur moyenne de la densité spectrale. Un exemple d'architecture optique d'un tel lidar cohérent à modulation de fréquence 13 selon l'invention est illustré figure 16. 15 Le lidar 13 selon l'invention est compatible de l'utilisation d'une détection équilibrée pour les ensembles de détection D1 et D2. Ainsi selon un mode de réalisation également illustré figure 14 pour le cas non limitatif d'un lidar à modulation de fréquence, le premier ensemble de détection D1 et/ou le deuxième ensemble de détection D2 sont des détecteurs balancés 20 comprenant chacun un premier détecteur, PD1 pour D1, PD1' pour D2, et un second détecteur, PD2 pour D1 et PD2 pour D2 configurés pour réaliser une détection balancée. Ainsi les premiers détecteurs PD1 et PD'1 reçoivent un signal proportionnel à la différence entre le signal d'oscillateur local SOL et le signal rétrodiffusé Sig, 25 et les deuxièmes détecteurs PD2,PD'2 reçoivent une somme du signal d'oscillateur local SOL et du signal rétrodiffusé Sig. Dans cette configuration, le premier signal de battement Sb1 est issu de la différence des intensités reçues par le premier détecteur PD1 et le deuxième détecteur PD2 du premier ensemble de détection D1, et le deuxième signal 30 de battement Sb2 est issu de la différence des intensités reçues par le premier détecteur PD1' et le deuxième détecteur PD2' du deuxième ensemble de détection D2. Chacun des modules de calcul que comporte le système selon l'invention et 35 plus particulièrement l'unité de traitement UT peut être réalisé sous forme 3061309 23 logicielle et/ou matérielle. Chaque module peut notamment consister en un processeur et une mémoire. 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- 5 Specific Integrated Circuit ») ou un réseau de portes programmables in situ (connu aussi sous le nom anglais de FPGA pour « Field-Programmable Gate Array )>). L'invention concerne également un produit programme d'ordinateur 10 comprenant des instructions de code permettant d'effectuer les étapes de la méthode de traitement selon l'invention. REVENDICATIONS1. Méthode de traitement (80) d'un signal issu d'un lidar cohérent comprenant une source cohérente (L), la méthode comprenant les étapes consistant à: -générer (100, 110) un premier signal de battement (Sb1) et un deuxième signal de battement (Sb2) respectivement à partir d'un premier ensemble de détection (D1) et d'un deuxième ensemble de détection (02), chaque signal de battement étant issu d'une interférence entre un signal d'oscillateur local (SOL) issu de la source cohérente et un signal (Sig) rétrodiffusé par une cible éclairée par le lidar, puis numériser ces signaux de battement, - pour une pluralité de n intervalles temporels (ti), déterminer (120) n valeurs respectives de densité spectrale (DSPi) à partir d'une transformée dans le domaine fréquentiel de la corrélation croisée entre le premier et le deuxième signal de battement, -déterminer (130) une valeur moyenne ( ) de la densité spectrale à 15 pardi- desdites n valeurs de densité spectrale (DSPi), -déterminer (140) une information de localisation de la cible à partir de la valeur moyenne de ladite densité spectrale. 2. Méthode selon la revendication 1 dans laquelle l'étape de détermination 20 (120) de n valeurs de la densité spectrale comprend les sous étapes consistant à: -déterminer une première valeur d'une transformée dans le domaine fréquentiel du premier signal de battement (TF1(Sb1)), -déterminer une deuxième valeur du conjugué d'une transformée dans le 25 domaine fréquentiel du deuxième signal de battement (TFi(Sb2)*), la valeur de la densité spectrale étant déterminée à partir du produit de la première et de la deuxième valeur. 3. Méthode de traitement selon l'une des revendications 1 ou 2 dans laquelle 30 la source cohérente (L) est modulée périodiquement en fréquence de sorte que le signal d'oscillateur local (SOL) présente une fréquence d'oscillateur local (foL(t)) constituée de la somme d'une valeur moyenne (f0) et d'une fréquence de modulation (fmc,d(t)) issue de la modulation de la source, la fréquence de modulation étant périodique selon une période de modulation 3061309 25 (TF0), et dans laquelle les intervalles temporels (ti) sont inférieurs ou égaux à la période de modulation, la méthode de traitement comprenant en outre une étape consistant à déterminer une information de distance (d) de la cible à partir de la valeur moyenne de la densité spectrale. 5 4. Méthode de traitement selon la revendication 3 dans laquelle chaque période de modulation (TF0) de la fréquence de modulation comprend n parties linéaires présentant respectivement n pentes de fréquence (ai), n étant supérieur ou égal à 2. 10 5. Système lidar cohérent (13) comprenant : -une source cohérente (L), -un dispositif d'émission (DE) d'un signal optique issu de la source cohérente et un dispositif de réception (DR) d'un signal rétrodiffusé par une cible (T) 15 illuminée par le lidar, -un premier (D1) et un deuxième (D2) ensemble de détection configurés pour générer respectivement un premier (Sb1) et un deuxième (Sb2) signal de battement, chaque signal de battement étant issu d'une interférence entre un signal d'oscillateur local (SOL) issu de la source cohérente et le signal 20 rétrodiffusé (Sig) par la cible (T), -une unité de traitement (UT) configurée pour numériser le premier et le deuxième signal de battement, et configurée pour déterminer pour une pluralité de n intervalles temporels (ti), n valeurs de densité spectrale (DSPi) correspondant à une transformée dans le domaine fréquentiel de la 25 corrélation croisée entre le premier et le deuxième signal de battement, -l'unité de traitement étant en outre configurée pour : *déterminer une valeur moyenne de la densité spectrale à partir des n valeurs de densité spectrale calculées, et *déterminer une information de localisation de la cible à partir de la 30 valeur moyenne de la densité spectrale. 6. Système lidar cohérent (13) selon la revendication 5 comprenant en outre un dispositif de modulation (WFC) synchronisé avec l'unité de traitement (UT) et configuré pour moduler périodiquement en fréquence la source 35 cohérente (L) de sorte que le signal d'oscillateur local (SOL) présente une 3061309 26 fréquence d'oscillateur local (foL(t)) constituée de la somme d'une valeur moyenne (f0) et d'une fréquence de modulation (f'd(t)) issue de la modulation de la source, la fréquence de modulation étant périodique selon une période de modulation (TF0), chaque période comprenant n parties 5 linéaires présentant respectivement n pentes de fréquence (ai), n étant supérieur ou égal à 2, et dans lequel l'unité de traitement (UT) est en outre configurée pour que les intervalles temporels ti soient inférieurs ou égaux à la période de modulation et pour déterminer une information de distance (d) de la cible à partir de la 10 valeur moyenne de la densité spectrale. 7. Système lidar selon l'une des revendications 5 ou 6 dans lequel le premier et/ou le deuxième ensemble de détection (D1, D2) sont des détecteurs balancés comprenant chacun un premier (PD1, PD1') et un second (PD2, 15 PD2') détecteur, les premiers détecteurs (PD1, PD11) recevant une différence entre le signal d'oscillateur local (SOL) et le signal rétrodiffusé (Sig), les deuxièmes détecteurs (PD2,PD12) recevant une somme du signal d'oscillateur local (SOL) et du signal rétrodiffusé (Sig), 20 les premier et deuxième signaux de battement (Sbl,Sb2) étant issus de la différence des intensités reçues par respectivement le premier (PD1) et le deuxième (PD2) détecteur du premier ensemble de détection (D1) et le premier (PD1') et le deuxième (PD2') détecteur du deuxième (D2) ensemble de détection. 25 8. Système lidar selon l'une des revendications 5 à 7 dans lequel le premier et le deuxième ensemble de détection (D1, D2) sont placés de sorte que la longueur des trajets empruntés par chacun des signaux jusqu'auxdits ensembles de détection soit sensiblement égale. 30 9. Un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes de la méthode de traitement selon l'une quelconque des revendications 1 à 4 lorsque ledit programme est exécuté sur un ordinateur.