计算机视觉实验报告Experiment3

1、Experiment 3 : Edge DetectionClass: 电子 1203班Student ID: 1210910322 Name: 王影I . AimThe aim of this laboratory sessi on is to lear n to deal with image data by Matlab. By the end of this sessi on, you should be able to perform image preprocess ing of edge detect ion in spatial doma in and freque ncy d

2、oma in.n . Knowledge required in the Experimenti .You are supposed to have learned the basic skills of using Matlab;ii .You need to review Matlab programming Ianguage and M-file format.iii . You should have studied edge detection methods.川.Experiment ContentsDema nd: Please show the figure on the le

3、ft and list the codes on the right respectively bellow eachquestion.(请将运行结果(图片)和程序代码贴在每题下方)i .Read “ar.jpg” file (to do this by imread function), convert the color image into grayscale image, and then perform edge detection using Roterts, Prewitt, Sobel operator separately in spatial domain and disp

4、lay the results in a Matlab win dow.程序：clear;im=imread( car.jpg);I=rgb2gray(im);subplot(3,2,1);imshow(I);title(Gray image);Y,X=size(l);im_edge=zeros(Y,X);T=30;for k=2:Y-1for kk=2:X-1im_edge(k,kk)=abs(l(k+1,kk+1)-l(k,kk)+abs(l(k,kk+1)-l(k+1,kk);if (im_edge(k,kk)T)im_edge(k,kk)=1;整理范本e-seimledge(k-kkH

5、qendendendsubp_2.(3-2-2)-imshowumledgeew 口 e 7 e y7Y X 77 77 7 卞 702255 卞 7 77 0 7 口 7 出s-e(-Robell image)-Y xHsNe (三imedgeHzeros(Y-xfor k2Y1for kk2x1imedge(k-kkHabs(-?1 -kk+1)(k?kk1)+abs(-(Kkk+1)m(Kkk1) +abs (一 (k+1 - kk+1 )(k+1 - kka)+ -abs (一 (k+1-kka(k?kka)+abs (一 (k+1- kk)(k?kk)+abs (一 (k+1-kk+

6、1)(k?kk+1)二 endendsubp_2.(3-2-3)-imshow(imedge-ss-e(-prewifimage)-Y xHsNe (三i mledge zeros (Y-xxfor k2Y1for kk2x1imledge(k-kkllabs(-(k?kk+v-?1 -kka)+2*abs(-(k-kk+1T-(Kkk1) +abs (- (k+1 - kk+1 )(k+1 - kk1)+ : abs(-(k+1-kk1T-(k?kka)+2*abs(-(k+1- kk)?1- kk)+abs(-(k+1-kk+1)?1-kk+1)_endendsubp-2.(3-2-4im

7、show(imledge-ss-e(-sobe-image)-mJlaMii. Read “ar.jpg” file (to do this by imread function), con vert the color image into grayscale image, then perform edge detect ion in freque ncy doma in using Gaussia n Highpass filter and display the result in a Matlab win dow.第二题程序：%频域边缘检测，利用Gassian高通滤波器进行滤波，进行

8、边缘检测%频域边缘检测，利用Gassian高通滤波器进行滤波，进行边缘检测clear; im=imread(car.jpg);I=rgb2gray(im); subplot(1,2,1); imshow(I); title(gray image);%shifting image (multiply the image by (_1)x+y)row,col=size(l);Y,X=meshgrid(1:col,1:row);ll=double(l).*(-1)A(X+Y);F=fft2(II);%subplot(2,2,2);%title(Fourier spectrum);%creat highpass filter%傅立叶变换中心D=zeros(row,col);u0=floor(row/2);v0=floor(col/2);D0=40;%截止频率n=2;for i=1:rowfor j=1:cold=(i-u0)A2+(j-v0)A2)A0.5;% D(i,j)=1心+(D0/d)A(2* n);D(i,j)=1-exp(-dA2)/(2*(D0)A2);endend%filteringG=F.*D;%Invert the result and s