图像重建MATLAB tell

qq_39215987 2019-04-28 07:24:16
clcclearI=imread('lena.bmp');I=rgb2gray(I);fftI=fft2(I);%二维离散傅立叶变换sfftI=fftshift(fftI);%直流分量移到频谱中心RR=real(sfftI);%取傅立叶变换的实部II=imag(sfftI); %取傅立叶变换的虚部A=sqrt(RR.^2+II.^2); %计算频谱幅值B=(A-min(min(A)))./(max(max(A))-min(min(A)))*225; %归一化fudupu=ifft2(abs(fftI)); %%幅度谱重建图像pinpu =ifft2(abs(fftI).*exp(j*(angle(fftI)))); %%频谱重建图像xiangweipu=ifft2(exp(j*(angle(fftI)))); %%相位谱重建figure%设定窗口subplot(2,3,1),imshow(I);title('原图');subplot(2,3,2),imshow(B);title('原图像的频谱'); %显示原图像的频谱subplot(2,3,3),imshow(log(B),[]);title('原图像取对数的频谱'); %显示原图像的频谱subplot(2,3,4),imshow(fudupu,[]);title('幅度谱重建图像');subplot(2,3,5),imshow(pinpu,[]);title('频谱重建图像'); %显示原图像的频谱subplot(2,3,6),imshow(xiangweipu,[]);title('相位谱重建'); %显示原图像的频谱
...全文
14 回复 打赏 收藏 转发到动态 举报
写回复
用AI写文章
回复
切换为时间正序
请发表友善的回复…
发表回复

27,579

社区成员

发帖
与我相关
我的任务
社区描述
MS-SQL Server 应用实例
社区管理员
  • 应用实例社区
加入社区
  • 近7日
  • 近30日
  • 至今
社区公告
暂无公告

试试用AI创作助手写篇文章吧