弱弱的问,图像频谱,怎么变成功率谱?
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);
页:
[1]