资源描述:
《matlab数字图像处理报告》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、昆明理工大学信息工程与自动化学院学生实验报告(2011—2012学年第1学期)课程名称:图形图像处理开课实验室:信自楼4442012年5月30日年级、专业、班计科091学号200910405116姓名李雅诚成绩实验项目名称图像分割指导教师王剑教师评语教师签名:年月日一、实验目的1.拉普拉斯算子、马尔算子(LOG)、Canny算子图像处理;2.边缘提取和跟踪3.简单图像的区域分割(峰谷法,类间最大距离法,最佳熵法)二.实验内容即步骤:1.马尔算子:I=imread('g:图片moon.tif');imshow(I);bw1=edge(I,
2、'roberts');%--------’roberts算子‘图像分割bw2=edge(I,'prewitt');%--------’prewitt算子‘图像分割bw3=edge(I,'sobel');%--------'sobel算子'图像分割bw4=edge(I,'log');%--------'log算子'(拉普拉斯-高斯)图像分割bw5=edge(I,'canny');%--------'canny算子'图像分割figure,subplot(2,3,1),imshow(bw1,[]);subplot(2,3,2),imshow(bw
3、2,[]);subplot(2,3,3),imshow(bw3,[]);subplot(2,3,4),imshow(bw4,[]);subplot(2,3,5),imshow(bw5,[]);h1=[01;-10];%--------’roberts算子‘矩阵h2=[-101;-101;-101];%--------’prewitt算子‘矩阵h3=[-10-1;-202;-101];%--------‘sobel算子'矩阵h4=[0-10;-141;0-10];%--------'log算子'矩阵J1=imfilter(I,h1);%---
4、-----采用’roberts算子‘矩阵化锐J2=imfilter(I,h2);%--------采用’prewitt算子‘矩阵化锐J3=imfilter(I,h3);%--------采用‘sobel算子'矩阵化锐J4=imfilter(I,h4);%---------采用'log算子'矩阵化锐figure,subplot(1,2,1),imshow(bw1,[]);title('roberts算子的分割');subplot(1,2,2),imshow(J1);title('roberts算子的锐化');figure,subplot(1,
5、2,1),imshow(bw2,[]);title('prewitt算子的分割');subplot(1,2,2),imshow(J2);title('prewitt算子的锐化');figure,subplot(1,2,1),imshow(bw3,[]);title('sobel算子的分割');subplot(1,2,2),imshow(J3);title('sobel算子的锐化');figure,subplot(1,2,1),imshow(bw4,[]);title('log算子的分割');subplot(1,2,2),imshow(J4)
6、;title('log算子的锐化');figure,imshow(bw5,[]);原始图像roberts算子、prewitt算子、sobel算子、log算、canny算子最终处理图像Canny算子代码:%log算子和canny边缘检测(拉普拉斯-高斯、canny)I=imread('g:图片testpat1.png');subplot(1,3,1);imshow(I);title('原始图像');I1=edge(I,'log');subplot(1,3,2);imshow(I1);title('log算子分割结果');I2=edge(I
7、,'canny');subplot(1,3,3);imshow(I2);title('canny算子分割结果');运行结果:2.边缘提取和跟踪边缘跟踪1:%functioncolordrawout(pix,n)%彩色图片轮廓线提取函数A=imread('h:图片peppers.png');%读取指定彩色图片n=20;B=A(:,:,1);%红色强度值矩阵C=A(:,:,2);%绿色强度值矩阵D=A(:,:,3);%蓝色强度值矩阵fori=1:3%依次从三个矩阵中提取轮廓线ifi==1%从红色矩阵提取E=B;elseifi==2%从绿色矩
8、阵提取E=C;elseE=D;%从蓝色矩阵提取end;end;H=double(E);%将选择的矩阵变为双精度矩阵%F=40*sin(1/255*H);%进行非线性变换F=H;[