机器视觉与图像处理

1、机器视觉与图像处理课程大作业专业：自动化班级：1402学号：201423020230姓名：杨坤翔2017.06.05大作业说明1 .要求每位同学独立完成课程大作业；2 .鼓励组成课程小组，集体讨论研究，课程总结内说明小组成员；3 .允许借鉴网络、书籍上相关代码资源，但一定要切合题目内容；4 .根据题目要求，原理解释部分若需要公式，使用公式编辑器编辑;代码部分保证完整、可运行；结果部分黏贴原图；5 .作业鼓励将个人调试经验、学习心得等个性化内容总结。题目1:高斯滤波器与拉普拉斯滤波器1.1给出高斯滤波器与拉普拉斯滤波器数学表达式;(1)给出高斯滤波器数学表达式x2G(x)=e。21G(x，y):

2、2；ex2_y2一2二2(2)给出拉普拉斯滤波器数学表达式：江卅冬3+冬必(1)rV2/=f(x+1,>)4/*_1J)+/(1J+1)+/(”-l)-4/(x,j)(2)g(x7y)=fx.y)+<f(x.y)(3)1.2使用Matlab绘制出高斯滤波器与拉普拉斯滤波器图形表达;图形表达:alf=3;n=7;%定义模板大小n1=floor(n+1)/2);%确定中心fori=1:na(i)=exp(-(i-n1).A2)/(2*alfA2);forj=1:nb(i,j)=exp(-(i-n1)A2+(j-n1)A2)/(4*alf)/(4*pi*alf);endendsubplo

3、t(121),plot(a),title('一维高斯函数)subplot(122),surf(b),title('二维高斯函数,)25260-106.jpg');1.3分别使用高斯滤波器和拉普拉斯滤波器对下列图片进行卷积运算操作，是否有快速方法进行拉普拉斯滤波器卷积运算?高斯滤波:I=imread('C:Usersdell-5000PictureslovewallpaperH=rgb2gray(I);Img=double(H);alf=3;n=10;n1=floor(n+1)/2);fori=1:nforj=1:nb(i,j)=exp(-(i-n1)A2+(j-

4、n1)A2)/(4*alf)/(4*pi*alf);DMendendImg_n=uint8(conv2(Img,b,'same');K=uint8(imfilter(Img,b);Img_n2=uint8(imfilter(Img,b,'conv');J=(Img_n2)-Img_n;原图')卷积运算图')相关运算图)flag=mean(J(:)subplot(131),imshow(I);title('subplot(132),imshow(Img_n);title('subplot(133),imshow(K);title(&

5、#39;国便卷利为舞用机美堆舞阻(2)拉普拉斯滤波：h1=fspecial('laplacian');25260-106.jpg');A=imread('C:Usersdell-5000PictureslovewallpaperK=rgb2gray(A);B=imfilter(K,h1);subplot(131),imshow(A);title('原图')subplot(132),imshow(K);title("灰度图')subplot(133),imshow(B);title('相关运算图')忆餐:曲海题目2

6、：使用Canny算子边缘检测2.1列写出Canny算子检测边缘算法原理；(1)图象边缘检测必须满足两个条件：一能有效地抑制噪声；二必须尽量精确确定边缘的位置。(2)根据对信噪比与定位乘积进行测度，得到最优化逼近算子。这就是Canny边缘检测算子。(3)类似与Marr(LoG)边缘检测方法，也属于先平滑后求导数的方法。(4)Canny边缘检测算法：stepl:用高斯滤波器平滑图象；step2:用一阶偏导的有限差分来计算梯度的幅值和方向；step3:对梯度幅值进行非极大值抑制；step4:用双阈值算法检测和连接边缘。2.2使用Canny算法对下图进行边缘检测，并对比其他边缘检测算法，如Sobel,

7、Roberts等；1 .canny算子：I=imread('dazuoye02.jpg');I=rgb2gray(I);imshow(I);title('原图')BW1=edge(I,'canny');figure,imshow(BW1);title('matlabcanny检测')rratab旬可氾2 .Roberts算子：I=imread('dazuoye02.jpg');K=rgb2gray(I);BW1=1,0;0,-1;BW2=0,1;-1,0;J1=filter2(BW1,K);J2=filter2(B

8、W2,K);K1=double(J1);K2=double(J2);M=(abs(K1)+abs(K2);figure,imshow(uint8(M)title('matlabRobert检测')mallabRobert.则3.Sobel算子:I=imread('dazuoye02.jpg');K=rgb2gray(I);BW1=-1,-2,-1;0,0,0;1,2,1;BW2=-1,0,1;-2,0,2;-1,0,1;J1=filter2(BW1,K);J2=filter2(BW2,K);K1=double(J1);K2=double(J2);检测')

9、M=(abs(K1)+abs(K2);figure,imshow(uint8(M)title('matlabsobel结论：Roberts算子：边缘定位准，但是对噪声敏感。适用于边缘明显且噪声较少的图像分割。Roberts边缘检测算子是一种利用局部差分算子寻找边缘的算子,Robert算子图像处理后结果边缘不是很平滑。经分析，由于Robert算子通常会在图像边缘附近的区域内产生较宽的响应，故采用上述算子检测的边缘图像常需做细化处理，边缘定位的精度不是很高。Sobel算子：其主要用于边缘检测，在技术上它是以离散型的差分算子，用来运算图像亮度函数的梯度的近似值，缺点是Sobel算子并没有将图

10、像的主题与背景严格地区分开来，换言之就是Sobel算子并没有基于图像灰度进行处理，由于Sobel算子并没有严格地模拟人的视觉生理特征，所以提取的图像轮廓有时并不能令人满意Canny算子：该算子功能比前面几种都要好，但是它实现起来较为麻烦，Canny算子是一个具有滤波，增强，检测的多阶段的优化算子，在进行处理前，Canny算子先利用高斯平滑滤波器来平滑图像以除去噪声，Canny分割算法采用一阶偏导的有限差分来计算梯度幅值和方向，在处理过程中，Canny算子还将经过一个非极大值抑制的过程，最后Canny算子还采用两个阈值来连接边缘。题目3角点与斑点检测3.1 使用Harris算法检测下图中角点（h

11、arris算法实现）Harris算法：filename,pathname,=uigetfile('dazuoye03.jpg');ifischar(filename)returnendstr=pathnamefilename;pic=imread(str);iflength(size(pic)=3img=rgb2gray(pic);endm,n=size(img);tmp=zeros(m+2,n+2);tmp(2:m+1,2:n+1)=img;Ix=zeros(m+2,n+2);Iy=zeros(m+2,n+2);Ix(:,2:n+1)=tmp(:,3:n+2)-tmp(:,1

12、:n);Iy(2:m+1,:)=tmp(3:m+2,:)-tmp(1:m,:);Ix2=Ix(2:m+1,2:n+1).A2;Iy2=Iy(2:m+1,2:n+1).A2;Ixy=Ix(2:m+1,2:n+1).*Iy(2:m+1,2:n+1);h=fspecial('gaussian',77,2);Ix2=filter2(h,Ix2);Iy2=filter2(h,Iy2);Ixy=filter2(h,Ixy);R=zeros(m,n);fori=1:mforj=1:nM=Ix2(i,j)Ixy(i,j);Ixy(i,j)Iy2(i,j);R(i,j)=det(M)-0.06*

13、(trace(M)A2;endendRmax=max(max(R);loc=;tmp(2:m+1,2:n+1)=R;fori=2:m+1forj=2:n+1iftmp(i,j)>0.01*Rmaxsq=tmp(i-1:i+1,j-1:j+1);sq=reshape(sq,1,9);sq=sq(1:4),sq(6:9);iftmp(i,j)>sqloc=loc;j-1,i-1;endendendendX=loc(:,1);Y=loc(:,2);subplot(1,2,1);imshow(pic);subplot(1,2,2);imshow(pic);holdonplot(X,Y,&#

14、39;*');holdoff3.2 使用Log算子检测下图中斑点(Matlab:iog_Biob)(1)构造LoG_Blob函数:9-LoG_Blob,rri算+I1EfunctionpointsoGJ101)(1mgjnujL.'blo'bs)2 百砌能；提取LoG时早，3 Sing输入图像?？4 Mum-需要检娜斑点蜀目？7-Spoinl精则出的斑点?？iung二double(iuig(;j;：j1)>|ifnargin=l%如果输入卷数仅有一个<ljie)?8-nmE沈；将将检测斑立数设置为12。及nm=nwii_blobs;年nd*设定LoG叁数?s

15、igma_begin=2:si£ina_end=15.sijma_step=：l:(2)构造draw函数:fS-E:nnatlabb4ndraw.rrdraw,m筑+1 functiondraw(imgjptjstr.)2 1焉中能;在图像中绘制特征点管3 Mm输入图像?？4 5*M一一特征点坐标”5 Sistr图上显示的名都?76 fifur("IFajie''?str.):7 -imshow(ling);8 holdon;9 -axisoff:10 switchsue(ptj2)11 -case212 -s=2;13 -fori=1:size(ptf1)

16、14 -rtctangl?CPosition"ptSJ-sptKLl)-&2尉&2*51,"Curvature.15j0,0,'EdgeColorJbLineTidth,2);16 -end17 -case3IS-fori=1:sizefpt,1)IS-rectangle(JPositiojiajpt(i,2)_pt(i,3),pttijl)-pt(i,3)j.?202*pt(if3),2*pt(i31；JCurvaturet1,11；JEdgeGolorJ,21fvf/LineWidth*,2)：(3)算法实现：img=imread('da

17、zuoye04.jpg');imshow(img);pt=LoG_Blob(rgb2gray(img);draw(img,pt,'LOG')个人体会：在做题时，通过查找网上的资料才做出来，我觉得难点在于函数的书写上，函数写出来了，题目就做出来了。题目4特征点匹配4.1完成下列两图中的特征点检测与匹配(1)算法实现：>>i1=imread('dazuoye05.jpg');i2=imread('dazuoye06.jpg');i11=rgb2gray(i1);i22=rgb2gray(i2);imwrite(i11,'v

18、1.jpg','quality',80);imwrite(i22,'v2.jpg','quality',80);match('v1.jpg','v2.jpg');Findingkeypoints.879keypointsfound.Findingkeypoints.2163keypointsfound.Found337matches.>>il=inueadCjlazuoyeOS.jpg)；i2=iJitreadCdazuoye06.jpg")：ill=rgb?gray(id;i2rcb

19、2tray(i2)：imwrite(ilLJvl.jpgquality"】SO)imwrile(i22j-v2.jpg'Jquality",8Q)matchCvkjpg?jv2,jpEJ)；Fmdingk?ypoin,ts.B79keypointsfound.Fmdmgkeypoints.,.2163ktrptintsfound.Found337jnatches.ft>I(2)算法函数及结果:-E:m3V3bbinsiftDenrioV45ift.rri&rft.m+JL16 甯Lred£ts:Thanksf心上initialversiono

20、fthisprogrsjhtD./17 %J.J.Guerrero；,UniversidaddeZaragazatuiodifiedt1319|-functzaninage,des匚rzupt凸匚*,Iocs二sift1imafeFile)2021 >Loadiiwag«22 ijiafe=imread(imafeFile?:232425262723293。-31券33-34话-35%IfyouhavetheImageProcessingloalbox,youcanuncaiiunent*Lmsst电allvinputQfcoLotimag巳whichvillb?%ifisrg

21、h(image)%imae-=r£bZgray(imaea):%endrovs,cols=size(linage);ftConvertiivtaPGMimagefil*?.,readablebykeypointse量吧cutwf=fopenC-tmp.pgm',:iff=-1errorCCouldnolcreg.tefiletmp.pgm-'.);endUFtfC)SREtmotI»t>bir»sit>omoV*trri£»tcK.mmutcMmLO£11213LALS一16LTISim2021222324*1fKincXionn3nmal:obi<xmsiee