c语言实现直方图均衡化

最新推荐文章于 2023-06-17 20:45:50 发布
最新推荐文章于 2023-06-17 20:45:50 发布
阅读量1.2w 收藏 68
点赞数 8
分类专栏： c语言图像处理 文章标签： 直方图均衡化 图像处理 c语言 灰度直方图
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u011450295/article/details/38228117
版权

直方图均衡化部分是用c语言写的，最后用opencv显示原图像，处理后图像以及原图和处理后图的灰度直方图。

虽然做出来了，均衡化效果还可以，但不知道为什么处理后图像中有三条白线，真心搞不懂，有看出来问题的大神麻烦留言告诉我，谢谢。

（终于知道哪出问题了，原来是每行字节数求错了，改为LineByte=(width*8/8+3)/4*4;即可。）

实现原理参考：http://wenku.baidu.com/link?url=jEEUyr14TNX1B88qLrB0IMSOKMA-S8hNidKU2BqxmnEqnEgva6d18IC38NjlPE0OvCc7fsk2Pa3Uld6pj5YF-zY8yoes31r6HQBiNmLvutq

下面是代码：

#include "stdafx.h"
#include<stdio.h>
#include<windows.h>

#include<opencv2\highgui\highgui.hpp>
#include<opencv2\core\core.hpp>
#include<cv.h>
  
int main(void)
{
	int width;//图像宽度
	int height;//图像高度
	RGBQUAD *pColorTable;
	unsigned char *pBmpBuf,*pBmpBuf1;
	BITMAPFILEHEADER bfhead;
	BITMAPINFOHEADER bihead;

	FILE *fp1=fopen("e:\\picture\\dog.bmp","rb");
	if(fp1==0)
		return 0;
	fread(&bfhead,14,1,fp1);
	fread(&bihead,40,1,fp1);
	width=bihead.biWidth;
	height=bihead.biHeight;

	pColorTable=new RGBQUAD[256];
	fread(pColorTable,4,256,fp1);
	int LineByte=0;
	LineByte=(width*1/4+1)*4;
<span style="white-space:pre">	</span>//LineByte=(width*8/8+3)/4*4;
	pBmpBuf = new unsigned char[LineByte*height];

	fread(pBmpBuf,LineByte*height,1,fp1);
	fclose(fp1);

	pBmpBuf1=new unsigned char[LineByte*height]; //用于存储均值化后的图像数据
	//统计每个灰度级像素点的个数
	int N[256]={0};
	for(int i=0;i<height;i++)
		for(int j=0;j<width;j++)
		{
			unsigned char *pb1,*pb2;
			pb1=pBmpBuf+i*LineByte+j;
			N[*pb1]++;
			pb2=pBmpBuf1+i*LineByte+j;
			*pb2=*pb1;	
		}
	/*for(int i=0;i<256;i++ )
		printf("%d  ",N[i]);*/


		//统计最小与最大灰度值
		int minGrayValue=255;    
		int maxGrayValue=0;
	for(int i=0;i<height;i++)
		for(int j=0;j<width;j++)
		{
			unsigned char *pb;
			pb=pBmpBuf+i*LineByte+j;
			if(*pb>maxGrayValue)
				maxGrayValue=*pb;
			else if(*pb<minGrayValue)
				minGrayValue=*pb;
		}
		printf("%d ,%d\n",minGrayValue,maxGrayValue);//输出最大与最小灰度值

		int x=maxGrayValue-minGrayValue+1;
		float *p;
		p=new float[x];
		
		for(int i=0;i<x;i++)
		{
			*(p+i)=(float)N[i]/(float)(width*height);   //*(p+i)中存放的是灰度级为i的像素在整幅图像中出现
		                                                //的概率（即*（p+i）i=0,1,2,3...中存放的就是这幅图像归一化后的直方图）	
		}
		
		
		float *c;
		c=new float[x];      //定义c，用来存放累积的归一化直方图
		for(int i=0;i<x;i++)  //对c进行初始化
		{
			*(c+i)=0;

		}

		for(int i=0;i<x;i++)
		{
			for(int j=0;j<=i;j++)
			{
				*(c+i)+=*(p+j);
			}
		}

	

		

		for(int i=0;i<height;i++)
			for(int j=0;j<width;j++)
			{
				unsigned char *pb;
				pb=pBmpBuf1+i*LineByte+j;
				*pb=*(c+*pb)*(maxGrayValue-minGrayValue)+minGrayValue;
			}

		FILE *fp2=fopen("junhenghua.bmp","wb");
		fwrite(&bfhead,14,1,fp2);
		fwrite(&bihead,40,1,fp2);
		fwrite(pColorTable,4,256,fp2);
		fwrite(pBmpBuf1,LineByte*height,1,fp2);
		fclose(fp2);
		
		
//显示原图与处理后的图像
		IplImage *src1=cvLoadImage("e:\\picture\\dog.bmp");
		IplImage *src2=cvLoadImage("junhenghua.bmp");

		cvNamedWindow("原图");
		cvNamedWindow("处理后图");

		cvShowImage("原图",src1);
		cvShowImage("处理后图",src2);




//显示原图像与处理后图像的灰度直方图
	int size=256;
	float range[]={0,255};
	float *ranges[]={range};
	CvHistogram *hist1=cvCreateHist(1,&size, CV_HIST_ARRAY,ranges,1);//创建一维直方图，
	CvHistogram *hist2=cvCreateHist(1,&size, CV_HIST_ARRAY,ranges,1);
	IplImage* gray1=cvCreateImage(cvGetSize(src1),8,1);
	IplImage* gray2=cvCreateImage(cvGetSize(src2),8,1);
	cvCvtColor(src1,gray1,CV_BGR2GRAY); 
	cvCvtColor(src2,gray2,CV_BGR2GRAY);
	//vCvtColor(...),是Opencv里的颜色空间转换函数，可以实现RGB颜色向HSV,HSI等颜色空间的转换，也可以转换为灰度图像。
	//参数CV_RGB2GRAY是RGB到gray,
    //参数CV_GRAY2RGB是gray到RGB

	cvCalcHist(&gray1,hist1,0,0);//统计图像在[0 255]像素的均匀分布，将统计结果存在结构体中
	cvCalcHist(&gray2,hist2,0,0);

	//draw histogram-----

	//统计直方图中的最大直方块
	float histMax1=0,histMax2=0;
	cvGetMinMaxHistValue(hist1,0,&histMax1,0);     
	cvGetMinMaxHistValue(hist2,0,&histMax2,0);
	//创建一张一维直方图的“图”，横坐标为灰度级，纵坐标为像素个数  
	IplImage *grayHist1=cvCreateImage(cvSize(256*2,64*2),8,1);
	IplImage *grayHist2=cvCreateImage(cvSize(256*2,64*2),8,1);
	cvZero(grayHist1);
	cvZero(grayHist2);

	//分别将每个直方块的值绘制到图中  
	for(int i=0;i<255;i++)
	{
		float histValue=cvQueryHistValue_1D(hist1,i);
		float nextValue=cvQueryHistValue_1D(hist1,i+1);

		//计算直方块4个点的值
		CvPoint pt1=cvPoint(i*2,64*2);
		CvPoint pt2=cvPoint((i+1)*2,64*2);
		CvPoint pt3=cvPoint((i+1)*2,(64-(nextValue/histMax1)*64)*2);
		//nextValue/histMax是将i级像素点个数归一到0~1，在*64是使其高对在0~64之间
		//由于opencv图像是以左上角为坐标原点，向右为x轴，向下时y轴，而显示的直方图是向上增长的，所以用64减，将其倒过来显示
		CvPoint pt4=cvPoint(i*2,   (64-(histValue/histMax1)*64)*2);

		int ptNum=5;
		CvPoint pt[5];
		pt[0]=pt1;
		pt[1]=pt2;
		pt[2]=pt3;
		pt[3]=pt4;
		pt[4]=pt1;

        cvFillConvexPoly(grayHist1,pt,ptNum,cvScalar(255)); //填充直方块
		   

	}

	for(int i=0;i<255;i++)
	{
		float histValue=cvQueryHistValue_1D(hist2,i);
		float nextValue=cvQueryHistValue_1D(hist2,i+1);

		//计算直方块4个点的值
		CvPoint pt1=cvPoint(i*2,64*2);
		CvPoint pt2=cvPoint((i+1)*2,64*2);
		CvPoint pt3=cvPoint((i+1)*2,(64-(nextValue/histMax2)*64)*2);
		//nextValue/histMax是将i级像素点个数归一到0~1，在*64是使其高对在0~64之间
		//由于opencv图像是以左上角为坐标原点，向右为x轴，向下时y轴，而显示的直方图是向上增长的，所以用64减，将其倒过来显示
		CvPoint pt4=cvPoint(i*2,   (64-(histValue/histMax2)*64)*2);

		int ptNum=5;
		CvPoint pt[5];
		pt[0]=pt1;
		pt[1]=pt2;
		pt[2]=pt3;
		pt[3]=pt4;
		pt[4]=pt1;

        cvFillConvexPoly(grayHist2,pt,ptNum,cvScalar(255)); //填充直方块
		   

	}

	cvNamedWindow("grayHistogram1");
	cvNamedWindow("grayHistogram2");
	cvShowImage("grayHistogram1",grayHist1);
	cvShowImage("grayHistogram2",grayHist2);


	

		cvWaitKey(0);

	system("pause");
	
return 0;  
}

原图：

处理后图：

原图直方图：


均衡化后直方图：


@leozhang
关注
  • 8
    点赞
  • 踩
  • 68
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
专栏目录
C语言实现直方图均衡化
wormtree的专栏
03-09 8574
#include #include #include using namespace std; int main(int argc, char** argv) { //载入图片 int i=0, j=0, temp=0; IplImage * img = cvLoadImage("D:\\2.jpg", CV_LOAD_IMAGE_GRAYSCALE);//图片路径 int hei
直方图均衡化——C语言实现
02-12
用C语言实现了直方图均衡化，内有测试图片。非常好用。
直方图均衡化的c语言实现
09-07
本例实现了直方图均衡化的算法，很简单，代码一共就80多行。代码就完全是压缩包里面文档的实现，对于初学者来说很有用。
用c语言实现图像的直方图均衡化_C# 直方图均衡
weixin_39669133的博客
11-23 615
概念直方图均衡化是图像处理领域中利用图像直方图对对比度进行调整的方法。通过这种方法，亮度可以更好地在直方图上分布。这样就可以用于增强局部的对比度而不影响整体的对比度，直方图均衡化通过有效地扩展常用的亮度来实现这种功能。直方图均衡_百度百科​baike.baidu.com介绍一个比较容易理解的链接内容https://blog.csdn.net/qq_39119348/article/details/...
直方图均衡化 C语言源代码
01-12
把原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布。对图像进行非线性拉伸，重新分配图像像素值，使一定灰度范围内的像素数量大致相同。把给定图像的直方图分布改变成“均匀”分布直方图分布。
直方图均衡化_c语言实现_源代码+实验报告+软件文档
06-04
本实验由北京交通大学计算机学院的阮秋琦教授指导，旨在让学生掌握C语言实现直方图均衡化的原理与方法。 直方图是描述图像中像素强度分布的一种统计图形，每个强度值对应的柱高表示该强度值出现的频次。在图像对比...
直方图均衡化
06-25
在C语言中实现直方图均衡化，通常包括以下几个步骤： 1. **计算直方图**：首先，我们需要遍历图像中的每一个像素，统计每个灰度级的像素数量，构建出原始的直方图。 2. **累积分布函数(CDF)**：直方图的累积分布...
用C实现直方图均衡化
10-10
在C语言中实现直方图均衡化，我们需要理解几个关键概念和技术。 首先，直方图是表示图像中像素亮度或颜色分布的统计图形。在灰度图像中，直方图由横坐标代表灰度级，纵坐标代表对应灰度级出现的频次。如果图像的...
c语言直方图均衡化程序,直方图均衡化程序 c++
weixin_28769989的博客
05-17 574
#include #include #include #include using namespace cv;using namespace std;Mat calchistcontrol(vector rgb);void draw(Mat hist, int histSize, Mat histImage, int bin_w, int hist_h, int rgb);Mat getHistI...
图像处理中的基础函数：直方图均衡化纯C代码
12-26
图像处理技术的基础函数，进行图像直方图均衡化的必要函数。代码质量很高，可作为C程序代码编写的规范模板。
数字图像处理-直方图均衡化c语言代码
04-24
基于VC++平台，带图片，绝度能运行，利用画图都方法，逐点描图，细致分析bmp结构及原理，全部利用代码进行手工描图。
图像处理C语言(读入图像，直方图均衡化，灰度线性变化，灰度拉伸)
03-26
用C语言写的图像处理程序，包括读入8位图像，对图像进行直方图均衡化，灰度线性变化和灰度线性拉伸。
灰度图像的直方图均衡化程序C语言版
05-26
将读入灰度图像进行直方图均衡化，使得原来比较暗的图像变的比较亮，就是对比度增加。
c语言，matlab语言实现直方图均衡化
11-25
见过最好的c语言，matlab语言实现直方图均衡化程序代码，完美通过 调试，希望对你理解有帮助
Intensity Transformations and HistogramEqualization
weixin_38825015的博客
01-15 291
1.1 对数变换方法进行图像增强的程序duishu.m function duishu(T,C) %对数变换对图像增强的效果 m=1;n=2; T=imread('Fig3.08(a).jpg'); C=35;%改变参变量C的值来探索最好的图像处理效果 double C; [x,y]=size(T); for i=1:x for j=1:y T1(i,j)=uint8(C*log(double
scratch lenet(3): 直方图均衡化的C语言实现
最新发布
baiyu33的博客
06-17 933
使用C语言实现直方图均衡化，并基于.pgm格式的单通道灰度图进行验证，比对原图和均衡化后的图像，比对原图的直方图和均衡化后的直方图。
c语言实现直方图均衡
G_peng的博客
10-27 958
#include<stdio.h> #include<malloc.h> #include<stdlib.h> #include <math.h> /* 位图头结构 */ #pragma pack(1) typedef struct tagBITMAPFILEHEADER { unsigned char bfType[2];//文件格式 unsigned long bfSize;//文件大小 unsigned short bfReser
c语言图像压缩算法实现_C语言实现BMP图像处理：直方图均衡化
weixin_39797532的博客
11-23 708
计算步骤：1)统计各灰度值的概率；2)计算了累积概率直方图(CDF)；3)取整扩展：Tk = int[ (L-1)*Tk];代码#include #include #include #include int main(int* argc, char** argv){FILE* fp = fopen("./01.bmp", "rb");if (fp == 0)return 0;BITMAPFILE...
如何用c语言实现直方图均衡化
06-01
直方图均衡化是一种图像增强的方法，可以使得图像的灰度变化更加均匀，从而提高图像的对比度和清晰度。下面是用C语言实现直方图均衡化的代码： ```c #include <stdio.h> #include <stdlib.h> #include <math.h> #define WIDTH 256 #define HEIGHT 256 int main(int argc, char *argv[]) { int histogram[WIDTH] = {0}; //输入图像的直方图 int cumulative_histogram[WIDTH] = {0}; //输入图像的累积直方图 int output_histogram[WIDTH] = {0}; //输出图像的直方图 int input_image[HEIGHT][WIDTH] = {0}; //输入图像 int output_image[HEIGHT][WIDTH] = {0}; //输出图像 int pixel_count = 0; int i, j, k; //读取输入图像 FILE *input_file = fopen("input_image.raw", "rb"); fread(input_image, sizeof(unsigned char), HEIGHT * WIDTH, input_file); fclose(input_file); //计算输入图像的直方图 for (i = 0; i < HEIGHT; i++) { for (j = 0; j < WIDTH; j++) { histogram[input_image[i][j]]++; } } //计算输入图像的累积直方图 cumulative_histogram[0] = histogram[0]; for (i = 1; i < WIDTH; i++) { cumulative_histogram[i] = cumulative_histogram[i - 1] + histogram[i]; } //计算输出图像的直方图 pixel_count = WIDTH * HEIGHT; for (i = 0; i < WIDTH; i++) { output_histogram[i] = round((float)(cumulative_histogram[i] * (WIDTH - 1)) / pixel_count); } //计算输出图像 for (i = 0; i < HEIGHT; i++) { for (j = 0; j < WIDTH; j++) { output_image[i][j] = output_histogram[input_image[i][j]]; } } //保存输出图像 FILE *output_file = fopen("output_image.raw", "wb"); fwrite(output_image, sizeof(unsigned char), HEIGHT * WIDTH, output_file); fclose(output_file); return 0; } ``` 在这个代码中，我们首先定义了一些常量，包括输入图像的大小和灰度级数。然后，我们读取输入图像，并计算输入图像的直方图和累积直方图。接着，我们根据累积直方图计算输出图像的直方图，并将其应用到输入图像上，生成输出图像。最后，我们将输出图像保存到文件中。 需要注意的是，这个代码中并没有进行图像预处理，如图像的归一化和平滑等。在实际应用中，这些预处理步骤可能是必要的。
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值