数字图像处理图像编码286.pdfVIP

数字图像处理上机实习报告 （DIP4DIP7 ） 学生姓名： 杜坤 班 级： 071123 学 号： 20121003699 指导老师： 傅 华 明 1 DIP-4 图像编码 一．题目要求 对图实施费诺-香农编码和解码，计算图像熵，平均码长和冗余度。 二．算法设计 1.测试脚本的程序框图 校对编码前后的数 开始 编码 计算图像的熵 据 读入图像的 根据编码的码字对 计算图像的 结束 数据为a 图像数据进行输出 平均码长 统计各个灰度 解码 编码的编码效率 值的概率 将解码后的数据 计算冗余度 data变行为8*8 将码字初始化 2.编码程序框图 读入图像的直方图，将图像的灰度值按照概率大小排序，按照香农编码的规 则编码。 香农编码将概率由大到小，由上到下排成一排，然后分为两组。是将大的一 组概率赋值为 0，概率小的一组赋值为 1，这是赋值的原则。然后依次的重复， 直到每组只有一种输入元素为止。 2 3.解码程序框图 三．实现代码 1.脚本文件 clear all load mat p = impr(a); %统计概率 code = FanoCodeInit(p); %Fano 编码初始化 3 code = FanoEncoder(code);%Fano 编码 outstream = FanoCodeStream(a,code); %输出 data = FanoDecoder(outstream,code);%解码 data = reshape(data,8,8); %恢复 8*8 的形状 data = data; %转置 I = abs(p.*log2(p)); disp(图像的熵为：); H = sum(I(:)) %计算熵 disp(图像的平局码长为：) B = FanoCodeLength(code); %求平均长度 disp(编码冗余度为： ); r = B/H - 1 %求冗余 disp(编码效率为： ) e = H/B %求编码效率 if isequal(a,data) msgbox( 解码后的数据和输入的数据完全吻合 ); end 2.统计灰度的概率 function [p]= impr(f) %概率统计 [m,n] = size(f); graymax = max(f(:)); %找出灰度最大值，划定统计范围 p = zeros(1,graymax + 1); for i = 1:m for j = 1:n x = f(i,j) + 1; p(x) = p(x) + 1;