实验四图像增强

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

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

2、all;I{1}=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}

3、{Index},[])endend执行结果：-word.zl--图1幂次变换增强结果实验结果分析：由实验结果可知，当r<1时，黑色区域被扩展，变的清晰；当r>1时，黑色区域被压缩，变的几乎不可见。2.直方图规定化处理程序代码：-word.zl--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_

4、image_cumulation=cumsum(Hist_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

5、(2*N-1-2*Index);Hist{2}=Hist{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=a

6、bs(Hist_image_cumulation(k)-Hist_cumulation{m});[Temp1,Project{m}(k)]=min(Temp);end%2.2变换后直方图fork=1:N-word.zl--Temp=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(['变换后的直方图',n

7、um2str(m)]);%2.3结果图Step=256/N;forK=1:NIndex=find(I>=Step*(k-1)&I