System for tracking an object

Abstract: A method for tracking an object (20) using energy fields (H 1 ,H 2 ,H 3 ), in the presence of interference due to introduction of an article (40) responsive to the fields, in the vicinity of the object (20). One or more desired energy fields are produced in the vicinity of the object (20), and a characteristic of energy fields induced due to introduction of the article is determined (40). Signals are generated responsive to the energy fields at a plurality of locations of the object after introduction of the article. Spatial coordinates of the object are determined, responsive to the generated signals and the determined characteristic.

Description: FIELD OF THE INVENTION The present invention relates generally to apparatus for generating and detecting electromagnetic fields, and specifically to non-contact, electromagnetic devices for tracking the position and orientation of an object, and associated disclosed method, not forming a part of the invention. BACKGROUND OF THE INVENTION Non-contact electromagnetic tracking systems are well known in the art, with a wide range of applications. For example, U.S. patent 4,054,881 describes a tracking system using three coils to generate electromagnetic fields in the vicinity of an object being tracked. The fields generated by these three coils are distinguished from one another by open loop multiplexing of time, frequency or phase. The signal currents flowing in three orthogonal sensor coils are used to determine the object's position, based on an iterative method of computation. U.S. patent 5,391,199, filed July 20, 1993, describes a system for generating three-dimensional location information regarding a medical probe or catheter. A sensor coil is placed in the catheter and generates signals in response to externally applied magnetic fields. The magnetic fields are generated by three radiator coils, fixed to an external reference frame in known, mutually spaced locations. The amplitudes of the signals generated in response to each of the radiator coil fields are detected and used to compute the location of the sensor coil. Each radiator coil is preferably driven by driver circuitry to generate a field at a known frequency, distinct from that of other radiator coils, so that the signals generated by the sensor coil may be separated by frequency into components corresponding to the different radiator coils. PCT patent publication WO/96/05768, filed January 24, 1995, describes a system that generates six-dimensional position and orientation information regarding the tip of a catheter and discloses an object tracking system of the type set forth in the preamble of the accompanying claim 1. This system uses a plurality of sensor coils adjacent to a locatable site in the catheter, for example near its distal end, and a plurality of radiator coils fixed in an external reference frame. These coils generate signals in response to magnetic fields generated by the radiator coils, which signals allow for the computation of six location and orientation coordinates. As in the case of the '539 patent application described above, the radiator coils preferably operate simultaneously at different frequencies, for example at 1000, 2000 and 3000 Hz, respectively. The above tracking systems rely on separation of position-responsive signals into components, most typically frequency components, wherein each such component is assumed to correspond uniquely to a single radiator coil, in a known position, radiating a magnetic field having a regular, well-defined spatial distribution. In practice, however, when a metal or other magnetically-responsive article is brought into the vicinity of the catheter or other object being tracked, the magnetic fields generated in this vicinity by the radiator coils are distorted. For example, the radiator coil's magnetic field may generate eddy currents in such an article, and the eddy currents will then cause a parasitic magnetic field to be radiated. Such parasitic fields and other types of distortion can lead to errors in determining the position of the object being tracked. In WO 97/29679, there is shown an endoscope with sensors that generate signals responsive to a magnetic field in order to determine positions and/or orientation of the endoscope. The sensors are located remote from interference causing structures in the endoscope. SUMMARY OF THE INVENTION It is an object of some aspects of the present invention to provide an object tracking system capable of performing a disclosed method, not forming a part of the invention, for improving the accuracy of position determination by a non-contact object tracking system. In one aspect of the disclosed method, the position determination is corrected so as to account for disturbances of an energy field used in tracking the object due to metal or other field-responsive articles in the vicinity of the object being tracked. In another aspect of the disclosed method, such disturbances are detected, so as to alert a user of the system that the accuracy of position determination may be compromised. Preferably, the energy field comprises a magnetic field, which causes position-responsive electrical signals to be generated in one or more coils associated with the object being tracked, for example a catheter or other medical probe. The present invention relies on the fact that parasitic magnetic fields, generated by metal or other field-responsive articles that receive and re-radiate energy from a radiator coil magnetic field are typically at the same frequency as the radiator coil field, but are shifted in phase relative thereto. The phase shift and the amplitudes of the parasitic fields generally depend on properties of the article, including dielectric constant, magnetic permeability and geometrical shape. Furthermore, when the article introduced into the magnet field is substantially magnetizable, i.e., its permeability is significantly different from the permeability of air, the lines of magnetic field in the vicinity of the article are generally distorted. Distortion of this type does not, however, substantially affect the phase of the magnetic field. In preferred embodiments of the present invention, an object tracking system comprises one or more sensor coils adjacent to a locatable point on an object being tracked, and one or more radiator coils, which generate magnetic fields when driven by electrical currents at respective driving frequencies in a vicinity of the object. Preferably each of the radiator coils has its own frequency, which is different from the frequencies of all the other radiator coils. Alternatively, however, two or more of the radiator coils may share a common frequency, and may be time-multiplexed so that only one of the two or more coils is driven to generate a magnetic field at any given time. The sensor coils generate electrical signals responsive to the magnetic fields, which signals are received by signal processing circuitry and analyzed by a computer or other processor. When a metal or other field-responsive article is in the vicinity of the object, the signals typically include position signal components responsive to the magnetic fields generated by the radiator coils at their respective driving frequencies, and parasitic signal components responsive to parasitic magnetic fields generated due to the article. The computer processes the signals to identify the parasitic components, preferably using a phase-sensitive method, as described below, and uses the position signal components to determine the position of the object. In some such preferred embodiments of the present invention, the signals due to the magnetic field generated by each of the radiator coils are first detected in the absence of any articles that could cause parasitic signal components to be generated. Baseline phases of the signals at each of the radiator coil frequencies, which are substantially independent of the position of the object relative to the radiator coils, are then determined relative to the phase of the current driving the respective radiator coil. When an article that generates parasitic magnetic fields is introduced into the vicinity of the object, at least some signal components generated by the sensor coils in response to the parasitic fields will generally be shifted in phase relative to the pre-determined baseline phase. The phase shift of the signals due to the parasitic signal components is detected and may be used to indicate that the position of the object determined by the computer may be inaccurate due to the presence of the article in the object's vicinity. Additionally or alternatively, while a metal or other field-responsive article is being moved into the vicinity of the object being tracked, the object is preferably held still. Changes in the amplitude and phase of the signals received from the sensor coils are then known to be associated with distortions of the magnetic field due to the article. These changes are detected and used in assessing the effect of the distortions on the determination of the object's position coordinates. In particular, the detected changes may be used to measure the phase shifts of the parasitic fields, relative to the radiator coil fields. In some cases, these phase shifts will be generally constant, independent of the position and orientation of the field-responsive article, for example, if the article is substantially symmetrical. Similarly, if the article comprises highly-conductive material, having negligible resistivity, the phase shift of the parasitic fields will be close to 90°, generally independent of the position and orientation of the article. Although not all articles will have such generally constant phase shifts, the dependence of the phase shift on position and orientation can be determined empirically in each case. In other preferred aspects of the disclosed method, a characteristic of the article, such as the distortion of the magnetic field that the article induces or the effect of the article on the determination of the object's coordinates, is known, for example based on prior measurement. The known characteristic may be used to correct the determination of the coordinates, as described below. In some preferred aspects of the disclosed method, the signals generated by the sensor coils are detected using a harmonic detection method, more preferably with the detection frequency synchronized with the frequencies of the radiator coils. Signal components that are out of phase with the baseline phase are eliminated from the sampled signal, for example using phase-sensitive signal correlation or other methods known in the art, thereby substantially removing at least some of the parasitic signal components and improving the accuracy of position determination in the presence of the metal or other field-responsive article. Preferably, if the phase shifts of the parasitic fields are known, having been measured as described above, for example, then the known phase shifts may be used in analyzing the signals so as to remove at least a portion of the parasitic signal components therefrom. By way of example, a method for tracking an object using an energy field, in the presence of interference due to introduction of an article responsive to the field, in the vicinity of the object, may include: producing a first energy field in the vicinity of the object; determining a characteristic of a second energy field induced responsive to the first field, due to introduction of the article; receiving signals responsive to the first and second energy fields generated at a plurality of locations of the object after introduction of the article; and determining spatial coordinates of the object responsive to the generated signals and the determined characteristic. Determining the characteristic includes determining a phase shift of the signals responsive to the second energy field with respect to the signals responsive to the first energy field. Preferably, the characteristic is not determined from measurements made at the location. Determining spatial coordinates includes determining an amplitude of the signals generated at a location of the object; processing the amplitude to find a corrected amplitude, using the determined characteristic; and calculating the coordinates based on the corrected amplitude. Preferably, determining the characteristic of the induced energy fields includes receiving a first signal responsive to the energy field in the absence of the field-responsive article; introducing the article into the vicinity of the object; receiving a second signal responsive to the energy field in the presence of the article; and processing the first and second signals to determine the characteristic. Preferably, processing the first and second signals includes finding a first amplitude and a first phase of the first signal; finding a second amplitude and a second phase of the second signal; and determining a phase shift associated with the induced field, based on the first and second amplitudes and the first and second phases. Further by way of example, determining spatial coordinates may include measuring an uncorrected amplitude and an uncorrected phase of the signals received at one of the plurality of locations of the object; finding a corrected amplitude, using the phase shift associated with the induced field; and calculating the coordinates based on the corrected amplitude. The disclosed method, not forming a part of the invention, includes detecting a malfunction in a system for tracking the object by finding a variability in one or more of the first and second phases. There is further disclosed a method including detecting the second energy field and notifying a user that the article has been introduced. Additionally or alternatively, determining the spatial coordinates of the object includes estimating an error bound on the coordinates in the presence of the article. Typically, producing the energy fields includes producing magnetic fields, and receiving the signals includes receiving electrical signals which are generated responsive to the magnetic field. There is further provided, in accordance with a preferred embodiment of the present invention, an object tracking system, including: a radiator, configured to generate a first energy field in the vicinity of the object; a sensor, fixed to the object, configured to generate signals responsive to the first energy field; and signal processing circuitry, configured to receive the signals from the sensor and configured to determine position coordinates of the object responsive thereto. The signal processing circuitry is configured to determine a characteristic of a second energy field induced responsive to the first field due to introduction of an article responsive to the first field in the vicinity of the object, configured to receive signals responsive to the first and second energy fields generated at a plurality of locations of the object after introduction of the article, and configured to determine spatial coordinates of the object responsive to the generated signals and the determined characteristic. The signal processing circuitry is configured to determine the characteristic by being configured to determine a phase shift of signals responsive to the second energy field with respect to the signals responsive to the first energy field. The signal processing circuitry is configured to determine spatial coordinates by being configured to determine an amplitude of the signals generated at a location of the object, process the amplitude to find a corrected amplitude, using the determined characteristic, and calculate the coordinates based on the corrected amplitude. There is disclosed the circuitry detecting a malfunction in the system by finding a variability in the phase shift. There is further disclosed the circuitry notifying a user of the system that the article has been introduced, based on the detected interference. Additionally or alternatively, the circuitry determines an error bound on the coordinates, based on the detected interference. Typically, the energy field includes a magnetic field, and the signals include electrical currents generated responsive to the magnetic field. A BRIEF DESCRIPTION OF THE DRAWINGS The present invention will be more fully understood from the following detailed description of the preferred embodiments thereof, taken together with the drawing in which: Fig. 1 is a schematic representation of an object tracking system in accordance with a preferred embodiment of the present invention; Fig. 2 is a flow chart illustrating a disclosed method, not forming a part of the present invention, for tracking an object; and Fig. 3 is a vector diagram useful in understanding principles of operation of the present invention DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS Reference is now made to Fig. 1, which schematically illustrates a system 10 for tracking a probe 20, such as a catheter for medical use, in accordance with preferred embodiment of the present invention. As described in the above-referenced U.S. patent 5,391,199 and PCT patent publication WO/96/05768, system 10 comprises a plurality of radiator coils 22, 24 and 26. These coils generate respective magnetic fields H 1 , H 2 and H 3 , at respective frequencies ω 1 , ω 2 and ω 3 , in the vicinity of probe 20. The probe includes sensor coils 27, 28 and 29, which generate electrical current signals in response to the magnetic fields. These signals comprise components at frequencies ω 1 , ω 2 and ω 3 , whose respective amplitudes are dependent on the position and orientation of probe 20. System 10 further comprises driver circuitry 30, 32 and 33, coupled to each of the radiator coils, which drives coils 22, 24 and 26 at the respective driving frequencies ω 1 , ω 2 and ω 3 . The signals generated by sensor coils 27, 28 and 29 are preferably received and processed by signal processing circuitry 34 and then used by computer 36 to calculate position and orientation coordinates of probe 20. Fig. 1 shows three radiator coils 22, 24 and 26 and three sensor coils 27, 28 and 29 in a probe 20. It will be understood, however, that the present invention is equally applicable to tracking systems comprising one, two, four or more radiator coils and one, two or more sensor coils. For six-dimensional tracking of probe 20 (three degrees of translation and three degrees of rotation), however, system 10 preferably include a total of six coils, counting both the radiator and sensor coils. Furthermore, although in preferred embodiments of the present invention are described here with reference to frequency-multiplexing of the radiator coils, the principles of the invention may similarly be applied to coils that are time-multiplexed or driven using other methods known in the art to distinguish their respective magnetic fields one from another. The present invention may be used in tracking other types of objects, as well. In the absence of parasitic effects, the signals generated by sensor coils 27, 28 and 29 at frequency ω 1 are proportional to the amplitude of the time derivative of the projection of field H 1 at probe 20 along the respective axes of the sensor coils. The signals generated at frequencies ω 2 and ω 3 are similarly proportional to the projections of H 2 and H 3 . Parasitic effects that may arise due to mutual inductance among the radiator coils are preferably substantially eliminated, as disclosed, for example, in PCT patent application no. PCT/IL/00100, filed March 18, 1997. Since the direction and amplitude of the magnetic field due to any one of radiator coils 22, 24 and 26 can be calculated easily using methods known in the art, the sensor coil signals due to the respective radiator coil field may be directly related to the sensor coil's distance from and orientation relative to the radiator coil. It will also be appreciated that in the absence of parasitic magnetic fields, such as will be described below, the phases of the signals at frequencies ω 1 , ω 2 and ω 3 are substantially constant relative to the phases of the magnetic fields generated by radiator coils 22, 24 and 26, independent of the position and orientation of sensor coils 27, 28 and 29. As shown in Fig. 1, however, when a metal or magnetic field-responsive article, for example surgical tool 40, is introduced into the vicinity of probe 20, the article will generally receive energy from fields H 1 , H 2 and H 3 , and re-radiate parasitic magnetic fields, H ' / 1 , H ' / 2 and H ' / 3 , at frequencies ω 1 , ω 2 and ω 3 . Generally the phases of the parasitic fields will be shifted relative to the radiator coil fields by phase angles  1 ',  2 ' and  3 ', respectively. The phases and amplitudes of the parasitic fields generally depend on properties of tool 40, including its dielectric constant, magnetic permeability, geometrical shape and orientation relative to the radiator coils. Although the phases may of the parasitic fields may generally vary as a function of the tool's position and orientation, in certain cases the phases will be substantially constant, for example, when tool 40 has suitable symmetry with respect to the fields generated by coils 22, 24 and 26, or when the tool comprises material of negligible resistivity. The parasitic fields produced by tool 40 will cause corresponding parasitic signal components to be generated in sensor coils 27, 28 and 29, so that the total signal I(t) received from any one of the sensor coils, including both position and parasitic signal components may generally be expressed as: where A i ,  i , A i ' and  i ' are the amplitude and phase of the position signal component and the parasitic signal component, respectively, at frequency ω i . It will be observed in equation (1) that for each of the signal frequency components I i , the superposition of the parasitic signal component will cause a phase shift in the total detected signal, relative to the signal phase in the absence of metal tool 40, given by: In preferred embodiments of the present invention, signal processing circuitry 34 and computer 36 detect and record a baseline phase of the signals received from sensor coils 27, 28 and 29 in the absence of any metal or other interfering magnetic field-responsive objects in the vicinity of probe 20. Alternatively, undisturbed phases of the position signal components may have been determined in advance for system 10 or are known based on the operation of the system. When metal tool 40 is introduced into the vicinity of probe 20, the phase shift due to the parasitic components engendered thereby in the signals is detected, and computer 36 may present a message to an operator of the probe, indicating that the determination of the probe's position may be inaccurate. Preferably, the amplitudes A i ' of the parasitic signal components are measured, for example as described below, and are then used to estimate and inform the operator of an error bound on the position coordinates determined by computer 36. Generally, the baseline phase of the signals received from sensor coils 27, 28 and 29 will be substantially constant, in the absence of interfering objects such as tool 40. If circuitry 34 receives a signal having a variable phase, particularly a randomly variable phase, such a signal is typically indicative of a malfunction in system 10, for example, a loose connection between one of the sensor coils and the circuitry. In such a case, probe 20 should be removed and repaired. Fig. 2 illustrates schematically a method for detecting interference due to tool 40 or due to another magnetic field-responsive article, and determining the correct position of probe 20 in the presence of such interference, in accordance with preferred embodiments of the present invention. Preferably, normalized baseline signals are first acquired by computer 36 at each of frequencies ω i in the absence of tool 40 or any other interfering article. Subsequently, after the tool is introduced into the vicinity of the probe 20, these baseline signals are compared with signals received from sensor coils 27, 28 and 29, using methods of signal correlation known in the art, for example, in order to measure the phase shifts  i '. Preferably, the phase shifts are measured for different positions and orientations of the tool, in order to ascertain the extent to which the phase shifts  i ' vary with tool position and orientation. If the measured phase shifts are found to be generally constant (or are known to be so from previous measurement), they may be used to calculate corrected values of A i , for example, based on formula (1). The position of probe 20 is then found with sufficient accuracy even in the presence of parasitic magnetic fields due to tool 40. Preferably, in order to measure phase shifts  ' / i , probe 20 is held immobile while tool 40 is being moved into its vicinity. Any changes in the amplitudes and phases of the signals received from sensor coils 27, 28 and 29 are then known to be associated with distortions of the magnetic field due to the article. These changes are detected and used in measuring the effect of the distortions on the determination of the probe's position coordinates. Specifically, the detected changes may be used to measure the phase shifts  ' / i of the parasitic fields, relative to the radiator coil fields. At low RF frequencies, for example up to 10 kHz, if tool 40 is sufficiently symmetrical and/or comprises material of negligible resistivity, the phase shifts  ' / i of the re-radiated fields will not change significantly as the tool is moved. Therefore, one measurement of  ' / i , made while probe 20 is held in any desired position and tool 40 is in any other suitable position in a vicinity of the probe, may be sufficient for calculating values of A i that are substantially corrected for the effect of the re-radiated fields, at all other positions of interest of probe 20 and tool 40. In the above preferred embodiment of the present invention, it is assumed that the re-radiated fields are the dominant source of magnetic field distortion due to tool 40. Distortion of the lines of magnetic field due to the magnetic permeability of tool 40 gives rise to parasitic magnetic field components that are generally in phase with fields H 1 , H 2 and H 3 , and are therefore more difficult to quantify and correct for. Fig. 3 is a vector diagram illustrating the method described above for separating the position and parasitic signal components. Signal vector 50, having amplitude |I i | and phase  i total , represents the signal received from one of the sensor coils at frequency ω i . Vector 50 is the vector sum of position signal component vector 52 and parasitic signal component vector 54. Position signal component 52 has amplitude A i and known, substantially constant phase  i , which is arbitrarily assigned the value zero, without loss of generality. If the phase  i ' (for the cases described above, in which  i ' can be considered to be substantially constant) of parasitic signal component 54 is also known, for example having been measured as described above, then a triangle defined by vectors 50, 52 and 54 is completely determined, since all its angles ( i ,  i ' and  i total ) and one of its sides (|I i |) are known. Respective amplitudes A i and A i ' are thus uniquely determined, and the coordinates of probe 20 may be found accurately based on the values of A i determined from the trigonometric formula: The method illustrated by Fig. 3 assumes that the phases of the parasitic signal components are shifted relative to the position signal components. However, when tool 40 or another magnetic field-responsive article in the vicinity of probe 20 has magnetic permeability, µ, that is significantly different from that of air, the lines of magnetic field may be substantially distorted, without significant phase shift. In this case, changes in the amplitudes of the signals from sensor coils 27, 28 and 29 are preferably detected as tool 40 is being introduced into the vicinity of probe 20, and are used to estimate an error bound on the position coordinates calculated by computer 36 in the presence of the tool. Generally the error will be greatest when probe 20 is closely adjacent to tool 40, and decrease as the distance between them increases. It will be understood that while the above methods have been described with reference to a system comprising three radiator coils and three sensor coils, they may equally be applied to other electromagnetic object tracking system, using greater or fewer numbers of coils or antennae. It will further be understood that the principles of the present invention may similarly be applied to object tracking systems based on other types of RF electromagnetic fields, wherein a sensor associated with the object generates signals whose amplitude and phase are dependent on the amplitude and phase of the radiation fields. Moreover, it will be appreciated that the preferred embodiments described above are cited by way of example, and the full scope of the invention is limited only by the claims.

Claims: An object tracking system (10), comprising: a radiator (22; 24; 26), configured to generate a first energy field (H 1 ; H 2 ; H 3 ) in the vicinity of the object (20) a sensor (27; 28; 29), fixed to the object (20), configured to generate signals (A i , (Φ i ) responsive to the first energy field; and signal processing circuitry (34), configured to receive the signals from the sensor (27; 28; 29) and configured to determine position coordinates of the object (20) responsive thereto, the signal processing circuitry (34) configured to: determine a characteristic of a second energy field (H 1 ; H 2 ; H 3 ) induced responsive to the,first field (H 1 ; H 2 ; H 3 ) due to introduction of an article (40) responsive to the first field, in the vicinity of the object (20); receive signals (I(t)) responsive to the first and second energy fields generated at a plurality of locations of the object (20) after introduction of the article (40); and characterised by the signal processing circuitry being configured to determine spatial coordinates of the object (20) responsive to the generated signals and the determined characteristic, wherein the signal processing circuitry (34) is configured to determine the characteristic by being configured to determine a phase shift (A i , Φ i ) of signals responsive to the second energy field with respect to the signals responsive to the first energy field, and wherein the signal processing circuitry (34) is configured to determine spatial coordinates by being configured to: determine an amplitude (|I i |, of the signals generated at a location of the object (20); process the amplitude to find a corrected amplitude (A 1 ), using the determined characteristic (Φ 1 ') ; and calculate the coordinates based on the corrected amplitude. A device in accordance with claim 1, wherein the signal processing circuitry (34) is configured to determine the characteristic (Φ i ') of the induced energy field (H' 1 ; H' 2; H' 3 ) by being configured to: receive a first signal (A i , Φ i ) responsive to the energy field (H 1 ; H 2 ; H 3 ) in the absence of the field-responsive article (40); receive a second signal (|I i |, Φ i total ) responsive to the energy field where the article (40) has been introduced into the vicinity of the object (20); and process the first and second signals to determine the characteristic (Φ i '). A device in accordance with claim 2, wherein the signal processing circuitry (34) is configured to process the first and second signals by being configured to: find a first amplitude (A i ) and a first phase (Φ i ) of the first signal; find a second amplitude (|I i |) and a second phase (Φ i total ) of the second signal; and determine a phase shift (Φ i '), associated with the induced field (H' 1 ; H' 2 ; H' 3 ), based on the first and second amplitudes (A I , |I i |) and the first and second phases (Φ i . Φ i total ). A device in accordance with claim 3, wherein the signal processing circuitry (34) is configured to determine the phase shift (Φ i ') by being configured to find values of (Φ i ') and of a third amplitude (A i ') associated with the induced field (H' 1 ; H' 2 ; H' 3 ), substantially as given by: A device in accordance with claim 4, wherein the signal processing circuitry is configured to determine spatial coordinates by being configured to: measure an uncorrected amplitude (|I i |) and an uncorrected phase (Φ i ) of the signals received at one of the plurality of locations of the object (20); find a corrected amplitude (A i ), using the phase shift (Φ i '); and calculate the coordinates based on the corrected amplitude. A device in accordance with claim 5, wherein the corrected amplitude is substantially as given by: A device in accordance with claim 6, further configured to detect a malfunction in a system (10) for tracking the object (20) by being configured to find a variability in one or more of the first and second phases (Φ i , Φ i total ).

Inventors: GOVARI ASSAF, OSADCHY DANIEL, GOVARI, ASSAF, OSADCHY, DANIEL
Applicants: BIOSENSE WEBSTER INC, BIOSENSE WEBSTER, INC.
Application year: 1999
Main classification: A61B5
IPC codes: A61N1, A61B5, A61M25, A61B19, G01B7