本文目录一览：
- 1、直方图均衡化步骤？
- 2、谁有部分子块重叠局部直方图均衡的程序，即poshe算法，谢谢，下面是poshe算法的步骤
- 3、用c语言编写一个程序,打印输入单词的长度的直方图?
- 4、求一个直方图均衡化的程序
- 5、用C语言编写灰度图像直方图变换增强
- 6、C语言绘制直方图(C程序设计语言习题1-13）
直方图均衡化步骤？
这个足够了 有问题再留言
clear all
%一，图像的预处理，读入彩色图像将其灰度化
PS=imread(‘s7.jpg’); %读入JPG彩色图像文件 ，注意路径
figure(1);subplot(2,2,1);imshow(PS);title(‘原图像灰度图’);
%二，绘制直方图
[m,n]=size(PS); %测量图像尺寸参数
GP=zeros(1,256); %预创建存放灰度出现概率的向量
for k=0:255
GP(k+1)=length(find(PS==k))/(m*n); %计算每级灰度出现的概率，将其存入GP中相应位置
end
figure(1);subplot(2,2,2);bar(0:255,GP,’g’) %绘制直方图
title(‘原图像直方图’)
xlabel(‘灰度值’)
ylabel(‘出现概率’)
%三，直方图均衡化
S1=zeros(1,256);
for i=1:256
for j=1:i
S1(i)=GP(j)+S1(i); %计算Sk
end
end
S2=round((S1*256)+0.5); %将Sk归到相近级的灰度
for i=1:256
GPeq(i)=sum(GP(find(S2==i))); %计算现有每个灰度级出现的概率
end
figure(1);subplot(2,2,4);bar(0:255,GPeq,’b’) %显示均衡化后的直方图
title(‘均衡化后的直方图’)
xlabel(‘灰度值’)
ylabel(‘出现概率’)
%四，图像均衡化
PA=PS;
for i=0:255
PA(find(PS==i))=S2(i+1); %将各个像素归一化后的灰度值赋给这个像素
end
figure(1);subplot(2,2,3);imshow(PA) %显示均衡化后的图像
title(‘均衡化后图像’)
imwrite(PA,’PicEqual.bmp’);
谁有部分子块重叠局部直方图均衡的程序，即poshe算法，谢谢，下面是poshe算法的步骤
clear;
I = imread(‘E:\ebook\lena.bmp’);
figure,imshow(I),title(‘original image’);
len = 5;
hstep = 2;
wstep = 2;
for k = 1:size(I,3)
img = I(:,:,k);
img = [img(:,len:-1:2) img img(:,end-1:-1:end-len)];
img = [img(len:-1:2,:) ;img; img(end-1:-1:end-len,:)];%%其实就是图像扩展，参考wextend函数，当然你的方法也可以
[h w] = size(img);
med_out(:,:,k) = int16(zeros(h,w));%%转化成int16型数据，原因下面会说明
counts = zeros(h,w);
for i =1:hstep:h
if h-i = len
for j = 1:wstep:w
if w-j = len
rect =int16(histeq(img(i:i+len-1,j:j+len-1)));%同样转化成int16类型
med_out(i:i+len-1,j:j+len-1,k) =med_out(i:i+len-1,j:j+len-1,k)+rect;%%rect和med_out转化成int16类型原因在于，你的med_out计算直方图均衡化并且把重复的像素的直方图均衡化结果进行累加，那么其累加的值很有可能超过255，所以你转化med_out为uint8型是不够的，它会把超过255的值强制为255。
counts(i:i+len-1,j:j+len-1) = counts(i:i+len-1,j:j+len-1)+1;
end
end
end
end
med_out1(:,:,k) = double(med_out(len:end-len,len:end-len,k));%%这里对图像两边各扩展了四个像素，所以起始位置应该是len，而不是len+1,同样结束就是end-len
med_out1(:,:,k) = med_out1(:,:,k)./counts(len:end-len,len:end-len);
out(:,:,k) = uint8(med_out1(:,:,k));%%转化成uint8型
end
figure,imshow(out),title(‘poshe result’);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
这样改就可以了，运行后可以出结果，图片就不贴了！！
编程时要考虑数据范围是否过界的问题。
用c语言编写一个程序,打印输入单词的长度的直方图?
/*
练习1.13；
编写一个程序，打印输入中长度的直方图。
*/
# include stdio.h
# define MAXHIST 15 //最大长度的直方图
# define MAXWORD 11 //一个单词的最大长度
# define IN 1 //文字内
# define OUT 0 //外一个单词
/*打印水平直方图*/
int main(void)
{
int c,i,nc,state;
int len; //每条的长度
int maxvalue; //最大者为
int ovflow; //数量的溢出的话
int wl[MAXWORD]; //单词长度计数器
state = OUT;
nc = 0; //在一个单词字符的数量
ovflow = 0; //字数 = MAXWORD
for(i = 0;i MAXWORD; ++i)
wl[i] = 0;
while ((c = getchar()) != EOF)
{
if (c == ‘ ‘ || c == ‘\n’ || c == ‘\t’)
{
state = OUT;
if(nc 0)
if(nc MAXWORD) //如果nc小于11则wl[nc]元素加一
++wl[nc];
else //否则ovflow加一
++ovflow;
nc = 0;
}
else if (state == OUT)
{
state = IN;
nc = 1; //开始一个新单词
}
else
++nc; //文字内
}
maxvalue = 0;
for (i = 1; i MAXWORD; ++i)
if(wl[i] maxvalue)
maxvalue = wl[i];
for (i = 1; i MAXWORD; ++i)
{
printf(“%5d – %5d : “,i,wl[i]);
if(wl[i] 0)
{
if ((len = wl[i] * MAXHIST / maxvalue) = 0)
len = 1;
}
else
len = 0;
while (len 0)
{
putchar(‘*’);
–len;
}
putchar(‘\n’);
}
if (ovflow 0)
printf(“There are %d words = %d\n”,ovflow,MAXWORD);
return 0;
}
求一个直方图均衡化的程序
图形处理中有一种对比度变换，像显示器就有对比度调节，PhotoShop也有图片的对比度修改，对比度的提高可以使图像细节清晰，相反，对比度的减小可以隐藏图像的细节，在一定程度上使图像柔和。
对比度变换其中一种比较简单的方法是直方图均衡化。
所谓直方图就是在某一灰度级的象素个数占整幅图像的象素比 h=nj/N,其中nj是灰度级在j的象素数，N是总象素数，扫描整幅图像得出的h的离散序列就是图像的直方图，h求和必然＝1，所以直方图可以看成是象素对于灰度的概率分布函数。
直方图是高低不齐的，因为象素灰度是随机变化的，直方图均衡化就是用一定的算法使直方图大致平和。
算法如下：
对于一个直方图
设 Pr(r)是原始图像直方图，Ps(s)是均衡化的直方图，
由于其是一个概率分布函数
所以有 Ps(s)ds=Pr(r)dr (编辑关系，ds,dr是积分变量)
因为要进行均衡化，令 Ps(s)=1,
得 ds=Pr(r)dr/1
两边积分得 s=F Pr(r)dr (因为编辑关系，左边F表示积分符号….-__-++)
数字图像是离散的，因此离散化上式得
sk=E{j=0,k}(nj/N) 左式k,j是离散量下标，因为编辑关系，E{0,k}表示下标0到k的连加符号，N是象素总数
由此得出每一象素的sk为均衡化后的正规化灰度（即灰度正规化到[0,1])，统计sk即可得出均衡化后的直方图。
在均衡化过程中可以对每一象素映射到新的实际灰度值sk*255，就实现了图像的变换
（严格理论中应该是灰度正规化到[0,1]区间，然后均衡化后的sk还要量化到原始的正规灰度以实现灰度合并，下面的BCB程序并没有量化，而且255是固定灰度级，因为256色BMP的彩色表就是256个表项)
现在开始实践
用BCB对一BMP灰度图像进行直方图均衡化处理，代码如下
//—————————-BCB6代码
#include vcl.h
#pragma hdrstop
#includestdio.h
#include “Unit1.h”
#include”File1.h”
#pragma pack(1)
//BMP文件头
struct BITMAPFILEHEADER_
{
short type;
int bfSize;
short re1,re2;
int Offbits;
};
//BMP信息头
struct BITMAPINFO_
{
long size;
long width,height;
short planes,bitCount;
long comp,sizeImg;
long xpels,ypels;
long used,important;
};
//BMP彩色表项
struct COLOR_
{
char blue,green,red,re;
};
//——将BMP彩色表的数据校正到BCB TColor的数据。
void SwitchColor(long c)
{
long blue=c 0x000000ff;
long green=c 0x0000ff00;
long red=c 0x00ff0000;
c=(blue16) | green | (red16);
}
void xxx()
{
FILE *f=fopen(“f:\\bbs_prev2.bmp”,”rb”);
if(f==NULL) /*判断文件是否打开成功*/
{
ShowMessage(“File open error”);
return;
}
fseek(f,0,0);//移动到开头
//———-读BMP文件头
BITMAPFILEHEADER_ *bmph=new BITMAPFILEHEADER_();
if(fread((char*)bmph,sizeof(BITMAPFILEHEADER_),1,f)==NULL)
{
ShowMessage(“File read error”);
return;
}
//———–读BMP信息头
BITMAPINFO_ *bmpi=new BITMAPINFO_();
if(fread((char*)bmpi,sizeof(BITMAPINFO_),1,f)==NULL)
{
ShowMessage(“File read error2”);
return;
}
//————–读彩色表
long *c=new long[bmph-Offbits-sizeof(BITMAPFILEHEADER_)-sizeof(BITMAPINFO_)];
fread((char*)c,bmph-Offbits-sizeof(BITMAPFILEHEADER_)-sizeof(BITMAPINFO_),1,f);
//———-显示一些信息
Form1-Edit1-Text=IntToStr(bmph-bfSize);
Form1-Edit2-Text=IntToStr(bmpi-width);
Form1-Edit3-Text=IntToStr(bmpi-height);
Form1-Edit4-Text=IntToStr(bmpi-comp);
Form1-Edit5-Text=IntToStr(bmpi-used);
int i,j,k,wc;
long N=bmph-bfSize- bmph-Offbits;//象素总数
unsigned char *image=new char[N]; //位图矩阵
unsigned char *newimage=new char[N];//变换后的位图矩阵
fread(image,N,1,f);//读入位图矩阵
//———直方图数列初始化
//———直方图数列用来存储正规化后的灰度
double *h=new double[255];//255个灰度级，保存原始图像正规化灰度直方图数据
for(i=0;i255;i++)
h[i]=0.0;
double *nh=new double[255];//255个灰度级，保存变换后的图像正规化灰度直方图
for(i=0;i255;i++)
nh[i]=0.0;
long *count=new long[255]; //每一灰度级的象素数量统计
for(i=0;i255;i++)
count[i]=0;
for(i=0;iN;i++)
{
count[image[i]]++;
}
//－－－－－正规化灰度概率统计
for(i=0;i255;i++)
{
h[i]=count[i]/(double)N;
}
//——正规化新灰度图
double hc;
for(i=0;iN;i++)
{
hc=0;
for(j=0;jimage[i];j++)
hc+=h[j];
nh[image[i]]+=hc; //保存新正规化灰度图
newimage[i]=hc*255; //保存新图像灰度索引
}
//———-显示直方图
for(i=0;i255;i++)
{
//原始直方图
Form1-Canvas-MoveTo(10+i,200);
Form1-Canvas-LineTo(10+i,200+h[i]*N);
//新直方图
Form1-Canvas-MoveTo(300+i,200);
Form1-Canvas-LineTo(300+i,200+nh[i]*255);
}
//——显示图形
TColor *tc;
if(bmpi-width%4==0)//———–因为BMP图像4字节对齐
wc=bmpi-width/4*4;
else
wc=(bmpi-width/4+1)*4;
long a;
long pos=0;
for( i=0;ibmpi-height;i++)
{
for(j=0;jwc;j++)
{
//—–原始图形
a= c[image[pos]];
SwitchColor(a);
Form1-Canvas-Pixels[10+j][600-i]=a;
//——新图形
a= c[newimage[pos]];
SwitchColor(a);
Form1-Canvas-Pixels[300+j][600-i]=a;
pos++;
}
}
fclose(f);
}
这个程序使用256色BMP文件，但程序代码是针对灰度图像的，用于彩色图像时得出一些古怪色彩配合而已。
在对灰度图像均衡化时
如果原始图像对比度本来就很高，如果再均衡化则灰度调和，对比度降低。
在泛白缓和的图像中，由于均衡化过程中会合并一些象素灰度，则会增大对比度，这里255灰度级太多，合并不明显。
用C语言编写灰度图像直方图变换增强
参考一下这个程序：
#define IMAGEWIDTH 80
#define IMAGEHEIGHT 80
void Histogram(unsigned char *pImage,int nWidth,int nHeight,float fHisto[256]);
void Enhance(unsigned char *pImage,unsigned char *pImage1,int nWidth,int nHeight,float fHisto[256],float fHisto1[256]);
unsigned char dbImage[IMAGEWIDTH*IMAGEHEIGHT];
unsigned char dbTargetImage[IMAGEWIDTH*IMAGEHEIGHT];
float fHistogram[256],fHistogram1[256];
unsigned char lut[256];
/* 直方图统计程序 */
int main()
{
… //准备一幅灰度图片到dbImage
Histogram(dbImage,IMAGEWIDTH,IMAGEHEIGHT,fHistogram); //做直方图统计
Enhance(dbImage,dbTargetImage,IMAGEWIDTH,IMAGEHEIGHT,fHistogram,fHistogram1); //直方图增强输出到dbTargetImage
…//完成后显示
}
void Histogram(unsigned char *pImage,int nWidth,int nHeight,float fHisto[256])
{
int i,j;
unsigned int uWork;
unsigned char *pWork;
for ( i=0;i256;i++ ) fHisto[i]=0.0f;
pWork=pImage;
for ( i=0;inHeight;i++ )
{
for ( j=0;jnWidth;j++,pWork++ )
{
uWork=(unsigned int)(*pWork);
fHisto[uWork]++;
}
}
uWork=nWidth*nHeight;
for ( i=0;i256;i++ )
{
fHisto[i]/=uWork;
fHisto[i]*=100;
}
}
void Enhance(unsigned char *pImage,unsigned char *pImage1,int nWidth,int nHeight,float fHisto[256],float fHisto1[256])
{
int i,j;
unsigned int uWork;
unsigned char *pWork,*pWork1;
for ( i=0;i256;i++ )
fHisto1[i]=fHisto[i]/100;
for ( i=1;i256;i++ )
fHisto1[i]+=fHisto1[i-1];
for ( i=0;i256;i++ )
lut[i]=fHisto1[i]*256;
for ( i=0;i256;i++ )
if ( lut[i]=256 )
lut[i]=255;
pWork=pImage; pWork1=pImage1;
for ( i=0;inHeight;i++ )
for ( j=0;jnWidth;j++,pWork++,pWork1++ )
(*pWork1)=lut[(*pWork)];
}
C语言绘制直方图(C程序设计语言习题1-13）
你好，大概看了一下，这个好像是统计一段文字中单词的长度并绘制直方图。这里wl就是用来记录单词长度的。wl[i] = j 代表长度为i的单词有j个。
比如最开始，wl数组全是0：
wl = [0, 0, 0, 0, …]
这个时候我们读入了一个长度为2的单词，于是wl[2]++，数组变成：（注意数组下标从0开始）
wl = [0, 0, 1, 0, …]
如此往复就记录了所有单词的长度统计。
原创文章，作者：小蓝，如若转载，请注明出处：https://www.506064.com/n/257199.html
微信扫一扫 
支付宝扫一扫 