哈弗曼编码 图像处理 压缩率

%以下为主程序mainp.m clc clear close all; %定义HufData/Len为全局变量的结构体 global HufData; global Len disp('计算机正在准备输出哈夫曼编码结果,请耐心等待……'); %原始码字的灰度 a=imread('Fig3.41(c).jpg'); %分区画出原始图像和灰度直方图 figure; subplot(1,2,1) imshow(a); %取消坐标轴和边框 axis off box off title('系统图像','fontsize',13); subplot(1,2,2); axis off box off imhist(a); title('图像灰度直方图','fontsize',13); %图像的灰度统计 GrayStatistics=imhist(a); GrayStatistics=GrayStatistics'; GrayRatioo=GrayStatistics/sum(GrayStatistics); GrayRatioNO=find(GrayRatioo~=0); Len=length(GrayRatioNO); %初始化灰度集，防止系统随即赋予其垃圾值 GrayRatio=ones(1,Len); for i=1:Len GrayRatio(i)=GrayRatioo(i); end GrayRatio=abs(sort(-GrayRatio)); %将图像灰度概率赋予结构体 for i=1:Len HufData(i).value=GrayRatio(i); end % 哈夫曼编码/霍夫曼编码 HuffmanCode(Len); %输出码字 zippedHuffman=1; for i=1:Len tmpData=HufData(i).code; str=''; for j=1:length(tmpData) str=strcat(str,num2str(tmpData(j))); zippedHuffman=zippedHuffman+1; end disp(strcat('a',num2str(i),'= ',str)) end i; %计算计算机一共输出多少个哈夫曼编码/霍夫曼编码 zippedHuffman; %计算在删去0灰度级压缩之前的原始图像字节容量 unzipped_delete=i*8; %计算压缩比率 ratio_delete=zippedHuffman/unzipped_delete; %计算图像的压缩比率 ad=num2str(ratio_delete*100); str2=strcat(ad,'%'); disp(strcat('哈夫曼编码压缩比率','= ',str2))