20陈细云作业3.doc
光电图像处理作业专业： 应用物理 任课老师： 作业名称:图像的均值处理完成时间：1、 一、作业要求：1,对一个图像用一个3*3均值板模糊处理,比较模糊后图像的直方图和原图像的直方图2 选用一幅自拍图像作为原图像,分别添加高斯噪声、椒盐噪声,生成高斯噪声污染图像和椒盐噪声污染图像；再分别使用3×3均值滤波器、中点滤波器与自适应中值滤波器对污染图像进行滤波,得到三幅复原图像,计算复原图像与原图像的结构相似度.2、 作业步骤与运行结果：作业一function hist_testf=imread<'cameraman.tif'>h=imhist<f>h1=h<1:10:256>horz=1:10:256;% bar<horz,h1>% stem<horz,h1,'fill'>% figure,plot<h>figure,imshow<f> figure,imhist<f>axis<-5 255 0 7000>set<gca,'xtick',-5:50:255>set<gca,'ytick',0:500:7000>xlabel<'pixels gray value','fontsize',12>ylabel<'pixels number','fontsize',12>title<'直方图'>% figure,imshow<f>% figure,imhist<f>% axis<0 255 0 20000>% set<gca,'xtick',0:50:255>% set<gca,'ytick',0:2000:20000>% ylim<'auto'>g=histeq<f,256>figure,imshow<g>figure,imhist<g>axis<0 255 0 7000>set<gca,'xtick',0:50:255>set<gca,'ytick',0:500:7000>% ylim<'auto'>hnorm=imhist<f>./numel<f>cdf=cumsum<hnorm>x=linspace<0,1,256>figure,plot<x,cdf>axis<0 1 0 1>set<gca,'xtick',0:.2:1>set<gca,'ytick',0:.2:1>xlabel<'Input intensity values','fontsize',9>ylabel<'Output intensity values','fontsize',9>text<0.18,0.5,'Transformation function','fontsize',9>s=0;end实验2function experience4_ssimf=imread<'chenxiyun.jpg'>f=mat2gray<im2double<f>>f=f<:,:,2>figure<1>,imshow<f,>m=0;sigma=0.025;g_gaussian=imnoise<f,'gaussian',m,sigma> noise=imnoise<zeros<size<f>>,'gaussian',m,sigma>-g_sp=imnoise<f,'salt & pepper',0.01>figure<2>,subplot<,2,1>,imshow<g_gaussian,>text<100,950,'高斯噪声污染图像'>subplot<1,2,2>,imshow<g_sp,>text<100,950,'椒盐噪声污染图像'> figure<3>,space_filter<f,g_gaussian>figure<4>,space_filter<f,g_sp>x=7;sita=100; PSF=fspecial<'motion',x,sita>g_b=imfilter<g_gaussian,PSF,'circular'> fr1=deconvwnr<g_b,PSF>Sn=abs<fft2<noise>>.2;nA=sum<Sn<:>>./prod<size<noise>>Sf=abs<fft2<f>>.2;fA=sum<Sf<:>>./prod<size<f>>R=nA/fA;fr2=deconvwnr<g_b,PSF,R> NCORR=fftshift<real<ifft2<Sn>>>ICORR=fftshift<real<ifft2<Sf>>>fr3=deconvwnr<g_b,PSF,NCORR,ICORR> M,N=size<f>noisepower=M*N*<sigma-m.2>fr4=deconvreg<g_b,PSF,noisepower>fr5=deconvreg<g_b,PSF,0.1*noisepower,1e-7 1e7>figure<5>,subplot<2,3,1>,imshow<g_b,>text<100,950,'退化图像'>subplot<2,3,2>,imshow<fr1,>text<100,950,'直接逆滤波图像'>subplot<2,3,3>,imshow<fr2,>text<100,950,'参数维纳滤波图像'>subplot<2,3,4>,imshow<fr3,>text<100,950,'相关系数维纳滤波图像'>subplot<2,3,5>,imshow<fr4,>text<100,950,'单参数正则滤波图像'>subplot<2,3,6>,imshow<fr5,>text<100,950,'2参数正则滤波图像'>endfunction space_filter<f,g_noise>g_noise=mat2gray<g_noise>fp_amean=spfilt<g_noise,'amean',3,3>fp_amean=mat2gray<im2double<fp_amean>>fp_gmean=spfilt<g_noise,'gmean',3,3>fp_gmean=mat2gray<im2double<fp_gmean>>fp_hmean=spfilt<g_noise,'hmean',3,3>fp_hmean=mat2gray<im2double<fp_hmean>>fp_chmean_plus=spfilt<g_noise,'chmean',3,3,1.5>fp_chmean_plus=mat2gray<im2double<fp_chmean_plus>>fp_chmean_minus=spfilt<g_noise,'chmean',3,3,-1.5>fp_chmean_minus=mat2gray<im2double<fp_chmean_minus>>fp_median=spfilt<g_noise,'median',3,3>fp_median=mat2gray<im2double<fp_median>>fp_max=spfilt<g_noise,'max',3,3>fp_max=mat2gray<im2double<fp_max>>fp_min=spfilt<g_noise,'min',3,3>fp_min=mat2gray<im2double<fp_min>>fp_midpoint=spfilt<g_noise,'midpoint',3,3>fp_midpoint=mat2gray<im2double<fp_midpoint>>fp_atrimmed=spfilt<g_noise,'atrimmed',3,3,2>fp_atrimmed=mat2gray<im2double<fp_atrimmed>>fp_adp=adpmedian<g_noise,7>fp_adp=mat2gray<im2double<fp_adp>>ssim_amean=spatial_structure_similarity<f,fp_amean>ssim_gmean=spatial_structure_similarity<f,fp_gmean>ssim_hmean=spatial_structure_similarity<f,fp_hmean>ssim_chmean_plus=spatial_structure_similarity<f,fp_chmean_plus>ssim_chmean_minus=spatial_structure_similarity<f,fp_chmean_minus>ssim_median=spatial_structure_similarity<f,fp_median>ssim_max=spatial_structure_similarity<f,fp_max>ssim_min=spatial_structure_similarity<f,fp_min>ssim_midpoint=spatial_structure_similarity<f,fp_midpoint>m_atrimmed=spatial_structure_similarity<f,fp_atrimmed>m_adp=spatial_structure_similarity<f,fp_adp>subplot<3,4,1>,imshow<g_noise>text<100,550,'噪声污染图像'>subplot<3,4,2>,imshow<fp_amean>text<100,550,'算术均值滤波图像'>subplot<3,4,3>,imshow<fp_gmean>text<100,550,'几何均值滤波图像'>subplot<3,4,4>,imshow<fp_hmean>text<100,550,'几何均值滤波图像'>subplot<3,4,5>,imshow<fp_chmean_plus>subplot<3,4,6>,imshow<fp_chmean_minus>subplot<3,4,7>,imshow<fp_median>subplot<3,4,8>,imshow<fp_max>subplot<3,4,9>,imshow<fp_min>subplot<3,4,10>,imshow<fp_midpoint>subplot<3,4,11>,imshow<fp_atrimmed>subplot<3,4,12>,imshow<fp_adp>ssim_amean = 0.5115ssim_amean =0.8702ssim_gmean = 0.2842ssim_gmean = 0.5775ssim_hmean = 0.3167ssim_hmean = 0.6679ssim_chmean_plus = 0.4485ssim_chmean_plus = 0.8150ssim_chmean_minus = 0.2412ssim_chmean_minus = 0.6864ssim_median = 0.3703ssim_median =0.9786ssim_max = 0.4136ssim_max = 0.6974ssim_min = 0.1779ssim_min =0.6590ssim_midpoint = 0.4092ssim_midpoint =0.6194m_atrimmed =0.4969m_atrimmed =0.9162m_adp = 0.1595m_adp =0.9895教师批阅 年 月 日
