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);
}
查看全文
相关阅读:
Apache httpd和JBoss构建高可用集群环境
Ubuntu 14.04下NFS安装配置
Ubuntu 14.04 安装 JDK 7.0
Docker第三方项目小结
Shipyard远程API
Linux软件管理——yum命令详解
Quartz集群原理及配置应用
Rsync原理介绍及配置应用
python构造wireshark可以解析的LTE空口数据
A Simple Web Server
原文地址:https://www.cnblogs.com/wqj1212/p/1009655.html
最新文章
HttpClient方式调用接口的实例
JAVA中创建线程的三种方法及比较
JAVA中发送电子邮件的方法
alibaba架包FastJson使用例子
MyEclipse中Lombok的安装及使用
消息发送函数OSMboxPostOpt()
NB-IOT技术学习问题记录
雨量传感器测试(大众凌渡高尔夫7雨量传感器)
STM32F405的内部ADC采集
海思NB-IOT模组在平台上注册
热门文章
NB-IOT技术 UP模式 和CP模式,用户面和控制面,数据面
英雄联盟测试静态IP(固态IP)和动态IP的网速测试
CRC校验的问题(c语言实现CRC校验和)
Centos7下卸载docker
Centos7下使用ELK(Elasticsearch + Logstash + Kibana)搭建日志集中分析平台
Linux设置FQDN
ubunt1204安装配置vsftp
阅读推荐——深入浅出Mesos
ffmpeg无法接收组播流问题处理
Centos7 安装配置NFS
Copyright © 2011-2022 走看看