<div dir="ltr"><div><div><div><div>Hello Matias,<br><br></div>I used the Euler3DTransform in a python terminal to see if I had problems setting the center and the angle, and everything seems to work as expected:<br><br>trans=itk.Euler3DTransform.D.New()<br>trans.SetCenter([0,0,0])<br></div>trans.SetRotation(10,0,0)<br></div>trans.TransformPoint([0,0,0])<br>Out[46]: itkPointD3 ([0, 0, 0]) <-expected as the transform center is [0,0,0]<br><br>trans.SetCenter([10,0,0])<br>trans.TransformPoint([0,0,0])<br>Out[48]: itkPointD3 ([5.16482, -7.96596, 3.62828]) <-rotation not around [0,0,0] but around [10,0,0]<br><br></div>To try to understand you problem, I would advise you to print the content of your transform to see if it is set as you expect. In my case, I get:<br><br><div>Euler3DTransform (0x18da8a0)<br> RTTI typeinfo: itk::Euler3DTransform<double><br> Reference Count: 1<br> Modified Time: 377<br> Debug: Off<br> Object Name: <br> Observers: <br> none<br> Matrix: <br> 0.483518 -0.362828 0.796596 <br> 0.796596 0.559609 -0.228631 <br> -0.362828 0.745113 0.559609 <br> Offset: [5.16482, -7.96596, 3.62828]<br> Center: [10, 0, 0]<br> Translation: [0, 0, 0]<br> Inverse: <br> 0.483518 0.796596 -0.362828 <br> -0.362828 0.559609 0.745113 <br> 0.796596 -0.228631 0.559609 <br> Singular: 0<br> Euler's angles: AngleX=40 AngleY=10 AngleZ=10<br> m_ComputeZYX = 0<br><br></div><div>Hope this helps,<br></div><div>Francois<br></div><div><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 17, 2017 at 7:11 PM, Matias <span dir="ltr"><<a href="mailto:matimontg@gmail.com" target="_blank">matimontg@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,
I'm trying to rotate a volume and use an user defined center but no matter what value I set, it seems the rotation is always performed using the same center (not sure which, maybe a default center)?
Here's the code I'm using.. am I missing something?
typedef itk::Euler3DTransform< double > TransformType; //Transform
TransformType::Pointer transform = TransformType::New();
double alfa, beta, gamma, centro_rotacion_X, centro_rotacion_Y, centro_rotacion_Z, origenX,origenY;
gamma = atof(argv[2]);
beta = atof(argv[3]);
alfa = atof(argv[4]);
centro_rotacion_X = atof(argv[5]);
centro_rotacion_Y = atof(argv[6]);
centro_rotacion_Z = atof(argv[7]);
double centro[3] = { centro_rotacion_X, centro_rotacion_Y, centro_rotacion_Z };
transform->SetCenter(centro);
transform->SetRotation(gamma, beta, alfa);
<br><hr width="300" align="left">
View this message in context: <a href="http://itk-users.7.n7.nabble.com/Euler3DTransform-not-taking-user-defined-Center-of-Rotation-tp38126.html" target="_blank">Euler3DTransform not taking user defined Center of Rotation?</a><br>
Sent from the <a href="http://itk-users.7.n7.nabble.com/" target="_blank">ITK - Users mailing list archive</a> at Nabble.com.<br><br>______________________________<wbr>_______<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>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" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>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" rel="noreferrer" target="_blank">http://www.itk.org/Wiki/ITK_<wbr>FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/insight-users" rel="noreferrer" target="_blank">http://public.kitware.com/<wbr>mailman/listinfo/insight-users</a><br>
<br></blockquote></div><br></div>