陈细云作业.docVIP

光电图像处理作业 专业： 应用物理 班级： 11级应用物理（2）班 任课老师： 余宏生 姓名： 陈细云 学号： 201141220220 作业名称: 图像的均值处理 完成时间：2014-5-17 一、作业要求： 1,对一个图像用一个3*3均值板模糊处理，比较模糊后图像的直方图和原图像的直方图 2 选用一幅自拍图像作为原图像，分别添加高斯噪声、椒盐噪声，生成高斯噪声污染图像和椒盐噪声污染图像；再分别使用3×3均值滤波器、中点滤波器及自适应中值滤波器对污染图像进行滤波，得到三幅复原图像，计算复原图像与原图像的结构相似度。 作业步骤及运行结果： 作业一 function hist_test f=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 实验2 function experience4_ssim f=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);