matlab图像傅里叶变换图像滤噪

电子1004黄佳卿201081163图像添加高斯噪声与去除代码如下：I=imread('C:\Users\jiaqing\Desktop\matlabhomework\zhen.jpg');%读入图像%给图像添加噪声grayI=rgb2gray(I);nI=imnoise(grayI,'gaussian',0.05);subplot(2,2,1);imshow(nI);title('加入高斯噪声后的图像');%加入噪声后图像的傅立叶变换FI=abs((fft2(nI)));NFI=255*mat2gray(FI);%归一化SFI=fftshift(NFI);subplot(2,2,2),imshow(SFI);title('加入噪声后的傅里叶变换');xd=wiener2(nI,[33]);subplot(2,2,3);imshow(xd);%去除噪声后图像的傅立叶变换FI=abs((fft2(xd)));NFI=255*mat2gray(FI);%归一化SFI=fftshift(NFI);subplot(2,2,4),imshow(SFI);title('去除噪声后的傅里叶变换');处理结果：滤波矩阵代码：I=imread('zhen.jpg');subplot(2,3,1);imshow(I);title('彩色原图显示');grayI=rgb2gray(I);subplot(2,3,2);imshow(grayI);title('RGB图像的显示');filter1=[100;010;000;];filter2=[000;111;111;];afI1=imfilter(grayI,filter1);subplot(2,3,3);imshow(afI1);title('矩阵1滤波后图像的显示');afI2=imfilter(grayI,filter2);subplot(2,3,4);imshow(afI2);title('矩阵2滤波后图像的显示');FI=abs((fft2(afI1)));NFI=255*mat2gray(FI);%归一化SFI=fftshift(NFI);subplot(2,3,5),imshow(SFI);title('矩阵1滤波后图像傅里叶变换');FI=abs((fft2(afI2)));NFI=255*mat2gray(FI);%归一化SFI=fftshift(NFI);subplot(2,3,6),imshow(SFI);title('矩阵2滤波后图像傅里叶变换');没太理解题意，不知道思路对不对，不调用系统函数，自己写代码没写成功。I=imread('zhen.jpg');subplot(2,3,1);imshow(I);title('彩色原图显示');grayI=rgb2gray(I);subplot(2,3,2);imshow(grayI);title('RGB图像的显示');[m,n]=size(I);filter=[11111;11111;11111;11111;11111;];bigI=blkdiag(grayI,zeros(2,1));%补零subplot(2,3,3);imshow(bigI);title('补零后图像的显示');manyI=mat2cell(bigI,ones(770/5,1)*5,ones(1025/5,1)*5);%分割??????RGB???????1??????????2??????????1????????????2??????????