桂电图像处理实验5.docVIP

主程序 global HufData; global Len; img=imread(E:\gd.jpg); subplot(1,2,1);imshow(img);title(彩图); I=rgb2gray(img);subplot(1,2,2);imshow(I);imhist(I);title(直方图); GrayStatistics=imhist(I); %统计 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 %计算计算机一共输出多少个哈夫曼编码 zippedHuffman; %计算在删去0灰度级压缩之前的原始图像直接容量 unzipped_delete=i*8; %计算压缩比率 ratio_delete=zippedHuffman/unzipped_delete; %计算图像压缩比率 ad=num2str(ratio_delete*100); str2=strcat(ad,%); disp(strcat(哈夫曼编码压缩比率,=,str2)); %排序 function HuffmanCode(OriginSize) global HufData; global Len; global SortList; for i=1:OriginSize HufData(i).left=1;%%霍夫曼编码树左边纪录为1 HufData(i).right=0;%%霍夫曼编码树右边纪录为0 HufData(i).code=[];%%输出码初始化为0 SortList(i).symbol=i;%%排序列表初始化 SortList(i).value=HufData(i).value; end newsymbol=OriginSize;%初始化原始消息数目 for n=OriginSize:-1:2 SortList=sortdata(SortList,n);%将N个消息进行排序 newsymbol=newsymbol+1;%将最后两个出现概率最小的消息合成一个消息 HufData(newsymbol).value=SortList(n-1).value+SortList(n).value; HufData(newsymbol).left=SortList(n-1).symbol; HufData(newsymbol).right=SortList(n).symbol; SortList(n-1).symbol=newsymbol;%将消息添加到列队的最后，为N-1个消息重新排序作好准备 SortList(n-1).value=HufData(newsymbol).value; end visit(newsymbol,Len,[]);%遍历霍夫曼树，获得霍夫曼编码/哈夫曼编码 end 遍历霍夫曼编码 function visit(node,n,ocode) global HufData if node=n %如果没有哈夫曼编码/霍夫曼编码树的子接点直接输出原始码，这里为空码（[]） HufData(node).code=ocode; else if(HufData(node).left0) %遍历左分支接点输出1，这里采用子函数嵌套调用 ocode1=[ocode 1]; v