matlab最小分类错误全局二值化算法
转自：http://download.csdn.net/detail/hupeng810/1511870
function imagBW = kittlerMet(imag)
% KITTLERMET binarizes a gray scale image 'imag' into a binary image
% Input:
% imag: the gray scale image, with black foreground(0), and white
% background(255).
% Output:
% imagBW: the binary image of the gray scale image 'imag', with kittler's
% minimum error thresholding algorithm.
% Reference:
% J. Kittler and J. Illingworth. Minimum Error Thresholding. Pattern
% Recognition. 1986. 19(1):41-47
MAXD = 100000;
imag = imag(:,:,1);
[counts, x] = imhist(imag); % counts are the histogram. x is the intensity level.
GradeI = length(x); % the resolusion of the intensity. i.e. 256 for uint8.
J_t = zeros(GradeI, 1); % criterion function
prob = counts ./ sum(counts); % Probability distribution
meanT = x' * prob; % Total mean level of the picture
% Initialization
w0 = prob(1); % Probability of the first class
miuK = 0; % First-order cumulative moments of the histogram up to the kth level.
J_t(1) = MAXD;
n = GradeI-1;
for i = 1 : n
w0 = w0 + prob(i+1);
miuK = miuK + i * prob(i+1); % first-order cumulative moment
if (w0 < eps) || (w0 > 1-eps)
J_t(i+1) = MAXD; % T = i
else
miu1 = miuK / w0;
miu2 = (meanT-miuK) / (1-w0);
var1 = (((0 : i)'-miu1).^2)' * prob(1 : i+1);
var1 = var1 / w0; % variance
var2 = (((i+1 : n)'-miu2).^2)' * prob(i+2 : n+1);
var2 = var2 / (1-w0);
if var1 > eps && var2 > eps % in case of var1=0 or var2 =0
J_t(i+1) = 1+w0 * log(var1)+(1-w0) * log(var2)-2*w0*log(w0)-2*(1-w0)*log(1-w0);
else
J_t(i+1) = MAXD;
end
end
end
minJ = min(J_t);
index = find(J_t == minJ);
th = mean(index);
th = (th-1)/n
imagBW = im2bw(imag, th);
% figure, imshow(imagBW), title('kittler binary');
相关文章：
XShell连接Deepin
为什么80%的码农都做不了架构师？>>> 先让deepin安装openssh sudo apt-get install openssh-serverchkconfig ssh on 转载于:https://my.oschina.net/enzo/blog/110518
第三届“达观杯”文本智能信息抽取挑战赛丰厚奖金，群英集结，等你来战！...
近日，第三届“达观杯”文本智能信息抽取挑战赛正式上线启动（点击阅读原文，跳转报名页面），6月28日至8月31日，面向所有参赛选手开放竞赛结果提交。本届“达观杯”的任务是信息抽取。“达观杯”大赛由国内文本…
Spline interpolation and Savitzki-Golay smoothing
转自：http://octave.1599824.n4.nabble.com/Spline-interpolation-and-Savitzki-Golay-smoothing-td1675136.html ## natural-cubic-spline interpolation ## usage: yspline spline(x,y,xspline) ## example: ## x 0:10; y sin(x); ## xspline 0:0.1:10; y…
SpringBoot实现热部署(修改class不需要重启)
热部署: devtools可以实现页面热部署(即页面修改后会立即生效， 这个可以直接在application.properties文件中配置spring.thymeleaf.cachefalse来实现) 实现类文件热部署(类文件修改后不会立即生效),实现对属性文件的热部署。 注意：因为采用的虚拟机机制&…
Oracle中查看表空间的使用率的脚本
如题： select f.tablespace_name tablespace_name, round((d.sumbytes / 1024 / 1024 / 1024), 2) total_g, round(f.sumbytes / 1024 / 1024 / 1024, 2) free_g, round((d.sumbytes - f.sumbytes) / 1024 / 1024 / 1024, 2) used_g, round((d.sumbytes - f.sumbyte…
vue实现多个元素或多个组件之间动画效果
2019独角兽企业重金招聘Python工程师标准>>> 多个元素的过渡 <style>.v-enter,.v-leave-to{opacity: 0;}.v-enter-acitve,.v-leave-active{opacity: opacity 1s;} </style> <div idapp><transition><div v-ifshow>hello world</di…
干货 | 20个教程，掌握时间序列的特征分析（附代码）
作者 | Selva Prabhakaran 译者 | Tianyu责编 | Jane出品 | AI科技大本营（ID: rgznai100）【导语】时间序列是指以固定时间为间隔的序列值。本篇教程将教大家用 Python 对时间序列进行特征分析。1、什么是时间序列？时间序列是指以固定时间为间…
利用OpenCV求取图像的重心
转自：http://blog.csdn.net/lxiaoxiaot/article/details/6539834 不规则区域的矩，表示把一个归一化的灰度级图像函数理解为一个二维随机变量的概率密度。 这个随机变量的属性可以用统计特征--矩（Moments）来描述。通过假设非零的像…
C++中Ansi、Unicode、UTF8字符串之间的转换和写入
转自: http://dark0729.blogbus.com/logs/51496111.html Ansi字符串我们最熟悉，英文占一个字节，汉字2个字节，以一个\0结尾，常用于txt文本文件 Unicode字符串，每个字符(汉字、英文字母)都占2个字节，以2个连…
MyBatis的扩展点(plugins)
2019独角兽企业重金招聘Python工程师标准>>> 1、mybatis扩展点plugins mybatis的扩展是通过拦截器Interceptor来实现的，本质上就是JDK的动态代理，所以它只能对接口进行拦截，mybatis可以对以下四个接口类型进行拦截，也就…
linux中使用CST时间
GMT(Greenwich Mean Time，格林威治标准时间): 是指位于英国伦敦郊区的格林尼治天文台的标准时间，因为本初子午线被定义在通过那里的经线。 UTC(Universal Time/Temps Cordonn 世界标准时间)CST(Central Standard Time 國家標準時間，一說中原標…
到「黄埔学院」去：打造AI首席架构师，第二期限量招募！
今年 1 月，百度联合“深度学习技术及应用国家工程实验室”成立黄埔学院，旨在为产业培养第一批“首席AI架构师”。黄埔学院一期学员历时半年的学习和交流，6 月 16 日，黄埔学院一期学员迎来了毕业典礼，并在 7 月百度 AI开…
linux守护进程的创建
下面的完成了这样一个功能，创建一个守护进程，每个一秒在/tmp目录下的文件peng.txt中记录当前系统时间。 一、守护进程 守护进程是linux中的后台服务进程，在系统启动时开始运行，在系统关闭时终止。Linux系统中的大多数服务进程都是…
tesseract3.01的训练和使用
相关源码、资源下载：http://code.google.com/p/tesseract-ocr/downloads/list 训练步骤： 1、 Generate Training Images：生成tif图像文件(简单的几个汉字)； 如：ABC.Roman.exp0.tif([lang].[fontname].exp[num].tif)…
旷视推出鼻纹识别，用AI寻找丢失宠物
来源 | 转载自旷视城市大脑（ID：MEGVII_CityBrain)导读：随着人工智能技术（AI）的大热，基于深度学习方法的人脸识别技术已成熟落地，在解锁、支付、认证、摄像等生活方方面面，各个大厂推…
Qt浅谈之一：内存泄露（总结）
一、简介 Qt内存管理机制：Qt 在内部能够维护对象的层次结构。对于可视元素，这种层次结构就是子组件与父组件的关系；对于非可视元素，则是一个对象与另一个对象的从属关系。在 Qt 中，在 Qt 中，删除父对…
LINUX新手入门-1.装系统
LINUX新手入门-1.装系统首先我们用虚拟机模拟 装linux系统，然后下一步下一步，然后完成后，编辑一些设置，把镜像放上面就可以了选第一项，安装系统，查看镜像是否能运行，直接跳过，选择语…
Log4cplus1.04的使用
首先，从http://sourceforge.net/projects/log4cplus/files/log4cplus-stable/下载最新的版本，解压缩，用vs2008打开msvc8文件夹下的log4cplus.sln，并按照提示转换。在Solution Configurations下拉列表框中，会有Debug、D…
FRVT赛程全纪录：格灵深瞳全球排名前五
作者 | 张德兵，格灵深瞳首席科学家&算法部负责人来源 | 转载自知乎张德兵最近两个月，格灵深瞳首席科学家&算法部负责人张德兵与算法团队参加了全球人脸识别算法测试(FRVT、Face Recognition Vendor Test)。虽然是第一次参加此比赛，格…
反转比特位(文章最后有干货)【转】
转自：https://blog.csdn.net/wuxianglonghaohao/article/details/21602305 http://www.newhottopic.com/2014/03/20/reverse-bits/ 把一个无符号整数的比特位反转顺序。有很多种方法来实现这个。我们这里给出一个算法：通过异或运算来交换，然后…
过关斩将打进Kaggle竞赛Top 0.3%，我是这样做的
作者 | Lavanya Shukla译者 | Monanfei责编 | 夕颜出品 | AI科技大本营（id：rgznai100）导读：刚开始接触数据竞赛时，我们可能会被一些高大上的技术吓到。各界大佬云集，各种技术令人眼花缭乱，新手们…
JavaBean规范
2019独角兽企业重金招聘Python工程师标准>>> （1）JavaBean 类必须是一个公共类，并将其访问属性设置为 public （2）JavaBean 类必须有一个空的构造函数：类中必须有一个不带参数的公用构造器&#x…
vigra1.8.0的使用
VIGRA stands for "Vision with Generic Algorithms". Its a novel computer vision library that puts its main emphasis oncustomizablealgorithms and data structures. 1、首先，从http://hci.iwr.uni-heidelberg.de/vigra/下载最新源代码࿰…
17个Python小窍门
python中相对不常见却很实用的小窍门。 空谈不如来码代码吧： 交换变量值 给列表元素创建新的分隔符 找列表中出现次数最多的元素 核对两个字符是否为回文 反向输出字符串 反向输出列表 转置2维数组 链式比较 我刚整理了一套2018最新的0基础入门和进阶教程࿰…
用产品思路建设中台，这走得通吗？| 白话中台
作者 | 王健，ThoughtWorks首席咨询师。 十多年国内外大型企业软件设计开发，团队组织转型经验。一直保持着对技术的热爱，热衷于技术分享。目前专注在企业平台化转型、中台战略规划，微服务架构与实施，大型遗留系统服务化…
利用cvMinAreaRect2求取轮廓最小外接矩形
转自：http://blog.csdn.net/mine1024/article/details/6044856 对给定的 2D 点集，寻找最小面积的包围矩形，使用函数： CvBox2D cvMinAreaRect2( const CvArr* points, CvMemStorage* storageNULL ); points 点序列或点集数组 …
电脑开机显示Invalidsystemdisk
开机或重启无法进入系统，并在屏幕上显示Invalidsystemdisk，Replacethediskandthenpressanykey或者diskerror之类的字样，这是怎么回事，该如何解决？今天u大师就为大家解决下。 出现这个原因是因为现在的电脑没有可以启…
Windows7 64位下vs2008配置OpenCV2.3.1
1、下载OpenCV2.3.1：http://www.opencv.org.cn/index.php/Download； 2、下载后解压缩：OpenCV-2.3.1-win-superpack.exe，生成一个opencv文件夹； 3、下载CMake：http://www.cmake.org/cmake/resources/softw…
腾讯拥抱开源：首次公布开源路线图，技术研发向共享、复用和开源迈进
整理 | 夕颜出品 | AI科技大本营（ID:rgznai100）导读：去年，知乎上一篇讨论腾讯技术的帖子异常火爆，讨论的主题是当下（2018 年）腾讯的技术建设是否处于落后同体量公司的状态，这篇帖子得…
Babylon.js 3.3发布：更强大的粒子系统和WebVR支持
Babylon.js 3.3版本利用微软混合现实工具包（MRTK）的功能来改进WebVR开发，并改进了其粒子系统控件。 MRTK提供了一系列脚本和组件来加速混合现实应用程序的开发。为了简化GUI VR构建，Bablyon.js利用3D体积网格来布局VR场景的界面&a…