《20陈细云作业3》由会员分享，可在线阅读，更多相关《20陈细云作业3（5页珍藏版）》请在金锄头文库上搜索。
1、光电图像处理作业专业： 应用物理 班级： 11级应用物理（2）班 任课老师： 余宏生 姓名： 陈细云 学号： 0 作业名称: 图像的均值处理 完成时间：2014-5-17 1、 一、作业要求：1,对一个图像用一个3*3均值板模糊处理，比较模糊后图像的直方图和原图像的直方图2 选用一幅自拍图像作为原图像，分别添加高斯噪声、椒盐噪声，生成高斯噪声污染图像和椒盐噪声污染图像；再分别使用33均值滤波器、中点滤波器及自适应中值滤波器对污染图像进行滤波，得到三幅复原图像，计算复原图像与原图像的结构相似度。2、 作业步骤及运行结果：作业一function hist_testf=imread(camerama
2、n.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);%
3、 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);fi
4、gure,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);fig
5、ure(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
6、,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(if
7、ft2(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,直接逆滤波图像)
8、;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,
9、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_m
10、inus=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);
11、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=s
12、patial_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_si
13、milarity(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_am
14、ean);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.2842
