图像处理第二组梁金玉 刘光旭 刘薇 赵野 刘阳.docVIP

6.4 层叠滤波技术的Matlab实现 在Matlab中利用现有函数并不能在多大程度上提高层叠运算的速度。但它为Matlab提供了一个很好的思路，那就是可以利用Matlab提供的二值图像操作去实现相应的灰度图像操作。实际上很多图像滤波方法都可以转化为相应的层叠滤波方法。 下面给出了利用层叠方法实现中值滤波和灰度形态滤波的两个例子。 6.4.1 层叠中值滤波的Matlab实现 在上一章，我们讨论了各种中值滤波及其Matlab实现，下面的源代码是在Matlab中利用层叠方法实现层叠滤波。 % 读取原图像 a0=imread(testpat1.png); figure(1) imshow(a0); %如图6-5所示 %加入高斯噪声 a1=imnoise(a0,gaussian,0,0.5); figure(2) imshow(a1) %如图6-6所示 %直接进行中值滤波 a2=medfilt2(a1); figure(3) imshow(a2) %如图6-7所示 c0=zeros(256); c1=zeros(256); c2=zeros(256); %按照阈值为1进行阈值分解, 进行层叠形态滤波 for j=1:255 b=a1; for i=1:65536 if b(i)=j b(i)=1; else b(i)=0; end end b=im2double(b); c0=c0+medfilt2(b); end figure(4) imshow(c0); %如图6-8所示 %按照阈值为2进行阈值分解, 进行层叠形态滤波 for j=1:2:255 b=a1; for i=1:65536 if b(i)=j b(i)=1; else b(i)=0; end end b=im2double(b); c1=c1+medfilt2(b); end figure(5) imshow(c1); %如图6-9所示 %按照阈值为8进行阈值分解, 进行层叠形态滤波 for j=1:8:255 b=a1; for i=1:65536 if b(i)=j b(i)=1; else b(i)=0; end end b=im2double(b); c2=c2+medfilt2(b); end figure(6) imshow(c2); %如图6-10所示 图6-5 图6-6 图6-7 图6-8 图6-9 图6-10 6.4.2 层叠秩排序滤波的Matlab实现 在下一章里我们将会讲到在Matlab中提供了三个形态滤波函数，实质上它们都属于秩排序统计滤波器，它们都是针对二值图像的，如前面所属利用层叠算法可以方便实现灰度图像的形态滤波。下面是在Matlab利用层叠方法实现秩排序统计滤波器的例子。 %读取源图像 a0=imread(rice.png); figure(1) imshow(a0); %如图6-11所示 %加入高斯噪声 a1=imnoise(a0,salt pepper, 0.02); figure(2) imshow(a1); %如图6-12所示 %直接进行秩排序膨胀操作 a2=ordfilt2(a1,9, ones(3)); figure(3) imshow(a2); %如图6-13所示 c0=zeros(256); c1=zeros(256); c2=zeros(256); c3=zeros(256); c4=zeros(256); c5=zeros(256); c6=zeros(256); c7=zeros(256); %按照阈值为1进行阈值分解, 进行层叠膨胀 for j=1:255 b=a1; for i=1:65536 if b(i)=j b(i)=1; else b(i)=0; end