FR 3142999 A1 20240614 FR 2213238 A 20221213 Procédé et dispositif de gestion d’une signature magnétique d’un navire, produit programme d'ordinateur et médium de stockage correspondants Procédé de gestion d’une signature magnétique d’un navire, exécuté par une machine de calcul et comprenant une étape (20) de modélisation magnétique du navire comprenant : obtention (21) de mesures de champs magnétiques associées chacune à des informations comprenant une position géographique, un cap magnétique du navire et une distance par rapport au navire ; et détermination (22) d’un modèle en moments magnétiques du navire, par résolution d’un premier problème inverse visant à déterminer des moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées. La résolution du premier problème inverse est réalisée avec une première méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. Le modèle en moments magnétiques du navire est utilisé dans une étape (23) de prédiction de la signature magnétique du navire et/ou dans une étape (24) de c ompensation de la signature magnétique du navire. Figure d’abrégé : figure 2 Procédé et dispositif de gestion d’une signature magnétique d’un navire, produit programme d'ordinateur et médium de stockage correspondants 1. DOMAINE TECHNIQUE Le domaine de l’invention est celui du magnétisme des navires (bâtiments navals). L’invention se rapporte plus particulièrement à un procédé et un dispositif de gestion d’une signature magnétique d’un navire. L’invention a ainsi de nombreuses applications, notamment, mais non exclusivement, la prédiction de la signature d’un navire et la compensation de la signature d’un navire. De façon générale, la signature magnétique d'un navire de surface correspond au champ total sur une ligne ou un plan de référence situé sous la quille à une distance égale à la largeur du navire. Elle dépend de la position du navire sur le globe terrestre et de son orientation. Elle représente la déformation du champ magnétique terrestre par l'aimantation du navire. La connaissance et la réduction de la signature magnétique d’un navire (notamment militaire) est cruciale pour assurer sa sécurité en approche des côtes (zone potentiellement minée) et, dans le cas particulier d’un sous-marin, c’est sa discrétion qui est directement impactée (détection par un aéronef/drone). Il est donc primordial d’avoir un moyen de mesurer cette signature magnétique. Cette mesure n’est cependant pas suffisante et un passage par une modélisation est indispensable pour pouvoir prédire la signature magnétique à un endroit (position géographique) différent du globe et/ou à une distance du navire différente et/ou à une attitude du navire différente, la signature magnétique dépendant entre autres du champ magnétique terrestre local et de l’orientation du navire dans celui-ci. L’attitude du navire est définie par le roulis, le tangage et le cap magnétique, qui correspondent à trois rotations autour des trois axes orthogonaux conventionnels définis pour le navire. Plus précisément, les coordonnées du navire sont rapportées à un trièdre orthogonal direct dont l’axe des x est orienté vers la proue, l’axe des y vers tribord et l’axe des z vers le bas. Dans un tel système, le roulis correspond à des rotations autour de l’axe des x (positives quand le navire gîte vers tribord), le tangage correspond à des rotations autour de l’axe des y (positives quand la proue pointe vers le haut) et les variations de cap (aussi appelées « lacets ») correspondent à des rotations autour de l’axe des z (positives dans le sens des aiguilles d’une montre). Le cap magnétique (aussi appelé cap géomagnétique) d’un navire est l’angle exprimé en degrés (de 0 à 360°), dans le sens des aiguilles d’une montre, entre la direction du nord magnétique et la ligne de foi du navire (c'est-à-dire la direction dans laquelle pointe l’étrave du navire, selon l’axe longitudinal du navire). Le cap magnétique d’un navire peut être obtenu soit directement (par exemple avec un magnétomètre), soit indirectement (par exemple par calcul à partir du cap compas ou du cap géographique). Le cap compas d’un navire est le cap indiqué par le compas, c'est-à-dire l’angle entre le nord du compas et la ligne de foi du navire. La différence angulaire entre le nord du compas et le nord magnétique est appelée déviation. Le cap géographique (aussi appelé cap vrai) d’un navire est l’angle entre le nord géographique (aussi appelé nord vrai) et la ligne de foi du navire. La différence angulaire entre le nord géographique et le nord magnétique est appelée déclinaison. Au-delà de la prédiction de la signature magnétique du navire, la connaissance de celle-ci permet de déterminer quelles actions sont à mener pour la réduire. Ceci peut se faire notamment par un système d’immunisation (aussi appelé « système de degaussing ») embarqué à bord du navire et constitué de boucles de courants. Pour assurer la réduction de la signature magnétique, il faut régler le système d’immunisation (en général ceci est effectué sur une station de mesure) et, là aussi, par modélisation, distribuer correctement les courants des boucles (aussi appelés « courants de degaussing ») en fonction de la signature prédite (fonction des limites de localisation (champ terrestre), de distance et d’attitude). 2. ARRIÈRE-PLAN TECHNOLOGIQUE Les techniques utilisées actuellement permettent de déterminer un modèle magnétique du navire en moments magnétiques, par résolution d’un problème inverse (c'est-à-dire un problème mathématique dit « inverse ») visant à déterminer des moments magnétiques du navire à partir de mesures de champs magnétiques et d’informations associées (position géographique du capteur effectuant la mesure, cap magnétique du navire et une distance par rapport au navire ). Une telle détermination d’un modèle magnétique du navire par résolution d’un problème inverse (aussi appelée « modélisation inverse » ou encore « modélisation selon le modèle inverse ») est décrite par exemple dans la demande de brevet WO2017064432A1. Mais les méthodes mathématiques employées (notamment la méthode des moindres carrés) ne mènent pas à une solution unique car le problème est mal posé (de petites variations dans les champs mesurés peuvent changer drastiquement la solution calculée). En effet, la philosophie de modélisation du navire et de réglage du système de degaussing et la logique d’interprétation, employées depuis des dizaines d’années par les experts dans le domaine du magnétisme du navire, n’ont jamais été explicitement traduites mathématiquement dans le domaine des moments magnétiques pour la résolution du problème de modélisation et de réglage du système de degaussing. En effet, le navire est communément modélisé par un plan de moments magnétiques (dipôles ; sources élémentaires de magnétisme) et les méthodes mathématiques employées tentent de trouver la valeur de chacun de ces moments. Plus précisément, dans l’état de l’art, les méthodes habituelles cherchent à trouver la valeur et la direction de chacune des sources également réparties en minimisant un critère entre la mesure réelle de champ magnétique et la mesure reconstruite, en considérant que toutes les sources potentielles ont une existence (c'est-à-dire une valeur non nulle). La solution trouvée est ainsi correcte mathématiquement mais a très peu de chance d’avoir une signification physique (dépend beaucoup du type de régularisation employée pour la résolution du problème inverse). Or, il n’y a qu’une solution qui satisfasse la physique et ce problème impose de repasser dans le domaine des champs magnétiques pour prédire la signature et régler le système de degaussing. Il y a donc un risque de mauvaise interprétation du champ magnétique (notamment, il y a un risque de compenser un moment magnétique selon un axe par des boucles de degaussing selon les autres axes) et ceci réserve l’exploitation de la prédiction et du réglage du système de degaussing à une poignée d’experts à travers le monde. Sans expert, et c’est notamment le cas pour exploiter une station de mesure portable embarquée à bord du navire, il n’est plus possible de correctement exploiter les mesures. En résumé, les modélisations actuelles souffrent de donner une solution mathématiquement correcte mais cette solution n’étant pas unique, elle n’a aucune raison d’être proche de la réalité physique et cela contraint la capacité à prédire la signature magnétique du navire et le risque de détection à une zone proche de la mesure initiale et, pour le réglage du système de degaussing, ne permet pas un réglage automatique simple ni un réglage en dehors de cette même zone. L’utilisation de méthodes mathématiques différentes (méthodes bayésiennes, notamment) a commencé depuis quelques années, mais pour le moment, ces méthodes venant de la statistique ont été exploitées juste pour le complément d’information qu’elles apportent (écart type, estimation de l’erreur,…). Les auteurs continuent à chercher les moments magnétiques également répartis dans le plan de modélisation et à repasser dans le domaine des champs magnétiques pour la prédiction de signature et le réglage du système de degaussing. Afin d’illustrer le fait que la modélisation, la prédiction de signature et le réglage des courants de degaussing passent aujourd’hui inexorablement par le domaine des champs magnétiques, on présente maintenant, en relation avec la , un procédé de gestion d’une signature magnétique d’un navire selon un mode de réalisation particulier de l'art antérieur. Dans une étape 11, pour un cap magnétique donné, le dispositif (exécutant le procédé) obtient des mesures de champs magnétiques associées chacune à des informations comprenant une position géographique, le cap magnétique donné du navire et une distance par rapport au navire. Dans une étape 12, pour le cap magnétique donné précité (voir étape 11), le dispositif détermine un modèle magnétique du navire en moments magnétiques, par résolution d’un problème inverse visant à déterminer des moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées. Dans une étape 13, pour le cap magnétique donné précité (voir étapes 11 et 12), le dispositif détermine une signature magnétique du navire, en fonction du modèle magnétique du navire. En d’autres termes, le dispositif repasse dans le domaine des champs magnétiques. La signature magnétique déterminée à l’étape 13 est valable uniquement pour le jeu d’informations associées aux mesures de champs magnétiques (même position géographique, même cap magnétique du navire et même distance par rapport au navire). Comme indiqué par le cadre référencé 10, les étapes 11, 12 et 13 sont itérées pour chacun d’une pluralité de N caps magnétiques (généralement on considère les quatre (N = 4) caps magnétiques particuliers suivants : cap magnétique nord (quand le cap magnétique est égal à 0°), cap magnétique sud (quand le cap magnétique est égal à 180°), cap magnétique est (quand le cap magnétique est égal à 90°) et cap magnétique ouest (quand le cap magnétique est égal à 270°)). On obtient ainsi une signature magnétique du navire pour chacun de cette pluralité de caps magnétiques. En d’autres termes, on obtient une « signature à différents caps magnétiques » (si N est supérieur à 1). Dans une étape 14, le dispositif effectue une prédiction de la signature magnétique du navire à une position géographique et/ou un cap magnétique du navire et/ou une distance par rapport au navire autre(s) que celle, celui ou ceux des mesures de champs magnétiques, en fonction de la signature magnétique (domaine des champs magnétiques) pour différents caps magnétiques (résultat des N itérations de l’étape 13 pour les N différents caps magnétiques) et à l’aide d’abaques et d’estimateurs construits par expérimentations. Dans une étape 15, le dispositif effectue une compensation de la signature magnétique du navire à l’aide du système de degaussing (système d’immunisation comprenant des boucles de courant installées à bord du navire). L’étape de compensation 15 comprend elle-même les étapes suivantes : détermination à l’étape 16 de l’effet magnétique, dans le domaine des champs magnétiques, des boucles de courant ; pour chacune des boucles de courant, détermination à l’étape 17 d’un courant de compensation en fonction de la signature magnétique du navire (celle déterminée à différents caps magnétiques à l’étape 13, ou bien celle prédite à l’étape 14 ; le second cas est symbolisé par la flèche référencée 19) et de l’effet magnétique (en champs magnétiques) des boucles de courant. En d’autres termes, l’étape 17 comprend des opérations d’analyse et de calcul effectuées dans le domaine des champs magnétiques ; et application à l’étape 18, sur les boucles de courant, des courants de compensation déterminés. Pour plus de détail concernant la compensation de la signature magnétique d’un navire à l’aide d’un système de degaussing, on pourra par exemple se reporter à la demande de brevet WO2017064432A1 déjà mentionnée plus haut, ainsi qu’aux brevets EP0597014B1 et US6965505B1. 3. RÉSUMÉ Dans un mode de réalisation de l’invention, il est proposé un procédé de gestion d’une signature magnétique d’un navire, exécuté par une machine de calcul et comprenant une étape de modélisation magnétique du navire comprenant elle-même : obtention de mesures de champs magnétiques associées chacune à des informations comprenant une position géographique, un cap magnétique du navire et une distance par rapport au navire ; et détermination d’un modèle en moments magnétiques du navire, par résolution d’un premier problème inverse visant à déterminer des moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées ; la résolution du premier problème inverse étant réalisée avec une première méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. Ainsi, la solution proposée propose une approche tout à fait nouvelle et inventive. En effet, comme indiqué plus haut, la modélisation, la prédiction de signature et le réglage des courants de degaussing passent toujours, dans les solutions connues de l’art antérieur, par le domaine des champs magnétiques. Or le champ magnétique est une conséquence et non la source magnétique, le champ magnétique étant ce que l’on peut mesurer, mais il introduit des complexités et des biais. L’objet de l’invention est justement de pouvoir avoir accès à un modèle du navire en sources magnétiques, c'est-à-dire en moments magnétiques, exploitable pour faire des opérations directement dans ce domaine des moments magnétiques. L’enjeu est en effet qu’un moment magnétique a une direction donnée alors que le champ magnétique qui en résulte est selon les trois axes en fonction de l’endroit où l’on se place pour mesurer. L’objectif recherché étant de correctement connaître la direction de chacune des sources magnétiques pour prédire correctement le champ magnétique créé par ces sources et ainsi effectuer correctement la prédiction de signature et/ou le réglage des courants de degaussing. La solution proposée propose donc d’utiliser une méthode de régularisation parcimonieuse pour déterminer le modèle en moments magnétiques du navire. Du fait qu’elle utilise un a priori favorisant la parcimonie, cette méthode de régularisation parcimonieuse permet de contraindre la résolution du premier problème inverse (celui visant à déterminer des moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées) de sorte que le modèle magnétique du navire, obtenu en moments magnétiques, soit directement exploitable dans un dispositif utilisable par n’importe quelle personne à bord du navire. En d’autres termes, l’idée de l’invention est que la solution la plus simple est aussi la plus physique : rien ne sert de chercher à mettre absolument dans le modèle des sources magnétiques partout, également réparties, il vaut mieux chercher un nombre plus faible (préférentiellement le minimum) de sources qui expliquent les mesures et donc s’autoriser à annuler certaines sources (ceci est d’ailleurs en cohérence avec le principe de conception (design) d’un système de degaussing, quelques boucles donc moments magnétiques suffisant à compenser un navire). C’est là une différence essentielle avec l’état de l’art, puisque l’utilisation d’une méthode dite parcimonieuse vise justement à modéliser le navire par un nombre restreint de sources magnétiques (moments magnétiques). Selon une caractéristique particulière, la première méthode de résolution basée sur une régularisation parcimonieuse s’appuie sur au moins une des normes et . L’avantage de ces normes et est de permettre de trouver un nombre restreint de moments magnétiques (propriété de parcimonie). Dans la présente description, la norme du vecteur (qui s’écrit ) compte le nombre d’éléments non nuls du vecteur . On notera qu’il s’agit d’un abus de langage, et qu’en toute rigueur l’écriture correcte est « norme » , car c’est en réalité une pseudo-norme. La norme d’un vecteur , …, ) s’écrit : . Plus généralement, on rappelle que pour , la norme d’un vecteur , …, ) s’écrit : . Selon certains modes de réalisation, le premier problème inverse s’écrit sous une des formes mathématiques suivantes : avec un vecteur des moments magnétiques du navire à estimer, un vecteur des mesures de champs magnétiques et une matrice de coefficients dépendant des positions associées aux mesures de champs magnétiques. Cette liste de formes mathématiques n’est pas exhaustive. Dans une implémentation particulière, la première méthode de résolution basée sur une régularisation parcimonieuse est une méthode bayésienne. Un avantage de la méthode bayésienne est qu’elle donne un bon compromis entre vitesse de convergence et qualité d’estimation du vecteur , et fournit également l’incertitude sur les moments magnétiques estimés. Cependant la présente invention n’est pas limitée à ce type de méthode. Selon certains modes de réalisation, le procédé comprend une étape de prédiction de la signature magnétique du navire à une position géographique et/ou un cap magnétique du navire et/ou une distance par rapport au navire autre(s) que celle, celui ou ceux des mesures de champs magnétiques, en fonction du modèle en moments magnétiques du navire. Ainsi, le modèle en moments magnétiques du navire (obtenu à l’étape de modélisation) permet de prédire la signature magnétique de navire dans des conditions différentes de celles de la mesure des champs magnétiques (c'est-à-dire à un cap magnétique donné et/ou à une distance donnée et/ou à une localisation donnée autre(s) que ceux de la mesure), et ainsi de donner aux marins à bord du navire une estimation de leur distance de détection par un dispositif ennemi (distance d’approche, profondeur minimale de navigation ou profondeur pour un navire et direction de plongée pour un sous-marin). On notera que, contrairement à l’état de l’art, la solution proposée permet une prédiction avec un changement de distance. En effet, du fait que, dans l’invention, le modèle en moments magnétiques du navire résulte de l’utilisation d’une méthode de régularisation parcimonieuse, il est suffisamment proche de la réalité physique (contrairement aux modèles de navire obtenus avec les solutions de l’état de l’art) pour avoir une validité en dehors de la zone de mesure initiale (dans la limite de validité d’un système de degaussing ou, autrement dit, en dehors de la zone dite de champ proche). Selon certains modes de réalisation, le procédé comprend une étape de compensatio n de la signature magnétique du navire à l’aide d’un système d’immunisation comprenant des boucles de courant installées à bord du navire, l’étape de compensation de la signature magnétique du navire comprenant : détermination d’un modèle en moment(s) magnétique(s) de chacune des boucles de courant ; pour chacune des boucles de courant, détermination d’un courant de compensation en fonction du modèle en moments magnétiques du navire et des modèles en moments magnétiques des boucles de courant ; et application, sur les boucles de courant, des courants de compensation déterminés. Ainsi, le modèle en moments magnétiques du navire (obtenu à l’étape de modélisation) permet d’effectuer une compensation de la signature magnétique du navire à l’aide d’un système d’immunisation. Comme déjà mentionné plus haut, dans l’invention, le modèle en moments magnétiques du navire résulte de l’utilisation d’une méthode de régularisation parcimonieuse, il est suffisamment proche de la réalité physique (contrairement aux modèles de navire obtenus avec les solutions de l’état de l’art) pour avoir une validité en dehors de la zone de mesure initiale. Les moments magnétiques étant très directement liés aux courants électriques, et donc aux effets de boucles, ils offrent l’avantage de pouvoir régler les boucles de courants plus facilement et surtout automatiquement, rendant un opérateur lambda autonome. La compensation de la signature magnétique du navire peut donc être effectuée sous la forme d’un réglage automatique des courants de degaussing, avec des quantités physiques en « ampères x mètres² » qui correspondent directement aux caractéristiques des boucles de courant (boucles de compensation magnétique). In fine , la méthodologie fournit à l’opérateur des courants (en ampères) à appliquer sur les boucles de compensation existantes. Selon une caractéristique particulière, le calcul des courants de compensation est réalisé par résolution d’un deuxième problème inverse, visant à déterminer les courants de compensation en fonction du modèle en moments magnétiques du navire et des modèles en moments magnétiques des boucles de courant, avec une deuxième méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. Ainsi, il est proposé d’utiliser une méthode de régularisation parcimonieuse également pour déterminer les courants de compensation (courants de degaussing). Du fait qu’elle utilise un a priori favorisant la parcimonie, une telle méthode de régularisation parcimonieuse permet de contraindre la résolution du deuxième problème inverse (celui visant à déterminer les courants de compensation en fonction du modèle en moments magnétiques du navire et des modèles en moments magnétiques des boucles de courant) de sorte que la solution est physiquement pertinente. Selon une caractéristique particulière, la deuxième méthode de résolution basée sur une régularisation parcimonieuse s’appuie sur au moins une des normes et . Comme déjà mentionné plus haut, l’avantage de ces normes et est de faciliter les calculs. Selon certains modes de réalisation, le deuxième problème inverse s’écrit sous une des formes mathématiques suivantes : avec un vecteur des courants de compensation à estimer, une matrice de coefficients dépendant des modèles magnétiques des boucles de courant et un courant de référence (appliqué à chacune des boucles de compensation lors de la mesure des effets de celles-ci). Cette liste de formes mathématiques n’est pas exhaustive. Dans une implémentation particulière, la deuxième méthode de résolution basée sur une régularisation parcimonieuse utilise un algorithme de type méthode de Monte-Carlo par chaînes de Markov. La présente invention n’est cependant pas limitée à ce type de méthode. Dans un autre mode de réalisation de l'invention, il est proposé un produit programme d'ordinateur qui comprend des instructions de code de programme pour la mise en œuvre du procédé précité (dans l'un quelconque de ses différents modes de réalisation), lorsque ledit programme est exécuté sur un ordinateur. Dans un autre mode de réalisation de l'invention, il est proposé un médium de stockage lisible par ordinateur et non transitoire, stockant un produit programme d'ordinateur comprenant un jeu d'instructions exécutables par un ordinateur pour mettre en œuvre le procédé précité (dans l'un quelconque de ses différents modes de réalisation). Dans un autre mode de réalisation de l'invention, il est proposé un dispositif de gestion d’une signature magnétique d’un navire, le dispositif comprenant une machine de calcul, dédiée ou reprogrammable, configurée pour effectuer une étape de modélisation magnétique du navire comprenant elle-même : obtention de mesures de champs magnétiques associées chacune à des informations comprenant une position géographique, un cap magnétique du navire et une distance par rapport au navire ; et détermination d’un modèle en moments magnétiques du navire, par résolution d’un premier problème inverse visant à déterminer des moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées ; la résolution du premier problème inverse étant réalisée avec une première méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. Dans des modes de réalisation du dispositif de gestion d’une signature magnétique d’un navire, la machine de calcul est également configurée pour la mise en œuvre d’autre(s) étape(s) (notamment l’étape de prédiction et/ou l’étape de compensation) et/ou caractéristique(s) du procédé précité, dans l'un quelconque de ses différents modes de réalisation. 4. LISTE DES FIGURES D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante, donnée à titre d'exemple illustratif et non limitatif, en relation avec les dessins annexés, dans lesquels : La , déjà décrite en relation avec l'art antérieur, présente un organigramme d'un procédé de gestion d’une signature magnétique d’un navire selon un mode de réalisation de l’art antérieur ; La présente un organigramme d'un procédé de gestion d’une signature magnétique d’un navire selon un mode de réalisation particulier de l'invention ; et La représente un exemple de structure d’un dispositif de gestion d’une signature magnétique d’un navire selon un mode de réalisation de l’invention, permettant la mise en œuvre des étapes du procédé de la . 5. DESCRIPTION DÉTAILLÉE Sur toutes les figures du présent document, les éléments et étapes identiques sont désignés par une même référence numérique. On présente désormais, en relation avec la , un procédé de gestion d’une signature magnétique d’un navire selon un mode de réalisation particulier de l'invention. Ce procédé est par exemple mis en œuvre par un dispositif 30, tel que décrit ci-dessous en relation avec la . Dans un mode de réalisation particulier, le dispositif 30 est par exemple intégré dans une station de mesure portable embarquée à bord du navire. Dans une étape 21, le dispositif 30 obtient des mesures de champs magnétiques associées chacune à des informations comprenant une position géographique, un cap magnétique du navire et une distance par rapport au navire. Comme déjà mentionné plus haut, l’information de cap magnétique du navire peut être obtenue soit directement (par exemple avec un magnétomètre), soit indirectement (par exemple par calcul à partir du cap compas ou du cap géographique). Dans une étape 22, le dispositif détermine un modèle en moments magnétiques du navire, par résolution d’un premier problème inverse visant à déterminer des moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées. Comme détaillé plus bas, la résolution du premier problème inverse est réalisée avec une première méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. Les étapes 21 et 22 font partie d’une étape 20 de modélisation magnétique du navire, à l’issue de laquelle on obtient un modèle en moments magnétiques du navire, à un cap magnétique. Comme indiqué par le cadre référencé 20, les étapes 21 et 22 sont itérées pour chacun d’une pluralité de N caps magnétiques (généralement on considère les quatre (N = 4) caps magnétiques particuliers suivants : cap magnétique nord (quand le cap magnétique est égal à 0°), cap magnétique sud (quand le cap magnétique est égal à 180°), cap magnétique est (quand le cap magnétique est égal à 90°) et cap magnétique ouest (quand le cap magnétique est égal à 270°)). On obtient ainsi un modèle en moments magnétiques du navire pour chacun de cette pluralité de caps magnétiques. Dans une étape 23, le dispositif effectue une prédiction de la signature magnétique du navire à une position géographique et/ou un cap magnétique du navire et/ou une distance par rapport au navire autre(s) que celle, celui ou ceux des mesures de champs magnétiques, en fonction du modèle en moments magnétiques du navire (domaine des moments magnétiques) déterminé à l’étape 22. Dans une étape 24, le dispositif effectue une compensation de la signature magnétique du navire à l’aide d’un système d’immunisation (système de degaussing) comprenant des boucles de courant installées à bord du navire. L’étape de compensation 24 comprend elle-même les étapes suivantes : détermination à l’étape 25 d’un modèle en moment(s) magnétique(s) de chacune des boucles de courant, par simple soustraction entre le modèle en moments magnétiques obtenu pour le navire avec une boucle activée et le modèle en moments magnétiques obtenu pour le navire sans boucle activée ; pour chacune des boucles de courant, détermination à l’étape 26 d’un courant de compensation en fonction du modèle en moments magnétiques du navire et des modèles en moments magnétiques des boucles de courant ; et application à l’étape 27, sur les boucles de courant, des courants de compensation déterminés. Comme détaillé plus bas, dans une implémentation particulière, le calcul des courants de compensation est réalisé par résolution d’un deuxième problème inverse, visant à déterminer les courants de compensation en fonction du modèle en moments magnétiques du navire et des modèles en moments magnétiques des boucles de courant, avec une deuxième méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. On présente maintenant plus en détail divers principes et caractéristiques permettant d’expliciter le procédé de gestion d’une signature magnétique d’un navire de la et les diverses étapes comprises dans celui-ci. Principe général (modélisation magnétique du navire, prédiction et compensation) Le navire est considéré comme un plan de dipôles/moments magnétiques. Chaque point de ce plan se décompose ainsi en trois axes et l’ensemble des moments magnétiques associés à certains points de ce plan sont notés comme suit : : vecteur des moments selon l’axe longitudinal du navire (nous considérerons un moment positif s’il est orienté vers l’avant) ; : vecteur des moments selon l’axe transversal du navire (nous considérerons un moment positif s’il est orienté vers tribord) ; et : vecteur des moments selon l’axe vertical du navire (nous considérerons un moment positif s’il est orienté vers le bas). La signature magnétique du navire dans le domaine statique est principalement due à deux contributions distinctes : l’aimantation permanente : liée aux caractéristiques magnétiques intrinsèques des matériaux ferromagnétiques constitutifs du navire et à leur histoire (rémanence). Cette contribution est considérée comme variant lentement sauf sous certaines conditions particulières. Pour la prédiction de signature et le réglage de degaussing, cette contribution est considérée comme constante ; et l’aimantation induite : liée à la présence de matériaux ferromagnétiques dans le champ terrestre. Cette contribution est donc variable selon la localisation du navire (champ terrestre local) et l’attitude du navire (cap magnétique, roulis, tangage). On définit ainsi : : vecteur des moments longitudinaux permanents : vecteur des moments longitudinaux induits : vecteur des moments transversaux permanents : vecteur des moments transversaux induits : vecteur des moments verticaux permanents : vecteur des moments verticaux induits. La modélisation magnétique du navire va donc consister à estimer les vecteurs de moments (longitudinaux, transversaux et verticaux) et à séparer les deux contributions (permanent et induit). Pour cela, dans un système classique, le navire passe sur une station de mesure sous-marine où des capteurs vectoriels tri-axiaux sont alignés et fixes (profondeur de mesure figée) et où le navire passe au-dessus de cette ligne de sorte que des capteurs sont à bâbord et d’autres à tribord. La trace de la signature magnétique sur chacun des capteurs est ainsi récupérée via les trois composantes du champ magnétique mesuré (Bx, By, Bz). Cette mesure est réalisée à chacun des quatre caps magnétiques particuliers précités (cap magnétique nord, cap magnétique sud, cap magnétique est et cap magnétique ouest) de sorte que le champ magnétique terrestre local se trouve, pour sa part horizontale, aligné avec un des axes horizontaux du navire (longitudinal ou transversal) pour pouvoir modéliser le navire à partir de ces composantes du champ magnétique : Cap magnétique nord : la part permanente est toujours présente et le champ terrestre va induire une aimantation dans l’axe vertical et l’axe longitudinal : ; Cap magnétique sud : la part permanente est toujours présente et le champ terrestre va induire une aimantation dans l’axe vertical et l’axe longitudinal : ; Cap magnétique ouest : la part permanente est toujours présente et le champ terrestre va induire une aimantation dans l’axe vertical et l’axe transversal : ; et Cap magnétique est : la part permanente est toujours présente et le champ terrestre va induire une aimantation dans l’axe vertical et l’axe transversal : Sans boucle de compensation du champ terrestre vertical, on ne peut pas éliminer . En obtenant les différents moments magnétiques du navire, on peut ainsi recalculer la signature magnétique du navire selon un autre champ terrestre (autre position géographique) et/ou à une orientation différente (autre cap magnétique). En revanche le changement de distance demande à avoir une modélisation proche de la réalité physique (ce que l’invention permet contrairement à l’état de l’art) pour avoir une validité en dehors de la zone de mesure initiale. Pour le système de degaussing, le principe repose sur des boucles de courant positionnées selon l’un des trois axes du navire : Boucle L : boucle créant un moment magnétique longitudinal ; Boucle A : boucle créant un moment magnétique transversal ; et Boucle M : boucle créant un moment magnétique vertical. Le nombre de boucles de chaque type (L, A et M) et les caractéristiques de celles-ci dépendent du niveau de réduction de signature magnétique (compensation) désiré. Pour réaliser le réglage du degaussing, on va donc, classiquement, repasser sur la station de mesure à un cap magnétique autant de fois que l’on a de boucles afin de mesurer l’effet magnétique de chacune des boucles à un courant connu Iref. On obtiendra ainsi dans la mesure l’effet magnétique du navire déjà mesuré seul et le champ magnétique créé par la boucle active. Schématiquement, un système de degaussing classique cherche donc à estimer les courants à partir de mesures de champs magnétiques acquises par une station de mesure dans les trois axes Bx, By, Bz pour les quatre caps magnétiques particuliers précités (cap magnétique nord, cap magnétique sud, cap magnétique est et cap magnétique ouest). Dans une implémentation particulière, la présente invention vise à régler les courants directement dans l’espace des moments magnétiques au lieu de raisonner en champ magnétique. Principe de fonctionnement d’un système de degaussing Comme indiqué ci-dessus, le système de degaussing est constitué de boucles L, A et M. Considérons par exemple quatre boucles de chaque type. Le réglage consiste donc à trouver les courants à injecter dans chacune des boucles en divisant ce courant en deux parties : : courant dans la boucle X qui est constant pour compenser une aimantation permanente ; et : courant dans la boucle X qui est variable pour compenser une aimantation induite. Le réglage est fait sur station de mesure à un champ terrestre local donné. Le système calculera ainsi, par règle de proportionnalité, le courant induit à injecter pour un champ terrestre et une attitude différente (projection du champ terrestre dans le repère navire). Le principe du degaussing est donc de trouver : , , , permettant de compenser ; , , , permettant de compenser ; , , , permettant de compenser ; , , , permettant de compenser ; , , , permettant de compenser ; et , , , permettant de compenser . Définition du problème pour la compensation (système de degaussing) On cherche le modèle magnétique du navire en moments magnétiques : et ; et ; et (pas séparable sauf si on a une boucle de compensation du champ terrestre vertical sur la station de mesure). Ce modèle est calculé à partir des mesures de champ magnétique (Bx, By, Bz) aux quatre caps magnétiques particuliers précités sur X capteurs. Pour le système de degaussing, la mesure de chaque effet de boucle est réalisé à un courant donné (ici qui est le même quelle que soit la boucle, dans la pratique, il peut être différent). On modélise chaque boucle X par un moment (ici un seul mais on pourra en trouver plusieurs, l’important étant qu’une boucle longitudinale (respectivement transversale ou verticale) compense un moment longitudinal (respectivement transversal ou vertical) du navire. On cherche les courants tels que : (1) (2) (3) (4) (5) Avec de dimension x*1 et l , t , v respectivement le nombre de moments dans le vecteur , et respectivement. En notant : On obtient sous forme matricielle : (6) Mesures de champs magnétiques du navire La mesure vaut : avec les vecteurs de coefficients kn dépendant de la position de la mesure. Ces coefficients ont donc une valeur différente pour chaque point de mesure. Mesures des effets des boucles La mesure va contenir les effets magnétiques du navire et de la boucle activée. Le champ magnétique étant régi par le principe de superposition, on pourra, connaissant la mesure de l’effet magnétique du navire, en déduire l’effet de la boucle seule. Les effets de boucles mesurés valent : (7) (8) (9) Relation entre champs magnétiques mesurés, effets de boucles et courants recherchés D’après les équations et en supposant un système parfait permettant de réaliser les mêmes points de mesure pour le navire et les effets de boucle (ou par modélisation et recalcul en un même point), on obtient : (10) En notant : (11) Ces équations montrent bien l’intérêt de la présente invention de pouvoir résoudre le problème directement dans le domaine des moments magnétiques (équation 6) plutôt que dans celui des champs magnétiques (équation 11) comme dans les méthodes classiques. D’une part on s’affranchit de la géométrie (termes kn) et d’autre part le problème à résoudre est de dimension plus faible. Propriétés des moments aux caps magnétiques Comme décrit plus haut, la mesure est réalisée à quatre caps magnétiques particuliers précités, de sorte que : Cap magnétique nord : avec : Cap magnétique sud : avec : Cap magnétique ouest : avec : Cap magnétique est : avec : D’où : (12) (13) (14) (15) (16) Champs magnétiques mesurés pour les passes aux quatre caps magnétiques particuliers précités (cap magnétique nord, cap magnétique sud, cap magnétique est et cap magnétique ouest) La mesure donnant le champ magnétique, on peut donc écrire que : (17) (18) (19) (20) On va ainsi, par modélisation, et avec les équations (12) à (16), trouver , , , et (étape 22 de la ). Effets de boucles mesurés Chaque boucle est mesurée à un cap magnétique quelconque avec un courant connu (ici ) puisque le moment magnétique, et donc le champ magnétique, créé par une boucle à courant constant ne dépend pas du cap magnétique. Là où les méthodes antérieures retiraient à la mesure le champ magnétique créé par le navire seul au même cap magnétique pour obtenir le champ magnétique créé par une boucle selon les équations suivantes : , et , la méthode proposée va elle directement déterminer le modèle en moments magnétiques pour chacune des mesures et déterminer les moments magnétiques créés par les boucles par différence entre modèles en moments magnétiques au même cap magnétique (on soustrait au modèle issu de la mesure comprenant le navire et une boucle active, le modèle issu de la mesure du navire seul au même cap magnétique). On va ainsi, par modélisation, trouver , , , , , , , , , , et (étape 25 de la ). Réglage des courants Nous cherchons à compenser l’effet magnétique du navire avec les boucles de degaussing. Le problème se décompose donc comme suit : nous cherchons la combinaison linéaire des boucles longitudinales qui compense le(s) moment(s) longitudinal(aux) induit(s) du navire : (21) nous cherchons la combinaison linéaire des boucles longitudinales qui compense le(s) moment(s) longitudinal(aux) permanent(s) du navire : (22) nous cherchons la combinaison linéaire des boucles transversales qui compense le(s) moment(s) transversal(aux) induit(s) du navire : (23) nous cherchons la combinaison linéaire des boucles transversales qui compense le(s) moment(s) transversal(aux) permanent(s) du navire : (24) nous cherchons la combinaison linéaire des boucles verticales qui compense le(s) moment(s) vertical(aux) du navire : (25) Le même type de calcul que pour les autres axes est possible si une boucle de compensation du champ terrestre est employée pour mesurer le navire. Sans cela, le réglage induit permanent vertical est réalisé par expérience. D’où le problème à résoudre : (26) (27) (28) (29) (30) Ces équations sont une forme simplifiée considérant les boucles comme parfaites et équivalentes à un moment unique à la distance de mesure. Ces équations sont complexifiables en fonction de la réalité des modèles obtenus (multi moments dans des directions différentes). Dans ce cas, il faut prendre en compte l’impact d’une boucle dans les autres axes (minime en principe, sauf géométrie particulière). Plus généralement (écriture sous forme générique), nous cherchons à compenser le navire avec p+q+r boucles de degaussing (p boucle(s) de type L, q boucle(s) de type A et r boucle(s) de type M). Le problème se décompose donc comme suit : nous cherchons la combinaison linéaire des boucles longitudinales qui compense le(s) moment(s) longitudinal(aux) induit(s) du navire : nous cherchons la combinaison linéaire des boucles longitudinales qui compense le(s) moment(s) longitudinal(aux) permanent(s) du navire : nous cherchons la combinaison linéaire des boucles transversales qui compense le(s) moment(s) transversal(aux) induit(s) du navire : nous cherchons la combinaison linéaire des boucles transversales qui compense le(s) moment(s) transversal(aux) permanent(s) du navire : nous cherchons la combinaison linéaire des boucles verticales qui compense le(s) moment(s) vertical(aux) du navire : Le même type de calcul que pour les autres axes est possible si une boucle de compensation du champ terrestre est employée pour mesurer le navire. Sans cela, le réglage induit permanent vertical est réalisé par expérience. D’où le problème à résoudre : Ces équations sont une forme simplifiée considérant les boucles comme parfaites et équivalentes à un moment unique à la distance de mesure. Ces équations sont complexifiables en fonction de la réalité des modèles obtenus (multi moments dans des directions différentes). Dans ce cas, il faut prendre en compte l’impact d’une boucle dans les autres axes (minime en principe). Le problème considéré à travers les équations (26’) à (30’) s’écrit sous la forme matricielle suivante : (31) avec le vecteur des moments magnétiques du navire, , la matrice des moments des boucles de dimension ( l + t + v )*( p + q + r ), et le vecteur des courants à estimer . Estimation des moments magnétiques du navire à partir des champs magnétiques mesurés (étape 22 de la ) L’estimation des moments magnétiques du navire (c'est-à-dire la détermination du modèle en moments magnétique du navire) à partir des champs magnétiques est un problème inverse bien connu en magnétisme du navire. L’angle pris dans cette invention est de faire l’hypothèse (qui est vérifiée en pratique) que les sources d’aimantation (ou moments dipolaires) qui expliquent les champs magnétiques observés en dehors de la zone de champ proche sont en nombre restreints (la nature des sources physiques d’aimantation sont connues des marins et ce principe est en cohérence avec le design du système de degaussing avec un nombre restreint de boucles et donc de moments magnétiques équivalents). Mathématiquement, ce type de problème inverse trouve une traduction pertinente dans la manière de régulariser le problème ; régulariser revient à imposer des contraintes sur les solutions attendues. Une régularisation du problème couramment utilisée dans l’état de l’art pour les problèmes inverses est celle de Tikhonov (aussi appelée en norme , notée et qui contient la somme des carrés de la variable régularisée), telle que : avec le vecteur de moments magnétiques du modèle en moments magnétiques du navire, le vecteur des mesures de champs magnétiques du navire, et la matrice (matrice de coefficients dépendant des positions associées aux mesures de champs magnétiques) telle que avec , x,y,z étant les coordonnées dans le repère navire et r=racine(x²+y²+z²). Celle-ci (et ses variantes) n’est toutefois pas pertinente au regard de la réalité physique car elle ne permet pas de prendre en compte que les sources d’aimantation sont localisées et en nombre restreint. Comme déjà mentionné plus haut, la présente invention repose sur l’utilisation d’une méthode de résolution basée sur une régularisation parcimonieuse (utilisant un a priori favorisant la parcimonie, aussi appelé « a priori parcimonieux » ou « sparse promoting priors » en anglais) pour résoudre le problème inverse visant à déterminer les moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées. L’usage d’une telle méthode de régularisation parcimonieuse dans l’étape 22 est justifié par la réalité physique : il est ainsi possible d’estimer automatiquement un nombre de sources limité expliquant les champs magnétiques mesurés. Dans des modes de réalisation de l’invention de l’étape 22, la méthode de résolution basée sur une régularisation parcimonieuse s’appuie sur au moins une des normes et . Ainsi, le problème inverse peut s’écrire sous plusieurs formes mathématiques, comme listé ci-après. Pour la norme : Pour la norme : Ou bien de manière équivalente : Il également possible de combiner les pénalités des normes et : Ou encore de combiner les pénalités des normes et : De nombreux algorithmes existent pour résoudre ces problèmes qu’on retrouve dans de nombreux domaines scientifiques. Un compromis en pratique doit être trouvé selon la puissance de calcul disponible, la facilité d’implantation, les performances attendues, l’estimation de l’incertitude sur le vecteur . Deux grandes familles d’algorithmes peuvent être convoquées pour ces problèmes : les algorithmes d’optimisation algébrique d’un côté et les algorithmes par simulation de l’autre. Dans un mode de réalisation particulier, pour l’estimation des moments magnétiques du navire, une méthode bayésienne et l’algorithme associé donnent un bon compromis entre vitesse de convergence et qualité d’estimation du vecteur , tout en fournissant les incertitudes pour les moments estimés. Estimation des courants à partir des moments magnétiques du navire (étape 26 de la ) Les moments magnétiques du navire (moments dipolaires) ayant été estimés (étape 22 de la ), il s’agit ensuite, dans l’étape 26 de la , d’estimer les courants de compensation à partir de ceux-ci. La même hypothèse de parcimonie de la solution du jeu d’équations (26) à (30) - ou plus généralement (26’) à (30’) - est adoptée et un a priori favorisant la parcimonie est choisi pour trouver cette solution. Autrement dit, certains courants peuvent être nuls dans telle ou telle situation en fonction de moments magnétiques estimés à l’étape 22. Dans des modes de réalisation de l’invention de l’étape 22, la méthode de résolution basée sur une régularisation parcimonieuse s’appuie sur au moins une des normes et . Ainsi, le problème inverse peut s’écrire sous plusieurs formes mathématiques, comme listé ci-après. Pour la norme : avec le vecteur des courants de compensation à estimer, la matrice des moments des boucles (c'est-à-dire une matrice de coefficients dépendant des modèles en moment magnétiques des boucles de courant) et un courant de référence (déjà discuté plus haut). Un exemple du vecteur est mentionné plus haut dans l’équation (6). Pour la norme : Ou bien de manière équivalente : Il également possible de combiner les pénalités des normes et : Ou encore de combiner les pénalités des normes et : Pour cette étape 26 du procédé, la méthode de résolution basée sur une régularisation parcimonieuse utilise par exemple un algorithme de type méthode de Monte-Carlo par chaînes de Markov. En d’autres termes, l’estimation du vecteur est par exemple assurée par un algorithme de type MCMC. Mais toute autre approche est envisageable. On présente désormais, en relation avec la un exemple de structure d’un dispositif 30 de gestion d’une signature magnétique d’un navire selon un mode de réalisation de l’invention, permettant la mise en œuvre des étapes du procédé de la . Le dispositif 30 comprend une mémoire vive 33 (par exemple une mémoire RAM), une unité de traitement 31 équipée par exemple d’un processeur, et pilotée par un programme d’ordinateur 320 stocké dans une mémoire morte 32 (par exemple une mémoire ROM ou un disque dur). A l’initialisation, les instructions de code du programme d’ordinateur sont par exemple chargées dans la mémoire vive 33 avant d’être exécutées par le processeur de l’unité de traitement 31. Cette illustre seulement une manière particulière, parmi plusieurs possibles, de réaliser le dispositif 30 afin qu’il effectue certaines étapes du procédé de gestion d’une signature magnétique d’un navire selon l’invention (selon l’un quelconque des modes de réalisation et/ou variantes décrit(e)s ci-dessus en relation avec la ). En effet, ces étapes peuvent être réalisées indifféremment sur une machine de calcul reprogrammable (un ordinateur PC, un processeur DSP ou un microcontrôleur) exécutant un programme comprenant une séquence d’instructions, ou sur une machine de calcul dédiée (par exemple un ensemble de portes logiques comme un FPGA ou un ASIC, ou tout autre module matériel). Dans le cas où le dispositif 30 est réalisé avec une machine de calcul reprogrammable, le programme correspondant (c’est-à-dire la séquence d’instructions) pourra être stocké dans un médium de stockage amovible (tel que par exemple un CD-ROM, un DVD-ROM, une clé USB) ou non, ce médium de stockage étant lisible partiellement ou totalement par un ordinateur ou un processeur. Procédé de gestion d’une signature magnétique d’un navire, exécuté par une machine de calcul et comprenant une étape (20) de modélisation magnétique du navire comprenant elle-même : obtention (21) de mesures de champs magnétiques associées chacune à des informations comprenant une position géographique, un cap magnétique du navire et une distance par rapport au navire ; et détermination (22) d’un modèle en moments magnétiques du navire, par résolution d’un premier problème inverse visant à déterminer des moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées ; caractérisé en ce que la résolution du premier problème inverse est réalisée avec une première méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. Procédé selon la revendication 1, dans lequel la première méthode de résolution basée sur une régularisation parcimonieuse s’appuie sur au moins une des normes et . Procédé selon la revendication 2, dans lequel le premier problème inverse s’écrit sous une des formes mathématiques suivantes : avec un vecteur des moments magnétiques du navire à estimer, un vecteur des mesures de champs magnétiques et une matrice de coefficients dépendant des positions associées aux mesures de champs magnétiques. Procédé selon l’une quelconque des revendications 1 à 3, dans lequel la première méthode de résolution basée sur une régularisation parcimonieuse est une méthode bayésienne. Procédé selon l’une quelconque des revendications 1 à 4, comprenant une étape (23) de prédiction de la signature magnétique du navire à une position géographique et/ou un cap magnétique du navire et/ou une distance par rapport au navire autre(s) que celle, celui ou ceux des mesures de champs magnétiques, en fonction du modèle en moments magnétiques du navire. Procédé selon l’une quelconque des revendications 1 à 4, comprenant une étape (24) de compensation de la signature magnétique du navire à l’aide d’un système d’immunisation comprenant des boucles de courant installées à bord du navire, l’étape de compensation de la signature magnétique du navire comprenant : détermination (25) d’un modèle en moment(s) magnétique(s) de chacune des boucles de courant ; pour chacune des boucles de courant, détermination (26) d’un courant de compensation en fonction du modèle en moments magnétiques du navire et des modèles en moments magnétiques des boucles de courant ; et application (27), sur les boucles de courant, des courants de compensation déterminés. Procédé selon la revendication 6, dans lequel le calcul des courants de compensation est réalisé par résolution d’un deuxième problème inverse, visant à déterminer les courants de compensation en fonction du modèle en moments magnétiques du navire et des modèles en moments magnétiques des boucles de courant, avec une deuxième méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. Procédé selon la revendication 7, dans lequel la deuxième méthode de résolution basée sur une régularisation parcimonieuse s’appuie sur au moins une des normes et . Procédé selon la revendication 8, dans lequel le deuxième problème inverse s’écrit sous une des formes mathématiques suivantes : avec un vecteur des courants de compensation à estimer, une matrice de coefficients dépendant des modèles en moments magnétiques des boucles de courant et un courant de référence. Procédé selon l’une quelconque des revendications 7 à 9, dans lequel la deuxième méthode de résolution basée sur une régularisation parcimonieuse utilise un algorithme de type méthode de Monte-Carlo par chaînes de Markov. Produit programme d'ordinateur (320), comprenant des instructions de code de programme pour la mise en œuvre du procédé selon au moins une des revendications 1 à 10, lorsque ledit programme est exécuté sur un ordinateur. Médium de stockage (32) lisible par ordinateur et non transitoire, stockant un produit programme d'ordinateur (320) selon la revendication 11. Dispositif (30) de gestion d’une signature magnétique d’un navire, le dispositif comprenant une machine de calcul (31, 32, 33), dédiée ou reprogrammable, configurée pour effectuer une étape de modélisation magnétique du navire comprenant elle-même : obtention de mesures de champs magnétiques associées chacune à des informations comprenant une position géographique, un cap magnétique du navire et une distance par rapport au navire ; et détermination d’un modèle en moments magnétiques du navire, par résolution d’un premier problème inverse visant à déterminer des moments magnétiques du navire à partir des mesures de champs magnétiques et des informations associées ; caractérisé en ce que la résolution du premier problème inverse est réalisée avec une première méthode de résolution basée sur une régularisation parcimonieuse utilisant un a priori favorisant la parcimonie. Dispositif selon la revendication 13, dans lequel la machine de calcul est également configurée pour effectuer une étape de prédiction de la signature magnétique du navire à une position géographique et/ou un cap magnétique du navire et/ou une distance par rapport au navire autre(s) que celle, celui ou ceux des mesures de champs magnétiques, en fonction du modèle en moments magnétiques du navire. Dispositif selon la revendication 13, dans lequel la machine de calcul est également configurée pour effectuer une étape de compensation de la signature magnétique du navire à l’aide d’un système d’immunisation comprenant des boucles de courant installées à bord du navire, l’étape de compensation de la signature magnétique du navire comprenant : détermination d’un modèle en moment(s) magnétique(s) de chacune des boucles de courant ; pour chacune des boucles de courant, détermination d’un courant de compensation en fonction du modèle en moments magnétiques du navire et des modèles en moments magnétiques des boucles de courant ; et application, sur les boucles de courant, des courants de compensation déterminés.