数字图像(实验四.docVIP

福建农林大学计算机与信息学院 实验报告 课程名称： 数字图像处理（双语） （DigitalImageProcessing） 姓名： 系： 计算机科学与技术 专业： 计算机科学与技术 年级： 指导教师： 王长缨 职称： 副教授 2014～2015学年下学期 福建农林大学计算机与信息学院实验报告 系： 计算机科学与技术 专业：计算机科学与技术 年级：2012级 姓名： 学号： 实验室号田家炳510 计算机号 实验时间：2015年5月 日指导教师签字： 成绩： ExperimentIVNoisereductionwithadaptivefiltering(validated) （2periods） Objective Establishthedatadrivenideafromadaptivefilterinbothspatialandfrequencydomain. Methodology ApplyspatialadaptivefilterandWienerfiltertoanoisedphotorespectivelyandmakeacommentbycomparisonoftwomethods. Environment Hardware:GeneralPC Operatingsystem:WindowsXP Calculatingplatform:MATLAB6.5orhigherversion Experimentindetail: 1．DesignanadaptivelocalnoisereductionfilterandapplyittoanimagewithGaussiannoise.Comparetheperformanceoftheadaptivelocalnoisereductionfilterwitharithmeticmeanandgeometricmeanfilter. 代码: rt=imread(D:\tupian1.jpg); gray=rgb2gray(rt); subplot(2,3,1);imshow(rt); title(Originalpicture); subplot(2,3,2);imshow(gray); title(Graypicture); rtg=im2double(gray); rtg=imnoise(rtg,gaussian,0,0.005)%加入均值为0，方差为0.005的高斯噪声 subplot(2,3,3);imshow(rtg); title(Gaussianpicture); [a,b]=size(rtg); n=3; smax=7; nrt=zeros(a+(smax-1),b+(smax-1)); fori=((smax-1)/2+1):(a+(smax-1)/2) forj=((smax-1)/2+1):(b+(smax-1)/2) nrt(i,j)=rtg(i-(smax-1)/2,j-(smax-1)/2); end end figure; imshow(nrt); title(扩充后的图像); nrt2=zeros(a,b); fori=n+1:a+n forj=n+1:b+n form1=3:2 m2=(m1-1)/2; c=nrt2(i-m2:i+m2,j-m2:j+m2);%使用7*7的滤波器 Zmed=median(median(c)); Zmin=min(min(c)); Zmax=max(max(c)); A1=Zmed-Zmin; A2=Zmed-Zmax; if(A10A20) B1=nrt2(i,j)-Zmin; B2=nrt2(i,j)-Zmax; if(B10B20) nrt2(i,j)=nrt2(i,j); else nrt2(i,j)=Zmed; end continue; end end end end nrt3=im2uint8(nrt2); figure; imshow(nrt3); title(adaptivemedianfilter); ImplementWienerfilterwith“wiener2”functionofMatLabtoanimagewithGaussiannoiseandcomparetheperformancewithadaptivelocalnoisereductionfilter. Experimentalresult Datainwhateverchart,imageandotherformisrequ