Get the world coordinates from VTK, then pass them to itk::Image
's method TransformPhysicalPointToIndex()
. That method will compute an itk::Index
of the voxel in question.
Significant difference between itk::Image
and vtkImageData
: VTK’s image does not have orientation, whereas ITK’s does.
A bit deeper discussion is in [SOLVED]TransformIndexToPhysicalPoint manually.