数字图像处理报告

退化结果');Рnoise=0.01*randn(size(I));РBn=imadd(B,im2uint8(noise));Рsubplot(2,3,3),imshow(Bn,[]);Рtitle('运动且加高斯噪声的图像');Рimg1=deconvwnr(B,PSF);Рsubplot(2,3,4),imshow(img1,[]);Рtitle('对运动模糊图像逆滤波');Р%对高斯噪声图像进行逆滤波复原Рimg2=deconvwnr(Bn,PSF);Рsubplot(2,3,5),imshow(img2,[]);Рtitle('对运动且加噪图像逆滤波');Р程序结果图:Р程序清单:Рclc;РI=imread('cameraman.tif');Р%转换图像矩阵为双精度型РI=im2double(I);РT=dctmtx(8);%产生二维DCT变换矩阵;Р%计算二维DCT,矩阵T及其转置是DCT函数P1*X*P2的参数РB=blkproc(I,[8,8],'P1*x*P2',T,T');Р%二值掩膜,只保留DCT系数中左上角10个Рmask=[1 1 1 1 0 0 0 0Р 1 1 1 0 0 0 0 0Р 1 1 0 0 0 0 0 0Р 1 0 0 0 0 0 0 0Р 0 0 0 0 0 0 0 0Р 0 0 0 0 0 0 0 0Р 0 0 0 0 0 0 0 0Р 0 0 0 0 0 0 0 0];РB2=blkproc(B,[8,8],'P1.*x',mask);%只保留DCT系数中左上角10个РI2=blkproc(B2,[8,8],'P1*x*P2',T,T');%重构图像Рsubplot(1,2,1),imshow(I);Рtitle('cameraman原图');Рsubplot(1,2,2),imshow(I2);Рtitle('压缩图像');