Installation tracking by measuring a perturbed oscillating magnetic field

Abstract

This method for tracking the pose of a receiver (24) relative to a transmitter (22) comprises obtaining a raw magnetic measurement, estimating the contribution of a disturbing element (90) to the raw magnetic measurement, calculating a refined magnetic measurement by subtracting the contribution of the disturbing element (90) from the raw magnetic measurement, and deducing a refined pose of the receiver (24) relative to the transmitter (22) from the refined magnetic measurement. The contribution of the disturbing element (90) is estimated from secondary current densities of eddy currents induced by the transmitter (22) in the disturbing element (90), said secondary current densities resulting from an interpolation of primary current density components of said eddy currents along different directions and at different points on a surface (94) of the disturbing element (90). Figure for the abstract: Fig. 1

Description

Title of the invention: Pose tracking by measuring a perturbed oscillating magnetic field. Field of the invention

[0001] The present invention relates to tracking the pose, that is, the position and orientation, of a first object relative to a second object, by measuring, at the level of the first object, an oscillating magnetic field produced by a transmitter attached to the second object. More particularly, the present invention relates to implementing this tracking when this oscillating magnetic field is disturbed by the presence, in the measurement environment, of a conductive element excited by the transmitter such that it emits a disturbing magnetic field that distorts the lines of the oscillating magnetic field. Technological background

[0002] To track the pose of an object, that is, its position and orientation, in a frame of reference, it is known to use magnetic technology pose tracking systems, commonly called EMTs (from the English "Electromagnetic tracker"). These systems generally comprise a transmitter fixed in the frame of reference, consisting of three coils arranged to form a fixed, substantially orthogonal trihedron, and a receiver attached to the object whose pose is to be determined. This receiver consists of three magnetic sensors arranged to form a movable, substantially orthogonal trihedron. Time-dependent electric currents flow in the coils of the transmitter and generate three magnetic fields, which are detected by the sensors of the receiver. Each sensor of the receiver measures, for each of the emitted magnetic fields, the projection of that field onto the direction in which the sensor is directed.This gives a total of nine components that allow us to move from a moving trihedron to a fixed trihedron. Indeed, these nine components depend on the position and orientation of the receiver relative to the transmitter. These nine components form what we will subsequently call a "magnetic matrix" in which the components of the same field are grouped within the same column of the matrix and the components measured by the same sensor are grouped within the same row of the matrix.

[0003] To enable the identification, on the receiver side, of the contribution of each coil to the captured magnetic field, it is known to excite the transmitter coils by means of alternating currents at their natural frequencies. The generated magnetic fields are thus alternating magnetic fields that can be distinguished from one another by demodulation using knowledge of the excitation frequency of the coils. Sensitivity to constant magnetic fields can thus be neglected and the sensors can be formed from simple coils across whose terminals the induced voltage is measured.

[0004] A drawback of magnetic technology pose tracking systems is their sensitivity to the presence of conductive objects in the environment. Consider a conductive object subjected to a single-frequency incident magnetic field (generated by a single coil) that varies over time and is not uniform throughout the object's volume. Eddy currents at the frequency of the incident field flow through this conductive object. These currents generate a reflected field that is superimposed on the incident field, distorting the lines of the total field at that frequency. This phenomenon disrupts the pose measurement by the magnetic technology.The magnetic matrix Mt at the actual pose (r, R) can therefore be decomposed into Mt = Me + Mc where Me is the matrix measured in the absence of a disturbance (which can easily be predicted at a given pose thanks to models of the source and sensor defects) and Mc is the contribution of the conducting object to the magnetic matrix Mt after demodulation (i.e. after projection of the measurement of the total magnetic field onto a "demodulation direction").

[0005] To circumvent this difficulty, different types of solutions have been proposed.

[0006] Solutions of a first type rely on exploiting the variation in the response of disturbances to incident fields of different frequencies. For example, methods are known from US6172499 and US6762600 in which the contribution of a disturbance is identified and removed by taking advantage of the variation in the phases of the disturbance's responses at different frequencies. A method is also known from US11187823 for obtaining a reference signal or pose by occasionally using very low frequencies (for example, around 100 Hz) at which conductive objects have a negligible impact. This reference signal / pose is then used to identify a correction for the fields at the main frequencies (for example, around 30 kHz), and this correction can be extrapolated to correct the fields at the main frequencies during subsequent poses.

[0007] Solutions of this first type have the advantage of being generic (they do not require dedicated work to study each disturbance whose impact must be compensated). However, they present several drawbacks. First, frequency changes are restrictive. Second, they require the sensor to be almost static when low frequencies are used.

[0008] Solutions of a second type consist of methods that perform corrections without a detailed physical model and employ a large number of correction parameters calibrated through the acquisition of a calibration trajectory (long and delicate) that meshes the target space. These parameters are typically either polynomial coefficients (or spherical harmonic development moments), or lookup tables (well known by the acronym LUT, from the English “Look-up table”) that associate a correction with an approximate exposure or field measurement, or a combination of both (a table that provides coefficients for a local correction calculation). Such an approach is known in particular from FR2664044, which describes the correction of disturbances caused by an extended conducting object (an aircraft cockpit) fixed relative to the source. A polynomial correction of the field measurement is proposed there. The correction parameters are constructed from a calibration trajectory acquired in the final operating environment.Several improvements have been proposed since, notably by US6400139, which suggests using "witness sensors" fixed in the source reference frame to make the correction model more adaptable by eliminating the need for impact mapping steps and robust to its potential movements. A similar approach is also described by H. Himberg in "Latency and Distortion Compensation in Augmented Environment Using Electromagnetic Trackers," PhD thesis, Virginia Commonwealth University, 2010.

[0009] However, solutions of this second type are only suitable for correcting disturbances in a very small and fairly stable space. Furthermore, the calibration phases are cumbersome: they require significant acquisition times and most often necessitate a reference system capable of providing the exact sensor pose at which the measurements are acquired by the EMT. Finally, such methods are not sustainable when the disturbance is significant and moves with the sensor in the source frame (indeed, whereas, when the disturbance is fixed in the source frame, the disturbance depends only on the sensor's position in a subspace V of p3 defining the desired range of the EMT, when the disturbance moves with the sensor, the disturbance changes at each pose, which increases the space to be traversed: it is then necessary to map all SO(3)xV, SO(3) being the set of rotation matrices in dimension 3).

[0010] A third type of solution consists of methods that attempt a detailed model of conducting objects to reconstruct their impact on the fields. In the case where the perturbator is simple (geometric properties that trivialize the problem), exact theoretical calculations of the impact of the perturbators can be used once their pose is known or estimated. This is what is proposed for infinite planes in US2023083524. In other cases, a detailed model of the perturbator with methods approaching finite elements can be carried out, the main limitation of this type of approach being the computation time, incompatible with use for real-time correction. (See H. Gietler, Object localization for) Autonomous Systems based on electromagnetic fields. PhD thesis, University of Klagenfurt, 2021 presents a method in which computation times are presented as acceptable thanks to a discretization of the perturbing volume and a resolution by integral calculus with simplifying assumptions on the elementary volumes, this method is not functional in practice (indeed the calculation of the necessary field integrals is not straightforward and is not even explained in the document) and is only effective in the case of a fixed perturbing element in the source frame (its adaptation to a fixed perturbing element in the sensor frame risks encountering an execution speed incompatible with a real-time application).

[0011] Thus, solutions of this third type are either limited to correcting disturbances produced by perturbators having a very simple geometric shape, or incompatible with a real-time application.

[0012] It should be noted that solutions of this third type can be combined with solutions of the second type, as described for example in FR2722299. Indeed, this document describes a process using a fine modeling of small perturbators (metallic parts of a pilot's helmet) fixed with respect to the sensor, combined with a calibration of multipolar moments of a spherical harmonic development centered on other perturbators.

[0013] Finally, a fourth type of solution combines EMT measurements with inertial measurements. Such solutions are known, in particular, from US20160356601, US 10746819, and K. O'Donoghue, A. Jaeger, Herman, and C. Murphy, Padraig. Sensor fusion hardware platform for robust electromagnetic navigation. IEEE International Instrumentation and Measurement Technology Conference (I2MTC), 2022. However, while these solutions offer high robustness to the occurrence of disturbances or noise, they do not improve the static accuracy of the system, which is defined by EMT alone.

[0014] Thus, there is a need for a solution to correct the disturbance produced by a conductive object of relatively complex shape present in the environment of an EMT and to improve the static accuracy of this EMT, when this conductive object is mobile with the sensor or when the desired range of the EMT is extended, said solution being suitable for a real-time application and not requiring quasi-static phases of the sensor. Description of the invention

[0015] One objective of the invention is to correct the disturbance produced by a relatively complex-shaped disturbing element present in the environment of an EMT, when this disturbing element is moving with the sensor and / or when the specified high-angle zone of the desired EMT is extended. Other objectives are to to allow this correction in real time, to improve the static accuracy of the EMT, and to be able to do without quasi-static phases of the EMT sensor.

[0016] To this end, the invention relates, according to a first aspect, to a method for tracking the positioning of a receiver relative to a transmitter, the transmitter comprising at least two magnetic generators, each oriented along its own direction and capable of generating an alternating magnetic field having a distinctive characteristic enabling the alternating magnetic field generated to be distinguished from the alternating magnetic field generated by each other magnetic generator, the magnetic generators being arranged relative to each other so that their directions are non-coplanar, the method being implemented by a data processing unit and comprising the following step: - determination of a position of the receiver relative to the transmitter at a determination instant, said determination comprising the following sub-steps: • a) obtaining a raw magnetic measurement composed of measurements, carried out at the time of determination, of local components, along at least two non-coplanar measurement directions attached to the receiver, of at least two resulting alternating magnetic fields constituting an ambient magnetic field, each of said resulting alternating magnetic fields having a distinctive characteristic identical to that of one of the generated alternating magnetic fields, • b) determination of an imprecise pose of the receiver relative to the emitter at the time of determination, • c) estimation, based on an approximate pose that is a function of the imprecise pose, of the contribution of a disturbing element to the raw magnetic measurement, • d) calculation of a refined magnetic measurement by subtracting the contribution of the disturbing element from the raw magnetic measurement, and • e) deduction of a refined pose of the receiver relative to the emitter at the time of determination from the refined magnetic measurement,

[0017] wherein the disturbing element is modeled as a surface mesh composed of a set of mesh points regularly distributed over a surface of the disturbing element and connected to each other by straight segments delimiting between them two-dimensional and parallelogrammatic mesh cells, and substep c) of estimating the contribution of the disturbing element comprises, for each of the generated alternating magnetic fields,

[0018] the estimation of local components, at the receiver level, for a receiver pose equal to the approximate pose, of a disturbing magnetic field produced by eddy currents induced in the disturbing element by said generated alternating magnetic field, by applying the Biot-Savart law to secondary current densities of these eddy currents, said secondary current densities resulting from an interpolation of primary current density components of said eddy currents, said primary current density components being a function of a pose of the disturbing element relative to the emitter and being composed of: - first components, along a first direction of the surface mesh, of first primary current densities in first internal segments of the surface mesh oriented along a second direction of the surface mesh, in particular in the middle of said first internal segments, and - of second components, along the second direction of the surface mesh, of second primary current densities in second internal segments of the surface mesh oriented along the first direction of the surface mesh, in particular in the middle of said second internal segments. According to particular embodiments of the invention, the tracking method also has one or more of the following characteristics, taken individually or in any technically possible combination(s): - the measurement directions are greater than or equal to three; - The primary current density components are obtained using the following formula or an equivalent formula:

[0019]

[0020] where jx is a vertical matrix of the set of the first components of the Jj first primary current densities, is a vertical matrix of the set Jj of the second components of the second secondary current densities, is a vertical matrix of normal components of the generated alternating magnetic field, orthogonal to the surface of the disturbing element, at each of the mesh points not belonging to a contour of the surface mesh, and jç1 is the inverse matrix of a matrix constructed from the local Maxwell equations discretized between the primary current density components and normal components of the generated alternating magnetic field; - the disturbing element is fixed relative to the receiver and substep c) of estimating the contribution of the disturbing element includes determining a pose of the disturbing element relative to the emitter, estimating the components of the primary current densities as a function of said pose of the disturbing element relative to the emitter, and determining the secondary current densities by interpolating the components of the primary current densities; - the matrix is ​​determined prior to the implementation of the step of determining the pose; the disturbing element is fixed relative to the emitter and the secondary current densities are determined prior to the implementation of the installation determination step; the pose determination step includes several successive iterations of substeps c), d) and e), the approximate pose being constituted, for the first iteration, by the imprecise pose and, for each subsequent iteration, by the refined pose deduced at the end of the previous iteration; the tracking process includes several iterations of the pose determination step to determine receiver poses at different determination times, the first iteration constituting a step of determining an initial receiver pose and subsequent iterations constituting steps of determining an updated pose, the imprecise pose determined during the initial pose determination step being deduced from the raw magnetic measurement; the imprecise pose determined during at least one of the steps of determining an updated pose is extrapolated from at least one previous pose determined during a previous iteration of the pose determination step; secondary current densities are determined at evaluation points regularly distributed within each mesh cell; The local components of the disturbing magnetic field are estimated using the following formula or an equivalent formula:

[0022] where bj(F) is a vector of said local components of the disturbing magnetic field, Pq is the magnetic permeability of free space, is a 3-dimensional vector of the components, in a frame attached to the disturbing element, of one of the secondary current densities, Si is the surface area of ​​an elementary cell associated with said estimated surface current density, is an average thickness of a portion of the surface of the disturbing element, within which the eddy currents are concentrated, at the level of the mesh cell comprising the elementary cell, is a vector of the coordinates of a center of said elementary cell in a frame attached to the disturbing element, r is a vector of the coordinates of the approximate pose in said frame attached to the disturbing element, X is the cross product operator and | • | is the 2-norm,

[0023] the elementary cell being constituted by one of the mesh cells or by a subdivision of a subdivided mesh cell, each subdivided mesh cell being made up of elementary cells, each associated with one of the secondary current densities, regularly distributed within said subdivided mesh cell and covering its entirety; and - the mesh cells are rectangular.

[0024] The invention also relates, according to a second aspect, to a method for calibrating an electromagnetic tracking system comprising a step of tracking a pose of a receiver of the electromagnetic tracking system relative to a transmitter of the electromagnetic tracking system by means of a tracking method according to the first aspect, on the basis of hypothetical values ​​of calibration parameters of the electromagnetic tracking system, and updating the hypothetical values ​​of the calibration parameters so as to minimize a sum of error criteria each representative of an error on a pose of the receiver relative to the transmitter determined during the tracking step.

[0025] The invention also relates, according to a third aspect, to a computer program product comprising code instructions for the implementation, by a processor, of a tracking method according to the first aspect or of a calibration method according to the second aspect.

[0026] The invention also relates, according to a fourth aspect, to a computer-readable recording medium on which a computer program product is stored according to the third aspect.

[0027] Finally, according to a fifth aspect, the invention relates to an assembly comprising a reference object, a moving object, and an electromagnetic tracking system for tracking the position of the moving object relative to the reference object, said electromagnetic tracking system comprising: - a transmitter attached to the reference object, said transmitter comprising at least two magnetic generators, each oriented along its own direction and capable of generating an alternating magnetic field having a distinctive characteristic enabling said alternating magnetic field to be distinguished from the alternating magnetic field generated by each other magnetic generator, the magnetic generators being arranged relative to each other so that their directions are non-coplanar, and - a receiving device comprising a receiver attached to the moving object, said receiving device being capable of measuring local components, along at least two non-coplanar directions, of at least two resulting alternating magnetic fields constituting an ambient magnetic field, each of said resulting alternating magnetic fields having a distinctive characteristic identical to that of one of the generated alternating magnetic fields, so as to obtain a raw magnetic measurement,

[0028] the assembly also comprising a disturbing element capable of impacting the measurement of the local components of the resulting alternating magnetic fields by the receiving device when the disturbing element is subjected to the alternating magnetic fields generated by the transmitter,

[0029] in which the electromagnetic tracking system includes a data processing unit configured for the implementation of a tracking method according to the first aspect or a calibration method according to the second aspect. Brief description of the Figures

[0030] Other features and advantages of the invention will become apparent from the following description, given solely by way of example and with reference to the accompanying drawings, in which: - [Fig. 1] is a schematic view of an assembly consisting of a first object and a second object comprising a tracking system according to an example of an embodiment of the invention, - [Fig.2] is a schematic view of an example of a surface mesh of a portion of the surface of a perturbing element of the whole of [Fig.1], - [Fig.3] is a diagram of an example of a determination process that can be implemented by the monitoring systems of [Fig.1], - [Fig.4] is a diagram illustrating a first variant of a step for providing calculation parameters for the process of [Fig.3], - [Fig.5] is a diagram illustrating a second variant of the step of providing calculation parameters from [Fig.4], - [Fig.6] is a diagram illustrating a step in determining an initial pose of the process of [Fig.3], - [Fig.7] is a diagram illustrating a first variant of a substep of the initial pose determination step of [Fig.6] for estimating the contribution of a disturbing element to a raw magnetic measurement produced by a receiver of the tracking system, - [Fig.8] is a diagram illustrating a second variant of the sub-step of [Fig.7], - [Fig. 8] is a diagram illustrating a step in determining an updated pose of the process of [Fig. 3], and - [Fig. 10] is a diagram illustrating a calibration process that can be implemented by the monitoring systems of [Fig. 1]. Detailed description of an example of implementation

[0031] The sets 10, 11 shown in Figures 1 and 2 each comprise a first object 12 and a second object 14, the second object 14 being mobile relative to the first object 12.

[0032] The first object 12 consists of a reference object whose pose in a reference frame, for example an inertial frame, is known. Here and subsequently, the pose is defined as consisting of all the position and orientation data of an object in space. In three-dimensional space, the pose therefore comprises six components, consisting of: - three position components, and - three orientation components.

[0033] The first object 12 is, for example, a fixed structure in the reference frame. Alternatively (not shown), the first object is itself mobile in the reference frame but is equipped with sensors enabling the movement of the first object to be tracked in the reference frame.

[0034] The second object 14 consists of a mobile object in the reference frame and whose movement relative to this reference frame we wish to follow.

[0035] The first object 12 and the second object 14 are typically arranged relative to each other so that there is no permanent line of sight between the two objects 12, 14. In other words, the second object 14 is likely to be hidden from the view of the first object 12 during its movement.

[0036] The speed of movement of the second object 14 relative to the first object 12 is typically less than 30 m / s and 10 rad / s.

[0037] For example, the second object 14 consists of a pen and the first object 12 consists of a housing that is stationary relative to the writing surface: it is thus possible to record the writing. Alternatively, the second object 14 consists of a helmet and the first object 12 of a structure fixed in the helmet's movement environment.

[0038] A primary frame (Ri) is attached to the first object 12. This primary frame (Ri) is a direct orthonormal frame with origin Oi formed of a triplet of axes, represented in [Fig. 1], comprising: - a first primary axis ef, - a second primary axis orthogonal to the first axis ef, and - a third primary axis orthogonal to the first and second axes And

[0039] A secondary frame (R2) is attached to the second object 12. This secondary frame (R2) is a direct orthonormal frame with origin O2 formed by a triplet of axes, shown in [Fig. 1], comprising: - a first secondary axis e^, - a second secondary axis orthogonal to the first axis ct - a third secondary axis orthogonal to the first and second axes ^2 and e£.

[0040] The position of the second object 14 relative to the first object 12 can be characterized by the set formed by the rotation matrix R[j2 from the primary frame (Ri) to the secondary frame (R2) and the vector aUant from the origin O1 of the primary frame (Ri) to the origin O2 of the secondary frame (R2). Here and subsequently, the notation convention for rotation matrices from a frame A to a frame B is to use matrices whose columns are the unit vectors of frame B expressed in frame A. Thus, in this specific case, the rotation matrix R from the primary frame (Ri) to the secondary frame (R2) is constituted by the matrix whose columns are the coordinates of the axes 62, ^2 of the secondary frame (R2) expressed in the system of coordinates of the primary frame (Ri).

[0041] The assembly 10 also includes a tracking system 20 for tracking the placement of the second object 14 relative to the first object 12. This tracking system 20 consists of an electromagnetic tracking system, or EMT. It includes a transmitting device 21 comprising a transmitter 22 attached to the first object 12 and a receiving device 23 comprising a receiver 24 attached to the second object 14.

[0042] Said transmitter 22 and receiver 24 are typically configured to be functional when separated by a distance of less than 10 m. This helps to limit the size of the transmitter 22 and the receiver 24.

[0043] The emitter 22 comprises at least two, here three, magnetic generators 30, 32, 34. Each magnetic generator 30, 32, 34 behaves essentially as a magnetic dipole having an alternating magnetic dipole moment of amplitude m and frequency f. It is oriented along a proper direction, respectively gb g2, g3, defined here and subsequently as the direction of said dipole moment.

[0044] The amplitude m of the dipole moment is preferably substantially equal for all generators 30, 32, 34. The frequency of the dipole moment, however, is specific to each generator 30, 32, 34, so that each generator 30, 32, 34 is thus capable of generating an alternating magnetic field having a distinctive characteristic that allows said alternating magnetic field to be distinguished from the alternating magnetic field generated by each other magnetic generator 30, 32, 34, said distinctive characteristic being here constituted by the frequency of the magnetic field. Alternatively, the frequency of the dipole moment is substantially equal for all generators 30, 32, 34, the generators 30, 32, 34 then having another distinctive characteristic that allows the production of magnetic fields that can be distinguished from one another. The frequency of the dipole moment is typically between 500 Hz and 500 kHz, preferably between 5 kHz and 20 kHz.

[0045] For this purpose, each magnetic generator 30, 32, 34 comprises a coil, respectively 36, 38, 40, connected to a coil excitation current generator, respectively 42, 44, 46. Each coil 36, 38, 40 is oriented along a proper direction, defined by the axis around which the coil extends, this direction constituting the orientation direction of the magnetic generator 30, 32, 34 to which it belongs.

[0046] Each current generator 42, 44, 46 belongs to the transmitting device 21. Each current generator 42, 44, 46 is, for example, as shown, integrated into the first object 12. Alternatively (not shown), at least part of the current generators 42, 44, 46 is offset relative to the first object 12.

[0047] Each current generator 42, 44, 46 is configured to generate the excitation current of the corresponding coil 36, 38, 40 at a natural frequency, this frequency constituting the frequency of the alternating magnetic field generated by the magnetic generator 30, 32, 34 to which the coil 36, 38, 40 connected to this current generator 42, 44, 46 belongs.

[0048] Each current generator 42, 44, 46 is in particular configured to generate the excitation current of the corresponding coil 36, 38, 40 at a frequency between 500 Hz and 500 kHz, preferably between 5 kHz and 20 kHz.

[0049] The magnetic generators 30, 32, 34 are arranged relative to each other such that each is oriented along a proper direction gb g2, g3, said directions gi, g2, g3 being non-coplanar. In particular, the magnetic generators 30, 32, 34 are arranged such that their directions gi, g2, g3 are such that: - a first generator direction gi is coincident with the first primary axis ej, - a second direction of generator g2 is included in a plane defined by the first and second primary axes and forms a first angle 0 between 0 and j radians (inclusive) with the second primary axis and - a third generator direction g3 is included in a plane passing through the third primary axis and forming a second angle q> between 0 and - 0 radians (exclusive) with the plane defined by the first and third primary axes ef, said third generator direction g3 forming a third angle rp between 0 and j radians (inclusive) with the third primary axis e|.

[0050] For example, the magnetic generators 30, 32, 34 are arranged so that their directions gb g2, g3 are substantially orthogonal to each other, each direction gb g2, g3 being typically, as shown, substantially collinear with one of the axes ef, ej, ef of the primary frame (RJ, i.e. that the first and second angles 0, q> are less than or equal to 0.05 radian.

[0051] Preferably, the magnetic generators 30, 32, 34 have substantially coincident centers, that is, they can be modeled as magnetic dipoles located substantially at the same point. To this end, the coils 36, 38, 40 have substantially coincident centers, that is, the centers of the coils are separated in pairs by a distance less than the average radius of the coils 36, 38, 40. The center of the magnetic generators 30, 32, 34 constitutes the origin Oi of the primary frame (Ri).

[0052] The transmitting device 21 further includes a clock 48 providing a reference signal representative of time. This reference signal is used by the magnetic generators 30, 32, 34 to vary over time the alternating magnetic fields that they generate.

[0053] The receiving device 23 is capable of measuring the local components, at the level of the receiver 24, along the three axes ©2, 62 of the secondary frame (R2), of three resulting alternating magnetic fields constituting an ambient magnetic field, each of said resulting alternating magnetic fields having a distinctive characteristic, allowing said alternating magnetic field to be distinguished resulting from the other resulting alternating magnetic fields, identical to that of one of the generated alternating magnetic fields

[0054] To this end, the receiver 24 comprises three magnetic sensors 50, 52, 54, each capable of measuring a component of the ambient magnetic field along a measurement direction cb c2, c3 specific to the sensor 50, 52, 54, these measurement directions cb c2, c3 being non-coplanar. These sensors 50, 52, 54 are preferably arranged so that their measurement directions cb c2, c3 are substantially orthogonal to each other, each measurement direction cb c2, c3 being typically, as shown, substantially collinear with one of the axes Cp e^, ef of the primary frame (Ri).

[0055] Each magnetic sensor 50, 52, 54 is here formed of a coil, respectively 56, 58, 60. Each coil 56, 58, 60 is oriented along a proper direction, defined by the axis around which the coil extends, this direction constituting the measurement direction of the magnetic sensor 50, 52, 54 to which it belongs.

[0056] Preferably, the magnetic sensors 50, 52, 54 have substantially coincident centers. To this end, the coils 56, 58, 60 have substantially coincident centers, that is to say, the centers of the coils are separated in pairs by a distance less than the average radius of the coils 56, 58, 60. The center of the magnetic sensors 50, 52, 54 constitutes the origin O2 of the secondary frame (R2).

[0057] The receiving device 23 also includes tensiometers, 62, 64, 66 respectively, each connected to one of the coils, 56, 58, 60 respectively, to measure the voltage across said coil. Each tensiometer 62, 64, 66 is, for example, as shown, integrated into the first object 12. Alternatively (not shown), at least part of the tensiometers 62, 64, 66 is offset relative to the first object 12.

[0058] Each tensiometer 62, 64, 66 is designed to produce a voltage signal representative of the voltage across the terminals of the associated coil 56, 58, 60 and therefore representative of the component of the ambient magnetic field along one of the measurement directions cb c2, c3.

[0059] For this purpose, each tensiometer 62, 64, 66 is typically made in the form of an amplifier coupled to an analog-to-digital converter (ADC).

[0060] The voltage signals produced by the tensiometers 62, 64, 66 together form a measurement signal. Since the measurement directions cb c2, c3 are not coplanar, this measurement signal is representative of the ambient magnetic field at the receiver 24.

[0061] The receiving device 23 further includes a processing unit 70 for processing the measurement signal so as to deduce the local components, at the level of the receiver 24, along the three axes of the secondary frame (R2), of the resulting alternating magnetic fields. Such treatment is known in itself and does not will not be detailed further. It typically involves demodulating the measurement signal with demodulation signals at frequencies equal to those of the resulting alternating magnetic fields.

[0062] The determination system 20 further includes a data processing unit 80. This data processing unit 80 is configured to deduce from the measurements made by the receiving device 23 the position of the receiver 24 relative to the transmitter 22. It is also configured to deduce from this position the position of the first object 12 relative to the second object 14.

[0063] For this purpose, the data processing unit 80 is typically a computer. It comprises a processor or CPU (Central Processing Unit) 82 and a RAM (Random Access Memory) and / or ROM (Read Only Memory) type memory 84. The processor 82 is configured to execute instructions loaded into the memory 84. When the processing unit 80 is powered on, the processor 82 is able to read instructions from the memory 84 and execute them. These instructions form a computer program causing the processor 82 to implement processes 1000 ([Fig. 3]) and 2000 ([Fig. 10]), which will be detailed below.

[0064] The data processing unit 80 also includes a buffer memory 86 for the temporary storage of information necessary for the implementation of processes 1000 and 2000.

[0065] Alternatively, the data processing unit 80 is constituted by a microcontroller.

[0066] In the example shown, the data processing unit 80 is integrated into the second object 14. Alternatively (not shown), at least part of the data processing unit 80 is remote, for example, in a mobile terminal (not shown) and / or in a remote server (not shown). In other words, at least part of the steps of process 1000 and / or 2000 is carried out by a mobile terminal and / or a remote server. The second object 14 then includes a communication system, typically a wireless communication system, configured to send the data from the receiving device 23 to the mobile terminal and / or the remote server.

[0067] The assembly 10 also includes a disturbing element 90 suitable for impacting the measurement of the local components of the resulting alternating magnetic fields by the receiving device 23 when the disturbing element 90 is subjected to the alternating magnetic fields generated by the transmitter 22.

[0068] This disturbing element 90 consists of an object made of current-conducting material. It is configured so that the eddy currents induced by the emitter 22 in the disturbing element 90 are concentrated in a portion of surface 92 ([Fig.2]) of the disturbing element 90 delimiting a surface 94 of the disturbing element 90. In other words, any eddy currents induced by the emitter 22 in the disturbing element 90 that do not propagate in said portion of surface 92 are negligible. This portion of surface 92 preferably has a thickness at every point less than 5 mm, for example less than 2.5 mm, while being less than one-tenth, for example less than one-hundredth, of each other dimension of said portion of surface 92. To this end, the disturbing element 90 is typically constituted by a so-called thin plate, that is to say a plate whose thickness at every point is less than 5 mm, for example less than 2.5 mm, while being less than one-tenth, for example less than one-hundredth, of each other dimension of said plate; the portion of surface 92 is then constituted by the entire plate.Alternatively, the disturbing element 90 has a conductivity such that the skin thickness of the disturbing element 90 at any point at the frequencies of the alternating magnetic fields generated by the emitter 22 is less than 1.67 mm, for example less than 0.85 mm, while being less than one-third of one-tenth, for example one-third of one-hundredth, of each dimension of a surface of the disturbing element 90 oriented towards the emitter 22; the portion of surface 92 is then constituted by a portion of the disturbing element 90 extending over three times the skin thickness.

[0069] The surface portion 92 can be of any shape. In particular, it can be solid or perforated, flat or curved, and have edges of any shape. To keep the explanation simple and understandable, we will henceforth consider a surface portion 92 that is flat, solid (i.e., without a hole), parallelepiped-shaped, and of uniform thickness and conductivity (i.e., identical at every point), the disturbing element 90 typically being constituted by said surface portion 92 (i.e., consisting of a thin, flat, solid, parallelepiped-shaped plate). This surface portion 92 has a length P, a width ly, a thickness lz, and a conductivity 0.

[0070] A tertiary frame (R3) is attached to the disturbing element 90. This tertiary frame (R3) is a direct orthonormal frame with origin O3 formed of a triplet of axes, shown in [Fig. 1], comprising: - a first tertiary axis 83, - a second tertiary axis eZ orthogonal to the first axis 83, and - a third tertiary axis 83 orthogonal to the first and second axes 83 and e3.

[0071] Here, the third tertiary axis 83 is chosen so as to be orthogonal to the surface 94 of the perturbing element 90 delimited by the portion of surface 92. The first tertiary axis 83 is preferably chosen so as to be parallel to a first edge 96 of said surface 94 and the second tertiary axis is advantageously chosen so as to be parallel to a second edge 98 of said surface 94.

[0072] The pose of the perturbing element 90 relative to the first object 12 can thus be characterized by the set formed by the rotation matrix R[j3 from the primary frame (Ri) to the tertiary frame (R3) and the vector (expressed in the primary frame (RJ) going from the origin Oi of the primary frame (Ri) to the origin O3 of the tertiary frame (R3) and, conversely, the pose of the first object 12 relative to the perturbing element 90 can be characterized by the set formed by the rotation matrix R3j[ from the tertiary frame (R3) to the primary frame (RJ (equal to the transpose of the matrix R1>3) and the vector (expressed in the tertiary frame (R3)) going from the origin O3 of the tertiary frame (R3) to the origin Oi of the primary frame (Ri).Similarly, the pose of the perturbing element 90 relative to the second object 14 can be characterized by the set formed by the rotation matrix R23 from the secondary frame (R2) to the tertiary frame (R3) and the vector (exPrimc in the secondary frame (R2)) going from the origin O2 of the secondary frame (R2) to the origin O3 of the tertiary frame (R3) and, conversely, the pose of the second object 14 relative to the perturbing element 90 can be characterized by the set formed by the rotation matrix R3 2 of the frame. tertiary (R3) to secondary frame (R2) (equal to the transpose of the matrix R2 3) and of the vector (expressed in the tertiary frame (R3)) going from the origin O3 of the tertiary frame (R3) to the origin O2 of the secondary frame (Ri).

[0073] In a first embodiment of the assembly 10, the perturbing element 90 is fixed to the first object 12. The position of the perturbing element 90 relative to the first object 12 (and conversely the position of the first object 12 relative to the perturbing element 90) is then constant. In other words, the rotation matrices Ri 3 and R3 i and the vectors p^ and are constant.

[0074] In a second embodiment of assembly 10, the disruptive element 90 is attached to the second object 14. The placement of the disturbing element 90 relatively the position of the second object 14 (and conversely the position of the second object 14 relative to the perturbing element 90) is then constant. In other words, the rotation matrices R23 and R3 2 and the vectors p^ and IL are constant.

[0075] With reference to [Fig.2], the perturbing element 90 is modeled as a surface mesh 100. This surface mesh 100 is composed of a set of mesh points 102 regularly distributed over the surface 94 of the perturbing element 90 delimited by the portion of surface 92. These mesh points 102 are connected to each other by straight segments 104 delimiting between them two-dimensional mesh cells 106.

[0076] The straight segments 104 include external segments 107, each bordering one of the mesh cells 106 and thus delimiting a contour of the mesh 100. The straight segments 104 also include internal segments 108, each bordering two of the mesh cells 106 and thus forming the boundary between these cells 106. The external segments 107 and internal segments 108 together constitute the entirety of the straight segments 104; in other words, there is no straight segment 104 that does not constitute an external segment 107 or an internal segment 108.

[0077] For the surface mesh 100, a first direction di and a second direction d2, different from the first direction dp, are identified. The internal segments 108 comprise first internal segments 109 oriented along this second direction d2 and second internal segments 110 oriented along this first direction dp. The first internal segments 109 are thus all parallel to each other, and the second internal segments 110 are all parallel to each other. The first and second internal segments 109 and 110 together constitute the entirety of the internal segments 108; in other words, there is no internal segment 108 that is not a first internal segment 109 or a second internal segment 110.

[0078] Advantageously, the first direction di is contained in the plane defined by the first and third tertiary axes C*. ef. Here, the first direction di is in particular substantially parallel to the first tertiary axis ©3.

[0079] Preferably, the second direction d2 lies in the plane defined by the second and third tertiary axes 63. Here, the second direction d2 is in particular substantially parallel to the second tertiary axis

[0080] The mesh cells 106 are parallelogrammatic, specifically rectangular, here square. Each is bounded along the first direction di by a pair of second segments 110 of the surface mesh 100, and along the second direction d2 by a pair of first segments 109 of the surface mesh 100. Each mesh cell 106 has a first dimension along the first direction di and a second dimension along the second direction d2. Here, since the mesh cells 106 are square, these first and second dimensions are equal. Moreover, the mesh cells 106 are of identical size, so that the first and second dimensions of each cell 106 are equal to the first and second dimensions of every other cell 106. We will denote h the value of these first and second dimensions in what follows.

[0081] The mesh cells 106 together cover substantially the entirety of a relevant region of the surface 94, that is to say, at least 90% of the relevant region of the surface 94 is covered by the mesh cells 106. By "relevant region", we mean here and thereafter a region which affects the shape of the global eddy currents in the portion of surface 92.

[0082] In order to keep the exposition simple and understandable, the mesh cells 106 in the detailed example here are nine in number. However, in practice, the number of mesh cells 106 will be much higher and will typically be between 100 and 10000 cells 106.

[0083] This model of the disturbing element 90 is typically stored in the memory 84 of the data processing unit 80.

[0084] A method 1000 for tracking the positioning of the receiver 24 relative to the transmitter 22, implemented by the data processing unit 80, will now be described, with reference to [Fig.3].

[0085] This method 1000 includes a first step 1100 of providing calculation parameters, followed by a step 1200 of determining an initial pose at an initial determination time, itself followed by a repeated step 1300 of determining an updated pose at a later determination time.

[0086] Step 1100 of providing calculation parameters has two variants, depending on whether one considers the first or the second embodiment described above (disruptive element 90 attached to the first or second object 12, 14). A first of these variants, corresponding to the first embodiment (disruptive element 90 attached to the first object 12), will now be described, with reference to [Fig.4].

[0087] In this first variant, step 1100 begins with a substep 1110 for providing parameters of the disturbing element 90. These parameters of the disturbing element 90 typically include the surface mesh 100, the dimensions lx, ly, and lz of the surface portion 92, and the conductivity θ. They also include the position of the disturbing element 90 relative to the object, here the first object 12, to which it is attached. These parameters are typically stored in the memory 84 of the data processing unit 80.

[0088] Substep 1110 is followed by substep 1120 for determining the normal components of each of the alternating magnetic fields generated by the emitter 22, orthogonal to the surface 94 of the disturbing element 90 bounded by the surface portion 92, at several of the mesh points 102, in particular at each of the mesh points 102 not belonging to the contour of the mesh 100. These normal components are typically determined from the pose of each mesh point 102 relative to the emitter 22 (known because the pose of the disturbing element 90 relative to the first object 12 is predetermined), from a theoretical model of the alternating magnetic fields generated by the emitter 22, and of a correction algorithm for this theoretical model aimed at correcting, for example, the defects of emitter 22, using techniques known to a person skilled in the art.

[0089] This substep 1120 is itself followed by a substep 1130 of estimation of primary current density components in the surface portion 92 of the disturbing element 90 at the level of the internal segments 108 of the surface mesh 100, said primary current densities being eddy current densities induced in the disturbing element 90 by the magnetic fields generated by the emitter 22.

[0090] The primary current density components thus estimated consist of first components, along the first direction di of the mesh 100, of first primary current densities at first primary evaluation points Of,..., (Figure 2), and second components, along the second direction d2 of the mesh 100, of second primary current densities n equal to the number of second internal segments 110) at second primary evaluation points Og ([Fig. 2]). They are estimated for each of the alternating magnetic fields generated by the emitter 22.

[0091] The first primary evaluation points Of,..., Og belong in particular to the first internal segments 109 and are typically constituted by the midpoints of these first internal segments 109. The second primary evaluation points Op..., Og belong in particular to the second internal segments 110 and are typically constituted by the midpoints of these second internal segments 110.

[0092] Primary current densities are typically volume current densities.

[0093] These primary current density components are typically obtained at the means of the following formula or an equivalent formula:

[0094]

[0095] where jx is a vertical matrix of the set of the first components of the The first primary current densities j^,..., jx for an alternating magnetic field j generated by the emitter 22, is a vertical matrix of the set of second components of the second primary current densities ^,..., f for the alternating magnetic field j, is a vertical matrix of the normal components of the alternating magnetic field j determined in substep 1120, and jçd is the inverse matrix of a matrix Kj constructed from the local Maxwell equations

[0096]

[0097]

[0098] discretized between the estimated primary current density components and the normal components of the alternating magnetic field j. In the detailed (and simplified) example given here, the matrix Ky is given by the following formula: (1 1 JW / \ ^C x + r^N x M x ôrCy + ï-JJ-NyMy Dx Dy t Or: - 1 is imaginary such that = _ ; - is the angular frequency (i.e., angular frequency) of the magnetic field alternative j; - Dx and Dy are sparse matrices, each containing the factors +1 and -1, such that the following equality is satisfied: jx → pj j? - q of JJ in order to reflect the result of applying the divergence operator to the two terms of the Maxwell-Ampère equation in quasi-stationary regime on each of the mesh cells 106; - Mx is a square matrix whose: • the diagonal terms (Mx) are equal to J where is a prismatic volume of the surface portion 92 which extends over the entire thickness of the surface portion 92 and whose base is a parallelogram, centered on one of the first primary evaluation points O,..., Oq, whose edges are parallel to the first and second directions dh d2 and which extends, along the first direction db from the center of one of the mesh cells 106 framing said evaluation point Oj to the center of the other mesh cell 106 framing said evaluation point Oj and, along the second direction d2, from one of the segments 104 framing said evaluation point Oj to the other segment 104 framing said evaluation point Oj, and • the other terms (Mx) (with k # 7) are equal to II yk II —L_ where I^x^xl is the norm 2 of a vector going from one of the first I^k^l I primary evaluation points Of,..., Oq to another Oj of the first primary evaluation points O y..., Oq and || y^|| is the volume of the prismatic volume yk centered on the evaluation point O^; My is a square matrix whose: • the diagonal terms (My) are equal to J, A where is a prismatic volume of the surface portion 92 which extends over the entire thickness of the surface portion 92 and whose base is a parallelogram, centered on one of the second primary evaluation points 0,---, Oq, whose edges are parallel to the first and second directions dh d2 and which extends, along the first direction db from one of the segments 104 framing said evaluation point to the other segment 104 framing said evaluation point o^, and, following the second direction d2, from the center of one of the mesh cells 106 framing said evaluation point to the center of the other mesh cell 106 framing said Op evaluation point and the other terms (Af v) J' 1k (with k^l) are equal to il y^ll__1_. 11 where is the norm 2 of a vector going from one C'y of the II K second primary evaluation points Op..., Oq to another of the second primary evaluation points Op..., Oq and || y^|| is the volume of the prismatic volume y^ centered on the evaluation point Qr, Cy, and Ny are hollow matrices, each gathering the factor +1 and -1 so that the following equality is satisfied: \ JJ f \ XJJ / / in order to reflect the result of applying the Maxwell-Faraday law in quasi-steady state at each of the mesh points 102, the term A y jx । j representing the (discretized) curl of the current density at each of the mesh points 102 and the term / kj ivr î1 jM M representing the self-induced magnetic field 4n.h ' ^y^yJj J by the perturbing element at each of the mesh points 102, resulting from the application of the rotational operator to the vector potentials of this self-induced magnetic field obtained by application of the Biot and Savart law at each of the evaluation points of a primary current density.

[0099] In particular, taking the ordering of the first primary evaluation points Of,..., Og and the second primary evaluation points Oq shown in Figure 2, the matrices Dx, Dy, Cx, Cy, Nx and Ny can be written:

[0100] 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 -1 0 0 1 0 0 D, = ^x 0 -1 0 0 1 0 0 0 -1 0 0 1 0 0 0 -1 0 0 0 0 0 0 -1 0

[0101] 1 0 0 0 0 0| -1 1 0 0 0 0 0 -1 0 0 0 0 0 0 1 0 0 0 Dy = 0 0 -1 1 0 0 0 0 0 -1 0 0 0 0 0 0 1 0 0 0 0 0 -1 1'

[0102] 1 1 -1 0 0 0 0 0 1 -1 0 0 0 Q = 0 0 0 1 -1 0 0 0 0 0 1 -1 /

[0103] -1 0 1 0 o en 0 -1 0 1 0 0 Cy = 0 0 -1 0 1 0 0 0 0 -1 0 1 /

[0104] H -1 0 0 0 0 0 1 -1 0 0 0 Ny = 0 0 0 1 -1 0 \0 0 0 0 1 -1 /

[0105] pi 0 1 0 0 0\ 0 -1 0 1 0 0 Ny = 0 0 -1 0 1 0 0 0 0 -1 0 1 /

[0106] It should be noted that the number of rows of the matrices Dx and Dy is not equal to the number of 106 mesh cells (8 rows instead of 9) due to the application of the divergence operator to the two terms of the Maxwell-Ampère equation in the regime quasi-stationary on each of the 106 mesh cells leads to redundant equations.

[0107] For further information on how to determine the Kj matrix and the underlying theory, the reader may refer to the articles J. Nagel, "Fast finite-difference calculation of eddy currents in thin metal sheets," Applied Computational Electromagnetics Society Journal, 33:575-584, June 2018, and J. Nagel, "Finite-difference simulation of eddy currents in nonmagnetic sheets via electric vector potential," IEEE Transactions on Magnetics, 55:1-8, December 2019, the contents of which are incorporated herein by reference.

[0108] It should be noted that an example of a formula equivalent to the formula

[0110] given above is the following formula: [YES]

[0112] where Tj is the vector of a potential such that: the first components of the first primary current densities constitute the variation along the second direction d2 of the mesh 100 at the mesh points 102 where the normal components of the alternating magnetic fields generated by the emitter 22 are determined, and the opposites of the second components of the second primary current densities constitute the variation along the first direction di at said mesh points 102.

[0113] Substep 1130 is followed by substep 1140 of determination of secondary current densities in the portion of surface 92 at secondary evaluation points sb ..., Si2 ([Fig.2]), said secondary current densities remaining densities of the eddy currents induced in the disturbing element 90 by the magnetic fields generated by the emitter 22.

[0114] These secondary evaluation points sb ..., Si2 are regularly distributed within each mesh cell 106. In other words, each mesh cell 106 comprises at least one secondary evaluation point sb ..., s[2 and, for each mesh cell 106, the secondary evaluation points sb ..., Si2 contained within said mesh cell 106 are regularly distributed within it. The distribution of said secondary evaluation points sb ..., s[2 may nevertheless vary from one mesh cell 106 to another, i.e., the number of secondary evaluation points sb ..., Si2 within a mesh cell 106 may be different from the number of secondary evaluation points sb ..., Si2 within of another mesh cell 106. There can thus be mesh cells 106 with a single secondary evaluation point sb ..Si2, constituted by the center of the mesh cell 106, and other mesh cells 106 containing several secondary evaluation points sb ..Si2; a particular example of this is given on [Fig.2], where all mesh cells 106 include a single secondary evaluation point sb ..., si2, with the exception of the central mesh cell 106, which includes four secondary evaluation points s5, s6, s7, s8.

[0115] Preferably, a secondary evaluation point density sb ..., Si2 will be chosen in the mesh cell(s) 106 closest to the receiver 24 and / or in the mesh cell(s) 106 where the current densities are highest.

[0116] Each secondary evaluation point sb ..., Si2 is associated with an elementary cell 120 centered on said secondary evaluation point sb ..., Si2. Each mesh cell 106 consists of the elementary cell(s) associated with the secondary evaluation point(s) sb ..., Si2 included by said mesh cell 106. Thus, for each mesh cell 106 comprising a single secondary evaluation point sb ..., s[2, the elementary cell 120 associated with this secondary evaluation point sb ..., s[2 is constituted by the mesh cell 106 and, for each mesh cell 106 comprising several secondary evaluation points sb ..., Si2, the elementary cells 120 associated with these secondary evaluation points sb ..., s H are each constituted by a subdivision of the mesh cell 106 and cover the entirety of it.

[0117] These secondary current densities are volume current densities represented by a 3-dimensional vector. They are obtained by interpolating the components of the primary current densities. The interpolation is, for example, a linear interpolation. We then have:

[0119] where: - js is a secondary current density, measured at a secondary evaluation point sb ..., Si2 of a 106 mesh cell, expressed in the tertiary frame (R3), - jx and ix are the first components of first densities of Jk Jk+1 primary current in segments 104 of the mesh 100 bordering the cell of mesh 106 along the first direction db these first components being constituted by the first components determined during substep 1130 for the first primary current densities evaluated in internal segments 108, and being zero for the first primary current densities evaluated in external segments 107, jy and jy are the second components of second primary current densities in segments 104 of the mesh 100 bordering the mesh cell 106 along the second direction d2, these second components being constituted by the second components determined during substep 1130 for the second primary current densities evaluated in internal segments 108, and being zero for the second primary current densities evaluated in external segments 107, and Ki, K2, K3 and K4 are interpolation coefficients.

[0120] For example:

[0122] with K and K' being interpolation coefficients, these interpolation coefficients K, K' preferably being barycentric weights such as: where 1253*5.| is the magnitude of the vector going from point 1UJ evaluation of one of the first components of the first densities of Jk primary current in segments 104 of the mesh 100 bordering the mesh cell 106 along the first direction di to the secondary evaluation point Si and | is the norm of the vector going from the evaluation point on the other hand îx of the first components of first current densities Jk+1 primary in segments 104 of the mesh 100 bordering the mesh cell 106 along the first direction di to the secondary evaluation point s;, and imf 7sJ , °ù I ^y^ I is the norm of the vector going from the point mj evaluation of one of the second components of second primary current densities in segments 104 of the mesh 100 bordering the mesh cell 106 along the second direction d2 at the secondary evaluation point s; and I 7 I is the magnitude of the vector going from the point 1^7+1^1 7+1 evaluation of the other jy of the second components of second J1+1 primary current densities in segments 104 of the mesh 100 bordering the mesh cell 106 along the second direction d2 at the secondary evaluation point S;.

[0123] The secondary current densities are thus determined prior to the implementation of steps 1200, 1300 of determining an initial pose and of determining an updated pose.

[0124] A second variant of step 1100, corresponding to the second embodiment (disruptive element 90 attached to the second object 14), will now be described, with reference to [Fig.5].

[0125] In this second variant, step 1100 always begins with step 1110, which provides parameters for the disturbing element 90. As in the first variant, these parameters of the disturbing element 90 typically include the surface mesh 100, the dimensions lx, ly, and lz of the surface portion 92, the conductivity θ, and the pose of the disturbing element 90 relative to the object to which it is attached. However, unlike the first variant, this pose is not relative to the first object 12, since the disturbing element 90 is not attached to it, but relative to the second object 14.

[0126] This second variant differs further from the first variant in that step 1100 does not include substeps 1120, 1130 and 1140. Indeed, since the position of the disturbing element 90 relative to the first object 12 is not known, it is not possible to determine the normal components of each of the alternating magnetic fields generated by the emitter 22 at several of the mesh points 102. It is a fortiori not possible to determine the current densities of the eddy currents induced by these magnetic fields.

[0127] In this second variant, step 1100 nevertheless includes a substep 1150 for determining each matrix (or each matrix K^) for calculating the primary current density components from the normal components of each of the alternating magnetic fields generated by the emitter 22 at each of the mesh points 102 not belonging to the contour of the mesh 100, and a substep 1160 for determining an interpolation matrix grouping the interpolation coefficients for calculating the secondary current densities from the primary current density components. It is indeed clear that the interpolation equations presented in substep 1140 can be written in matrix form so as to give all the secondary current densities from all the primary current density components in a single operation. It is also clear that: - the Kj matrices depend only on the frequencies of the magnetic fields generated by the emitter 22 and on the parameters of the disturbing element 90 (in particular its surface mesh 100) and can therefore be calculated (as well as their inverses) without knowing the position of the disturbing element 90 relative to the emitter 22, - the interpolation coefficients depend exclusively on parameters of the surface mesh 100 (in particular the distances between the primary and secondary evaluation points) and can therefore be calculated (as well as the matrix grouping them) without knowing the pose of the disturbing element 90 relative to the emitter 22.

[0128] It also includes a substep 1170 for determining a calculation matrix for disturbing magnetic fields at the receiver 24. This matrix is ​​configured to allow the calculation, by a simple matrix operation, of local components, at the receiver 24, of disturbing magnetic fields produced by eddy currents induced in the disturbing element 90 by the magnetic fields generated by the emitter 22, when the secondary densities of these eddy currents are known. Indeed, it is known that eddy currents propagating in a current-conducting material such as that of the disturbing element 90 produce magnetic fields at the same frequency as said currents.In accordance with Biot and Savart's law, the local components of such magnetic fields at a measurement point are a function of these secondary current densities, the areas of the elementary cells 120 associated with the secondary evaluation points sb ..., Si2 of these secondary current densities, the average thicknesses of the surface portion 92 of the disturbing element 90 at the level of these elementary cells 120, and the vectors linking said measurement point to the secondary evaluation points sb ..., Si2. The position of the disturbing element 90 relative to the second object 14 (and therefore relative to the receiver 24) being known, as well as the dimensions of the disturbing element 90 and the parameters of the surface mesh 100, it is therefore possible to establish a matrix function linking the local components of such a magnetic field at the receiver 24 to the secondary eddy current densities induced in the disturbing element 90.

[0129] Thus, in substep 1170, a computation matrix £2 is determined such that:

[0130] B(r)=QJs

[0131] where: - J. is the following matrix: J s Asj JW j3^ JW

[0132]

[0133] j\sN) j\sN) where each ;7 / c 1 is a 2-dimensional vector of components, according to the JX^i) first and second directions dh d2, of one of the secondary current densities, evaluated at a secondary evaluation point si5 of the eddy currents induced in the disturbing element 90 by an alternating magnetic field j, and N is the total number of secondary assessment points Si, and B(r) is a matrix of local components, at the level of receptor 24, disruptive magnetic fields produced by these eddy currents.

[0134]

[0135]

[0136]

[0137] The matrix £2 is typically given by the following formula: O = — 12 4n with Or : 0 -(rS / )3 ( r - s f)2 ( r ' s ù3 ' 0 Pq is the magnetic permeability of free space. is a vector of the coordinates of a secondary evaluation point s; in the tertiary frame (R3), S7 is the surface area of ​​elementary cell 120 associated with the secondary evaluation point Si, is an average thickness of the surface portion 92 of the perturbing element 90 at the level of the mesh cell 106 comprising the elementary cell 120, TV is the total number of secondary assessment points si5 r is a vector of the coordinates of the position of receiver 24 in the coordinate system tertiary (R3), | • | is standard 2, and is the k component of the vector shown in parentheses.

[0138] These substeps 1150, 1160 and 1170 can be implemented in parallel or sequentially.

[0139] Step 1200 of determining an initial pose of the second object 14 relative to the first object 12 will now be described, with reference to [Fig.6].

[0140] This step 1200 consists of a first iteration of a step 1201 of determining a precise pose of the second object 12.

[0141] This step of determining a precise pose 1201 begins with a substep 1210 of obtaining a raw magnetic measurement Mb produced by the receiving device 23 at the initial determination time. This raw magnetic measurement Mb is composed of measurements, made by the receiving device 23 at the initial determination time, of the local components, at the receiver 24, along the three axes 6, 62 of the secondary frame (R2), of the three resulting alternating magnetic fields constituting the ambient magnetic field. During this acquisition substep 1210, said raw magnetic measurement Mb is typically transmitted by the receiving device 23 and received by the data processing unit 80.

[0142] Substep 1210 is followed by substep 1220 for determining an imprecise pose p0 of the receiver 24 relative to the transmitter 22 at the initial determination time. During this substep 1220, said imprecise pose p0 is advantageously deduced, for this first iteration, from the raw magnetic measurement Mb, typically by means of a theoretical model of the alternating magnetic fields generated by the transmitter 22 and an algorithm for correcting this theoretical model aimed at correcting, for example, the defects of the transmitter 22 and those of the receiver 24, by techniques known to those skilled in the art.

[0143] Substep 1220 is followed by several substeps 1230, 1240, 1250, 1260, 1270 which will be repeated.

[0144] A first of these substeps 1230, 1240, 1250, 1260, 1270 is a substep 1230 of determining an approximate pose pijk, where k=l,.. .,m is an index incremented at each iteration of substep 1230. This approximate pose pi>k is a function of the imprecise pose p0. During the first iteration of substep 1230, this approximate pose, denoted p;> b, is constituted by the imprecise pose p0.

[0145] This substep 1230 is followed by a substep 1240 for estimating the contribution of the disturbing element 90 to the raw magnetic measurement Mb. Indeed, each disturbing magnetic field produced by the eddy currents induced in the disturbing element 90 by an incident magnetic field such as one of those generated by the emitter 22 has the same eigencharacteristics, and in particular the same differentiating characteristics, as this incident magnetic field. Consequently, each of the resulting alternating magnetic fields measured at the receiver 24 is formed not only from one of the alternating magnetic fields generated by the emitter 22, but also from the corresponding disturbing magnetic field re-emitted by the disturbing element 90. Therefore, the raw magnetic measurement Mb can be written Mp = Me+ Mp, where Me is the contribution of the emitter 22, i.e. the alternating magnetic fields generated by the emitter 22, to the raw magnetic measurement Mb, and Mp is the contribution of the disturbing element 90, i.e. the alternating disturbing magnetic fields re-emitted by the disturbing element 90, to the raw magnetic measurement Mb.

[0146] Substep 1240 has two variants, depending on whether one considers the first or second embodiment described above (disruptive element 90 attached to the first or second object 12, 14). A first of these variants, corresponding to the first embodiment (disruptive element 90 attached to the first object 12), will now be described, with reference to [Fig. 7].

[0147] In this first embodiment, substep 1240 begins with a substep 1241 for estimating local components, at the receiver 24, of disturbing magnetic fields produced by the disturbing element 90. These local components are estimated from the approximate pose pijk and the secondary current densities induced in the disturbing element 90 by the magnetic fields generated by the emitter 22. Since the approximate pose pi>k was determined during the previous substep 1230 and the secondary current densities induced in the disturbing element 90 by the magnetic fields generated by the emitter 22 were determined prior to the implementation of the initial pose determination step 1200, during the calculation parameter provision step 1100, substep 1241 can indeed be implemented without any other prior step.

[0148] In this substep 1241, the data processing unit 80 estimates, for each of the alternating magnetic fields generated by the emitter 22, the local components, at the receiver 24, of the disturbing magnetic field produced by the eddy currents induced in the disturbing element 90 by said generated alternating magnetic field, taking the approximate pose pijk as the pose of the receiver 24 in the primary frame (Ri). To this end, the data processing unit 80 applies the Biot-Savart law to the secondary current densities of these eddy currents determined in substep 1140 of step 1100, which provides the calculation parameters. These local components of the disturbing magnetic field are thus typically estimated using the following formula or an equivalent formula:

[0149] h V ç ~

[0150] where: bj(r) is a vector of said local components of the disturbing magnetic field, here expressed in the tertiary frame (R3), .Pq is the magnetic permeability of free space, each is a 3-dimensional vector of the components, in the tertiary frame (R3), of one of the secondary current densities, evaluated at a secondary evaluation point si5 Sj is the area of ​​the elementary cell 120 associated with the secondary evaluation point Si, is an average thickness of the surface portion 92 of the perturbing element 90 at the level of the mesh cell 106 comprising the elementary cell 120, Sj is a vector of the coordinates of the secondary evaluation point Sj in

[0151]

[0152]

[0153] the tertiary reference point (R3), - r is a vector of the coordinates of the approximate pose in the tertiary frame (R3) (such a vector can be easily determined by means of the knowledge of the approximate pose and the pose of the perturbing element 90 in the primary frame (Ri)), - is the cross product operator, and - | • | is standard 2. Substep 1241 is followed by substep 1242 of deduction of the contribution of the disturbing element 90 to the raw magnetic measurement Mb from the local components of the disturbing magnetic fields thus estimated. During this substep 1242, the data processing unit 80 projects these local components onto the measurement directions cb c2, c3 of the receiver 24, taking into account the phase shifts between these local components and the demodulation signals used by the receiving device 23 to demodulate the measurement signal produced by the receiver 24. The projection of the local components onto the measurement directions cb, c2, c3 of receiver 24 is typically obtained by applying to said local components a rotation matrix defining the orientation of the disturbing element 90 relative to the approximate pose, constituted by the matrix R2 3 (assuming the pose of the secondary frame (R2) in the primary frame (RJ) equals the approximate pose). This matrix R2 3 can be easily obtained by the product , where r is the transpose of the rotation matrix defining the approximate pose in the primary frame (Ri).

[0154] The phase shifts between these local components and the demodulation signals are typically taken into account by temporal projection, onto the demodulation direction, of the spatial projection of said local components onto the measurement directions cb c2, c3. In other words, for each spatial projection of local components on one of the measurement directions cb c2, c3 of receiver 24, the ,erme Bfcos^ où : - is the amplitude of the projection of the local components of a field magnetic disturbance on a measurement direction cb and - (pJ is the phase shift, relative to the demodulation signals used by the receiving device 23 for demodulating the measurement signal produced by the receiver 24, of the projection of the local components of the disturbing magnetic field jgJ onto the measurement direction C;.

[0155] It should be noted that the amplitude can be intrinsically deduced from the projection of the local components of a disturbing magnetic field along a direction of measure C; As for the phase shift (pj, it is equal to the sum is a first phase shift between, on the one hand, the projection of the local components of the disturbing magnetic field ÿj onto the measurement direction c; and, on the other hand, the alternating magnetic field generated by the emitter 22 at the origin of the eddy currents inducing this field, which first phase shift can be intrinsically deduced from the projection of the local components of the disturbing magnetic field onto the direction of CB measurement and - ÔCPj is a second phase shift between, on the one hand, the demodulation signals and, on the other hand, the alternating magnetic field generated by the transmitter 22 which causes the eddy currents inducing the disturbing magnetic field pk which second phase shift is known thanks to a model of the transmitter 22 and, if necessary, thanks to a synchronization system of the receiving device 23 with the transmitter 22.

[0156] It should also be noted that, as an alternative, one can first make a time projection of the local components of the disturbing magnetic fields onto the demodulation direction, and then transport the results thus obtained into the secondary frame (R2) by applying the matrix R2 3 to the matrix of these results.

[0157] This yields an estimation matrix of a theoretical contribution of the disturbing element 90 to the raw magnetic measurement, the coefficients of which are made up of the terms A correction is then applied to this theoretical matrix of way to obtain a matrix jVf estimation of the actual contribution of the element The 90 interference factor in the raw magnetic measurement is typically adapted to simulate measurement errors due in particular to defects in receiver 24, using techniques known to those skilled in the art.

[0158] A second variant of substep 1240, corresponding to the second embodiment (disruptive element 90 attached to the second object 14), will now be described, with reference to [Fig.8].

[0159] In this second variant, substep 1240 comprises the same substeps 1241 and 1242 as described previously. These substeps 1241 and 1242 differ from those described above only in that the local components at the receiver 24 of the disturbing magnetic fields produced by the disturbing element 90 are estimated using the calculation matrix £2 determined during substep 1170 of step 1100, which provides the calculation parameters. The pre-existence of this calculation matrix £2 allows for a faster calculation of these local components once the secondary current densities of the eddy currents are known.

[0160] However, unlike the first variant described above, these secondary current densities were not determined prior to the implementation of substep 1240. Substep 1240 therefore also includes substeps 1243, 1244, 1245 and 1246 prior to the implementation of substep 1241 for the determination of said secondary current densities.

[0161] These preliminary substeps 1243, 1244, 1245, 1246 include a first substep 1243 of determining a pose of the disturbing element 90 relative to the emitter 22. This pose is determined from the approximate pose of the second object 14 relative to the first object 12 and the (known) pose of the disturbing element 90 relative to the second object 14.

[0162] This substep 1243 is followed by a substep 1244 for determining the normal components of each of the alternating magnetic fields generated by the emitter 22, orthogonal to the surface 94 of the disturbing element 90 delimited by the portion of surface 92, at several of the mesh points 102, in particular at each of the mesh points 102 not belonging to the contour of the mesh 100. These normal components are typically determined from the pose of each mesh point 102 relative to the emitter 22 (known from the fact that a pose of the disturbing element 90 relative to the first object 12 was determined in substep 1243), a theoretical model of the alternating magnetic fields generated by the emitter 22, and a correction algorithm for this theoretical model aimed at to correct, for example, the defects of emitter 22, by techniques known to a person skilled in the art.

[0163] This substep 1244 is itself followed by a substep 1245 of estimating components of the primary current densities of the eddy current densities induced in the disturbing element 90 by the magnetic fields generated by the emitter 22.

[0164]

[0165]

[0166] These primary current density components are typically obtained at the means of the following formula or an equivalent formula: .x Jj jy \JJ (J_z\ 0 J or Ji is a vertical matrix of the set of the first components of the The first primary current densities ix,..., jx for an alternating magnetic field j generated by the emitter 22, is a vertical matrix of the set of second components of the second primary current densities j^,..., jy for the alternating magnetic field j, is a vertical matrix of the normal components of the alternating magnetic field j determined during substep 1244, and jçl is the inverse matrix of the matrix Kj which was determined prior to the setting in work of the initial pose determination step 1200, during sub-step 1150 of step 1100 of providing calculation parameters.

[0167] The fact that the matrix jçl is precalculated allows a substantial time saving in the execution of this substep 1245.

[0168] Substep 1245 is followed by substep 1246 for determining the secondary current densities of the eddy currents induced in the disturbing element 90 by the magnetic fields generated by the emitter 22. As described for substep 1140, these secondary current densities are obtained by interpolating the components of the primary current densities. However, unlike substep 1140, this interpolation is performed using the interpolation matrix determined prior to the implementation of the initial pose determination step 1200, during substep 1160 of step 1100, which provides calculation parameters.

[0169] Again, the fact that this interpolation matrix is ​​pre-calculated allows a substantial time saving in the execution of this substep 1246.

[0170] Substep 1246 is followed by substeps 1241, then 1242.

[0171] Returning to Figure 6, substep 1240 is followed by substep 1250 for calculating a refined magnetic measurement. Put simply, in this substep 1250, the processing unit 80 subtracts the contribution of the disturbing element 90 from the raw magnetic measurement. In other words, it performs the following calculation: Ma = °where Ma is the refined magnetic measurement, is The raw magnetic measurement is the actual contribution of the disturbing element 90 estimated during substep 1240. It should be noted that this estimated contribution jÇfp differs from the actual contribution Mp because the pose of the second object 14 relative to the first object 12 is not precisely known. Thus, although constituting a better estimate of the contribution Me of the emitter 22 than the raw measurement M jr, the refined measurement Ma is not equal to this contribution Me, and a residual Me-Ma remains.

[0172] Substep 1250 is followed by substep 1260 of deduction of a refined pose pfjk, where k=l,...,m is an index incremented at each iteration of substep 1260. This refined pose pf>k is deduced from the refined magnetic measurement Ma, typically by means of a theoretical model of the alternating magnetic fields generated by the emitter 22 and an algorithm for correcting this theoretical model aimed at correcting, for example, the defects of the emitter 22 and those of the receiver 24, by techniques known to the person skilled in the art.

[0173] Substep 1260 is followed by substep 1270 of calculating a difference between the approximate pose pijk and the refined pose pfjk, here by means of a predefined function f, and of comparing this difference to a threshold, typically a predetermined threshold.

[0174] If this difference is greater than the threshold, then substeps 1230, 1240, 1250, 1260 and 1270 are repeated, the approximate pose pi>k+i being constituted, for the new iteration of these substeps, by the refined pose pf>k deduced at the end of the previous iteration.

[0175] If the difference is less than the threshold, then substep 1270 is followed by substep 1280 in which the precise pose (which constitutes, for this first iteration, the initial pose) is determined to be constituted by the refined pose pfjk deduced at the end of the last iteration of substep 1260.

[0176] Step 1300 of determining an updated pose will now be described, with reference to [Fig.9].

[0177] This step 1300 consists of a new iteration of the step of determining a precise pose 1201. Step 1300 therefore includes the same substeps 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280. Reference is made here to the description of these substeps given above.

[0178] This new iteration of the step for determining a precise pose 1201 differs from the first iteration only in the following characteristics: - the imprecise pose p0 determined during substep 1220 is preferably extrapolated from at least one previous pose determined during a previous iteration of the precise pose determination step 1201, rather than deduced from the raw magnetic measurement Mb (although such an implementation of substep 1220 remains possible, in a non-preferred embodiment of the invention), and - Logically, the precise position determined during sub-step 1280 constitutes the updated position and not the initial position.

[0179] A calibration method 2000 for the tracking system 20, implemented by the data processing unit 80, will now be described, with reference to [Fig. 10].

[0180] This method 2000 includes a first step 2100 of providing initial hypothetical values ​​for calibration parameters of the tracking system 20. These calibration parameters typically include one or more of the following parameters: - the positioning parameters of the disturbing element 90 relative to those of objects 12, 14 with respect to which it is fixed, - geometric parameters of the perturbing element 90, such as dimensional parameters of its surface portion 92, parameters of concavity or convexity, parameters of holes formed in the surface portion 92, etc. - physical parameters of the disturbing element 90, such as the conductivity (possibly non-uniform and anisotropic) of the surface portion 92.

[0181] These hypothetical values ​​of the calibration parameters are for example stored in memory 84 of the data processing unit 80.

[0182] Other parameters of the tracking system 20 are, on the other hand, predetermined and are not subject to calibration.

[0183] Step 2100 is followed by a pose tracking step 2200 based on hypothetical values. During this step 2200, the pose of the second object 14 relative to the first object 12 is tracked using the pose tracking method 1000, for a limited number of preselected poses, using the hypothetical values ​​of the calibration parameters as the parameters of the disturbing element 90. For this implementation of the pose tracking method 1000, it is preferable that the imprecise pose p0 determined during each iteration of substep 1220 be systematically deduced from the raw magnetic measurement Mb or, alternatively, deduced from a measurement supplied by an external sensor (not shown), for example an imaging motion capture system.

[0184] Step 2200 is followed by a step 2300 of calculation, for each determined pose, of a criterion representative of an error on the determined pose.

[0185] This error criterion is preferably chosen so as to be an increasing function of the likelihood of an error on the determined pose.

[0186] According to a first variant of this step 2300, this criterion is a function of the diagonal entries of a diagonal matrix S of positive real values ​​Su, S22 and S33 ordered such that > ​​S22 — S33. This diagonal matrix results from a singular value decomposition of a corrected matrix Mc deduced from the matrix of the refined measure Ma obtained during the last iteration of substep 1250, said corrected matrix Mc being obtained by means of a correction algorithm aimed at correcting the errors of the refined measure Ma related, for example, to defects in the transmitter 22 and those in the receiver 24 by techniques known to those skilled in the art. This singular value decomposition is such that the diagonal matrix S satisfies the following relation:

[0187] Mc = PSQT

[0188] where P is a first matrix belonging to the SO(3) group of rotation matrices in dimension 3 and QT is the transpose of a second matrix Q also belonging to the SO(3) group.

[0189] The error criterion is thus, for example, determined by the following formula: [°190] _ / \1 / 2 1 S?1+4S12+4S23 /

[0191] According to a second variant of this step 2300, the representative criterion of an error on the determined pose is the variation of the gap between the pose determined for the receiver 24 and the pose determined, using the same tracking method 1000 and with the same hypothetical values ​​of the calibration parameters, for another receiver (not shown), fixed relative to the receiver 24.

[0192] Step 2300 is followed by a step 2400 of calculation of a sum of said error criteria on all the poses determined during substep 2200.

[0193] If this is the first iteration of step 2400, it is followed by a step 2700 of updating the hypothetical values, then steps 2200, 2300 and 2400 are repeated, the follow-up of the placement of step 2200 being carried out on the basis of the hypothetical values ​​thus updated.

[0194] If it is a subsequent iteration of step 2400, this is followed by a substep 2500 of calculating a gradient of said sum with respect to the variation of the hypothetical values, that is to say, a difference between said sum and the calculated sum During the previous iteration of substep 2400, the change in hypothetical values ​​between the two iterations was reported. This gradient is then compared to a threshold during step 2600. If it is above the threshold, steps 2700 and 2200 to 2600 are repeated. If it is below the threshold, step 2600 is followed by step 2800, which validates the hypothetical values ​​provided during the last iteration of step 2700.

[0195] The hypothetical discounted values ​​provided in step 2700 are preferably established on the basis of a gradient descent algorithm for which the function to be minimized is the sum determined in step 2400.

[0196] Thanks to the invention described above, it is possible to correct in time In reality, for example at a frequency of 40 Hz, measurements from an electromagnetic tracking system 20 are affected by a complex-shaped disturbing element 90, even for receiver 24 poses close to the disturbing element 90. The calculations required for this correction are indeed simple, allowing for rapid execution, yet still enable a relatively accurate modeling of the disturbing magnetic fields re-emitted by the disturbing element 90, thus giving good accuracy to the corrected measurement. Consequently, the accuracy of the electromagnetic tracking system 20 is improved.

[0197] This correction also makes it possible to improve the accuracy of the static measurements of the electromagnetic tracking system 20, since it directly corrects the measurements of the receiver 24. It also makes it possible to do without quasi-static phases of the receiver 24, since it does not require variation of the frequency of the alternating magnetic fields generated by the transmitter 22.

[0198] Finally, the invention allows for autonomous calibration of the electromagnetic tracking system 20 without requiring the use of a reference system. This calibration is also relatively fast, especially when compared to that of a mapping model.

Claims

1. Demands Method (1000) for tracking the placement of a receiver (24) relative to a transmitter (22), the transmitter (22) comprising at least two magnetic generators (30, 32, 34) each oriented along a proper direction (gb g2, g3) and capable of generating an alternating magnetic field having a distinctive characteristic enabling the alternating magnetic field to be distinguished from the alternating magnetic field generated by each other magnetic generator (30, 32, 34), the magnetic generators (30, 32, 34) being arranged relative to each other such that their directions (gb g2, g3) are non-coplanar, the method (1000) being implemented by a data processing unit (80) and comprising the following step: - determination (1201) of a pose of the receiver (24) relative to the transmitter (22) at a determination instant, said determination (1201) comprising the following sub-steps: • a) obtaining (1210) a raw magnetic measurement composed of measurements, carried out at the time of determination, of local components, along at least two non-coplanar measurement directions (cb c2, c3) attached to the receiver (24), of at least two resulting alternating magnetic fields constituting an ambient magnetic field, each of said resulting alternating magnetic fields having a distinctive characteristic identical to that of one of the generated alternating magnetic fields, • b) determination (1220) of an imprecise pose of the receiver (24) relative to the emitter (22) at the time of determination, • c) estimation (1240), based on an approximate pose that is a function of the imprecise pose, of a contribution from a disturbing element (90) to the raw magnetic measurement, • d) calculation (1250) of a refined magnetic measurement by subtracting the contribution of the disturbing element (90) from the raw magnetic measurement, and • e) deduction (1260) of a refined pose of the receiver (24) relative to the emitter (22) at the time of determination from the refined magnetic measurement, in which the disturbing element (90) is modeled as a surface mesh (100) composed of a set of mesh points (102) regularly distributed over a surface (94) of the disturbing element (90) and connected to each other by straight segments (104) delimiting between them two-dimensional and parallelogrammatic mesh cells (106), and substep c) of estimation (1240) of the contribution of the disturbing element (90) includes, for each of the magnetic fields alternating currents generated, the estimation (1241) of local components, at the level of the receiver (24), for a pose of the receiver (24) equal to the approximate pose of a disturbing magnetic field produced by eddy currents induced in the disturbing element (90) by said generated alternating magnetic field, by applying the Biot-Savart law to secondary current densities of these eddy currents, said secondary current densities resulting from an interpolation of primary current density components of said eddy currents, said primary current density components being a function of a pose of the disturbing element (90) with respect to the emitter (22) and being constituted: - first components, along a first direction (di) of the surface mesh (100), of first primary current densities in first internal segments (109) of the surface mesh (100) oriented along a second direction (d2) of the surface mesh (100), in particular in the middle of said first internal segments (109), and - of second components, along the second direction (d2) of the surface mesh (100), of second primary current densities in second internal segments (110) of the surface mesh (100) oriented along the first direction (di) of the surface mesh (100), in particular in the middle of said second internal segments (110).

2. A tracking method (1000) according to claim 1, wherein the primary current density components are obtained by means of the following formula or an equivalent formula: where jX is a vertical matrix of the set of first Jj components of the first primary current densities, is a vertical matrix of all second components of the second secondary current densities, hz- is a matrix J vertical of normal components of the alternating magnetic field generated, orthogonal to the surface (94) of the perturbing element (90), in each of the mesh points (102) not belonging to a contour of the surface mesh (100), and jçl is the inverse matrix of a matrix constructed from the local Maxwell equations discretized between the primary current density components and the normal components of the generated alternating magnetic field.

3. A tracking method (1000) according to claim 1 or 2, wherein the disturbing element (90) is fixed relative to the receiver (24) and substep c) of estimating (1240) the contribution of the disturbing element (90) comprises determining (1243) a pose of the disturbing element (90) relative to the emitter (22), estimating (1245) the components of the primary current densities as a function of said pose of the disturbing element (90) relative to the emitter (22), and determining (1246) the secondary current densities by interpolating the components of the primary current densities.

4. A tracking method (1000) according to claims 2 and 3 taken together, wherein the matrix jçd is determined prior to the implementation of the pose determination step (1201).

5. A tracking method (1000) according to claim 1 or 2, wherein the disturbing element (90) is fixed relative to the emitter (22) and the secondary current densities are determined prior to the implementation of the pose determination step (1201).

6. A tracking method (1000) according to any one of the preceding claims, wherein the pose determination step (1201) comprises several successive iterations of substeps c), d) and e), the approximate pose being constituted, for the first iteration, by the imprecise pose and, for each subsequent iteration, by the refined pose deduced at the end of the previous iteration.

7. A tracking method (1000) according to any one of the preceding claims, comprising several iterations of the pose determination step (1201) for determining poses of the receiver (24) at different determination times, the first iteration constituting a step (1200) for determining an initial pose of the receiver (24) and subsequent iterations constituting steps (1300) for determining an updated pose, the pose imprecise determined during the initial pose determination step (1200) being deduced from the raw magnetic measurement.

8. A tracking method (1000) according to claim 7, wherein the imprecise pose determined during at least one of the steps for determining an updated pose (1300) is extrapolated from at least one prior pose determined during a prior iteration of the pose determination step (1201).

9. A tracking method (1000) according to any one of the preceding claims, wherein the secondary current densities are determined at evaluation points (sb ..., Si2) regularly distributed within each mesh cell (106).

10. A tracking method (1000) according to any one of the preceding claims, wherein the local components of the disturbing magnetic field are estimated by means of the following formula or an equivalent formula: where bj(r) is a vector of said local components of the disturbing magnetic field, Pq is the magnetic permeability of free space, j(Si) is a 3-dimensional vector of the components, in a frame attached to the disturbing element (90), of one of the secondary current densities, Sj is the area of ​​an elementary cell (120) associated with said estimated surface current density, Q is an average thickness of a surface portion (92) of the disturbing element (90), within which the eddy currents are concentrated, at the level of the mesh cell (106) comprising the elementary cell (120), and is a vector of the coordinates of a center of said elementary cell (120) in a frame attached to the element disruptor (90),r is a vector of the coordinates of the approximate pose in said frame attached to the disturbing element (90), X is the cross product operator and | • | is the norm 2, the elementary cell (120) being constituted by one of the mesh cells (106) or by a subdivision of a subdivided mesh cell (106), each subdivided mesh cell (106) being made up of elementary cells (120), each associated with one of the secondary current densities, regularly distributed at,

11.

12.

13.

14. the interior of said mesh cell (120) subdivided and covering the entirety of it. Calibration method (2000) of an electromagnetic tracking system (20) comprising a step (2200) of tracking a pose of a receiver (24) of the electromagnetic tracking system (20) relative to a transmitter (22) of the electromagnetic tracking system (20) by means of a tracking method (1000) according to any one of the preceding claims, based on hypothetical values ​​of calibration parameters of the electromagnetic tracking system (20), and updating (2700) the hypothetical values ​​of the calibration parameters so as to minimize a sum of error criteria each representative of an error on a pose of the receiver (24) relative to the transmitter (22) determined during the tracking step (2200). Product computer program comprising code instructions for the implementation, by a processor, of a tracking method (1000) or a calibration method (2000) according to any one of the preceding claims. A computer-readable recording medium on which a computer program product according to claim 12 is stored. An assembly (10) comprising a reference object (12), a moving object (14), and an electromagnetic tracking system (20) for tracking the position of the moving object (14) relative to the reference object (12), said electromagnetic tracking system (20) comprising: - an emitter (22) integral with the reference object (12), said emitter (22) comprising at least two magnetic generators (30, 32, 34), each oriented along a proper direction (gb g2, g3) and capable of generating an alternating magnetic field having a distinctive characteristic enabling said alternating magnetic field to be distinguished from the alternating magnetic field generated by each other magnetic generator (30, 32, 34), the magnetic generators (30, 32, 34) being arranged relative to each other such that their directions (gb g2, g3) are non-coplanar, and - a receiving device (23) comprising a receiver (24) attached to the moving object (14), said receiving device (23) being capable of measuring local components, along at least two non-coplanar directions (cb c2, c3), of at least two resulting alternating magnetic fields constituting an ambient magnetic field, each of said resulting alternating magnetic fields having a distinctive characteristic identical to that of one of the generated alternating magnetic fields, so as to obtain a raw magnetic measurement, the assembly (10) also comprising a disturbing element (90) capable of impacting the measurement of the local components of the resulting alternating magnetic fields by the receiving device (23) when the disturbing element (90) is subjected to the alternating magnetic fields generated from the emitter (22), in which the electromagnetic tracking system (20) includes a data processing unit (80) configured for the implementation of a tracking method (1000) or a calibration method (2000) according to any one of claims 1 to 11.

Images