用MATLAB实现直方图均衡化和中值滤波.docVIP

直方图均衡化直方图均衡化处理的“中心思想”是把原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布。直方图均衡化就是对图像进行非线性拉伸，重新分配图像像素值，使一定灰度范围内的像素数量大致相同。直方图均衡化就是把给定图像的直方图分布改变成“均匀”分布直方图分布。 2．直方图均衡化程序 clear all ; close all ; tuu=imread(picture1.JPG); tu=rgb2gray(tuu); %将彩色图片转换为灰度图 graydis=zeros(1,256); %设置矩阵大小 graydispro=zeros(1,256); new_graydis=zeros(1,256); new_graydispro=zeros(1,256); [h, w]=size(tu); new_tu=zeros(h,w); % 计算原始直方图各灰度级像素个数graydis for x=1:h for y=1:w graydis(1,(tu(x,y)+1))=graydis(1,(tu(x,y)+1))+1; end end %计算原始直方图graydispro graydispro=graydis./sum(graydis); subplot(1,2,1); plot(graydispro); title(灰度直方图); xlabel(灰度值); ylabel(像素的概率密度); %计算原始累计直方图 for i=2:256 graydispro(1,i)=graydispro(1,i)+graydispro(1,i-1); end %计算和原始灰度对应的新的灰度t[]，建立映射关系 for i=1:256 t(1,i)=floor(255*graydispro(1,i)+0.5); end %统计新直方图各灰度级像素个数new_graydis for i=1:256 new_graydis(1,t(1,i)+1)=new_graydis(1,t(1,i)+1)+graydis(1,i); end %计算新的灰度直方图new_graydispro new_graydispro=new_graydis./sum(new_graydis); subplot(1,2,2); plot(new_graydispro); title(均衡化后的灰度直方图); xlabel(灰度值);ylabel(像素的概率密度); %计算直方图均衡后的新图new_tu for x=1:h for y=1:w new_tu(x,y)=t(1,tu(x,y)+1); end end figure,imshow(tu,[]); title(原图); figure,imshow(new_tu,[]); title(直方图均衡化后的图); 3．程序运行结果 图1 原始图像 图2 直方图均衡化后的图 图3 原始灰度直方图和均衡化后的灰度直方图比较 中值滤波中值滤波法是一种非线性平滑技术，它将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值 中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术，中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替，让周围的像素值接近的真实值，从而消除孤立的噪声点。方法是用某种结构的二维滑动模板，将板内像素按照像素值的大小进行排序，生成单调上升（或下降）的为二维数据序列。二维滤波输出为g（x,y）=med{f(x-k,y-l),(k,lW)} ，其中，f(x,y)，g(x,y)分别为原始图像和处理后图像。W为二维模板，通常为2*2，3*3区域，也可以是不同的的形状，如线状，圆形，十字形，圆环形等。 clear all close all tuu=imread(picture1.JPG); tu=rgb2gray(tuu); J=imnoise(tu,salt pepper,0.02); imshow(J,[]); title(加入椒盐噪声后的图像); [h,w]=size(J); new_tu=J; k=zeros(1,9); for i=1:(h-2) for j=1:(w-2) k=[J(i,j) J(i+1,j) J(i+2,j) J(i,j+1) J(i+1,j+1) J(i+2,j+1) J(i,j+2) J(i+1,j+2) J(i+2,j+2)]; t=9;s=1; while t1 while t1