华工电信数学实验轮廓线地提取.docx

题目：轮廓线的提取 学 院： 电子与信息学院 班 级： 2010级信工五班 学 生： 陈耿涛 序 号： 3号 一．实验内容 1．任意选取一幅灰度图像和一幅彩色图像，对算法中若干关键语句中进行调整，得出不同的实验结果，对这些结果进行分析，并与MATLAB自带的边缘检测做对比。 2．提出其它的轮廓线提取方法，与简单阈值法进行比较分析。 二．实验过程 1.（1）代码 ①灰度图轮廓线提取函数 function linegray(pix,n) %灰度图的轮廓线提取 figure; A=imread(pix); %读取指定的灰度图 T=linecalculate(A,n); %计算矩阵轮廓线 subplot(2,1,1); %将窗口分割为两行一列，下图显示于第一行 imshow(A); %显示原图像 title(灰度图原图); %图释 axis image; %保持图片显示比例 subplot(2,1,2); %下图显示于第二行 imshow(T); %显示提取轮廓线后的图片 title(提取轮廓线); %图释 axis image; %保持图片显示比例 ②计算矩阵轮廓线函数（利用sin函数进行非线性变换） function tlab=linecalculate(piclab,n) %计算矩阵轮廓线 [a,b]=size(piclab); %a,b分别等于矩阵的行数和列数 B=double(piclab); %将矩阵变为双精度矩阵 D=40*sin((pi/2)*(1/255)*B); %将矩阵进行非线性变换 tlab=piclab; %新建同等大小矩阵 for p=2:a-1 %计算矩阵轮廓线 for q=2:b-1 if abs((D(p,q)-D(p,q+1)))n||abs((D(p,q)-D(p,q-1)))n||abs((D(p,q)-D(p+1,q)))n||abs((D(p,q)-D(p-1,q)))n||abs((D(p,q)-D(p-1,q+1)))n||abs((D(p,q)-D(p+1,q-1)))n||abs((D(p,q)-D(p-1,q-1)))n||abs((D(p,q)-D(p+1,q+1)))n tlab(p,q)=0; %置轮廓线点为黑色 else tlab(p,q)=255; %置非轮廓线点为白色 end; end; end; 所得图像： 1、linegray(E:\pixx.jpg,1) 2、linegray(E:\pixx.jpg,3.5) ③计算矩阵轮廓线函数（利用log函数进行非线性变换） function tlab=linecalculate1(piclab,n) %计算矩阵轮廓线 [a,b]=size(piclab); %a,b分别等于矩阵的行数和列数 B=double(piclab); %将矩阵变为双精度矩阵 D=40*log(B./255+1)%将矩阵进行非线性变换 tlab=piclab; %新建同等大小矩阵 for p=2:a-1 %计算矩阵轮廓线 for q=2:b-1 if abs((D(p,q)-D(p,q+1)))n||abs((D(p,q)-D(p,q-1)))n||abs((D(p,q)-D(p+1,q)))n||abs((D(p,q)-D(p-1,q)))n||abs((D(p,q)-D(p-1,q+1)))n||abs((D(p,q)-D(p+1,q-1)))n||abs((D(p,q)-D(p-1,q-1)))n||abs((D(p,q)-D(p+1,q+1)))n tlab(p,q)=0; %置轮廓线点为黑色 else