湖南大学脑与认知科学概论实验报告3详解

1、课 程 名 称脑与认知科学概论实验课时4实 验 项 目 名 称 和 编 号CS0502503同组者姓 名实 验 目 的1进一步理解边沿检测的基本原理。2 掌握对图像边沿检测的基本方法。3学习利用Matlab图像工具箱对图像进行边沿检测。实 验 环 境硬件：联想 Y50计算机软件：WINDOWS操作系统应用软件：MATLAB2014A实 验 内 容 和 原 理实验要求对边缘检测的要求：使用Matlab图像处理工具箱中的不冋方法对下图的边缘进行提取；注意观察不冋操作对垂直方向、水平方向、斜方向的提取效果有何区别；注意观察提取后的边界是否连续，若 不连续可采用什么方法使其连续。实验步骤1打开计算机，

2、启动 MATLA龍序；2. 调入数字图像1.bmp,利用MATLABi像工具箱中已有函数进行编程以实现图像的边缘检测；重点理 解edge函数的使用方法，参数的含义，具体如下：canny算子（sy3_1.m）:不同取值范围的效果。重点理解canny的意思。sobel 算子（sy3_2.m）: log拉普拉斯算法（sy3_3.m） imfilter 函数的使用（sy3_4.m）3. 调入数字图像 2.jpg，检测其边缘（sy3_5.m）4. 显示原图和处理过的图像。5 记录和整理实验报告实验原理调用matlab下的help edge函数可以看到有关 edge的使用说明：BW= edge(l) 采用

3、1作为它的输入，并返回一个与1相同大小的二值化图像BW在函数检测到边缘的地方为1，其他地方为0。BW = edge(l,sobel)自动选择阈值用 Sobel算子进行边缘检测。BW = edge(l,sobel,thresh)根据所指定的敏感度阈值thresh，用Sobel算子进行边缘检测，它忽略了所有小于阈值的边缘。当thresh为空时，自动选择阈值。BW = edge(I,sobel,thresh,direction)根据所指定的敏感度阈值thresh，在所指定的方向direction 上，用Sobel算子进行边缘检测。Direction 可取的字符串值为horizontal( 水平方向)

4、、vertical(垂直方向)或both(两个方向)。BW,thresh = edge(l,sobel,.)返回阈值以上标红的语句是 edge函数的完整参数实现，而我们在做实验的过程中只需要填写一个或几个即可完 成边缘检测。下面，考虑几个算子的使用方法：Canny算子，Canny的目标是找到一个最优的边缘检测算法，最优边缘检测的含义是： 好的检测-算法能够尽可能多地标识出图像中的实际边缘。好的定位-标识出的边缘要尽可能与实际图像中的实际边缘尽可能接近。最小响应-图像中的边缘只能标识一次，并且可能存在的图像噪声不应标识为边缘。为了满足这些要求 Canny使用了变分法，这是一种寻找满足特定功能的函

5、数的方法。最优检测使用四 个指数函数项的和表示，但是它非常近似于高斯函数的一阶导数。边氽检测是一种比较新的辿缰检则算子，具有很好的边塢监测性能，在圉樣处理中得到了越来拯广泛的应用。它依 据圏像边煨检测最优;隹则设计。込泌燼楡测算法：（1）首先用烈高斯遥逋複扳魅行卷稅以消除噪靑（2）利用导教空艸到图慢灰度地沿着两个右向的偏导教（0斗*q ）,弄求出梯摩的大小,冏=阳心利用（2）的结果i+H曲第度的方向&二Af e ,ilu 一且知道了瞒的方向，就可以把边逢的棵廣方向大致分为四种；水平、姿直、4艰方向、1节播方向沪通过梯 度的方向，就可以找到这个橡章梯度方向的邻接像養。（5） 遍历圏橡，若某个像耒

6、的页度值与其梯度方向上前后两个愎薰的换度值相比不是最大的，那么这个橡臺值養 拘山即不呈边煖v（B）使用累计宜方图计算两个阈值，大于裔阈值的一定是边曝，小于低阈值的一定不是边壕，介于之间的，看返个像素的邻接像素中肓没有超过葛阈值的迩家橡素，如果有的话那么它就是边壕了，若则它就不是谨康。可以看到，在调用时，thresh阈值可以自己设定，在接下来的实验中，我会自己尝试设定thresh值,进行边缘检测。Sobel算子，是计算机视觉领域的一种重要处理方法。主要用于获得数字图像的一阶梯度，常见的应用和物理意义是边缘检测。在技术上，它是一个离散的一阶差分算子，用来计算图像亮度函数的一阶梯 度之近似值。在图像

7、的任何一点使用此算子，将会产生该点对应的梯度矢量或是其法矢量。SotaM子是憑波算子的形式来提取迤“X，丫方向各用一个檯板，两个榄板组合起来构感1个梯度茸子Q K方向摸丽寸垂直边爆嶷响最夫，丫方向模板对水平边塚融向最犬b0 -1-1LoG算子也就是Laplace of Gaussian function （高斯拉普拉斯函数）。常用于数字图像的边缘提取和 二值化。LoG算子源于D.Marr计算视觉理论中提出的边缘提取思想，即首先对原始图像进行最佳平滑 处理，最大程度地抑制噪声，再对平滑后的图像求取边缘。由于噪声点(灰度与周围点相差很大的像素点)对边缘检测有一定的影响，所以效果更好的边缘检测器是L

8、oG算子，也就是Laplacian-Gauss 算子。它把的 Gauss平滑滤波器和Laplacian锐化滤波器结 合了起来，先平滑掉噪声，再进行边缘检测，所以效果会更好。常用的LoG算子是5X 5的模板。拉普谊斯高斯算子是一种二阶导数算子，将在逖杲处产生一个陡耀的窸交罠。前面介錯的几种梯度法具有方向性苹能对 各种走向的逖泉鄒具育相同的增强議果、fBLapliC19nM子是各向同性的，能对任何未向的霁绒和线築进行锐北，无方向 性亠这是拉晉也斯聲子区别于其也算法的最犬优点。对一个连续,它在位畫(/& fQ1 f/ 询啦普扭斯算子走文如下：/ = 4- + 4- mch*在鬥橡边媒檎测中，为了运箕

9、方便，函埶的拉晋竝斯高斯箕子也是借助檯板来实现的其桓扳有一个基本要求：檯板中心的至 数为正，其余相邻系数为负，所有轴的和应该为寒*-1-10-1-1-1Imfilter 函数：用法：B = imfilter(A,H)B = imfilter(A,H,optio n1,optio n2,.)或写作 g = imfilter(f, w, filtering_mode, boundary_options, size_options)其中，f 为输入图像，w为滤波掩模，g为滤波后图像。filtering_mode用于指定在滤波过程中是使用“相关”还是“卷积”。boundary options 用于处理

10、边界充零问题，边界的大小由滤波器的大小确定。1) canny算子边缘检测： close all;clear all;%author I1=imread(1.bmp); %读取图像%考虑canny的参数，本程序通过改变检测阈值进行调试 11= rgb2gray(l1);实 验 步 骤方法关键代码subplot(2,3,1),imshow(l1),title( I2=edge(I1,ca nn y); subplot(2,3,2),imshow(I2),title(ca nny I3=edge(I1,ca nn y,0.002,0.2,1);subplot(2,3,3),imshow(I3),ti

11、tle( I4=edge(I1,ca nn y,0,0.1,1);subplot(2,3,4),imshow(l4),title( I5=edge(I1,ca nn y,0,0.99,1);subplot(2,3,5),imshow(l5),title( I6=edge(I1,ca nn y,0,0.01,1);subplot(2,3,6),imshow(l6),title(2)sobel算子边缘检测close all;clear all;%authorI1=imread(1.bmp); %读取图像I1= rgb2gray(I1); subplot(2,3,1),imshow(I1),titl

12、e( I2=edge(I1,sobel);%考虑sobel的参数，有阈值和检测方向 subplot(2,3,2),imshow(I2),title(sobel I4=edge(I1,sobel,horizo ntal);subplot(2,3,3),imshow(I4),title( I5=edge(I1,sobel,vertical);subplot(2,3,4),imshow(I5),title( I6=edge(I1,sobel,0.01);subplot(2,3,5),imshow(I6),title( I7=edge(I1,sobel,0);subplot(2,3,6),imshow

13、(I7),title(原始图象1);边缘检测方法-默认阈值);阈值范围0.002,0.2);阈值范围0,0.1);阈值范围0,0.99);阈值范围0,0.01);原始图象1);下面更改两种方式进行实验 边缘检测方法-默认方向);仅水平方向);仅竖直方向);阈值0.1的sobel算子);阈值0的sobel算子);3)log算子边缘检测close all;clear all;%authorI1=imread(1.bmp); %读取图像11= rgb2gray(l1);subplot(2,3,1),imshow(l1),title(原始图象 1);I2=edge(l1,log);%考虑log的参数，

14、有阈值和检测方向，下面更改两种方式进行实验subplot(2,3,2),imshow(I2),title(log边缘检测方法-默认方向);I4=edge(I1,log,horizo ntal);subplot(2,3,3),imshow(I4),title(仅水平方向);I5=edge(I1,log,vertical);subplot(2,3,4),imshow(l5),title(仅竖直方向);I6=edge(l1,log,0.01);subplot(2,3,5),imshow(I6),title(阈值 0.01 的 log 算子);I7=edge(I1,log,0);subplot(2,3

15、,6),imshow(I7),title(阈值 0 的 log 算子);4) matlab在edge函数中给出了很多算子，下面进行对2.jpg进行测试%厂图展示6个算子边缘检测close all;clear all;%authorI1=imread(2.jpg); %读取图像I1= rgb2gray(I1);I2=edge(I1,sobel); %Sobel算子subplot(2,3,1),imshow(l2),title(Sobel边缘检测方法);I3=edge(I1,prewitt); %prewitt算子subplot(2,3,2),imshow(I3),title(prewitt边缘检

16、测方法);I4=edge(I1,roberts); %roberts算子subplot(2,3,3),imshow(I4),title(roberts边缘检测方法);I5=edge(I1,log); %Laplacia n-Gaussia n算子subplot(2,3,4),imshow(l5),title(Laplacian-Gaussian边缘检测方法);I6=edge(I1,zerocross); %zerecross算子subplot(2,3,5),imshow(I6),title(过零点边缘检测方法);I7=edge(I1,ca nny); %canny算子subplot(2,3,6

17、),imshow(I7),title(canny边缘检测方法);5) 不使用算子对2.jpg的检测clear all;a=imread(2.jpg);A=rgb2gray(a);%authorh1=1,1,1;0,0,0;-1,-1,-1;%水平方向h2=0,1,1;-1,0,1;-1,-1,0; %45度万向h3=1 0 -1;1 0 -1;1 0 -1;%竖直方向h4=0,-1,-1;1,0,-1;1,1,0;%135j仁imfilter(A,h1);%对数组滤波，j2=imfilter(A,h2);j3=imfilter(A,h3); j4=imfilter(A,h4); x=uint8

18、(j1); %转化成 uint8 型y=ui nt8(j2);z=ui nt8(j3); m=u in t8(j4); D=imadd(x,y);D仁imadd(D,z);%进行叠加D2=imadd(D1,m);D3=imadd(D1,m); subplot(321),imshow(A);title(origi n) subplot(322),imshow(j1);title( subplot(323),imshow(j2);title(45 subplot(324),imshow(j3);title( subplot(325),imshow(j4);title(135 subplot(326

19、),imshow(D2);title(度方向F同水平方向) 度方向)竖直方向) 度方向) 叠加后的图像)1)canny算子多种检测结果:测 试 记 录 分 析 结 论匣射国眾I聞建圖fwzqzIfl 范因ID.D IrtiilSLirooi词皆克曲仙2)sobel算子多种检测结果:噩1E：列曰Sii 1 的m -1?3 log算子多种检测结果:個静:平角问硯（1.州旳网干SttwKiBoWr4）对2.jpg的多种算子检测5）不使用算子对2.jpg的检测结果unyri常$厘方向墜方向7分析与总结：1） canny算子的目标是找到一个最优的边缘检测算法，在对一个陌生图像进行边缘检测时，canny算 子是最好的选择。2） 关于canny算子的参数设置，阈值范围设定，A,B下限决定了检测的细节是否清晰，上限决定了是否会将大量的次要边缘包括在内，考虑canny算子的参数sigma即标准差，我对sigma进行了第二次设定后（选定值为 0.5 ）进行了又一次检