[ITK] [ITK-users] Inverse of Versor Rigid Transform and RIRE

Gabriel A. Giménez gabrielgimenez85 at gmail.com
Mon Apr 6 22:24:08 EDT 2015


Hello everyone again, hope you are well!

I am use RIRE project, specifically  CT (movig) and MR_PD (fixed) images.
Basically, I hava a set of point (in millimeters) of the CT image to which
apply the trasform result of the registration and updaload this results in
the web for the evaluation. Example of set of points and his "ground truth"
:

Point      x          y          z        new_x       new_y       new_z

  1       0.0000     0.0000     0.0000   -7.5731    -41.2534    -27.3093
  2     333.9870     0.0000     0.0000  324.8722    -72.8159    -32.9063
  3       0.0000   333.9870     0.0000   24.1607    291.0393    -16.2727
  4     333.9870   333.9870     0.0000  356.6060    259.4768    -21.8697
  5       0.0000     0.0000   112.0000   -6.0554    -45.1157     84.6137
  6     333.9870     0.0000   112.0000  326.3899    -76.6782     79.0168
  7       0.0000   333.9870   112.0000   25.6784    287.1769     95.6503
  8     333.9870   333.9870   112.0000  358.1237    255.6145     90.0534 

So, the first I need is the transformation to apply, for that I do the
following : 

   //get the inverse transform
   TransformType::Pointer inverseTransform = TransformType::New();
   inverseTransform->SetCenter( finalTransform->GetCenter() );
   bool response = finalTransform->GetInverse(inverseTransform);

It makes sense to use the same center in the inverse transform?. A
"quaternion" define an "axis" (right part) of rotation and an angle to use
for rotate the image about this axis...why use a center of rotation...? 

Second, apply this transform...as follows: 

NewPoint = RotationMatrix * OriginalPoint + Offset

The rotation matrix and the offset are obtained from the inverse transforme
objetc. Found something wrong? something that is not taking into account ?
The results do not appear to be correct...the calculated error is too big
and does not correspond with the visual result.

Thanks in advance ...








--
View this message in context: http://itk-insight-users.2283740.n2.nabble.com/Inverse-of-Versor-Rigid-Transform-and-RIRE-tp7587127.html
Sent from the ITK Insight Users mailing list archive at Nabble.com.
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.php

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/insight-users


More information about the Community mailing list