学号：111204510    姓名：邹龙   
      1、对一幅图像进行平移，显示原始图像与处理后图像，分别对其进行傅里叶变换，显示变换后结果，分析原图的傅里叶谱与平移后傅里叶频谱的对应关系。
s=imread('d:\1.jpg');%读入原图像
i=rgb2gray(s)
i=double(i)
j=fft2(i);                        %傅里叶变换   
k=fftshift(j);                  % 直流分量移到频谱中心
l=log(abs(k));                    %对数变换
m=fftshift(j);                    %直流分量移到频谱中心
RR=real(m);                        %取傅里叶变换的实部
II=imag(m);                        %取傅里叶变换的虚部
A=sqrt(RR.^2+II.^2);                %计算频谱府幅值
A=(A-min(min(A)))/(max(max(A)))*255;  % 归一化
b=circshift(s,[800 450]);  %对图像矩阵im中的数据进行移位操作matlab求傅里叶变换
b=rgb2gray(b)
b=double(b)
c=fft2(b);                        %傅里叶变换   
e=fftshift(c);                  % 直流分量移到频谱中心
l=log(abs(e));                    %对数变换
f=fftshift(c);                    %直流分量移到频谱中心
WW=real(f);                        %取傅里叶变换的实部B
ZZ=imag(f);                        %取傅里叶变换的虚部
B=sqrt(WW.^2+ZZ.^2);                %计算频谱府幅值
B=(B-min(min(B)))/(max(max(B)))*255;  % 归一化
subplot(2,2,1);imshow(s);title('原图像')
subplot(2,2,2);imshow(uint8(b));;title('平移图像')
subplot(2,2,3);imshow(A);title('离散傅里叶频谱');
subplot(2,2,4);imshow(B);title('平移图像离散傅里叶频谱')
2、对图像进行除法运算
moon=imread('moon.tif');
subplot(1,2,1);imshow(moon);
I=double(moon);
J=I*0.43+80; moon2=uint8(J);
Ip=imdivide(moon,moon2);
subplot(1,2,2);imshow(Ip,[]);

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。