I am trying to further understand how itk handles sort.
It seems like the foundation is calculating a normal vector in:
cosines = ImageHelper::GetDirectionCosinesValue( **it ); // You only have to do this once for all slices in the volume. Next, // for each slice, calculate the distance along the slice normal // using the IPP ("Image Position Patient") tag. // ("dist" is initialized to zero before reading the first slice) : normal = cosines*cosines - cosines*cosines; normal = cosines*cosines - cosines*cosines; normal = cosines*cosines - cosi*emphasized text*nes*cosines;
Is there documentation in dicom docs or somewhere else on how this is working? I have seen many people ask a similar question on sorting so I would love to understand this solution as it seems robust against all orientations including obliques.
Thanks for any insights!