文档章节

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

t
 tigerBin
发布于 2017/07/25 15:25
字数 395
阅读 223
收藏 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
粉丝 0
博文 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

没有更多内容

加载失败,请刷新页面

加载更多

区块链入门教程以太坊源码分析ethdb源码分析

兄弟连区块链入门教程以太坊源码分析ethdb源码分析,2018年下半年,区块链行业正逐渐褪去发展之初的浮躁、回归理性,表面上看相关人才需求与身价似乎正在回落。但事实上,正是初期泡沫的渐退...

兄弟连区块链入门教程
16分钟前
1
0
Eclipse下忽略掉node_modules目录相关配置

https://blog.csdn.net/yzf913214/article/details/72872523

张宏亮1982
17分钟前
1
0
在Python-dataframe中如何把出生日期转化为年龄?

我们在做数据挖掘项目或大数据竞赛时,如果个体是人的时候,获得的数据中可能有出生日期的Series,举个简单例子,比如这样的一些数: # -*- coding: utf-8 -*- import pandas as pd from ...

Mr_zebra
20分钟前
1
0
Android入门—activity生命周期

activity生命周期如下所示: onCreate()方法 创建activity时使用onStart()方法 当activity界面变为用户可见时调用onResume()方法 当activity界面获取到焦点时调用(界面按钮可点击,...

haoyuehong
21分钟前
1
0
使用C++ Builder 6编译产生错误

错误提示: [C++ Error] E2075 Incorrect project override option: Files\Borland\CBuilder6\lib\vcl60.csm[C++ Error] E2075 Incorrect project override option: Files\Borland\CBuilde......

simpower
24分钟前
1
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部