Locally rigid, vessel-based registration for laparoscopic liver surgery
Identifieur interne : 000819 ( Pmc/Curation ); précédent : 000818; suivant : 000820Locally rigid, vessel-based registration for laparoscopic liver surgery
Auteurs : Yi Song [Royaume-Uni] ; Johannes Totz [Royaume-Uni] ; Steve Thompson [Royaume-Uni] ; Stian Johnsen [Royaume-Uni] ; Dean Barratt [Royaume-Uni] ; Crispin Schneider [Royaume-Uni] ; Kurinchi Gurusamy [Royaume-Uni] ; Brian Davidson [Royaume-Uni] ; Sébastien Ourselin [Royaume-Uni] ; David Hawkes [Royaume-Uni] ; Matthew J. Clarkson [Royaume-Uni]Source :
- International Journal of Computer Assisted Radiology and Surgery [ 1861-6410 ] ; 2015.
Abstract
Laparoscopic liver resection has significant advantages over open surgery due to less patient trauma and faster recovery times, yet is difficult for most lesions due to the restricted field of view and lack of haptic feedback. Image guidance provides a potential solution but is challenging in a soft deforming organ such as the liver. In this paper, we therefore propose a laparoscopic ultrasound (LUS) image guidance system and study the feasibility of a locally rigid registration for laparoscopic liver surgery.
We developed a real-time segmentation method to extract vessel centre points from calibrated, freehand, electromagnetically tracked, 2D LUS images. Using landmark-based initial registration and an optional iterative closest point (ICP) point-to-line registration, a vessel centre-line model extracted from preoperative computed tomography (CT) is registered to the ultrasound data during surgery.
Using the locally rigid ICP method, the RMS residual error when registering to a phantom was 0.7 mm, and the mean target registration error (TRE) for two in vivo porcine studies was 3.58 and 2.99 mm, respectively. Using the locally rigid landmark-based registration method gave a mean TRE of 4.23 mm using vessel centre lines derived from CT scans taken with pneumoperitoneum and 6.57 mm without pneumoperitoneum.
In this paper we propose a practical image-guided surgery system based on locally rigid registration of a CT-derived model to vascular structures located with LUS. In a physical phantom and during porcine laparoscopic liver resection, we demonstrate accuracy of target location commensurate with surgical requirements. We conclude that locally rigid registration could be sufficient for practically useful image guidance in the near future.
Url:
DOI: 10.1007/s11548-015-1236-8
PubMed: 26092658
PubMed Central: 4642598
Links toward previous steps (curation, corpus...)
- to stream Pmc, to step Corpus: Pour aller vers cette notice dans l'étape Curation :000819
Links to Exploration step
PMC:4642598Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Locally rigid, vessel-based registration for laparoscopic liver surgery</title>
<author><name sortKey="Song, Yi" sort="Song, Yi" uniqKey="Song Y" first="Yi" last="Song">Yi Song</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Totz, Johannes" sort="Totz, Johannes" uniqKey="Totz J" first="Johannes" last="Totz">Johannes Totz</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Thompson, Steve" sort="Thompson, Steve" uniqKey="Thompson S" first="Steve" last="Thompson">Steve Thompson</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Johnsen, Stian" sort="Johnsen, Stian" uniqKey="Johnsen S" first="Stian" last="Johnsen">Stian Johnsen</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Barratt, Dean" sort="Barratt, Dean" uniqKey="Barratt D" first="Dean" last="Barratt">Dean Barratt</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Schneider, Crispin" sort="Schneider, Crispin" uniqKey="Schneider C" first="Crispin" last="Schneider">Crispin Schneider</name>
<affiliation wicri:level="1"><nlm:aff id="Aff2">Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Gurusamy, Kurinchi" sort="Gurusamy, Kurinchi" uniqKey="Gurusamy K" first="Kurinchi" last="Gurusamy">Kurinchi Gurusamy</name>
<affiliation wicri:level="1"><nlm:aff id="Aff2">Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Davidson, Brian" sort="Davidson, Brian" uniqKey="Davidson B" first="Brian" last="Davidson">Brian Davidson</name>
<affiliation wicri:level="1"><nlm:aff id="Aff2">Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Ourselin, Sebastien" sort="Ourselin, Sebastien" uniqKey="Ourselin S" first="Sébastien" last="Ourselin">Sébastien Ourselin</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Hawkes, David" sort="Hawkes, David" uniqKey="Hawkes D" first="David" last="Hawkes">David Hawkes</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Clarkson, Matthew J" sort="Clarkson, Matthew J" uniqKey="Clarkson M" first="Matthew J." last="Clarkson">Matthew J. Clarkson</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">26092658</idno>
<idno type="pmc">4642598</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4642598</idno>
<idno type="RBID">PMC:4642598</idno>
<idno type="doi">10.1007/s11548-015-1236-8</idno>
<date when="2015">2015</date>
<idno type="wicri:Area/Pmc/Corpus">000819</idno>
<idno type="wicri:Area/Pmc/Curation">000819</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Locally rigid, vessel-based registration for laparoscopic liver surgery</title>
<author><name sortKey="Song, Yi" sort="Song, Yi" uniqKey="Song Y" first="Yi" last="Song">Yi Song</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Totz, Johannes" sort="Totz, Johannes" uniqKey="Totz J" first="Johannes" last="Totz">Johannes Totz</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Thompson, Steve" sort="Thompson, Steve" uniqKey="Thompson S" first="Steve" last="Thompson">Steve Thompson</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Johnsen, Stian" sort="Johnsen, Stian" uniqKey="Johnsen S" first="Stian" last="Johnsen">Stian Johnsen</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Barratt, Dean" sort="Barratt, Dean" uniqKey="Barratt D" first="Dean" last="Barratt">Dean Barratt</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Schneider, Crispin" sort="Schneider, Crispin" uniqKey="Schneider C" first="Crispin" last="Schneider">Crispin Schneider</name>
<affiliation wicri:level="1"><nlm:aff id="Aff2">Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Gurusamy, Kurinchi" sort="Gurusamy, Kurinchi" uniqKey="Gurusamy K" first="Kurinchi" last="Gurusamy">Kurinchi Gurusamy</name>
<affiliation wicri:level="1"><nlm:aff id="Aff2">Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Davidson, Brian" sort="Davidson, Brian" uniqKey="Davidson B" first="Brian" last="Davidson">Brian Davidson</name>
<affiliation wicri:level="1"><nlm:aff id="Aff2">Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Ourselin, Sebastien" sort="Ourselin, Sebastien" uniqKey="Ourselin S" first="Sébastien" last="Ourselin">Sébastien Ourselin</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Hawkes, David" sort="Hawkes, David" uniqKey="Hawkes D" first="David" last="Hawkes">David Hawkes</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
<author><name sortKey="Clarkson, Matthew J" sort="Clarkson, Matthew J" uniqKey="Clarkson M" first="Matthew J." last="Clarkson">Matthew J. Clarkson</name>
<affiliation wicri:level="1"><nlm:aff id="Aff1">Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</nlm:aff>
<country xml:lang="fr">Royaume-Uni</country>
<wicri:regionArea>Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London</wicri:regionArea>
</affiliation>
</author>
</analytic>
<series><title level="j">International Journal of Computer Assisted Radiology and Surgery</title>
<idno type="ISSN">1861-6410</idno>
<idno type="eISSN">1861-6429</idno>
<imprint><date when="2015">2015</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><sec><title>Purpose</title>
<p>Laparoscopic liver resection has significant advantages over open surgery due to less patient trauma and faster recovery times, yet is difficult for most lesions due to the restricted field of view and lack of haptic feedback. Image guidance provides a potential solution but is challenging in a soft deforming organ such as the liver. In this paper, we therefore propose a laparoscopic ultrasound (LUS) image guidance system and study the feasibility of a locally rigid registration for laparoscopic liver surgery.</p>
</sec>
<sec><title>Methods</title>
<p>We developed a real-time segmentation method to extract vessel centre points from calibrated, freehand, electromagnetically tracked, 2D LUS images. Using landmark-based initial registration and an optional iterative closest point (ICP) point-to-line registration, a vessel centre-line model extracted from preoperative computed tomography (CT) is registered to the ultrasound data during surgery.</p>
</sec>
<sec><title>Results</title>
<p>Using the locally rigid ICP method, the RMS residual error when registering to a phantom was 0.7 mm, and the mean target registration error (TRE) for two in vivo porcine studies was 3.58 and 2.99 mm, respectively. Using the locally rigid landmark-based registration method gave a mean TRE of 4.23 mm using vessel centre lines derived from CT scans taken with pneumoperitoneum and 6.57 mm without pneumoperitoneum.</p>
</sec>
<sec><title>Conclusion</title>
<p>In this paper we propose a practical image-guided surgery system based on locally rigid registration of a CT-derived model to vascular structures located with LUS. In a physical phantom and during porcine laparoscopic liver resection, we demonstrate accuracy of target location commensurate with surgical requirements. We conclude that locally rigid registration could be sufficient for practically useful image guidance in the near future.</p>
</sec>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Barratt, Dc" uniqKey="Barratt D">DC Barratt</name>
</author>
<author><name sortKey="Davies, Ah" uniqKey="Davies A">AH Davies</name>
</author>
<author><name sortKey="Hughes, Ad" uniqKey="Hughes A">AD Hughes</name>
</author>
<author><name sortKey="Thom, Sa" uniqKey="Thom S">SA Thom</name>
</author>
<author><name sortKey="Humphries, Kn" uniqKey="Humphries K">KN Humphries</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Clarkson, M" uniqKey="Clarkson M">M Clarkson</name>
</author>
<author><name sortKey="Zombori, G" uniqKey="Zombori G">G Zombori</name>
</author>
<author><name sortKey="Thompson, S" uniqKey="Thompson S">S Thompson</name>
</author>
<author><name sortKey="Totz, J" uniqKey="Totz J">J Totz</name>
</author>
<author><name sortKey="Song, Y" uniqKey="Song Y">Y Song</name>
</author>
<author><name sortKey="Espak, M" uniqKey="Espak M">M Espak</name>
</author>
<author><name sortKey="Johnsen, S" uniqKey="Johnsen S">S Johnsen</name>
</author>
<author><name sortKey="Hawkes, D" uniqKey="Hawkes D">D Hawkes</name>
</author>
<author><name sortKey="Ourselin, S" uniqKey="Ourselin S">S Ourselin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Croome, Kp" uniqKey="Croome K">KP Croome</name>
</author>
<author><name sortKey="Yamashita, Mh" uniqKey="Yamashita M">MH Yamashita</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Feuerstein, M" uniqKey="Feuerstein M">M Feuerstein</name>
</author>
<author><name sortKey="Reichl, T" uniqKey="Reichl T">T Reichl</name>
</author>
<author><name sortKey="Vogel, J" uniqKey="Vogel J">J Vogel</name>
</author>
<author><name sortKey="Traub, J" uniqKey="Traub J">J Traub</name>
</author>
<author><name sortKey="Navab, N" uniqKey="Navab N">N Navab</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Franz, A" uniqKey="Franz A">A Franz</name>
</author>
<author><name sortKey="Haidegger, T" uniqKey="Haidegger T">T Haidegger</name>
</author>
<author><name sortKey="Birkfellner, W" uniqKey="Birkfellner W">W Birkfellner</name>
</author>
<author><name sortKey="Cleary, K" uniqKey="Cleary K">K Cleary</name>
</author>
<author><name sortKey="Peters, T" uniqKey="Peters T">T Peters</name>
</author>
<author><name sortKey="Maier Hein, L" uniqKey="Maier Hein L">L Maier-Hein</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Guerrero, J" uniqKey="Guerrero J">J Guerrero</name>
</author>
<author><name sortKey="Salcudean, S" uniqKey="Salcudean S">S Salcudean</name>
</author>
<author><name sortKey="Mcewen, J" uniqKey="Mcewen J">J McEwen</name>
</author>
<author><name sortKey="Masri, B" uniqKey="Masri B">B Masri</name>
</author>
<author><name sortKey="Nicolaou, S" uniqKey="Nicolaou S">S Nicolaou</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Kingham, Tp" uniqKey="Kingham T">TP Kingham</name>
</author>
<author><name sortKey="Jayaraman, S" uniqKey="Jayaraman S">S Jayaraman</name>
</author>
<author><name sortKey="Clements, Lw" uniqKey="Clements L">LW Clements</name>
</author>
<author><name sortKey="Scherer, Ma" uniqKey="Scherer M">MA Scherer</name>
</author>
<author><name sortKey="Stefansic, Jd" uniqKey="Stefansic J">JD Stefansic</name>
</author>
<author><name sortKey="Jarnagin, Wr" uniqKey="Jarnagin W">WR Jarnagin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lange, T" uniqKey="Lange T">T Lange</name>
</author>
<author><name sortKey="Eulenstein, S" uniqKey="Eulenstein S">S Eulenstein</name>
</author>
<author><name sortKey="Hunerbein, M" uniqKey="Hunerbein M">M Hünerbein</name>
</author>
<author><name sortKey="Schlag, Pm" uniqKey="Schlag P">PM Schlag</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lange, T" uniqKey="Lange T">T Lange</name>
</author>
<author><name sortKey="Papenberg, N" uniqKey="Papenberg N">N Papenberg</name>
</author>
<author><name sortKey="Heldmann, S" uniqKey="Heldmann S">S Heldmann</name>
</author>
<author><name sortKey="Modersitzki, J" uniqKey="Modersitzki J">J Modersitzki</name>
</author>
<author><name sortKey="Fischer, B" uniqKey="Fischer B">B Fischer</name>
</author>
<author><name sortKey="Lamecker, H" uniqKey="Lamecker H">H Lamecker</name>
</author>
<author><name sortKey="Schlag, Pm" uniqKey="Schlag P">PM Schlag</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mercier, L" uniqKey="Mercier L">L Mercier</name>
</author>
<author><name sortKey="Lang, T" uniqKey="Lang T">T Lang</name>
</author>
<author><name sortKey="Lindseth, F" uniqKey="Lindseth F">F Lindseth</name>
</author>
<author><name sortKey="Collins, Ld" uniqKey="Collins L">LD Collins</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Nicolau, S" uniqKey="Nicolau S">S Nicolau</name>
</author>
<author><name sortKey="Soler, L" uniqKey="Soler L">L Soler</name>
</author>
<author><name sortKey="Mutter, D" uniqKey="Mutter D">D Mutter</name>
</author>
<author><name sortKey="Marescaux, J" uniqKey="Marescaux J">J Marescaux</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Noble, J" uniqKey="Noble J">J Noble</name>
</author>
<author><name sortKey="Boukerroui, D" uniqKey="Boukerroui D">D Boukerroui</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Penney, Gp" uniqKey="Penney G">GP Penney</name>
</author>
<author><name sortKey="Blackall, Jm" uniqKey="Blackall J">JM Blackall</name>
</author>
<author><name sortKey="Hamady, M" uniqKey="Hamady M">M Hamady</name>
</author>
<author><name sortKey="Sabharwal, T" uniqKey="Sabharwal T">T Sabharwal</name>
</author>
<author><name sortKey="Adam, A" uniqKey="Adam A">A Adam</name>
</author>
<author><name sortKey="Hawkes, Dj" uniqKey="Hawkes D">DJ Hawkes</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Suwelack, S" uniqKey="Suwelack S">S Suwelack</name>
</author>
<author><name sortKey="Rhl, S" uniqKey="Rhl S">S Rhl</name>
</author>
<author><name sortKey="Bodenstedt, S" uniqKey="Bodenstedt S">S Bodenstedt</name>
</author>
<author><name sortKey="Reichard, D" uniqKey="Reichard D">D Reichard</name>
</author>
<author><name sortKey="Dillmann, R" uniqKey="Dillmann R">R Dillmann</name>
</author>
<author><name sortKey="Dos Santos, T" uniqKey="Dos Santos T">T dos Santos</name>
</author>
<author><name sortKey="Maier Hein, L" uniqKey="Maier Hein L">L Maier-Hein</name>
</author>
<author><name sortKey="Wagner, M" uniqKey="Wagner M">M Wagner</name>
</author>
<author><name sortKey="Wnscher, J" uniqKey="Wnscher J">J Wnscher</name>
</author>
<author><name sortKey="Kenngott, H" uniqKey="Kenngott H">H Kenngott</name>
</author>
<author><name sortKey="Mller, Bp" uniqKey="Mller B">BP Mller</name>
</author>
<author><name sortKey="Speidel, S" uniqKey="Speidel S">S Speidel</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></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">Int J Comput Assist Radiol Surg</journal-id>
<journal-id journal-id-type="iso-abbrev">Int J Comput Assist Radiol Surg</journal-id>
<journal-title-group><journal-title>International Journal of Computer Assisted Radiology and Surgery</journal-title>
</journal-title-group>
<issn pub-type="ppub">1861-6410</issn>
<issn pub-type="epub">1861-6429</issn>
<publisher><publisher-name>Springer Berlin Heidelberg</publisher-name>
<publisher-loc>Berlin/Heidelberg</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">26092658</article-id>
<article-id pub-id-type="pmc">4642598</article-id>
<article-id pub-id-type="publisher-id">1236</article-id>
<article-id pub-id-type="doi">10.1007/s11548-015-1236-8</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Original Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Locally rigid, vessel-based registration for laparoscopic liver surgery</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" corresp="yes"><name><surname>Song</surname>
<given-names>Yi</given-names>
</name>
<address><email>yi.song@ucl.ac.uk</email>
</address>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Totz</surname>
<given-names>Johannes</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Thompson</surname>
<given-names>Steve</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Johnsen</surname>
<given-names>Stian</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Barratt</surname>
<given-names>Dean</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Schneider</surname>
<given-names>Crispin</given-names>
</name>
<xref ref-type="aff" rid="Aff2"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Gurusamy</surname>
<given-names>Kurinchi</given-names>
</name>
<xref ref-type="aff" rid="Aff2"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Davidson</surname>
<given-names>Brian</given-names>
</name>
<xref ref-type="aff" rid="Aff2"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Ourselin</surname>
<given-names>Sébastien</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Hawkes</surname>
<given-names>David</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Clarkson</surname>
<given-names>Matthew J.</given-names>
</name>
<address><email>m.clarkson@ucl.ac.uk</email>
</address>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<aff id="Aff1"><label></label>
Centre For Medical Image Computing, Engineering Front Building, University College London, Malet Place, London, UK</aff>
<aff id="Aff2"><label></label>
Royal Free Campus, 9th Floor, Royal Free Hospital, UCL Medical School, Rowland Hill Street, London, UK</aff>
</contrib-group>
<pub-date pub-type="epub"><day>20</day>
<month>6</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>20</day>
<month>6</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="ppub"><year>2015</year>
</pub-date>
<volume>10</volume>
<issue>12</issue>
<fpage>1951</fpage>
<lpage>1961</lpage>
<history><date date-type="received"><day>24</day>
<month>3</month>
<year>2015</year>
</date>
<date date-type="accepted"><day>30</day>
<month>5</month>
<year>2015</year>
</date>
</history>
<permissions><copyright-statement>© The Author(s) 2015</copyright-statement>
<license license-type="OpenAccess"><license-p><bold>Open Access</bold>
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.</license-p>
</license>
</permissions>
<abstract id="Abs1"><sec><title>Purpose</title>
<p>Laparoscopic liver resection has significant advantages over open surgery due to less patient trauma and faster recovery times, yet is difficult for most lesions due to the restricted field of view and lack of haptic feedback. Image guidance provides a potential solution but is challenging in a soft deforming organ such as the liver. In this paper, we therefore propose a laparoscopic ultrasound (LUS) image guidance system and study the feasibility of a locally rigid registration for laparoscopic liver surgery.</p>
</sec>
<sec><title>Methods</title>
<p>We developed a real-time segmentation method to extract vessel centre points from calibrated, freehand, electromagnetically tracked, 2D LUS images. Using landmark-based initial registration and an optional iterative closest point (ICP) point-to-line registration, a vessel centre-line model extracted from preoperative computed tomography (CT) is registered to the ultrasound data during surgery.</p>
</sec>
<sec><title>Results</title>
<p>Using the locally rigid ICP method, the RMS residual error when registering to a phantom was 0.7 mm, and the mean target registration error (TRE) for two in vivo porcine studies was 3.58 and 2.99 mm, respectively. Using the locally rigid landmark-based registration method gave a mean TRE of 4.23 mm using vessel centre lines derived from CT scans taken with pneumoperitoneum and 6.57 mm without pneumoperitoneum.</p>
</sec>
<sec><title>Conclusion</title>
<p>In this paper we propose a practical image-guided surgery system based on locally rigid registration of a CT-derived model to vascular structures located with LUS. In a physical phantom and during porcine laparoscopic liver resection, we demonstrate accuracy of target location commensurate with surgical requirements. We conclude that locally rigid registration could be sufficient for practically useful image guidance in the near future.</p>
</sec>
</abstract>
<kwd-group xml:lang="en"><title>Keywords</title>
<kwd>Registration</kwd>
<kwd>Ultrasound</kwd>
<kwd>Laparoscopic</kwd>
<kwd>Liver surgery</kwd>
</kwd-group>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© CARS 2015</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1"><title>Introduction</title>
<p>In the UK, approximately 1800 liver resections are performed annually for primary or metastatic cancer. Liver cancer is a major global health problem, and 150,000 patients per year could benefit from liver resection. Currently, approximately 10 % of patients are considered suitable for laparoscopic liver resection, mainly those with small cancers on the periphery of the liver. Potentially, laparoscopic resection has significant benefits in reduced pain and cost savings due to shorter hospital stays [<xref ref-type="bibr" rid="CR7">7</xref>
]. Larger lesions and those close to major vascular/biliary structures are generally considered high risk for the laparoscopic approach mainly due to the restricted field of view and lack of haptic feedback.</p>
<p>We have developed a system that provides wider spatial context and potentially greater accuracy by aligning a preoperative plan derived from magnetic resonance (MR) or CT scans with the laparoscopic view. In this paper, we describe a freehand laparoscopic ultrasound (LUS)-based system that registers liver vessels in ultrasound (US) with MR/CT data. Specifically, we evaluate whether within a small region of interest a locally rigid registration is sufficiently accurate for surgical guidance.</p>
<sec id="Sec2"><title>Background</title>
<p>Previously reported commercial systems register using either surfaces of the liver reconstructed using a dragged pointer [<xref ref-type="bibr" rid="CR14">14</xref>
] or manual identification of four points (CAS-ONE<xref ref-type="fn" rid="Fn1">1</xref>
). The former will lead to errors due to direct contact with a soft tissue, while both are limited to a global rigid registration which is clearly unrealistic with the abdominal insufflation needed in laparoscopy. We have previously developed a system [<xref ref-type="bibr" rid="CR24">24</xref>
] for laparoscopic guidance based on dense stereo surface reconstruction [<xref ref-type="bibr" rid="CR25">25</xref>
] and an iterative closest point (ICP) [<xref ref-type="bibr" rid="CR5">5</xref>
]-based alignment to a surface derived from a preoperative CT model. However, the research literature suggests that deformable registration is a necessity for image guidance [<xref ref-type="bibr" rid="CR13">13</xref>
, <xref ref-type="bibr" rid="CR23">23</xref>
]. But, deformable models are difficult to validate [<xref ref-type="bibr" rid="CR19">19</xref>
] and may have multiple plausible solutions. It is also essential that a surgeon understands the registration accuracy while operating. We therefore propose a system based on locally rigid registration and test whether such a system is sufficiently accurate for surgical guidance. In this paper, a LUS probe is used to scan a local region of interest and update the global rigid registration based on the alignment of vessels in the preoperative CT and intraoperative ultrasound data within the region of interest.</p>
<p>In the literature, Aylward proposed rigid body registration of 3D B-mode ultrasound to preoperative CT for radio frequency ablation, based on a feature-to-image metric [<xref ref-type="bibr" rid="CR2">2</xref>
]. Lange, however, used a feature-to-feature method by extracting vessel centre lines from CT and 3D power Doppler ultrasound and then used ICP followed by multi-level B-splines for non-rigid alignment [<xref ref-type="bibr" rid="CR15">15</xref>
]. This was subsequently extended to incorporate vessel branch points as registration constraints [<xref ref-type="bibr" rid="CR16">16</xref>
]. The branch points were automatically identified in advance of surgery in the CT data, but selected manually in the ultrasound.</p>
<p>Accurate segmentation is a critical prerequisite for feature-based registration, and ultrasound image segmentation is itself a challenging problem, due to the low signal-to-noise ratio. Noble provides a thorough review [<xref ref-type="bibr" rid="CR20">20</xref>
]. Subsequently, Guerrero used an ellipse model to constrain an edge detection algorithm [<xref ref-type="bibr" rid="CR12">12</xref>
], thereby extracting vessels from ultrasound data for assessment of deep vein thrombosis. Later, Schneider used power Doppler ultrasound to initialize and guide vessel segmentation in B-mode images [<xref ref-type="bibr" rid="CR22">22</xref>
], replacing the previously required [<xref ref-type="bibr" rid="CR12">12</xref>
] manual initialisation of vessel centres. A scale-space blob detection approach has been used by Dagon et al. [<xref ref-type="bibr" rid="CR8">8</xref>
] and Anderegg et al. [<xref ref-type="bibr" rid="CR1">1</xref>
] to initialise vessel regions and approximate vessel walls using an ellipse model.</p>
<p>An alternative approach to feature-to-feature registration is image-to-image registration. Penney et al. [<xref ref-type="bibr" rid="CR21">21</xref>
] transformed a sparse set of freehand ultrasound slices to probability maps and registered with resampled and preprocessed CT data. Subsequently, Wein et al. [<xref ref-type="bibr" rid="CR26">26</xref>
] used a magnetic tracker to perform freehand 3D ultrasound registration of a sweep of data to preprocessed CT, using a <italic>semi-affine</italic>
(rotations, translations, 2 scaling, 1 skew) transformation. This work was extended to non-rigid deformation using B-splines and tested in a neurosurgical application [<xref ref-type="bibr" rid="CR27">27</xref>
].</p>
<p>Currently, there still exists challenges that are specific to the use of freehand LUS in surgical applications. The methods of Aylward et al. [<xref ref-type="bibr" rid="CR2">2</xref>
] and Lange et al. [<xref ref-type="bibr" rid="CR16">16</xref>
] are based on a 3D percutaneous probe. The probe is held stationary while a mechanical motor sweeps the ultrasound transducer in a predictable arc. Unfortunately, there are currently no commercially available laparoscopic 3D ultrasound probes. Wein’s work is based on a percutaneous probe, swept through a volume collecting a dense set of slices [<xref ref-type="bibr" rid="CR26">26</xref>
], and Penney’s work collects a sparse set of slices [<xref ref-type="bibr" rid="CR21">21</xref>
]. However, in a freehand laparoscopic setting, port positions and positioning of the LUS probe are often restrictive, and control of the motion during a sweep of data is often difficult, resulting in jerky motion. Moreover, the relatively small field of view makes the context difficult to interpret, and in our experience, it is often difficult to obtain elliptical vessel outlines.
</p>
</sec>
<sec id="Sec3"><title>Contribution of this paper</title>
<p>In this paper, we describe a registration system to align models derived from preoperative MR/CT data to intraoperative freehand ultrasound data taken using a 2D LUS probe. The method has similarities to the preceding literature in that we extract vessel centre lines as in [<xref ref-type="bibr" rid="CR1">1</xref>
, <xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR22">22</xref>
] and use an ICP registration as in [<xref ref-type="bibr" rid="CR15">15</xref>
, <xref ref-type="bibr" rid="CR22">22</xref>
]. In addition, to the best of our knowledge, while globally rigid [<xref ref-type="bibr" rid="CR22">22</xref>
] and additionally deformable [<xref ref-type="bibr" rid="CR1">1</xref>
, <xref ref-type="bibr" rid="CR8">8</xref>
] registration of vessel models from CT and US data have been proposed, an evaluation of registration accuracy in vivo in a clinically usable laparoscopic ultrasound system has not been reported. The hypothesis of this paper is that local rigid registration within a small region of interest is sufficient for image guidance without deformable modelling, and the specific contributions of this paper are the delivery of a system to achieve that, and a thorough evaluation of errors using a phantom and during porcine laparoscopic liver resection.</p>
</sec>
</sec>
<sec id="Sec4"><title>Methods</title>
<p>Figure <xref rid="Fig1" ref-type="fig">1</xref>
shows an overview of the registration process. Vessel centre points <italic>P</italic>
are detected in 2D ultrasound images and converted into 3D space via the ultrasound calibration and tracking transformations. The preoperative CT scan is preprocessed to extract a graph <italic>G</italic>
representing vessel centre lines. The ultrasound- derived data <italic>P</italic>
and CT-derived data <italic>G</italic>
are registered using manually picked landmarks and/or the ICP algorithm. The rigid body registration transformation <inline-formula id="IEq2"><alternatives><tex-math id="M1">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^GT_{P}$$\end{document}</tex-math>
<mml:math id="M2"><mml:mrow><mml:msup><mml:mrow></mml:mrow>
<mml:mi>G</mml:mi>
</mml:msup>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq2.gif"></inline-graphic>
</alternatives>
</inline-formula>
enables the preoperative data to be visualised relative to the live ultrasound imaging plane, as shown in Fig. <xref rid="Fig2" ref-type="fig">2</xref>
.<fig id="Fig1"><label>Fig. 1</label>
<caption><p>Overview of the registration process. Vessel centre points <italic>P</italic>
from ultrasound data are registered to a vessel centre-line graph <italic>G</italic>
giving rigid body transformation <inline-formula id="IEq1"><alternatives><tex-math id="M3">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^GT_{P}$$\end{document}</tex-math>
<mml:math id="M4"><mml:mrow><mml:msup><mml:mrow></mml:mrow>
<mml:mi>G</mml:mi>
</mml:msup>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq1.gif"></inline-graphic>
</alternatives>
</inline-formula>
</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig1_HTML" id="MO1"></graphic>
</fig>
<fig id="Fig2"><label>Fig. 2</label>
<caption><p>Applying the registration transformation to anatomical models derived from preoperative CT data enables live visualisation of CT data, within the context of live laparoscopic video and ultrasound data</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig2_HTML" id="MO2"></graphic>
</fig>
</p>
<sec id="Sec5"><title>Preprocessing preoperative data</title>
<p>A standard clinical tri-phase abdominal CT scan is obtained and segmented<xref ref-type="fn" rid="Fn2">2</xref>
to represent important structures such as the liver, tumours, arteries, hepatic vein, portal vein, gall bladder. Centre lines are extracted using the Vascular Modelling Tool Kit (VMTK).<xref ref-type="fn" rid="Fn3">3</xref>
This yields a vessel graph <italic>G</italic>
, which can be trivially processed to identify vessel bifurcation points.
</p>
</sec>
<sec id="Sec6"><title>Real-time Ultrasound Segmentation</title>
<p>Previous works on 2D ultrasound vessel segmentation use an ellipse model to constrain the edge detection process [<xref ref-type="bibr" rid="CR1">1</xref>
, <xref ref-type="bibr" rid="CR8">8</xref>
, <xref ref-type="bibr" rid="CR12">12</xref>
]. This approach assumes that vessels are imaged approximately perpendicular to the vessel centre line, which is not practical for laparoscopic use where movement is often restricted by the position of a trocar. Moreover, it is not clear how this approach handles topological changes of the external contours of vessels in the 2D US images. Therefore, we propose a flexible segmentation method that is not limited to cross-sectional scans and can cope with topology changes during the course of scanning. An example is shown in Fig. <xref rid="Fig3" ref-type="fig">3</xref>
.<fig id="Fig3"><label>Fig. 3</label>
<caption><p>Vessel segmentation: <bold>a</bold>
ultrasound B-mode image. <bold>b</bold>
Vessel-enhanced image. <bold>c</bold>
Thresholded vessel-enhanced image. <bold>d</bold>
<italic>Dip</italic>
image [<xref ref-type="bibr" rid="CR21">21</xref>
]. <bold>e</bold>
Thresholded <italic>Dip</italic>
image. <bold>f</bold>
The candidate seeds of vessels after thresholded vessel- enhance image is masked with the thresholded <italic>Dip</italic>
image. <bold>g</bold>
Vessel contours are depicted in <italic>red</italic>
, fitted ellipses, and centres are in <italic>green</italic>
</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig3_HTML" id="MO3"></graphic>
</fig>
</p>
<sec id="Sec7"><title>Vessel enhancement image</title>
<p>The standard B-mode ultrasound images have a low signal-to-noise ratio (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
a), so vessel structures are first enhanced for more reliable vessel segmentation. The multi-scale vessel enhancement filter [<xref ref-type="bibr" rid="CR10">10</xref>
] is used, which is based on eigenvalue analysis of the Hessian. The eigen values are ordered as <inline-formula id="IEq5"><alternatives><tex-math id="M5">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$| \lambda _1 | < | \lambda _2 |$$\end{document}</tex-math>
<mml:math id="M6"><mml:mrow><mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:msub><mml:mi mathvariant="italic">λ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
<mml:mo><</mml:mo>
<mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
<mml:msub><mml:mi mathvariant="italic">λ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow><mml:mo stretchy="false">|</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq5.gif"></inline-graphic>
</alternatives>
</inline-formula>
. The 2D “vesselness” of a pixel is measured by<disp-formula id="Equ1"><label>1</label>
<alternatives><tex-math id="M7">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} v_0(s) = {\left\{ \begin{array}{ll} 0 &{} \text {if } \lambda _2 < 0 \\ \hbox {e}^{\left( - \frac{{R_B}^2}{2\beta ^2} \right) } \left( 1 - \hbox {e}^{\left( - \frac{S^2}{2c^2}\right) } \right) &{} \end{array}\right. } \end{aligned}$$\end{document}</tex-math>
<mml:math id="M8" display="block"><mml:mrow><mml:mtable columnspacing="0.5ex"><mml:mtr><mml:mtd columnalign="right"><mml:mrow><mml:msub><mml:mi>v</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:mi>s</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mfenced close="" open="{" separators=""><mml:mrow><mml:mtable columnspacing="0.5ex"><mml:mtr><mml:mtd columnalign="left"><mml:mn>0</mml:mn>
</mml:mtd>
<mml:mtd columnalign="left"><mml:mrow><mml:mrow></mml:mrow>
<mml:mtext>if</mml:mtext>
<mml:mspace width="0.333333em"></mml:mspace>
<mml:msub><mml:mi mathvariant="italic">λ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo><</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd columnalign="left"><mml:mrow><mml:mrow></mml:mrow>
<mml:msup><mml:mtext>e</mml:mtext>
<mml:mfenced close=")" open="(" separators=""><mml:mo>-</mml:mo>
<mml:mfrac><mml:msup><mml:mrow><mml:msub><mml:mi>R</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow><mml:mn>2</mml:mn>
<mml:msup><mml:mi mathvariant="italic">β</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mfenced>
</mml:msup>
<mml:mfenced close=")" open="(" separators=""><mml:mn>1</mml:mn>
<mml:mo>-</mml:mo>
<mml:msup><mml:mtext>e</mml:mtext>
<mml:mfenced close=")" open="(" separators=""><mml:mo>-</mml:mo>
<mml:mfrac><mml:msup><mml:mi>S</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow><mml:mn>2</mml:mn>
<mml:msup><mml:mi>c</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mfenced>
</mml:msup>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
<mml:mtd columnalign="left"><mml:mrow></mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<graphic xlink:href="11548_2015_1236_Article_Equ1.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
where<disp-formula id="Equ2"><label>2</label>
<alternatives><tex-math id="M9">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&R_B = \lambda _1 / \lambda _2 \end{aligned}$$\end{document}</tex-math>
<mml:math id="M10" display="block"><mml:mrow><mml:mtable columnspacing="0.5ex"><mml:mtr><mml:mtd></mml:mtd>
<mml:mtd columnalign="left"><mml:mrow><mml:msub><mml:mi>R</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub><mml:mi mathvariant="italic">λ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo stretchy="false">/</mml:mo>
<mml:msub><mml:mi mathvariant="italic">λ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<graphic xlink:href="11548_2015_1236_Article_Equ2.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<disp-formula id="Equ3"><label>3</label>
<alternatives><tex-math id="M11">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned}&S = \sqrt{{\lambda _1}^2 + {\lambda _2}^2} \end{aligned}$$\end{document}</tex-math>
<mml:math id="M12" display="block"><mml:mrow><mml:mtable columnspacing="0.5ex"><mml:mtr><mml:mtd></mml:mtd>
<mml:mtd columnalign="left"><mml:mrow><mml:mi>S</mml:mi>
<mml:mo>=</mml:mo>
<mml:msqrt><mml:mrow><mml:msup><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msup><mml:mrow><mml:msub><mml:mi mathvariant="italic">λ</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:math>
<graphic xlink:href="11548_2015_1236_Article_Equ3.gif" position="anchor"></graphic>
</alternatives>
</disp-formula>
<inline-formula id="IEq6"><alternatives><tex-math id="M13">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\beta = 1$$\end{document}</tex-math>
<mml:math id="M14"><mml:mrow><mml:mi mathvariant="italic">β</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq6.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq7"><alternatives><tex-math id="M15">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$c = 10$$\end{document}</tex-math>
<mml:math id="M16"><mml:mrow><mml:mi>c</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>10</mml:mn>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq7.gif"></inline-graphic>
</alternatives>
</inline-formula>
are thresholds which control the sensitivity of the line filter to the measures <inline-formula id="IEq8"><alternatives><tex-math id="M17">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$${R_B}$$\end{document}</tex-math>
<mml:math id="M18"><mml:msub><mml:mi>R</mml:mi>
<mml:mi>B</mml:mi>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq8.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <italic>S</italic>
. In Fig. <xref rid="Fig3" ref-type="fig">3</xref>
b, it can be seen that some common artefacts on the ultrasound images, e.g. shadows, are wrongly picked up by the enhancement filter. For many cases, using only the prior knowledge of the vessel intensity distributions is not sufficient to exclude those non-vessel regions. To improve robustness, we adopt the <italic>Dip</italic>
image as proposed by Penny et al. [<xref ref-type="bibr" rid="CR21">21</xref>
].</p>
</sec>
<sec id="Sec8"><title>Creation of the Dip image</title>
<p>The <italic>Dip</italic>
image (<inline-formula id="IEq9"><alternatives><tex-math id="M19">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$I_{dip})$$\end{document}</tex-math>
<mml:math id="M20"><mml:mrow><mml:msub><mml:mi>I</mml:mi>
<mml:mrow><mml:mi>d</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow><mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq9.gif"></inline-graphic>
</alternatives>
</inline-formula>
was originally designed to produce vessel probability maps via a training data set. In this paper, we only use the intensity differences (i.e. intensity dips) between regions of interest. The size of a region is determined by the diameter of vessels. No additional artefact removal step is required, except for a Gaussian filter over the US image. Since we currently target the left liver lobe for surgical guidance, we set the search range of vessel diameters from 9 to 3 mm (roughly equal to 100–40 pixels on the LUS image) as a porcine left lobe features relatively large vessels.</p>
<p>The <italic>Dip</italic>
image is computed along the beam direction. As we use a linear LUS probe, the beam directions can be modelled as image columns. Figure <xref rid="Fig4" ref-type="fig">4</xref>
illustrates the calculation of three mean intensity values <italic>a</italic>
, <italic>b</italic>
and <italic>c</italic>
, within regions <inline-formula id="IEq10"><alternatives><tex-math id="M21">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$[x+v/2 , x+v], [x-v , x-v/2], [x-v/2 , x+v/2]$$\end{document}</tex-math>
<mml:math id="M22"><mml:mrow><mml:mo stretchy="false">[</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">]</mml:mo>
<mml:mo>,</mml:mo>
<mml:mo stretchy="false">[</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo stretchy="false">]</mml:mo>
<mml:mo>,</mml:mo>
<mml:mo stretchy="false">[</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq10.gif"></inline-graphic>
</alternatives>
</inline-formula>
, respectively, with <italic>x</italic>
being a pixel at the <italic>i</italic>
th column and <italic>v</italic>
the vessel width. If <inline-formula id="IEq11"><alternatives><tex-math id="M23">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$c < b$$\end{document}</tex-math>
<mml:math id="M24"><mml:mrow><mml:mi>c</mml:mi>
<mml:mo><</mml:mo>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq11.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq12"><alternatives><tex-math id="M25">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$c < a$$\end{document}</tex-math>
<mml:math id="M26"><mml:mrow><mml:mi>c</mml:mi>
<mml:mo><</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq12.gif"></inline-graphic>
</alternatives>
</inline-formula>
, every pixel in <inline-formula id="IEq13"><alternatives><tex-math id="M27">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$[x-v/2, x+v/2]$$\end{document}</tex-math>
<mml:math id="M28"><mml:mrow><mml:mo stretchy="false">[</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq13.gif"></inline-graphic>
</alternatives>
</inline-formula>
on the <italic>Dip</italic>
image will have the value <inline-formula id="IEq14"><alternatives><tex-math id="M29">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$b_v = \hbox {min} (a-c, b-c)$$\end{document}</tex-math>
<mml:math id="M30"><mml:mrow><mml:msub><mml:mi>b</mml:mi>
<mml:mi>v</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:mi>a</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>c</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>b</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>c</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq14.gif"></inline-graphic>
</alternatives>
</inline-formula>
. This process is repeated for each v in <inline-formula id="IEq15"><alternatives><tex-math id="M31">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$[v_\mathrm{min}, v_\mathrm{max}]$$\end{document}</tex-math>
<mml:math id="M32"><mml:mrow><mml:mo stretchy="false">[</mml:mo>
<mml:msub><mml:mi>v</mml:mi>
<mml:mi mathvariant="normal">min</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>v</mml:mi>
<mml:mi mathvariant="normal">max</mml:mi>
</mml:msub>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq15.gif"></inline-graphic>
</alternatives>
</inline-formula>
. The final pixel values at position <inline-formula id="IEq16"><alternatives><tex-math id="M33">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$[x-v/2, x+v/2]$$\end{document}</tex-math>
<mml:math id="M34"><mml:mrow><mml:mo stretchy="false">[</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>-</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>x</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq16.gif"></inline-graphic>
</alternatives>
</inline-formula>
will be <inline-formula id="IEq17"><alternatives><tex-math id="M35">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\hbox {max}(b_v)$$\end{document}</tex-math>
<mml:math id="M36"><mml:mrow><mml:mtext>max</mml:mtext>
<mml:mo stretchy="false">(</mml:mo>
<mml:msub><mml:mi>b</mml:mi>
<mml:mi>v</mml:mi>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq17.gif"></inline-graphic>
</alternatives>
</inline-formula>
. The steps above are repeated for every column of the US image and all pixels along that column. This can be parallelised easily as each column is processed independently of others. To reduce the search range of vessel diameters, a coarse-to-fine pyramidal approach is proposed to speed up the process further.<fig id="Fig4"><label>Fig. 4</label>
<caption><p>Creation of the Dip image. <italic>Left</italic>
Gaussian blurred ultrasound image <inline-formula id="IEq3"><alternatives><tex-math id="M37">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(I_{US})$$\end{document}</tex-math>
<mml:math id="M38"><mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:msub><mml:mi>I</mml:mi>
<mml:mrow><mml:mi>U</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq3.gif"></inline-graphic>
</alternatives>
</inline-formula>
. <italic>Centre</italic>
The intensity profile along line <inline-formula id="IEq4"><alternatives><tex-math id="M39">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(x_0, x_n)$$\end{document}</tex-math>
<mml:math id="M40"><mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:msub><mml:mi>x</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>x</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq4.gif"></inline-graphic>
</alternatives>
</inline-formula>
. The location and size of image regions which gives the values <italic>a</italic>
, <italic>b</italic>
and <italic>c</italic>
. <italic>Right</italic>
The resulting <italic>Dip</italic>
image</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig4_HTML" id="MO4"></graphic>
</fig>
</p>
</sec>
<sec id="Sec9"><title>Segmentation and reconstruction</title>
<p>The vessel-enhanced image is thresholded at <inline-formula id="IEq18"><alternatives><tex-math id="M41">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$T_\mathrm{e}$$\end{document}</tex-math>
<mml:math id="M42"><mml:msub><mml:mi>T</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq18.gif"></inline-graphic>
</alternatives>
</inline-formula>
to eliminate background noise; Fig. <xref rid="Fig3" ref-type="fig">3</xref>
c. We create a mask image <inline-formula id="IEq19"><alternatives><tex-math id="M43">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(I_\mathrm{mask})$$\end{document}</tex-math>
<mml:math id="M44"><mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:msub><mml:mi>I</mml:mi>
<mml:mi mathvariant="normal">mask</mml:mi>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq19.gif"></inline-graphic>
</alternatives>
</inline-formula>
by applying a threshold <inline-formula id="IEq20"><alternatives><tex-math id="M45">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(T_\mathrm{d})$$\end{document}</tex-math>
<mml:math id="M46"><mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq20.gif"></inline-graphic>
</alternatives>
</inline-formula>
to the <italic>Dip</italic>
image which is set as half of maximum value of the <italic>Dip</italic>
image; Fig. <xref rid="Fig3" ref-type="fig">3</xref>
e. These two thresholds are set for the given B-mode imaging parameters, e.g. gain, power, map.</p>
<p>The de-noised vessel-enhanced image is masked with <inline-formula id="IEq21"><alternatives><tex-math id="M47">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$I_\mathrm{mask}$$\end{document}</tex-math>
<mml:math id="M48"><mml:msub><mml:mi>I</mml:mi>
<mml:mi mathvariant="normal">mask</mml:mi>
</mml:msub>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq21.gif"></inline-graphic>
</alternatives>
</inline-formula>
. Regions appearing on both images are kept (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
f). The intensity distribution of those regions will be further compared against the prior knowledge of vessel intensity and will be removed if they are not matching, i.e. falling out of the vessel intensity range. The remaining pixels are candidate vessel seeds. The regions in the de-noised vessel enhancement image containing such candidate seeds are identified as vessels and their contours are detected.</p>
<p>Since vessel centre points are employed for registration in this paper, ellipses are fitted to those contours to derive centre points in each ultrasound image (Fig. <xref rid="Fig3" ref-type="fig">3</xref>
g). Outliers can be excluded by defining minimal and maximal ellipse axes ratio and length, as demonstrated in Fig. <xref rid="Fig5" ref-type="fig">5</xref>
. For example, when an image is scanned near parallel to a vessel centre-line direction, it results in large ellipse axes. This can be removed by constraining the short axis length to the pre-defined vessel diameter range <inline-formula id="IEq22"><alternatives><tex-math id="M49">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$[v_\mathrm{min}, v_\mathrm{max}]$$\end{document}</tex-math>
<mml:math id="M50"><mml:mrow><mml:mo stretchy="false">[</mml:mo>
<mml:msub><mml:mi>v</mml:mi>
<mml:mi mathvariant="normal">min</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub><mml:mi>v</mml:mi>
<mml:mi mathvariant="normal">max</mml:mi>
</mml:msub>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq22.gif"></inline-graphic>
</alternatives>
</inline-formula>
, as described in section “Creation of the Dip image”. An additional criterion is that the axes ratio should be larger than 0.5. Otherwise, the vessel could be scanned less than 30<inline-formula id="IEq23"><alternatives><tex-math id="M51">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^\circ $$\end{document}</tex-math>
<mml:math id="M52"><mml:msup><mml:mrow></mml:mrow>
<mml:mo>∘</mml:mo>
</mml:msup>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq23.gif"></inline-graphic>
</alternatives>
</inline-formula>
to its centre-line direction, which often does not produce reliable ellipse centres.<fig id="Fig5"><label>Fig. 5</label>
<caption><p>Example of outlier rejection. The ellipse is fitted to the vessel outline but the detected centre is rejected by the ellipse axes criteria</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig5_HTML" id="MO8"></graphic>
</fig>
</p>
<p>The vessel centres, in 2D pixel coordinates are multiplied by the ultrasound calibration, and the probe tracking transformation and hence converted into 3D data points (<italic>P</italic>
), which are used to register the preoperative CT data to the patient in the operation room.</p>
</sec>
</sec>
<sec id="Sec10"><title>Registration</title>
<p>Figure <xref rid="Fig6" ref-type="fig">6</xref>
a, b illustrates the landmarks and vectors used for initial alignment. A landmark <italic>L</italic>
and two vectors, <inline-formula id="IEq24"><alternatives><tex-math id="M53">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{u}$$\end{document}</tex-math>
<mml:math id="M54"><mml:mi mathvariant="bold">u</mml:mi>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq24.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq25"><alternatives><tex-math id="M55">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbf{v}$$\end{document}</tex-math>
<mml:math id="M56"><mml:mi mathvariant="bold">v</mml:mi>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq25.gif"></inline-graphic>
</alternatives>
</inline-formula>
, are defined on the preoperative centre-line model <italic>G</italic>
, along with their correspondences <inline-formula id="IEq26"><alternatives><tex-math id="M57">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$L', \mathbf{u', v'}$$\end{document}</tex-math>
<mml:math id="M58"><mml:mrow><mml:msup><mml:mi>L</mml:mi>
<mml:mo>′</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mrow><mml:msup><mml:mi mathvariant="bold">u</mml:mi>
<mml:mo>′</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:msup><mml:mi mathvariant="bold">v</mml:mi>
<mml:mo>′</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq26.gif"></inline-graphic>
</alternatives>
</inline-formula>
in the derived centre points <italic>P</italic>
. Currently, this was done manually. An initial rigid registration is obtained by the alignment of landmarks <inline-formula id="IEq27"><alternatives><tex-math id="M59">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\{L, L'\}$$\end{document}</tex-math>
<mml:math id="M60"><mml:mrow><mml:mo stretchy="false">{</mml:mo>
<mml:mi>L</mml:mi>
<mml:mo>,</mml:mo>
<mml:msup><mml:mi>L</mml:mi>
<mml:mo>′</mml:mo>
</mml:msup>
<mml:mo stretchy="false">}</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq27.gif"></inline-graphic>
</alternatives>
</inline-formula>
which gives the translation, and vectors <inline-formula id="IEq28"><alternatives><tex-math id="M61">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\{\mathbf{u, u'}\}$$\end{document}</tex-math>
<mml:math id="M62"><mml:mrow><mml:mo stretchy="false">{</mml:mo>
<mml:mrow><mml:mi mathvariant="bold">u</mml:mi>
<mml:mo>,</mml:mo>
<mml:msup><mml:mi mathvariant="bold">u</mml:mi>
<mml:mo>′</mml:mo>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">}</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq28.gif"></inline-graphic>
</alternatives>
</inline-formula>
and <inline-formula id="IEq29"><alternatives><tex-math id="M63">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\{\mathbf{v, v'}\}$$\end{document}</tex-math>
<mml:math id="M64"><mml:mrow><mml:mo stretchy="false">{</mml:mo>
<mml:mrow><mml:mi mathvariant="bold">v</mml:mi>
<mml:mo>,</mml:mo>
<mml:msup><mml:mi mathvariant="bold">v</mml:mi>
<mml:mo>′</mml:mo>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">}</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq29.gif"></inline-graphic>
</alternatives>
</inline-formula>
which computes the rotation. After the initial alignment, the ICP algorithm [<xref ref-type="bibr" rid="CR5">5</xref>
] is applied to further refine the registration of preoperative data <italic>G</italic>
to the intraoperative data <italic>P</italic>
(Fig. <xref rid="Fig6" ref-type="fig">6</xref>
c).<fig id="Fig6"><label>Fig. 6</label>
<caption><p>Example of corresponding landmarks and vectors in the hepatic vein, as used for initial alignment. <bold>a</bold>
Intraoperative centre points <italic>P</italic>
. <bold>b</bold>
Preoperative centre-line model <italic>G</italic>
. <bold>c</bold>
Preoperative centre-line model <italic>G</italic>
is aligned to intraoperative centre points <italic>P</italic>
using ICP</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig6_HTML" id="MO9"></graphic>
</fig>
</p>
</sec>
</sec>
<sec id="Sec11"><title>Experiments and results</title>
<p>Experiments were performed to determine the overall registration accuracy of the system, and to identify sources of error from various component parts (sections “Ultrasound calibration error” and “Vessel segmentation error”). Our system uses an electromagnetic (EM) tracker, which is known to display tracking inaccuracies due to magnetic field inhomogeneities [<xref ref-type="bibr" rid="CR11">11</xref>
]. Various works have tried to mitigate against EM tracking inaccuracies by calibration [<xref ref-type="bibr" rid="CR18">18</xref>
] and combination with optical trackers [<xref ref-type="bibr" rid="CR9">9</xref>
]. The focus of this paper is the practicalities of intraoperative registration, so we refer to the manufacturer-claimed position accuracy of 1.4 mm RMS and orientation accuracy of 0.5<inline-formula id="IEq30"><alternatives><tex-math id="M65">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^{\circ }$$\end{document}</tex-math>
<mml:math id="M66"><mml:msup><mml:mrow></mml:mrow>
<mml:mo>∘</mml:mo>
</mml:msup>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq30.gif"></inline-graphic>
</alternatives>
</inline-formula>
RMS. A fundamental point for surgical navigation is that while the presented algorithm determines the registration transformation <inline-formula id="IEq31"><alternatives><tex-math id="M67">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$^{P}T_{G}$$\end{document}</tex-math>
<mml:math id="M68"><mml:mrow><mml:msup><mml:mrow></mml:mrow>
<mml:mi>P</mml:mi>
</mml:msup>
<mml:msub><mml:mi>T</mml:mi>
<mml:mi>G</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq31.gif"></inline-graphic>
</alternatives>
</inline-formula>
from preoperative data <italic>G</italic>
to intraoperative data <italic>P</italic>
, the actual navigation accuracy will be the composition of the registration accuracy, the EM tracking accuracy as the probe moves, the US calibration accuracy and the deformation of the liver due to the US probe itself. For this reason, separate data are used to assess registration accuracy (section “Registration accuracy: in vivo”), and navigation accuracy (section “Navigation accuracy: in vivo”). In experiments “Registration accuracy: in vivo” and “Navigation accuracy: in vivo”, we use vessel models derived from CT scans taken using pneumoperitoneum (insufflated), which are not available clinically. So in section “Comparison of insufflated versus non-insufflated models”, we specifically compare registration and navigation accuracy when registering to CT-derived vessel models using pneumoperitoneum (insufflated) and without pneumoperitoneum (non-insufflated). US images were collected under controlled breathing (Boyles apparatus), which is discussed later.</p>
<sec id="Sec12"><title>Experimental set-up</title>
<p>Our data acquisition system is built upon the NifTK platform [<xref ref-type="bibr" rid="CR6">6</xref>
]. Live LUS images are acquired at 25 frames per second (fps). We used an Analogic<xref ref-type="fn" rid="Fn4">4</xref>
SonixMDP, and a Vermon<xref ref-type="fn" rid="Fn5">5</xref>
LP7 linear probe. An Ascension<xref ref-type="fn" rid="Fn6">6</xref>
3D Guidance medSafe mid-range electromagnetic (EM) tracker was used to track the LUS probe at 60 fps via a six-degrees-of-freedom (6-DOF) sensor (Model 180) attached to the articulated tip.</p>
</sec>
<sec id="Sec13"><title>Ultrasound calibration error</title>
<p>The LUS probe was calibrated at a scanning depth of 45 mm before surgery using an invariant point method [<xref ref-type="bibr" rid="CR17">17</xref>
]. The scanning depth of the LUS probe was not changed throughout our experiments. The validation phantom is shown in Fig. <xref rid="Fig7" ref-type="fig">7</xref>
a, and described further in [<xref ref-type="bibr" rid="CR4">4</xref>
]. Eight pins on the phantom were scanned in turn using the LUS probe. The pin heads were manually segmented from the US images, Fig. <xref rid="Fig7" ref-type="fig">7</xref>
b. 100 frames were collected at each pin to minimise the impact of manual segmentation error. Their 3D positions in the EM coordinate system were computed by multiplying the 2D pixel location by the calibration transformation and then the EM tracking transformation, Fig. <xref rid="Fig7" ref-type="fig">7</xref>
c. The accuracy of these computed 3D positions were validated based on two ground truths. The first ground truth is the known geometry of the 8-pin phantom, where the pins are arranged on a <inline-formula id="IEq32"><alternatives><tex-math id="M69">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$4 \times 2$$\end{document}</tex-math>
<mml:math id="M70"><mml:mrow><mml:mn>4</mml:mn>
<mml:mo>×</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq32.gif"></inline-graphic>
</alternatives>
</inline-formula>
grid, with each side being 25 mm in length. The resulting mean edge length was 24.62 mm. The second ground truth is the physical positions of the eight phantom pins in the EM coordinate system, which are measured by using another EM sensor tracked by the same EM transmitter, Fig. <xref rid="Fig7" ref-type="fig">7</xref>
c. The distance between each reconstructed pin and its ground truth position is listed in Table <xref rid="Tab1" ref-type="table">1</xref>
.<fig id="Fig7"><label>Fig. 7</label>
<caption><p>Evaluation of ultrasound calibration using an eight- point phantom. <bold>a</bold>
Eight-point phantom. <bold>b</bold>
A LUS B-mode scan of pins on phantom. <bold>c</bold>
3D positions of eight pins obtained from tracked LUS scans are depicted in <italic>yellow</italic>
. The ground truth position of eight pins is depicted in <italic>green</italic>
</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig7_HTML" id="MO10"></graphic>
</fig>
<table-wrap id="Tab1"><label>Table 1</label>
<caption><p>Error measures for each reconstructed pin position</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th align="left" rowspan="2"></th>
<th align="left" colspan="8">Pin number</th>
</tr>
<tr><th align="left">1</th>
<th align="left">2</th>
<th align="left">3</th>
<th align="left">4</th>
<th align="left">5</th>
<th align="left">6</th>
<th align="left">7</th>
<th align="left">8</th>
</tr>
</thead>
<tbody><tr><td align="left">RMS error (mm)</td>
<td align="left">2.89</td>
<td align="left">3.40</td>
<td align="left">1.28</td>
<td align="left">0.81</td>
<td align="left">2.35</td>
<td align="left">1.59</td>
<td align="left">2.20</td>
<td align="left">2.82</td>
</tr>
</tbody>
</table>
<table-wrap-foot><p>See Fig. <xref rid="Fig7" ref-type="fig">7</xref>
c</p>
</table-wrap-foot>
</table-wrap>
</p>
</sec>
<sec id="Sec14"><title>Vessel segmentation error</title>
<p>LUS images were acquired from a phantom made from Agar. The phantom contained tubular structures filled with water, as shown in Fig. <xref rid="Fig8" ref-type="fig">8</xref>
a, b. The ground truth is the diameter of the tubular structures manufactured as 6.5 mm. One hundred and sixty images (<inline-formula id="IEq33"><alternatives><tex-math id="M71">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$640\times 480$$\end{document}</tex-math>
<mml:math id="M72"><mml:mrow><mml:mn>640</mml:mn>
<mml:mo>×</mml:mo>
<mml:mn>480</mml:mn>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq33.gif"></inline-graphic>
</alternatives>
</inline-formula>
pixels) were collected. The contours of the tubular structures were automatically segmented from the US images and fitted with ellipses, so that the short ellipse axis approximated the diameter of the tubular structures, Fig. <xref rid="Fig8" ref-type="fig">8</xref>
c–e. The resulting mean (standard deviation) diameter of the segmented contours was <inline-formula id="IEq34"><alternatives><tex-math id="M73">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$6.4\;(0.17)$$\end{document}</tex-math>
<mml:math id="M74"><mml:mrow><mml:mn>6.4</mml:mn>
<mml:mspace width="0.277778em"></mml:mspace>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>0.17</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq34.gif"></inline-graphic>
</alternatives>
</inline-formula>
mm. The average time of the image processing on one US image was 100 ms.<fig id="Fig8"><label>Fig. 8</label>
<caption><p>Validation of vessel segmentation on a phantom. <bold>a</bold>
Phantom design. The rods will be removed after filling the <italic>box</italic>
with Agar. <bold>b</bold>
LUS probe sweeping across the surface of the phantom which is made from Agar. An EM sensor is attached to the probe and tracked. <bold>c</bold>
–<bold>e</bold>
LUS images of the tubular structures. The outlines are depicted in <italic>red</italic>
. The ellipses fitted to the outlines are depicted in <italic>green</italic>
. The extracted ellipse centres are depicted in <italic>green</italic>
</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig8_HTML" id="MO11"></graphic>
</fig>
</p>
</sec>
<sec id="Sec15"><title>Registration accuracy: phantom</title>
<p>The registration accuracy was assessed on the same phantom as section “Vessel segmentation error”, Fig. <xref rid="Fig8" ref-type="fig">8</xref>
. Using the presented algorithm, the tubular structures were automatically segmented, the centre points extracted, and converted to EM coordinates by multiplication with the US calibration matrix and EM tracker matrix. These reconstructed points were rigidly registered to the centre lines of the phantom tubular structures using the ICP method. Figure <xref rid="Fig9" ref-type="fig">9</xref>
illustrates the registration of reconstructed points to the phantom model. The RMS residual error given by the ICP method was 0.7 mm.<fig id="Fig9"><label>Fig. 9</label>
<caption><p>Validation of vessel registration on a phantom: the reconstructed contours from the ultrasound data (<italic>yellow rings</italic>
) were rigidly registered to the phantom using ICP</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig9_HTML" id="MO12"></graphic>
</fig>
</p>
</sec>
<sec id="Sec16"><title>Registration accuracy: in vivo</title>
<p>The overall registration accuracy was evaluated during porcine laparoscopic liver resection using two studies of the same subject. The LUS images were acquired from the left lobe of the liver, before and after a significant repositioning of the lobe. The surgeon swept the liver surface in a steady way to make sure vessel centre points were densely sampled and gently so as not to cause significant deformation of the liver surface. The US imaging parameters for brightness, contrast and gain control were preset values and not changed during scanning. About 10 images per second were segmented. In the first study, in total 370 images (640 <inline-formula id="IEq35"><alternatives><tex-math id="M75">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\times $$\end{document}</tex-math>
<mml:math id="M76"><mml:mo>×</mml:mo>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq35.gif"></inline-graphic>
</alternatives>
</inline-formula>
480 pixels) were processed. In the second study, 340 images were processed. The detected vessel centres were converted into 3D data points <italic>P</italic>
. Two tri-phase clinical CT scans had been obtained a week earlier, one with insufflation (12 mm Hg) and one without. Vessel centre lines were extracted using the model derived from the insufflated CT scan. The registration method of section “Registration” was applied, registering the preoperative centre-line model <italic>G</italic>
to the intraoperative data set <italic>P</italic>
. For the first study, eight bifurcations were manually identified and labelled in both the US images and the CT data, to be used for anatomical targets, as shown in Fig. <xref rid="Fig10" ref-type="fig">10</xref>
b. The mean TRE was 3.58 mm, and the maximum TRE was 5.76 mm. For the second study, three bifurcations (i.e. number 1, 2, 4 in Fig. <xref rid="Fig10" ref-type="fig">10</xref>
) were identified, as only the middle part of the left lobe of the liver was scanned. The mean TRE was 2.99 mm and the maximum TRE was 4.37 mm.<fig id="Fig10"><label>Fig. 10</label>
<caption><p>Hepatic vein landmark positions. <bold>a</bold>
Eight bifurcation landmarks on the centre-line model were used to measure TRE in the first study. <bold>b</bold>
Three bifurcation landmarks on the centre-line model were used to measure TRE in the second study</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig10_HTML" id="MO13"></graphic>
</fig>
</p>
</sec>
<sec id="Sec17"><title>Navigation accuracy: in vivo</title>
<p>To evaluate the navigation accuracy, the surgeon scanned another LUS image sequence for each study (giving four US data sets in total), again using minimal force on the LUS probe to avoid deformation. Using the same bifurcation landmarks as in the registration experiment (section “Registration accuracy: in vivo”), the corresponding landmarks on LUS images were manually identified. For the first study, the mean TRE was 4.48 mm and the maximum TRE was 7.18 mm. For the second study, the mean TRE was 3.71 mm and the maximum TRE was 4.40 mm.</p>
</sec>
<sec id="Sec18"><title>Comparison of insufflated versus non-insufflated models</title>
<p>In sections “Registration accuracy: in vivo” and “Navigation accuracy: in vivo”, we used the insufflated CT model to evaluate the registration and navigation accuracy. In clinical practice, the patient would be scanned without insufflation, so in this section we used vessel centre lines derived from both insufflated and non-insufflated CT data. From the first study, landmarks 1, 2, 4, 5 (see Fig. <xref rid="Fig10" ref-type="fig">10</xref>
a) were manually identified and labelled in both the US images and the CT data. From the second study, landmarks 1, 2, 4 (see Fig. <xref rid="Fig10" ref-type="fig">10</xref>
b) were used. Using each landmark a registration was performed, registering the CT data to the US using the manual registration method (a landmark and two vectors, illustrated in Fig. <xref rid="Fig6" ref-type="fig">6</xref>
a, b). For each registration the TRE was evaluated as in section “Registration accuracy: in vivo” using the eight bifurcations for the first study and the three bifurcations for the second study. The measures of TRE are presented graphically in Fig. <xref rid="Fig11" ref-type="fig">11</xref>
. Similarly the navigation error is measured on the second LUS sequence for each study for each locally rigid registration. The measures of navigation error are illustrated in Fig. <xref rid="Fig12" ref-type="fig">12</xref>
.<fig id="Fig11"><label>Fig. 11</label>
<caption><p>Evaluation of registration accuracy with locally rigid registration. The errors are shown as a function of distance from the landmark used to register. Within 35-mm distance to the reference points, 76 % landmarks have TRE smaller or equal to 10 mm with the insufflated CT model; 72 % for the non-insufflated CT model</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig11_HTML" id="MO14"></graphic>
</fig>
<fig id="Fig12"><label>Fig. 12</label>
<caption><p>Evaluation of navigation accuracy with locally rigid registration. The errors are shown as a function of distance from the reference landmarks. Within 35-mm distance to the reference points, 74 % landmarks have TRE smaller or equal to 10 mm with the insufflated CT model; 71 % for the non-insufflated CT model</p>
</caption>
<graphic xlink:href="11548_2015_1236_Fig12_HTML" id="MO15"></graphic>
</fig>
</p>
</sec>
</sec>
<sec id="Sec19"><title>Discussion</title>
<p>In this paper, we describe and evaluate a practical laparoscopic image guidance system based on a fast and accurate vessel centre-point reconstruction coupled with a locally rigid registration to the preoperative model using vascular features visible in LUS.</p>
<p>In section “Ultrasound calibration error”, we checked the accuracy of the invariant point calibration method. The mean edge length between pins in the 8-pin phantom was 24.62 mm compared with a manufactured edge length of 25 mm. Table <xref rid="Tab1" ref-type="table">1</xref>
shows the reconstructed physical position errors between 0.81 and 3.40 mm, and an average of 2.17 mm, and this includes errors in measuring the gold standard itself. We concluded that although simple, this was comparable to other methods [<xref ref-type="bibr" rid="CR17">17</xref>
]. Future work could try improve calibration accuracy, specifically for LUS probes.</p>
<p>We subsequently checked the segmentation accuracy on a phantom (section “Vessel segmentation error”). The phantom was constructed via 3D printing a computer-aided design (CAD) model. So, the plastic phantom had known geometry with a tolerance of 0.1 mm. The reconstructed size of the internal diameter of the tubes was 6.4 mm compared with the diameter in the CAD model of 6.5 mm and was deemed within tolerance. Furthermore, in section “Registration accuracy: phantom” we see that the ICP-based registration of the point cloud resulting from the US segmentation to the CAD model itself gave a RMS error of 0.7 mm.</p>
<p>In section “Registration accuracy: in vivo”, we evaluated the registration accuracy in two in vivo studies. The mean TRE was 3.58 and 2.99 mm, measured at eight and three identifiable landmarks,respectively. This represents a best-case scenario for rigid registration, as we used an insufflated CT model, and a large region of interest (left temporal lobe). However, it does include movement due to respiration and cardiac pulsatile motion. The controlled breathing means that most of the time is spent near maximum exhale. We collected data for around 40 seconds, over several breathing cycles. Thus we assume that for the ICP-based methods, over a large region of interest, the data will be somewhat noisy, but the registration will average over the noise. For the manual landmark- based registration, future work will consider breath-holding techniques, faster software or a footswitch synchronised to the breathing. During the cardiac cycle, vessels pulsate and change size. We mitigated against this problem by using vessel centre lines which should be more reliable than vessel external contours.</p>
<p>From the initial registration, a second test data set was used to evaluate navigation accuracy. This incorporates the error due to registration, additional nonlinear EM tracking errors and errors due to further liver deformation via the US probe. Comparing the TRE errors of the corresponding data set in sections “Registration accuracy: in vivo” and “Navigation accuracy: in vivo”, the navigation accuracy is only slightly worse than the registration accuracy, if the surgeon performed the US scans in a consistent way. This also suggests the EM tracking error may not be a major problem, although further work is needed here.</p>
<p>In clinical practice, the patient will not be CT scanned while insufflated. The preoperative, non-insufflated CT will have a significantly different shape to that seen during surgery. So we compared registration of both insufflated and non-insufflated CT. It was difficult to identify corresponding landmarks in both CT scans. So rather than having eight landmarks in study 1, we could only identify landmarks labelled as 1, 2, 4 and 5 in Fig. <xref rid="Fig10" ref-type="fig">10</xref>
a consistently in both insufflated and non-insufflated CT models. If a large region of interest was scanned using the US probe, the ICP-based registration to non-insufflated CT models was unreliable, due to the significantly different shape. If a small region of interest was scanned, then the smaller the structure, the more likely it was to be featureless, e.g. more closely resembling a line. So, to directly compare insufflated with non-insufflated registration, the manual landmark- based method (section “Registration”) was used around individual bifurcations, so as to be consistent across the two studies. Comparing Figs. <xref rid="Fig11" ref-type="fig">11</xref>
and <xref rid="Fig12" ref-type="fig">12</xref>
, we can see that there are similar errors when using non-insufflated or insufflated errors. But an acceptable level (<inline-formula id="IEq36"><alternatives><tex-math id="M77">\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$<$$\end{document}</tex-math>
<mml:math id="M78"><mml:mo><</mml:mo>
</mml:math>
<inline-graphic xlink:href="11548_2015_1236_Article_IEq36.gif"></inline-graphic>
</alternatives>
</inline-formula>
5 mm) is achievable only near to a registration point. Interestingly, the navigation errors are similar. We tested locally rigid registrations on both insufflated and non-insufflated CT models which gave mean (standard deviation) errors of 4.23 (2.18) mm and 6.57 (3.41) mm, when measured at target landmarks located within 10 mm of the landmark used to register. When measured within 35 mm to the reference points, over 70 % of the target landmarks have errors smaller or equal to 10 mm for both model. Figures <xref rid="Fig11" ref-type="fig">11</xref>
and <xref rid="Fig12" ref-type="fig">12</xref>
show that if TREs are assessed away from the reference points, then errors do indeed increase.</p>
<p>Considering the state of the art in finite element methods that do attempt to compensate for tissue deformation, Suwelack et al. [<xref ref-type="bibr" rid="CR23">23</xref>
] measured errors of 5.05 mm and 8.7 mm on a liver phantom, Haouchine et al. [<xref ref-type="bibr" rid="CR13">13</xref>
] measure registration accuracy at two points as 2.2 and 5.3, in an ex-vivo trial while Bano et al. [<xref ref-type="bibr" rid="CR3">3</xref>
] measure 4 mm error at the liver surface but 10 mm error at structures internal to the liver. Although deformable models based on understanding of the biomechanics of tissue deformation are developing rapidly [<xref ref-type="bibr" rid="CR3">3</xref>
, <xref ref-type="bibr" rid="CR13">13</xref>
, <xref ref-type="bibr" rid="CR23">23</xref>
], there remain significant issues of validation in a surgical environment. We anticipate that it will be a long time before the surgeon has sufficient faith in a deforming model alone to guide surgical decisions during resection itself. However we do propose that this locally rigid registration system is practical and could relatively easily be automated with minimal user intervention. Local regions could also be used to drive and validate a deformable model.</p>
</sec>
<sec id="Sec20"><title>Conclusion</title>
<p>While this work is preliminary and performed on two studies from a single porcine experiment, we believe that the results provide preliminary evidence that our method is sufficiently accurate to be further developed and validated using animal models and clinically. Given a simple interface and a sufficiently close initial estimate, the liver could be scanned round the target lesion and nearby vessel bifurcations, and then it may be possible to obtain registration errors of the order of 4–6 mm with no deformable modelling. Our proposed method is both practical and provides guidance to the surgical target. It also implicitly includes information on the location of nearby vasculature structures which are the same structures that the surgeon needs to be aware of when undertaking laparoscopic resection. It may be that such a system has advantages over open surgery and haptics where the surgeon still remains blind to the precise location of these structures. Future work, considering either the vessel centre lines and deformable registration or combining ultrasound vessel centre lines into a deformable model will likely provide either better accuracy, more robustness or both.</p>
</sec>
</body>
<back><fn-group><fn id="Fn1"><label>1</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.cascination.com">http://www.cascination.com</ext-link>
.</p>
</fn>
<fn id="Fn2"><label>2</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.visiblepatient.com">http://www.visiblepatient.com</ext-link>
.</p>
</fn>
<fn id="Fn3"><label>3</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.vmtk.org/tutorials/Centerlines.html">http://www.vmtk.org/tutorials/Centerlines.html</ext-link>
.</p>
</fn>
<fn id="Fn4"><label>4</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.analogicultrasound.com">http://www.analogicultrasound.com</ext-link>
.</p>
</fn>
<fn id="Fn5"><label>5</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.vermon.com">http://www.vermon.com</ext-link>
.</p>
</fn>
<fn id="Fn6"><label>6</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.ascension-tech.com">http://www.ascension-tech.com</ext-link>
.</p>
</fn>
</fn-group>
<ack><p>This publication presents independent research funded by the Health Innovation Challenge Fund (HICF-T4-317), a parallel funding partnership between the Wellcome Trust and the Department of Health. The views expressed in this publication are those of the author(s) and not necessarily those of the Wellcome Trust or the Department of Health. DB and DJH received funding from EPSRC EP/F025750/1. SO and DJH receive funding from EPSRC EP/H046410/1 and the National Institute for Health Research (NIHR) University College London Hospitals Biomedical Research Centre (BRC) High Impact Initiative. We would like to thank NVidia Corporation for the donation of the Quadro K5000 and SDI capture cards used in this research.</p>
</ack>
<notes notes-type="conflict-interest"><sec id="FPar1"><title>Conflict of interest</title>
<p>The authors declare that they have no conflict of interest.</p>
</sec>
<sec id="FPar2"><title>Ethical standard</title>
<p>All applicable international, national and/or institutional guidelines for the care and use of animals were followed. All procedures performed in studies involving animals were in accordance with the ethical standards of the institution or practice at which the studies were conducted.</p>
</sec>
</notes>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<mixed-citation publication-type="other">Anderegg S, Peterhans M, Weber S (2010) Ultrasound segmentation in navigated liver surgery. <ext-link ext-link-type="uri" xlink:href="http://www.cascination.com/information/publications/">http://www.cascination.com/information/publications/</ext-link>
</mixed-citation>
</ref>
<ref id="CR2"><label>2.</label>
<mixed-citation publication-type="other">Aylward SR, Jomier J, Guyon JP, Weeks S (2002) Intra-operative 3D ultrasound augmentation. In: Proceedings, 2002 IEEE international symposium on biomedical imaging, pp 421–424. IEEE. doi:10.1109/ISBI.2002.1029284</mixed-citation>
</ref>
<ref id="CR3"><label>3.</label>
<mixed-citation publication-type="other">Bano J, Nicolau S, Hostettler A, Doignon C, Marescaux J, Soler L (2013) Registration of preoperative liver model for laparoscopic surgery from intraoperative 3d acquisition. In: Liao H, Linte C, Masamune K, Peters T, Zheng G (eds) Augmented reality environments for medical imaging and computer-assisted interventions. Lecture notes in computer science, vol 8090. Springer, Berlin, pp 201–210. doi:10.1007/978-3-642-40843-4_22</mixed-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Barratt</surname>
<given-names>DC</given-names>
</name>
<name><surname>Davies</surname>
<given-names>AH</given-names>
</name>
<name><surname>Hughes</surname>
<given-names>AD</given-names>
</name>
<name><surname>Thom</surname>
<given-names>SA</given-names>
</name>
<name><surname>Humphries</surname>
<given-names>KN</given-names>
</name>
</person-group>
<article-title>Accuracy of an electromagnetic three-dimensional ultrasound system for carotid artery imaging</article-title>
<source>Ultrasound Med Biol</source>
<year>2001</year>
<volume>27</volume>
<issue>10</issue>
<fpage>1421</fpage>
<lpage>1425</lpage>
<pub-id pub-id-type="doi">10.1016/S0301-5629(01)00447-1</pub-id>
<pub-id pub-id-type="pmid">11731055</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<mixed-citation publication-type="other">Besl PJ, McKay ND (1992) Method for registration of 3-D shapes. In: Robotics-DL tentative. International Society for Optics and Photonics, pp 586–606. doi:10.1117/12.57955</mixed-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Clarkson</surname>
<given-names>M</given-names>
</name>
<name><surname>Zombori</surname>
<given-names>G</given-names>
</name>
<name><surname>Thompson</surname>
<given-names>S</given-names>
</name>
<name><surname>Totz</surname>
<given-names>J</given-names>
</name>
<name><surname>Song</surname>
<given-names>Y</given-names>
</name>
<name><surname>Espak</surname>
<given-names>M</given-names>
</name>
<name><surname>Johnsen</surname>
<given-names>S</given-names>
</name>
<name><surname>Hawkes</surname>
<given-names>D</given-names>
</name>
<name><surname>Ourselin</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>The NifTK software platform for image-guided interventions: platform overview and NiftyLink messaging</article-title>
<source>Int J Comput Assist Radiol Surg</source>
<year>2015</year>
<volume>10</volume>
<issue>3</issue>
<fpage>301</fpage>
<lpage>316</lpage>
<pub-id pub-id-type="doi">10.1007/s11548-014-1124-7</pub-id>
<pub-id pub-id-type="pmid">25408304</pub-id>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Croome</surname>
<given-names>KP</given-names>
</name>
<name><surname>Yamashita</surname>
<given-names>MH</given-names>
</name>
</person-group>
<article-title>Laparoscopic vs open hepatic resection for benign and malignant tumors: an updated meta-analysis</article-title>
<source>Arch Surg</source>
<year>2010</year>
<volume>145</volume>
<issue>11</issue>
<fpage>1109</fpage>
<lpage>1118</lpage>
<pub-id pub-id-type="doi">10.1001/archsurg.2010.227</pub-id>
<pub-id pub-id-type="pmid">21079101</pub-id>
</element-citation>
</ref>
<ref id="CR8"><label>8.</label>
<mixed-citation publication-type="other">Dagon B, Baur C, Bettschart V (2008) Real-time update of 3D deformable models for computer aided liver surgery. In: 19th international conference on pattern recognition (ICPR 2008), pp. 1–4. IEEE. doi:10.1109/ICPR.2008.4761741</mixed-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Feuerstein</surname>
<given-names>M</given-names>
</name>
<name><surname>Reichl</surname>
<given-names>T</given-names>
</name>
<name><surname>Vogel</surname>
<given-names>J</given-names>
</name>
<name><surname>Traub</surname>
<given-names>J</given-names>
</name>
<name><surname>Navab</surname>
<given-names>N</given-names>
</name>
</person-group>
<article-title>Magneto-optical tracking of flexible laparoscopic ultrasound: model-based online detection and correction of magnetic tracking errors</article-title>
<source>IEEE Trans Med Imaging</source>
<year>2009</year>
<volume>28</volume>
<issue>6</issue>
<fpage>951</fpage>
<lpage>967</lpage>
<pub-id pub-id-type="doi">10.1109/TMI.2008.2008954</pub-id>
<pub-id pub-id-type="pmid">19211352</pub-id>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<mixed-citation publication-type="other">Frangi A, Niessen W, Vincken K, Viergever M (1998) Multiscale vessel enhancement filtering. In: Wells W, Colchester A, Delp S (eds) Medical image computing and computer-assisted interventation MICCAI98. Lecture notes in computer science, vol 1496. Springer, Berlin, pp 130–137. doi:10.1007/BFb0056195</mixed-citation>
</ref>
<ref id="CR11"><label>11.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Franz</surname>
<given-names>A</given-names>
</name>
<name><surname>Haidegger</surname>
<given-names>T</given-names>
</name>
<name><surname>Birkfellner</surname>
<given-names>W</given-names>
</name>
<name><surname>Cleary</surname>
<given-names>K</given-names>
</name>
<name><surname>Peters</surname>
<given-names>T</given-names>
</name>
<name><surname>Maier-Hein</surname>
<given-names>L</given-names>
</name>
</person-group>
<article-title>Electromagnetic tracking in medicine 2014: a review of technology, validation, and applications</article-title>
<source>IEEE Trans Med Imaging</source>
<year>2014</year>
<volume>33</volume>
<issue>8</issue>
<fpage>1702</fpage>
<lpage>1725</lpage>
<pub-id pub-id-type="doi">10.1109/TMI.2014.2321777</pub-id>
<pub-id pub-id-type="pmid">24816547</pub-id>
</element-citation>
</ref>
<ref id="CR12"><label>12.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guerrero</surname>
<given-names>J</given-names>
</name>
<name><surname>Salcudean</surname>
<given-names>S</given-names>
</name>
<name><surname>McEwen</surname>
<given-names>J</given-names>
</name>
<name><surname>Masri</surname>
<given-names>B</given-names>
</name>
<name><surname>Nicolaou</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Real-time vessel segmentation and tracking for ultrasound imaging applications</article-title>
<source>IEEE Trans Med Imaging</source>
<year>2007</year>
<volume>26</volume>
<issue>8</issue>
<fpage>1079</fpage>
<lpage>1090</lpage>
<pub-id pub-id-type="doi">10.1109/TMI.2007.899180</pub-id>
<pub-id pub-id-type="pmid">17695128</pub-id>
</element-citation>
</ref>
<ref id="CR13"><label>13.</label>
<mixed-citation publication-type="other">Haouchine N, Dequidt J, Peterlik I, Kerrien E, Berger MO, Cotin S (2013) Image-guided simulation of heterogeneous tissue deformation for augmented reality during hepatic surgery. In: 2013 IEEE international symposium on mixed and augmented reality (ISMAR), pp 199–208. doi:10.1109/ISMAR.2013.6671780</mixed-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kingham</surname>
<given-names>TP</given-names>
</name>
<name><surname>Jayaraman</surname>
<given-names>S</given-names>
</name>
<name><surname>Clements</surname>
<given-names>LW</given-names>
</name>
<name><surname>Scherer</surname>
<given-names>MA</given-names>
</name>
<name><surname>Stefansic</surname>
<given-names>JD</given-names>
</name>
<name><surname>Jarnagin</surname>
<given-names>WR</given-names>
</name>
</person-group>
<article-title>Evolution of image-guided liver surgery: transition from open to laparoscopic procedures</article-title>
<source>J Gastrointest Surg</source>
<year>2013</year>
<volume>17</volume>
<issue>7</issue>
<fpage>1274</fpage>
<lpage>1282</lpage>
<pub-id pub-id-type="doi">10.1007/s11605-013-2214-5</pub-id>
<pub-id pub-id-type="pmid">23645420</pub-id>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lange</surname>
<given-names>T</given-names>
</name>
<name><surname>Eulenstein</surname>
<given-names>S</given-names>
</name>
<name><surname>Hünerbein</surname>
<given-names>M</given-names>
</name>
<name><surname>Schlag</surname>
<given-names>PM</given-names>
</name>
</person-group>
<article-title>Vessel-based non-rigid registration of MR/CT and 3D ultrasound for navigation in liver surgery</article-title>
<source>Comput Aided Surg</source>
<year>2003</year>
<volume>8</volume>
<issue>5</issue>
<fpage>228</fpage>
<lpage>240</lpage>
<pub-id pub-id-type="doi">10.3109/10929080309146058</pub-id>
<pub-id pub-id-type="pmid">15529952</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lange</surname>
<given-names>T</given-names>
</name>
<name><surname>Papenberg</surname>
<given-names>N</given-names>
</name>
<name><surname>Heldmann</surname>
<given-names>S</given-names>
</name>
<name><surname>Modersitzki</surname>
<given-names>J</given-names>
</name>
<name><surname>Fischer</surname>
<given-names>B</given-names>
</name>
<name><surname>Lamecker</surname>
<given-names>H</given-names>
</name>
<name><surname>Schlag</surname>
<given-names>PM</given-names>
</name>
</person-group>
<article-title>3D ultrasound-CT registration of the liver using combined landmark-intensity information</article-title>
<source>Int J Comput Assist Radiol Surg</source>
<year>2009</year>
<volume>4</volume>
<issue>1</issue>
<fpage>79</fpage>
<lpage>88</lpage>
<pub-id pub-id-type="doi">10.1007/s11548-008-0270-1</pub-id>
<pub-id pub-id-type="pmid">20033605</pub-id>
</element-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mercier</surname>
<given-names>L</given-names>
</name>
<name><surname>Lang</surname>
<given-names>T</given-names>
</name>
<name><surname>Lindseth</surname>
<given-names>F</given-names>
</name>
<name><surname>Collins</surname>
<given-names>LD</given-names>
</name>
</person-group>
<article-title>A review of calibration techniques for freehand 3-d ultrasound systems</article-title>
<source>Ultrasound Med Biol</source>
<year>2005</year>
<volume>31</volume>
<issue>2</issue>
<fpage>143</fpage>
<lpage>165</lpage>
<pub-id pub-id-type="doi">10.1016/j.ultrasmedbio.2004.11.001</pub-id>
<pub-id pub-id-type="pmid">15708453</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<mixed-citation publication-type="other">Nakada K, Nakamoto M, Sato Y, Konishi K, Hashizume M, Tamura S (2003) A rapid method for magnetic tracker calibration using a magneto-optic hybrid tracker. In: Ellis R, Peters T (eds) Medical image computing and computer-assisted intervention–MICCAI 2003. Lecture notes in computer science, vol 2879. Springer, Berlin, pp 285–293. doi:10.1007/978-3-540-39903-2_36</mixed-citation>
</ref>
<ref id="CR19"><label>19.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nicolau</surname>
<given-names>S</given-names>
</name>
<name><surname>Soler</surname>
<given-names>L</given-names>
</name>
<name><surname>Mutter</surname>
<given-names>D</given-names>
</name>
<name><surname>Marescaux</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Augmented reality in laparoscopic surgical oncology</article-title>
<source>Surg Oncol</source>
<year>2011</year>
<volume>20</volume>
<issue>3</issue>
<fpage>189</fpage>
<lpage>201</lpage>
<pub-id pub-id-type="doi">10.1016/j.suronc.2011.07.002</pub-id>
<pub-id pub-id-type="pmid">21802281</pub-id>
</element-citation>
</ref>
<ref id="CR20"><label>20.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Noble</surname>
<given-names>J</given-names>
</name>
<name><surname>Boukerroui</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Ultrasound image segmentation: a survey</article-title>
<source>IEEE Trans Med Imaging</source>
<year>2006</year>
<volume>25</volume>
<issue>8</issue>
<fpage>987</fpage>
<lpage>1010</lpage>
<pub-id pub-id-type="doi">10.1109/TMI.2006.877092</pub-id>
<pub-id pub-id-type="pmid">16894993</pub-id>
</element-citation>
</ref>
<ref id="CR21"><label>21.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Penney</surname>
<given-names>GP</given-names>
</name>
<name><surname>Blackall</surname>
<given-names>JM</given-names>
</name>
<name><surname>Hamady</surname>
<given-names>M</given-names>
</name>
<name><surname>Sabharwal</surname>
<given-names>T</given-names>
</name>
<name><surname>Adam</surname>
<given-names>A</given-names>
</name>
<name><surname>Hawkes</surname>
<given-names>DJ</given-names>
</name>
</person-group>
<article-title>Registration of freehand 3D ultrasound and magnetic resonance liver images</article-title>
<source>Med Image Anal</source>
<year>2004</year>
<volume>8</volume>
<issue>1</issue>
<fpage>81</fpage>
<lpage>91</lpage>
<pub-id pub-id-type="doi">10.1016/j.media.2003.07.003</pub-id>
<pub-id pub-id-type="pmid">14644148</pub-id>
</element-citation>
</ref>
<ref id="CR22"><label>22.</label>
<mixed-citation publication-type="other">Schneider C, Guerrero J, Nguan C, Rohling R, Salcudean S (2011) Intra-operative pick-up ultrasound for robot assisted surgery with vessel extraction and registration: a feasibility study. In: Taylor R, Yang GZ (eds) Information processing in computer-assisted interventions. Lecture notes in computer science, vol 6689. Springer, Berlin, pp 122–132. doi:10.1007/978-3-642-21504-9_12</mixed-citation>
</ref>
<ref id="CR23"><label>23.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Suwelack</surname>
<given-names>S</given-names>
</name>
<name><surname>Rhl</surname>
<given-names>S</given-names>
</name>
<name><surname>Bodenstedt</surname>
<given-names>S</given-names>
</name>
<name><surname>Reichard</surname>
<given-names>D</given-names>
</name>
<name><surname>Dillmann</surname>
<given-names>R</given-names>
</name>
<name><surname>dos Santos</surname>
<given-names>T</given-names>
</name>
<name><surname>Maier-Hein</surname>
<given-names>L</given-names>
</name>
<name><surname>Wagner</surname>
<given-names>M</given-names>
</name>
<name><surname>Wnscher</surname>
<given-names>J</given-names>
</name>
<name><surname>Kenngott</surname>
<given-names>H</given-names>
</name>
<name><surname>Mller</surname>
<given-names>BP</given-names>
</name>
<name><surname>Speidel</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Physics-based shape matching for intraoperative image guidance</article-title>
<source>Med Phys</source>
<year>2014</year>
<volume>41</volume>
<issue>11</issue>
<fpage>111901</fpage>
<pub-id pub-id-type="doi">10.1118/1.4896021</pub-id>
<pub-id pub-id-type="pmid">25370634</pub-id>
</element-citation>
</ref>
<ref id="CR24"><label>24.</label>
<mixed-citation publication-type="other">Thompson S, Totz J, Song Y, Stoyanov D, Ourselin S, Hawkes DJ, Clarkson MJ (2015) Accuracy validation of an imageguided laparoscopy system for liver resection. In: Proceedings of SPIE medical imaging</mixed-citation>
</ref>
<ref id="CR25"><label>25.</label>
<mixed-citation publication-type="other">Totz J, Thompson S, Stoyanov D, Gurusamy K, Davidson B, Hawkes DJ, Clarkson MJ (2014) Fast semi-dense surface reconstruction from stereoscopic video in laparoscopic surgery. In: Stoyanov D, Collins D, Sakuma I, Abolmaesumi P, Jannin P (eds) Information processing in computer-assisted interventions. Lecture notes in computer science, vol 8498. Springer, pp 206–215. doi:10.1007/978-3-319-07521-1_22</mixed-citation>
</ref>
<ref id="CR26"><label>26.</label>
<mixed-citation publication-type="other">Wein W, Brunke S, Khamene A, Callstrom MR, Navab N (2008) Automatic CT-ultrasound registration for diagnostic imaging and image-guided intervention. Med Image Anal 12(5):577–585. doi:10.1016/j.media.2008.06.006</mixed-citation>
</ref>
<ref id="CR27"><label>27.</label>
<mixed-citation publication-type="other">Wein W, Ladikos A, Fuerst B, Shah A, Sharma K, Navab N (2013) Global registration of ultrasound to mri using the LC2 metric for enabling neurosurgical guidance. In: Mori K, Sakuma I, Sato Y, Barillot C, Navab N (eds) Medical image computing and computer-assisted intervention-MICCAI 2013. Springer, Berlin, Heidelberg, pp 34–41. doi:10.1007/978-3-642-40811-3_5</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/HapticV1/Data/Pmc/Curation
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000819 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Curation/biblio.hfd -nk 000819 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Ticri/CIDE |area= HapticV1 |flux= Pmc |étape= Curation |type= RBID |clé= PMC:4642598 |texte= Locally rigid, vessel-based registration for laparoscopic liver surgery }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Curation/RBID.i -Sk "pubmed:26092658" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Curation/biblio.hfd \ | NlmPubMed2Wicri -a HapticV1
![]() | This area was generated with Dilib version V0.6.23. | ![]() |