計(jì)算機(jī)視覺實(shí)驗(yàn)報(bào)告Experiment3

1、Experiment 3: Edge DetectionClass:電子 1203 班 Student ID: 1210910322 Name:王影AimThe aim of this laboratory session is to learn to deal with image data by Matlab. By the end of this session, you should be able to perform image preprocessing of edge detection in spatial domain and frequency domain.Knowle

2、dge required in the Experiment.You are supposed to have learned the basic skills of using Matlab;.You need to review Matlab programming language and M-file format.iii. You should have studied edge detection methods.Experiment ContentsDemand: Please show the figure on the left and list the codes on t

3、he right respectively bellow each question.(請將運(yùn)行結(jié)果(圖片)和程序代碼貼在每題下方)i . Read “car.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 display the results in a Ma

4、tlab window.程序：clear;im=imread(car.jpg);I=rgb2gray(im);subplot(3,2,1);imshow(I);title(Gray image);Y,X=size(I);im_edge=zeros(Y,X);T=30;for k=2:Y-1for kk=2:X-1im_edge(k,kk)=abs(I(k+1,kk+1)-I(k,kk)+abs(I(k,kk+1)-I(k+1,kk);if (im_edge(k,kk)T) im_edge(k,kk)=1;elseim_edge(k,kk)=0;endendendsubplot(3,2,2);i

5、mshow(im_edge,);%EAEyYxOE6-Ap0255pA-IAUjtitle(Robert image);Y X=size(I);imedge=zeros(Y,X);for k=2:Y-1for kk=2:X-1imedge(k,kk)=abs(I(k-1,kk+1)-I(k-1,kkT)+abs(I(k,kk+1)-im(k,kkT) + abs(I(k+1,kk+1)-I(k+1,kk-1)+.abs(I(k+1,kk -1)-I(k-1,kk-1)+abs(I(k+1, kk)-I(k-1, kk)+abs(I(k+1,kk+1)-I(k-1,kk+1);end endsu

6、bplot(3,2,3);imshow(imedge,);title(Prewit image);Y X=size(I);im_edge=zeros(Y,X);for k=2:Y-1for kk=2:X-1im_edge(k,kk)=abs(I(k-1,kk+1)-I(k-1,kkT)+2火abs(I(k,kk+1)-I(k,kk-1) + abs(I(k+1,kk+1)-I(k+1,kk-1)+.abs(I(k+1,kk -1)-I(k-1,kk-1)+2火abs(I(k+1, kk)-I(k-1, kk)+abs(I(k+1,kk+1)-I(k-1,kk+1);end endsubplot

7、(3,2,4);imshow(im_edge,);title(Sobel image);圖像如下:ii. Read “car.jpg” file (to do this by imread function), convert the color image into grayscale image, then perform edge detection in frequency domain using Gaussian Highpass filter and display the result in a Matlab window.第二題程序：%頻域邊緣檢測，利用Gassian高通濾波

8、器進(jìn)行濾波，進(jìn)行邊緣檢測%頻域邊緣檢測，利用Gassian高通濾波器進(jìn)行濾波，進(jìn)行邊緣檢測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(I);Y,X=meshgrid(1:col,1:row);II=double(I).*(-1).A(X+Y);F=fft2(II);%subplot(2,2,2);%title(Fourier spectrum);%creat highpass filterD=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/(1+(D0/d)A(2*n);D(i,j)=1-exp(-dA2)/(2*(D0)A2);endend%filteringG=F.*D;%Invert the result and shift