图像工程-灰度线性化-均衡化

《图像工程-灰度线性化-均衡化》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、一、主程序class_ex.m(一)流程图(二)程序代码clearMyImg_origMyImg_antiMyImg_lineMyImg_adj;%clearall%原始图像MyImg_origMyImg_orig=imread('pout.tif');MyImg_adj=MyImg_orig;%像素点行列[m,n]=size(MyImg_orig);%图像反色，得到MyImg_countMyImg_anti=256-MyImg_orig;%分段函数step_scale处理每个像素点的灰度值，得到MyImg_

2、2fori=1:mforj=1:nMyImg_line(i,j)=step_scale(MyImg_orig(i,j));endend%灰度均衡化处理，合并某些灰度%求累计概率,个数hNum,概率pLevel,累计分布sLevel,图像灰度级L[hNum,hLevel]=imhist(MyImg_orig);L=length(hLevel);pLevel=hNum/(m*n);sLevel(1)=pLevel(1);fori=2:LsLevel(i)=sLevel(i-1)+pLevel(i);end%计算新

3、的灰度级newLevel%用公式S=int[(L-1)*sLevel+0.5]得到新的灰度值new_hLevel=fix((L-1)*sLevel+0.5);%根据需要合并的灰度值，调整图像的灰度值foriL=1:Lifnew_hLevel(iL)~=hLevel(iL)MyImg_adj(find(MyImg_orig==hLevel(iL)))=new_hLevel(iL);endend%对比处理结果,显示对应的灰度直方图subplot(2,4,1);imshow(MyImg_orig);title('原

4、始图像');subplot(2,4,5);imhist(MyImg_orig);title('灰度直方图');subplot(2,4,2);imshow(MyImg_anti);title('反色图像');subplot(2,4,6);imhist(MyImg_anti);title('灰度直方图');subplot(2,4,3);imshow(MyImg_line);title('分段线性化');subplot(2,4,7);imhist(MyImg_line);title('灰度直方图');subplot

5、(2,4,4);imshow(MyImg_adj);title('灰度均衡化');subplot(2,4,8);imhist(MyImg_adj);title('灰度直方图');一、分段函数子程序step_scale.m(一)函数曲线(二)程序代码functionnew_a=step_scale(a)%分段函数：%0~64:x1.5%65~191:x1%192~255:x0.5ifa>=192new_a=fix((a+255)*0.5);elseifa>=65new_a=fix(64*1.5+(a-65)*1

6、);elsenew_a=fix(a*1.5);end二、结果