code.txt

离散余弦新型编码算法

function [ output_args ] = Untitled1( input_args )
%UNTITLED1 Summary of this function goes here
%  Detailed explanation goes here
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%载入原始图像
I=imread('lena.bmp');
%装入原始图像
J=dct2(I) ;
%对I进行二维DCT
Q=log(abs(J));
imshow(Q);
%%%%上边两行作了改动，是添加上的%%
%%%%%%imshow(log(abs(J)),[]);
figure(1);
colormap(jet(64)),colorbar; 
figure(2);
%进行 DCT变 换后的二维变换谱
J(abs(J)<100)=0; 
% 其中‘10’ 可以改为其它值
%将DCT变换值小于10 的元素设为0
K=idct2(J)/255; 
%对逆DCT变换值归一化
imshow(I);
figure(3);
title('图像K ');
imshow(K);
figure(4);
%显示原始图像 I和压缩图像 K
error=J(:,2)-K(:,2);
%下面要画出误差函数
plot (error) ;
title('误差函数');
figure(5);
MSE=sum(error(:,1))/prod(size(K))
%计算归一化图像的均方误差值