Brevet FR3043788A1 — publié le 20170519

DOMAINE TECHNIQUE AUQUEL SE RAPPORTE L'INVENTION La présente invention concerne de manière générale le domaine des systèmes et méthodes de localisation de fond marin basés sur la réflexion d'ondes acoustiques ou lumineuses. Elle concerne plus particulièrement un système bathymétrique ou de mesure de profondeur du fond marin ayant une grande précision de mesure. ARRIERE-PLAN TECHNOLOGIQUE Une grande partie de l'activité des services hydrographiques de différents pays consiste à déterminer la topographie des fonds marins grâce à des mesures de profondeur et de relief de ces fonds. Ce type de mesure, appelées mesures de bathymétrie, est également effectué par de nombreuses entreprises travaillant dans les domaines de la construction marine. Il existe différents systèmes permettant d'effectuer des mesures de bathymétrie basés sur la propagation et la réflexion d'ondes acoustiques, les sonars, ou d'ondes optiques, les lidars. Un sonar ou un lidar peut être tracté par un bateau ou fixé à la coque d'un véhicule porteur, par exemple un bateau ou un sous-marin. Des systèmes sonar utilisant la réflexion des ondes acoustiques sont couramment utilisés pour obtenir des mesures de bathymétrie et déterminer la distance entre le fond marin et le véhicule porteur du sonar. On connaît en particulier des systèmes sonar dits sonars multifaisceaux. Un sonar multifaisceaux comporte au moins un émetteur acoustique et plusieurs récepteurs acoustiques, et fonctionne par récurrence acoustique, où chaque récurrence acoustique consiste en l'émission d'un signal acoustique incident vers un fond marin suivi de la réception de plusieurs signaux réfléchis sur le fond marin sur les différents détecteurs acoustiques. Un sonar multifaisceaux de type classique fournit ainsi à chaque récurrence acoustique, la mesure de la profondeur du fond marin dans plusieurs directions (correspondant aux n faisceaux acoustiques réfléchis détectés) formant une fauchée, généralement perpendiculaire à l'axe du véhicule porteur. La combinaison de plusieurs récurrences acoustiques d'un sonar multifaisceaux permet de déterminer la bathymétrie des fonds marins sur un couloir le long de la route du véhicule porteur. 3043788 2 Lorsque les conditions le permettent et en particulier lorsque le fond marin est peu profond (moins d'environ une dizaine de mètres), un système LIDAR, basé sur la réflexion d'un faisceau laser, peut aussi être utilisé pour obtenir des mesures de bathymétrie et déterminer la distance entre le fond marin 5 et le véhicule porteur de ce système LIDAR. On connait en particulier des systèmes de LIDAR à balayage de faisceau laser qui permettent, à chaque récurrence laser, de mesurer la profondeur du fond marin dans n directions (n réflexions de faisceau laser suivant n directions) formant une fauchée, en général perpendiculaire à l'axe du véhicule porteur. La combinaison de plusieurs 10 récurrences lidar permet ainsi de déterminer la bathymétrie d'un fond marin sur un couloir le long de la route du véhicule porteur. Pour effectuer une topographie de fond marin à partir des mesures provenant d'un système de type sonar multi-faisceaux ou lidar, il est nécessaire de prendre en compte la position et l'orientation du véhicule porteur du système de 15 mesure, afin de positionner correctement les points de mesures correspondant aux différents faisceaux, appelés sondes, sur une carte géographique. Ainsi, les systèmes de mesure de la topographie des fonds marins combinent en général un système SONAR ou LIDAR avec un système de navigation inertielle (INS) lui-même fréquemment couplé à un récepteur de 20 géolocalisation et de navigation par satellite (GPS). L'utilisation d'un système de navigation inertielle permet de déterminer la position (x, y, z) et l'orientation (cap, roulis, tangage) du véhicule porteur et ainsi d'en tenir compte dans le positionnement géographique des sondes. La méthode consistant à associer un système de mesure de la distance 25 au fond avec un système inertiel est la méthode actuellement utilisée par l'ensemble des opérateurs de systèmes de bathymétrie dans le monde. Cette méthode présente néanmoins deux inconvénients majeurs. D'une part le coût des systèmes de navigation inertiels requis, incluant fréquemment un récepteur GPS précis de type RTK, est relativement élevé. D'autre part, tandis 30 que les systèmes inertiels permettent de déterminer précisément l'attitude (roulis, tangage), la position (x, y) dans un plan horizontal (avec des erreurs relatives de l'ordre de quelques centimètres sur des durées de quelques dizaines de secondes) et, pour les plus élaborés d'entre eux, le cap du véhicule porteur, ces systèmes inertiels ont pour défaut d'être relativement imprécis dans la 3043788 3 détermination de l'altitude (z) du véhicule porteur suivant un axe vertical (avec des erreurs relatives atteignant rapidement plusieurs dizaines de centimètres voir quelques mètres sur des durées de quelques dizaines de secondes). La détermination de l'altitude (z) du véhicule porteur dépend donc de façon très 5 importante de la disponibilité d'une information GPS de très grande qualité, avec une précision de l'ordre de quelques centimètres. Or cette information GPS de position verticale précise n'est en général pas disponible de façon continue dans l'espace et en particulier n'est pas disponible au large des côtes. Certes, il est en général possible de déterminer l'altitude moyenne du 10 véhicule porteur sur une période de temps de quelques minutes, typiquement en moyennant les valeurs fournies par un capteur imprécis. Il est par contre particulièrement difficile de déterminer très précisément les variations d'erreur d'altitude à court terme (quelques secondes à quelques dizaines de secondes). Or les variations d'erreur d'altitude à court terme ont un fort impact sur la qualité de la 15 topographie obtenue. Les techniques actuellement disponibles pour résoudre ce problème de détermination des variations d'erreur d'altitude du véhicule porteur à court terme se basent en général sur des hypothèses concernant soit le mouvement du véhicule porteur soit la topographie du sol soit la combinaison des deux. Ainsi, de 20 nombreux systèmes inertiels, équipant les bateaux hydrographiques, utilisent des filtres de houles. Par construction de ces filtres de houle, on suppose que les fréquences des mouvements en altitude sont situées dans une gamme de relativement hautes fréquences ce qui permet d'éliminer les fréquences plus basses par filtrage numérique et donc de réduire l'erreur intrinsèque d'altitude des 25 systèmes inertiels. Certaines méthodes de traitement purement acoustiques, telle la méthode décrite dans le document de brevet US2013/0016584, supposent également que les fréquences des mouvements du véhicule porteur sont plus élevées que les fréquences de variation de l'altitude du fond marin le long du parcours du véhicule porteur. 30 Malheureusement, les hypothèses de ces méthodes ne sont en général pas parfaitement satisfaites ce qui entraine des imperfections des cartes topographiques déterminées en les appliquant. Un des buts de l'invention est de proposer un système et une méthode de bathymétrie permettant de corriger les mouvements du véhicule porteur, et en 3043788 4 particulier les mouvements verticaux, sans faire d'hypothèse sur la topographie du fond marin ou sur les mouvements du véhicule porteur. OBJET DE L'INVENTION 5 Afin de remédier à l'inconvénient précité de l'état de la technique, la présente invention propose un système bathymétrique comprenant un calculateur et un dispositif d'émission et de réception multi-faisceaux attaché à un véhicule porteur, le dispositif d'émission et de réception multi-faisceaux étant configuré pour émettre, à une première récurrence, un premier faisceau incident vers le fond 10 marin suivant une première surface d'émission et pour détecter, à la première récurrence, l'amplitude d'une première pluralité de faisceaux réfléchis suivant une première pluralité de directions par réflexion du premier faisceau incident sur le fond marin, le calculateur étant configuré pour déterminer, à partir de l'amplitude de la première pluralité de faisceaux réfléchis détectés, les coordonnées spatiales 15 tridimensionnelles d'une première pluralité de points de sonde situés à l'intersection du fond marin et de la première surface d'émission, ladite première pluralité de points de sonde formant une première fauchée à la première récurrence. Plus particulièrement, on propose selon l'invention un système 20 bathymétrique dans lequel le dispositif d'émission et de réception multi-faisceaux est configuré pour émettre, à une deuxième récurrence, un deuxième faisceau incident vers le fond marin suivant une deuxième surface d'émission, la deuxième surface d'émission étant non parallèle à la première surface d'émission, et les premier et deuxième faisceaux incidents non interférents entre eux, autrement dit 25 les premier et deuxième faisceaux incidents étant décorrélés temporellement et/ou en fréquence et/ou en phase de telle façon que leur inter-corrélation à la réception soit nulle, le dispositif d'émission et de réception multi-faisceaux étant configuré pour détecter l'amplitude, à la deuxième récurrence, d'une deuxième pluralité de faisceaux réfléchis, suivant une deuxième pluralité de directions, par réflexion du 30 deuxième faisceau incident sur le fond marin, le calculateur étant configuré pour déterminer, à partir de l'amplitude de la deuxième pluralité de faisceaux réfléchis détectés, les coordonnées spatiales tridimensionnelles d'une deuxième pluralité de points de sonde situés à l'intersection du fond marin et de la deuxième surface d'émission, ladite deuxième pluralité de points de sonde formant une deuxième 3043788 5 fauchée, le dispositif d'émission et de réception multi-faisceaux étant configuré pour que la première fauchée et la deuxième fauchée comportent au moins un point de sonde redondant ayant des coordonnées spatiales bidimensionnelles dans un plan horizontal identiques à la première récurrence et à la deuxième 5 récurrence, le calculateur étant configuré pour calculer une variation d'erreur d'altitude du véhicule porteur à partir des variations de coordonnées spatiales en altitude d'au moins un point de sonde redondant entre la première récurrence et la deuxième récurrence, et le calculateur étant configuré pour corriger les coordonnées spatiales tridimensionnelles de la première fauchée et/ou de la 10 deuxième fauchée, en fonction de la variation d'erreur d'altitude du véhicule porteur calculée entre la première récurrence et la deuxième récurrence. L'invention est ainsi mise en oeuvre en combinant, d'une part, un dispositif physique configuré pour effectuer à chaque instant la mesure des sondes de telle façon que des points de sondes effectuées à deux instants de 15 récurrence éloignés temporellement de plusieurs secondes se superposent spatialement en partie, et, d'autre part, un algorithme permettant de combiner des mesures effectuées à des instants de récurrence différents afin de déterminer la variation d'altitude du véhicule porteur entre ces instants de récurrence. L'invention exploite la superposition spatiale partielle de différentes 20 fauchées effectuées à différentes récurrences pour déterminer précisément les variations d'erreur d'altitude entre ces différentes récurrences. L'invention permet ainsi de corriger les variations d'erreur d'altitude entre deux récurrences non corrélées entre elles, sans faire d'hypothèse sur les mouvements du véhicule porteur ou sur la topographie du fond marin. 25 Le système permet de corriger les mesures de bathymétrie en altitude entre deux récurrences, que ces récurrences soient consécutives ou non, dès qu'il existe au moins un point de sonde redondant appartenant à deux fauchées associées à ces deux récurrences. Ce système permet, de manière récursive, d'améliorer progressivement la précision des mesures de bathymétrie au fur et à 30 mesure des récurrences. Le système est configuré pour sonder le fond marin suivant au moins deux surfaces incidentes non parallèles entre elles, et pour collecter les ondes réfléchies dans plusieurs directions pour chaque surface incidente sondée. Ce système de correction des variations d'erreur d'altitude du véhicule porteur peut 3043788 6 être adapté aisément à des systèmes bathymétriques existants, ayant par exemple une antenne multi-faisceaux à deux dimensions, ou deux antennes multifaisceaux disposées transversalement l'une par rapport à l'autre. D'autres caractéristiques non limitatives et avantageuses du système 5 bathymétrique conforme à l'invention, prises individuellement ou selon toutes les combinaisons techniquement possibles, sont les suivantes : - le dispositif d'émission et de réception multi-faisceaux comporte un sonar multi-faisceaux ayant une première et une deuxième antennes linéaires d'émission et de réception multi-faisceaux, la première antenne linéaire étant 10 disposée transversalement par rapport à la deuxième antenne linéaire ; - d'une part, la première antenne est configurée pour émettre le premier faisceau incident lors de la première récurrence, et la deuxième antenne est configurée pour détecter la première fauchée lors de la première récurrence ; et, d'autre part, la deuxième antenne est configurée pour émettre le deuxième 15 faisceau incident lors de la deuxième récurrence, et la première antenne est configurée pour détecter la deuxième fauchée lors de la deuxième récurrence ; - la première antenne et la deuxième antenne sont disposées respectivement à environ 45 degrés d'une ligne de foi du véhicule porteur ; - le véhicule porteur comporte une ligne de foi et la première antenne est 20 alignée sur cette ligne de foi; - le dispositif d'émission et de réception multi-faisceaux comporte un sonar multi-faisceaux ayant au moins une antenne de réception multi-capteurs à deux dimensions ; - le système bathymétrique comprend deux antennes d'émission et une 25 antenne de réception linéaire multi-capteurs ; - le dispositif d'émission et de réception multi-faisceaux comporte un système lidar configuré pour émettre un faisceau laser vers le fond marin suivant une première et respectivement une deuxième surface d'émission, et pour détecter une pluralité de faisceaux laser réfléchis sur le fond marin, formant la 30 première fauchée et respectivement la deuxième fauchée ; - le système lidar comporte une source laser à balayage bidirectionnel de faisceau ou deux sources laser à balayage unidirectionnel de faisceau. Selon un mode de réalisation particulier, le système bathymétrique comporte en outre au moins un capteur externe choisi parmi un système de 3043788 7 géolocalisation par satellite (GPS), un capteur de profondeur ou une centrale de navigation inertielle, ce capteur externe étant configuré pour fournir au moins une mesure complémentaire de position et/ou d'orientation du dispositif d'émission et de réception multi-faisceaux, et le calculateur est configuré pour calculer une 5 variation d'erreur d'altitude du véhicule porteur à partir des variations de coordonnées spatiales en altitude du au moins un point de sonde redondant entre la première récurrence et la deuxième récurrence et en fonction de ladite au moins une mesure complémentaire de position et/ou d'orientation. L'invention propose également une méthode de bathymétrie comprenant 10 les étapes suivantes : a) Emission, à une première récurrence, par un dispositif d'émission et de réception multi-faisceaux, d'un premier faisceau incident apte à se propager dans un milieu sous-marin vers un fond marin suivant une première surface d'émission ; 15 b) Détection par le dispositif d'émission et de réception multi-faisceaux, à la première récurrence, de l'amplitude d'une première pluralité de faisceaux réfléchis suivant une première pluralité de directions, par réflexion du premier faisceau incident sur le fond marin ; c) calcul, à partir de l'amplitude de la première pluralité de faisceaux 20 réfléchis détectés, des coordonnées spatiales tridimensionnelles d'une première pluralité de points de sonde situés à l'intersection du fond marin et de la première surface d'émission, ladite première pluralité de points de sonde formant une première fauchée à la première récurrence ; 25 d) émission, à une deuxième récurrence, par le dispositif d'émission et de réception multi-faisceaux, d'un deuxième faisceau incident apte à se propager dans le milieu sous-marin vers le fond marin suivant une deuxième surface d'émission, la deuxième surface d'émission étant non parallèle à la première surface d'émission, et le deuxième faisceau 30 incident et le premier faisceau incident étant non interférents entre eux, autrement dit les premier et deuxième faisceaux incidents étant décorrélés temporellement et/ou en fréquence et/ou en phase de telle façon que leur inter-corrélation à la réception soit nulle; e) détection, par le dispositif d'émission et de réception multi-faisceaux, à 3043788 8 la deuxième récurrence, de l'amplitude d'une deuxième pluralité de faisceaux réfléchis suivant une deuxième pluralité de directions, par réflexion du deuxième faisceau incident sur le fond marin ; f) calcul, à partir de l'amplitude de la deuxième pluralité de faisceaux 5 réfléchis détectés, des coordonnées spatiales tridimensionnelles d'une deuxième pluralité de points de sonde situés à l'intersection du fond marin et de la deuxième surface d'émission, ladite deuxième pluralité de points de sonde formant une deuxième fauchée à la deuxième récurrence, 10 g) détermination d'au moins un point de sonde redondant appartenant aux première et deuxième fauchées, un point de sonde redondant ayant des coordonnées spatiales bidimensionnelles dans un plan horizontal identiques aux première et deuxième récurrences, h) calcul d'une variation d'erreur d'altitude du véhicule porteur en fonction 15 des coordonnées spatiales en altitude du au moins un point de sonde redondant calculées à l'étape c) pour la première récurrence et respectivement à l'étape f) pour la deuxième récurrence, et i) correction des coordonnées spatiales tridimensionnelles desdites première, et respectivement deuxième, fauchées, en fonction de la 20 variation d'erreur d'altitude du véhicule porteur calculée à l'étape précédente. DESCRIPTION DETAILLEE D'UN EXEMPLE DE REALISATION La description qui va suivre en regard des dessins annexés, donnés à titre d'exemples non limitatifs, fera bien comprendre en quoi consiste l'invention et 25 comment elle peut être réalisée. Sur les dessins annexés : - la figure 1 représente en vue de dessus un système bathymétrique, de type sonar multi-faisceaux ou lidar, de l'art antérieur formant, à deux récurrences distinctes, des fauchées perpendiculaires à l'axe de déplacement horizontal du 30 véhicule porteur ; - la figure 2 représente en vue de dessus un système bathymétrique, de type sonar multi-faisceaux ou lidar, selon un mode de réalisation de l'invention ; - la figure 3 représente une autre vue de dessus du système bathymétrique, de type sonar multi-faisceaux ou lidar, selon un mode de 3043788 9 réalisation de l'invention ; - la figure 4 illustre schématiquement en vue de dessus les zones mesurées par le système bathymétrique, de type sonar multi-faisceaux ou lidar, selon un mode de réalisation de l'invention. 5 Dispositif Sur les figures 1 à 4, on a représenté un repère orthonormé (X, Y, Z) local où le plan XY est un plan horizontal et l'axe Z un axe vertical. Sur la figure 1, on a représenté, en projection sur un plan horizontal, le fonctionnement d'un système bathymétrique classique de type sonar 10 multifaisceaux ou lidar selon l'art antérieur. A titre d'exemple, on détaille la structure et le fonctionnement d'un système bathymétrique de type sonar multifaisceaux. Toutefois, les mêmes principes s'appliquent de manière analogue à un système bathymétrique de type lidar. 15 Le système bathymétrique comprend au moins un émetteur acoustique et une antenne de récepteurs acoustiques, par exemple de type linéaire, disposée transversalement à la ligne de foi 5 du bateau 1. L'émetteur acoustique et l'antenne de récepteurs acoustiques sont attachés à un véhicule porteur, par exemple fixés sous la coque d'un bateau 1 ayant une ligne de foi 5, c'est à dire 20 une ligne de référence pointant vers la proue du bateau et correspondant à l'axe longitudinal du bateau. Lors d'une première récurrence, le bateau 1 est à une première position repérée A, l'émetteur acoustique émet un signal acoustique vers le fond marin suivant une première surface d'émission et l'antenne de récepteurs acoustiques 25 détecte l'amplitude des signaux acoustiques réfléchis dans des directions correspondant à chaque récepteur acoustique. La surface d'émission intersecte le fond marin suivant une figure géométrique, représentée de manière schématique par un segment de droite 10 sur la figure 1. A titre d'exemple, le segment de droite 10 a une longueur L et une largeur WL. La largeur WL correspond en général à la 30 résolution longitudinale le long de la ligne de foi 5 du bateau. La longueur L et la largeur WL sont fonction de multiples paramètres, entre autres de l'ouverture du faisceau incident, du nombre de récepteurs acoustiques et de la profondeur du fond marin. Lors de la première récurrence, l'antenne de récepteurs acoustiques multi-faisceaux permet de mesurer l'amplitude d'une première pluralité de 3043 788 10 faisceaux acoustiques réfléchis suivant n directions. Le calculateur détermine, à partir de l'amplitude de la première pluralité de faisceaux acoustiques réfléchis détectés, les coordonnées spatiales tridimensionnelles d'un ensemble de points de sonde sur le segment de droite, avec une résolution transverse WT. On obtient 5 ainsi les coordonnées spatiales tridimensionnelles (x, y, z) de M points de sonde 1001, 1002, 1003, ..., 1099 sur le segment de droite, où M un nombre entier naturel est égal à la longueur L divisée par la résolution spatiale transverse WT. M est en général choisi égal à 64, 128 ou 256 selon l'antenne acoustique. Les M points de sonde du fond marin forment une première fauchée 10, à la première 10 récurrence. On collecte ainsi une première fauchée 10 dans l'axe de l'antenne linéaire de récepteurs acoustiques. Le bateau se déplace d'une distance C suivant l'axe de déplacement 2 et se trouve à une deuxième position repérée B à un deuxième instant de récurrence. On effectue de manière analogue, lors de la deuxième récurrence, 15 une deuxième fauchée 20, constituée de M points de sonde 2001, 2002, ..., 2099. La deuxième fauchée 20 a généralement la même résolution spatiale longitudinale et respectivement transverse que la première fauchée 11. En projection dans un plan horizontal, la deuxième fauchée 20 est ici parallèle à la première fauchée 10. On constate que, sous réserve d'un avancement significatif du bateau, la 20 distance C étant supérieure à la résolution spatiale longitudinale WL, les fauchées 10 et 20 effectuées à deux instants de récurrence différents n'ont aucun point de sonde en commun. Les informations de ces deux fauchées 10, 20 étant indépendantes, il n'est généralement pas possible de les combiner pour estimer une variation d'altitude du véhicule porteur entre ces deux instants de récurrence 25 sans faire d'hypothèses sur les fréquences des mouvements du véhicule porteur par rapport aux fréquences de variation de l'altitude du fond marin le long du parcours du véhicule porteur. Un aspect de l'invention repose sur l'acquisition de mesure des sondes de telle façon que les fauchées effectuées à des instants de récurrence différents 30 se superposent en partie spatialement, ce qui permet de déterminer précisément les variations d'erreur d'altitude du véhicule porteur entre ces différents instants de récurrence. Sur la figure 2, on représente schématiquement en vue de dessus le fonctionnement d'un système bathymétrique selon un premier mode de 3043 788 11 réalisation. Le système comporte ici une première antenne linéaire et une deuxième antenne linéaire, orientées de manière à former un angle ALPHA, respectivement BETA avec la ligne de foi 5 du véhicule porteur 1. De préférence, les angles ALPHA et BETA sont opposés, de manière à ce que les deux antennes 5 linéaires soient disposées symétriquement par rapport à la ligne de foi 5 du véhicule porteur 1. Par exemple, la première antenne est alignée sur l'axe X et la deuxième antenne alignée sur l'axe Y du repère orthonormé (X, Y, Z). De façon avantageuse, les deux antennes linéaires ont la même longueur, le même nombre de récepteurs acoustiques et fournissent la même résolution spatiale longitudinale 10 WL, respectivement la même résolution spatiale transverse WT. A un premier instant de récurrence, le véhicule porteur 1 étant situé au point A, on collecte une première fauchée 11 suivant l'axe X de la première antenne. A titre d'exemple non limitatif, cette première fauchée est obtenue par l'émission d'un signal acoustique incident S1 sur la deuxième antenne et la 15 formation des voies en réception sur la première antenne, en utilisant les méthodes classiques de formation de voies et un filtre adapté au signal S. Le procédé de formation de voies, classique pour tous les sonars multi-faisceaux, permet de collecter M points de sondes le long de la fauchée déterminée par l'antenne d'émission avec en général M=64, 128 ou 256 en fonction des 20 spécifications de l'antenne. La première fauchée 11 comporte ainsi M points de sonde s'étendant suivant l'axe X. A un deuxième instant de récurrence, le véhicule porteur 1 étant situé au point B, on collecte une deuxième fauchée 22 suivant l'axe Y de la deuxième antenne. De manière symétrique, cette deuxième fauchée 22 est obtenue par 25 l'émission d'un autre signal acoustique incident S2 sur la première antenne et la formation des voies en réception sur la deuxième antenne. On obtient ainsi une deuxième fauchée 22 comportant aussi M points de sonde dans l'axe Y de la deuxième antenne. Le déplacement C du bateau 1 du point A au point B est tel que la 30 première fauchée intersecte la deuxième fauchée en un point de sonde redondant 31. Plus précisément, le point de sonde redondant 31 a les mêmes coordonnées en X et en Y dans un plan horizontal (XY) à la résolution des faisceaux acoustique près. Les faisceaux ayant une certaine largeur, chaque point de sonde correspond en fait à une petite zone de sonde de dimension WL*WT. Pour que deux points de 3043788 12 sonde soient redondants, il suffit que les deux zones correspondante s'intersectent. Toutefois, la première fauchée est obtenue à un premier instant de récurrence et la deuxième fauchée à un deuxième instant de récurrence, qui peut être séparé de plusieurs dizaines de secondes du premier instant de récurrence. 5 Par contre, la mesure de profondeur du point de sonde redondant 31 faisant partie de la première fauchée 11 au premier instant de récurrence est différente et indépendante de la mesure de profondeur du point de sonde redondant 31 faisant partie de la deuxième fauchée 22 au deuxième instant de récurrence, du fait des variations d'altitude et d'erreur d'altitude du véhicule porteur 1 entre la première et 10 la deuxième récurrence. Le calculateur utilise la différence de profondeur mesurée suivant l'axe Z de ce point de sonde redondant 31 entre le premier et le deuxième instant de récurrence pour estimer la variation d'altitude du véhicule porteur 1 entre ces deux instants de récurrence. On fait l'hypothèse que la première fauchée est à une 15 position de référence en Z. Le calculateur compense cette variation d'altitude pour toutes les mesures de points de sonde de la deuxième fauchée. On obtient ainsi une topographie du fond marin de meilleure précision en altitude. En fonction de la configuration des émissions-réceptions du sonar multifaisceaux, on peut former des fauchées de telle façon à avoir un ou plusieurs 20 recouvrements partiels entre des fauchés obtenues à des instants de récurrence séparés temporellement de plusieurs secondes. Différentes variantes sont considérées ici pour déterminer des fauchées ayant un ou plusieurs points de sonde redondants à des instants de récurrence séparés. 25 Selon une variante, de manière alternée, aux récurrences paires, on émet un signal acoustique sur la deuxième antenne et on reçoit des signaux acoustiques réfléchis en formant des voies sur la première antenne, et, aux récurrences impaires, on émet un autre signal acoustique sur la première antenne et on reçoit d'autres signaux acoustiques réfléchis en formant des voies sur la 30 deuxième antenne. Selon une autre variante, lors d'une récurrence sur K récurrences, où K est un nombre entier naturel, on émet un signal acoustique sur la deuxième antenne et on reçoit des signaux acoustiques réfléchis en formant des voies sur la première antenne, et aux K-1 récurrences suivantes, on émet un autre signal 3043788 13 acoustique sur la première antenne et on reçoit des signaux acoustiques réfléchis en formant des voies sur la deuxième antenne. De préférence dans ce cas, la valeur de K est limitée, par exemple K compris entre 3 et 10 afin de conserver les propriétés de redondance entre fauchées non consécutives dans le temps. 5 Selon encore une autre variante, illustrée sur la figure 2, lors de la première récurrence, on émet un signal acoustique S2 sur la deuxième antenne et on reçoit des signaux acoustiques réfléchis par réflexion du signal S2 sur le fond marin en formant des voies sur la première antenne, et simultanément, lors de la première récurrence, on émet un autre signal acoustique 51 sur la première 10 antenne et on reçoit des signaux acoustiques réfléchis par réflexion du signal 51 sur le fond marin en formant des voies sur la deuxième antenne. On procède de même lors de chaque récurrence. On choisit de préférence les signaux acoustiques 51 et S2 de telle façon à ce qu'ils présentent une inter-corrélation proche de zéro ce qui permet en réception de les distinguer puisque le filtre 15 adapté d'un des signaux appliqué à l'autre donne une valeur de sortie proche de zéro. Par exemple 51 et S2 occupent des bandes de fréquences acoustiques distinctes. Dans ce cas, on acquiert, à la première récurrence, une première fauchée 11 suivant l'axe X et une première fauchée 12 suivant l'axe Y, puis on acquiert, à la deuxième récurrence, une deuxième fauchée 21 suivant l'axe X et 20 une deuxième fauchée 22 suivant l'axe Y. La première fauchée 11 suivant l'axe X présente un point de sonde redondant 31 avec la deuxième fauchée 22 suivant l'axe Y. De plus, la première fauchée 12 suivant l'axe Y présente un autre point de sonde redondant 32 avec la deuxième fauchée 21 suivant l'axe X. On obtient ainsi deux points de sonde redondants 31, 32 distincts entre ces deux instants de 25 récurrence. Par contre, les autres points de sonde de ces deux fauchées ne sont pas redondants entre eux. Dans un exemple de réalisation d'un système de bathymétrie illustré sur la figure 3, les antennes sont orientées à 45 degrés de la ligne de foi du bateau 1, l'émission et la réception des signaux acoustiques sont simultanées sur les deux 30 antennes, l'ouverture des faisceaux acoustiques est de ±60 degrés, et la distance du fond environ de 150 mètres. L'empreinte des faisceaux acoustiques émis à un instant donné sur le fond marin est composée de deux segments orthogonaux ayant chacun une longueur L de 260 mètres et se croisant à la verticale du bateau. En calculant les 3043788 14 distances, comme illustré sur la figure 3, on observe que les fauchées 21, 22 obtenues à un deuxième instant de récurrence continuent à intersecter aux points de sonde 31, 32 situés aux extrémités des fauchées 11, 12 émises obtenues au premier instant de récurrence, lorsque le déplacement C du bateau reste inférieur 5 ou égal à 180 m= V2.2602 entre ces deux instants de récurrence. Avec une vitesse de déplacement du bateau de l'ordre de 4 noeuds, on en déduit que les fauchées 21, 22 intersectent avec les fauchées 11, 12 pendant 90 s. Cette durée est largement suffisante pour résoudre le problème d'estimation des variations d'erreur d'altitude du véhicule porteur à court terme. 10 On constate que ces propriétés de recouvrement sont également valables dans les deux autres schémas d'émission consistant à alterner l'utilisation des antennes en fonction des récurrences, avec alternance de récurrences paire et impaire ou avec alternance d'une récurrence sur K. Le choix du schéma d'émission est alors déterminé en fonction d'autres spécifications de 15 l'opération du système bathymétrique. La configuration dans laquelle les antennes forment un angle de 45 degrés avec la ligne de foi du bateau est la plus intéressante en termes de redondance des points de sonde, ce qui est utile en particulier pour estimer les autres paramètres du mouvement. Néanmoins d'autres configurations peuvent 20 être utilisées. Par exemple, on peut utiliser une configuration classique dans laquelle une des antennes est alignée avec l'axe longitudinal du bateau et l'autre antenne est disposée transversalement à cet axe. La durée du recouvrement spatial obtenue en croisant les fauchées est ainsi bien supérieure à la durée du recouvrement de deux fauchées parallèles 25 (comme illustré sur la figure 1), qui est limitée par la largeur WL d'une fauchée. Selon une variante, il peut être utile de déphaser les signaux d'émission acoustique des différents transducteurs de chaque antenne afin de dépointer le signal émis par rapport au plan orthogonal à l'antenne. Dans ce cas, les fauchées du sonar ne sont plus des segments de droites mais des segments d'hyperboles, 30 cependant on conserve en général la propriété de redondance partielle des points de sonde. La fréquence des signaux acoustiques utilisés est généralement comprise 5 kHz et 1 MHz. La fréquence acoustique peut être adaptée en fonction de la portée et de la résolution souhaitées. La fréquence acoustique détermine 3043788 15 principalement la portée et la résolution du système. L'invention s'applique indifféremment à n'importe quelle fréquence acoustique réfléchie par le fond marin et détectée par l'antenne multi-faisceaux. L'invention ne se limite nullement à la configuration à deux antennes telle 5 que décrite en lien avec les figures 2-4. De nombreuses autres configurations d'antenne acoustique permettent également d'assurer la propriété de redondance des points de sonde entre des instants de récurrence séparés temporellement. Par exemple, une antenne multi-capteurs bidimensionnelle permet aussi de collecter les positions de points de sonde redondants à des instants de 10 récurrence distincts. Une antenne bidimensionnelle est cependant plus coûteuse que deux antennes linéaires. Dans un autre exemple de réalisation, le système comporte une seule antenne de réception linéaire multi-capteurs et deux antennes d'émission monolithiques, une antenne d'émission émettant sous le bateau et l'autre antenne 15 d'émission vers l'avant du bateau. Cependant, cette configuration est moins avantageuse en termes de redondance. Comme indiqué plus haut, l'invention ne se limite pas à l'émission et la réception de signaux acoustiques. L'invention s'applique de manière analogue à l'émission et la réception de signaux optiques. Par exemple, on obtiendrait les 20 mêmes propriétés de redondance de points de sonde en combinant deux LIDAR à balayage de faisceau, chaque lidar comprenant une source laser, et un système de balayage de faisceau suivant un axe, de manière à balayer chaque faisceau laser dans un plan, les deux plans étant disposés orthogonalement l'un par rapport à l'autre. Les faisceaux émis par les deux sources laser n'interfèrent pas 25 optiquement entre eux du fait de leur géométrie. Alternativement, on considère aussi un système bathymétrique basé sur un LIDAR comprenant une seule source laser et un dispositif de balayage de faisceau suivant deux axes transverses. Un tel système lidar permet de collecter des points de sonde redondants à des instants de récurrence distincts et de calculer les variations d'erreur d'altitude du 30 véhicule porteur entre ces instants de récurrence. Procédé Dans la suite, nous détaillons un exemple de procédé de traitement des signaux entre récurrences distinctes, pour en extraire au moins un point redondant, et corriger les mesures de bathymétrie des mouvements du véhicule 3043788 16 porteur. Supposons que préalablement à la récurrence N on dispose d'une carte bathymétrique partielle du fond marin réalisée à partir des récurrences 1 à N-1 d'un système bathymétrique, basé sur un sonar multi-faisceaux ou sur un lidar tel 5 que décrit plus haut, par exemple en lien avec les figures 2 et 3. A la récurrence N, on obtient de nouveaux points de sondes. En première approximation, les coordonnées de ces points de sondes sont correctes à un décalage ou offset commun près, noté DZN, sur l'axe vertical. Etant donnée la configuration générale du système bathymétrique décrit plus haut, un certain 10 nombre M>=1 des points de sonde obtenus à la récurrence N appartiennent à la zone déjà cartographiée et forment donc M points de sonde redondants. On dispose donc de M équations linéaires sur DZN, c'est-à-dire une équation par point de sonde redondant avec un point de sonde de la zone déjà cartographiée. Il suffit donc de trouver la valeur de DZN qui approche le mieux la solution de ces M 15 équations. Par exemple on pose DZN égale à la moyenne des M' valeurs données par le sous-ensemble des M équations obtenu en retirant les points aberrants éventuels. On observe que si M est égal à 1, on ne peut pas retirer de point aberrant. Par contre plus M est grand, plus il est possible de retirer un nombre 20 important de points aberrants en utilisant une méthode statistique appropriée, par exemple une méthode de vraisemblance. En pratique, le nombre M est en général très significativement supérieur à 1. Généralement, dans un sonar multi-faisceaux, un segment de fauchée correspond à de nombreux faisceaux acoustiques, allant de quelques dizaines à quelques centaines de faisceaux acoustiques comme 25 indiqué plus haut, et donc à de nombreuses sondes. Si la configuration d'antenne est bien choisie, le nombre de points de sondes redondants est très important. Sur la figure 4, on a représenté un sonar comportant deux antennes d'émission-réception disposées symétriquement à 45 degrés de la ligne de foi 5 du bateau 1, comme décrit en lien avec les figures 2-3, à un instant de récurrence N. Le sonar 30 est par exemple configuré pour obtenir à chaque récurrence une fauchée s'étendant le long de l'axe X et une autre fauchée s'étendant le long de l'axe Y. Plus précisément, la fauchée suivant l'axe X comporte une moitié de segment de fauchée 13 et une autre moitié de segment de fauchée 11. De manière analogue, la fauchée suivant l'axe Y comporte, de part et d'autre d'un plan vertical passant 3043788 17 par la ligne de foi du bateau, une moitié de segment de fauchée 12 et une autre moitié de segment de fauchée 14. Le bateau se déplace suivant l'axe de déplacement 2. A la récurrence N, la zone 60 du fond marin située à l'avant du bateau n'a pas encore été sondée. Les zones 41; 42 et 50 situées le long de la 5 trace de déplacement du bateau, ont déjà été sondées lors des récurrences 1 à N1. Plus précisément, la zone 41 a été cartographiée par l'antenne s'étendant suivant l'axe X, la zone 42 a été cartographiée par l'antenne s'étendant suivant l'axe Y et la zone 50, située à l'arrière du bateau, a déjà été cartographiée par ces deux antennes à différentes récurrences. Ainsi, à la récurrence N, la moitié de 10 segment de fauchée 13 suivant l'axe X forme des points de sonde potentiellement redondants avec des points de sonde obtenus à une récurrence antérieure sur une moitié de l'antenne s'étendant suivant l'axe Y. De manière analogue, à la récurrence N, la moitié de segment de fauchée 14 suivant l'axe Y forme des points de sonde potentiellement redondants avec des points de sonde obtenus à une 15 récurrence antérieure sur une moitié de l'antenne s'étendant suivant l'axe X. A la récurrence N, les points redondants peuvent former jusqu'à la moitié du nombre total de points de sonde obtenus, ce qui fournit un nombre de points de sonde redondants de l'ordre de 30 à 250, à la récurrence N. Une fois la valeur de DZN estimée, on complète la carte bathymétrique 20 du fond marin en utilisant les points de sonde obtenus à la récurrence N, corrigés de l'erreur d'altitude DZN. On observe qu'une fois obtenue une cartographie partielle du fond marin, il est facile de la compléter en éliminant les effets des variations d'erreur d'altitude du véhicule porteur du système de mesure. 25 Le problème de recalage se ramène donc à obtenir une carte initiale avant de mettre en oeuvre le procédé décrit ci-dessus. Il existe pour cela de nombreuses possibilités plus ou moins complexes suivant la qualité désirée pour l'estimation initiale de l'altitude. La plus simple consiste à prendre comme carte initiale la carte, très partielle mais suffisante, obtenue après la première sonde à 30 partir des données d'altitude disponibles à la première récurrence. Plus généralement on peut prendre comme carte initiale la carte obtenue à partir des K premières récurrences en utilisant des estimations approchées de l'altitude du porteur. Dans le cas où le véhicule porteur est un bateau, on connait son altitude approchée puisque le bateau est par définition à la surface de la mer. 3043788 18 Selon une variante, le système bathymétrique comporte en outre un capteur d'altitude externe. La méthode de correction des mouvements en altitude du véhicule porteur décrite ci-dessus se combine avantageusement avec l'utilisation d'un 5 capteur d'altitude externe connecté au système bathymétrique. Par exemple, pour un bateau, un système GPS fournit l'altitude du bateau. Autre exemple, pour un système sous-marin, un capteur de pression fournit une mesure de la profondeur d'immersion du système d'émission-réception multi-faisceaux. Cette combinaison d'un système bathymétrique multi-faisceaux et d'un capteur externe permet 10 d'améliorer progressivement la qualité de la carte bathymétrique en profitant de la bonne précision moyenne du capteur externe sur le long terme. Pour combiner efficacement les mesures provenant du système bathymétrique multi-faisceaux et du capteur externe, on adapte l'étape d'estimation de la valeur du décalage DZN de la façon suivante. 15 A chaque récurrence, on pose DZN égal au barycentre de la moyenne des M' points de sondes redondants avec la carte et de la valeur ZN fournie par le capteur externe, les coefficients barycentriques étant déterminé à partir des incertitudes du capteur d'altitude externe d'une part, des points de sonde de la carte d'autre part. 20 A titre d'exemple, on pose : DZN = 1/ (SIGMA2_CARTE / M' + SIGMA2_ZN) * ((SIGMA2_CARTE / M') * ZN + (SIGMA2_ZN) * MOYENNE) où SIGMA2_ZN est la variance de l'erreur de mesure du capteur d'altitude externe, SIGMA2_CARTE est la variance de l'erreur d'altitude des points de la 25 carte et MOYENNE est la moyenne des M' valeurs d'erreurs d'altitude données par le sous-ensemble des M équations de sondes redondantes obtenu en retirant les points aberrants. Selon une variante du procédé, une estimation des autres paramètres du mouvement en translation et/ou en rotation du véhicule porteur est déduite des 30 mesures de points de sonde redondants à des récurrences distinctes. Dans ce cas, on peut utiliser la même approche en modifiant légèrement l'algorithme pour estimer les variations des autres paramètres du mouvement entre différents instants de mesure des points de sonde. En particulier on peut utiliser la même approche pour estimer les variations d'orientation (cap, roulis, 3043788 19 tangage) du système bathymétrique multi-faisceaux. En effet chaque variation d'un de ces paramètres entre deux récurrences se traduit par une formule d'erreur définie sur la position géographique estimée des points de sonde en fonction de leur position relative par rapport au système de mesure bathymétrique multi- 5 faisceaux. Or la position relative du système de mesure bathymétrique est connue précisément par le principe de mesure du sonar multifaisceaux. Par exemple, une erreur de roulis de R radians entraine pour un point à une distance D par rapport à l'axe principal du système de mesure, une erreur d'altitude en Z de la sonde égale à R * D. 10 En profitant des redondances multiples entre les points de sondes calculées à la récurrence N, avant l'étape de correction, et de la carte déterminée à partir des récurrences précédentes, on peut non seulement déterminer la valeur de DZN mais aussi celle de R. Cette estimation se fait par exemple par une méthode de type moindres 15 carrés dans laquelle on cherche à minimiser les distances entre les nouveaux points de sondes redondants et la carte préétablie (les méthodes de moindres carrées sont bien connues de l'homme de l'art ; les algorithmes donnés plus haut relèvent d'ailleurs également de la méthode des moindres carrés mais dans des cas plus simples). 20 Un procédé de correction analogue est applicable aux variations de cap et/ou de tangage. On observe que la qualité de ces estimations dépend fortement du recouvrement spatial qui est obtenu entre des fauchées obtenues à des instants de récurrence différents. En pratique, la configuration avec les deux antennes à 45 25 degrés de l'axe longitudinal du bateau donne les meilleures conditions d'observabilité des erreurs et donc les meilleures estimations. Là encore, de façon similaire à ce qui est indiqué plus haut concernant un capteur d'altitude externe, la méthode de correction peut être combinée de façon optimale avec l'utilisation d'un capteur externe, par exemple de type système de 30 navigation inertielle, fournissant des estimations des variations de cap, roulis, tangage et déplacement. REVENDICATIONS1. Système bathymétrique comprenant un calculateur et un dispositif d'émission et de réception multi-faisceaux attaché à un véhicule porteur (1), le dispositif d'émission et de réception multi-faisceaux étant configuré à une première récurrence pour émettre un premier faisceau incident vers le fond marin suivant une première surface d'émission et pour mesurer l'amplitude d'une première pluralité de faisceaux réfléchis par réflexion du premier faisceau incident sur le fond marin, le calculateur étant configuré pour déterminer, à partir de l'amplitude de la première pluralité de faisceaux réfléchis détectés, les coordonnées spatiales tridimensionnelles d'une première pluralité de points de sonde (1101, 1102, ..., 1199) situés à l'intersection du fond marin et de la première surface d'émission, ladite première pluralité de points de sonde (1101, 1102, ..., 1199) formant une première fauchée (11) à la première récurrence, caractérisé en ce que : - le dispositif d'émission et de réception multi-faisceaux est configuré pour émettre, à une deuxième récurrence, un deuxième faisceau incident vers le fond marin suivant une deuxième surface d'émission, la deuxième surface d'émission étant non parallèle à la première surface d'émission, et les premier et deuxième faisceaux incidents étant non-interférents entre eux, - le dispositif d'émission et de réception multi-faisceaux est configuré pour mesurer l'amplitude, à la deuxième récurrence, d'une deuxième pluralité de faisceaux réfléchis par réflexion du deuxième faisceau incident sur le fond marin, - le calculateur étant configuré pour déterminer, à partir de l'amplitude de la deuxième pluralité de faisceaux réfléchis détectés, les coordonnées spatiales tridimensionnelles d'une deuxième pluralité de points de sonde (2201, 2202, ..., 2299) situés à l'intersection du fond marin et de la deuxième surface d'émission, ladite deuxième pluralité de points de sonde formant une deuxième fauchée (22), de manière à ce que au moins un point de sonde de la première fauchée (11) soit redondant avec un point de sonde de la deuxième fauchée (22), un point de sonde redondant (31) ayant des coordonnées spatiales bidimensionnelles dans un plan horizontal identiques aux première et deuxième récurrence, et - le calculateur étant configuré pour calculer une variation d'erreur d'altitude du véhicule porteur à partir des variations de coordonnées spatiales en altitude du au 3043788 21 moins un point de sonde redondant (31) entre la première récurrence et la deuxième récurrence, et - le calculateur étant configuré pour corriger les coordonnées spatiales tridimensionnelles desdites première, et respectivement deuxième, fauchées (11, 5 22), en fonction de la variation d'erreur d'altitude du véhicule porteur (1) calculée entre lesdites récurrences. 2. Système bathymétrique selon la revendication 1, dans lequel le dispositif d'émission et de réception multi-faisceaux comporte un sonar multifaisceaux ayant une première et une deuxième antennes linéaires d'émission et 10 de réception multi-faisceaux, la première antenne linéaire étant disposée transversalement par rapport à la deuxième antenne linéaire. 3. Système bathymétrique selon la revendication 2, dans lequel la première antenne est configurée pour émettre le premier faisceau incident, la deuxième antenne est configurée pour détecter la première fauchée à la première 15 récurrence, et dans lequel la deuxième antenne est configurée pour émettre le deuxième faisceau incident, et la première antenne est configurée pour détecter la deuxième fauchée à la deuxième récurrence. 4. Système bathymétrique selon l'une des revendications 2 ou 3, dans lequel la première antenne linéaire et la deuxième antenne linéaire sont disposées 20 symétriquement à environ 45 degrés d'une ligne de foi (5) du véhicule porteur (1). 5. Système bathymétrique selon l'une des revendications 2 ou 3, dans lequel la première antenne linéaire est alignée sur une ligne de foi du véhicule porteur. 6. Système bathymétrique selon la revendication 1, dans lequel le 25 dispositif d'émission et de réception multi-faisceaux comporte un sonar multi- faisceaux ayant au moins une antenne de réception multi-capteurs à deux dimensions. 7. Système bathymétrique selon la revendication 6, comprenant deux antennes d'émission et une antenne de réception linéaire multi-capteurs. 30 8. Système bathymétrique selon la revendication 1, dans lequel le dispositif d'émission et de réception multi-faisceaux comporte un système lidar configuré pour émettre un faisceau laser vers le fond marin suivant une première et respectivement une deuxième surface d'émission, et pour détecter une pluralité de faisceaux laser réfléchis sur le fond marin, formant la première fauchée (11) et 3043788 22 respectivement la deuxième fauchée (22). 9. Système bathymétrique selon la revendication 8, dans lequel le système lidar comporte une source laser à balayage bidirectionnel de faisceau ou deux sources laser à balayage unidirectionnel de faisceau. 5 10. Système bathymétrique selon l'une des revendications 1 à 9, comportant en outre au moins un capteur externe choisi parmi un système de géolocalisation par satellite, un capteur de profondeur ou une centrale de navigation inertielle, et dans lequel le capteur externe est configuré pour fournir au moins une mesure complémentaire de position et/ou d'orientation du dispositif 10 d'émission et de réception multi-faisceaux, et dans lequel le calculateur est configuré pour calculer une variation d'altitude du véhicule porteur à partir des variations de coordonnées spatiales en altitude du au moins un point de sonde redondant (31) entre la première récurrence et la deuxième récurrence et en fonction de ladite au moins une mesure complémentaire de position et/ou 15 d'orientation. 11. Méthode de bathymétrie comprenant les étapes suivantes : a) émission, à une première récurrence, par un dispositif d'émission et de réception multi-faisceaux, d'un premier faisceau incident apte à se propager dans un milieu sous-marin vers un fond marin 20 suivant une première surface d'émission, b) détection par le dispositif d'émission et de réception multifaisceaux, à la première récurrence, de l'amplitude d'une première pluralité de faisceaux réfléchis suivant une première pluralité de directions par réflexion du premier faisceau incident 25 sur le fond marin, c) calcul, à partir de l'amplitude de la première pluralité de faisceaux réfléchis détectés, des coordonnées spatiales tridimensionnelles d'une première pluralité de points de sonde situés à l'intersection du fond marin et de la première surface d'émission, ladite 30 première pluralité de points de sonde formant une première fauchée à la première récurrence, d) émission, à une deuxième récurrence, par le dispositif d'émission et de réception multi-faisceaux, d'un deuxième faisceau incident apte à se propager dans le milieu sous-marin vers le fond marin 3043788 23 suivant une deuxième surface d'émission, la deuxième surface d'émission étant non parallèle à la première surface d'émission, le premier faisceau incident et le deuxième faisceau incident étant non interférents entre eux, 5 e) détection, par le dispositif d'émission et de réception multi- faisceaux, à la deuxième récurrence, de l'amplitude d'une deuxième pluralité de faisceaux réfléchis, suivant une deuxième pluralité de directions, par réflexion du deuxième faisceau incident sur le fond marin, 10 f) calcul, à partir de l'amplitude de la deuxième pluralité de faisceaux réfléchis détectés, des coordonnées spatiales tridimensionnelles d'une deuxième pluralité de points de sonde situés à l'intersection du fond marin et de la deuxième surface d'émission, ladite deuxième pluralité de points de sonde formant 15 une deuxième fauchée à la deuxième récurrence, g) détermination d'au moins un point de sonde redondant appartenant aux première et deuxième fauchées, un point de sonde redondant ayant des coordonnées spatiales bidimensionnelles dans un plan horizontal identiques aux 20 première et deuxième récurrences, h) calcul d'une variation d'erreur d'altitude du véhicule porteur en fonction des mesures d'altitude du au moins un point de sonde redondant aux première et deuxième récurrences, et i) correction des coordonnées spatiales tridimensionnelles desdites 25 première, et respectivement deuxième, fauchées, en fonction de la variation d'erreur d'altitude du véhicule porteur calculée à l'étape précédente.