实验四图像增强

《实验四图像增强》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、.信息工程学院实验报告成绩：指导老师(签名)：课程名称：数字图像处理实验项目名称：实验四图像增强实验时间：2016.11.08班级：姓名：学号：一、实验目的1.了解图像增强的目的及意义，加深对图像增强的感性认识，巩固所学理论知识。2.掌握图像空域增强算法的基本原理。3.掌握图像空域增强的实际应用及MATLAB实现。4.掌握频域滤波的概念及方法。5.熟练掌握频域空间的各类滤波器。6.掌握怎样利用傅立叶变换进行频域滤波。7.掌握图像频域增强增强的实际应用及MATLAB实现。二、实验步骤及结果分析1.基于幂次变换的图像增强程序代码：..clearall;closeall;I{1}

2、=double(imread('fig534b.tif'));I{1}=I{1}/255;figure,subplot(2,4,1);imshow(I{1},[]);holdonI{2}=double(imread('room.tif'));I{2}=I{2}/255;subplot(2,4,5);imshow(I{2},[]);holdonform=1:2Index=0;forlemta=[0.55]Index=Index+1;F{m}{Index}=I{m}.^lemta;subplot(2,4,(m-1)*4+Index+1),imshow(F{m}{Index},[

3、])endend执行结果：..图1幂次变换增强结果实验结果分析：由实验结果可知，当r<1时，黑色区域被扩展，变的清晰；当r>1时，黑色区域被压缩，变的几乎不可见。2.直方图规定化处理程序代码：..clearallclccloseall%0.读图像I=double(imread('lena.tiff'));subplot(2,4,1);imshow(I,[]);title('原图')N=32;Hist_image=hist(I(:),N);Hist_image=Hist_image/sum(Hist_image);Hist_image_cumulation=cumsum(Hi

4、st_image);%累计直方图subplot(245);stem(0:N-1,Hist_image);title('原直方图');%1.设计目标直方图Index=0:N-1;%正态分布直方图Hist{1}=exp(-(Index-N/2).^2/N);Hist{1}=Hist{1}/sum(Hist{1});Hist_cumulation{1}=cumsum(Hist{1});subplot(242);stem([0:N-1],Hist{1});title('规定化直方图1');%倒三角形状直方图Hist{2}=abs(2*N-1-2*Index);Hist{2}=His

5、t{2}/sum(Hist{2});Hist_cumulation{2}=cumsum(Hist{2});subplot(246);stem(0:N-1,Hist{2});title('规定化直方图2');%2.规定化处理Project{1}=zeros(N);Project{2}=zeros(N);Hist_result{1}=zeros(N);Hist_result{2}=zeros(N);form=1:2Image=I;%SML处理(SML,SingleMappingLaw单映射规则fork=1:NTemp=abs(Hist_image_cumulation(k)-H

6、ist_cumulation{m});..[Temp1,Project{m}(k)]=min(Temp);end%2.2变换后直方图fork=1:NTemp=find(Project{m}==k);ifisempty(Temp)Hist_result{m}(k)=0;elseHist_result{m}(k)=sum(Hist_image(Temp));endendsubplot(2,4,(m-1)*4+3);stem(0:N-1,Hist_result{m});title(['变换后的直方图',num2str(m)]);%2.3结果图Step=256/N;forK=1:N

7、Index=find(I>=Step*(k-1)&I