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>