实验二信息熵与图像熵.docVIP

信电学院 信息论与编码 实验报告书 （2013/2014学年第二学期） 实验名称 ： 信息熵与图像熵计算 专业班级 ： 通信三班 学生姓名 ： 包博文 学 号： 120310323 指导教师 ： 张龙 设计成绩 ： 2014年 5月10日 实验二 信息熵与图像熵计算 一、实验目的 1 通过本实验复习MATLAB的基本命令，熟悉MATLAB下的基本函数 2 复习信息熵基本定义，能够自学图像熵定义和基本概念 二、实验要求 1 能够写出MATLAB源代码，求信源的信息熵 2 根据图像熵基本知识，设计出MATLAB程序，求出给定图像的图像熵 三、实验步骤 1 求解信息熵过程： ① 输入一个离散信源，并检查该信源是否是完备集。 ② 去除信源中符号分布概率为零的元素。 ③ 根据平均信息量公式，求出离散信源的熵。 2 图像熵计算过程： ① 输入一幅图像，并将其转换成灰度图像。 ② 统计出图像中每个灰度阶象素概率。 ③ 统计出图像中相邻两象素的灰度阶联合分布矩阵。 ④ 根据图像熵和二阶熵公式，计算出一幅图像的熵。 四、程序的流程图 1 信息熵 2 图像熵 五、试验运行结果记录 1、信息熵 （1）、A=[0.3 0.5 0.2 0.1]; H=xinxishang(A) ???Error using==xinxishang 不是完备集。 （2）、A=[0.3 0.05 0.35 0.1 0.2 0]; H=xinxishang(A) 去除概率为零的元素 A =0.3000 0.0500 0.3500 0.1000 0.2000 H =2.0639； （3）、A=[0.05 0.1 0.3 0.1 0.35]; H=xinxishang(A) H=2.0639; 2、图像熵 [H1,H2] = tuxiangshang(111.jpg) H1 =7.2250; H2 =5.3577; 六、编写的程序 1、信息熵 function H=xinxishang(A) if sum(A)~=1 error(不是完备集) end if find(A==0) A(find(A==0))=[]; disp(去除概率为零的元素);A end H=-sum(A.*log2(A)); 2、图像熵 function [H1,H2]=tuxiangshang(img) I=imread(img); img=rgb2gray(I); imview(I), imview(img); [ix,iy]=size(img); P1=imhist(img)/(ix*iy); temp=double(img); temp=[temp,temp(:,1)]; A=zeros(256,256); for x=1:ix for y=1:iy i=temp(x,y);j=temp(x,y+1); A(i+1,j+1)=A(i+1,j+1)+1; end end P2=A./(ix*iy);H1=0;H2=0; for i=1:256 if P1(i)~=0 H1=H1-P1(i)*log2(P1(i)); end for j=1:256 if P2(i,j)~=0 H2=H2-P2(i,j)*log2(P2(i,j)); end end end H2=H2/2 七、实验小结 通过本次实验熟悉了Matlab软件编写程序环境和一些函数的功能及使用，掌握了信息熵，图像熵的计算方法。 开始 判断是否 为完备集 查找去除为 零的元素 计算信息熵 结束 输出错误 否 是 开始 读取图像 转化为灰度图像 计算图像大小并获取数据直方图 判断是否循环 统计灰度阶联合分布矩阵 判断是否循环 计算图像熵 结束 假