<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap:break-word; color:rgb(0,0,0); font-size:14px; font-family:Calibri,sans-serif">
<div>Bill,</div>
<div><br>
</div>
<div>In ITKv4 the "Centered" transforms are likely to be deprecated.</div>
<div><br>
</div>
<div>Brian Avants and Luis Ibanez provided the following:</div>
<div><br>
</div>
<div>
<blockquote class="gmail_quote" style="margin-top:0px; margin-right:0px; margin-bottom:0px; margin-left:0.8ex; border-left-width:1px; border-left-color:rgb(204,204,204); border-left-style:solid; padding-left:1ex; font-family:Calibri; font-size:medium">
<div style="word-wrap:break-word; color:rgb(0,0,0); font-size:14px; font-family:Calibri,sans-serif">
<ol>
<li>Rationale for removing the CenteredTransforms</li></ol>
</div>
</blockquote>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium">The goal of the centered transforms was to attempt to optimize both the center of rotation and the transform parameters. This leads to a poorly posed optimization problem because, as one changes the center
of rotation, the optimization landscape also changes leading to unstable performance. </div>
<div style="font-family:Calibri; font-size:medium"> </div>
<blockquote class="gmail_quote" style="margin-top:0px; margin-right:0px; margin-bottom:0px; margin-left:0.8ex; border-left-width:1px; border-left-color:rgb(204,204,204); border-left-style:solid; padding-left:1ex; font-family:Calibri; font-size:medium">
<div style="word-wrap:break-word; color:rgb(0,0,0); font-size:14px; font-family:Calibri,sans-serif">
<ol>
<li>A recommended ITKv4 upgrade path for those using the centered transform</li></ol>
</div>
</blockquote>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium">If you were using the centered transforms, we recommend replacing the centered transform with its non-centered equivalent. Then, where one previously wrote:</div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium">transform->SetCenter(center) </div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium">one would now write </div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium">transform->SetFixedParameters(fixed_parameters) </div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium">where the fixed_parameters will depend on which transform is being used. For matrix transforms, this will be the center of rotation (check this) and may be equal to the image's origin. For dense mappings
( BSpline, displacement fields ) , this will be the image spacing, direction, origin. </div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
</div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium">Hans</div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<div style="font-family:Calibri; font-size:medium"><br>
</div>
<span id="OLK_SRC_BODY_SECTION">
<div style="font-family:Calibri; font-size:11pt; text-align:left; color:black; border-bottom:medium none; border-left:medium none; padding-bottom:0in; padding-left:0in; padding-right:0in; border-top:#b5c4df 1pt solid; border-right:medium none; padding-top:3pt">
<span style="font-weight:bold">From: </span>Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>><br>
<span style="font-weight:bold">Date: </span>Sun, 23 Oct 2011 11:11:32 -0400<br>
<span style="font-weight:bold">To: </span>Dean Inglis <<a href="mailto:dean.inglis@camris.ca">dean.inglis@camris.ca</a>><br>
<span style="font-weight:bold">Cc: </span><<a href="mailto:insight-users@itk.org">insight-users@itk.org</a>><br>
<span style="font-weight:bold">Subject: </span>Re: [Insight-users] initialize VersorRigid3DTransform<br>
</div>
<div><br>
</div>
Is this class useful?<br>
./Registration/Common/include/itkCenteredVersorTransformInitializer<br>
<br>
<br>
<div class="gmail_quote">On Sat, Oct 22, 2011 at 12:44 PM, Dean Inglis <span dir="ltr">
<<a href="mailto:dean.inglis@camris.ca">dean.inglis@camris.ca</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex; border-left:1px solid rgb(204,204,204); padding-left:1ex">
how can I initialize a VersorRigid3DTransform for a point set<br>
to image registration problem using the following information:<br>
1) center of the moving points, c<br>
2) principal axes of the moving points, e1, e2, e3<br>
3) center of the image calculated using the ImageMomentsCalculator, cofg<br>
4) principal axes of the image calculated using the ImageMomentsCalculator, p1, p2, p3<br>
<br>
where all variables above are double arrays of size 3? Bascally, I would like<br>
to be able to initialize the transform:<br>
ParametersType parameters( transform-><u></u>GetNumberOfParameters() );<br>
<br>
// Versor type<br>
typedef TransformType::VersorType VersorType;<br>
VersorType versor;<br>
<br>
parameters[0] = versor.GetX(); // Rotations<br>
parameters[1] = versor.GetY();<br>
parameters[2] = versor.GetZ();<br>
parameters[3] = 0.0; // Translations<br>
parameters[4] = 0.0;<br>
parameters[5] = 0.0;<br>
<br>
<br>
transform->SetParameters( parameters );<br>
<br>
but I dont know how the rotations should be initialized. I could calculate a 3x3 rotation<br>
matrix :<br>
A = e * p<br>
where e = {e1,e2,e3} and p = {p1,p2,p3} but how could that be used to initialize the<br>
transform?<br>
<br>
Dean <br>
______________________________<u></u>_______<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/<u></u>opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/<u></u>products/protraining.html</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_<u></u>FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/<u></u>listinfo/insight-users</a><br>
</blockquote>
</div>
<br>
_____________________________________ Powered by www.kitware.com Visit other Kitware open-source projects at
<a href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a> Kitware offers ITK Training Courses, for more information visit:
<a href="http://www.kitware.com/products/protraining.html">http://www.kitware.com/products/protraining.html</a> Please keep messages on-topic and check the ITK FAQ at:
<a href="http://www.itk.org/Wiki/ITK_FAQ">http://www.itk.org/Wiki/ITK_FAQ</a> Follow this link to subscribe/unsubscribe:
<a href="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users</a>
</span><br>
<br>
<hr>
Notice: This UI Health Care e-mail (including attachments) is covered by the Electronic Communications Privacy Act, 18 U.S.C. 2510-2521, is confidential and may be legally privileged. If you are not the intended recipient, you are hereby notified that any
retention, dissemination, distribution, or copying of this communication is strictly prohibited. Please reply to the sender that you have received the message in error, then delete it. Thank you.
<hr>
</body>
</html>