I did a registration in two steps, Affine -> B-Spline, and now I want to know what is the correspondent index in the new space for a given index in the original space. Here is the code I’m using:
#Receiving the voxel coordinates print("Please, insert the voxel coordinates (X, Y, Z): ") voxel = [int(x) for x in input().split()] physical_voxel = moving_volume.TransformIndexToPhysicalPoint(voxel) print(physical_voxel) #Applying the transformation transformed_point = BS_transformation.TransformPoint(Affine_transformation.TransformPoint(physical_voxel)) print(transformed_point) fixed_index = Fixed_volume.TransformPhysicalPointToIndex(transformed_point) print(fixed_index) label = fixed_atlas.GetPixel(fixed_index)
Voxel is the index in the original space.
The problem is: I do not get the correct index in the registered space but for me everything is right. So, How to solve this?