图像处理与模式识别作业4图像与滤波.docVIP

图像处理与模式识别作业4 ----- 图像 滤波 姓名： 班级：自动化 学院：电信学院 联系方式： 日期：2012/4/6【运行环境：vs平台opencv】 分别用高斯滤波器和中值滤波器去平滑测试图像test1和2，模板大小分别是3x3 ， 5x5 ，7x7 ;1). matlab高斯滤波器程序代码： si=3; ss=(si+1)/2; sig2=1.5*1.5; mins=exp(-2.0*ss*ss/(2*sig2));%令最小项系数为1 he=0;%整数系数和 guass1=ones(si); for i=1:si for j=1:si temp=(i-ss)*(i-ss)+(j-ss)*(j-ss); guass1(i,j)=round(exp(-temp/(2*sig2))/mins); he=he+guass1(i,j); end end I=imread(F:/图像处理作业/第四次作业/第四次作业/test2.tif); [h,w,x]=size(I); J=I; for i=ss:h+1-ss for j=ss:w+1-ss tempp=double(I(i-ss+1:i+ss-1,j-ss+1:j+ss-1)); z=tempp.*guass1; J(i,j)=round(sum(sum(z))/he); end end figure; subplot(121); imshow(I); title(原图像) subplot(122); imshow(J); title(高斯滤波后) 执行结果： 3x3: 5x5: 以3x3高斯滤波器为例， test1.pgm: test2.tif: 2). c中值滤波器程序代码： si=3; ss=(si+1)/2; I=imread(F:/图像处理作业/第四次作业/第四次作业/test2.tif); [h,w,x]=size(I); J=I; for i=ss:h+1-ss for j=ss:w+1-ss tempp=[]; for k=1:si tempp=[tempp,(I(i-ss+1:i+ss-1,j-ss+k:j-ss+k))]; end tempp=sort(tempp); J(i,j)=tempp(round((si*si+1)/2)); end end figure; subplot(121); imshow(I); title(原图像) subplot(122); imshow(J); title(中值滤波后) 执行结果： test1.pgm（3x3中值滤波器）: test1.pgm（5x5中值滤波器）: test1.pgm（7x7中值滤波器）: test2.tif（3x3中值滤波器）: test2.tif（5x5中值滤波器）: test2.tif（7x7中值滤波器）: 2.利用高通滤波器滤波测试图像test3,4：包括unsharp masking, Sobel edge detector, and Laplace edge detection ;  1).unsharp masking Matlab程序代码： si=3; ss=(si+1)/2; I=imread(F:/图像处理作业/第四次作业/第四次作业/test4 copy.bmp); [h,w,x]=size(I); J=I; for i=ss:h+1-ss for j=ss:w+1-ss tempp=double(I(i-ss+1:i+ss-1,j-ss+1:j+ss-1)); z=tempp.*[1 1 1;1 -9 1;1 1 1]; J(i,j)=round(0.1*sum(sum(z)))+I(i,j); end end figure; subplot(121); imshow(I); title(原图像) subplot(122); imshow(J); title(unsharp masking) test3_corrupt.pgm: test4.tif: 2). Sobel edge detector(阈值取100): Matlab 程序代码： si=3; ss=(si+1)/2; I=imread(F:/图像处理作业/第四次作业/第四次作业/test