olisp:

but I should not permute the axes as the data uses entirely different coordinate system ?


Could you explain this statement a bit?
Does the depth dimension represent the number of slices or what do you mean by “different coordinate system”?
What wouldn’t it be possible to permute the data?