直方图均衡实验报告.pdf

《直方图均衡实验报告.pdf》由会员分享，可在线阅读，更多相关《直方图均衡实验报告.pdf（26页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数字图像处理实验报告 实验二 图像直方图均衡 姓名：*学号：*专业：*一 实验目的 学习并掌握图像直方图均衡的基本理论，并通过分析均衡前后的图像验证课堂教学内容，总结直方图均衡的特点。二 实验内容 对图像 img2 进行直方图均衡（img2 为 X 光片图像）1 对比均衡前后图像的直方图及特点，1 图形显示其直方图以及灰度映射关系 2 计算以下参数：概率非零灰度数，概率非零灰度中最大概率、最小概率、最大最小概率之比。3 统计图像中概率大于平均概率的灰度级数；4 计算非零概率的平均值和方差 5 您认为哪些参数能够表现图像直方图分布的均匀程度？试提取相关参数进行测试。2找一到两幅图像重复上述实验。

2、3通过实验结果对比，能得出什么结论？三、实验报告要求 1源程序（C 或 Matlab）：（1）.画直方图的 Matlab 程序 filname=E:课件大三下图像处理实验二1.bmp;/打开文件 A=imread(filname);/读取文件信息 imhist(A);/画出图像文件的直方图 （2）C 程序#include#include#include#include BITMAPFILEHEADER bfh;BITMAPINFOHEADER bih;typedef struct PP unsigned char GRAY;PIXEL;/画出灰度映射关系图/void draw(double s

3、)FILE*fout,*fin;int i,j;PIXEL p256256;BITMAPFILEHEADER m_bfh;BITMAPINFOHEADER m_bih;RGBQUAD colorPanel256;fin=fopen(1.bmp,rb);fread(&m_bfh,1,sizeof(BITMAPFILEHEADER),fin);fread(&m_bih,1,sizeof(BITMAPINFOHEADER),fin);fread(colorPanel,256,sizeof(RGBQUAD),fin);fclose(fin);fout=fopen(Gray mapping relati

4、onship.bmp,wb);m_bih.biWidth=256;m_bih.biHeight=256;m_bih.biSizeImage=256*256;m_bfh.bfSize=m_bfh.bfOffBits+m_bih.biSizeImage;for(i=0;i256;i+)for(j=0;j256;j+)pij.GRAY=255;for(i=0;i256;i+)j=(int)si;pij.GRAY=0;fwrite(&m_bfh,1,sizeof(BITMAPFILEHEADER),fout);fwrite(&m_bih,1,sizeof(BITMAPINFOHEADER),fout)

5、;fwrite(colorPanel,256,sizeof(RGBQUAD),fout);fwrite(p,256*256,1,fout);fclose(fout);/统计概率非零灰度数子函数/double nonzero(double p)int i;double sum=0;for(i=0;i256;i+)if(pi!=0)sum=sum+1;printf(概率非零灰度数:%.2fn,sum);return sum;/找到最大概率子函数/double lmax(double p)int i;double max;for(i=0;i256;i+)/对 max 初始化 if(pi!=0)max

6、=pi;break;for(i=0;imax)max=pi;/找出概率非零最大值 printf(概率非零最大值:%lfn,max);return max;/找到最小概率子函数/double lmin(double p)int i;double min;for(i=0;i256;i+)/对 min 初始化 if(pi!=0)min=pi;break;for(i=0;i256;i+)if(pi!=0&pimin)min=pi;/找出概率非零最小值 printf(概率非零最小值:%lfn,min);return min;/求最大最小概率之比子函数/double rate(double max,dou

7、ble min)double r;r=max/min;printf(最大最小概率之比:%lfn,r);return r;/统计图像中概率大于平均概率的灰度级数子函数/double lager(double p)int i;double sum=0;for(i=0;i1.0/256)sum=sum+1;printf(图像中概率大于平均概率的灰度级数:%lfn,sum);return sum;/求非零概率的平均值的子函数/double average(double sum)double aver;aver=1.0/sum;printf(非零概率的平均值:%lfn,aver);return aver

8、;/求非零概率的方差的子函数/double variance(double sum,double p,double aver)int i;double asum=0,var;for(i=0;i256;i+)if(pi!=0)asum=asum+pow(pi-aver),2);/printf(%lfn,asum);var=asum/sum;printf(非零概率的方差:%lfn,var);return var;/主函数/int main()FILE*fin,*fout;RGBQUAD colorPanel256;PIXEL*bmpData;long i,j;double num256=0,p25

9、6=0,s256=0;double max,min,r,sum=0,asum=0,aver,psum=0,var;fin=fopen(1.bmp,rb);fread(&bfh,1,sizeof(BITMAPFILEHEADER),fin);fread(&bih,1,sizeof(BITMAPINFOHEADER),fin);fread(colorPanel,256,sizeof(RGBQUAD),fin);bmpData=(PIXEL*)malloc(bih.biHeight*sizeof(PIXEL*);for(i=0;ibih.biHeight;i+)bmpDatai=(PIXEL*)ma

10、lloc(bih.biWidth*sizeof(PIXEL);fread(bmpDatai,bih.biWidth,sizeof(PIXEL),fin);fclose(fin);/对图像进行直方图均衡/for(i=0;ibih.biHeight;i+)/二重循环遍历图象数据 for(j=0;jbih.biWidth;j+)numbmpDataij.GRAY+;/当前灰度级象素/总象素 for(i=0;i256;i+)pi=numi/bih.biWidth/bih.biHeight;/计算灰度概率 s0=p0;for(i=1;i256;i+)si=si-1+pi;/计算累积概率（低灰度概率加到高

11、灰度概率上）for(i=0;i256;i+)/01 映射到 0255 si=si*255;for(j=0;j=j-0.5&sij+0.5)si=j;for(i=0;ibih.biHeight;i+)/对原图像进行直方图均衡 for(j=0;jbih.biWidth;j+)bmpDataij.GRAY=(unsigned char)sbmpDataij.GRAY;/计算各个参数/sum=nonzero(p);/概率非零灰度数 max=lmax(p);/概率非零灰度中最大概率 min=lmin(p);/概率非零灰度中最小概率 r=rate(max,min);/最大最小概率之比 asum=lager

12、(p);/图像中概率大于平均概率的灰度级数 aver=average(sum);/非零概率的平均值 var=variance(sum,p,aver);/非零概率的方差 draw(s);/画灰度映射关系曲线 fout=fopen(3.bmp,wb);fwrite(&bfh,1,sizeof(BITMAPFILEHEADER),fout);fwrite(&bih,1,sizeof(BITMAPINFOHEADER),fout);fwrite(colorPanel,256,sizeof(RGBQUAD),fout);for(i=0;ibih.biHeight;i+)fwrite(bmpDatai,b

13、ih.biWidth,sizeof(PIXEL),fout);fclose(fout);return 0;2.程序流程图：3实验结果图像及对应参数：（1）图一：均衡前：均衡前直方图：灰度映射：对应参数：均衡后：均衡后的直方图：对应参数：均衡前后灰度映射关系：横轴表示均衡前的灰度，纵轴表示均衡后的灰度。（2）图二：均衡前：均衡前直方图：对应的参数：均衡后：均衡后直方图：对应参数：均衡前后灰度映射关系：横轴表示均衡前的灰度，纵轴表示均衡后的灰度。（3）图三：均衡前：均衡前直方图：对应参数：均衡后：均衡后直方图：对应参数：均衡前后灰度映射关系：横轴表示均衡前的灰度，纵轴表示均衡后的灰度。（4）图四：

14、均衡前：均衡前直方图：对应参数：均衡后：均衡后直方图：对应参数：均衡前后灰度映射关系：横轴表示均衡前的灰度，纵轴表示均衡后的灰度。（5）图五：均衡前：均衡前直方图：对应参数：均衡后：均衡后直方图：对应参数：均衡前后灰度映射关系：横轴表示均衡前的灰度，纵轴表示均衡后的灰度。4实验结果分析及结论 由实验结果的参数分析，直方图均衡前后有几项参数变化较大。均衡后概率非零的灰度数减少了，表示图像灰度分布得更均匀了。同时均衡后概率非零的最大值和最小值的比明显减小了，从这项数据中可以看出图像数据确实是被均衡了。而且非零概率的方差和均值都变大了，也就是亮暗区域之间的变化大了，使得图像更清晰了。直方图均衡能够使图像暗区变得更暗，亮区变得更亮，因此可以使一些看上去较暗的图变得更加容易辨认，视觉效果比较好。并且通过比较种子图，可以看出，直方图均衡对那些整体画面比较暗或者看上去颜色比较浅的图效果比较好，可以使图像看上去变得更亮，对本身就亮暗较分明的图像效果不明显，由 X 光图片可以看出，直方图均衡可以使图像更好的显示骨骼结构中的。