图像处理(第六周备课).pptVIP

3.2 图像的邻域操作 邻域操作 指在图像处理时，输入要处理像素的某邻域内各个像素值，输出处理后像素的新值，处理结果与邻域内各像素值有关。 【例3-10】读入荷花图像，并化为灰度图像矩阵。 A0=imread(D:\hehua.jpg); A1=rgb2gray(A0) A2=double(A0); A3=floor((A0(:,:,1)+A0(:,:,2)+A0(:,:,3))/3) subplot(1,3,1) imshow(A0) subplot(1,3,2) imshow(A1) subplot(1,3,3) imshow(A3) Matlab邻域操作函数 nlfilter函数： 该函数被调用的一般形式为：B =nlfilter(A,[M N],FUN)，其中A是被处理的图像，[M N]表示滑动块的高为M、宽为N。FUN是要对滑动覆盖区域进行操作的函数名，该函数可以是各种滤波算子，可以是任意定义的矩阵。 nlfilter函数 nlfilter函数 【例3-13】使用nlfilter函数对图像进行均值平滑邻域操作。 除了排序滤波是非线性滤波，还可以定义一个非线性函数来进行滤波操作。 随着（x,y）的改变， A表示的算子在图像上滑动。 在数学意义上完成了二阶偏倒数的计算； 在图像处理意义上完成了基于拉普拉期算子的滤波操作。 算子：A= fspecial(laplacian,ALPHA)产生的滤波器 非线性滤波也是基于邻域操作的，但是不象线性滤波那样，滤波器与对应元素相乘然后相加。 非线性滤波不再用线性组合的方法得到新的元素，使用的是统计的方法或非线性组合的方法。 fun = inline(median(x(:))); B = nlfilter(A1,[3 3],fun); median(x(:))，这种排序取中值操作就是一种非线性滤波操作。 3、非线性空间滤波 MATLAB还提供了函数medfilt2()专门进行二维中值滤波。 【例3-19】使用函数medfilt2进行二维中值滤波。 设计如下程序： A= imread(‘D:\005.jpg); A1=rgb2gray(A); B=medfilt2(A1); subplot(1,2,1) imshow(A) subplot(1,2,2) imshow(B) 从上图可以看出，中值滤波器具有去掉噪音的功能，，但是同时也把窄的边线也滤掉了。中值滤波的效果比较好是因为噪音的颜色比较突出，多数不是一个小区域内各个颜色值的中间值。 从【例3-19】结果可以看出，去除噪音的效果不是太好，下面可以使用一个通用的统计排序滤波器函数ordfilt2,该函数的调用方式为: B=ordfilt2(A,order,domain) A:为要处理的图像矩阵。 domain:为一个矩阵，大小为滑动操作时邻域的大小。 【例3-20】使用函数ordfilt2进行排序滤波。 A= imread(D:\005.jpg); A1=rgb2gray(A); order=1; domain=[1 0 0 1; 0 0 0 0;0 0 0 0;1 0 0 1]; B=ordfilt2(A1,order,domain); subplot(1,2,1);imshow(A) subplot(1,2,2);imshow(B) 上例程序：B=ordfilt2(A1,order,domain)中： domain：进行滑动操作时，当算子滑动到一个新的位置，就把domain中有1的位置对应的（图像上）元素的值取出排序。 order：为排序后选取第几个作为输出。 例如：domain为[4 4]大小，四个角上的值为1，其他位置为0。当滑动到一个新的位置，就把[4 4]大小的被覆盖区域中的4个顶点的值取出排序。order=1，所以取4个顶点的最大值。 【例3-21】使用colfilt函数进行非线性滤波操作。 function v=gmean(A) [m n]=size(A) v=prod(A,1).^(1/(m+n)) B=imread(d:\005.jpg) C=rgb2gray(B) D=double(C) E=padarray(D,[3 3],replicate) F=colfilt(E,[3 3],sliding,@gmean) subplot(1,2,1) imshow(C) subplot(1,2,2) imshow(F) * * 图 像 处 理 邻域操作使用的方法 使用算子模版逐行（或逐列）地图像上滑动，滑动时对模版经过区域进行运算，把计算的结果作为区域中心像素的新值，所以有时也称为滑动邻域操作。 荷花图像是方形的，高与宽都为480。A1(:,:,1) 、A1(:,:,2) 与