zoukankan
html css js c++ java
OTSU求阈值程序
writen by
wqj1212@yahoo.com.cn
int
otsu (IplImage
*
image,
int
rows,
int
cols,
int
x0,
int
y0,
int
dx,
int
dy,
int
vvv)
{
unsigned
char
*
np;
//
图像指针
int
thresholdValue
=
1
;
//
阈值
int
ihist[
256
];
//
图像直方图,256个点
int
i, j, k;
//
various counters
int
n, n1, n2, gmin, gmax;
double
m1, m2, sum, csum, fmax, sb;
//
对直方图置零
memset(ihist,
0
,
sizeof
(ihist));
gmin
=
255
; gmax
=
0
;
//
生成直方图
/**/
/*
for (i = y0 + 1; i < y0 + dy - 1; i++) {
np = &image[i*cols+x0+1];
for (j = x0 + 1; j < x0 + dx - 1; j++) {
ihist[*np]++;
if(*np > gmax) gmax=*np;
if(*np < gmin) gmin=*np;
np++; /* next pixel
}
}
*/
for
(j
=
y0;j
<
dy;j
++
)
{
for
(i
=
0
;i
<
dx;i
++
)
{
unsigned
char
temp
=
CV_IMAGE_ELEM(image,uchar,j,i);
ihist[temp]
++
;
}
}
//
set up everything
sum
=
csum
=
0.0
;
n
=
0
;
for
(k
=
0
; k
<=
255
; k
++
)
{
sum
+=
(
double
) k
*
(
double
) ihist[k];
/**/
/*
x*f(x) 质量矩
*/
n
+=
ihist[k];
/**/
/*
f(x) 质量
*/
}
if
(
!
n)
{
//
if n has no value, there is problems
fprintf (stderr,
"
NOT NORMAL thresholdValue = 160\n
"
);
return
(
160
);
}
//
do the otsu global thresholding method
fmax
=
-
1.0
;
n1
=
0
;
for
(k
=
0
; k
<
255
; k
++
)
{
n1
+=
ihist[k];
if
(
!
n1)
{
continue
; }
n2
=
n
-
n1;
if
(n2
==
0
)
{
break
; }
csum
+=
(
double
) k
*
ihist[k];
m1
=
csum
/
n1;
m2
=
(sum
-
csum)
/
n2;
sb
=
(
double
) n1
*
(
double
) n2
*
(m1
-
m2)
*
(m1
-
m2);
/**/
/*
bbg: note: can be optimized.
*/
if
(sb
>
fmax)
{
fmax
=
sb;
thresholdValue
=
k;
}
}
//
at this point we have our thresholding value
//
debug code to display thresholding values
if
( vvv
&
1
)
fprintf(stderr,
"
# OTSU: thresholdValue = %d gmin=%d gmax=%d\n
"
,thresholdValue, gmin, gmax);
return
(thresholdValue);
}
查看全文
相关阅读:
在 Docker 搭建 Maven 私有库
Maven:mirror和repository 区别
ubuntu DEBIAN_FRONTEND环境变量用法
Redis常见面试题
Error:(1, 1) java: 非法字符: ‘ufeff’
jpa 查询方法和sql查询语句对应关系
net.sf.json.JSONObject对时间戳的格式化处理
美团Leaf——全局序列生成器
Logstash
Kafka和SpringBoot
原文地址:https://www.cnblogs.com/wqj1212/p/1009655.html
最新文章
python redis 的基本操作指令
PGSQL常用操作
图解 | 原来这就是网络
Python Flask web入门例子
基于TypeScript装饰器定义Express RESTful 服务
TypeScript装饰器(decorators)
TypeScript和Node模块解析策略
webpack打包非模块化js
TypeScript笔记 6--接口
TypeScript笔记 5--变量声明(解构和展开)
热门文章
TypeScript笔记 4--变量声明
TypeScript笔记 3--基础类型
TypeScript笔记 2--代码调试
Ubuntu Server 19.04配置静态IP
Nginx+uWSGI+Django+Python在Linux上的部署
MySQL 数据备份与还原
阿里云的maven仓库
Linux下is not in the sudoers file解决方法
修改 Ubuntu SSH 登录后的欢迎信息
ubuntu14.04设置静态ip
Copyright © 2011-2022 走看看