Electromagnetic tracking system with reduced distortion using quadratic excitation
Identifieur interne : 000124 ( Pmc/Corpus ); précédent : 000123; suivant : 000125Electromagnetic tracking system with reduced distortion using quadratic excitation
Auteurs : Tomasz Bien ; Mengfei Li ; Zein Salah ; Georg RoseSource :
- International Journal of Computer Assisted Radiology and Surgery [ 1861-6410 ] ; 2013.
Abstract
Electromagnetic tracking systems, frequently used in minimally invasive surgery, are affected by conductive distorters. The influence of conductive distorters on electromagnetic tracking system accuracy can be reduced through magnetic field modifications. This approach was developed and tested.
The voltage induced directly by the emitting coil in the sensing coil without additional influence by the conductive distorter depends on the first derivative of the voltage on the emitting coil. The voltage which is induced indirectly by the emitting coil across the conductive distorter in the sensing coil, however, depends on the second derivative of the voltage on the emitting coil. The electromagnetic tracking system takes advantage of this difference by supplying the emitting coil with a quadratic excitation voltage. The method is adaptive relative to the amount of distortion cause by the conductive distorters. This approach is evaluated with an experimental setup of the electromagnetic tracking system.
In vitro testing showed that the maximal error decreased from 10.9 to 3.8 mm when the quadratic voltage was used to excite the emitting coil instead of the sinusoidal voltage. Furthermore, the root mean square error in the proximity of the aluminum disk used as a conductive distorter was reduced from 3.5 to 1.6 mm when the electromagnetic tracking system used the quadratic instead of sinusoidal excitation.
Electromagnetic tracking with quadratic excitation is immune to the effects of a conductive distorter, especially compared with sinusoidal excitation of the emitting coil. Quadratic excitation of electromagnetic tracking for computer-assisted surgery is promising for clinical applications.
Url:
DOI: 10.1007/s11548-013-0925-4
PubMed: 23918006
PubMed Central: 3955490
Links to Exploration step
PMC:3955490Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Electromagnetic tracking system with reduced distortion using quadratic excitation</title>
<author><name sortKey="Bien, Tomasz" sort="Bien, Tomasz" uniqKey="Bien T" first="Tomasz" last="Bien">Tomasz Bien</name>
<affiliation><nlm:aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Li, Mengfei" sort="Li, Mengfei" uniqKey="Li M" first="Mengfei" last="Li">Mengfei Li</name>
<affiliation><nlm:aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Salah, Zein" sort="Salah, Zein" uniqKey="Salah Z" first="Zein" last="Salah">Zein Salah</name>
<affiliation><nlm:aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Rose, Georg" sort="Rose, Georg" uniqKey="Rose G" first="Georg" last="Rose">Georg Rose</name>
<affiliation><nlm:aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">23918006</idno>
<idno type="pmc">3955490</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3955490</idno>
<idno type="RBID">PMC:3955490</idno>
<idno type="doi">10.1007/s11548-013-0925-4</idno>
<date when="2013">2013</date>
<idno type="wicri:Area/Pmc/Corpus">000124</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000124</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Electromagnetic tracking system with reduced distortion using quadratic excitation</title>
<author><name sortKey="Bien, Tomasz" sort="Bien, Tomasz" uniqKey="Bien T" first="Tomasz" last="Bien">Tomasz Bien</name>
<affiliation><nlm:aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Li, Mengfei" sort="Li, Mengfei" uniqKey="Li M" first="Mengfei" last="Li">Mengfei Li</name>
<affiliation><nlm:aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Salah, Zein" sort="Salah, Zein" uniqKey="Salah Z" first="Zein" last="Salah">Zein Salah</name>
<affiliation><nlm:aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Rose, Georg" sort="Rose, Georg" uniqKey="Rose G" first="Georg" last="Rose">Georg Rose</name>
<affiliation><nlm:aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">International Journal of Computer Assisted Radiology and Surgery</title>
<idno type="ISSN">1861-6410</idno>
<idno type="eISSN">1861-6429</idno>
<imprint><date when="2013">2013</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Purpose</title>
<p> Electromagnetic tracking systems, frequently used in minimally invasive surgery, are affected by conductive distorters. The influence of conductive distorters on electromagnetic tracking system accuracy can be reduced through magnetic field modifications. This approach was developed and tested.</p>
</sec>
<sec><title>Methods</title>
<p> The voltage induced directly by the emitting coil in the sensing coil without additional influence by the conductive distorter depends on the first derivative of the voltage on the emitting coil. The voltage which is induced indirectly by the emitting coil across the conductive distorter in the sensing coil, however, depends on the second derivative of the voltage on the emitting coil. The electromagnetic tracking system takes advantage of this difference by supplying the emitting coil with a quadratic excitation voltage. The method is adaptive relative to the amount of distortion cause by the conductive distorters. This approach is evaluated with an experimental setup of the electromagnetic tracking system.</p>
</sec>
<sec><title>Results</title>
<p> In vitro testing showed that the maximal error decreased from 10.9 to 3.8 mm when the quadratic voltage was used to excite the emitting coil instead of the sinusoidal voltage. Furthermore, the root mean square error in the proximity of the aluminum disk used as a conductive distorter was reduced from 3.5 to 1.6 mm when the electromagnetic tracking system used the quadratic instead of sinusoidal excitation.</p>
</sec>
<sec><title>Conclusions</title>
<p> Electromagnetic tracking with quadratic excitation is immune to the effects of a conductive distorter, especially compared with sinusoidal excitation of the emitting coil. Quadratic excitation of electromagnetic tracking for computer-assisted surgery is promising for clinical applications.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Birkfellner, W" uniqKey="Birkfellner W">W Birkfellner</name>
</author>
<author><name sortKey="Hummel, J" uniqKey="Hummel J">J Hummel</name>
</author>
<author><name sortKey="Wilson, E" uniqKey="Wilson E">E Wilson</name>
</author>
<author><name sortKey="Cleary, K" uniqKey="Cleary K">K Cleary</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Weiser, T" uniqKey="Weiser T">T Weiser</name>
</author>
<author><name sortKey="Hyman, K" uniqKey="Hyman K">K Hyman</name>
</author>
<author><name sortKey="Yun, J" uniqKey="Yun J">J Yun</name>
</author>
<author><name sortKey="Litle, V" uniqKey="Litle V">V Litle</name>
</author>
<author><name sortKey="Chin, C" uniqKey="Chin C">C Chin</name>
</author>
<author><name sortKey="Swanson, S" uniqKey="Swanson S">S Swanson</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Lascalza, S" uniqKey="Lascalza S">S LaScalza</name>
</author>
<author><name sortKey="Arico, J" uniqKey="Arico J">J Arico</name>
</author>
<author><name sortKey="Hughes, R" uniqKey="Hughes R">R Hughes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Birkfellner, W" uniqKey="Birkfellner W">W Birkfellner</name>
</author>
<author><name sortKey="Watzinger, F" uniqKey="Watzinger F">F Watzinger</name>
</author>
<author><name sortKey="Wanschitz, F" uniqKey="Wanschitz F">F Wanschitz</name>
</author>
<author><name sortKey="Ewers, R" uniqKey="Ewers R">R Ewers</name>
</author>
<author><name sortKey="Bergmann, H" uniqKey="Bergmann H">H Bergmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schuler, N" uniqKey="Schuler N">N Schuler</name>
</author>
<author><name sortKey="Bey, M" uniqKey="Bey M">M Bey</name>
</author>
<author><name sortKey="Shearn, J" uniqKey="Shearn J">J Shearn</name>
</author>
<author><name sortKey="Butler, D" uniqKey="Butler D">D Butler</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Poulin, F" uniqKey="Poulin F">F Poulin</name>
</author>
<author><name sortKey="Amiot, L" uniqKey="Amiot L">L Amiot</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hummel, J" uniqKey="Hummel J">J Hummel</name>
</author>
<author><name sortKey="Figl, M" uniqKey="Figl M">M Figl</name>
</author>
<author><name sortKey="Birkfellner, W" uniqKey="Birkfellner W">W Birkfellner</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Kindratenko, V" uniqKey="Kindratenko V">V Kindratenko</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Plotkin, A" uniqKey="Plotkin A">A Plotkin</name>
</author>
<author><name sortKey="Kucher, V" uniqKey="Kucher V">V Kucher</name>
</author>
<author><name sortKey="Horen, Y" uniqKey="Horen Y">Y Horen</name>
</author>
<author><name sortKey="Paperno, E" uniqKey="Paperno E">E Paperno</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fitzpatrick, Jm" uniqKey="Fitzpatrick J">JM Fitzpatrick</name>
</author>
<author><name sortKey="Hill, Dl" uniqKey="Hill D">DL Hill</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Int J Comput Assist Radiol Surg</journal-id>
<journal-id journal-id-type="iso-abbrev">Int J Comput Assist Radiol Surg</journal-id>
<journal-title-group><journal-title>International Journal of Computer Assisted Radiology and Surgery</journal-title>
</journal-title-group>
<issn pub-type="ppub">1861-6410</issn>
<issn pub-type="epub">1861-6429</issn>
<publisher><publisher-name>Springer Berlin Heidelberg</publisher-name>
<publisher-loc>Berlin/Heidelberg</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">23918006</article-id>
<article-id pub-id-type="pmc">3955490</article-id>
<article-id pub-id-type="publisher-id">925</article-id>
<article-id pub-id-type="doi">10.1007/s11548-013-0925-4</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Original Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Electromagnetic tracking system with reduced distortion using quadratic excitation</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" corresp="yes"><name><surname>Bien</surname>
<given-names>Tomasz</given-names>
</name>
<address><phone>+49-391-6720105</phone>
<fax>+49-391-6720051</fax>
<email>Tomasz.Bien@ovgu.de</email>
</address>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Li</surname>
<given-names>Mengfei</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Salah</surname>
<given-names>Zein</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Rose</surname>
<given-names>Georg</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<aff id="Aff1">Chair for Healthcare Telematics and Medical Engineering, Otto-von-Guericke University, Magdeburg, Germany</aff>
</contrib-group>
<pub-date pub-type="epub"><day>6</day>
<month>8</month>
<year>2013</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>6</day>
<month>8</month>
<year>2013</year>
</pub-date>
<pub-date pub-type="ppub"><year>2014</year>
</pub-date>
<volume>9</volume>
<issue>2</issue>
<fpage>323</fpage>
<lpage>332</lpage>
<history><date date-type="received"><day>10</day>
<month>1</month>
<year>2013</year>
</date>
<date date-type="accepted"><day>4</day>
<month>7</month>
<year>2013</year>
</date>
</history>
<permissions><copyright-statement>© The Author(s) 2013</copyright-statement>
<license license-type="OpenAccess"><license-p><bold>Open Access</bold>
This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.</license-p>
</license>
</permissions>
<abstract id="Abs1"><sec><title>Purpose</title>
<p> Electromagnetic tracking systems, frequently used in minimally invasive surgery, are affected by conductive distorters. The influence of conductive distorters on electromagnetic tracking system accuracy can be reduced through magnetic field modifications. This approach was developed and tested.</p>
</sec>
<sec><title>Methods</title>
<p> The voltage induced directly by the emitting coil in the sensing coil without additional influence by the conductive distorter depends on the first derivative of the voltage on the emitting coil. The voltage which is induced indirectly by the emitting coil across the conductive distorter in the sensing coil, however, depends on the second derivative of the voltage on the emitting coil. The electromagnetic tracking system takes advantage of this difference by supplying the emitting coil with a quadratic excitation voltage. The method is adaptive relative to the amount of distortion cause by the conductive distorters. This approach is evaluated with an experimental setup of the electromagnetic tracking system.</p>
</sec>
<sec><title>Results</title>
<p> In vitro testing showed that the maximal error decreased from 10.9 to 3.8 mm when the quadratic voltage was used to excite the emitting coil instead of the sinusoidal voltage. Furthermore, the root mean square error in the proximity of the aluminum disk used as a conductive distorter was reduced from 3.5 to 1.6 mm when the electromagnetic tracking system used the quadratic instead of sinusoidal excitation.</p>
</sec>
<sec><title>Conclusions</title>
<p> Electromagnetic tracking with quadratic excitation is immune to the effects of a conductive distorter, especially compared with sinusoidal excitation of the emitting coil. Quadratic excitation of electromagnetic tracking for computer-assisted surgery is promising for clinical applications.</p>
</sec>
</abstract>
<kwd-group xml:lang="en"><title>Keywords</title>
<kwd>Electromagnetic tracking</kwd>
<kwd> Computer-assisted surgery</kwd>
<kwd>Minimally invasive surgery</kwd>
<kwd>Compensation of conductive distortion</kwd>
</kwd-group>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© CARS 2014</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1" sec-type="intro"><title>Introduction</title>
<p>An electromagnetic tracking system is a navigation system that is based on the measurement of electromagnetic fields with a known distribution. Based on those measurements, the position and orientation of the field sensor can be calculated. Electromagnetic tracking systems are not dependent on line-of-sight or rigid connections like in optical and mechanical tracking systems. It is used in computer-assisted surgery to track the position of the instrument relative to the patient’s anatomy [<xref ref-type="bibr" rid="CR1">1</xref>
]. It has been used successfully for navigation in ear, nose and throat surgery (ENT) as well as in cranial and spine applications. The availability of microsensors has also allowed for the use of EMTS in tightly controlled cardiac mapping and pulmonary applications. EMTS has also shown its benefit in laparoscopic ultrasonography, needle aspiration biopsy, ultrasound freehand tracking, brachytherapy, ultrasound bone registration, and has also been introduced in endoscopic and bronchoscopic procedures [<xref ref-type="bibr" rid="CR1">1</xref>
–<xref ref-type="bibr" rid="CR7">7</xref>
].</p>
<p>EMTS is often applied in environments that typically contain metallic objects (conductive, para- or ferromagnetic) or electronic devices that may affect the magnetic reference fields. Ferromagnetic materials alter the shape of the reference fields and consequently distort both the direct current (DC) and the alternating current (AC) EMTS [<xref ref-type="bibr" rid="CR8">8</xref>
–<xref ref-type="bibr" rid="CR10">10</xref>
]. The eddy currents induced by conductive materials generate secondary AC magnetic fields, which distort the reference field pattern [<xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR11">11</xref>
–<xref ref-type="bibr" rid="CR13">13</xref>
]. The resulting fields are difficult to determine or characterize analytically [<xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR14">14</xref>
]. Stray fields produced by electronic devices [<xref ref-type="bibr" rid="CR4">4</xref>
, <xref ref-type="bibr" rid="CR14">14</xref>
] or power lines [<xref ref-type="bibr" rid="CR9">9</xref>
] also distort the reference field. Internal electronic components of EMTS may also have an influence on the accuracy of position measurements [<xref ref-type="bibr" rid="CR14">14</xref>
]. Distortions originating from electronic devices and eddy currents additively overlay the reference field. Errors caused by nearby metallic objects are the most serious problem in clinical applications [<xref ref-type="bibr" rid="CR14">14</xref>
]. In a minimally invasive surgery operating table, LCD displays, light, metal studded walls, C-arms, operating microscopes and surgical instruments placed in proximity to the field generator affect the accuracy of EMTS [<xref ref-type="bibr" rid="CR9">9</xref>
, <xref ref-type="bibr" rid="CR14">14</xref>
–<xref ref-type="bibr" rid="CR16">16</xref>
].</p>
<p>Distortions in EMTS could be eliminated by placing potential sources of distortions far away from the region of interest (ROI). This is typically not possible in clinical applications without changing the medical workflow. There are many methods presented in literature for the compensation of distortions caused by metal devices fixed in position and orientation. A survey of such methods of calibration has been published [<xref ref-type="bibr" rid="CR17">17</xref>
].</p>
<p>Calibration techniques, however, cannot compensate dynamic distortions originating from mobile objects such as surgical tools or movable imaging equipment. In such circumstances, measurements could be averaged or filtered (e.g., with a Kalman filter) [<xref ref-type="bibr" rid="CR18">18</xref>
]. Another method utilizing reverse tracking has been claimed, whereby a magnetic field source is integrated within the medical device, i.e., the catheter to be tracked contains the magnetic field source while a sensor coil array is distributed outside the patient body [<xref ref-type="bibr" rid="CR19">19</xref>
]. This approach is less sensitive to external distortions by conductive objects but suffers from very weak magnetic fields which could be generated within a small volume inside the patient’s body. Another approach is the shielding of at least one side of the tracking system [<xref ref-type="bibr" rid="CR20">20</xref>
]. Aside from strategies to avoid distortions, there are also approaches for the detection and correction of the errors caused by conductive objects.</p>
<p>The detection of field distortion can be accomplished using two field sensors attached rigidly to the surgical tool. The known distance between the sensors is compared with the measured value. If the deviations exceed a given tolerance, an error message may be presented [<xref ref-type="bibr" rid="CR9">9</xref>
].</p>
<p>The correction of metallic disturbances can be based on redundant measurements [<xref ref-type="bibr" rid="CR21">21</xref>
]. In this manner, another approach measured the position using two frequencies [<xref ref-type="bibr" rid="CR22">22</xref>
]. The measurement with lower frequency was interpreted as less distorted and was used to correct the measurement with higher frequency. In order to decrease its influence on the accuracy of the tracking system, the medical instrument in use could be made from laminated conductive materials [<xref ref-type="bibr" rid="CR23">23</xref>
] similar to transformer cores. In the patent [<xref ref-type="bibr" rid="CR24">24</xref>
], the authors model distortion with RL circuits and claim to correct the measured voltage with the calculated parameter of the conductive distorter. The author of the patent [<xref ref-type="bibr" rid="CR25">25</xref>
] proposes to represent the conductive distorter as lag or lead–lag networks. The excitation of the emitting coil with at least two frequencies leads to the calculation of the parameters of this lead–lag network. The influence of the conductive distorters is removed by adjusting the magnetic fields based on the estimated parameters.</p>
<p>In contrast, the method of distortion compensation presented in this paper does not require modeling of the parameter of the conductive distorter. In proposed approach, the emitting coil is supplied with a quadratic excitation, which therefore allows the voltage induced in the sensing coil by the emitting coil to be distinguished from the voltage induced by the conductive distorter. Only the undistorted voltage induced in the sensing coil directly by the emitting coil is used for the position and orientation calculation.</p>
</sec>
<sec id="Sec2" sec-type="methods"><title>Methods</title>
<sec id="Sec3"><title>Equivalent circuit</title>
<p>The presented method of the distortion minimization is based on a model of the tracking system and on a model of a conductive distorter. AC EMTS typically uses a sensing coil as a field sensor, which measures the magnetic field of the emitting coil for further position estimation. The equivalent circuit of the magnetically coupled emitting and sensing coil is presented in Fig. <xref rid="Fig1" ref-type="fig">1</xref>
. In this equivalent circuit, the emitting coil is supplied with a voltage <inline-formula id="IEq1"><alternatives><tex-math id="M1">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${U}_{1}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq1.gif"></inline-graphic>
</alternatives>
</inline-formula>
and the voltage <inline-formula id="IEq2"><alternatives><tex-math id="M2">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${U}_{2}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq2.gif"></inline-graphic>
</alternatives>
</inline-formula>
is measured on the sensing coil. The amplitude of the voltage on the sensing coil is used for the position estimation. <inline-formula id="IEq3"><alternatives><tex-math id="M3">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${R}_{L1}, {R}_{L2}, {L}_{1}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq3.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq4"><alternatives><tex-math id="M4">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${L}_{2}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq4.gif"></inline-graphic>
</alternatives>
</inline-formula>
are resistances and inductances of the emitting coil and of the sensing coil, respectively. The currents <inline-formula id="IEq5"><alternatives><tex-math id="M5">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${I}_{1}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq5.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq6"><alternatives><tex-math id="M6">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${I}_{2}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq6.gif"></inline-graphic>
</alternatives>
</inline-formula>
are flowing in the emitting and in the sensing coil, respectively. <inline-formula id="IEq7"><alternatives><tex-math id="M7">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M_{12}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq7.gif"></inline-graphic>
</alternatives>
</inline-formula>
is the mutual inductance between the emitting and the sensing coil.<fig id="Fig1"><label>Fig. 1</label>
<caption><p>Equivalent circuit of the coupling between the emitting and the sensing coil</p>
</caption>
<graphic xlink:href="11548_2013_925_Fig1_HTML" id="MO1"></graphic>
</fig>
</p>
<p>The output of the model may be considered as an open loop due to the high impedance of the AD card. The current flow in the sensing coil <inline-formula id="IEq8"><alternatives><tex-math id="M8">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${I}_{2}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq8.gif"></inline-graphic>
</alternatives>
</inline-formula>
is therefore approximately equal to zero. Such a system may be also considered as an air transformer with a very low efficiency.</p>
<p>The resulting transmittance <inline-formula id="IEq9"><alternatives><tex-math id="M9">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{1{-}2}(s)$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq9.gif"></inline-graphic>
</alternatives>
</inline-formula>
linking the input voltage <inline-formula id="IEq10"><alternatives><tex-math id="M10">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${U}_{1}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq10.gif"></inline-graphic>
</alternatives>
</inline-formula>
and the output voltage <inline-formula id="IEq11"><alternatives><tex-math id="M11">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${U}_{2}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq11.gif"></inline-graphic>
</alternatives>
</inline-formula>
is given in the Laplace domain in Eq. (<xref rid="Equ1" ref-type="">1</xref>
), where <inline-formula id="IEq12"><alternatives><tex-math id="M12">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\tau }_{1}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq12.gif"></inline-graphic>
</alternatives>
</inline-formula>
is a time constant depending on the inductance and on the resistance of the emitting coil (<xref rid="Equ3" ref-type="">3</xref>
).<disp-formula id="Equ1"><label>1</label>
<alternatives><tex-math id="M13">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&{G}_{1{-}2}({s})=\frac{{U}_2 ({s})}{{U}_1 ({s})}=\frac{M_{12}^{\prime } \cdot {s}}{\left( \frac{1}{{\tau }_{1}}+{s}\right) },\end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ1.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ2"><label>2</label>
<alternatives><tex-math id="M14">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&M_{12}^{\prime }=\frac{M_{12}}{{L}_{1}},\end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ2.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ3"><label>3</label>
<alternatives><tex-math id="M15">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&{\tau }_1 =\frac{{L}_{1}}{{R}_{L1}}, \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ3.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
According to the transmittance, the voltage on the sensing coil is proportional to the derivative of the voltage on the emitting coil. The value of the mutual inductance <inline-formula id="IEq13"><alternatives><tex-math id="M16">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M_{12}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq13.gif"></inline-graphic>
</alternatives>
</inline-formula>
is used for the calculation of the magnetic field of the emitting coil and for solving the position and orientation algorithm.</p>
<p>Conductive distorters in the EMTS may be modeled as an additional RL circuit [<xref ref-type="bibr" rid="CR24">24</xref>
, <xref ref-type="bibr" rid="CR25">25</xref>
]. The equivalent circuit of two magnetically coupled coils with an additional conductive distorter is presented in Fig. <xref rid="Fig2" ref-type="fig">2</xref>
. Along with the coefficients common with the previous model, Fig. <xref rid="Fig2" ref-type="fig">2</xref>
comprises the current flow within the conductive distorter <inline-formula id="IEq14"><alternatives><tex-math id="M17">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${I}_{\mathrm{D}}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq14.gif"></inline-graphic>
</alternatives>
</inline-formula>
, inductance <inline-formula id="IEq15"><alternatives><tex-math id="M18">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${L}_{\mathrm{D}}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq15.gif"></inline-graphic>
</alternatives>
</inline-formula>
and resistance <inline-formula id="IEq16"><alternatives><tex-math id="M19">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${R}_{\mathrm{D}}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq16.gif"></inline-graphic>
</alternatives>
</inline-formula>
of the conductive distorter. The model also includes the mutual inductance between the emitting coil and the distorter <inline-formula id="IEq17"><alternatives><tex-math id="M20">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M_{1{\mathrm{D}}}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq17.gif"></inline-graphic>
</alternatives>
</inline-formula>
and the mutual inductance between the distorter and the sensing coil <inline-formula id="IEq18"><alternatives><tex-math id="M21">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M_{2{\mathrm{D}}}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq18.gif"></inline-graphic>
</alternatives>
</inline-formula>
.<fig id="Fig2"><label>Fig. 2</label>
<caption><p>Equivalent circuit of two magnetically coupled coils with conductive distorter modeled as RL circuit</p>
</caption>
<graphic xlink:href="11548_2013_925_Fig2_HTML" id="MO5"></graphic>
</fig>
</p>
<p>The resulting transmittance <inline-formula id="IEq19"><alternatives><tex-math id="M22">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_{1-D-2}(s)$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq19.gif"></inline-graphic>
</alternatives>
</inline-formula>
linking the input voltage <inline-formula id="IEq20"><alternatives><tex-math id="M23">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${U}_{1}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq20.gif"></inline-graphic>
</alternatives>
</inline-formula>
and the output voltage <inline-formula id="IEq21"><alternatives><tex-math id="M24">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${U}_{2}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq21.gif"></inline-graphic>
</alternatives>
</inline-formula>
of the two coils coupled indirectly across the RL circuit is presented in the Laplace domain in Eq. (<xref rid="Equ4" ref-type="">4</xref>
). The value of the time constant <inline-formula id="IEq22"><alternatives><tex-math id="M25">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\tau }_\mathrm{D}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq22.gif"></inline-graphic>
</alternatives>
</inline-formula>
of the conductive distorter (<xref rid="Equ6" ref-type="">6</xref>
) depends on the size, conductivity, form and the orientation of the conductive object relative to the magnetic field of the emitting coil. The voltage on the sensing coil is proportional to the second derivative of the voltage on the emitting coil.<disp-formula id="Equ4"><label>4</label>
<alternatives><tex-math id="M26">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&{G}_{1-D-2} ({s})=\frac{{U}_2 (s)}{{U}_1 (s)}=\frac{\frac{M_{1\mathrm{D}}}{{L}_1}\cdot {s}}{\left( {\frac{1}{{\tau }_1 }+{s}} \right) }\cdot \frac{\frac{M_{\mathrm{D}2}}{{L}_\mathrm{D}}\cdot {s}}{\left( {\frac{1}{{\tau }_\mathrm{D}}+{s}}\right) }\nonumber \\&\qquad \qquad \qquad \quad =\frac{M_{1\mathrm{D}2}^{\prime } \cdot {s}^{2}}{\left( {\frac{1}{{\tau }_1}+{s}} \right) \cdot \left( \frac{1}{{\tau }_{\mathrm{D}} +{s}} \right) },\end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ4.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ5"><label>5</label>
<alternatives><tex-math id="M27">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&M_{1\mathrm{D}2}^{\prime } =\frac{M_{1\mathrm{D}}}{{L}_{1}}\cdot \frac{M_{1\mathrm{D}}}{{L}_\mathrm{D} },\end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ5.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ6"><label>6</label>
<alternatives><tex-math id="M28">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&{\tau }_\mathrm{D} =\frac{{L}_\mathrm{D}}{{R}_\mathrm{D} }, \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ6.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
The transmittance <inline-formula id="IEq23"><alternatives><tex-math id="M29">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${G}_\mathrm{E}(s)$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq23.gif"></inline-graphic>
</alternatives>
</inline-formula>
linking the input voltage <inline-formula id="IEq24"><alternatives><tex-math id="M30">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${U}_{1}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq24.gif"></inline-graphic>
</alternatives>
</inline-formula>
and the output voltage <inline-formula id="IEq25"><alternatives><tex-math id="M31">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${U}_{2}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq25.gif"></inline-graphic>
</alternatives>
</inline-formula>
including both the directly and indirectly induced voltages (<xref rid="Equ7" ref-type="">7</xref>
) is a sum of both transmittances (<xref rid="Equ1" ref-type="">1</xref>
) and (<xref rid="Equ4" ref-type="">4</xref>
).<disp-formula id="Equ7"><label>7</label>
<alternatives><tex-math id="M32">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {G}_\mathrm{E} (s)=\frac{{U}_2 (s)}{{U}_1 (s)}=\frac{M_{12}^{\prime } \cdot {s}}{\left( {\frac{1}{{\tau }_1}+{s}} \right) }+\frac{M_{1\mathrm{D}2}^{\prime } \cdot {s}^{2}}{\left( {\frac{1}{{\tau }_1 }+{s}} \right) \cdot \left( \frac{1}{{\tau }_\mathrm{D}}+{s} \right) }, \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ7.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
The proposed model of the conductive distorter exhibits properties of a high-pass filter and therefore contributes additional filtered voltage to the voltage induced directly by the emitting coil. The transfer function is similar to the transfer function proposed by Schneider [<xref ref-type="bibr" rid="CR25">25</xref>
]. However, the model presented in this work links the voltage on the emitting coil with the voltage on the sensing coil. This difference is crucial for understanding the approach proposed in the following section.</p>
<p>The response of the system to a sinusoidal excitation without distortions is presented in the Laplace domain (<xref rid="Equ8" ref-type="">8</xref>
) and in the time domain (<xref rid="Equ9" ref-type="">9</xref>
), respectively.<disp-formula id="Equ8"><label>8</label>
<alternatives><tex-math id="M33">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_{\mathrm{WS}} (s)&= {U}_1 (s)\cdot {G}(s)=\frac{{U}\cdot {\omega } }{\left( {{\varpi } ^{2}+{s}^{2}} \right) }\cdot \frac{{s}\cdot M_{12}^{\prime } }{\left( {\frac{1}{{\tau }_1 }+{s}} \right) }\nonumber \\&= \frac{{U}\cdot {\omega } \cdot {s}\cdot M_{12}^{\prime } }{\left( {{\varpi }^{2}+{s}^{2}} \right) \cdot \left( {\frac{1}{{\tau }_1 }+{s}} \right) },\end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ8.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ9"><label>9</label>
<alternatives><tex-math id="M34">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_{\mathrm{WS}} ({t})&= {U}\cdot M_{12}^{\prime } \nonumber \\&\quad \cdot \frac{{\omega } \cdot {\tau } _{1}\cdot \left( {\cos \left( {{\varpi } {t}} \right) \!-\!\mathrm{e}^{-\frac{{t}}{{\tau }_1 }}\!+\!{\omega } \cdot {\tau }_1 \cdot \sin ({\varpi } {t})} \right) }{\left( {{\varpi }^{2}\cdot {\tau }_1^2 +1} \right) },\nonumber \\ \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ9.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
The response may be distinguished between a steady-state response comprising a sum of sinusoidal functions and a transient response comprising an exponential function. After a time period depending on the time constant <inline-formula id="IEq26"><alternatives><tex-math id="M35">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${\tau }_{1}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq26.gif"></inline-graphic>
</alternatives>
</inline-formula>
, the transient response decays to zero and the response comprises only the sum of sinusoidal voltages.</p>
<p>The response of the system to a sinusoidal excitation caused by distortions only is presented in the Laplace domain (<xref rid="Equ10" ref-type="">10</xref>
) and in the time domain (<xref rid="Equ11" ref-type="">11</xref>
), respectively, where U is the amplitude of the sinusoidal excitation.<disp-formula id="Equ10"><label>10</label>
<alternatives><tex-math id="M36">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_{\mathrm{DS}} (s)&= {U}_1 (s)\cdot {G}_{1{-}2} (s)=\frac{{U}\cdot {\omega }}{\left( {{\varpi }^{2}+{s}^{2}} \right) }\nonumber \\&\quad \cdot \frac{{s}^{2}\cdot M_{1\mathrm{D}2}^{\prime } }{\left( {\frac{1}{{\tau }_1}+{s}} \right) \cdot \left( {\frac{1}{{\tau }_\mathrm{D}}+{s}} \right) }\nonumber \\&= \frac{{U}\cdot {\omega } \cdot {s}^{2}\cdot M_{1\mathrm{D}2}^{\prime } }{\left( {{\varpi } ^{2}+s^{2}} \right) \cdot \left( {\frac{1}{{\tau }_1}+{s}} \right) \cdot \left( {\frac{1}{{\tau }_\mathrm{D}}+{s}}\right) }, \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ10.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ11"><label>11</label>
<alternatives><tex-math id="M37">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&{U}_{\mathrm{DS}} ({t})={U}\cdot M_{1\mathrm{D}2}^{\prime } \cdot \frac{{\omega } \cdot {\tau }_1\cdot {\tau }_\mathrm{D} \cdot \mathrm{e}^{-\frac{{t}}{{\tau }_1}}}{\left( {{\varpi }^{2}\cdot {\tau } _1^2 \!+\!1}\right) \cdot \left( {{\tau }_1 \!-\!{\tau }_\mathrm{D}}\right) }\nonumber \\&\!\quad -{U}\cdot M_{1\mathrm{D}2}^{\prime } \cdot \frac{{\omega } \cdot {\tau }_{1} \cdot {\tau }_\mathrm{D} \cdot \mathrm{e}^{-\frac{{t}}{{\tau }_\mathrm{D}}}}{\left( {{\varpi } ^{2}\cdot {\tau }_\mathrm{D}^2 +1} \right) \cdot \left( {{\tau }_1 -{\tau }_\mathrm{D}} \right) }\nonumber \\&\!\quad + {U}\cdot M_{1\mathrm{D}2}^{\prime }\cdot \frac{{\omega }^{2}\cdot {\tau } _1 \cdot {\tau }_\mathrm{D} \cdot \left( {\omega } \cdot {\tau }_{1} \cdot \cos \left( {\varpi } {t} \right) -\sin ({\varpi } {t})\right) }{\left( {\varpi }^{2}\cdot {\tau }_{1}^{2} +1\right) \cdot \left( {\varpi }^{2}\cdot {\tau }_\mathrm{D}^2 +1\right) }\nonumber \\&\!\quad +{U} \cdot M_{1\mathrm{D}2}^{\prime } \cdot \frac{{\omega }^{2} \cdot {\tau }_{1} \cdot {\tau }_\mathrm{D} \cdot \left( {\omega } \cdot {\tau }_\mathrm{D} \cdot \cos \left( {\varpi } {t} \right) \!+\!\omega ^{2}\cdot {\tau }_{1} \cdot {\tau }_\mathrm{D} \cdot \sin \left( {\varpi } {t} \right) \right) }{\left( {\varpi }^{2}\cdot {\tau }_{1}^{2} \!+\!1\right) \cdot \left( {\varpi }^{2} \cdot {\tau }_\mathrm{D}^2 \!+\!1\right) },\nonumber \\ \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ11.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
The sum of the voltage induced directly (<xref rid="Equ9" ref-type="">9</xref>
) and indirectly (<xref rid="Equ11" ref-type="">11</xref>
) in the sensing coil builds a resulting voltage. The voltage induced by the conductive distorter has the same frequency as the voltage induced directly by the emitting coil, and therefore, the distortion due to the conductive distorter cannot be filtered out of the measured sum of signals. The only feature which distinguishes the voltage induced in the sensing coil by the emitting coil from the voltage induced by the conductive distorter is a slightly different phase shift.</p>
</sec>
<sec id="Sec4"><title>Quadratic excitation</title>
<p>As a solution to the problem of distinguishing the voltage induced in the sensing coil directly by the emitting coil and indirectly by the conductive distorter, an EMTS with a quadratic excitation on the emitting coil (<xref rid="Equ12" ref-type="">12</xref>
) is proposed.<disp-formula id="Equ12"><label>12</label>
<alternatives><tex-math id="M38">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_\mathrm{Q} ({t})={U}\cdot {t}^{2}, \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ12.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
The response of the sensing coil on the quadratic excitation of the emitting coil without conductive distorter is presented in the Laplace domain (<xref rid="Equ13" ref-type="">13</xref>
) and in the time domain (<xref rid="Equ14" ref-type="">14</xref>
), respectively.<disp-formula id="Equ13"><label>13</label>
<alternatives><tex-math id="M39">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_{\mathrm{WQ}} (s)&= X(s)\cdot {G}_{1{-}2} (s)=\frac{2\cdot {U}}{{s}^{3}}\cdot \frac{{s}\cdot M_{12}^{\prime } }{\left( {\frac{1}{{\tau } _{1}}+{s}}\right) }\nonumber \\&= \frac{2\cdot {U}\cdot {s}\cdot M_{12}^{\prime } }{{s}^{3}\cdot \left( {\frac{1}{{\tau }_1 }+{s}} \right) },\end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ13.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ14"><label>14</label>
<alternatives><tex-math id="M40">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_{\mathrm{WQ}} ({t})&= 2\cdot {U}\cdot {\tau }_{1} \cdot M_{12}^{\prime } \cdot \left( {t}-{\tau }_{1} +{\tau }_{1} \cdot \mathrm{e}^{-\frac{{t}}{{\tau }_{1}}} \right) , \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ14.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
The system response on the sensing coil to the quadratic excitation caused by the conductive distorter only is presented in the Laplace domain (<xref rid="Equ15" ref-type="">15</xref>
) and in the time domain (<xref rid="Equ16" ref-type="">16</xref>
), respectively.<disp-formula id="Equ15"><label>15</label>
<alternatives><tex-math id="M41">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_{\mathrm{DQ}} (s)&= X(s)\cdot {G}_{1{-}2} (s)\nonumber \\&= \frac{2\cdot {U}}{{s}^{3}}\cdot \frac{{s}^{2}\cdot M_{\mathrm{1D2}}^{\prime } }{\left( {\frac{1}{{\tau }_{1}}+{s}} \right) \cdot \left( \frac{1}{{\tau }_\mathrm{D}}+{s} \right) }\nonumber \\&= \frac{2\cdot {U}\cdot {s}^{2}\cdot M_{\mathrm{1D2}}^{\prime } }{{s}^{3}\cdot \left( {\frac{1}{{\tau }_{1}}+{s}} \right) \cdot \left( {\frac{1}{{\tau }_\mathrm{D}}+{s}}\right) }, \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ15.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ16"><label>16</label>
<alternatives><tex-math id="M42">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&{U}_{\mathrm{DQ}} ({t})=2\cdot {U}\cdot M_{\mathrm{1D2}}^{\prime }\nonumber \\&\quad \quad \cdot \left( \frac{1}{{\tau }_1 \cdot {\tau }_\mathrm{D}}\!+\!\frac{1}{{\tau }_1 \cdot \left( {{\tau }_{1} -{\tau }_\mathrm{D}}\right) }\cdot \mathrm{e}^{-\frac{{t}}{{\tau }_{1}}}-\frac{1}{{\tau }_{\mathrm{D}} \cdot \left( {\tau }_{1} -{\tau }_\mathrm{D} \right) }\cdot \mathrm{e}^{-\frac{{t}}{{\tau }_{\mathrm{D}}}} \right) ,\nonumber \\ \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ16.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
</p>
<p>The directly induced voltage (<xref rid="Equ14" ref-type="">14</xref>
) comprises a saw tooth, a rectangular pulse and an exponential transient response. The voltage induced indirectly (<xref rid="Equ16" ref-type="">16</xref>
) by the distorter comprises a rectangular pulse and an exponential transient response only. This difference is caused by a double derivative of the signal induced by the emitting coil induced in the sensing coil indirectly across a conductive distorter (<xref rid="Equ4" ref-type="">4</xref>
). According to the model, a saw tooth signal form of the voltage is induced in the conductive distorter. After the second derivative, the signal form is changed into a rectangular wave with an additional exponential transient response. The missing saw tooth signal form distinguishes the signal induced in the sensing coil directly by the emitting coil from the signal induced indirectly across a conductive distorter in the sensing coil. It also provides a measurement of the amplitude of the signal induced directly by the emitting coil while neglecting the voltage induced by the source of distortion. However, the measured signal also includes the exponential transient response, which overlays the steady-state response. After a time period dependent on the properties of the conductive distorter, the transient response decays to zero and the signal comprises the saw tooth and the rectangular pulse only. The measurement of the saw tooth wave allows calculation of the mutual inductance between the emitting and the sensing coil. This value of the mutual inductance allows estimation of the position of the sensor without an influence from the conductive distorter.</p>
<p>When the emitting coil is excited by quadratic excitation, the transient response in the form of the sum of exponential functions overlays the steady-state response comprising the saw tooth function. This problem could be solved by increasing the period of the quadratic excitation. However, the amplitude of the signal in the sensing coil depends on the time derivative of the current in the emitting coil. The increase in the length of the period would decrease the amplitude of the measured voltage in the sensing coil and therefore significantly reduce the working volume of the electromagnetic tracking system. In the presented paper, the transient response of the system is estimated in the first step using the response on the excitation by an additional rectangle pulse (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
c) supplied into the emitting coil shortly after supplying it with the quadratic excitation (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
a). The resulting response of the system (<xref rid="Equ17" ref-type="">17</xref>
) and of the conductive distorters (<xref rid="Equ18" ref-type="">18</xref>
) on a rectangle input voltage is as follows:<disp-formula id="Equ17"><label>17</label>
<alternatives><tex-math id="M43">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_{\mathrm{WR}} ({t})&= {U}\cdot M_{\mathrm{\mathrm{1D2}}}^{\prime } \cdot \mathrm{e}^{-\frac{{t}}{{\tau }_{1}}},\end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ17.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ18"><label>18</label>
<alternatives><tex-math id="M44">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} {U}_{\mathrm{DR}} ({t})&= {U}\cdot M_{\mathrm{1D2}}^{\prime } \nonumber \\&\quad \cdot \left( {\frac{{\tau }_{1}}{{\tau }_{1} -{\tau }_{\mathrm{D}}}\cdot \mathrm{e}^{-\frac{{t}}{{\tau }_1 }}-\frac{{\tau }_\mathrm{D} }{{\tau }_{1} -{\tau }_\mathrm{D}}\cdot \mathrm{e}^{-\frac{{t}}{{\tau }_\mathrm{D}}}} \right) , \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ18.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
This response to a rectangular excitation contains information about the tracking system and the conductive distorter and comprises the transient response only (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
d). The emitting coil will be supplied with a rectangular signal directly after supplying it with the quadratic signal form. Therefore, it is assumed that the position of the sensing coil and the position of the conductive distorter relative to the emitting coils are constant over the entire position and orientation estimation. The response of the system to a rectangular input signal leads to the calculation of the transient response of the system to a quadratic input signal. This calculation is based on a relationship between the rectangular and the quadratic signal, that is, the second time derivative of the quadratic signal produces a rectangular form of the signal.<fig id="Fig3"><label>Fig. 3</label>
<caption><p>Extraction of the steady-state response from the response on the quadratic excitation of the emitting coil measured by ADC: <bold>a</bold>
quadratic excitation fed into the emitting coil measured by ADC, <bold>b</bold>
response on the sensing coil to the quadratic excitation of the emitting coil measured by ADC, <bold>c</bold>
rectangular excitation fed into the emitting coil measured by ADC, <bold>d</bold>
response on the sensing coil to the rectangular excitation of the emitting coil measured by ADC, <bold>e</bold>
transient response of the system on the rectangular excitation estimated with double integration of the rectangular excitation of the emitting coil, <bold>f</bold>
response on the sensing coil to the quadratic excitation with subtracted transient response</p>
</caption>
<graphic xlink:href="11548_2013_925_Fig3_HTML" id="MO19"></graphic>
</fig>
</p>
<p>Therefore, the second numerical integration of the response on the rectangular wave with suitable amplitude is equal to the transient response of the system to the quadratic signal when the constant resulting from the numerical integration is removed. After the calculation, the estimation with the described method of the transient response of the system on the quadratic excitation (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
e) will be subtracted from the voltage measured on the sensing coil (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
f).</p>
<p>The transient response on the quadratic excitation is calculated with double integration of the response on the rectangle excitation:<disp-formula id="Equ19"><label>19</label>
<alternatives><tex-math id="M45">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\int \int {{U}_{\mathrm{WR}} ({t})} \mathrm{dt}^{2}={U}\cdot M_{12}^{\prime } \cdot {\tau } _{1}^{2} \cdot \mathrm{e}^{-\frac{{t}}{{\tau }_{1}}},\end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ19.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ20"><label>20</label>
<alternatives><tex-math id="M46">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&\int \int {{U}_{\mathrm{DR}} ({t})} \mathrm{dt}^{2}={U}\cdot M_{\mathrm{1D2}}^{\prime } \nonumber \\&\quad \cdot \left( \frac{1}{{\tau }_{1} \cdot \left( {\tau }_{1} -{\tau }_{\mathrm{D}} \right) }\cdot \mathrm{e}^{-\frac{{t}}{{\tau }_{1}}}-\frac{1}{{\tau }_{\mathrm{D}} \cdot \left( {\tau }_{1} -{\tau } _\mathrm{D}\right) }\cdot \mathrm{e}^{-\frac{{t}}{{\tau }_\mathrm{D}}}\! \right) \!, \end{aligned}$$\end{document}</tex-math>
<graphic xlink:href="11548_2013_925_Article_Equ20.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
The estimation of the position of the sensing coil requires calculation of the mutual inductance between every coupled emitting and sensing coil. The known dimensions of the emitting and sensing coil and the mutual inductance between the emitting and the sensing coil lead to the estimation of the position of the sensing coil. When the emitting coil is fed with the sinusoidal excitation, the mutual inductance is calculated based on the voltage measured in the sensing coil and the known amplitude of the voltage in the emitting coil. On the other hand, when the emitting coil is fed with the quadratic excitation, the mutual inductance between the emitting and the sensing coil will be calculated based on the gradient of the saw tooth signal induced in the sensing coil. Its curve is fitted with the least squares method. The gradient of the saw tooth signal leads to the calculation of the mutual inductance that consequently leads to the estimation of the position of the sensing coil.</p>
</sec>
</sec>
<sec id="Sec5"><title>Experimental evaluation</title>
<sec id="Sec6"><title>Experimental setup of the tracking system</title>
<p>For an experimental evaluation of the proposed method, an experimental setup of the tracking system was used (Fig. <xref rid="Fig4" ref-type="fig">4</xref>
). It includes 8 emitting coils each with a diameter of 40 mm supplied with a voltage generated by an analog-to-digital converter NI PXI 6281 (ADC) and amplified by an amplifier. The experimental setup utilizes the NDI field sensor with a diameter of 0.3 mm as a sensing coil. The voltage induced in the sensing coil is amplified by an additional amplifier and measured by ADC. The resulting position of the sensing coil is estimated in Matlab by solving the position and orientation algorithm [<xref ref-type="bibr" rid="CR26">26</xref>
].<fig id="Fig4"><label>Fig. 4</label>
<caption><p>Experimental setup of the electromagnetic tracking with attached optical marker for the evaluation of the EMTS with quadratic excitation</p>
</caption>
<graphic xlink:href="11548_2013_925_Fig4_HTML" id="MO24"></graphic>
</fig>
</p>
</sec>
<sec id="Sec7"><title>Measurement setup</title>
<p>In multiple interventions in computer-assisted surgery, the sensing coil of the EMTS is attached to surgical instruments which are made of conductive materials. In the experimental validation of the proposed method, an aluminum disk with a radius of 0.05 m and a thickness of 0.03 m has been utilized as such a conductive distorter. This disk represents the surgical tools used in minimal invasive surgery, e.g., hammer, trocar, endoscope or ultrasonic probe. Those instruments are usually made of less electrically conductive and non-ferromagnetic steel. Therefore, in this experiment, the more electrically conductive aluminum disk may be considered as a worst-case scenario of distortion.
</p>
<p>The goal of the performed experiment was to examine the influence of the conductive distorter on the accuracy of the tracking system with the quadratic excitation and compare it to the accuracy of the tracking system with the sinusoidal excitation. The sensing coil was moved within the working surface by a non-conductive plastic robotic arm in 200 random positions in proximity to the conductive distorter. The working surface had a rectangular shape (<inline-formula id="IEq27"><alternatives><tex-math id="M47">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$0.1\hbox {m}\times 0.1\hbox {m}$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq27.gif"></inline-graphic>
</alternatives>
</inline-formula>
), where the middle of the rectangle is placed in the middle of the coordinate system parallel to the surface of the field generator at a height of 0.15 m. The working surface was chosen to be relatively small because the experiment was designed in order to focus on the influence of the conductive distorter on the accuracy of the EMTS. The measurements were therefore performed in direct proximity to the conductive distorter in order to achieve the highest influence on the accuracy of the electromagnetic tracker. A larger distance between the sensing coil and emitting coils reduces the accuracy of the EMTS. Furthermore, the voltage on the sensing coil was measured for 8 s for the sinusoidal excitation (1 kHz) and 32 s for the quadratic excitation and averaged in order to increase the accuracy of the experimental setup of the tracking system.</p>
<p>After each movement of the robotic arm, the position of the sensing coil was measured by the electromagnetic tracking system and by the optical tracking system Polaris Spectra. Accuracy of this optical tracking system (root mean square error <inline-formula id="IEq28"><alternatives><tex-math id="M48">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$<$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq28.gif"></inline-graphic>
</alternatives>
</inline-formula>
0.25 mm) exceeds the expected accuracy of the EM tracking and therefore was used as a reference for the evaluation of the accuracy of the position estimation by the EMTS. The EMTS uses sinusoidal excitation in the first position measurement and shortly after, while in the same position, uses quadratic excitation. After 200 random measurements in the presence of the aluminum disk, the entire experiment was repeated without the presence of any conductive distorter.</p>
<p>In order to evaluate the electromagnetic tracking system, the coordinate system of the EMTS was registered with the coordinate system of the optical tracking system prior to the accuracy evaluation. For the performed registration, 27 additional measurements within a working volume of 0.2 m <inline-formula id="IEq29"><alternatives><tex-math id="M49">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\times $$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq29.gif"></inline-graphic>
</alternatives>
</inline-formula>
0.2 m <inline-formula id="IEq30"><alternatives><tex-math id="M50">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\times $$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq30.gif"></inline-graphic>
</alternatives>
</inline-formula>
0.2 m with both the electromagnetic and the optical tracking system were performed. The registration of both coordinate systems causes additional inaccuracy that increases the mean root square error of the position estimation with the optical tracking system. The registration error was calculated according to the work by Fitzpatrick [<xref ref-type="bibr" rid="CR27">27</xref>
]. The registration was carried out by fiducial registration via the software 3D Slicer.</p>
</sec>
</sec>
<sec id="Sec8"><title>Results and discussion</title>
<sec id="Sec9"><title>Results</title>
<p>The registration error introduced by registering the optical and electromagnetic tracking system is equal to 0.7 mm. The experimental errors in the position estimation described in the previous chapter are presented in Fig. <xref rid="Fig5" ref-type="fig">5</xref>
. The errors were calculated in four different scenarios. Figure <xref rid="Fig5" ref-type="fig">5</xref>
a, b presents the errors in the position estimation without the aluminum disk with the quadratic and sinusoidal excitation, respectively. Furthermore, Fig. <xref rid="Fig5" ref-type="fig">5</xref>
c, d presents the errors in the position estimation with the aluminum disk with the quadratic and sinusoidal excitation, respectively.<fig id="Fig5"><label>Fig. 5</label>
<caption><p>The influence of the aluminum disk on the accuracy of the electromagnetic tracking system as a function of the distance between aluminum disk and the sensing coil: <bold>a</bold>
tracking with quadratic excitation without the aluminum disk, <bold>b</bold>
tracking with sinusoidal excitation without the aluminum disk, <bold>c</bold>
tracking with quadratic excitation nearby to the aluminum disk, <bold>d</bold>
tracking with sinusoidal excitation nearby to the aluminum disk</p>
</caption>
<graphic xlink:href="11548_2013_925_Fig5_HTML" id="MO25"></graphic>
</fig>
</p>
<p>The maximal error, the root mean square error (RMSE) and the standard deviation (std) for four different scenarios within the performed experiment are presented in Table <xref rid="Tab1" ref-type="table">1</xref>
.<table-wrap id="Tab1"><label>Table 1</label>
<caption><p>Accuracy of the electromagnetic tracker with and without the aluminum disk in the proximity</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th align="left"></th>
<th align="left">Max error (mm)</th>
<th align="left">RMSE (mm)</th>
<th align="left">Std (mm)</th>
</tr>
</thead>
<tbody><tr><td align="left">Sinusoidal input voltage without conductive distorter</td>
<td align="center">1.5</td>
<td align="center">0.9</td>
<td align="left">0.3</td>
</tr>
<tr><td align="left">Quadratic input voltage without conductive distorter</td>
<td align="center">3.4</td>
<td align="center">1.7</td>
<td align="left">0.6</td>
</tr>
<tr><td align="left">Sinusoidal input voltage with conductive distorter</td>
<td align="center">10.9</td>
<td align="center">2.4</td>
<td align="left">1.7</td>
</tr>
<tr><td align="left">Quadratic input voltage with conductive distorter</td>
<td align="center">3.8</td>
<td align="center">2.0</td>
<td align="left">0.6</td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
</sec>
<sec id="Sec10"><title>Discussion</title>
<p>According to Fig. <xref rid="Fig5" ref-type="fig">5</xref>
, when the aluminum disk is placed in proximity of the sensing coil, the maximal error of the position estimation in the scenario with the sinusoidal excitation increased from 1.5 to 10.9 mm. The root mean square error of the position estimation increased from 0.9 to 2.4 mm. This shows a high influence of the conductive aluminum disk on the accuracy of the electromagnetic tracker. Despite the high increase in the maximal error of the position estimation in the sinusoidal excitation scenario with the conductive distorter in proximity to the sensing coil, the maximal error in the scenario with quadratic excitation changed only slightly from 3.4 to 3.8 mm. The RMSE with quadratic excitation changed from 1.7 mm without any distorters to 2.0 mm in the presence of conductive distorters. This increase in the RMSE in the presence of the distorter is smaller than the standard deviation of the measurements.</p>
<p>Additionally, Fig. <xref rid="Fig5" ref-type="fig">5</xref>
c does not indicate the existence of a trend in the error measurement as the distance between the surface of the aluminum disk and the sensing coil increased for the quadratic excitation scenario. A declining trend could be seen, however, in Fig. <xref rid="Fig5" ref-type="fig">5</xref>
d, where the sinusoidal excitation has been used. The results indicate very low influence of the conductive distorter on the accuracy of the electromagnetic tracker with quadratic excitation.</p>
<p>The root mean square error of the position estimation without the additional distorter for the EMTS with quadratic excitation (1.7 mm) is higher than with the sinusoidal excitation (0.9 mm). The lower accuracy in the position estimation is caused by the noise originating from, e.g., the amplifier, the discretization, the electromagnetic field generated by surrounding electronic devices or noise with a frequency of 50 Hz induced by the power supply. Due to the analysis of the system response on the quadratic excitation in the time domain, the response of the system was not filtered out. Only a low-pass filter of the ADC with a cutoff frequency equal to 80 kHz was used. In the proposed method, the emitting coils have been supplied with both the quadratic and rectangular signal form, which has a negative effect on the measuring rate.</p>
</sec>
</sec>
<sec id="Sec11" sec-type="conclusions"><title>Conclusion</title>
<p>In this paper, a novel method of electromagnetic tracking is proposed, which by using quadratic excitation for position estimation reduces the error in the position estimation in the presence of the conductive distorter. The EMTS is able to compensate for dynamic distortions even when previous calibrations have not been made in the presence of the fixed distorting object. The proposed method is based on the difference in the transmittance between the voltage induced directly in the sensing coil and the voltage induced indirectly in the sensing coil across a conductive distorter. The undistorted voltage in the emitting coil depends on the first derivative of the voltage in the emitting coil. In contrast, the distorted voltage depends on its second derivative.</p>
<p>In many applications of computer-assisted surgery, the sensing coil is placed in direct proximity of conductive surgical instruments or even attached to such instruments. The proposed method was evaluated using an aluminum disk to represent those surgical instruments used in minimal invasive surgery like hammer, scalpel or ultrasonic probe.</p>
<p>In the performed experiment, the influence of the conductive distorter on the accuracy of the electromagnetic tracking system was measured simultaneously with a quadratic and a sinusoidal excitation. The measurements in the distance between the sensing coil and aluminum plate smaller than 20 mm show a non-significant impact of the conductive distorter on the accuracy of the electromagnetic tracking system when supplied with a quadratic excitation (Mann–Whitney <inline-formula id="IEq31"><alternatives><tex-math id="M51">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$U$$\end{document}</tex-math>
<inline-graphic xlink:href="11548_2013_925_Article_IEq31.gif"></inline-graphic>
</alternatives>
</inline-formula>
-test with the significance level 0.01). The same experiments performed with the sinusoidal excitation show a significant influence of the conductive distorter on its accuracy. On the other hand, it should be noted that the position estimation with the quadratic excitation is less accurate without additional distortions due to the conductive objects.</p>
<p>The methods proposed in the literature require parameter calculations of the conductive distorter equivalent circuit. These calculations are complicated and lead to high inaccuracies, especially when there are multiple distorters to be modeled as multiple RL circuits. In the proposed method, it is not necessary to calculate these parameters of the conductive distorter equivalent circuit, and therefore, the proposed method is applicable when multiple distorters with complicated geometry are present—a common situation in surgical and radiological suites.</p>
<p>The proposed method focuses on the distortions due to the conductive distortions and does not take into consideration distortions due to ferromagnetic distorters.</p>
</sec>
</body>
<back><ack><title>Conflict of interest</title>
<p>The Tomasz Bien, Mengfei Li, Zein Salah and Georg Rose declare that they have no conflict of interest.</p>
</ack>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Birkfellner</surname>
<given-names>W</given-names>
</name>
<name><surname>Hummel</surname>
<given-names>J</given-names>
</name>
<name><surname>Wilson</surname>
<given-names>E</given-names>
</name>
<name><surname>Cleary</surname>
<given-names>K</given-names>
</name>
</person-group>
<person-group person-group-type="editor"><name><surname>Peters</surname>
<given-names>T</given-names>
</name>
<name><surname>Cleary</surname>
<given-names>K</given-names>
</name>
</person-group>
<article-title>Tracking devices</article-title>
<source>Image-guided interventions</source>
<year>2008</year>
<publisher-loc>US, New York</publisher-loc>
<publisher-name>Springer</publisher-name>
</element-citation>
</ref>
<ref id="CR2"><label>2.</label>
<mixed-citation publication-type="other">Nafis C, Jensen V, Jako R (2008) Method for evaluating compatibility of commercial electromagnetic (EM) micro sensor tracking systems with surgical and imaging tables. Proc SPIE. doi:10.1117/12.769513</mixed-citation>
</ref>
<ref id="CR3"><label>3.</label>
<mixed-citation publication-type="other">Nafis C, Jensen V, Beauregard L, Anderson P (2006) Method for estimating dynamic EM tracking accuracy of surgical navigation Tools. Proc SPIE. doi:10.1117/12.653448</mixed-citation>
</ref>
<ref id="CR4"><label>4.</label>
<mixed-citation publication-type="other">Traub J, Kaur S, Kneschaurek P, Navab N (2007), Evaluation of electromagnetic error correction methods—correcting distortion fields for appliance in the radiation therapy room. Bildverarbeitung für die Medizin. Available: <ext-link ext-link-type="uri" xlink:href="http://citeseerx.ist.psu.edu/viewdoc/download?">http://citeseerx.ist.psu.edu/viewdoc/download?</ext-link>
. doi:10.1.1.141.936&rep=rep1&type=pdf</mixed-citation>
</ref>
<ref id="CR5"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Weiser</surname>
<given-names>T</given-names>
</name>
<name><surname>Hyman</surname>
<given-names>K</given-names>
</name>
<name><surname>Yun</surname>
<given-names>J</given-names>
</name>
<name><surname>Litle</surname>
<given-names>V</given-names>
</name>
<name><surname>Chin</surname>
<given-names>C</given-names>
</name>
<name><surname>Swanson</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Electromagnetic navigational bronchoscopy: a surgeon’s perspective</article-title>
<source>Ann Thorac Surg</source>
<year>2008</year>
<volume>85</volume>
<fpage>797</fpage>
<lpage>801</lpage>
<pub-id pub-id-type="doi">10.1016/j.athoracsur.2007.11.052</pub-id>
</element-citation>
</ref>
<ref id="CR6"><label>6.</label>
<mixed-citation publication-type="other">Levy EB, Tang J, Lindisch D, Glossop N, Banovac F, Cleary K (2007) Implementation of an electromagnetic tracking system for accurate intrahepatic puncture needle guidance: accuracy results in an in vitro model. Acad Radiol 14(3):344–354</mixed-citation>
</ref>
<ref id="CR7"><label>7.</label>
<mixed-citation publication-type="other">Zhou J, Sebastian E, Mangona V, Yan D (2013) Real-time catheter tracking for high-dose-rate prostate brachytherapy using an electromagnetic 3D-guidance device: A preliminary performance study. Med Phys 40(2):021716. <ext-link ext-link-type="uri" xlink:href="http://online.medphys.org/resource/1/mphya6/v40/i2/p021716_s1?isAuthorized=no">http://online.medphys.org/resource/1/mphya6/v40/i2/p021716_s1?isAuthorized=no</ext-link>
</mixed-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>LaScalza</surname>
<given-names>S</given-names>
</name>
<name><surname>Arico</surname>
<given-names>J</given-names>
</name>
<name><surname>Hughes</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Effect of metal and sampling rate on accuracy of Flock of Birds electromagnetic tracking system</article-title>
<source>J Biomech</source>
<year>2003</year>
<volume>36</volume>
<issue>1</issue>
<fpage>141</fpage>
<lpage>144</lpage>
<pub-id pub-id-type="doi">10.1016/S0021-9290(02)00322-6</pub-id>
<pub-id pub-id-type="pmid">12485650</pub-id>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Birkfellner</surname>
<given-names>W</given-names>
</name>
<name><surname>Watzinger</surname>
<given-names>F</given-names>
</name>
<name><surname>Wanschitz</surname>
<given-names>F</given-names>
</name>
<name><surname>Ewers</surname>
<given-names>R</given-names>
</name>
<name><surname>Bergmann</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Calibration of tracking systems in a surgical environment</article-title>
<source>IEEE Trans Med Imag</source>
<year>1998</year>
<volume>17</volume>
<issue>5</issue>
<fpage>737</fpage>
<lpage>742</lpage>
<pub-id pub-id-type="doi">10.1109/42.736028</pub-id>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schuler</surname>
<given-names>N</given-names>
</name>
<name><surname>Bey</surname>
<given-names>M</given-names>
</name>
<name><surname>Shearn</surname>
<given-names>J</given-names>
</name>
<name><surname>Butler</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Evaluation of an electromagnetic position tracking for measuring in vivo, joint dynamic kinematics</article-title>
<source>J Biomech</source>
<year>2005</year>
<volume>38</volume>
<issue>10</issue>
<fpage>2113</fpage>
<lpage>2117</lpage>
<pub-id pub-id-type="doi">10.1016/j.jbiomech.2004.09.015</pub-id>
<pub-id pub-id-type="pmid">16084212</pub-id>
</element-citation>
</ref>
<ref id="CR11"><label>11.</label>
<mixed-citation publication-type="other">Raab F, Blood E, Steiner T, Jones H (1979) Magnetic position and orientation tracking systems. IEEE Trans Aerosp Electron Syst 15(5):709–718</mixed-citation>
</ref>
<ref id="CR12"><label>12.</label>
<mixed-citation publication-type="other">Seiler P, Blattmann H, Kirsch S, Schilling Ch (2000) A novel tracking technique for the continuous precise measurement of tumour position in conformal radiotherapy. Phys Med Biol 45: 103–110</mixed-citation>
</ref>
<ref id="CR13"><label>13.</label>
<mixed-citation publication-type="other">Schneider M, Stevens Ch (2007) Development and testing of a new magnetic-tracking device for image guidance. Proc SPIE. doi:10.1117/12.713249</mixed-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Poulin</surname>
<given-names>F</given-names>
</name>
<name><surname>Amiot</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>Interference during the use of an electromagnetic tracking system under OR conditions</article-title>
<source>J Biomech</source>
<year>2002</year>
<volume>35</volume>
<issue>6</issue>
<fpage>733</fpage>
<lpage>737</lpage>
<pub-id pub-id-type="doi">10.1016/S0021-9290(02)00036-2</pub-id>
<pub-id pub-id-type="pmid">12020992</pub-id>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hummel</surname>
<given-names>J</given-names>
</name>
<name><surname>Figl</surname>
<given-names>M</given-names>
</name>
<name><surname>Birkfellner</surname>
<given-names>W</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Evaluation of a new electromagnetic tracking system using a standardized assessment protocol</article-title>
<source>Phys Med Biol</source>
<year>2006</year>
<volume>51</volume>
<fpage>205</fpage>
<lpage>210</lpage>
<pub-id pub-id-type="doi">10.1088/0031-9155/51/10/N01</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<mixed-citation publication-type="other">Wilson E, Yaniv Z, Lindisch D, Cleary K (2008) A buyer’s guide to electromagnetic tracking systems for clinical applications. Proc SPIE. doi:10.1117/12.770509</mixed-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kindratenko</surname>
<given-names>V</given-names>
</name>
</person-group>
<article-title>A survey of electromagnetic position tracker calibration techniques</article-title>
<source>Virtual Real</source>
<year>2000</year>
<volume>5</volume>
<issue>3</issue>
<fpage>169</fpage>
<lpage>182</lpage>
<pub-id pub-id-type="doi">10.1007/BF01409422</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<mixed-citation publication-type="other">Mucha D (2009) An error compensating electromagnetic navigation for head surgery. Technische Universität Berlin, Dissertation</mixed-citation>
</ref>
<ref id="CR19"><label>19.</label>
<mixed-citation publication-type="other">Rodgers A, Higgins R, Hanf K, Farr J (2006), Distributed array magnetic tracking. Patent Application US 2006/0170417</mixed-citation>
</ref>
<ref id="CR20"><label>20.</label>
<mixed-citation publication-type="other">Scully J, Schneider M (2007) DC magnetic-based position and orientation monitoring system for tracking medical instruments. Patent Application US 2007/0078334A1</mixed-citation>
</ref>
<ref id="CR21"><label>21.</label>
<mixed-citation publication-type="other">Bar-Tal M, Govari A (2006) Distortion-immune position tracking using redundant measurements. European patent application, EP 1887309A1</mixed-citation>
</ref>
<ref id="CR22"><label>22.</label>
<mixed-citation publication-type="other">Lewandowski R, Wier E (1998) Metal immune magnetic tracker. US Patent 6,154,024</mixed-citation>
</ref>
<ref id="CR23"><label>23.</label>
<mixed-citation publication-type="other">Shalgi A, Nitzan Y, Tal M, Yaron U (2006) Reduced field distortion in medical tools. Patent Application US 2007/0255132A1</mixed-citation>
</ref>
<ref id="CR24"><label>24.</label>
<mixed-citation publication-type="other">Nieminen J, Kirsch S (2006) Eddy current detection and compensation. US Patent 7,353,125B2</mixed-citation>
</ref>
<ref id="CR25"><label>25.</label>
<mixed-citation publication-type="other">Schneider M (2010) System and method for measuring position and orientation using distortion compensated magnetic fields. US Patent 7,788,060B2 </mixed-citation>
</ref>
<ref id="CR26"><label>26.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Plotkin</surname>
<given-names>A</given-names>
</name>
<name><surname>Kucher</surname>
<given-names>V</given-names>
</name>
<name><surname>Horen</surname>
<given-names>Y</given-names>
</name>
<name><surname>Paperno</surname>
<given-names>E</given-names>
</name>
</person-group>
<article-title>A new calibration procedure for magnetic tracking systems</article-title>
<source>IEEE Trans Magn</source>
<year>2008</year>
<volume>44</volume>
<issue>11</issue>
<fpage>4525</fpage>
<lpage>4528</lpage>
<pub-id pub-id-type="doi">10.1109/TMAG.2008.2003056</pub-id>
</element-citation>
</ref>
<ref id="CR27"><label>27.</label>
<element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Fitzpatrick</surname>
<given-names>JM</given-names>
</name>
<name><surname>Hill</surname>
<given-names>DL</given-names>
</name>
</person-group>
<person-group person-group-type="editor"><name><surname>Sonka</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Image registration</article-title>
<source>Medical image processing and analysis</source>
<year>2000</year>
<publisher-loc>Washington</publisher-loc>
<publisher-name>Bellingham</publisher-name>
<fpage>488</fpage>
<lpage>496</lpage>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/TelematiV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000124 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000124 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Ticri/CIDE |area= TelematiV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:3955490 |texte= Electromagnetic tracking system with reduced distortion using quadratic excitation }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:23918006" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a TelematiV1
![]() | This area was generated with Dilib version V0.6.31. | ![]() |