六 、图像复原.ppt

顺序统计滤波器 中值滤波器 用滤波器确定的该像素领域像素的 灰度中值来代替该像素的值 从小到大排列，取中间值 中值滤波器 中值滤波器应用非常广泛，因为 对许多随机噪声，他都有很好的 去噪作用，且在相同尺寸下比起 线性平滑滤波器引起的模糊较少 中值滤波器尤其对椒盐噪声非常有效 B = medfilt2(A,[m n]) B = medfilt2(A) 中值滤波在matlab中的实现 A:待处理的图像 [m,n]:滤波器的尺寸， 若省略，默认为[3 3] f1=medfilt2(g1); 最大值最小值滤波器 最大值滤波器 从小到大排列，取最大值 可以很好消除图像中“胡椒”噪声 1 1 1 1 1 1 1 1 1 0 1 0 1 1 1 0 1 0 最大值最小值滤波器 最小值滤波器 从小到大排列，取最小值 可以很好消除图像中“盐”噪声 1 1 1 1 1 1 1 1 1 在matlab中的实现 B = ordfilt2(A,order,domain) 作用：使用由domain确定的邻域 的一组排序元素中的第order个元素 来代替A中的每个元素 domain：一个由0和1组成的m*n的 矩阵，他指定了将在计算中使用的 邻域中像素的位置。计算中不使用 对应于domain中是0位置的像素 排序滤波器 排序滤波器 例：实现最小值滤波器 1 1 1 1 1 1 1 1 1 g=ordfilt2(f,1, domain); 0 1 0 1 1 1 0 1 0 domain=[0 1 0;1 1 1;0 1 0]; g=ordfilt2(f,1,domain); domain=ones(3,3); 最大值滤波器如何实现？ g=ordfilt2(f,m*n, domain); 最小值滤波器 最大 值滤波器 中点滤波器 中点滤波器是在滤波器涉及范围内 计算最大值和最小值之间的中点 找到最大值：212 最小值：198 对高斯和均匀分布 的噪声有较好的效果 Matlab中的实现 f1=ordfilt2(g,1,ones(m,n),’symmetric’); f2=ordfilt2(g,m*n,ones(m,n),’symmetric’); f=imlincomb（0.5,f1,0.5,f2); 修正后的阿尔法均值滤波器 在模版所确定的邻域中去掉d/2个最大值和d/2个最小值，剩下的像素的平均值形成的滤波器称为修正后的阿尔法均值滤波器 从小到大排列，去掉最大值最小值，取剩余值的平均值 求平均 d=2 修正后的阿尔法均值滤波器 修正后的阿尔法均值滤波器在包含多种噪声的情况下非常适用，例如高斯噪声和椒盐噪声混合的情况下 Matlab实现 function f=alphatrim(g,m,n,d) inclass =class(g); g=im2double(g); f=imfilter(g,ones(m,n),symmetric) for k=1:d/2 f=imsubstract(f,ordfilter2(g,k,ones(m,n),symmetric); end for k=(m*n-d/2+1):m*n f=imsubstract(f,ordfilter2(g,k,ones(m,n),symmetric); end f=f/(m*n-d); f=changclass(inclass,f); 被均匀噪声和椒盐噪声污染的图像 均值滤波后图像 几何均值滤波后图像 修正后阿尔法均值滤波 中值滤波后图像 自适应中值空间滤波器 令： zmin 表示Sxy中的最小亮度 zmax 表示Sxy中的最大亮度 zmed 表示Sxy中亮度中值 zxy 表示坐标(x,y)处的亮度 这个自适应中值滤波器算法工作在 两个层面，表示为level A、level B level A: 若zminzmedzmax则转向level B 否则增加窗口尺寸 若窗口尺寸=Smax，重复levelA 否则输出Zmed level B: 若zminzxyzmax,则输出zxy 否则输出zmed 在ftp上找到adpmedian.m，下载到 自己的工作路径，即可使用本函数 function f=adpmedian(g,Smax) Smax:滤波器的最大尺寸 中值滤波 5*5 7*7 自适应中值滤波器 f=imread(flowersh.tif); g3=imnoise(f,salt pepper,0.3); figure,imshow(g3); f3=medfilt2(g3); figur