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.