[code]彩色图像直方图均衡化 histogram_rgb

最新推荐文章于 2023-12-09 20:15:57 发布

eaglediao123

最新推荐文章于 2023-12-09 20:15:57 发布

阅读量1k

点赞数

分类专栏： c/cpp

本文链接：https://blog.csdn.net/eaglediao123/article/details/26506071

版权

c/cpp 专栏收录该内容

14 篇文章 0 订阅

订阅专栏

//2013.9 eageldiao
#ifdef HISTOGRAM_RGB
unsigned int lut[256];
unsigned intncount[256]={0},ncount1[256]={0},ncount2[256]={0};
int nTemp;
//b
for(y=0;yheight;y++)
{
unsigned char *srcrow= (unsignedchar*)(src->imageData+y*src->widthStep);
for (x=0;xwidth;x++)
{   
ncount[srcrow[3*x]]++;               //统计灰度级数量
}
}
for (int i=0;i<256;i++)             
{
nTemp=0;
for (int j=0;j<=i;j++)
{
nTemp+=ncount[j];
}
lut[i]=nTemp*255/src->width/src->height;//确定变换函数
}
for(y=0;yheight;y++)               //均衡化
{
unsigned char *srcrow= (unsignedchar*)(src->imageData+y*src->widthStep);
for (x=0;xwidth;x++)
{  
srcrow[3*x]=lut[srcrow[3*x]];
}
}
//g
for(y=0;yheight;y++)                //统计灰度级 数量
{
unsigned char *srcrow= (unsignedchar*)(src->imageData+y*src->widthStep);
for (x=0;xwidth;x++)
{   
ncount1[srcrow[3*x+1]]++;
}
}

for (int i=0;i<256;i++)             
{
nTemp=0;
for (int j=0;j<=i;j++)
{
nTemp+=ncount1[j];
}
lut[i]=nTemp*255/src->width/src->height; //确定变换函数
}

for(y=0;yheight;y++)                 //均衡化
{
unsigned char *srcrow= (unsignedchar*)(src->imageData+y*src->widthStep);
for (x=0;xwidth;x++)
{  
srcrow[3*x+1]=lut[srcrow[3*x+1]];
}
}

//r
for(y=0;yheight;y++)                  //统计灰度级 数量
{
unsigned char *srcrow= (unsignedchar*)(src->imageData+y*src->widthStep);
for (x=0;xwidth;x++)
{   
ncount2[srcrow[3*x+2]]++;
}
}
for (int i=0;i<256;i++)             
{
nTemp=0;
for (int j=0;j<=i;j++)
{
nTemp+=ncount2[j];
}
lut[i]=nTemp*255/src->width/src->height;//确定变换函数
}
for(y=0;yheight;y++)                //均衡化
{
unsigned char *srcrow= (unsignedchar*)(src->imageData+y*src->widthStep);
for (x=0;xwidth;x++)
{  
srcrow[3*x+2]=lut[srcrow[3*x+2]];
}
}
#endif