So how can we help Dean?<br><br><div class="gmail_quote">On Sun, Oct 23, 2011 at 11:18 AM, Johnson, Hans J <span dir="ltr"><<a href="mailto:hans-johnson@uiowa.edu">hans-johnson@uiowa.edu</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;">
<div 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: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204); 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: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204); 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>
<div style="font-family: Calibri; font-size: 11pt; text-align: left; color: black; border-width: 1pt medium medium; border-style: solid none none; border-color: rgb(181, 196, 223) -moz-use-text-color -moz-use-text-color; padding: 3pt 0in 0in;">
<span style="font-weight: bold;">From: </span>Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com" target="_blank">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" target="_blank">dean.inglis@camris.ca</a>><br>
<span style="font-weight: bold;">Cc: </span><<a href="mailto:insight-users@itk.org" target="_blank">insight-users@itk.org</a>><br>
<span style="font-weight: bold;">Subject: </span>Re: [Insight-users] initialize VersorRigid3DTransform<br>
</div><div><div></div><div class="h5">
<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" target="_blank">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 <a href="http://www.kitware.com" target="_blank">www.kitware.com</a> Visit other Kitware open-source projects at
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">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" target="_blank">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" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a> Follow this link to subscribe/unsubscribe:
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a>
</div></div></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>
</div>
</blockquote></div><br>