Thanks for providing more details to reproduce it.
Please post the minimal set of matlab commands you are using to read the image, apply the fftn, and write the output image, with the switch of axis you commented.
The output image generated by a fft has positive and negative frequencies, and also Nyquist frequencies depending on the oddness in size of the axis, so after switching of axis, the output has to be interpreted with care.