分别利用canny、robert、和Laplace边缘提取算子提取lenna图像的边界并比较结果的空域与频域。

数字图像处理实验报告

姓名: 学号: 日期：

一、实验要求

（1）分别利用canny、robert、和Laplace边缘提取算子提取lenna图像的边界并比较结果的空域与频域。

（2）用matlab读取和显示

二、实验代码

I=imread('D:\1.jpg'); %读入图像文件

if isrgb(I)

I=rgb2gray(I);

end

%Laplace算法

[m,n]=size(I);

f=abs(I);

B=zeros(m,n);

for i=2:m-1

for j=2:n-1

x(i,j)=f(i+1,j)+f(i-1,j)-2*f(i,j);

y(i,j)=f(i,j+1)+f(i,j-1)-2*f(i,j);

B(i,j)=x(i,j)+y(i,j);

end

end

subplot(3,2,1),imshow(uint8(B));

title('Laplace算法');

A1=fft2(B);

C1=real(120*A1./abs(A1));

subplot(3,2,2),imshow(uint8(C1));

title('频谱');

%Roberts算法

C=zeros(m,n);

for i=2:m-1

for j=2:n-1

x(i,j)=f(i,j)-f(i+1,j+1);

y(i,j)=f(i+1,j)+-f(i,j+1);

C(i,j)=x(i,j)+y(i,j);

end

end

subplot(323),imshow(uint8(C));

title('Roberts算法');

A2=fft2(C);

C2=real(120*A2./abs(A2));

subplot(324),imshow(uint8(C2));

title('频谱');

%canny算法

D=zeros(m,n);

for i=2:m-1

for j=2:n-1

x(i,j)=(f(i+1,j)-f(i,j)+f(i+1,j+1)-f(i,j+1))/2;

y(i,j)=(f(i,j+1)-f(i,j)+f(i+1,j+1)-f(i+1,j))/2;

D(i,j)=x(i,j)+y(i,j);

end

end

subplot(325),imshow(uint8(D));

title('Roberts算法');

A3=fft2(D);

C3=real(120*A3./abs(A3));

subplot(326),imshow(uint8(C3));

title('频谱');

三、实验结果截图并做分析

分析：

小结：