资源描述：
《直方图均衡化和直方图规定化-实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、一、实验目的掌握直方图均衡化和直方图规定化的图像增强方法掌握图像平滑滤波和锐化滤波的模板计算方法二、实验内容：1.使用IPT中imhist，histeq函数进行直方图的均衡化和规定化，并显示结果2.编写myhisteq函数实现直方图均衡化，与1中结果进行对比3.给读取的图像叠加椒盐噪音(imnoise)，分别使用均值滤波和中值滤波进行去噪，并对比图像处理的结果(使用IPT函数）4.自定义3*3模板矩阵F，编写myfilter函数实现模板和图像的卷积运算，设计模板矩阵，实现图像的平滑和锐化。三、实验代码及结果
2、(1)直方图的均衡化和规定化clc;clear;I=imread('H:imagejpgflower.jpg');I=rgb2gray(I);%将图像转换为灰度图像J=histeq(I);%对I进行直方图均衡化figure,subplot(2,2,1),imshow(I),title('原始图像');subplot(2,2,2),imshow(J),title('直方图均衡化后的图像');subplot(2,2,3),imhist(I,64),title('原始的直方图');subplot(2,2,4
3、),imhist(J,64),title('均衡化后的直方图');clc;clear;I=imread('H:imagejpgflower.jpg');I=rgb2gray(I);%将图像转换为灰度图像h=0:255;h=1-h/255;J=histeq(I,h);figure,subplot(2,3,1),imshow(I),title('原图像');subplot(2,3,2),imhist(I,64),title('原图像的直方图');subplot(2,3,3),stem(h),title('
4、目标直方图');subplot(2,3,4),imshow(I,64),title('规定化后的图像');subplot(2,3,5),imhist(J,64),title('规定化后的直方图');二、myhisteq函数实现直方图均衡化I=imread('j:imagejpgflower.jpg');I=rgb2gray(I);[height,width]=size(I);figuresubplot(2,2,1)imshow(I)%显示原始图像title('原图像');subplot(2,2,2)i
5、mhist(I)%显示原始图像直方图title('原图像直方图');%进行像素灰度统计;s=zeros(1,256);%统计各灰度数目，共256个灰度级fori=1:heightforj=1:widths(I(i,j)+1)=s(I(i,j)+1)+1;%对应灰度值像素点数量增加一endend%计算灰度分布密度p=zeros(1,256);fori=1:256p(i)=s(i)/(height*width*1.0);end%计算累计直方图分布c=zeros(1,256);c(1)=p(1);fori=2:2
6、56c(i)=c(i-1)+p(i);end%累计分布取整,将其数值归一化为1~256c=uint8(255.*c+0.5);%对图像进行均衡化fori=1:heightforj=1:widthI(i,j)=c(I(i,j)+1);endendsubplot(2,2,3)imshow(I)%显示均衡化后的图像title('均衡化后图像');subplot(2,2,4)imhist(I)%显显示均衡化后的图像的直方图title('均衡化后图像的直方图');三、使用均值滤波和中值滤波进行去噪I=imread('
7、j:imagejpgflower.jpg');I=rgb2gray(I);%将图像转换为灰度图像I1=imnoise(I,'salt&pepper',0.02);%0.02是噪声强度，其值越大噪声越多h=fspecial('average',5);J1=filter2(h,I1,'valid');J2=imfilter(I1,h,'full');figure,subplot(3,3,1),imshow(I1);title('原始椒盐噪声图像图1');subplot(3,3,2),imshow(uint
8、8(J1));title('filter2均值滤波图2');subplot(3,3,3),imshow(J2);title('imfilter均值滤波图3');J3=medfilt2(I1,[5,5]);subplot(3,3,4),imshow(J3),title('中值滤波效果图4');G1=histeq(I1);subplot(3,3,5),imhist(G1,64),title('原图均衡化后的直方图')
