<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:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",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;}
span.EmailStyle17
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
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"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">It seems that you want to do a simple landmark registration (you touch exactly the same points with the stylus, in the same order, as you defined in the image).
If this is the case, then you should not use ICP but </span><span style="font-size:9.0pt;font-family:Consolas;color:#333333;background:white">itk::LandmarkBasedTransformInitializer</span><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">
(with optional preprocessing step to determine the scale), as it is done in this Slicer module:<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><a href="https://github.com/Slicer/Slicer/blob/master/Modules/CLI/FiducialRegistration/FiducialRegistration.cxx">https://github.com/Slicer/Slicer/blob/master/Modules/CLI/FiducialRegistration/FiducialRegistration.cxx</a><o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Andras<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> Community [mailto:community-bounces@itk.org]
<b>On Behalf Of </b>Matias Montroull<br>
<b>Sent:</b> Wednesday, November 11, 2015 4:58 PM<br>
<b>To:</b> Matt McCormick <matt.mccormick@kitware.com><br>
<b>Cc:</b> insight-users@itk.org Users <insight-users@itk.org><br>
<b>Subject:</b> Re: [ITK] [ITK-users] Pointset registration<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Yes, here it goes..<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">you will notice there are points that I had to repeat, the original dataset was 4 points on each file and I was getting warnings: <i><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#262626"> Number of unknowns(12) greater
than number of data (4)</span></i><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">"puntosfijos.txt" is the fixed pointsets<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">"puntosflotantes.txt" is the moving pointset..<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks for your help.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Matias.<o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">El mié., 11 de nov. de 2015 a la(s) 6:43 p. m., Matt McCormick <<a href="mailto:matt.mccormick@kitware.com">matt.mccormick@kitware.com</a>> escribió:<o:p></o:p></p>
</div>
<blockquote style="border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-top:5.0pt;margin-right:0in;margin-bottom:5.0pt">
<p class="MsoNormal">Hi Matias,<br>
<br>
Can the datasets be shared?<br>
<br>
Thanks,<br>
Matt<br>
<br>
On Wed, Nov 11, 2015 at 4:38 PM, Matias Montroull <<a href="mailto:matimontg@gmail.com" target="_blank">matimontg@gmail.com</a>> wrote:<br>
> Thanks.<br>
><br>
> Here's another question, it seems I need 12 points to perform the<br>
> registration and obtain the matrix, now, the results I'm getting are not<br>
> very good. We have the algorithm programmed manually (not using ITK) and<br>
> when we compare the results are different (way too different).<br>
> I just pulled the example and modified to 3 dimensions and changed the<br>
> translationtransform to an affine transform. Is there an example using<br>
> Affine that I can check? I'm thinking is not just changing the transform,<br>
> there may be other variables I need to tweak..<br>
> Matias.<br>
><br>
> El mié., 11 de nov. de 2015 a la(s) 5:51 p. m., Dženan Zukić<br>
> <<a href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>> escribió:<br>
>><br>
>> Nowhere, which means it's the >> operator's standard white-space<br>
>> separators (spaces, tabs, newlines).<br>
>><br>
>> On Wed, Nov 11, 2015 at 3:03 PM, Matias Montroull <<a href="mailto:matimontg@gmail.com" target="_blank">matimontg@gmail.com</a>><br>
>> wrote:<br>
>>><br>
>>> by the way, where in the code is the space separator specified?<br>
>>><br>
>>> El mié., 11 de nov. de 2015 a la(s) 5:02 p. m., Matias Montroull<br>
>>> <<a href="mailto:matimontg@gmail.com" target="_blank">matimontg@gmail.com</a>> escribió:<br>
>>>><br>
>>>> that worked! Thanks for the advise, I was under the impression that<br>
>>>> points needed to be separated by ";"<br>
>>>> Thanks again<br>
>>>><br>
>>>> El mié., 11 de nov. de 2015 a la(s) 5:00 p. m., Dženan Zukić<br>
>>>> <<a href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>> escribió:<br>
>>>>><br>
>>>>> By looking at the code, the numbers should be separated by space, not<br>
>>>>> semicolon. Can you try that?<br>
>>>>><br>
>>>>> HTH<br>
>>>>><br>
>>>>> On Wed, Nov 11, 2015 at 2:51 PM, Matias Montroull <<a href="mailto:matimontg@gmail.com" target="_blank">matimontg@gmail.com</a>><br>
>>>>> wrote:<br>
>>>>>><br>
>>>>>> Hi Dženan, thanks for the suggestion, I tried that example and I think<br>
>>>>>> I'm missign something as my disk drive became crazy.<br>
>>>>>> I basically created 2 txt files with a set of 4 points each in this<br>
>>>>>> format:<br>
>>>>>> -0.585938;92.406328;38.200<br>
>>>>>> 44.7654;52.73442;61.000<br>
>>>>>> -53.5;49.218792;68.200<br>
>>>>>> -1.757814;81.445382;70.600<br>
>>>>>><br>
>>>>>> second file has the same format but different numbers of course.<br>
>>>>>><br>
>>>>>> Is there anything wrong with how I entered the point sets in both<br>
>>>>>> files?<br>
>>>>>><br>
>>>>>> As you can see it is 3 points per line separated by this character ";"<br>
>>>>>><br>
>>>>>> Thanks,<br>
>>>>>><br>
>>>>>> Matias.<br>
>>>>>><br>
>>>>>><br>
>>>>>> El mar., 10 de nov. de 2015 a la(s) 6:16 p. m., Dženan Zukić<br>
>>>>>> <<a href="mailto:dzenanz@gmail.com" target="_blank">dzenanz@gmail.com</a>> escribió:<br>
>>>>>>><br>
>>>>>>> Have you seen these two examples:<br>
>>>>>>><br>
>>>>>>> <a href="http://www.itk.org/Doxygen/html/Examples_2RegistrationITKv3_2IterativeClosestPoint3_8cxx-example.html" target="_blank">
http://www.itk.org/Doxygen/html/Examples_2RegistrationITKv3_2IterativeClosestPoint3_8cxx-example.html</a><br>
>>>>>>><br>
>>>>>>> <a href="http://www.itk.org/Doxygen/html/Examples_2RegistrationITKv4_2IterativeClosestPoint3_8cxx-example.html" target="_blank">
http://www.itk.org/Doxygen/html/Examples_2RegistrationITKv4_2IterativeClosestPoint3_8cxx-example.html</a><br>
>>>>>>><br>
>>>>>>> HTH<br>
>>>>>>><br>
>>>>>>> On Tue, Nov 10, 2015 at 3:18 PM, Matias Montroull<br>
>>>>>>> <<a href="mailto:matimontg@gmail.com" target="_blank">matimontg@gmail.com</a>> wrote:<br>
>>>>>>>><br>
>>>>>>>> Hi,<br>
>>>>>>>><br>
>>>>>>>> I have 2 sets of points (one is points on an image and another one<br>
>>>>>>>> is points on a NDI tracker) and I need to obtain the transformation matrix.<br>
>>>>>>>> What would be the best method? I'm currently using a Solver but I'd like to<br>
>>>>>>>> implement this using ITK and better if it is SimpleITK.<br>
>>>>>>>><br>
>>>>>>>> Thank you!<br>
>>>>>>>><br>
>>>>>>>> Matias.<br>
>>>>>>>> --<br>
>>>>>>>> Matias<br>
>>>>>>>><br>
>>>>>>>> _____________________________________<br>
>>>>>>>> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
>>>>>>>><br>
>>>>>>>> Visit other Kitware open-source projects at<br>
>>>>>>>> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">
http://www.kitware.com/opensource/opensource.html</a><br>
>>>>>>>><br>
>>>>>>>> Kitware offers ITK Training Courses, for more information visit:<br>
>>>>>>>> <a href="http://www.kitware.com/products/protraining.php" 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="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
>>>>>>>><br>
>>>>>>>> Follow this link to subscribe/unsubscribe:<br>
>>>>>>>> <a href="http://public.kitware.com/mailman/listinfo/insight-users" target="_blank">
http://public.kitware.com/mailman/listinfo/insight-users</a><br>
>>>>>>>><br>
>>>>>>><br>
>>>>>> --<br>
>>>>>> Matias<br>
>>>>><br>
>>>>><br>
>>>> --<br>
>>>> Matias<br>
>>><br>
>>> --<br>
>>> Matias<br>
>><br>
>><br>
> --<br>
> Matias<br>
><br>
> _____________________________________<br>
> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at<br>
> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Kitware offers ITK Training Courses, for more information visit:<br>
> <a href="http://www.kitware.com/products/protraining.php" 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="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://public.kitware.com/mailman/listinfo/insight-users" target="_blank">
http://public.kitware.com/mailman/listinfo/insight-users</a><br>
><o:p></o:p></p>
</blockquote>
</div>
</div>
<div>
<p class="MsoNormal">-- <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Matias<o:p></o:p></p>
</div>
</div>
</body>
</html>