data = load_nii(‘ImageName.nii’);
fftImg = fftn(data.img);
realPart = make_nii(real(fftmg),[1 1 1]);
imagPart = make_nii(imag(fftImg),[1 1 1]);
save_nii(realPart,‘fftRealPart.nii’);
save_nii(imagPart,‘fftImagPart.nii’);
flip_lrtd('fftRealPart.nii",“fftRealPart_flipped_lrtd”);
flip_lrtd('fftImagPart.nii",“fftImagPart_flipped_lrtd”);
Here are the infos on the toolbox.
Last information that might be useful, ITK transformation matrix (displayed by slicer in the volume module) is [-1 0 0; 0 -1 0; 0 0 1]. The matlab transformation matrix generated before the flip is [1 0 0; 0 1 0; 0 0 1].