文档章节

图像二值化(直方图法,C语言实现)

t
 tigerBin
发布于 2017/07/25 15:25
字数 395
阅读 330
收藏 0

1.首先构造像素的频率直方图

2. 寻找直方图中的两个最大的波峰

3. 寻找这两个波峰之间的最小的波谷

4.波谷的index(像素值)为 K 值(阈值)


#include "JpegDecoder.h"
#include <stdio.h>
#include <opencv2/highgui.hpp>
#include <math.h>

using namespace JpegCodec;

static cv::Mat ConvertToMat(Matrix &mat)
{
	int channel = CV_8UC3;
	if (mat.channal == 1) channel = CV_8UC1;

	cv::Mat img(mat.rows, mat.cols, channel);  // create a new matrix

	for (int i = 0; i < mat.rows * mat.cols * mat.channal; i++)
	{
		img.data[i] = mat.data[i];
	}

	return img;
}


typedef struct
{
	int index;
	int freq; // frequence
}Point;

void Binary(Matrix &mat)
{
	Point peaks[256]; // 波峰
	int histogram[256];     //直方图

	// 构造直方图
	for (int i = 0; i < 255; i++) histogram[i] = 0;
	int cnt = mat.rows * mat.cols * mat.channal;
	for (int i = 0; i < cnt; i++)
	{
		histogram[mat.data[i]] ++;
	}

	// 统计波峰
	int size = 0;
	for (int i = 1; i < 254; i++)
	{
		if (histogram[i] > histogram[i + 1] && histogram[i] < histogram[i - 1])
		{
			peaks[size].index = i;
			peaks[size].freq = histogram[i];
			size++;
		}
	}

	// 寻找最大的两个波峰
	Point max = peaks[0], snd = peaks[0];
	for (int i = 0; i < size; i++)
	{
		if (peaks[i].freq > max.freq)
		{
			snd = max;
			max = peaks[i];
		}
		else if (peaks[i].freq > snd.freq)
		{
			snd = peaks[i];
		}
	}

	// 寻找波谷
	int i = (max.index > snd.index ? snd.index : max.index);
	size = (max.index > snd.index ? max.index : snd.index);
	Point K = max;
	for (; i < size; i++)
	{
		if (histogram[i] < K.freq)
		{
			K.freq = histogram[i];
			K.index = i;
		}
	}

	//二值化
	for (i = 0; i < cnt; i++)
	{
		if (mat.data[i] < K.index) mat.data[i] = 0;
		else mat.data[i] = 255;
	}
}

void ShowImage(Matrix &mat)
{
	cv::Mat img = ConvertToMat(mat);
	cv::imshow("Bitmap", img);
}


void Gray(Matrix &dst, Matrix &src)
{
	dst.Create(src.rows, src.cols, 1);

	for (int i = 0; i < src.rows; i++)
	{
		for (int j = 0; j < src.cols; j++)
		{
			int idx = (i * src.rows + j) * 3;
			dst.data[idx / 3] = (src.data[idx] + src.data[idx + 1] + src.data[idx + 2]) / 3;
		}
	}
}


int main(int argc, char *arrv[])
{
	JpegDecoder decoder("01.jpg");
	Matrix mat, dst;
	decoder.Decoder(mat);

	Gray(dst, mat);
	Binary(dst);
	ShowImage(dst);

	cvWaitKey(0);
	return 0;
}

JpegCodec : https://github.com/lzb-cc/JpegCodecs

运行示例

© 著作权归作者所有

共有 人打赏支持
t
粉丝 1
博文 65
码字总数 17440
作品 0
武汉
程序员
私信 提问
(三)matlab数字图像处理实验-图像灰度变换处理

前两节都是熟悉一下怎么在matlab底下对图片做一些操作,并没有什么卵用,这一节稍微有点卵用,灰度变换一般是图像处理的第一步。 数字图像处理实验1-9点击下列链接有源码和链接: matlab数字...

aninstein
01/06
0
0
论文学习:Contrast Limited Adaptive Histogram Equalization

目录 一、背景 1、对比度和直方图均衡HE 2、HE的问题 3、AHE 4、底噪问题 二、CLAHE 1、效果展示 2、算法格式和细节 一、背景 1、对比度和直方图均衡HE “对比度contrast ratio”这一概念,类...

RyanXing
07/28
0
0
图像的灰度图 和 二值化的概念

原文地址 在图像处理中,用RGB三个分量(R:Red,G:Green,B:Blue),即红、绿、蓝三原色来表示真彩色,R分量,G分量,B分量的取值范围均为0~255,比如电脑屏幕上的一个红色的像素点的三个...

liuhongyi0104
03/23
0
0
灰度图像的阈值化

通过对灰度图像二值化处理,能够凸现出感兴趣目标的轮廓. 灰度图像的二值化处理就是讲图像上的点的灰度置为0或255,也就是讲整个图像呈现出明显的黑白效果。即将256个亮度等级的灰度图像通过适...

文艺小青年
2017/07/06
0
0
图像直方图与直方图均衡化

图像直方图(英语:Image Histogram)是用以表示数字图像中亮度分布的直方图,标绘了图像中每个亮度值的像素数。可以借助观察该直方图了解需要如何调整亮度分布。这种直方图中,横坐标的左侧...

Tony沈哲
2017/05/22
0
0

没有更多内容

加载失败,请刷新页面

加载更多

java框架学习日志-2

上篇文章(java框架学习日志-1)虽然跟着写了例子,也理解为什么这么写,但是有个疑问,为什么叫控制反转?控制的是什么?反转又是什么? 控制其实就是控制对象的创建。 反转与正转对应,正转...

白话
32分钟前
0
0
Integer使用双等号比较会发生什么

话不多说,根据以下程序运行,打印的结果为什么不同? Integer a = 100;Integer b = 100;System.out.println(a == b);//print : trueInteger a = 200;Integer b = 200;System.out.pr...

兜兜毛毛
昨天
8
0
CockroachDB

百度云上的CockroachDB 云数据库 帮助文档 > 产品文档 > CockroachDB 云数据库 > 产品描述 开源NewSQL – CockroachDB在百度内部的应用与实践 嘉宾演讲视频及PPT回顾:http://suo.im/5bnORh ...

miaojiangmin
昨天
5
0
I2C EEPROM驱动实例分析

上篇分析了Linux Kernel中的I2C驱动框架,本篇举一个具体的I2C设备驱动(eeprom)来对I2C设备驱动有个实际的认识。 s3c24xx系列集成了一个基于I2C的eeprom设备at24cxx系列。at24cxx系列芯片包...

yepanl
昨天
6
0
设计模式之工厂模式

本篇博文主要翻译这篇文章: https://www.journaldev.com/1392/factory-design-pattern-in-java 由于翻译水平有限,自认为许多地方翻译不恰当,欢迎各位给出宝贵的建议,建议大家去阅读原文。...

firepation
昨天
15
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部