Serveur d'exploration sur les dispositifs haptiques

Attention, ce site est en cours de développement !
Attention, site généré par des moyens informatiques à partir de corpus bruts.
Les informations ne sont donc pas validées.

Towards Active Tracking of Beating Heart Motion in the Presence of Arrhythmia for Robotic Assisted Beating Heart Surgery

Identifieur interne : 000856 ( Pmc/Checkpoint ); précédent : 000855; suivant : 000857

Towards Active Tracking of Beating Heart Motion in the Presence of Arrhythmia for Robotic Assisted Beating Heart Surgery

Auteurs : E. Erdem Tuna [États-Unis] ; Jamshid H. Karimov [États-Unis] ; Taoming Liu [États-Unis] ; Özkan Bebek [Turquie] ; Kiyotaka Fukamachi [États-Unis] ; M. Cenk Çavu O Lu [États-Unis]

Source :

RBID : PMC:4105597

Abstract

In robotic assisted beating heart surgery, the control architecture for heart motion tracking has stringent requirements in terms of bandwidth of the motion that needs to be tracked. In order to achieve sufficient tracking accuracy, feed-forward control algorithms, which rely on estimations of upcoming heart motion, have been proposed in the literature. However, performance of these feed-forward motion control algorithms under heart rhythm variations is an important concern. In their past work, the authors have demonstrated the effectiveness of a receding horizon model predictive control-based algorithm, which used generalized adaptive predictors, under constant and slowly varying heart rate conditions. This paper extends these studies to the case when the heart motion statistics change abruptly and significantly, such as during arrhythmias. A feasibility study is carried out to assess the motion tracking capabilities of the adaptive algorithms in the occurrence of arrhythmia during beating heart surgery. Specifically, the tracking performance of the algorithms is evaluated on prerecorded motion data, which is collected in vivo and includes heart rhythm irregularities. The algorithms are tested using both simulations and bench experiments on a three degree-of-freedom robotic test bed. They are also compared with a position-plus-derivative controller as well as a receding horizon model predictive controller that employs an extended Kalman filter algorithm for predicting future heart motion.


Url:
DOI: 10.1371/journal.pone.0102877
PubMed: 25048462
PubMed Central: 4105597


Affiliations:


Links toward previous steps (curation, corpus...)


Links to Exploration step

PMC:4105597

Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">Towards Active Tracking of Beating Heart Motion in the Presence of Arrhythmia for Robotic Assisted Beating Heart Surgery</title>
<author>
<name sortKey="Tuna, E Erdem" sort="Tuna, E Erdem" uniqKey="Tuna E" first="E. Erdem" last="Tuna">E. Erdem Tuna</name>
<affiliation wicri:level="2">
<nlm:aff id="aff1">
<addr-line>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Karimov, Jamshid H" sort="Karimov, Jamshid H" uniqKey="Karimov J" first="Jamshid H." last="Karimov">Jamshid H. Karimov</name>
<affiliation wicri:level="2">
<nlm:aff id="aff2">
<addr-line>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Liu, Taoming" sort="Liu, Taoming" uniqKey="Liu T" first="Taoming" last="Liu">Taoming Liu</name>
<affiliation wicri:level="2">
<nlm:aff id="aff1">
<addr-line>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Bebek, Ozkan" sort="Bebek, Ozkan" uniqKey="Bebek O" first="Özkan" last="Bebek">Özkan Bebek</name>
<affiliation wicri:level="1">
<nlm:aff id="aff3">
<addr-line>Department of Mechanical Engineering, Özyeğin University, Istanbul, Turkey</addr-line>
</nlm:aff>
<country xml:lang="fr">Turquie</country>
<wicri:regionArea>Department of Mechanical Engineering, Özyeğin University, Istanbul</wicri:regionArea>
<wicri:noRegion>Istanbul</wicri:noRegion>
</affiliation>
</author>
<author>
<name sortKey="Fukamachi, Kiyotaka" sort="Fukamachi, Kiyotaka" uniqKey="Fukamachi K" first="Kiyotaka" last="Fukamachi">Kiyotaka Fukamachi</name>
<affiliation wicri:level="2">
<nlm:aff id="aff2">
<addr-line>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Cavu O Lu, M Cenk" sort="Cavu O Lu, M Cenk" uniqKey="Cavu O Lu M" first="M. Cenk" last="Çavu O Lu">M. Cenk Çavu O Lu</name>
<affiliation wicri:level="2">
<nlm:aff id="aff1">
<addr-line>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">25048462</idno>
<idno type="pmc">4105597</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4105597</idno>
<idno type="RBID">PMC:4105597</idno>
<idno type="doi">10.1371/journal.pone.0102877</idno>
<date when="2014">2014</date>
<idno type="wicri:Area/Pmc/Corpus">002364</idno>
<idno type="wicri:Area/Pmc/Curation">002364</idno>
<idno type="wicri:Area/Pmc/Checkpoint">000856</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">Towards Active Tracking of Beating Heart Motion in the Presence of Arrhythmia for Robotic Assisted Beating Heart Surgery</title>
<author>
<name sortKey="Tuna, E Erdem" sort="Tuna, E Erdem" uniqKey="Tuna E" first="E. Erdem" last="Tuna">E. Erdem Tuna</name>
<affiliation wicri:level="2">
<nlm:aff id="aff1">
<addr-line>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Karimov, Jamshid H" sort="Karimov, Jamshid H" uniqKey="Karimov J" first="Jamshid H." last="Karimov">Jamshid H. Karimov</name>
<affiliation wicri:level="2">
<nlm:aff id="aff2">
<addr-line>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Liu, Taoming" sort="Liu, Taoming" uniqKey="Liu T" first="Taoming" last="Liu">Taoming Liu</name>
<affiliation wicri:level="2">
<nlm:aff id="aff1">
<addr-line>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Bebek, Ozkan" sort="Bebek, Ozkan" uniqKey="Bebek O" first="Özkan" last="Bebek">Özkan Bebek</name>
<affiliation wicri:level="1">
<nlm:aff id="aff3">
<addr-line>Department of Mechanical Engineering, Özyeğin University, Istanbul, Turkey</addr-line>
</nlm:aff>
<country xml:lang="fr">Turquie</country>
<wicri:regionArea>Department of Mechanical Engineering, Özyeğin University, Istanbul</wicri:regionArea>
<wicri:noRegion>Istanbul</wicri:noRegion>
</affiliation>
</author>
<author>
<name sortKey="Fukamachi, Kiyotaka" sort="Fukamachi, Kiyotaka" uniqKey="Fukamachi K" first="Kiyotaka" last="Fukamachi">Kiyotaka Fukamachi</name>
<affiliation wicri:level="2">
<nlm:aff id="aff2">
<addr-line>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
<author>
<name sortKey="Cavu O Lu, M Cenk" sort="Cavu O Lu, M Cenk" uniqKey="Cavu O Lu M" first="M. Cenk" last="Çavu O Lu">M. Cenk Çavu O Lu</name>
<affiliation wicri:level="2">
<nlm:aff id="aff1">
<addr-line>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio, United States of America</addr-line>
</nlm:aff>
<country xml:lang="fr">États-Unis</country>
<wicri:regionArea>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio</wicri:regionArea>
<placeName>
<region type="state">Ohio</region>
</placeName>
</affiliation>
</author>
</analytic>
<series>
<title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint>
<date when="2014">2014</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<p>In robotic assisted beating heart surgery, the control architecture for heart motion tracking has stringent requirements in terms of bandwidth of the motion that needs to be tracked. In order to achieve sufficient tracking accuracy, feed-forward control algorithms, which rely on estimations of upcoming heart motion, have been proposed in the literature. However, performance of these feed-forward motion control algorithms under heart rhythm variations is an important concern. In their past work, the authors have demonstrated the effectiveness of a receding horizon model predictive control-based algorithm, which used generalized adaptive predictors, under constant and slowly varying heart rate conditions. This paper extends these studies to the case when the heart motion statistics change abruptly and significantly, such as during arrhythmias. A feasibility study is carried out to assess the motion tracking capabilities of the adaptive algorithms in the occurrence of arrhythmia during beating heart surgery. Specifically, the tracking performance of the algorithms is evaluated on prerecorded motion data, which is collected
<italic>in vivo</italic>
and includes heart rhythm irregularities. The algorithms are tested using both simulations and bench experiments on a three degree-of-freedom robotic test bed. They are also compared with a position-plus-derivative controller as well as a receding horizon model predictive controller that employs an extended Kalman filter algorithm for predicting future heart motion.</p>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct>
<analytic>
<author>
<name sortKey="Ginhoux, R" uniqKey="Ginhoux R">R Ginhoux</name>
</author>
<author>
<name sortKey="Gangloff, Ja" uniqKey="Gangloff J">JA Gangloff</name>
</author>
<author>
<name sortKey="Demathelin, Mf" uniqKey="Demathelin M">MF DeMathelin</name>
</author>
<author>
<name sortKey="Soler, L" uniqKey="Soler L">L Soler</name>
</author>
<author>
<name sortKey="Leroy, J" uniqKey="Leroy J">J Leroy</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Richa, R" uniqKey="Richa R">R Richa</name>
</author>
<author>
<name sortKey="B, Ap" uniqKey="B A">AP Bó</name>
</author>
<author>
<name sortKey="Poignet, P" uniqKey="Poignet P">P Poignet</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Yuen, Sg" uniqKey="Yuen S">SG Yuen</name>
</author>
<author>
<name sortKey="Kettler, Dt" uniqKey="Kettler D">DT Kettler</name>
</author>
<author>
<name sortKey="Novotny, Pm" uniqKey="Novotny P">PM Novotny</name>
</author>
<author>
<name sortKey="Plowes, Rd" uniqKey="Plowes R">RD Plowes</name>
</author>
<author>
<name sortKey="Howe, Rd" uniqKey="Howe R">RD Howe</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bachta, W" uniqKey="Bachta W">W Bachta</name>
</author>
<author>
<name sortKey="Renaud, P" uniqKey="Renaud P">P Renaud</name>
</author>
<author>
<name sortKey="Cuvillon, L" uniqKey="Cuvillon L">L Cuvillon</name>
</author>
<author>
<name sortKey="Laroche, E" uniqKey="Laroche E">E Laroche</name>
</author>
<author>
<name sortKey="Forgione, A" uniqKey="Forgione A">A Forgione</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bogatyrenko, E" uniqKey="Bogatyrenko E">E Bogatyrenko</name>
</author>
<author>
<name sortKey="Pompey, P" uniqKey="Pompey P">P Pompey</name>
</author>
<author>
<name sortKey="Hanebeck, U" uniqKey="Hanebeck U">U Hanebeck</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Ortmaier, T" uniqKey="Ortmaier T">T Ortmaier</name>
</author>
<author>
<name sortKey="Groeger, M" uniqKey="Groeger M">M Groeger</name>
</author>
<author>
<name sortKey="Boehm, Dh" uniqKey="Boehm D">DH Boehm</name>
</author>
<author>
<name sortKey="Falk, V" uniqKey="Falk V">V Falk</name>
</author>
<author>
<name sortKey="Hirzinger, G" uniqKey="Hirzinger G">G Hirzinger</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bebek, O" uniqKey="Bebek O">O Bebek</name>
</author>
<author>
<name sortKey="Cavusoglu, Mc" uniqKey="Cavusoglu M">MC Cavusoglu</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Tuna, Ee" uniqKey="Tuna E">EE Tuna</name>
</author>
<author>
<name sortKey="Franke, Tj" uniqKey="Franke T">TJ Franke</name>
</author>
<author>
<name sortKey="Bebek, O" uniqKey="Bebek O">O Bebek</name>
</author>
<author>
<name sortKey="Shiose, A" uniqKey="Shiose A">A Shiose</name>
</author>
<author>
<name sortKey="Fukamachi, K" uniqKey="Fukamachi K">K Fukamachi</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Newman, Mf" uniqKey="Newman M">MF Newman</name>
</author>
<author>
<name sortKey="Kirchner, Jl" uniqKey="Kirchner J">JL Kirchner</name>
</author>
<author>
<name sortKey="Phillips Bute, B" uniqKey="Phillips Bute B">B Phillips-Bute</name>
</author>
<author>
<name sortKey="Gaver, V" uniqKey="Gaver V">V Gaver</name>
</author>
<author>
<name sortKey="Grocott, H" uniqKey="Grocott H">H Grocott</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Roka, A" uniqKey="Roka A">A Roka</name>
</author>
<author>
<name sortKey="Toth, E" uniqKey="Toth E">E Toth</name>
</author>
<author>
<name sortKey="Szilagyi, S" uniqKey="Szilagyi S">S Szilagyi</name>
</author>
<author>
<name sortKey="Merkely, B" uniqKey="Merkely B">B Merkely</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Takagaki, M" uniqKey="Takagaki M">M Takagaki</name>
</author>
<author>
<name sortKey="Mccarthy, Pm" uniqKey="Mccarthy P">PM McCarthy</name>
</author>
<author>
<name sortKey="Chung, M" uniqKey="Chung M">M Chung</name>
</author>
<author>
<name sortKey="Connor, J" uniqKey="Connor J">J Connor</name>
</author>
<author>
<name sortKey="Dessoffy, R" uniqKey="Dessoffy R">R Dessoffy</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Takagaki, M" uniqKey="Takagaki M">M Takagaki</name>
</author>
<author>
<name sortKey="Mccarthy, P" uniqKey="Mccarthy P">P McCarthy</name>
</author>
<author>
<name sortKey="Inoue, M" uniqKey="Inoue M">M Inoue</name>
</author>
<author>
<name sortKey="Chung, M" uniqKey="Chung M">M Chung</name>
</author>
<author>
<name sortKey="Connor, J" uniqKey="Connor J">J Connor</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Clarke, Dw" uniqKey="Clarke D">DW Clarke</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Garcia, Ce" uniqKey="Garcia C">CE Garcia</name>
</author>
<author>
<name sortKey="Prett, Dm" uniqKey="Prett D">DM Prett</name>
</author>
<author>
<name sortKey="Morari, M" uniqKey="Morari M">M Morari</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Cavusoglu, Mc" uniqKey="Cavusoglu M">MC Cavusoglu</name>
</author>
<author>
<name sortKey="Feygin, D" uniqKey="Feygin D">D Feygin</name>
</author>
<author>
<name sortKey="Tendick, F" uniqKey="Tendick F">F Tendick</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Ratcliffe, Mb" uniqKey="Ratcliffe M">MB Ratcliffe</name>
</author>
<author>
<name sortKey="Gupta, Kb" uniqKey="Gupta K">KB Gupta</name>
</author>
<author>
<name sortKey="Streicher, Jt" uniqKey="Streicher J">JT Streicher</name>
</author>
<author>
<name sortKey="Savage, Eb" uniqKey="Savage E">EB Savage</name>
</author>
<author>
<name sortKey="Bogen, Dk" uniqKey="Bogen D">DK Bogen</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bowthorpe, M" uniqKey="Bowthorpe M">M Bowthorpe</name>
</author>
<author>
<name sortKey="Tavakoli, M" uniqKey="Tavakoli M">M Tavakoli</name>
</author>
<author>
<name sortKey="Becher, H" uniqKey="Becher H">H Becher</name>
</author>
<author>
<name sortKey="Howe, R" uniqKey="Howe R">R Howe</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Liu, H" uniqKey="Liu H">H Liu</name>
</author>
<author>
<name sortKey="Shah, S" uniqKey="Shah S">S Shah</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Pearson, Rk" uniqKey="Pearson R">RK Pearson</name>
</author>
</analytic>
</biblStruct>
<biblStruct></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">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group>
<journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">25048462</article-id>
<article-id pub-id-type="pmc">4105597</article-id>
<article-id pub-id-type="publisher-id">PONE-D-14-11482</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0102877</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v2">
<subject>Engineering and Technology</subject>
<subj-group>
<subject>Mechanical Engineering</subject>
<subj-group>
<subject>Robotics</subject>
<subj-group>
<subject>Robotics Systems</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v2">
<subject>Medicine and Health Sciences</subject>
<subj-group>
<subject>Cardiology</subject>
<subj-group>
<subject>Cardiovascular Diseases</subject>
<subj-group>
<subject>Arrhythmia</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group>
<subject>Surgical and Invasive Medical Procedures</subject>
<subj-group>
<subject>Cardiovascular Procedures</subject>
<subj-group>
<subject>Cardiac Surgery</subject>
</subj-group>
</subj-group>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Towards Active Tracking of Beating Heart Motion in the Presence of Arrhythmia for Robotic Assisted Beating Heart Surgery</article-title>
<alt-title alt-title-type="running-head">Robotic Tracking of Heart Motion during Arrhythmia</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Tuna</surname>
<given-names>E. Erdem</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="cor1">
<sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Karimov</surname>
<given-names>Jamshid H.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Taoming</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Bebek</surname>
<given-names>Özkan</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fukamachi</surname>
<given-names>Kiyotaka</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Çavuşoğlu</surname>
<given-names>M. Cenk</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<label>1</label>
<addr-line>Department of Electrical Engineering and Computer Science, Case Western Reserve University, Cleveland, Ohio, United States of America</addr-line>
</aff>
<aff id="aff2">
<label>2</label>
<addr-line>Department of Biomedical Engineering, Lerner Research Institute, Cleveland Clinic, Cleveland, Ohio, United States of America</addr-line>
</aff>
<aff id="aff3">
<label>3</label>
<addr-line>Department of Mechanical Engineering, Özyeğin University, Istanbul, Turkey</addr-line>
</aff>
<contrib-group>
<contrib contrib-type="editor">
<name>
<surname>Talkachova</surname>
<given-names>Alena</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1">
<addr-line>University of Minnesota, United States of America</addr-line>
</aff>
<author-notes>
<corresp id="cor1">* E-mail:
<email>eet12@case.edu</email>
</corresp>
<fn fn-type="conflict">
<p>
<bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con">
<p>Conceived and designed the experiments: JHK OB KF MCC. Performed the experiments: EET JHK TL OB KF MCC. Analyzed the data: EET. Contributed reagents/materials/analysis tools: EET JHK TL OB KF MCC. Wrote the paper: EET JHK TL OB KF MCC.</p>
</fn>
</author-notes>
<pub-date pub-type="collection">
<year>2014</year>
</pub-date>
<pub-date pub-type="epub">
<day>21</day>
<month>7</month>
<year>2014</year>
</pub-date>
<volume>9</volume>
<issue>7</issue>
<elocation-id>e102877</elocation-id>
<history>
<date date-type="received">
<day>13</day>
<month>3</month>
<year>2014</year>
</date>
<date date-type="accepted">
<day>24</day>
<month>6</month>
<year>2014</year>
</date>
</history>
<permissions>
<copyright-year>2014</copyright-year>
<copyright-holder>Tuna et al</copyright-holder>
<license>
<license-p>This is an open-access article distributed under the terms of the
<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License</ext-link>
, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<abstract>
<p>In robotic assisted beating heart surgery, the control architecture for heart motion tracking has stringent requirements in terms of bandwidth of the motion that needs to be tracked. In order to achieve sufficient tracking accuracy, feed-forward control algorithms, which rely on estimations of upcoming heart motion, have been proposed in the literature. However, performance of these feed-forward motion control algorithms under heart rhythm variations is an important concern. In their past work, the authors have demonstrated the effectiveness of a receding horizon model predictive control-based algorithm, which used generalized adaptive predictors, under constant and slowly varying heart rate conditions. This paper extends these studies to the case when the heart motion statistics change abruptly and significantly, such as during arrhythmias. A feasibility study is carried out to assess the motion tracking capabilities of the adaptive algorithms in the occurrence of arrhythmia during beating heart surgery. Specifically, the tracking performance of the algorithms is evaluated on prerecorded motion data, which is collected
<italic>in vivo</italic>
and includes heart rhythm irregularities. The algorithms are tested using both simulations and bench experiments on a three degree-of-freedom robotic test bed. They are also compared with a position-plus-derivative controller as well as a receding horizon model predictive controller that employs an extended Kalman filter algorithm for predicting future heart motion.</p>
</abstract>
<funding-group>
<funding-statement>This work was supported in part by NSF under grants IIS-0905344 and CNS-1035602, and National Institutes of Health under grants R21 HL096941 and R01 EB018108. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts>
<page-count count="8"></page-count>
</counts>
</article-meta>
</front>
<body>
<sec id="s1">
<title>Introduction</title>
<p>The present work is focused on a previously unexplored aspect of robotic-assisted beating heart surgery; namely, evaluating performance of robotic active relative motion cancellation (ARMC) control algorithms for heart motion tracking under atrial fibrillation (AF) induced arrhythmia conditions. The occurrence of arrhythmia during cardiac procedures cannot be ruled out, hence robots must be able to handle this case to be clinically useful. This paper is the first to acknowledge this essential problem of tracking heart motion in the presence of heart rhythm irregularities during arrhythmia.</p>
<p>There is substantial amount of literature in heart motion tracking studies for robotic-assisted beating heart surgery
<xref rid="pone.0102877-Ginhoux1" ref-type="bibr">[1]</xref>
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
(see
<italic>Background</italic>
for further details). However, the studies under slow heart rate variations are limited
<xref rid="pone.0102877-Yuen1" ref-type="bibr">[4]</xref>
,
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
and none of the heart motion tracking algorithms are evaluated under arrhythmia conditions.</p>
<p>The purpose of this study is to evaluate the feasibility of the ARMC control algorithms introduced previously in
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
under arrhythmia conditions. Specifically, a novel set of
<italic>in-vivo</italic>
heart motion data is acquired from a bovine model, where the data includes artificially induced atrial fibrillation rhythm irregularities. Then, the ARMC control algorithms are tested on this data using simulations and bench-top experiments on a hardware test-bed. The algotihms are further compared with an extended Kalman filter (EKF) based control algorithm and a position-plus-derivative (PD) control algorithm. Although the control algorithms are formerly described in
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
and the proposed approach for inducing arrhythmia might be a limited model for the clinical situation, the presented results provide novel contribution to the literature of heart motion tracking for robotic assisted beating heart surgery.</p>
<sec id="s1a">
<title>Background</title>
<p>Coronary artery bypass graft (CABG) surgery requires surgeons to operate on blood vessels that move with high bandwidth. This rapid motion of heart makes it difficult to track these arteries by hand effectively
<xref rid="pone.0102877-Trejos1" ref-type="bibr">[15]</xref>
. Contemporary techniques either stop the heart and use a cardio-pulmonary bypass machine,
<italic>on-pump</italic>
, or passively restrain the beating heart with mechanical stabilizers,
<italic>off-pump</italic>
, in order to cancel the biological motion of the heart during CABG surgery. However, using on-pump CABG surgery might cause significant complications that might occur during or after surgery, which includes long-term cognitive loss
<xref rid="pone.0102877-Newman1" ref-type="bibr">[16]</xref>
. Off-pump CABG surgery is mostly limited to the front surface of the heart and significant residual motion is observed during stabilization
<xref rid="pone.0102877-Lemma1" ref-type="bibr">[17]</xref>
.</p>
<p>Robotic-assisted beating heart surgery has been proposed to enable high dexterity surgical manipulations to be performed on a beating heart without mechanical stabilization. In robotic-assisted beating heart surgery conventional surgical tools are replaced with robotic instruments, which are directly controlled by the surgeon through teleoperation. In this system, the surgeon views the surgical site through a camera mounted on a robotic arm, which follows the heart motion, giving a stabilized view. The robotic surgical instruments also track the heart motion, canceling the relative motion between the surgical site and the instruments. As a result, the surgeon operates on the heart as if it was motionless, while the robotic system actively compensates the relative motion of the heart
<xref rid="pone.0102877-Nakamura1" ref-type="bibr">[18]</xref>
.</p>
<p>Due to the high bandwidth of the heart motion, it is necessary to employ model predictive feed-forward control algorithms, which rely on estimation of the future motion of the point-of-interest (POI) to achieve sufficient tracking accuracy
<xref rid="pone.0102877-Ginhoux1" ref-type="bibr">[1]</xref>
,
<xref rid="pone.0102877-Yuen1" ref-type="bibr">[4]</xref>
,
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
. However, the performance of these algorithms under changes in heart rhythm is a valid and important concern.</p>
<p>Several approaches have been proposed in the literature for modeling and predicting the motion of a POI on the heart, including algorithms using harmonics models
<xref rid="pone.0102877-Ginhoux1" ref-type="bibr">[1]</xref>
, extended Kalman filters
<xref rid="pone.0102877-Richa1" ref-type="bibr">[2]</xref>
<xref rid="pone.0102877-Yuen1" ref-type="bibr">[4]</xref>
, coupled breathing and heartbeat motion models
<xref rid="pone.0102877-Bachta1" ref-type="bibr">[5]</xref>
,
<xref rid="pone.0102877-Cuvillon1" ref-type="bibr">[6]</xref>
, surface deformation models
<xref rid="pone.0102877-Bader1" ref-type="bibr">[7]</xref>
<xref rid="pone.0102877-Bogatyrenko2" ref-type="bibr">[9]</xref>
, ECG synchronized periodic models
<xref rid="pone.0102877-Ortmaier1" ref-type="bibr">[10]</xref>
,
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
, and generalized adaptive filters
<xref rid="pone.0102877-Franke1" ref-type="bibr">[12]</xref>
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
. However, heart motion tracking performance using these algorithms under heart rhythm irregularities has not been addressed in the literature.</p>
<p>In our previous work two least-squares-based prediction algorithms, namely one-step and generalized, using an adaptive filter to generate future position estimates, were introduced
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
. These algorithms were studied by using a range of prerecorded
<italic>in vivo</italic>
constant and slowly varying heart rate motion data. Our previous report showed that the tracking of POI motion is no longer the bottleneck since the necessary amount of root-mean-square (RMS) tracking error on the order of 100–250
<italic> µm</italic>
is achieved. Furthermore, if the heart behavior slowly changes, then adaptive predictors are able to adjust to these changes quickly enough and yield good tracking results.</p>
<p>This paper extends our previous efforts on active tracking of the beating heart to the case of arrhythmia presence, where heart behavior changes significantly and abruptly. Here, a feasibility study is performed to determine motion tracking capabilities of these algorithms in the occurrence of arrhythmia. We evaluate the tracking performance of one-step and generalized predictors presented in
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
used in conjunction with the receding horizon model predictive control algorithm (RHMPC) presented in
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
with a range of prerecorded
<italic>in vivo</italic>
arrhythmia data by simulations and on a three degree-of-freedom (3-DOF) robotic test bed. Additionally, tracking results of RHMPC with an EKF (adapted from the implementation presented in
<xref rid="pone.0102877-Yuen1" ref-type="bibr">[4]</xref>
) and of a PD controller are given for comparison purposes.</p>
<p>The rest of this paper is organized as follows. First, experimental setup and analysis of the arrhythmia data is described. Then, the adaptive prediction methods are briefly revisited, and the results and the comparison of the algorithms are presented. Finally, conclusions and possible extensions are given. The adaptive prediction algorithms are proposed in
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
; the applications of these algorithms to arrhythmia tracking is the original contribution and has not been published previously.</p>
</sec>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>Materials and Methods</title>
<sec id="s2a">
<title>Analysis of Heart Data</title>
<p>The heart motion data were collected
<italic>in vivo</italic>
from three calves using a sonomicrometer and all of the bench-top experiments were performed with these prerecorded data. The size of the hearts were approximately 16
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e001.jpg"></inline-graphic>
</inline-formula>
2 cm in length and 12
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e002.jpg"></inline-graphic>
</inline-formula>
1 cm in width. In the experimental set-up for measurement of heart motion, two sonomicrometry crystals were placed on the epicardial surface. One crystal was sutured to the left side of the left anterior descending artery (LAD), whose location is referred to as “Anterior” in the rest of the paper. The second crystal was sutured on the right side of the LAD, whose location is referred as “Lateral” in the paper. Eight other crystals were asymmetrically mounted upon a rigid plastic base of diameter 60 mm, on a circle of diameter 50 mm, forming a reference coordinate frame. This rigid plastic sensor base was placed in a rubber latex balloon, which was filled with a 9.5% glycerine solution. The reason of using such a set-up was to ensure a continuous line of sight between the base crystals and the crystal on the heart surface through a liquid medium for proper operation of the sonomicrometry sensor system.</p>
<p>Two pacemaker electrodes were sutured to the right atrium and connected to a pacemaker to set a desired cardiac rhythm and induce arrhythmias. The method of inducing atrial fibrillation (AF) by electrically stimulating right atrium has been used in the literature for easy and reliable modeling of AF
<xref rid="pone.0102877-Roka1" ref-type="bibr">[19]</xref>
<xref rid="pone.0102877-Takagaki2" ref-type="bibr">[21]</xref>
. The AF was induced by stimulating the right atrium via two electrodes connected to an Electric Transformer with 2.5V AC @60 Hz. The ECG pattern of the AF (i.e. the absence of P waves, disorganized electrical activity in their place, and irregular R-R intervals due to irregular conduction iimparted to the ventricles) induced intraoperatively in our experimental setting was similar to that found clinically in patients
<xref rid="pone.0102877-Roka1" ref-type="bibr">[19]</xref>
<xref rid="pone.0102877-Takagaki2" ref-type="bibr">[21]</xref>
. This set-up for arrhythmia generation was chosen for motion control and tracking purposes only and not for the functional aspect of rhythm disturbances. The formulation of the biochemical and the pathophysiological aspects of the cardiac muscle motion, and the related neurostimulation are beyond the scope of the present study.</p>
<p>Data were processed offline and only filtering performed on the data (using the proprietary software provided with the system) was very limited removal of the outliers, which occasionally occur as a result of ultrasound echoing effects. Description of the sonomicrometer crystal locations on heart with respect to the left anterior descending artery (LAD) are given in
<xref ref-type="table" rid="pone-0102877-t001">Table 1</xref>
.
<xref ref-type="fig" rid="pone-0102877-g001">Fig. 1</xref>
shows the experimental setup for data collection. The sonomicrometer crystals, the sonomicrometer base, and the pacemaker leads are visible. A Fourier transform (
<xref ref-type="fig" rid="pone-0102877-g002">Fig. 2</xref>
) of the 183 s heart signal data from animal 1 reveals the inherently non-periodic nature of heart motion during arrhythmia. The abundance of intermittent frequencies in the power spectral density (PSD) with the absence of tall, narrow peaks indicates this point.</p>
<fig id="pone-0102877-g001" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.g001</object-id>
<label>Figure 1</label>
<caption>
<title>Experimental setup for the measurement of the heart motion.</title>
<p>Two sonomicrometer crystals that are sutured on the anterior and posterior surfaces of the heart are used for data collection. Pacemaker leads and sonomicrometer base are also visible in the image.</p>
</caption>
<graphic xlink:href="pone.0102877.g001"></graphic>
</fig>
<fig id="pone-0102877-g002" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.g002</object-id>
<label>Figure 2</label>
<caption>
<title>Power spectral density of the heart motion in the z-direction.</title>
<p>Heart motion modes are inseparable. The frequency axis is set to 12-s arrhythmia data from animal 1 is shown. The spectrum corresponds to motion of the POI located at 0.5 cm on the right side of LAD.</p>
</caption>
<graphic xlink:href="pone.0102877.g002"></graphic>
</fig>
<table-wrap id="pone-0102877-t001" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.t001</object-id>
<label>Table 1</label>
<caption>
<title>Arrhythmia Data.</title>
</caption>
<alternatives>
<graphic id="pone-0102877-t001-1" xlink:href="pone.0102877.t001"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">DataSet (Body weight)</td>
<td colspan="2" align="left" rowspan="1">Animal 1 (75.8 kg)</td>
<td colspan="2" align="left" rowspan="1">Animal 2 (68.7 kg)</td>
<td colspan="2" align="left" rowspan="1">Animal 3 (70.2 kg)</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">
<bold>Location of crystal on cardiac surface</bold>
</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<bold>Position relative to LAD</bold>
</td>
<td align="left" rowspan="1" colspan="1">2 cm left</td>
<td align="left" rowspan="1" colspan="1">0.5 cm right</td>
<td align="left" rowspan="1" colspan="1">1 cm left</td>
<td align="left" rowspan="1" colspan="1">5 cm right</td>
<td align="left" rowspan="1" colspan="1">1 cm left</td>
<td align="left" rowspan="1" colspan="1">5 cm right</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="nt101">
<label></label>
<p>Data from animal 1 have a sampling rate of 180 Hz with durations 85 s and 183 s respectively for “Anterior” and “Lateral” locations. Data from animal 2 have a sampling rate of 404 Hz with a duration of 85 s for both locations. Data from animal 3 have a sampling rate of 404 Hz with durations respectively 215 s and 105 s for “Anterior” and “Lateral” locations.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s2b">
<title>Ethics Statement</title>
<p>The
<italic>in vivo</italic>
study protocol was approved by the Cleveland Clinics Institutional Animal Care and Use Committee, and all animals received humane care in compliance with the Guide for the Care and Use of Laboratory Animals (Institute of Laboratory Animal Resources, Commission on Life Sciences, National Research Council, National Academy Press, Washington, DC, 2011) and institutional guidelines.</p>
</sec>
<sec id="s2c">
<title>Adaptive Motion Estimation Algorithms</title>
<p>The most important aspects of robotic surgical instruments in robotic assisted CABG surgery are accurately measuring and predicting the heart motion as they are instrumental in tracking and canceling the relative motion between the heart surface and surgical tools attached to the robotic manipulators. The rapid motion of the heartbeat component possesses demanding requirements on the control architecture of the robotic system in terms of the bandwidth of the motion that needs to be tracked. Feed-forward control algorithms, which rely on estimations of upcoming heart motion, need to be utilized to achieve sufficient tracking accuracy
<xref rid="pone.0102877-Ginhoux1" ref-type="bibr">[1]</xref>
,
<xref rid="pone.0102877-Yuen1" ref-type="bibr">[4]</xref>
,
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
.</p>
<p>Two adaptive filter-based predictors used in this study provide an estimate of the immediate future of the POI motion over a prediction horizon to the feed-forward control algorithm. The best estimate is defined to be the one that minimizes the square of the estimation error, where the estimation error is the difference between the prediction and the observed value at that time. Once a method is established to predict the next observations, a sequence of future observations can be estimated.
<xref ref-type="fig" rid="pone-0102877-g003">Fig. 3</xref>
provides a graphical schematic of the prediction problem.</p>
<fig id="pone-0102877-g003" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.g003</object-id>
<label>Figure 3</label>
<caption>
<title>A schematic of the heart motion prediction problem.</title>
<p>The circles represent past observations, now in memory, the
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e003.jpg"></inline-graphic>
</inline-formula>
is the current observation, and the short curve originating from there is the horizon estimate. The predictor takes the past observations and produces the horizon estimate from past observations.</p>
</caption>
<graphic xlink:href="pone.0102877.g003"></graphic>
</fig>
<p>The implementations of the predictors parameterize a linear system to predict POI motion and rely on recursive least-squares (RLS) adaptive filter algorithms. The one-step predictor assumes a linear system relation between the consecutive samples in the prediction horizon, whereas the generalized method performs the parametrization of the linear system independently for each point throughout the horizon.</p>
<p>The estimated trajectory of the POI on the heart is used as a feed-forward control signal in a model predictive controller
<xref rid="pone.0102877-Clarke1" ref-type="bibr">[22]</xref>
<xref rid="pone.0102877-Camacho1" ref-type="bibr">[25]</xref>
. The RHMPC algorithm combines linear quadratic optimal control strategy with prediction. The robot model and the estimated trajectory of the POI, which extends into the future, provide the prediction. At each time step, the control action is calculated by solving a finite horizon linear quadratic optimal control problem, which compares the predicted plant signals with the provided desired trajectory and the control objectives for the given time horizon.</p>
<p>The details of the RHMPC algorithm can be found in
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
. The motion estimation problem and adaptive estimation algorithms are comprehensively explained in
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
. They have been omitted here due to space constraints.</p>
<p>In the authors' previous work
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
the effectiveness of the RHMPC algorithm, which employs adaptive predictors to estimate upcoming heart motion, was studied during regular and slowly-changing heart rate conditions. Building on this work, this study focuses on the problem of heart motion tracking in the presence of highly irregular heart rhythm conditions. Specifically, the tracking performance of the algorithms is evaluated with arrhythmia data. To the best of the authors' knowledge, there are no previous studies in the literature to address this problem.</p>
</sec>
</sec>
<sec id="s3">
<title>Results</title>
<p>It is necessary to evaluate algorithms on some type of hardware testbeds. The algorithms were tested on a PHANToM Premium 1.5A haptic device, which is a 3-DOF robotic system. The PHANToM provides insight into the effectiveness of the algorithms on a desired system. Its lightweight frame and drive system allow for sufficient motion and speed to attempt to track the heartbeat signal. With its high bandwidth and low inertia features, the PHANToM robot possesses characteristics similar to those desired in a surgical robot that would be capable of heart motion compensation. Such a surgical robot would need to have lightweight links, low inertia design, and low-friction actuation system for tracking the heart with sufficient motion and speed. Such a prototype is currently being developed
<xref rid="pone.0102877-Liu1" ref-type="bibr">[26]</xref>
.</p>
<p>In the experimental setup, control algorithms were executed in xPC Target and run in a real-time kernel with a sampling time of 0.5 ms on a Intel Xeon 2.33 GHz Core PC. The nonlinearities of the test-bed system (i.e., gravitational effects, joint frictions, and Coriolis and centrifugal forces) were canceled independently from the controller (more details can be found in
<xref rid="pone.0102877-Cavusoglu1" ref-type="bibr">[27]</xref>
). The robot was commanded to follow the combined motion of heartbeat and breathing. The system used online streaming position data in place of real-time measurements.</p>
<p>The robot model was controlled using RHMPC. The encoder positions on the PHANToM were recorded and these positions were transformed into end effector positions. The reported RMS errors, which represent the tracking performance of the algorithms, were calculated from the difference between the prerecorded target point and the actual end effector position calculated from joint angles.</p>
<p>In both simulations and experiments, the same control methods and reference data were used. During the trials, a
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e004.jpg"></inline-graphic>
</inline-formula>
order correlated signal one-step estimator and a
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e005.jpg"></inline-graphic>
</inline-formula>
order generalized estimator were used. The predictors were downsampled by a factor of 15, processing observations that were 7.5 ms apart
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
and used to predict a 25 ms horizon. The 25 ms horizon corresponds to 50 control samples into the future. The length of the control horizon and estimator orders were chosen for the optimum error/performance under real-time computation requirements
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
,
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
.</p>
<p>As sonomicrometer was employed to acquire heart motion data used in this study, the sensor measurement delays were not considered, while determining the 25 ms prediction horizon. The sonomicrtometry sensor operates at high update rates (
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e006.jpg"></inline-graphic>
</inline-formula>
125 Hz) and has a high spatial resolution (24
<italic>µm</italic>
)
<xref rid="pone.0102877-Dione1" ref-type="bibr">[28]</xref>
. The primary sources of measurement error include sonomicrometry crystal geometry, ultrasound echoes, and either weak or missing signals
<xref rid="pone.0102877-Ratcliffe1" ref-type="bibr">[29]</xref>
. In this sense, sonomicrometry system architecture does not cause significant data acquisition and processing delays, unlike the vision-based and the ultrasound-based sensory systems, where there are inherent to image acquisition and processing delays
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
,
<xref rid="pone.0102877-Bowthorpe1" ref-type="bibr">[30]</xref>
. When such image-based sensor modalities are employed as in
<xref rid="pone.0102877-Yuen1" ref-type="bibr">[4]</xref>
, the non-negligible delays introduced due to image acquisition and processing must be compensated. Thus, 25 ms estimation horizon is almost exclusively used for feed-forward compensation of the robot dynamics for improved control performance in this study.</p>
<p>For each case, experiments on the PHANToM robot were run 10 times with the estimation algorithms and again with the actual heart motion data as future signal reference for the prediction horizon. The latter case represents a ‘perfect’ estimation, providing a performance base of the robotic system's capability. It was noted that the deviation between the trials had been very small. Among these results, the maximum values for the
<italic>End-effector RMS and Maximum Position Errors in millimeters</italic>
in 3D and
<italic>RMS Control Effort in millinewton meters</italic>
are summarized respectively in
<xref ref-type="table" rid="pone-0102877-t002">Table 2</xref>
and
<xref ref-type="table" rid="pone-0102877-t003">Table 3</xref>
for the simulations and respectively in
<xref ref-type="table" rid="pone-0102877-t004">Table 4</xref>
and
<xref ref-type="table" rid="pone-0102877-t005">Table 5</xref>
for experiments to project the worst cases. The
<italic>RMS magnitude</italic>
and
<italic>peak-to-peak amplitude</italic>
of the heart motion for each data set are also stated to highlight the importance of using motion prediction for decreasing the tracking error.</p>
<table-wrap id="pone-0102877-t002" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.t002</object-id>
<label>Table 2</label>
<caption>
<title>Simulation Results for End-Effector Tracking: RMS End-Effector Error and MAX Position Error for the Control Algorithms.</title>
</caption>
<alternatives>
<graphic id="pone-0102877-t002-2" xlink:href="pone.0102877.t002"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">End-effector Tracking Results</td>
<td colspan="6" align="left" rowspan="1">RMS Position Error [mm]</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td colspan="6" align="left" rowspan="1">(Maximum Position Error [mm])</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">DataSet</td>
<td colspan="2" align="left" rowspan="1">Animal 1</td>
<td colspan="2" align="left" rowspan="1">Animal 2</td>
<td colspan="2" align="left" rowspan="1">Animal 3</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Crystal Location</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">
<bold>P-P amp of POI motion [mm]</bold>
</td>
<td align="left" rowspan="1" colspan="1">11.811</td>
<td align="left" rowspan="1" colspan="1">20.946</td>
<td align="left" rowspan="1" colspan="1">7.619</td>
<td align="left" rowspan="1" colspan="1">10.399</td>
<td align="left" rowspan="1" colspan="1">13.796</td>
<td align="left" rowspan="1" colspan="1">18.676</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<bold>RMS mag of POI motion [mm]</bold>
</td>
<td align="left" rowspan="1" colspan="1">4.488</td>
<td align="left" rowspan="1" colspan="1">4.419</td>
<td align="left" rowspan="1" colspan="1">3.238</td>
<td align="left" rowspan="1" colspan="1">3.063</td>
<td align="left" rowspan="1" colspan="1">4.819</td>
<td align="left" rowspan="1" colspan="1">4.058</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Exact Reference Information</td>
<td align="left" rowspan="1" colspan="1">0.201</td>
<td align="left" rowspan="1" colspan="1">0.283</td>
<td align="left" rowspan="1" colspan="1">0.188</td>
<td align="left" rowspan="1" colspan="1">0.224</td>
<td align="left" rowspan="1" colspan="1">0.392</td>
<td align="left" rowspan="1" colspan="1">0.310</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(1.066)</td>
<td align="left" rowspan="1" colspan="1">(1.689)</td>
<td align="left" rowspan="1" colspan="1">(1.132)</td>
<td align="left" rowspan="1" colspan="1">(1.062)</td>
<td align="left" rowspan="1" colspan="1">(1.934)</td>
<td align="left" rowspan="1" colspan="1">(4.479)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with One-Step Adaptive Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">0.204</td>
<td align="left" rowspan="1" colspan="1">0.300</td>
<td align="left" rowspan="1" colspan="1">0.193</td>
<td align="left" rowspan="1" colspan="1">0.239</td>
<td align="left" rowspan="1" colspan="1">0.388</td>
<td align="left" rowspan="1" colspan="1">0.334</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(1.603)</td>
<td align="left" rowspan="1" colspan="1">(2.431)</td>
<td align="left" rowspan="1" colspan="1">(1.325)</td>
<td align="left" rowspan="1" colspan="1">(1.554)</td>
<td align="left" rowspan="1" colspan="1">(1.970)</td>
<td align="left" rowspan="1" colspan="1">(6.223)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Generalized Adaptive Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.192</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.262</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.185</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.221</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.380</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.299</bold>
</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(1.402)</td>
<td align="left" rowspan="1" colspan="1">(1.751)</td>
<td align="left" rowspan="1" colspan="1">(1.054)</td>
<td align="left" rowspan="1" colspan="1">(1.012)</td>
<td align="left" rowspan="1" colspan="1">(1.792)</td>
<td align="left" rowspan="1" colspan="1">(3.671)</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<table-wrap id="pone-0102877-t003" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.t003</object-id>
<label>Table 3</label>
<caption>
<title>Simulation Results for End-Effector Tracking: RMS Control Effort for the Control Algorithms.</title>
</caption>
<alternatives>
<graphic id="pone-0102877-t003-3" xlink:href="pone.0102877.t003"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">End-effector Tracking Results</td>
<td colspan="6" align="left" rowspan="1">Control Effort [mNm]</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">DataSet</td>
<td colspan="2" align="left" rowspan="1">Animal 1</td>
<td colspan="2" align="left" rowspan="1">Animal 2</td>
<td colspan="2" align="left" rowspan="1">Animal 3</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Crystal Location</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Exact Reference Information</td>
<td align="left" rowspan="1" colspan="1">12.965</td>
<td align="left" rowspan="1" colspan="1">15.649</td>
<td align="left" rowspan="1" colspan="1">10.943</td>
<td align="left" rowspan="1" colspan="1">13.894</td>
<td align="left" rowspan="1" colspan="1">17.425</td>
<td align="left" rowspan="1" colspan="1">15.761</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with One-Step Adaptive Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">19.827</td>
<td align="left" rowspan="1" colspan="1">34.248</td>
<td align="left" rowspan="1" colspan="1">16.723</td>
<td align="left" rowspan="1" colspan="1">22.823</td>
<td align="left" rowspan="1" colspan="1">33.500</td>
<td align="left" rowspan="1" colspan="1">27.034</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Generalized Adaptive Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">15.294</td>
<td align="left" rowspan="1" colspan="1">21.475</td>
<td align="left" rowspan="1" colspan="1">11.861</td>
<td align="left" rowspan="1" colspan="1">15.095</td>
<td align="left" rowspan="1" colspan="1">23.003</td>
<td align="left" rowspan="1" colspan="1">20.195</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<table-wrap id="pone-0102877-t004" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.t004</object-id>
<label>Table 4</label>
<caption>
<title>Experimental Results for End-Effector Tracking: RMS End-Effector Error and MAX Position Error for the Control Algorithms.</title>
</caption>
<alternatives>
<graphic id="pone-0102877-t004-4" xlink:href="pone.0102877.t004"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">End-effector Tracking Results</td>
<td colspan="6" align="left" rowspan="1">RMS Position Error [mm]</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td colspan="6" align="left" rowspan="1">(Maximum Position Error [mm])</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">DataSet</td>
<td colspan="2" align="left" rowspan="1">Animal 1</td>
<td colspan="2" align="left" rowspan="1">Animal 2</td>
<td colspan="2" align="left" rowspan="1">Animal 3</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Crystal Location</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">
<bold>P-P amp of POI motion [mm]</bold>
</td>
<td align="left" rowspan="1" colspan="1">11.811</td>
<td align="left" rowspan="1" colspan="1">20.946</td>
<td align="left" rowspan="1" colspan="1">7.619</td>
<td align="left" rowspan="1" colspan="1">10.399</td>
<td align="left" rowspan="1" colspan="1">13.796</td>
<td align="left" rowspan="1" colspan="1">18.676</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">
<bold>RMS mag of POI motion [mm]</bold>
</td>
<td align="left" rowspan="1" colspan="1">4.488</td>
<td align="left" rowspan="1" colspan="1">4.419</td>
<td align="left" rowspan="1" colspan="1">3.238</td>
<td align="left" rowspan="1" colspan="1">3.063</td>
<td align="left" rowspan="1" colspan="1">4.819</td>
<td align="left" rowspan="1" colspan="1">4.058</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Exact Reference Information</td>
<td align="left" rowspan="1" colspan="1">0.260</td>
<td align="left" rowspan="1" colspan="1">0.301</td>
<td align="left" rowspan="1" colspan="1">0.274</td>
<td align="left" rowspan="1" colspan="1">0.248</td>
<td align="left" rowspan="1" colspan="1">0.303</td>
<td align="left" rowspan="1" colspan="1">0.309</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(0.809)</td>
<td align="left" rowspan="1" colspan="1">(1.933)</td>
<td align="left" rowspan="1" colspan="1">(1.412)</td>
<td align="left" rowspan="1" colspan="1">(0.779)</td>
<td align="left" rowspan="1" colspan="1">(1.276)</td>
<td align="left" rowspan="1" colspan="1">(9.191)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with One-Step Adaptive Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">0.293</td>
<td align="left" rowspan="1" colspan="1">0.403</td>
<td align="left" rowspan="1" colspan="1">0.278</td>
<td align="left" rowspan="1" colspan="1">0.263</td>
<td align="left" rowspan="1" colspan="1">0.405</td>
<td align="left" rowspan="1" colspan="1">0.388</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(1.459)</td>
<td align="left" rowspan="1" colspan="1">(5.911)</td>
<td align="left" rowspan="1" colspan="1">(1.643)</td>
<td align="left" rowspan="1" colspan="1">(1.513)</td>
<td align="left" rowspan="1" colspan="1">(3.383)</td>
<td align="left" rowspan="1" colspan="1">(8.289)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Generalized Adaptive Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.282</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.340</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.276</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.247</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.319</bold>
</td>
<td align="left" rowspan="1" colspan="1">
<bold>0.330</bold>
</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(1.038)</td>
<td align="left" rowspan="1" colspan="1">(2.703)</td>
<td align="left" rowspan="1" colspan="1">(1.405)</td>
<td align="left" rowspan="1" colspan="1">(1.174)</td>
<td align="left" rowspan="1" colspan="1">(1.442)</td>
<td align="left" rowspan="1" colspan="1">(8.620)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Extended Kalman Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">0.570</td>
<td align="left" rowspan="1" colspan="1">0.839</td>
<td align="left" rowspan="1" colspan="1">0.546</td>
<td align="left" rowspan="1" colspan="1">0.561</td>
<td align="left" rowspan="1" colspan="1">1.179</td>
<td align="left" rowspan="1" colspan="1">0.913</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(5.033)</td>
<td align="left" rowspan="1" colspan="1">(6.702)</td>
<td align="left" rowspan="1" colspan="1">(2.630)</td>
<td align="left" rowspan="1" colspan="1">(3.342)</td>
<td align="left" rowspan="1" colspan="1">(6.499)</td>
<td align="left" rowspan="1" colspan="1">(12.430)</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">PD Controller</td>
<td align="left" rowspan="1" colspan="1">0.465</td>
<td align="left" rowspan="1" colspan="1">0.775</td>
<td align="left" rowspan="1" colspan="1">0.419</td>
<td align="left" rowspan="1" colspan="1">0.484</td>
<td align="left" rowspan="1" colspan="1">1.130</td>
<td align="left" rowspan="1" colspan="1">0.886</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">(3.674)</td>
<td align="left" rowspan="1" colspan="1">(6.026)</td>
<td align="left" rowspan="1" colspan="1">(2.207)</td>
<td align="left" rowspan="1" colspan="1">(4.610)</td>
<td align="left" rowspan="1" colspan="1">(7.227)</td>
<td align="left" rowspan="1" colspan="1">(12.231)</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<table-wrap id="pone-0102877-t005" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.t005</object-id>
<label>Table 5</label>
<caption>
<title>Experimental Results for End-Effector Tracking: RMS Control Effort for the Control Algorithms.</title>
</caption>
<alternatives>
<graphic id="pone-0102877-t005-5" xlink:href="pone.0102877.t005"></graphic>
<table frame="hsides" rules="groups">
<colgroup span="1">
<col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead>
<tr>
<td align="left" rowspan="1" colspan="1">End-effector Tracking Results</td>
<td colspan="6" align="left" rowspan="1">Control Effort [mNm]</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">DataSet</td>
<td colspan="2" align="left" rowspan="1">Animal 1</td>
<td colspan="2" align="left" rowspan="1">Animal 2</td>
<td colspan="2" align="left" rowspan="1">Animal 3</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">Crystal Location</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
<td align="left" rowspan="1" colspan="1">Anterior</td>
<td align="left" rowspan="1" colspan="1">Lateral</td>
</tr>
</thead>
<tbody>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Exact Reference Information</td>
<td align="left" rowspan="1" colspan="1">29.284</td>
<td align="left" rowspan="1" colspan="1">38.286</td>
<td align="left" rowspan="1" colspan="1">28.734</td>
<td align="left" rowspan="1" colspan="1">32.217</td>
<td align="left" rowspan="1" colspan="1">48.327</td>
<td align="left" rowspan="1" colspan="1">43.298</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with One-Step Adaptive Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">34.416</td>
<td align="left" rowspan="1" colspan="1">51.048</td>
<td align="left" rowspan="1" colspan="1">29.813</td>
<td align="left" rowspan="1" colspan="1">31.574</td>
<td align="left" rowspan="1" colspan="1">54.779</td>
<td align="left" rowspan="1" colspan="1">43.241</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Generalized Adaptive Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">32.289</td>
<td align="left" rowspan="1" colspan="1">43.109</td>
<td align="left" rowspan="1" colspan="1">29.358</td>
<td align="left" rowspan="1" colspan="1">35.092</td>
<td align="left" rowspan="1" colspan="1">50.758</td>
<td align="left" rowspan="1" colspan="1">45.348</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">RHMPC with Extended Kalman Filter Estimation</td>
<td align="left" rowspan="1" colspan="1">35.092</td>
<td align="left" rowspan="1" colspan="1">48.341</td>
<td align="left" rowspan="1" colspan="1">32.947</td>
<td align="left" rowspan="1" colspan="1">32.657</td>
<td align="left" rowspan="1" colspan="1">67.176</td>
<td align="left" rowspan="1" colspan="1">55.756</td>
</tr>
<tr>
<td align="left" rowspan="1" colspan="1">PD Controller</td>
<td align="left" rowspan="1" colspan="1">49.938</td>
<td align="left" rowspan="1" colspan="1">83.809</td>
<td align="left" rowspan="1" colspan="1">42.894</td>
<td align="left" rowspan="1" colspan="1">57.258</td>
<td align="left" rowspan="1" colspan="1">140.301</td>
<td align="left" rowspan="1" colspan="1">99.273</td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>Tracking results of 183-s arrhythmia data from animal 1 for the generalized predictor are shown in
<xref ref-type="fig" rid="pone-0102877-g004">Fig. 4</xref>
for each PHANToM axis. Experimental results of the RHMPC with EKF predictor and PD controller are also given in
<xref ref-type="table" rid="pone-0102877-t004">Table 4</xref>
and
<xref ref-type="table" rid="pone-0102877-t005">Table 5</xref>
for comparison purposes (EKF predictor was implemented as described in
<xref rid="pone.0102877-Yuen1" ref-type="bibr">[4]</xref>
with the experimental parameters presented in
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
. More details on the PD controller can be found in
<xref rid="pone.0102877-Bebek2" ref-type="bibr">[31]</xref>
. Implementation details are omitted here due to space constraints). Tracking result for Axis 1 of PD controller is shown in
<xref ref-type="fig" rid="pone-0102877-g005">Fig. 5</xref>
.</p>
<fig id="pone-0102877-g004" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.g004</object-id>
<label>Figure 4</label>
<caption>
<title>Tracking results of 183-s arrhythmia data (only a part of the data is presented) from animal 1 for the generalized predictor.</title>
<p>Reference and PHANToM positions, RMS position error and MPC control effort are shown (A) Axis 1 results. (B) Axis 2 results. (C) Axis 3 results.</p>
</caption>
<graphic xlink:href="pone.0102877.g004"></graphic>
</fig>
<fig id="pone-0102877-g005" orientation="portrait" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0102877.g005</object-id>
<label>Figure 5</label>
<caption>
<title>Tracking results of 183-s arrhythmia data (only a part of the data is presented) from animal 1 for the PD Controller.</title>
<p>Axis 1 results are shown only.</p>
</caption>
<graphic xlink:href="pone.0102877.g005"></graphic>
</fig>
<p>The weighting parameters of the optimal index in RHMPC algorithms and the position and derivative constants in PD control were tuned to minimize RMS tracking error in both simulations and experiments. The tuning was performed individually for simulations and experiments. In experiments, additional tuning was performed to avoid the high frequency resonances so that no high frequency vibrations would be reflected to the structure.</p>
<p>In simulations, the generalized adaptive estimator gave better results than the exact heart signal in terms of RMS end-effector error (
<xref ref-type="table" rid="pone-0102877-t002">Table 2</xref>
). This is likely due to the combination of two factors. First, the simulation model is a linearized and reduced order model of the actual hardware, hence does not completely reflect the influence of the nonlinearities of the hardware on the tracking performance. Second, the adaptive estimator has a robustness characteristic that makes its output (i.e. predicted trajectory) less susceptible to the noise existent in the raw motion data
<xref rid="pone.0102877-Franke1" ref-type="bibr">[12]</xref>
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
. Accordingly, RHMPC with generalized estimator yields better results in the linear case. However, when the experiment is performed on the hardware, the effects of the nonlinearities have become apparent and the performance of the estimator-driven controller decreases (
<xref ref-type="table" rid="pone-0102877-t004">Table 4</xref>
). It should be noted that although the simulation provides valuable insight about the effectiveness of the controller, the experimental trials are the best indicator of performance.</p>
<p>The experimental results reveal that adaptive predictors provide satisfactory tracking performance for the arrhythmia trajectories considered (
<xref ref-type="table" rid="pone-0102877-t004">Table 4</xref>
). When the tracking results of the adaptive predictors are compared with each other, the generalized predictor gave better results than the one-step predictor in all experiments and yields close results to the the controller with exact heart signal reference. The RHMPC with adaptive estimation algorithms also outperformed RHMPC with EKF algorithm and PD controller.</p>
<p>In
<xref rid="pone.0102877-Yuen1" ref-type="bibr">[4]</xref>
, Yuen et al. used an EKF algorithm with a quasi-periodic motion model to predict the path of mitral valve motion in order to compensate the time delay resulting from the 3-D ultrasound (3DUS) measurements. They concluded that since the EKF explicitly models the quasi-periodic motion of the heart, it can adjust to normal variations in heart rhythm. In
<xref rid="pone.0102877-Yuen2" ref-type="bibr">[32]</xref>
, they noted that this does not extend to cases of arrhythmia, in which the motion of the heart is inherently non-periodic. Tracking results of RHMPC with EKF algorithm in
<xref ref-type="table" rid="pone-0102877-t004">Table 4</xref>
agree with this remark.</p>
<p>It is important to note that, RHMPC heavily relies on the estimate of the immediate future of the POI motion and the tracking problem is reduced to predicting the estimated reference heartbeat signal, when RHMPC has high enough precision to perform the necessary tracking
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
. Since EKF yielded poor predictions of the future POI position due to the unpredictable nature of heart during arrhythmia, the RHMPC with EKF algorithm was not able to track heart motion accurately. As a result, the PD controller yielded better performance than RHMPC with EKF algorithm.</p>
<p>The means of the RMS position errors and the standard deviations for RHMPC with generalized predictor, RHMPC with EKF, and PD controller are respectively computed as: 0.299 mm
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e007.jpg"></inline-graphic>
</inline-formula>
0.036 mm, 0.768 mm
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e008.jpg"></inline-graphic>
</inline-formula>
0.255 mm, 0.693 mm
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e009.jpg"></inline-graphic>
</inline-formula>
0.285 mm. These results further emphasize the better performance and robustness of RHMPC with generalized predictor on heart motion tracking during irregular beating. The effects of using RHMPC with generalized predictor on the arrhythmia tracking were also tested for statistical significance by paired t-test. For the results presented in
<xref ref-type="table" rid="pone-0102877-t004">Table 4</xref>
, RHMPC with generalized predictor has led to improvements in tracking when compared to RHMPC with EKF (P
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e010.jpg"></inline-graphic>
</inline-formula>
0.004) and PD controller (P
<inline-formula>
<inline-graphic xlink:href="pone.0102877.e011.jpg"></inline-graphic>
</inline-formula>
0.009).</p>
<p>We conjecture that the maximum error values are affected from the noise in the sensor data. In particular, the sonomicrometer measurements are susceptible to non-Gaussian noise resulting form echo effects and interruptions in ultrasound transmission paths that appear as large spikes or jump discontinuities. It is unlikely that the POI on the heart is capable of moving as much as 12 mm in a few milliseconds as it can occasionally be observed in the experimental data. More aggressive filtering can be performed to eliminate such high frequency motions (e.g.
<xref rid="pone.0102877-Liu2" ref-type="bibr">[33]</xref>
<xref rid="pone.0102877-Menold1" ref-type="bibr">[35]</xref>
), but the data has been kept as-is without applying any filtering to eliminate these jumps, as currently we do not have an independent set of sensor measurements (such as from a vision sensor) that would validate this conjecture. Though, it should be noted that such noise occurs only occasionally, as it can be observed through the significantly lower RMS error values, relative to maximum error values. At these high jumps, RHMPC algorithms with adaptive predictors outperformed PD controller, because with the PD controller these jumps initiated system oscillations (
<xref ref-type="fig" rid="pone-0102877-g005">Fig. 5</xref>
). Also, when the RMS control efforts are compared, the PD controller performed poorly (
<xref ref-type="table" rid="pone-0102877-t005">Table 5</xref>
).</p>
</sec>
<sec id="s4">
<title>Discussions</title>
<p>In this paper, heart motion tracking with adaptive estimation algorithms in the presence of arrhythmia is presented. Performance of the algorithms are evaluated with a range of data. These algorithms were previously explored under slowly varying heart rate conditions in
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
. To the best of our knowledge, no previous POI motion tracking work (see
<xref ref-type="sec" rid="s1">
<italic>Introduction</italic>
</xref>
) has presented results in the presence of arrhythmia. The experimental RMS errors on the order of 0.250–0.340 mm obtained using the generalized estimator (in comparison to the RMS tracking errors on the order of 0.170–0.350 mm in constant heart rate and 0.160–0.180 mm in slowly varying heart rate conditions, as reported in
<xref rid="pone.0102877-Tuna1" ref-type="bibr">[14]</xref>
) represent satisfactory tracking performance during arrhythmia. Results show that if the behavior of the heart changes abruptly, the predictors are able to adapt the new heart behavior and can track the ideal time-varying solution. Although the resulting tracking errors are above the desired specification for effectively performing anastomosis, the system still achieved sufficient tracking accuracy to maintain safe tracking, until the system can potentially switch to a safe mode of operation. The results of the study should also be validated
<italic>in vivo</italic>
. An
<italic>in vivo</italic>
validation study would be valuable to verify that the proposed scheme will be effective under practical constraints of an operating room setting.</p>
<p>It is important to note that a certain level of consistency was necessary during the experimental setting of this study. For this reason, atrial fibrillation and its electrocardiographic pattern have been chosen as an arrhythmic scenario to provide a comparable and repeatable type of cardiac rhythm disturbance. The validity of the presented approach on other types of the peri-operative rhythm disturbance like tachycardia or sinus bradycardia will be explored in future studies. In addition, atrial fibrillation is a complex type of arrhythmia with high ventricular rates. Accordingly, the baseline heart rate (HR) is an important variable that was continuously monitored throughout the study. Yet, since the HR values may vary in each and every given setting, the high variability of the baseline heart rate was not considered to be a drawback of the experimental approach in choosing a bovine model. Similar to humans, the HR may vary per baseline and intra-operative hemodynamic conditions, but no interference and/or vagal tone dominance was observed during the studies. However, if the native cardiac pacing is altered, then the vagal tone and correlated effects may become more significant in heart failure model.</p>
<p>We also note that the presented study is not intended to evaluate the tracking performance of the algorithms on each and every type of atrio-ventricular arrhythmia. Rather, by simulating one particular type of arrhythmia from a wide variety of patterns, we aim to show feasibility of employing adaptive estimation algorithms for heart motion tracking in the occurrence of arrhythmia during beating heart surgery. The underlying electro-physiological mechanisms of arrhythmia are very complex. Atrio-ventricular contraction is affected by multitude of factors, such as, changes in the function of ion channels in atrial cells. The complex pattern of electrical impulse propagation within atrio-ventricular tissues is also determined by the frequency of stimulation and summation of waveforms. As such, constructing
<italic>in vivo</italic>
models of arbitrary arrhythmia types is not practical, if not impossible. A more broad
<italic>in vivo</italic>
study would overcome the limitations of the employed method of artificially inducing arrhythmias, and would cover more diverse situations in CABG surgeries and so will evaluate the performance and robustness of the proposed method in a wide range of conditions. Additionally, the applicability of the proposed system to specific arrhythmia patterns would be considered in further studies as a part of the development and safety control evaluation.</p>
<p>One way to improve tracking quality is to incorporate other types of data into the estimation scheme. One such possibility is to include the electrocardiogram (ECG) signal into the observations. In this way, the predictor is able to use the electrical signals that activate heart contraction in order to improve the prediction as in
<xref rid="pone.0102877-Bebek1" ref-type="bibr">[11]</xref>
. The ECG signal is very suitable for period-to-period synchronization with sufficient lead time for the feed-forward control, and identification of arrhythmias. If an arrhythmia is detected, feed-forward component of the controller can be turned off and switched to a further fail-safe mode if necessary, which in turn may improve performance during heart rhythm abnormalities.</p>
<p>Future works will also include multi sensor fusion where complementary and redundant sensors will be used for superior performance and safety, e.g., a vision based sensor system could be used as a secondary sensor for the
<italic>in vivo</italic>
validation of the proposed concept. Merging the sensor data from multiple position sources would increase accuracy of the motion estimation and improve tracking results. Adding more mechanical sensors that measure heart motion would improve the measurement precision and help to resolve calibration problems of the sonomicrometry system.</p>
</sec>
</body>
<back>
<ack>
<p>We would like to thank Dr. Akira Shiose for his help during collection of the experimental data.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="pone.0102877-Ginhoux1">
<label>1</label>
<mixed-citation publication-type="journal">
<name>
<surname>Ginhoux</surname>
<given-names>R</given-names>
</name>
,
<name>
<surname>Gangloff</surname>
<given-names>JA</given-names>
</name>
,
<name>
<surname>DeMathelin</surname>
<given-names>MF</given-names>
</name>
,
<name>
<surname>Soler</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Leroy</surname>
<given-names>J</given-names>
</name>
,
<etal>et al</etal>
(
<year>2005</year>
)
<article-title>Active filtering of physiological motion in robotized surgery using predictive control</article-title>
.
<source>IEEE Trans Robot</source>
<volume>21</volume>
:
<fpage>67</fpage>
<lpage>79</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Richa1">
<label>2</label>
<mixed-citation publication-type="other">Richa R, Bo APL, Poignet P (2008) Motion prediction for tracking the beating heart. In: Proc. of the Annual International Conference of the Engineering in Medicine and Biology Society (EMBC). Vancouver, British Columbia, Canada, pp. 3261–3264.</mixed-citation>
</ref>
<ref id="pone.0102877-Richa2">
<label>3</label>
<mixed-citation publication-type="journal">
<name>
<surname>Richa</surname>
<given-names>R</given-names>
</name>
,
<name>
<surname></surname>
<given-names>AP</given-names>
</name>
,
<name>
<surname>Poignet</surname>
<given-names>P</given-names>
</name>
(
<year>2011</year>
)
<article-title>Towards robust 3d visual tracking for motion compensation in beating heart surgery</article-title>
.
<source>Medical Image Analysis</source>
<volume>15</volume>
:
<fpage>302</fpage>
<lpage>315</lpage>
<pub-id pub-id-type="pmid">21277821</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Yuen1">
<label>4</label>
<mixed-citation publication-type="journal">
<name>
<surname>Yuen</surname>
<given-names>SG</given-names>
</name>
,
<name>
<surname>Kettler</surname>
<given-names>DT</given-names>
</name>
,
<name>
<surname>Novotny</surname>
<given-names>PM</given-names>
</name>
,
<name>
<surname>Plowes</surname>
<given-names>RD</given-names>
</name>
,
<name>
<surname>Howe</surname>
<given-names>RD</given-names>
</name>
(
<year>2009</year>
)
<article-title>Robotic motion compensation for beating heart intracardiac surgery</article-title>
.
<source>International Journal of Robotics Research</source>
<volume>28</volume>
:
<fpage>1355</fpage>
<lpage>1372</lpage>
<pub-id pub-id-type="pmid">20436927</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Bachta1">
<label>5</label>
<mixed-citation publication-type="journal">
<name>
<surname>Bachta</surname>
<given-names>W</given-names>
</name>
,
<name>
<surname>Renaud</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Cuvillon</surname>
<given-names>L</given-names>
</name>
,
<name>
<surname>Laroche</surname>
<given-names>E</given-names>
</name>
,
<name>
<surname>Forgione</surname>
<given-names>A</given-names>
</name>
,
<etal>et al</etal>
(
<year>2009</year>
)
<article-title>Motion prediction for computer-assisted beating heart surgery</article-title>
.
<source>IEEE Trans Biomed Eng</source>
<volume>56</volume>
:
<fpage>2551</fpage>
<lpage>63</lpage>
<pub-id pub-id-type="pmid">19567337</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Cuvillon1">
<label>6</label>
<mixed-citation publication-type="other">Cuvillon L, Gangloff J, DeMathelin M, Forgione A (2005) Toward robotized beating heart TECABG: Assessment of the heart dynamics using high-speed vision. In: Proc. of Medical Image Computing and Computer-Assisted Intervention (MICCAI). Palm Springs, USA, volume 2, pp. 551–558.</mixed-citation>
</ref>
<ref id="pone.0102877-Bader1">
<label>7</label>
<mixed-citation publication-type="other">Bader T, Wiedemann A, Roberts K, Hanebeck UD (2007) Model-based motion estimation of elastic surfaces for minimally invasive cardiac surgery. In: Proc. of the International Conference on Robotics and Automation (ICRA). Rome, Italy, pp. 2261–2266.</mixed-citation>
</ref>
<ref id="pone.0102877-Bogatyrenko1">
<label>8</label>
<mixed-citation publication-type="other">Bogatyrenko E, Hanebeck U, Szabo G (2009) Heart surface motion estimation framework for robotic surgery employing meshless methods. In: Intelligent Robots and Systems, 2009. IROS 2009. IEEE/RSJ International Conference on. pp. 67–74.</mixed-citation>
</ref>
<ref id="pone.0102877-Bogatyrenko2">
<label>9</label>
<mixed-citation publication-type="journal">
<name>
<surname>Bogatyrenko</surname>
<given-names>E</given-names>
</name>
,
<name>
<surname>Pompey</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Hanebeck</surname>
<given-names>U</given-names>
</name>
(
<year>2011</year>
)
<article-title>Efficient physics-based tracking of heart surface motion for beating heart surgery robotic systems</article-title>
.
<source>International Journal of Computer Assisted Radiology and Surgery</source>
<volume>6</volume>
:
<fpage>387</fpage>
<lpage>399</lpage>
<pub-id pub-id-type="pmid">20694522</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Ortmaier1">
<label>10</label>
<mixed-citation publication-type="journal">
<name>
<surname>Ortmaier</surname>
<given-names>T</given-names>
</name>
,
<name>
<surname>Groeger</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Boehm</surname>
<given-names>DH</given-names>
</name>
,
<name>
<surname>Falk</surname>
<given-names>V</given-names>
</name>
,
<name>
<surname>Hirzinger</surname>
<given-names>G</given-names>
</name>
(
<year>2005</year>
)
<article-title>Motion estimation in beating heart surgery</article-title>
.
<source>IEEE Trans Biomed Eng</source>
<volume>52</volume>
:
<fpage>1729</fpage>
<lpage>1740</lpage>
<pub-id pub-id-type="pmid">16235658</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Bebek1">
<label>11</label>
<mixed-citation publication-type="journal">
<name>
<surname>Bebek</surname>
<given-names>O</given-names>
</name>
,
<name>
<surname>Cavusoglu</surname>
<given-names>MC</given-names>
</name>
(
<year>2007</year>
)
<article-title>Intelligent control algorithms for robotic-assisted beating heart surgery</article-title>
.
<source>IEEE Trans Robot</source>
<volume>23</volume>
:
<fpage>468</fpage>
<lpage>480</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Franke1">
<label>12</label>
<mixed-citation publication-type="other">Franke T, Bebek O, Cavusoglu MC (2007) Improved prediction of heart motion using an adaptive filter for robot assisted beating heart surgery. In: Proc. of IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS). San Diego, CA, USA, pp. 509–515.</mixed-citation>
</ref>
<ref id="pone.0102877-Franke2">
<label>13</label>
<mixed-citation publication-type="other">Franke T, Bebek O, Cavusoglu MC (2008) Prediction of heartbeat motion with a generalized adaptive filter. In: Proc. of the International Conference on Robotics and Automation (ICRA). Pasadena, CA, USA, pp. 2916–2921.</mixed-citation>
</ref>
<ref id="pone.0102877-Tuna1">
<label>14</label>
<mixed-citation publication-type="journal">
<name>
<surname>Tuna</surname>
<given-names>EE</given-names>
</name>
,
<name>
<surname>Franke</surname>
<given-names>TJ</given-names>
</name>
,
<name>
<surname>Bebek</surname>
<given-names>O</given-names>
</name>
,
<name>
<surname>Shiose</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Fukamachi</surname>
<given-names>K</given-names>
</name>
,
<etal>et al</etal>
(
<year>2013</year>
)
<article-title>Heart motion prediction based on adaptive estimation algorithms for robotic-assisted beating heart surgery</article-title>
.
<source>IEEE Trans Robot</source>
<volume>29</volume>
:
<fpage>261</fpage>
<lpage>276</lpage>
<pub-id pub-id-type="pmid">23976889</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Trejos1">
<label>15</label>
<mixed-citation publication-type="other">Trejos AL, Salcudean SE, Sassani F, Lichtenstein S (1999) On the feasibility of a moving support for surgery on the beating heart. In: Proc. of Medical Image Computing and Computer-Assisted Interventions (MICCAI). Cambridge, UK, pp. 1088–1097.</mixed-citation>
</ref>
<ref id="pone.0102877-Newman1">
<label>16</label>
<mixed-citation publication-type="journal">
<name>
<surname>Newman</surname>
<given-names>MF</given-names>
</name>
,
<name>
<surname>Kirchner</surname>
<given-names>JL</given-names>
</name>
,
<name>
<surname>Phillips-Bute</surname>
<given-names>B</given-names>
</name>
,
<name>
<surname>Gaver</surname>
<given-names>V</given-names>
</name>
,
<name>
<surname>Grocott</surname>
<given-names>H</given-names>
</name>
,
<etal>et al</etal>
(
<year>2001</year>
)
<article-title>Longitudinal assessment of neurocognitive function after coronary-artery bypass surgery</article-title>
.
<source>New England Journal of Medicine</source>
<volume>344</volume>
:
<fpage>395</fpage>
<lpage>402</lpage>
<pub-id pub-id-type="pmid">11172175</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Lemma1">
<label>17</label>
<mixed-citation publication-type="other">Lemma M, Mangini A, Redaelli A, Acocella F (2005) Do cardiac stabilizers really stabilize? experimental quantitative analysis of mechanical stabilization. Interactive CardioVascular and Thoracic Surgery: 222–226.</mixed-citation>
</ref>
<ref id="pone.0102877-Nakamura1">
<label>18</label>
<mixed-citation publication-type="other">Nakamura Y, Kishi K, Kawakami H (2001) Heartbeat synchronization for robotic cardiac surgery. In: Proc. of IEEE International Conference on Robotics and Automation (ICRA). Seoul, Korea, volume 2, pp. 2014–2019.</mixed-citation>
</ref>
<ref id="pone.0102877-Roka1">
<label>19</label>
<mixed-citation publication-type="journal">
<name>
<surname>Roka</surname>
<given-names>A</given-names>
</name>
,
<name>
<surname>Toth</surname>
<given-names>E</given-names>
</name>
,
<name>
<surname>Szilagyi</surname>
<given-names>S</given-names>
</name>
,
<name>
<surname>Merkely</surname>
<given-names>B</given-names>
</name>
(
<year>2008</year>
)
<article-title>Electrical atrial fibrillation induction affects the characteristics of induced arrhythmia</article-title>
.
<source>Journal of Electrocardiology</source>
<volume>41</volume>
:
<fpage>131</fpage>
<lpage>137</lpage>
<pub-id pub-id-type="pmid">17631892</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Takagaki1">
<label>20</label>
<mixed-citation publication-type="journal">
<name>
<surname>Takagaki</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>McCarthy</surname>
<given-names>PM</given-names>
</name>
,
<name>
<surname>Chung</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Connor</surname>
<given-names>J</given-names>
</name>
,
<name>
<surname>Dessoffy</surname>
<given-names>R</given-names>
</name>
,
<etal>et al</etal>
(
<year>2002</year>
)
<article-title>Preload-adjusted maximal power: a novel index of left ventricular contractility in atrial fibrillation</article-title>
.
<source>Heart</source>
<volume>88</volume>
:
<fpage>170</fpage>
<lpage>176</lpage>
<pub-id pub-id-type="pmid">12117848</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Takagaki2">
<label>21</label>
<mixed-citation publication-type="journal">
<name>
<surname>Takagaki</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>McCarthy</surname>
<given-names>P</given-names>
</name>
,
<name>
<surname>Inoue</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Chung</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Connor</surname>
<given-names>J</given-names>
</name>
,
<etal>et al</etal>
(
<year>2005</year>
)
<article-title>Myocardial compliance was not altered after acute induction of atrial fibrillation in sheep</article-title>
.
<source>Medical Science Monitor</source>
<volume>41</volume>
:
<fpage>147</fpage>
<lpage>153</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Clarke1">
<label>22</label>
<mixed-citation publication-type="journal">
<name>
<surname>Clarke</surname>
<given-names>DW</given-names>
</name>
(
<year>1988</year>
)
<article-title>Application of generalized predictive control to industrial processes</article-title>
.
<source>IEEE Control Syst Mag</source>
<volume>8</volume>
:
<fpage>49</fpage>
<lpage>55</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Garcia1">
<label>23</label>
<mixed-citation publication-type="journal">
<name>
<surname>Garcia</surname>
<given-names>CE</given-names>
</name>
,
<name>
<surname>Prett</surname>
<given-names>DM</given-names>
</name>
,
<name>
<surname>Morari</surname>
<given-names>M</given-names>
</name>
(
<year>1989</year>
)
<article-title>Model predictive control: Theory and practice-a survey</article-title>
.
<source>Automatica</source>
<volume>25</volume>
:
<fpage>335</fpage>
<lpage>348</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Anderson1">
<label>24</label>
<mixed-citation publication-type="other">Anderson BD, Moore JB (1990) Optimal Control Linear Quadratic Methods. Prentice Hall.</mixed-citation>
</ref>
<ref id="pone.0102877-Camacho1">
<label>25</label>
<mixed-citation publication-type="other">Camacho E, Bordons C (1999) Model Predictive Control. Springer.</mixed-citation>
</ref>
<ref id="pone.0102877-Liu1">
<label>26</label>
<mixed-citation publication-type="other">Liu T (2011) Design and Prototyping of a Three Degrees of Freedom Robotic Wrist Mechanism for a Robotic Surgery System. M.s. thesis, Case Western Reserve University, Cleveland, OH, USA. Available:
<ext-link ext-link-type="uri" xlink:href="http://etd.ohiolink.edu/view.cgi?acc_num=case1283538593">http://etd.ohiolink.edu/view.cgi?acc_num=case1283538593</ext-link>
Advisor: M. Cenk Cavusoglu.</mixed-citation>
</ref>
<ref id="pone.0102877-Cavusoglu1">
<label>27</label>
<mixed-citation publication-type="journal">
<name>
<surname>Cavusoglu</surname>
<given-names>MC</given-names>
</name>
,
<name>
<surname>Feygin</surname>
<given-names>D</given-names>
</name>
,
<name>
<surname>Tendick</surname>
<given-names>F</given-names>
</name>
(
<year>2002</year>
)
<article-title>A critical study of the mechanical and electrical properties of the phantom haptic interface and improvements for highperformance control</article-title>
.
<source>Presence: Teleoperators and Virtual Environments</source>
<volume>11</volume>
:
<fpage>555</fpage>
<lpage>568</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Dione1">
<label>28</label>
<mixed-citation publication-type="other">Dione D, Shi P, Smith W, DeMan P, Soares J,
<etal>et al</etal>
. (1997) Three-dimensional regional left ventricular deformation from digital sonomicrometry. In: Proc. of the Annual International Conference of the Engineering in Medicine and Biology Society (EMBC). pp. 848–851.</mixed-citation>
</ref>
<ref id="pone.0102877-Ratcliffe1">
<label>29</label>
<mixed-citation publication-type="journal">
<name>
<surname>Ratcliffe</surname>
<given-names>MB</given-names>
</name>
,
<name>
<surname>Gupta</surname>
<given-names>KB</given-names>
</name>
,
<name>
<surname>Streicher</surname>
<given-names>JT</given-names>
</name>
,
<name>
<surname>Savage</surname>
<given-names>EB</given-names>
</name>
,
<name>
<surname>Bogen</surname>
<given-names>DK</given-names>
</name>
,
<etal>et al</etal>
(
<year>1995</year>
)
<article-title>Use of sonomicrometry and multidimensional scaling to determine the three-dimensional coordinates of multiple cardiac locations: Feasibility and initial implementation</article-title>
.
<source>IEEE Trans Biomed Eng</source>
<volume>42</volume>
:
<fpage>587</fpage>
<lpage>598</lpage>
<pub-id pub-id-type="pmid">7790015</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0102877-Bowthorpe1">
<label>30</label>
<mixed-citation publication-type="journal">
<name>
<surname>Bowthorpe</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Tavakoli</surname>
<given-names>M</given-names>
</name>
,
<name>
<surname>Becher</surname>
<given-names>H</given-names>
</name>
,
<name>
<surname>Howe</surname>
<given-names>R</given-names>
</name>
(
<year>2014</year>
)
<article-title>Smith predictor-based robot control for ultrasound-guided teleoperated beating-heart surgery</article-title>
.
<source>Biomedical and Health Informatics, IEEE Journal of</source>
<volume>18</volume>
:
<fpage>157</fpage>
<lpage>166</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Bebek2">
<label>31</label>
<mixed-citation publication-type="other">Bebek O (2008) Robotic-Assisted Beating Heart Surgery. Ph.D. Dissertation, Case Western Reserve University, Cleveland, OH, USA. Available:
<ext-link ext-link-type="uri" xlink:href="http://rave.ohiolink.edu/etdc/view?acc_num=case1201289393">http://rave.ohiolink.edu/etdc/view?acc_num=case1201289393</ext-link>
Advisor: M. Cenk Cavusoglu.</mixed-citation>
</ref>
<ref id="pone.0102877-Yuen2">
<label>32</label>
<mixed-citation publication-type="other">Yuen SG, Vasilyev NV, del Nido PJ, Howe RD (2010) Robotic tissue tracking for beating heart mitral valve surgery. Medical Image Analysis.</mixed-citation>
</ref>
<ref id="pone.0102877-Liu2">
<label>33</label>
<mixed-citation publication-type="journal">
<name>
<surname>Liu</surname>
<given-names>H</given-names>
</name>
,
<name>
<surname>Shah</surname>
<given-names>S</given-names>
</name>
(
<year>2004</year>
)
<collab>JiangW</collab>
(
<year>2004</year>
)
<article-title>Online outlier detection and data cleaning</article-title>
.
<source>Computers and Chemical Engineering</source>
<volume>28</volume>
:
<fpage>1635</fpage>
<lpage>1647</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Pearson1">
<label>34</label>
<mixed-citation publication-type="journal">
<name>
<surname>Pearson</surname>
<given-names>RK</given-names>
</name>
(
<year>2002</year>
)
<article-title>Outliers in process modeling and identification</article-title>
.
<source>IEEE Trans Contr Sys Techn</source>
<volume>10</volume>
:
<fpage>55</fpage>
<lpage>63</lpage>
</mixed-citation>
</ref>
<ref id="pone.0102877-Menold1">
<label>35</label>
<mixed-citation publication-type="other">Menold P, Pearson R, Allgöwere F (1999) Online outlier detection and removal. In: Proc. of the 7th Mediterranean Conference on Control and Automation. pp. 1110–114.</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
<affiliations>
<list>
<country>
<li>Turquie</li>
<li>États-Unis</li>
</country>
<region>
<li>Ohio</li>
</region>
</list>
<tree>
<country name="États-Unis">
<region name="Ohio">
<name sortKey="Tuna, E Erdem" sort="Tuna, E Erdem" uniqKey="Tuna E" first="E. Erdem" last="Tuna">E. Erdem Tuna</name>
</region>
<name sortKey="Cavu O Lu, M Cenk" sort="Cavu O Lu, M Cenk" uniqKey="Cavu O Lu M" first="M. Cenk" last="Çavu O Lu">M. Cenk Çavu O Lu</name>
<name sortKey="Fukamachi, Kiyotaka" sort="Fukamachi, Kiyotaka" uniqKey="Fukamachi K" first="Kiyotaka" last="Fukamachi">Kiyotaka Fukamachi</name>
<name sortKey="Karimov, Jamshid H" sort="Karimov, Jamshid H" uniqKey="Karimov J" first="Jamshid H." last="Karimov">Jamshid H. Karimov</name>
<name sortKey="Liu, Taoming" sort="Liu, Taoming" uniqKey="Liu T" first="Taoming" last="Liu">Taoming Liu</name>
</country>
<country name="Turquie">
<noRegion>
<name sortKey="Bebek, Ozkan" sort="Bebek, Ozkan" uniqKey="Bebek O" first="Özkan" last="Bebek">Özkan Bebek</name>
</noRegion>
</country>
</tree>
</affiliations>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/HapticV1/Data/Pmc/Checkpoint
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000856 | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Checkpoint/biblio.hfd -nk 000856 | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Ticri/CIDE
   |area=    HapticV1
   |flux=    Pmc
   |étape=   Checkpoint
   |type=    RBID
   |clé=     PMC:4105597
   |texte=   Towards Active Tracking of Beating Heart Motion in the Presence of Arrhythmia for Robotic Assisted Beating Heart Surgery
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Checkpoint/RBID.i   -Sk "pubmed:25048462" \
       | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Checkpoint/biblio.hfd   \
       | NlmPubMed2Wicri -a HapticV1 

Wicri

This area was generated with Dilib version V0.6.23.
Data generation: Mon Jun 13 01:09:46 2016. Site generation: Wed Mar 6 09:54:07 2024