<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Verdana;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">You may also give elastix registration toolbox a try (<a href="http://elastix.isi.uu.nl/">http://elastix.isi.uu.nl/</a>). It is based on ITK but it adds a number of registration classes. My experience is that with BRAINS you often need
to tune parameters to get meaningful results, while elastix gives great results with default settings for all kinds of input data. Elastix is available in 3D Slicer, too (by installing SlicerElastix extension).<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Andras<o:p></o:p></p>
<p class="MsoNormal"><a name="_MailEndCompose"><o:p> </o:p></a></p>
<span style="mso-bookmark:_MailEndCompose"></span>
<p class="MsoNormal"><b>From:</b> Community [mailto:community-bounces@itk.org] <b>
On Behalf Of </b>Dženan Zukic<br>
<b>Sent:</b> Thursday, July 27, 2017 10:22 AM<br>
<b>To:</b> thanos thanos <thanosxania@gmail.com><br>
<b>Cc:</b> Insight-users <insight-users@itk.org><br>
<b>Subject:</b> Re: [ITK] [ITK-users] Cardiac Deformation using Segmentation and Registration of Ultrasound Images<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif">Hi Thanos,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif">segmentation is not a trivial task, and usually requires a lot of user intervention/interaction to achieve a satisfactory result.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif">Completely segmenting one case, and then registering other cases to the segmented one and transferring labels is called atlas segmentation.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif">To figure out proper registration parameters, you could play with BRAINSFit module of 3D Slicer (<a href="https://na01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.slicer.org%2Fwiki%2FModules%3ABRAINSFit&data=02%7C01%7Classo%40queensu.ca%7C247a47ea85324d5a764c08d4d4fae8f1%7Cd61ecb3b38b142d582c4efb2838b925c%7C1%7C0%7C636367621488662816&sdata=k94xFtCvTDijsGhFfuDeRSSkeFeQFobaTtVJAOElxzE%3D&reserved=0">docs</a>),
which is built on top of ITK and provides a nice user interface with significant number of options to tweak.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif">For B-Splines, common number of grid points are up to maybe 13x13x13. Since you have a high frame rate, you would probably benefit from initializing registration of time point
n+1 by resulting transform from time point n.<o:p></o:p></span></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif">Regards,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Verdana",sans-serif">Dženan Zukić, PhD, Senior R&D Engineer, Kitware (Carrboro, N.C.)<o:p></o:p></span></p>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">On Thu, Jul 27, 2017 at 8:08 AM, thanos thanos <<a href="mailto:thanosxania@gmail.com" target="_blank">thanosxania@gmail.com</a>> wrote:<o:p></o:p></p>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm">
<div>
<p class="MsoNormal">Hello everyone,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">The past months I have been working on tasks of Segmentation and Registration of 4D ultrasound images so I can eventually do a proper Cardiac Deformation analysis. Unfortunately I am struggling on my own to solve these problems so I still
haven't found a proper workflow and this is why I am posting here in case someone has a suggestion or a recommendation.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">For the acquisition of my data, a pig heart was placed inside a water tank with the 4D transducer on the bottom of the tank (the apex of the heart was about 5 cm from the transducer) and with around 50 volumes per second I obtained images
of 120x90x80, that after masking I ended up with volumes of 75x70x52 (x,y,z).<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">So far, I use Curvature Anisotropic Diffusion using the ITK libraries in order to perform an automatic segmentation using Active Contours (with Matlab). I filtered my data, mostly to assist my 2D segmentation algorithm, which I perform
slice by slice on the z axis. I obtain a 3D model using the isosurfaces.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Then I tried both the Demons registration and the Bsplines registration on the 3D volumes. I tried both on the original (unsegmented) data and two segmented volumes, but I can not say that I am quite satisfied with any of the results.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">So my questions are the following : <br>
1) Is it necessary to perform segmentation before the registration since my heart is placed in a water tank, so there are no other organs around?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">2) Should the registration be done on the filtered or the unfiltered data?<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">3) I am not satisfied at all with my segmentation algorithm since active contours are quite sensitive to noise and because of the nature of the data I have to use different parameters for every 4-5 slices. Therefore I was thinking to perform
an almost manual segmentation on one volume and then use registration for the segmentation of the rest (if that is necessary). So I am now studying and trying to say if something like that is possible (and especially on ITK). If someone knows something about
it, I would be happy to hear.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">4) What would be a proper number of grid points for my 75x70x52 images, when using the B-Splines registration?<br>
<br>
I am sorry for the long post but I have been struggling for too long with this problem and I would really appreciate any kind of help. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Best regards,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Thanos<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
_____________________________________<br>
Powered by <a href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.kitware.com&data=02%7C01%7Classo%40queensu.ca%7C247a47ea85324d5a764c08d4d4fae8f1%7Cd61ecb3b38b142d582c4efb2838b925c%7C1%7C0%7C636367621488662816&sdata=Sk2Z82q8%2FgZKcqZjhOhaDmbQ8%2BJNOO5TmFcBR14k8ok%3D&reserved=0" target="_blank">
www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.kitware.com%2Fopensource%2Fopensource.html&data=02%7C01%7Classo%40queensu.ca%7C247a47ea85324d5a764c08d4d4fae8f1%7Cd61ecb3b38b142d582c4efb2838b925c%7C1%7C0%7C636367621488662816&sdata=40oO6LOBFfUf94oulLX9ghuWI7pqR5Lzqb855h5Z%2F84%3D&reserved=0" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.kitware.com%2Fproducts%2Fprotraining.php&data=02%7C01%7Classo%40queensu.ca%7C247a47ea85324d5a764c08d4d4fae8f1%7Cd61ecb3b38b142d582c4efb2838b925c%7C1%7C0%7C636367621488662816&sdata=Vy4lHKRJBZlEW0%2BHX8Jvv7es1ESBRLmYsv%2BIooNEs48%3D&reserved=0" target="_blank">http://www.kitware.com/products/protraining.php</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.itk.org%2FWiki%2FITK_FAQ&data=02%7C01%7Classo%40queensu.ca%7C247a47ea85324d5a764c08d4d4fae8f1%7Cd61ecb3b38b142d582c4efb2838b925c%7C1%7C0%7C636367621488662816&sdata=npBMPAyC79LtcBrkLcBrElHeWQ1awEhV9IrfMiFmxqw%3D&reserved=0" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fpublic.kitware.com%2Fmailman%2Flistinfo%2Finsight-users&data=02%7C01%7Classo%40queensu.ca%7C247a47ea85324d5a764c08d4d4fae8f1%7Cd61ecb3b38b142d582c4efb2838b925c%7C1%7C0%7C636367621488662816&sdata=XKhOVF5hK3yQzPcjREpUntiv%2BZzrnAU91galWezU3Bg%3D&reserved=0" target="_blank">http://public.kitware.com/mailman/listinfo/insight-users</a><o:p></o:p></p>
</blockquote>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</body>
</html>