弱弱的问，图像频谱，怎么变成功率谱？

tjucruiser

I_sub=imread(demo.bmp');%图像是2048*2048的

%在进行傅立叶变换前，需要设置什么参数吗？我记得一维的傅立叶变换是需要设置采样频率之类的

I_FFT=fftshift(fft2(I_sub));%这里得到的是图像的频谱，没问题吧？

PS=abs(IFFT)*abs(I_FFT);%这么做是得到图像功率谱了吗？这个I_FFT竟然也是2048*2048的（复数矩阵），难道这个结果（即2048*2048)不受图像像素点大小什么的影响吗

Two_D=log(1+NPS);%对数功率谱
imshow(Two_D,[]);
colormap(jet(256)),colorbar  
%下面我想得到一维的功率谱随着频率的变化，即在二维图的中心位置取一个剖面。
One_D=Two_D(1000:1030,1024:2048);
One_D=mean(One_D,1);