MATLAB图像邻域运算,去除噪声,模板运算中值滤波,均值滤波

1、数字图像处理实验报告 #领域运算 领域运算列出常用的卷积模板基于3X3的模板，编写均值滤波的处理程序，处理含有加性高斯噪声和椒盐噪声的图像，观察处理结果编写中值滤波程序，处理相同的图像与均值滤波进行比较；改变模板尺寸观察处理结果编程实现利用一阶微分算子和二阶拉普拉斯算子进行图像锐化的程序对比不同的邻域运算结果，体会图像锐化与图像平滑的区别均值滤波处理含有椒盐噪声图像程序代码：模板运算f=imread(lenasp.bmp);f=double(f);row,col=size(f);r=1;fori=2:row-1forj=2:col-1g(i,j)=(f(i_1,j_1)+f(i_1,j)+f(

2、i_1,j+1)+f(i,j_1)+f(i,j)+f(i,j+1)+f(i+1,j-1)+f(i+1,j)+f(i+1,j+1)/9;g=uint8(g);endendsubplot(121);imshow(uint8(f);title(椒盐噪声图像);subplot(122);imshow(g);title(均值滤波处理后的图像)数字图像处理实验报告椒盐噪声图像均值滤波处理后的图像数字图像处理实验报告椒盐噪声图像均值滤波处理后的图像 领域运算 #领域运算：mean函数均值滤波f=imread(lenasp.bmp);f=double(f);row,col=size(f);g=size(f);

3、fori=2:row-1forj=2:col-1t=f(i-1:i+1,j-1:j+1);t=double(t);g(i,j)=mean(mean(t);g(i,j)=uint8(g(i,j);endendsubplot(121);imshow(uint8(f);title(椒盐噪声图像);subplot(122);imshow(uint8(g);title(均值滤波处理后的图像)数字图像处理实验报告均值滤波处理后的图像数字图像处理实验报告均值滤波处理后的图像 #领域运算 领域运算可变模板处理f=imread(lenagauss.bmp);f=double(f);row,col=size(f)

4、;r=1;g=size(f);fori=2:row-1forj=2:col-1s=0;s=double(s);form=-r:rforn=-r:rs=s+f(i+m,j+n);endendg(i,j)=s/power(r+2),2);endendsubplot(121);imshow(uint8(f);title(椒盐噪声图像);subplot(122);imshow(uint8(g);title(均值滤波处理后的图像)数字图像处理实验报告椒盐噪声图像均值滤波处理后的图像数字图像处理实验报告椒盐噪声图像均值滤波处理后的图像 领域运算 #领域运算均值滤波处理含有加性高斯噪声图像程序代码:模板运算

5、f=imread(lenagauss.bmp);f=double(f);row,col=size(f);r=1;fori=2:row-1forj=2:col-1g(i,j)=(f(i_1,j_1)+f(i_1,j)+f(i_1,j+1)+f(i,j_1)+f(i,j)+f(i,j+1)+f(i+1,j-1)+f(i+1,j)+f(i+1,j+1)/9;g=uint8(g);endendsubplot(121);imshow(uint8(f);title(加性高斯噪声图像);subplot(122);imshow(g);title(均值滤波处理后的图像)数字图像处理实验报告加性高斯噪声图像均值滤

6、波处理后的图像数字图像处理实验报告加性高斯噪声图像均值滤波处理后的图像 #领域运算 领域运算：mean函数均值滤波f=imread(lenagauss.bmp);f=double(f);row,col=size(f);r=1;g=size(f);fori=2:row-1forj=2:col-1t=f(i-1:i+1,j-1:j+1);t=double(t);g(i,j)=mean(mean(t);g(i,j)=uint8(g(i,j);endendsubplot(121);imshow(uint8(f);title(加性高斯噪声图像);subplot(122);imshow(uint8(g);

7、title(均值滤波处理后的图像)数字图像处理实验报告均值滤波处理后的图像数字图像处理实验报告均值滤波处理后的图像 领域运算 #领域运算:可变模板处理f=imread(lenagauss.bmp);f=double(f);row,col=size(f);r=1;g=size(f);fori=2:row-1forj=2:col-1s=0;s=double(s);form=-r:rforn=-r:rs=s+f(i+m,j+n);endendg(i,j)=s/power(r+2),endendsubplot(121);imshow(uint8(f);title(加性高斯噪声图像);subplot(1

8、22);imshow(uint8(g);数字图像处理实验报告数字图像处理实验报告 领域运算 领域运算title（均值滤波处理后的图像）椒盐噪声图像均值滤波处理后的图像中值滤波(椒盐噪声图像处理)f=imread(lenasp.bmp);f=double(f);row,col=size(f);g=size(f);r=1;fori=r+1:row-rforj=r+1:col-rt=f(i-r:i+r,j-r:j+r);id=0;fortm=1:5m=0;fork=1:9if(mt(k)m=t(k);id=k;endendt(id)=0;endg(i,j)=m;endendsubplot(121);

9、imshow(uint8(f);title(原图像);subplot(122);imshow(uint8(g);title(中值滤波处理后的图像)R=1时的图像原图像中值滤波处理后的图像R=5原图像中值滤波处理后的图像中值滤波(加性高斯噪声)f=imread(lenagauss.bmp);f=double(f);row,col=size(f);g=size(f);r=1;fori=r+1:row-rforj=r+1:col-rt=f(i-r:i+r,j-r:j+r);id=0;fortm=1:5m=0;fork=1:9if(mt(k)m=t(k);id=k;endendt(id)=0;endg

10、(i,j)=m;endendsubplot(121);imshow(uint8(f)title(原图像)subplot(122)imshow(uint8(g)title(中值滤波处理后的图像)原图像中值滤波处理后的图像R=5原图像中值滤波处理后的图像数字图像处理实验报告均值滤波处理后的图像数字图像处理实验报告 领域运算 领域运算不同领域运算结果比较r=2f=imread(lenasp.bmp);f=double(f);row,col=size(f);r=2;g=size(f);fori=r+1:row-rforj=r+1:col-rs=0;s=double(s);form=-r:rforn=-r:rs=s+f(i+m,j+n);end