%此函数来提取特征线链码
%-------------------------
%作者：孙寅乐
%时间：2008.12.3
%-------------------------
%im为输入图像。注意，该图像必须是经过sobel,log,canny,prewitt等算子进行过边缘提取后的二值图像
%result为返回值，返回一个矩阵。矩阵的每行前两个数记录每一条特征线的起始坐标x1,y1值，终点坐标x2,y2。每行其余值记录每条特征线的链码。
%--------方向表------
% 1 2 3
% 8 4
% 7 6 5
%--------------------
function result=listcode(im)
im2=im;
[m n]=size(im);
count=1; %用来记录特征线的数量
%遍历原图，搜索链码起始点
for i=1:m
for j=1:n
if(im2(i,j)==1)
xy=searchtip(im2,i,j);
[m n]=size(im2);
x0=xy(1); %用points记录所有链码起始点的位置
y0=xy(2);
next=1; %用来判断是否结束对当前链码的搜索
number=1; %用来记录当前链码以搜索点的数量
list=[];
list(number)=0;
x=xy(1);
y=xy(2);
%每当搜索到一个起始点，搜索其周围点，追寻链码轨迹
while(next==1)
if(x<=0|x>m|y<=0|y>n) %如果链码当前点已到达边界，则停止搜索
next=0;
else
im2(x,y)=0; %把已搜索过的点值置零，以免重复搜索
%按照方向表中从1到8的顺序检索当前点周围是否存在非0点
if(im2(x,y+1)==1)
list(number)=4;
x=x;
y=y+1;
x1=x;
y1=y;
number=1+number;
if(x==x0&y==y0) %判断是否回到起始点，若是则当前特征线链码提取结束
next=0;
end
elseif(im2(x+1,y)==1)
list(number)=6;
x=x+1;
y=y;
x1=x;
y1=y;
number=1+number;
if(x==x0&y==y0)
next=0;
end
elseif(im2(x,y-1)==1)
list(number)=8;
x=x;
y=y-1;
x1=x;
y1=y;
number=number+1;
if(x==x0&y==y0)
next=0;
end
elseif(im2(x-1,y)==1)
list(number)=2;
x=x-1;
y=y;
x1=x;
y1=y;
number=1+number;
if(x==x0&y==y0)
next=0;
end
elseif(im2(x+1,y+1)==1)
list(number)=5;
x=x+1;
y=y+1;
x1=x;
y1=y;
number=1+number;
if(x==x0&y==y0)
next=0;
end
elseif(im2(x+1,y-1)==1)
list(number)=7;
x=x+1;
y=y-1;
x1=x;
y1=y;
number=1+number;
if(x==x0&y==y0)
next=0;
end
elseif(im2(x-1,y-1)==1)
list(number)=1;
x=x-1;
y=y-1;
x1=x;
y1=y;
number=1+number;
if(x==x0&y==y0)
next=0;
end
elseif(im2(x-1,y+1)==1)
list(number)=3;
x=x-1;
y=y+1;
x1=x;
y1=y;
number=1+number;
if(x==x0&y==y0)
next=0;
end
else
next=0;
end
end
end
list1=[x0 y0 x1 y1 list(:)'];
[p q]=size(list1);
list2(count,1:q)=list1;
count=count+1;
end
end
end
%figure:imshow(im);
%figure:imshow(im2);
result=list2;
end
