博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
matlab练习程序(Renyi熵)
阅读量:6036 次
发布时间:2019-06-20

本文共 658 字,大约阅读时间需要 2 分钟。

Renyi熵是对通常的香农熵的扩展,算是q阶广义熵。公式如下:

其中P和香农熵公式中的P一样,是概率。当q=1时公式退化为香农熵公式。(如何证明?)

有用此公式寻找图像最佳二值化阈值的。

首先定义前景区域A,背景区域B。

那么前景与背景区域像素相应的Renyi熵就如下定义:

其中k是当前取的灰度级,P(A)是像素在A区域的概率,P(B)类似。当然,这里说区域,不是指空间区域,是像素灰度级区域。

最后图像Renyi熵求最佳阈值定义如下:

这里得到的K就是分割阈值。

分割效果如下:

原图:

分割后:

在编程时还需要确定阶数q,我取的是2。

matlab代码如下:

clear all;close all;clc;img=imread('lena.jpg');[m n]=size(img);imshow(img)Hist=imhist(img);q=2;H=[];for k=2:256    PA=sum(Hist(1:k-1));    PB=sum(Hist(k:255));        Pa=Hist(1:k-1)/PA;    Pb=Hist(k:256)/PB;        HA=(1/1-q)*log(sum(Pa.^q));    HB=(1/1-q)*log(sum(Pb.^q));        H=[H HA+HB];    end[junk level]=max(H);imgn=im2bw(mat2gray(img),level/256);figure;imshow(imgn)

 

转载地址:http://chohx.baihongyu.com/

你可能感兴趣的文章
Spring Security4实战与原理分析视频课程( 扩展+自定义)
查看>>
第一周博客作业
查看>>
thinkpython2
查看>>
String、StringBuffer和StringBuilder的区别
查看>>
oracle recyclebin与flashback drop
查看>>
svmlight使用说明
查看>>
Swing 和AWT之间的关系
查看>>
Mysql设置自增长主键的初始值
查看>>
Android计时器正确应用方式解析
查看>>
获取post传输参数
查看>>
ASP生成静态页面的方法
查看>>
HDU 1325 Is It A Tree? 判断是否为一棵树
查看>>
Shell命令-文件压缩解压缩之gzip、zip
查看>>
个人总结
查看>>
uva 673 Parentheses Balance
查看>>
Bzoj 2252: [2010Beijing wc]矩阵距离 广搜
查看>>
css 禁止选中文本
查看>>
bzoj2165
查看>>
tomcat 配置首页
查看>>
算术运算表达式正则及分析
查看>>