In case it can help someone later, I was able to recover proper orientation of the volume saved by ITK with the following 2 operators:
scanArray = np.swapaxes(scanArray,0,2).copy()
scanArray = np.rollaxis(scanArray,2,1).copy()
3 Likes