[Insight-users] QuaternionRigidTransform SetRotation

[Luis Ibanez]

Hi Tony,

As we pointed out in a previous email:

You *DON'T NEED TO INVERSE THE TRANSFORM*.

When you are doing PointSet to Image registration
the transform that you compute is the one that maps
the points from the fixed point set into the coordinate
system of the image.

Therefore, once you finish the registration you simply
need to pass the same transform to the TransformMeshFilter.

The output of the filter will be a point set that is
registered to the image coordinate system.


    Regards,


        Luis


----------------------
tony hakki wrote:
> yes I am doing this for initializing the transform before registration. 
> By the way at the end of the Pointsettoimageregistration  i get 
> finalparameters,shouldn't I calculate these parameters opposite,because 
> these parameters belong to image ,i mean they indicates that how far 
> image will transform and how far it will rotate. But I would like to 
> apply these finalparameters to my pointset. For this reason I use 
> TransformMeshFilter.  When I get the inverse of the final transform (  
> *finalTransform2->GetInverse(finalTransform2);*)the new  translation 
> parameters  are calculated,but rotation parameters not renew? If I  am 
> right,is there a way to find the opposite rotation patameters.
>  
> 
> typedef itk::TransformMeshFilter<
> 
> FMeshType,
> 
> FMeshType,
> 
> TransformType > FilterType;
> 
>  
> 
> FilterType::Pointer filter = FilterType::New();
> 
> filter->SetInput( outputMesh);
> 
> filter->SetTransform(finalTransform );
> 
> filter->Update();
> 
> FMeshType::Pointer outputMesh2 = filter->GetOutput();
> 
>  
> 
> 
> 
> ----- Original Message ----
> From: Luis Ibanez <luis.ibanez at kitware.com>
> To: tony hakki <tony2007vtk at yahoo.com>
> Cc: insight-users at itk.org
> Sent: Saturday, March 31, 2007 7:33:28 PM
> Subject: Re: [Insight-users] QuaternionRigidTransform SetRotation
> 
> Hi Tony,
> 
> Are you doing this for Initializing the transform
> before starting the registration process ?
> 
> 
> Please look at the documentation of this class:
> http://www.itk.org/Insight/Doxygen/html/classitk_1_1QuaternionRigidTransform.html
> 
> You can call the SetRotationMatrix() method.
> 
> 
> 
> 
>     Regards,
> 
> 
>        Luis
> 
> 
> 
> --------------------
> tony hakki wrote:
>  > Hello;
>  > I am trying to implement pointsettoimageregistration and I use
>  > QuaternionRigidTransform class. How should I set the rotation of this
>  > transform?for example I set the translation like that  
>  > *finalTransform->SetTranslation( centerMovingimage - centerFixedpointset
>  > );*  translation works fine ,but i also have to set rotation  . how can
>  > I calculate the rotation matrix? I hope somebody will help me.
>  > Thank you very but very much
>  > Tony
>  >
>  > ------------------------------------------------------------------------
>  >
>  > _______________________________________________
>  > Insight-users mailing list
>  > Insight-users at itk.org
>  > http://www.itk.org/mailman/listinfo/insight-users
> 
> 
> ------------------------------------------------------------------------
> Never miss an email again!
> Yahoo! Toolbar 
> <http://us.rd.yahoo.com/evt=49938/*http://tools.search.yahoo.com/toolbar/features/mail/> 
> alerts you the instant new Mail arrives. Check it out. < 
> http://us.rd.yahoo.com/evt=49937/*http://tools.search.yahoo.com/toolbar/features/mail/>