对二维数据进行fft2和fftshift后得到了频率波数图，但怎么确定横纵坐标的f，k值呢？

问题相关代码

全代码如下
filename = 'test1_inter_0.3-0.dat';
fid = fopen(filename, 'r');
int3 = fread (fid, [800 300],'float');
fclose(fid);
fftint3 = fft2(int3,800,300);
figure(4);
imagesc(abs(fftshift(fftint3)));