[ITK-users] Rotation of anisotropic 3D image

g2 antoine.letouzey at gmail.com
Wed Sep 13 07:16:44 EDT 2017


Hello Francois, Dženan,


Francois Budin-3 wrote
> I think you are missing a step in your computation:
> 1) You need to trnasform your point before rotation to find its position
> after transformation. You can directly use the affine transform for that:
> affineTransform->TransformPoint(my_point)
> 2) Compute the index of the transformed point:
>  filter->GetOutput()->TransformPhysicalPointToIndex(Pp, Pp_pix);
> 
> I think you forgot to compute 1)

I did not use the TransformPoint function but I did transform the points
using my own matrix multiplication method before using
transformPhysicalPointToIndex. And it gives the same results.



While cleaning up the code I discovered that the input 3D point I was
getting were actually not in proper physical space, they did not take image
origin into account, only pixel position * spacing. I'll try to see the
implication of that in my code and I'll come back to you.

Cheers,
A.




--
Sent from: http://itk-users.7.n7.nabble.com/


More information about the Insight-users mailing list