实验三_图像边缘检测new.docVIP

广西工学院信计系 上机实验报告 课　　程：数字图像处理 第 1 页 / 共 4 页 专业班级：电科082　　　　　　　 实验日期： 2011 年 6 月 10日 姓　　名：梁芳梅　　　　　　　　　　　　 学号：　200800902044 实验三 图像的边缘检测 一、实验目的 1、了解边缘检测的基本原理及方法； 2、掌握用MATLAB语言进行图象的边缘检测的方法。 二、实验要求 1、分别编写用soble、prewitt、LOG（高斯拉普拉斯算子）进行边缘检测的程序，对lena.bmp进行边缘检测，观察边缘检测的结果，给出你的分析； 2、选择其他不同类型（边缘较多与较少）的图像，再进行边缘检测，分析结果 三. 程序源代码： 1、（1） %边缘检测log.m f=imread(cameraman.tif); %读入灰度图像； [width,height]=size(f); %计算灰度图像的尺寸； f=double(f); %将整形转变为双精度形； %h=[-2 -4 -4 -4 -2; % -4 0 8 0 -4; % -4 8 24 8 -4; % -4 0 8 0 -4; % -2 -4 -4 -4 -2]; % h=[ 0 0 -1 0 0; 0 -1 -2 -1 0; -1 -2 16 -2 -1; 0 -1 -2 -1 0; 0 0 -1 0 0]; % g=zeros(width,height); %建立一个width*height的全零矩阵； for i=3:1:width-2 %变量i从1到width-4开始循环； for j=3:1:height-2 %变量j从1到height-4开始循环； sum=0; for m=1:1:5 for n=1:1:5 sum=sum+f(i+m-3,j+n-3)*h(m,n); end end g(i,j)=sum; end %变量j循环结束； end %变量i循环结束； for i=1:1:height %变量i从1到height开始循环； g(1,i)=0; % g(2,i)=0; % g(width,i)=0; % g(width-1,i)=0; % end %变量i循环结束； for i=1:1:width %变量i从1到height开始循环； g(i,1)=0; % g(i,2)=0; % g(i,height)=0; % g(i,height-1)=0; % end %变量i循环结束； subplot(1,2,1),imshow(uint8(f)); %将当前窗口分成1*2的区域，原图像显示在第一区； title(原图像); subplot(1,2,2),imshow(uint8(g),[]); %将当前窗口分成1*2的区域，滤波后的图像显示在第二区； title(边缘检测结果); （2） %边缘检测prewitt.m f=imread(cameraman.tif); %读入灰度图像； [width,height]=size(f); %计算灰度图像的尺寸； f=double(f); %将整形转变为双精度形； %prewitt %h1=[-1 0 1;-1 0 1; -1 0 1]; % %h2=[-1 -