The instant application claims the benefit of French patent application number 15/56827, filed Jul. 20, 2015, the teachings of which are incorporated herein by this reference.
The present application relates to electromagnetic position tracking systems or electromagnetic locating systems. It more specifically aims at electromagnetic position tracking systems adapted to the medical field, for example, to locate one or a plurality of tools used by a surgeon during a medical operation.
Position tracking systems of the type comprising the following elements are here more particularly considered:
an electromagnetic field emission device comprising one of a plurality of emitter coils; and
a device for receiving this field comprising one or a plurality of receiver coils.
In such a system, the conjoint analysis of the fields emitted by the emitter coils and of the fields received by the receiver coils enables to determine the position and/or the orientation of the receiver device with respect to the emission device.
An example of such a system is for example described in patent application US2005/0165297.
An electromagnetic position tracking system having an improved performance over known systems, particularly in terms of location accuracy, is needed.
Thus, an embodiment provides an electromagnetic position tracking system comprising: an electromagnetic field emission device comprising at least one emitter coil and, connected to the emitter coil, a first generator of an electric signal for exciting the emitter coil; a device for receiving the electromagnetic field emitted by the emission device, comprising at least one receiver coil, and, connected to the receiver coil, a circuit for reading an electric signal induced in the receiver coil; and a system for measuring at least one parameter of the reception device, comprising, connected to the receiver coil, a second generator of an electric signal for exciting the receiver coil.
According to an embodiment, the parameter is a parameter from the group comprising the series resistance and the self-inductance of the receiver coil.
According to an embodiment, the reception device comprises at least two receiver coils, and the parameter is a parameter from the group comprising the series resistance of each receiver coil, the self-inductance of each receiver coil, and the coefficient of mutual inductance between two receiver coils.
According to an embodiment, the parameter comprises the series resistance of each receiver coil, and the measurement system is capable, to measure the series resistance, of applying to one end of the receiver coil a DC voltage or a sinusoidal voltage having a frequency at least 20 times lower than the cut-off frequency of the coil.
According to an embodiment, the parameter comprises the self-inductance of each receiver coil, and the measurement system is capable, to measure the self-inductance, of applying on one end of the receiver coil a sinusoidal voltage having a frequency greater than the cut-off frequency of the coil.
According to an embodiment, the parameter comprises the coefficient of mutual inductance between the two receiver coils, and the measurement system is capable, to measure the coefficient, of applying on one end of one of the two receiver coils a sinusoidal voltage having a frequency greater than the cut-off frequency of the coil.
According to an embodiment, the read circuit comprises, associated with each receiver coil, a sense amplification stage comprising an operational amplifier having its inverting input coupled to the output via a feedback resistor and having its non-inverting input coupled to a node of application of a reference potential, a first end of the receiver coil being connected to the inverting input of the operational amplifier.
According to an embodiment, the second generator comprises, associated with each receiver coil, a controllable voltage source connected between the second end of the receiver coil and the node of application of a reference potential.
According to an embodiment, the second generator comprises a digital frequency synthesizer.
According to an embodiment, the emitter coil comprises three emitter coils oriented along different axes, and the receiver coil comprises three receiver coils oriented along different axes.
The foregoing and other features and advantages will be discussed in detail in the following non-limiting description of specific embodiments in connection with the accompanying drawings.
The same elements have been designated with the same reference numerals in the different drawings. For clarity, only those elements which are useful to the understanding of the described embodiments have been shown and are detailed. In particular, the calculation algorithms enabling, based on the values of the fields emitted by the emitter coils and received by the receiver coils, to determine the position and/or the orientation of the reception device with respect to the emission device, are not detailed. Indeed, the described embodiments are compatible with usual algorithms for determining the position and/or orientation information of the reception device based on the measured field values. Unless otherwise specified, expressions “substantially” and “in the order of” mean to within 10%, preferably to within 5%. Further, in the present description, term “connected” will be used to designate a direct electric connection, with no intermediate electronic component, for example, by means of one or a plurality of conductive tracks or conductive wires, and term “coupled” or term “linked” is used to designate either a direct electric connection (then meaning “connected”) or a connection via one or a plurality of intermediate components (resistor, capacitor, etc.).
The system of
The system of
The system of
The locating system of
Assembly 103 of emitter coils of emission device 101 defines a reference frame for the location of assembly 113 of receiver coils of reception device 111. Assembly 113 of receiver coils may be secured to an object to be located. As an example, in the case of an application to the locating of a surgical tool, assembly 103 of emitter coils may be secured to a patient's body or on a structure external to the patient's body used as a reference frame. The assembly 113 of receiver coils may be secured to a surgical tool which is desired to be located during the operation.
During a phase of tracking the position of the object to be located, emitter coils BE1, BE2 and BE3 of emission device 101, controlled by generator 105, each emit an AC electromagnetic field. As an example, generator 105 is configured to simultaneously apply excitation signals of different frequencies to the different coils of assembly 103, in which case the coils of assembly 103 simultaneously emit AC electromagnetic fields having different frequencies. This enables, at the level of the reception device, to be able to spectrally discriminate the fields originating from the different emitter coils. As a variation, the generator is configured to sequentially apply to the different coils of assembly 103 excitation signals of same frequency, so that at a given time, only one of the coils of assembly 103 emits an electromagnetic field. This enables, at the level of the reception device, to be able to discriminate in time the fields originating from the different emitter coils. The emission frequency or frequencies of device 101 are for example in the range from 1 to 50 kHz. Signals (for example digital) representative of the excitation signals applied by generator 105 to coils BE1, BE2 and BE3, are transmitted to processing unit 121.
At the level of reception device 111, the electric signals are induced in receiver coils BR1, BR2 and BR3 under the effect of the electromagnetic fields emitted by emitter coils BE1, BE2 and BE3. These signals are measured by read device 115 and transmitted (for example, in digital form) to processing unit 121.
Processing device 121 comprises a synchronous demodulation and filtering unit (not shown in the drawing) capable, based on signals transmitted by generator 105 and read device 115, of determining, for each of the receiver coils of assembly 113, the intensity of the field captured by the coil from each of the emitter coils of assembly 103. Based on the intensity of the fields captured by the receiver coils, processing device 121 determines the position and/or the orientation of the assembly of receiver coils 113 with respect to the assembly of emitter coils 103. More particularly, in the example of
The accuracy of a locating system of the type described in relation with
The case where read device 115 of reception device 111 comprises a transimpedance-type amplification device is here more particularly considered. An example of such an amplification device is illustrated in
In the example of
Coil BRi has its two ends connected, for example, via a coaxial cable 205i, respectively to the inverting input (−) and to the non-inverting input (+) of operational amplifier 203i.
In such an assembly, coil BRi is shorted due to the virtual ground present on the inverting input (−) of operational amplifier 203i. Coil BRi conducts a current iBRi proportional to the intensity of the magnetic field received by the coil.
In operation, operational amplifier 203i delivers on its output node a voltage Vouti (referenced with respect to node GND) which is an amplified image of current iBRi.
Read device 115 may further comprise, associated with each of receiver coils BRi, a sampling and digitization circuit (not shown) capable of sampling and digitizing voltage Vouti, the digital output of this circuit being coupled to processing device 121.
The transfer function or magnetic field/voltage conversion gain of such an assembly may be expressed as follows:
where B designates the intensity of the magnetic field received by coil BRi, j designates the imaginary unit, ω designates the pulse of field B (ω=2πf, f being the frequency of field B), MBRi designates the magnetic moment per current unit of coil BRi, in A·m2/A, LBRi designates the inductance of coil BRi, RsBRi designates the series resistance of the circuit connected between the inverting input (−) and the non-inverting input (+) of operational amplifier 203i, particularly comprising, in this example, the resistance of coil BRi and the resistance of cable 205i, and KBRI(jω) designates the field/voltage proportionality coefficient at pulse ω.
The frequency response of conversion gain KBRi is a response of high-pass type, with a gain equal to (MBRi/LBRi)*Rci, and a cut-off frequency equal to RsBRi/(2πLBRi).
To accurately know the magnetic field/voltage transfer function of such an assembly, one needs to accurately know parameters LBRi and RsBRi of the assembly.
The accurate knowledge of the value of resistor RsBRi raises a specific issue since this value is likely to vary by relatively large proportions under the effect of the temperature variations to which receiver coils BRi may be submitted.
The value of inductance LBRi, although it is less dependent on temperature variations, however remains capable of drifting (particularly depending on the materials used, such as for example, the use or not of a ferromagnetic core). It may further be useful to accurately know this value in a phase of calibration of the electromagnetic system to constrain certain coefficients used by the calculation algorithm for determining the position and/or the orientation.
Further, in the case of a reception device with a plurality of coils of the type shown in
where Mli,i′ designates the coefficient of mutual inductance between coils BRi and BRi′.
The values of the mutual inductance coefficients are little dependent on temperature variations, and fluctuate little on the short term. However, in the same way as for inductance values LBRi, it may be useful to know the mutual inductance coefficients in the electromagnetic calibration phase to constrain certain coefficients used by the calculation algorithm for determining the position and/or the orientation.
The system of
The system of
The parameter measurement system comprises a controllable device 331 for generating electric signals (voltages or currents) for exciting the receiver coils of assembly 113. Generator 331 is connected to at least one terminal of each of receiver coils BR1, BR2 and BR3.
Generator 331 is coupled to processing device 121 by a wire or wireless connection. Generator 331 is capable of transmitting to processing device 121 signals (for example, digital) representative of the excitation signals that it applies to the terminals of receiver coils BR1, BR2 and BR3.
In the example of
In the example of
The other end of coil BRi is connected to the inverting input (−) of operational amplifier 203i.
For simplification, only one coil BRi, one amplification device 201i associated with this coil, and one controllable voltage source 433i associated with this coil have been shown in
As an example, a cable 205i, for example, a coaxial cable, couples the two ends of coil BRi respectively to the inverting input (−) of operational amplifier 203i and to the output of controllable voltage source 433i.
In this assembly, coil BRi is not directly connected to the non-inverting input (+) of operational amplifier 203i, but is coupled to the non-inverting input (+) via controllable current source 433i and reference node GND.
The operation of the electromagnetic locating system of
The operation of the system for measuring parameters of reception device 101 of the electromagnetic locating system of
A) measurement of series resistance RsBRi
To measure the series resistance RsBRi associated with a coil BRi, generator 331 controls the application, via controllable voltage source 433i, of a DC or sinusoidal voltage Vmi on the terminal of coil BRi connected to voltage source 433i.
In the case of a sinusoidal voltage Vmi, the latter may be expressed as follows:
Vmi(t)=Vref cos(2πfit)
where Vref designates the amplitude of the sinusoidal voltage and fi designates the frequency thereof.
Voltage Vmi is amplified by transimpedance amplification stage 201i, the amplified image of this voltage forming a component of output voltage Vouti of amplification stage 201i. Frequency fi of voltage Vmi is preferably selected to be much smaller than the cut-off frequency of the coil (equal to RsBRi/2πLBRi), so that the impedance of the inductance of coil BRi can be neglected as compared with the series resistance to be measured. As an example, frequency fi is selected to be at least 20 times smaller than the cut-off frequency of the coil. As an example, frequency fi may be in the range from 1 to 20 Hz for a cut-off frequency of the coil substantially equal to 1 kHz. Amplitude Vref of the sinusoidal voltage is selected to avoid saturating operational amplifier 203i, that is, so that output voltage Vouti of the amplification stage remains lower than the saturation output voltage of amplifier 203i.
Processing device 121 comprises a synchronous detection unit (not shown in the drawing) capable of extracting peak value Voutifi of the frequency component fi of signal Vouti. Based on peak value Voutifi, processing device 121 may determine the value of series resistance RsBRi by the following formula:
In the case of a DC voltage Vmi of value Vref, only the processing implemented by the detection unit differs from what has just been described. This processing may be a simple low-pass filtering enabling to directly extract DC component Voutic from voltage Vouti. Processing device 121 can then determine the values of series resistance RsBRi by the following formula:
B) measurement of self-inductance LBRi
To measure the inductance LBRi of coil BRi, generator 331 controls the application, via controllable voltage source 433i, of a sinusoidal AC voltage Vmi on the terminal of coil BRi connected to voltage source 433i.
Voltage Vmi may, as in case A), be expressed as follows:
Vmi(t)=Vref cos(2πfit)
This time, frequency fi of voltage Vmi is selected to be close to or greater than the cut-off frequency of the coil (equal to RsBRi/2πLBRi), so that the impedance of the inductance of coil BRi is non-negligible as compared with the series resistance of coil RsBRi. As an example, to measure inductance LBRi, frequency fi is selected to be from 1 to 50 times the cut-off frequency of the coil. As an example, to measure inductance LBRi, frequency fi may be in the range from 1 to 20 kHz for a cut-off frequency in the order of 1 kHz. Preferably, frequency fi is different from the emission frequency or frequencies of emission device 101. As an example, frequency fi is smaller than the lowest emission frequency of emission device 101, or greater than the largest emission frequency of emission device 101. Amplitude Vref of voltage Vmi is selected to avoid saturating operational amplifier 203i.
Voltage Vmi is amplified by transimpedance amplification stage 201i. The amplified image of voltage Vmi thus forms a component of frequency fi of output voltage Vouti of amplification stage 201i, expressed as follows:
The detection unit of processing device 121 is capable of extracting peak value Voutifi of frequency component fi of signal Vouti, expressed as follows:
Based on this peak value, processing device 121 may determine the value of inductance LBRi by the following formula:
C) measurement of mutual inductance Mli,i′
To measure the coefficient Mli,i′ of mutual inductance between two receiver coils BRi and BRi′, generator 331 controls the application, via controllable voltage source 433i, of a sinusoidal AC voltage Vmi on the terminal of coil BRi connected to voltage source 433i. Voltage Vmi may, as in cases A) and B), be expressed as follows:
Vmi(t)=Vref cos(2πfit)
As in case B), frequency fi of voltage Vmi is selected to be close to or greater than the cut-off frequency of the coil, RsBRi/2πLBRi, so that the impedance of the inductance of coil BRi is non-negligible as compared with series resistance RsBRi. As an example, to measure the mutual inductance coefficient, frequency fi is selected to be from 1 to 50 times the cut-off frequency of the coil. As an example, to measure mutual inductance coefficient frequency fi may be in the range from 1 to 20 kHz for a cut-off frequency of the coil in the order of 1 kHz. Further, as in case B), frequency fi is preferably different from the emission frequency or frequencies of emission device 101. As an example, frequency fi is smaller than the lowest emission frequency of emission device 101, or greater than the largest emission frequency of emission device 101. Amplitude Vref of voltage Vmi is selected to avoid saturating operational amplifier 203i.
The voltage applied, via controllable current source 433i′, to the terminal of coil Biti′ connected to voltage source 433i′, is for example zero. Anyway, this voltage comprises no component of frequency fi.
Current iBRi crossing coil BRi is coupled to output voltage Vouti of transimpedance amplifier 201i by the following relation:
Under the effect of this sinusoidal current, coil BRi emits an AC electromagnetic field of frequency fi. This field is captured by coil BRi′, inducing across coil Biti′ a voltage vBRi′ expressed as follows:
vBRi′(t)=−jωMli,i′iBRi(t)
This voltage is amplified by transimpedance amplification stage 201i′, which outputs a voltage Vouti′ expressed as follows:
Processing device 121 is capable, based on the signals measured at the output of amplification stages 201i to 201i′, of determining coefficient Mli,i′ which can be expressed as follows:
As an example, the different mutual inductance coefficients of assembly 113 of receiver coils, that is, coefficients Ml1,2, Ml1,3 and Ml1,3 in the example of three receiver coils, may be determined simultaneously by using different excitation frequencies fi. As a variation, the different mutual inductance coefficients may be determined sequentially by using a same excitation frequency fi.
The measurement of parameters RsBRi, LBRi and/or Mli,i′ may be performed during a calibration phase prior to a position tracking phase, and/or during the actual position tracking phase. For example the measurements of parameters RsBRi, LBRi and/or Mli,i′ may be performed continuously in parallel with the measurements of the electromagnetic fields from the emission device, provided for the excitation frequency or frequencies applied by the parameter measurement system on the receiver coils to be different from the emission frequencies of the emission device. Further, parameters RsBRi, LBRi and/or Mli,i′ may be measured simultaneously as soon as the excitation frequencies used to measure them are different from one another.
As an example, generator 331 may comprise a digital frequency synthesizer (not shown) capable of generating one or a plurality of sinusoidal signals having controllable frequencies and, possibly, of adding these signals. The controllable voltage source 433i associated with each coil for example comprises a digital-to-analog converter (not shown) having its input coupled to a digital output of the frequency synthesizer, and having its output coupled to coil BRi. As a variation, a resistive voltage dividing bridge may be provided between the output of the digital-to-analog converter and the end of coil BRi connected to voltage source 433i, to adjust the level of the excitation voltage applied to coil BRi.
An advantage of the locating system described in relation with
Specific embodiments have been described. Various alterations, modifications, and improvements will readily occur to those skilled in the art.
In particular, the described embodiments are not limited to the above-mentioned application to the locating of surgical tools, but may apply to other fields requiring accurately locating an object in a reference frame.
Further, the described embodiments are not limited to the above-mentioned example of a locating system with three emitter coils and three receiver coils, but more generally apply to any electromagnetic locating system comprising at least one emitter coil and at least one receiver coil. It should be noted that in the case of a system with a single receiver coil, no measurement of the mutual inductance coefficient is necessary.
More generally, it should be noted that although a system enabling to measure possible drifts of parameters RsBRi, LBRi and/or Mli,i′ has been described, it may be chosen, according to the needs of the application, to only monitor a single one of these parameters, for example, parameter RsBRi which is generally the most critical, or a subset of these parameters.
Further, although an embodiment where the operations carried out by processing device 121 are implemented in the digital domain has been described hereabove, the described embodiments are not limited to this specific case. As a variation, the operations performed by device 121 may be totally or partly implemented in the analog domain. As an example, the operations of synchronous demodulation and of filtering and/or of synchronous detection may be implemented in the analog domain, and the operations of position calculation based on the measured field values may be implemented in the digital domain.
Such alterations, modifications, and improvements are intended to be part of this disclosure, and are intended to be within the spirit and the scope of the present invention. Accordingly, the foregoing description is by way of example only and is not intended to be limiting. The present invention is limited only as defined in the following claims and the equivalents thereto.
Number | Date | Country | Kind |
---|---|---|---|
1556827 | Jul 2015 | FR | national |
Number | Name | Date | Kind |
---|---|---|---|
6172499 | Ashe | Jan 2001 | B1 |
7321228 | Govari | Jan 2008 | B2 |
7373271 | Schneider | May 2008 | B1 |
20030184285 | Anderson et al. | Oct 2003 | A1 |
20040207389 | Nieminen et al. | Oct 2004 | A1 |
20050165297 | Anderson et al. | Jul 2005 | A1 |
20050246122 | Jones, Jr. et al. | Nov 2005 | A1 |
20130296651 | Ito et al. | Nov 2013 | A1 |
20140002063 | Ashe | Jan 2014 | A1 |
Number | Date | Country |
---|---|---|
2003102497 | Dec 2003 | WO |
Entry |
---|
Preliminary Search Report in French Patent Application No. 1556827, dated May 11, 2016, 2 pages. |
Number | Date | Country | |
---|---|---|---|
20170023381 A1 | Jan 2017 | US |