FR 3055424 A1 20180302 FR 1601251 A 20160823 1 Méthode de traitement d'un signal issu d'un lidar cohérent et système lidar associé DOMAINE DE L'INVENTION L'invention concerne le domaine des lidars cohérents à modulation de 5 fréquence utilisés par exemple pour la détection de cible à longue distance. ETAT DE LA TECHNIQUE 10 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 15 de l'onde lumineuse rétrodiffusée par une cible T. Le décalage de fréquence Doppler vpop de l'onde rétrodiffusée est fonction de la vitesse radiale v de la cible T. A la réception on effectue un mélange entre l'onde lumineuse rétrodiffusée reçue S de fréquence fs et une partie de l'onde émise dénommée onde OL 20 pour « oscillateur local ». 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 l'information de 25 vitesse v de la cible T. Préférentiellement l'unité de traitement filtre électroniquement le signal de battement Sb dans une bande étroite centrée sur la fréquence nulle, en l'absence de décalage de fréquence (voir plus loin). Pour les lidars cohérents, les dispositifs d'émission et de réception utilisent 30 préférentiellement une même optique (lidar monostatique). Cette caractéristique permet d'obtenir une bonne stabilité mécanique et de réduire 3055424 2 l'influence des turbulences atmosphériques à longue distance, les trajets de propagation des ondes incidente et rétrodiffusée étant confondus. Une solution pour la télémétrie/vélocimétrie par lidar consiste à réaliser un 5 système à modulation de fréquence. Cette technique, classique en radar, trouve un intérêt particulier actuellement compte tenu des progrès des sources lasers fibrées. Grâce à la modulation de fréquence, une analyse temps/fréquence permet de remonter à la distance d de la cible et à sa vitesse v. Ce type de lidar permet également de réaliser une fonction 10 d'anémométrie laser. Un exemple d'architecture optique d'un lidar à modulation de fréquence 20 est décrit figure 2. La source cohérente est modulée en fréquence de manière à ce que la fréquence de l'oscillateur local soit modulée selon une fonction prédéterminée dénommée forme d'onde dont le contrôle est assuré 15 par le module WFC, synchronisé avec l'unité de traitement UT. Le signal optique à l'émission est amplifié par un amplificateur EDFA, l'émission et la réception utilisent une même optique O 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 20 préférentiellement positionné avant l'amplificateur EDFA mais peut également être positionné sur le trajet de l'oscillateur local. Dans ce cas le filtrage électronique dans l'unité de traitement est effectué autour de la fréquence de décalage. Une ligne à retard LR permet d'égaliser les chemins optiques de l'oscillateur local et du signal d'émission de manière à filtrer, 25 dans le domaine RF, les défauts des composants optiques placés après l'amplificateur EDFA (défaut de cross-talk du circulateur C, imperfections des traitements anti-reflet de l'optique d'émission/réception 0,...). Un exemple d'un lidar cohérent à modulation de fréquence est décrit dans le 30 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 and Imaging; and Biological and Medical Applications of Photonics Sensing and Imaging, 819202 (August 19, 2011). La figure 3 décrit le principe de 35 fonctionnement de ce lidar. 3055424 3 On se place dans la description ci-dessous dans le cas où la fréquence optique d'émission et celle de l'oscillateur local ne sont pas décalées à l'aide d'un modulateur acousto-optique. La fréquence de l'oscillateur local foL est modulée linéairement selon deux pentes de fréquence ao et al de manière périodique de période TFO. Cette fréquence optique foL peut s'écrire comme la somme d'une fréquence optique f0 constante (ici la fréquence initiale du laser) et d'une fréquence de modulation dans le domaine radiofréquence dépendante du temps fmod(t) issue de la modulation de la source laser : foL(t) = f0+ fmod(t) La figure 3 illustre la variation au cours du temps des fréquences foL(t) et fs(t), la fréquence optique f0 ayant été soustraite pour plus de clarté. Comme illustré sur la figure 3a, le signal rétrodiffusé de fréquence fs(t) est décalé temporellement d'un temps 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). Le signal de battement détecté Sb présente une composante de fréquence fs-foL. La figure 3b illustre l'évolution au cours du temps de fs -foL. 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 de la cible D et à sa vitesse v radiale par les équations : ao = 2v 2a°D et v 2v 2aiD En mesurant ces deux fréquences caractéristiques vao et val du signal de battement Sb, par exemple en réalisant une transformée de Fourier de celui-ci, on remonte à d et y. Cependant, lorsque la distance à la cible conduit à un temps de vol supérieur à la durée de la forme d'onde TFO normalisée par le nombre de pentes de fréquence (2 dans l'exemple), l'analyse directe par transformée de Fourier est insuffisante. En effet, le mélange de l'oscillateur local et du signal rétrodiffusé conduit à la disparition des plateaux et à une fréquence instantanée constamment variable, qui après analyse par transformée de Fourier, ne présentera aucun pic. 3055424 4 Un exemple de cet effet est illustré figure 4, pour une modulation de l'oscillateur local selon deux pentes de fréquence ac, = 2 MHz/ps et al = - 2 MHz/ps, et une cible se déplaçant à la vitesse de 30 m/s. La figure 4a illustre la variation temporelle de fs par rapport à fol_ et la 5 composante de fréquence de Sb fs-foi pour une distance d de 1800m, la figure 4b pour une distance d de 14000 m et la figure 4c pour une distance d de 20000 m. Dans ce cas la portée du lidar est donc limitée par le traitement du signal quelque soit la puissance du laser. Il est théoriquement possible d'allonger 10 la période de modulation TFO de la forme de d'onde, mais la plage de modulation de certains lasers étant limitée, cet allongement ne permet pas d'atteindre simultanément une résolution élevée à longue distance. En effet, compte tenu de la bande de modulation limitée du laser, il est possible d'accroître la période TFO en réduisant les pentes de fréquences pour couvrir 15 la même bande de modulation. Dans ce cas, les plateaux de fréquence existeront à plus longue distance mais, pour une durée de transformée de Fourier TFFT constante et inférieure à la fréquence de modulation TFO, la bande de modulation couverte pendant TFFT sera moindre et donc, la résolution longitudinale, proportionnelle à cette bande, sera dégradée. 20 Un but de la présente invention est de remédier aux inconvénients précités en proposant une méthode de traitement du signal de battement permettant de s'affranchir de cette limitation en retrouvant un signal présentant les plateaux de fréquences caractéristiques. 25 DESCRIPTION DE L'INVENTION La présente invention a pour objet une méthode de traitement d'un signal issu 30 d'un lidar cohérent comprenant une source cohérente modulée périodiquement en fréquence, -un signal de battement étant généré par un photodétecteur à partir de l'interférence entre un signal optique dénommé oscillateur local présentant une fréquence d'oscillateur local et un signal optique rétrodiffusé par une cible 35 illuminée par le lidar, ledit signal de battement étant numérisé, 3055424 5 -la fréquence d'oscillateur local étant 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 5 respectivement n pentes de fréquence, n étant supérieur ou égal à 2, la méthode comprenant les étapes consistant à : -moduler de manière complexe le signal de battement par la fréquence de modulation pour obtenir un signal modulé, -démoduler de manière complexe le signal modulé par n fréquences de 10 démodulation présentant chacune une pente unique égale à la pente de fréquence respective de la fréquence de modulation, pour obtenir n signaux démodulés, -déterminer n densités spectrales des n signaux démodulés, -déterminer n fréquences caractéristiques correspondant respectivement au 15 maximum des n densités spectrales, -déterminer une information de vitesse et une information de distance de la cible à partir desdites n fréquences caractéristiques. Selon un mode de réalisation, l'étape de détermination de chaque densité 20 spectrale comprend les sous étapes consistant à : -déterminer une pluralité de densités spectrales élémentaires pour une pluralité d'intervalles de temps inférieurs ou égaux à la période de modulation, -déterminer ladite densité spectrale à partir de la somme de la pluralité de 25 densités spectrales élémentaires. Préférentiellement, chaque densité spectrale élémentaire est déterminée par transformée de Fourier rapide ou FFT, et dans laquelle la densité spectrale est égale à une moyenne des densités spectrales élémentaires. Avantageusement chaque fréquence de démodulation est périodique selon 30 la période de modulation, Avantageusement les pentes de fréquence sont indicées par un indice i variant de 0 à n-1 et dans laquelle chaque fréquence de démodulation présentant une pente d'indice i est décalée temporellement par rapport à la fréquence de modulation d'un temps de décalage dépendant de i, de n et de 35 la période de modulation. 3055424 6 Selon une variante, la forme d'onde comprend 4 pentes a0, al, a2, a3 avec : al = - a0 et a3 = - a2 L'invention concerne également un système lidar cohérent comprenant : 5 - une source cohérente modulée périodiquement en fréquence, -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 lidar, - un photodétecteur configuré pour générer un signal de battement à partir de 10 l'interférence entre un signal optique dénommé oscillateur local présentant une fréquence d'oscillateur local et le signal optique rétrodiffusé, la fréquence d'oscillateur local étant 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, 15 chaque période comprenant n parties linéaires présentant respectivement n pentes de fréquence, n étant supérieur ou égal à 2, -une unité de traitement configurée pour : *numériser le signal de battement, * moduler de manière complexe le signal de battement par la fréquence 20 de modulation pour obtenir un signal modulé, *démoduler de manière complexe le signal modulé par n fréquences de démodulations présentant chacune une pente unique égale à la pente de fréquence respective de la fréquence de modulation pour obtenir n signaux démodulés, 25 *déterminer n densités spectrales des n signaux démodulés, *déterminer n fréquences caractéristiques correspondant respectivement au maximum des n densités spectrales, *déterminer une information de vitesse et une information de distance de la cible à partir desdites n fréquences caractéristiques. 30 Préférentiellement l'unité de traitement est en outre configurée pour déterminer, pour chaque densité spectrale, une pluralité de densités spectrales élémentaires pour une pluralité d'intervalles de temps inférieurs ou égaux à la période de modulation, ladite densité spectrale étant déterminée à partir de la somme de la pluralité de densités spectrales 35 élémentaires. 3055424 7 Avantageusement chaque densité spectrale élémentaire est déterminée par transformée de Fourier rapide, et dans laquelle la densité spectrale est égale à une moyenne des densités spectrales élémentaires. Avantageusement l'unité de traitement comprend n voies, une voie par 5 pente, chaque voie fonctionnant en parallèle des autres et étant configurée pour déterminer la fréquence associée. 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 10 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. La figure 2 déjà citée décrit un exemple d'architecture optique d'un lidar à modulation de fréquence. 15 La figure 3a déjà citée illustre la variation au cours du temps des fréquences foL(t) et fs(t). La figure 3b déjà citée illustre l'évolution au cours du temps de fs -foL. La figure 4a déjà citée illustre la variation temporelle de fs par rapport à foL et la composante de fréquence de Sb fs-foL pour une distance d de 1800m, la 20 figure 4b pour une distance d de 14000 m et la figure 4c pour une distance d de 20000 m. La figure 5 illustre la méthode de traitement du signal issu d'un lidar cohérent selon l'invention. La figure 6 schématise la forme d'onde périodique d'une fréquence de 25 modulation fmod(t) en fonction du temps, constituée d'une suite de 4 pentes a0, a1, a2 et a3. La figure 7 illustre le cas d'une cible se déplaçant avec une vitesse de 40 m/s à une distance de 12km, la forme d'onde du signal fmod comprenant 2 pentes 30 a0 = 0.2 MHz/ps et al = -0.2 MHz/ps pour une fréquence moyenne du laser de 1.55 prn. La figure 7a illustre la variation en fonction du temps de la fréquence de l'oscillateur local foL(t) et de la fréquence signal fs(t). La figure 7b illustre les deux composantes fréquentielles du signal de battement lidar Sb, fs-foL et 35 foL-fs. La figure 7c illustre la variation en fonction du temps de la fréquence 3055424 8 du signal modulé obtenu Smod. Les figures 7d et 7e illustrent respectivement la variation en fonction du temps des fréquences de démodulation fmod(0) (pente a0) et fmod(1) (pente al ). Les figures 7f et 7g illustrent respectivement la variation en fonction du temps du signal démodulé Sdemod(0) et Sdemod(1) . 5 La figure 8 illustre les densités spectrales SP(0) (figure 8a) et SP(1) (figure 8b) des spectres déterminés à partir des signaux Sdemod(0) et Sdemod(1). La figure 9 est équivalente à la figure 7 mais pour une cible située à une 10 distance de 18 km. La figure 9a illustre la variation en fonction du temps de la fréquence de l'oscillateur local foL(t) et de la fréquence signal fs(t). La figure 9b illustre les deux composantes fréquentielles du signal de battement lidar Sb, fs-foi et foL-fs. La figure 9c illustre la variation en fonction du temps de la fréquence du signal modulé obtenu Smod. Les figures 9d et 9e illustrent 15 respectivement la variation en fonction du temps des fréquences de démodulation fmod(0) (pente a0) et fmod(1) (pente al). Les figures 9f et 9g illustrent respectivement la variation en fonction du temps du signal démodulé Sdemod(0) et Sdemod(1). 20 La figure 10a illustre la variation au cours du temps des fréquences foL(t) et fs(t) et la figure 10b illustre l'évolution au cours du temps de fs -foi pour le cas où fmod présente 4 pentes. La figure 11a illustre la variation au cours du temps des fréquences foL(t) et 25 fs(t), la figure 11 b illustre les deux composantes fréquentielles du signal de battement lidar Sb, fs-foi et foL-fs, la figure 11c illustre la variation en fonction du temps de la fréquence du signal modulé obtenu Smod pour une forme d'onde comprenant 4 pentes. 30 Les figures 12a, 12b, 12c et 12d illustrent respectivement la variation en fonction du temps des fréquences de démodulation fmod(0) (pente a0) et fmod(1) (pente al), fmod(2) (pente a2) et fmod(3) (pente a3), et les figures 12e, 12f, 12g et 12h illustrent respectivement la variation en fonction du temps du signal démodulé Sdemod(0), Sdemod(1), Sdemod(2) et Sdemod(3)- 35 3055424 9 La figure 13 schématise un système lidar selon l'invention. La figure 14 illustre un exemple d'implémentation d'une architecture parallèle de l'unité de traitement du lidar selon l'invention. 5 DESCRIPTION DETAILLEE DE L'INVENTION La méthode de traitement 50 du signal issu d'un lidar cohérent selon l'invention est décrite figure 5. Le lidar cohérent comprend une source 10 cohérente L modulée périodiquement selon un signal RF. La modulation RF peut être réalisée directement par le courant d'injection du laser ou par un composant extérieur. Par RF on entend une onde présentant une fréquence comprise entre 1 Hz et 10 GHz, et préférentiellement entre 0.1 kHz et 10 MHz. 15 Un signal de battement Sb est généré par un photodétecteur D à partir de l'interférence entre un signal optique, dénommé oscillateur local OL présentant une fréquence oscillateur local foL(t), et un signal optique fs(t) rétrodiffusé par une cible T illuminée par le lidar. Le signal de battement est numérisé pour être traité. 20 La fréquence oscillateur local foL(t) est constituée de la somme d'une valeur moyenne f0 et d'une fréquence de modulation fmod(t) issue de la modulation de la source. foL(t) = f0 + fmod(t) Lorsque qu'aucun modulateur acoustique de décalage n'est utilisé, la 25 fréquence f0 est égale à la fréquence optique initiale de la source L. Lorsque le signal OL est décalé en fréquence par un modulateur acousto optique, la fréquence f0 est égale à la fréquence optique de la source décalée. La fréquence de modulation fmod(t) est périodique selon une période de modulation TFO, et provient de la modulation RF périodique de la source, 30 mais n'est pas équivalente à cette dernière, du fait du comportement non linéaire du laser. Typiquement la période TFO est comprise entre 1 ns et une seconde, et préférentiellement entre 100 ns et 10 ms. Pour que la méthode selon l'invention fonctionne correctement, la fréquence de modulation fmod(t) doit être telle que chaque période comprenne n parties 35 linéaires, soit n pentes de fréquence ai, avec i indice variant de 0 à n-1, qui 3055424 10 se joignent selon des sommets. Le nombre de pentes n est supérieur ou égal à 2. Avantageusement, n est pair, car, comme précisé ci-dessous, cela permet d'alterner les signes des pentes ai et de simplifier ainsi le traitement de 5 signal. En pratique compte tenu de la bande de fréquence de modulation accessible pour les lasers actuels, il est difficile d'obtenir des angles aigus pour ces sommets, et ceux-ci sont généralement arrondis, comme illustré figure 6 pour une fréquence fmod(t) constituée d'une suite de 4 pentes a0, a1, a2 et a3. La 10 forme du signal de modulation de fmod sur une période TFO est dénommée forme d'onde. Préférentiellement, la pente d'indice i+1 ai+1 présente un signe opposé à la pente d'indice i ai . Cela permet de réduire la bande de fréquence couverte en maintenant la même fraction de la période TFO pour chaque pente (et 15 donc le même ordre de grandeur d'intensité des raies pour chaque pente de fréquence). Préférentiellement, les pentes d'indices impairs sont égales à l'opposé des pentes d'indice pair. Pour un signal fmod à deux pentes al = -a0 20 Pour un signal fmod à quatre pentes al = -a0 et a3 = -a2 Dans ce dernier cas, la forme d'onde peut être divisée en quatre parties égales (conduisant à quatre raies d'intensité proche) sans avoir recours à des discontinuités de fréquence. Préférentiellement, les pentes ai sont comprises entre 0.1 MHz/ps et 25 quelque centaines de MHz/ps . A noter qu'il n'est pas aisé d'obtenir une fréquence optique oscillateur local modulée selon une suite de pentes linéaires prédéterminées tel qu'illustré sur la figure 6. Pour cela il convient de pré-corriger le signal RF de modulation de 30 la source, comme par exemple décrit dans la demande de brevet FR N° 1500603. Pour l'application de l'invention, la forme d'onde de fmod est supposé connue avec une bonne précision. Avant de décrire les étapes de la méthode 50 selon l'invention, nous allons 35 préciser la terminologie employée. 3055424 11 On dénomme modulation l'opération consistant à ajouter une fréquence à un signal initial et démodulation l'opération consistant à soustraire une fréquence du signal initial. Ainsi moduler à +f est équivalent à démoduler à -f et réciproquement. 5 Dans l'espace des temps, la modulation ou la démodulation consistent à multiplier un signal temporel initial S0(t) par un nombre, qui est un nombre réel pour une modulation/démodulation réelle (un cosinus) et un nombre complexe pour une modulation/démodulation complexe. Par exemple moduler de manière complexe par une fréquence f équivaut à 10 multiplier S0(t) par exp(2jnft). De même démoduler de manière complexe par une fréquence f équivaut à multiplier S0(t) par exp(-2jnft). Lorsque la fréquence f(t) est fonction du temps, il convient de multiplier par - t t exp 2p-cf.!' (u)du pour une modulation et par exp - 2.»-t-f f (u)du pour une 0 o - 0 démodulation. 15 La méthode 50 selon l'invention consiste en un traitement numérique spécifique d'un signal issu d'un lidar cohérent, pour la détermination des informations en vitesse et en distance d'une cible illuminée par le lidar. Plus particulièrement la méthode s'applique au traitement du signal de battement 20 lidar Sb. Les premières étapes de la méthode sont illustrées figure 7 pour le cas d'un signal fmod présentant deux pentes a0 et al . La méthode 50 selon l'invention comprend une première étape 501 consistant à moduler de manière complexe le signal de battement Sb par la 25 fréquence de modulation fmod pour obtenir un signal modulé Smod. La figure 7a illustre la variation en fonction du temps de la fréquence de l'oscillateur local foL(t) et de la fréquence signal fs(t), la fréquence optique moyenne f0 ayant été soustraite pour plus de clarté. La figure 7b illustre les deux composantes fréquentielles du signal de 30 battement lidar Sb, fs-foL et foL-fs. En effet celui-ci étant réel, il présente une composante de fréquence positive et une composante de fréquence négative. La figure 7c illustre la variation en fonction du temps de la fréquence du signal modulé obtenu Smod. Le signal réel de battement lidar Sb est modulé 3055424 12 en complexe numériquement par une fréquence associée à la forme d'onde, c'est-à-dire foi-f0 , avec f0 fréquence moyenne de la source laser L. On reconstruit alors pour le signal modulé Smod une fréquence instantanée correspondant à : 5 fs-foL (fi_040) = fs-f0 fco- fs+ (foi-f0) = 2foL - f0 - fs Ensuite une étape 502 consiste à démoduler de manière complexe le signal modulé Smod par n fréquences de démodulation f -demod(i) présentant une pente 10 unique respectivement égale à une pente de fréquence ai de la fréquence de modulation fmod, pour obtenir n signaux démodulés Sdemod(i). Ainsi n démodulations complexes sont appliquées à l'aide de n signaux numériques fdemod(i) de pente unique ai. Pour prendre en compte la périodicité de la forme d'onde, il convient de 15 revenir régulièrement à zéro. Les fréquences de démodulation fdemod(i) sont préférentiellement périodiques selon un multiple de TFO, et préférentiellement présentent une période égale à TFO. Cette égalité permet de faire coïncider les plateaux de fréquence (et donc les raies, après analyse spectrale) des différentes périodes de forme d'onde analysées : pour chaque pente de 20 fréquence ai, la raie associée apparaitra à la même fréquence va; et, par conséquent, l'énergie associé à un signal de cible sera concentrée dans la même raie après analyse temps-fréquence. Les figures 7d et 7e illustrent respectivement la variation en fonction du temps des fréquences de démodulation fmod(0) (pente a0) et fmod(1) (pente 25 al ). Afin de recaler les différentes fréquences, la fréquence de démodulation d'indice i (correspondant à une pente ai) est décalée d'un temps de décalage tdi dépendant de i, de n et de la période de modulation TFO. Préférentiellement, le temps de décalage est égal à : 30 tdi = Un * TFO Ainsi pour 2 pentes fmod(0) n'est pas décalée et fmod(1) est décalée de TFO/2 (voir figure 7e). Les figures 7f et 7g illustrent respectivement la variation en fonction du temps du signal démodulé Sdemod(0) et -demod(1) 3055424 13 Chaque démodulation correspond à la recherche du signal d'intérêt dans l'ensemble des cases distances. On retrouve alors dans le signal démodulé d'indice i un plateau de fréquence caractéristique vo; . Pour le cas de 2 pentes, Sdemod(0) permet de déterminer voo tandis que Sdemod(1) permet de 5 déterminer vol. La fréquence vo; correspond à l'écart, mesuré à un temps pour lequel fs(t)40 présente une pente de fréquence ai, entre la fréquence de démodulation fmod(i) et la fréquence fs(t)-f0, elle-même ayant été reconstituée à l'aide de la modulation du signal de battement par la fréquence foi-f0. 10 Chaque fréquence vol correspond à l'écart, mesuré à un temps pour lequel fs(t)40 présente une pente de fréquence ai, entre la fréquence de démodulation frood(i) et la fréquence fs(t)-f0, elle-même ayant été reconstituée à l'aide de la modulation du signal de battement par la fréquence foi-f0. Dans la figure 7f, le plateau le plus large correspond à +voo tandis que le 15 plateau plus faible correspond à -voo. En effet on adapte la fonction de démodulation à la fréquence d'intérêt, ici +voo (vitesse radiale de la cible positive). De la même manière dans la figure 7g, le plateau le plus large correspond à +vol tandis que le plateau plus faible correspond à -val. 20 Afin de mesurer ces fréquences caractéristiques, la méthode 50 selon l'invention comprend également une étape 503 de détermination de n densités spectrales SP(i) des n signaux démodulés Sdemod(i)- Il s'agit de réaliser une analyse temps/fréquence, soit transformée fréquentielle du signal Sdemod(i) (t), pour faire apparaître la fréquence caractéristique vo; sous 25 forme de pics. Avantageusement, il est possible d'inclure un fenêtrage temporel qui dépend de la plage de distance d'analyse et de la pente de fréquence d'analyse. La figure 8 illustre les densités spectrales SP(0) (figure 8a) et SP(1) (figure 8b) des spectres déterminés à partir des signaux Sdemod(0) et Sdemod(1). La 30 fréquence caractéristique recherchée présente la densité spectrale la plus élevée. On retrouve un pic plus faible à la fréquence opposée. Un pic négatif à la fréquence zéro est dû à un filtrage du signal pour les basses fréquences. 3055424 14 Puis la méthode selon l'invention comprend une étape 504 de détermination des n fréquences caractéristiques va; correspondant respectivement au maximum des n densités spectrales SP(i). En effet la fréquence présentant le plateau le plus large dans le signal Sdemodii(t), qui correspond à la fréquence 5 caractéristique recherchée, est celle présentant la densité spectrale la plus élevée. Un deuxième plateau de durée moins importante (et donc conduisant à une raie moins intense après analyse spectrale) est également présent mais la fréquence correspondante présente une densité spectrale plus faible que 10 celle de la fréquence caractéristique. Ce signal provient des modulations et démodulations décrites précédemment sur l'autre composante du signal de battement générée par la détection réelle (composante de fréquence négative si le signal de cible correspond à une fréquence positive ou, inversement, composante de fréquence positive si le signal de cible 15 correspond à une fréquence négative) Enfin, la méthode 50 comprend une étape 505 de détermination d'une information de vitesse v et d'une information de distance D de la cible T à partir desdites n fréquences caractéristiques va;, à partir de la formule : 20 2 v 2aiD = ' Pour 2 pentes de fréquences : 25 v = ao 2 v 2a0D et V a1 - _ 2v 2a1D A noter que les formules ci-dessus sont valables lorsque la fréquence du laser n'est pas décalée par un modulateur acousto-optique. Lorsque tel est le cas, avec fMAo fréquence de décalage, les fréquences caractéristiques se 30 calculent avec la formule : 2v 2cx,D C + J mA0 3055424 15 L'invention est bien entendu compatible d'un tel décalage en adaptant l'étape 505 de détermination de d et v à partir des valeurs des fréquences caractéristiques 5 La figure 7 correspond à une cible se déplaçant avec une vitesse de 40 m/s à une distance de 12km, et la forme d'onde du signal fmad comprend 2 pentes a0 = 0.2 MHz/ps et al = -0.2 MHz/ps pour une fréquence moyenne du laser de 1.55 pm, soit 193.41 THz. La période TFO est égale à 532 tas. Les fréquences caractéristiques détectées sont vao = 35.6 MHz et val = 67.6 10 MHz. Pour Smod(0), un pic plus faible reste à -35.6 MHz et pour Smad(1) à 67.6 MHz correspondant au plateau le plus petit. La figure 9 est équivalente à la figure 7 mais pour une cible située à une distance de 18 km. 15 La figure9a illustre la variation en fonction du temps de la fréquence de l'oscillateur local foL(t) et de la fréquence signal fs(t), La figure 9b illustre les deux composantes fréquentielles du signal de battement lidar Sb, fs-foL et foL-fs. La figure 9c illustre la variation en fonction du temps de la fréquence du 20 signal modulé obtenu Smod. Les figures 9d et 9e illustrent respectivement la variation en fonction du temps des fréquences de démodulation fmod(0) (pente a0) et fmod(1) (pente al). Les figures 9f et 9g illustrent respectivement la variation en fonction du temps 25 du signal démodulé Sdemad(0) et -demod(1) . On constate que les plateaux réapparaissent, même à plus longue distance. Les fréquences caractéristiques détectées sont vao = 27.6 MHz et val = 75.6 MHz. Il ne reste presque plus de pics à - 27.6 MHz et -75.6 MHz 30 Ainsi la méthode proposée évite de tester toutes les cases distance (solution couteuse en calcul) et permet, par une simple opération modulation/démodulation, de retrouver la distance de la cible, tant que la puissance du laser reste suffisante. Les pics issus du signal rétrodiffusé réapparaissent, permettant ainsi d'obtenir une méthode qui n'est plus limitée 35 par le traitement du signal, mais uniquement par la puissance du laser. 3055424 16 Le calcul s'effectue à partir du signal de battement Sb(t) numérisé au fur et à mesure de l'écoulement du temps. Mathématiquement, l'étape 501 de modulation par la fréquence 5 fmod(t) = foL(t) -f0 revient à multiplier le signal Sb(t) par un nombre complexe C(t), également numérisé, égal à : C = exp[2»rf (foL (u) - fo )du = exp 2»-cf fmod(u)dui o 10 Soit Smod(t) = C* Sb(t) f0 : fréquence du laser sans modulation foi : fréquence de l'oscillateur local Puis dans l'étape 502 de démodulation, chaque démodulation revient à 15 multiplier le signal Smod(t) par un nombre complexe Ci(t) défini comme suit : C; = exp -2/.4 u- -ng,(u)+ floorr i+1\\ *Ti F° \ 2 ' 2 r Avec i indice de la pente ai, avec i varaint de 0 à n-1, TF0 : période de la forme d'onde r 4u i gi(u)-= floor 20 floor étant la fonction arrondi inférieure (par exemple floor(2.6) =2 et floor(- 3.2)=-4) Soit au final : nTF0 ni 25 S de modh(t) S b(t). exp 2 jn- foL(u)- ( ( -ngi(u)+ floor( T o u * FO 2 2 ) 2 S de lm& (t) = S b(t). exp2 jn-f { o fmod (u)- ai( u -( -n g ,(u) + floor( i ± if * TF0iclul 2 2 ) 2 La partie aiu correspond à la partie linéaire, la partie n/2.gi(u).TF0/2 exprime le retour à zéro régulier et le décalage temporel, et la partie floor(i+1/2)* 3055424 17 TFO/2 correspond à un décalage en fréquence permettant de ramener à fréquence nulle une situation pour laquelle la vitesse et la distance de la cible sont nulles. Ce dernier décalage en fréquence compense un effet parasite généré par le décalage temporel associé à la fonction gi(u). 5 A noter que si les sommets de la forme d'onde sont arrondis, ces équations restent valables car cet arrondi est pris en compte dans la définition de Smod(t). L'étape 503 d'obtention des densités spectrales SP(i) s'effectue typiquement 10 par transformée fréquentielle, en prenant le carré du module de la transformée de Fourier du signal temporel Sdemodti(t) : u - n - g ,(4 )+ floor +1 *FO du 2 ), 2 _ 2 } 7 {S inod, i(t)}i 2 = FFT Sb(t). exp 2J 2r f fnod (u)- ai o 15 Préférentiellement, l'étape 503 de détermination de chaque densité spectrale comprend les sous étapes consistant à : -déterminer une pluralité de densités spectrales élémentaires pour une pluralité d'intervalles de temps ôt inférieurs ou égaux à la période de 20 modulation TFO, -déterminer chaque densité spectrale d'indice i SP(i) à partir de la somme de la pluralité de densités spectrales élémentaires. Préférentiellement, chaque densité spectrale élémentaire est déterminée par transformée de Fourier rapide FFT. 25 En effet pour simplifier le traitement les transformées de Fourier réalisées pendant la période de la forme d'onde peuvent être directement sommées (en puissance). On réalise donc une accumulation incohérente de densités spectrales élémentaires, qui sont ensuite moyennées. 30 Cette opération permet de réaliser des calculs rapides, chaque densité spectrale élémentaire étant calculée sur un temps ôt court. Par exemple, pour une fréquence d'échantillonnage de 125 MHz et une période TFO de 500 ps, réaliser plusieurs calculs de FFT sur un ôt de 30ps 3055424 18 (correspondant à 4000 points) est beaucoup plus performant que d'effectuer un calcul sur la durée totale de TFO (trop de points). De plus, réaliser une moyenne sur un certain nombre de FFT pendant une période TFO permet d'améliorer le rapport signal sur bruit de SPi(v) sans 5 perte d'information en choisissant judicieusement les instants où l'on accumule le signal. En effet, le bruit est généralement limité par le bruit de photon. Le signal et le bruit présentent une statistique de chi2 et, par conséquent, le rapport signal sur bruit diminue en 1/sqrt(N) où N désigne le nombre de densités spectrales moyennées. Les figures 8a et 8b 10 correspondent à une moyenne des densités spectrales SP(0) et SP(1) effectuée sur plusieurs centaines de FFT (N= 864). Le signal décrit par les fréquences instantanées entre les plateaux a une puissance proportionnelle à la puissance du signal concentré dans les plateaux de fréquence, mais il est réparti sur un nombre de canaux spectraux 15 nettement plus importants. Après analyse temps/fréquence, ce signal est donc dilué dans la bande d'analyse et conduit : -à courte distance à un bruit additionnel qui réduit le rapport signal sur bruit (SNR). Cette diminution n'est cependant pas importante puisqu'à courte distance, le SNR est élevé. Selon un mode de réalisation, si l'on souhaite 20 éviter cette diminution, une étape de recherche des plages temporelles conduisant à un plateau de fréquence est ajoutée pour garantir un SNR maximal. -à longue distance (pour un SNR plus faible), le bruit additionnel l'effet reste inférieur aux bruits de détection (notamment le bruit de photon de l'oscillateur 25 local), mais réduire le temps d'accumulation aux seuls instants où le signal est présent permet de réduire le bruit de détection. Par ailleurs, réaliser une moyenne sur un certain nombre de FFT pendant une période TFO permet de caler la durée d'une transformée de Fourier sur le 30 temps de cohérence de la cible (qui dépend notamment des mouvements de cette cible) ce qui optimise également le rapport signal sur bruit. La densité spectrale calculée est préférentiellement égale à la moyenne des densités spectrales élémentaires, afin de toujours obtenir des valeurs 35 numériques normalisées. 3055424 19 D'un point de vue pratique, les calculs de modulation/démodulation, puis de FFT et du carré du module s'effectuent au fur et à mesure de la numérisation du signal de battement, en temps réel. Puis au bout d'un 5 certain temps d'accumulation, les densités spectrales SP(i) sont obtenues en effectuant la moyenne des densités spectrales élémentaires accumulées (voir plus loin figure 14). L'invention s'applique pour toute valeur de n supérieure ou égale à 2. La 10 figure 7 illustre la méthode appliquée pour n=2. 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, 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. Ceci permet d'obtenir une redondance et donc d'utiliser 15 une des équations 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é 20 de manière itérative (Iteratively reweighted least squares - IRLS). De même que pour la figure 7 pour le cas d'une forme d'onde à 2 pentes, la figure 10a illustre la variation au cours du temps des fréquences foL(t) et fs(t), la fréquence optique moyenne f0 ayant été soustraite pour plus de clarté. La 25 figure 10b illustre l'évolution au cours du temps de fs -foi pour le cas à 4 pentes. On voit sur la figure 10b que cette variation de fréquence au cours du temps présente 4 plateaux correspondant aux 4 fréquences caractéristiques. De même que pour la figure 7 pour le cas d'une forme d'onde à 2 pentes, la 30 figure 11 b illustre les deux composantes fréquentielles du signal de battement lidar Sb, fs-foi et foi-fs, la figure 11c illustre la variation en fonction du temps de la fréquence du signal modulé obtenu Smod pour une forme d'onde comprenant 4 pentes. Les figures 12a, 12b, 12c et 12d illustrent respectivement la variation en fonction du temps des fréquences de 35 démodulation fmod(0) (pente a0) et fmod(1) (pente al ), fmod(2) (pente a2) et 3055424 20 fmod(3) (pente a3), et les figures 12e, 12f, 12g et 12h illustrent respectivement la variation en fonction du temps du signal démodulé Sdemod(0), Sdemod(1), Sdemod(2) et Sdemod(3). Pour 4 pentes, fmod(0) n'est pas décalée (voir figure 12a), fmod(1) est décalée 5 de TFO/4 (voir figure 12b), fmod(2) est décalée de TF0/2 (voir figure 12c), fmod(3) est décalée de 3/4.1-F0 (voir figure 12d). Les figures 10 et 11 correspondent au cas d'une cible située à 12 km, se déplaçant à 40 m/s, la fréquence fmod présentant les valeurs de pentes 10 suivantes (laser de fréquence optique f0 = 1.55 pm): a0 = 0.2 MHz/ps al = -0.2 MHz/ps a2 = 0.3 MHz/ps a3 = -0.3 MHz/ps 15 Par transformée dans l'espace fréquentiel, on détecte les fréquences caractéristiques (plateaux les plus longs) : 35.6 MHz (ao), 67.6 MHz (ai), 27.6 MHz (a2) et 75.6 MHz (a3). Il reste également des pics plus faibles aux fréquences opposées. 20 L'invention concerne également un système lidar cohérent illustré figure 13 et comprenant : - une source cohérente L modulée périodiquement en fréquence, - un dispositif d'émission DE d'un signal optique issu de la source cohérente 25 et un dispositif de réception DR d'un signal rétrodiffusé par une cible T illuminée par le lidar, - un photodétecteur D configuré pour générer un signal de battement Sb à partir de l'interférence entre un signal optique dénommé oscillateur local présentant une fréquence d'oscillateur local foL(t) et le signal optique 30 rétrodiffusé, la fréquence d'oscillateur local foL(t) étant constituée de la somme d'une valeur moyenne f0 et d'une fréquence de modulation fmod(t) issue de la modulation de la source, la fréquence de modulation étant périodique selon une période de modulation TFO, chaque période comprenant n parties linéaires présentant respectivement n pentes de fréquence ai, n 35 étant supérieur ou égal à 2, i variant de 0 à n-1, 3055424 21 -une unité de traitement UT configurée pour : *numériser le signal de battement, *moduler de manière complexe le signal de battement Sb par la fréquence de modulation fmod pour obtenir un signal modulé Smod, 5 *démoduler de manière complexe le signal modulé Smod par n fréquences de démodulations f -demod(i) présentant chacune une pente unique égale à la pente de fréquence ai respective de la fréquence de modulation pour obtenir n signaux démodulés Sdemod(i), *déterminer n densités spectrales SP(i) des n signaux démodulés, 10 *déterminer n fréquences caractéristiques va; correspondant respectivement au maximum des n densités spectrales SP(i), *déterminer une information de vitesse v et une information de distance d de la cible T à partir desdites n fréquences caractéristiques va;. 15 Avantageusement, l'unité de traitement UT est en outre configurée pour déterminer, pour chaque densité spectrale, une pluralité de densités spectrales élémentaires pour une pluralité d'intervalles de temps inférieurs ou égaux à la période de modulation TFO, la densité spectrale SP(i) étant déterminée à partir de la somme de la pluralité de densités spectrales 20 élémentaires. Préférentiellement chaque densité spectrale élémentaire est déterminée par transformée de Fourier rapide FFT. Préférentiellement la densité spectrale est égale à une moyenne des densités spectrales élémentaires. 25 Préférentiellement, l'unité de traitement UT comprend n voies, une voie par pente, chaque voie fonctionnant en parallèle des autres et étant configurée pour déterminer la fréquence associée. En effet, la modulation et la démodulation peuvent être conduites simultanément, conduisant ainsi à un coût de calcul réduit (constitué d'une seule multiplication complexe). 30 Un exemple d'implémentation d'une architecture parallèle à 4 voies (4 pentes) de l'unité de traitement UT est illustré figure 14. Le signal de battement Sb est numérisé à l'aide d'un convertisseur analogique /digital ADC (par exemple un convertisseur 14 bits, 125 MHz) puis éventuellement filtré par un filtre fréquentiel F. Le signal numérisé et 35 filtré est ensuite réparti sur les 4 voies. Chaque voie fonctionne en parallèle 3055424 22 des autres et effectue la même chaîne de traitement. Seule la valeur de la fréquence de démodulation f -clemod(i) (et son décalage temporel) est différente d'une chaîne à l'autre. Le module 2 permet de définir l'amplitude et la phase des fonctions de 5 modulation et démodulation C et fmod(i). Le produit de ces fonctions est ensuite évalué dans le module 3. Le module 4 permet de réaliser la multiplication complexe du signal de battement Sb numérisé par la fonction calculée dans le module 3 (produit de la fonction de modulation C par la fonction de démodulation frood(0). 10 Le module 5 réalise les transformées de Fourier rapide (FFT) complexes. Le module 6 calcule la norme au carré des transformées de Fourier. Le module 7 réalise la somme des densités spectrales de puissance pendant un temps déterminé par les caractéristiques fournies par le module 12 (durée, cadence...). Ce résultat est transféré dans un buffer 8 avant d'être 15 transféré via un serveur TCP 9 et exploité, dans une seconde partie du traitement de signal qui peut être réalisé plus lentement. Cette seconde partie, module 11 sur la figure 15, permet de réaliser la détection de pics, l'évaluation des fréquences et effectue le calcul de v et d en prenant en compte l'ensemble de ces fréquences caractéristiques. Cette étape peut, par 20 exemple, être réalisée par la technique des moindres carrés ou, par « Iteratively reweighted least squares » (IRLS), connue dans la littérature. L'invention concerne également un produit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes de la 25 méthode de traitement selon l'invention. Dans les différentes variantes de réalisation du système selon l'invention, les modules de calcul peuvent être agencés selon différentes architectures, en particulier chaque étape du procédé peut être implémentée par un module distinct ou au contraire l'ensemble des étapes peuvent être regroupées au 30 sein d'un module de calcul unique. Chacun des modules de calcul que comporte le système selon l'invention peut être réalisé sous forme 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 35 circuit intégré propre à une application (connu aussi sous le nom anglais 3055424 23 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 Gate Array »). 5 REVENDICATIONS1. Méthode de traitement (50) d'un signal issu d'un lidar cohérent comprenant une source cohérente (L) modulée périodiquement en fréquence, -un signal de battement (Sb) étant généré par un photodétecteur (D) à partir 5 de l'interférence entre un signal optique dénommé oscillateur local présentant une fréquence d'oscillateur local (foL(t)) et un signal optique rétrodiffusé par une cible (T) illuminée par le lidar, ledit signal de battement étant numérisé, -la fréquence d'oscillateur local (foL(t)) étant constituée de la somme d'une valeur moyenne (f0) et d'une fréquence de modulation (fmod(t)) issue de la 10 modulation de la source, la fréquence de modulation étant périodique selon une période de modulation (TFO), chaque période comprenant n parties linéaires présentant respectivement n pentes de fréquence (ai), n étant supérieur ou égal à 2, la méthode comprenant les étapes consistant à : 15 -moduler (501) de manière complexe le signal de battement (Sb) par la fréquence de modulation (fmod) pour obtenir un signal modulé (Smod), -démoduler (502) de manière complexe le signal modulé (Smod) par n fréquences de démodulation (fdemod(i)) présentant chacune une pente unique égale à la pente de fréquence (ai) respective de la fréquence de modulation 20 (fmod), pour obtenir n signaux démodulés (Sdemod(I)), -déterminer (503) n densités spectrales (SP(i)) des n signaux démodulés (Sdemod(I)), -déterminer (504) n fréquences caractéristiques (va;) correspondant respectivement au maximum des n densités spectrales (SP(i)), 25 -déterminer (505) une information de vitesse (y) et une information de distance (d) de la cible (T) à partir desdites n fréquences caractéristiques (vai). 2. Méthode selon la revendication 1 dans laquelle l'étape de détermination 30 de chaque densité spectrale comprend les sous étapes consistant à : -déterminer une pluralité de densités spectrales élémentaires pour une pluralité d'intervalles de temps inférieurs ou égaux à la période de modulation ((TFO), 3055424 25 -déterminer ladite densité spectrale (SP(i)) à partir de la somme de la pluralité de densités spectrales élémentaires. 3. Méthode selon la revendication 2 dans laquelle chaque densité spectrale 5 élémentaire est déterminée par transformée de Fourier rapide (FFT), et dans laquelle la densité spectrale est égale à une moyenne des densités spectrales élémentaires. 4. Méthode selon l'une des revendications 1 ou 2 dans laquelle chaque (0) 10 fréquence de démodulation (fdemodest périodique selon la période de modulation (TFo). 5. Méthode selon la revendication 4 dans laquelle les pentes de fréquence (ai) sont indicées par un indice i variant de 0 à n-1 et dans laquelle la chaque 15 fréquence de démodulation (fdemod( i)) présentant une pente d'indice i est décalée temporellement par rapport à la fréquence de modulation (fmod) d'un temps de décalage (tdi) dépendant de i, de n et de la période de modulation (TFO). 20 6. Méthode selon l'une des revendications précédentes dans laquelle la forme d'onde comprend 4 pentes a0, a1, a2, a3 avec : al = - a0 et a3 = - a2 7. Système lidar cohérent comprenant : 25 - une source cohérente (L) modulée périodiquement en fréquence, -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 par le lidar, -un photodétecteur (D) configuré pour générer un signal de battement (Sb) à 30 partir de l'interférence entre un signal optique dénommé oscillateur local présentant une fréquence d'oscillateur local (foL(t)) et le signal optique rétrodiffusé, la fréquence d'oscillateur local (foL(t)) étant constituée de la somme d'une valeur moyenne (f0) et d'une fréquence de modulation (fmod(t)) issue de la modulation de la source, la fréquence de modulation étant 35 périodique selon une période de modulation (TF0), chaque période 3055424 26 comprenant n parties linéaires présentant respectivement n pentes de fréquence (ai), n étant supérieur ou égal à 2, -une unité de traitement (UT) configurée pour : *numériser le signal de battement, *moduler de manière complexe le signal de battement (Sb) par la fréquence de modulation (fmod) pour obtenir un signal modulé (Sn-KA), *démoduler de manière complexe le signal modulé (Smog) par n fréquences de démodulations (fdemod(i)) présentant chacune une pente unique égale à la pente de fréquence (ai) respective de la fréquence de modulation pour obtenir n signaux démodulés (Sdemod(I)), *déterminer n densités spectrales (SP(i)) des n signaux démodulés, *déterminer n fréquences caractéristiques (va;) correspondant respectivement au maximum des n densités spectrales (SP(i)), *déterminer une information de vitesse (v) et une information de distance (d) de la cible (T) à partir desdites n fréquences caractéristiques (va;). 8. Système lidar selon la revendication 7 dans laquelle l'unité de traitement (UT) est en outre configurée pour déterminer, pour chaque densité spectrale, une pluralité de densités spectrales élémentaires pour une pluralité d'intervalles de temps inférieurs ou égaux à la période de modulation (TFo), ladite densité spectrale (SP(i)) étant déterminée à partir de la somme de la pluralité de densités spectrales élémentaires. 9. Système lidar selon la revendication 8 dans laquelle chaque densité spectrale élémentaire est déterminée par transformée de Fourier rapide (FFT), et dans laquelle la densité spectrale est égale à une moyenne des densités spectrales élémentaires. 10. Système lidar selon l'une des revendications 7 à 9 dans laquelle l'unité de traitement (UT) comprend n voies, une voie par pente, chaque voie fonctionnant en parallèle des autres et étant configurée pour déterminer la fréquence associée. 3055424 27 11. Un produit programme d'ordinateur, ledit programme d'ordinateur comprenant des instructions de code permettant d'effectuer les étapes de la méthode de traitement selon l'une quelconque des revendications 1 à 6.