%% 学习目标：离散余弦变换图像去噪
X=imread('01.png'); %读取图像
[m,n]=size(X); %读取图像尺寸
Xnoised=imnoise(X,'speckle',0.01); %给图像加噪
subplot(1,2,1),imshow(Xnoised),title('加噪后图像'); %输出加噪图像
Y=dct2(Xnoised); %DCT变换
I=zeros(m,n);
I(1:m/3,1:n/3)=1; %高频屏蔽
Ydct=Y.*I;
Y=uint8(idct2(Ydct)); %逆DCT变换
subplot(1,2,2),imshow(Y),title('DCT去噪后图像'); %结果输出
%% 大仙QQ：1960009019
%% 在线教育微信公众号：大仙一品堂
