求出叶子边界周长 求出8方向链码 计算周长CI.doc

求出叶子边界周长 求出8方向链码 计算周长CIR clc; clear all; close all; I=imread(D:\我的文档\MATLAB\yezi.jpg); I=im2bw(I); [m,n]=size(I); %计算图像的大小 figure,imshow(I); R=1; %二值图像的背景像素 I0=zeros(m+2,n+2)+R; %原图像+边框，避免边界查找出错 I0(2:(m+1),2:(n+1))=I; x0=0; y0=0; fsp=0; %是否找到起始点 %%%%%%%%%%%%寻找目标区域的一个边界点%%%%%%%%%%%%%%%%%%%%%%%% for i=1:m for j=1:n if (I0(i,j)~=R) x0=i; %记录起始点的纵坐标i y0=j; %记录起始点的横坐标j fsp=1; break; end end if fsp==1 break; end end row=x0; col=y0; %%%%%%%%%%%%%%%%%%%%按某一方向查找边界坐标并且记录%%%%%%%%%%%%%%%%%%%%%% A=[]; %存储边界点坐标的矩阵 B=zeros(m,n); %边界跟踪的图像矩阵 ended=0; direction=4; %由起始点的左边开始查找 while(ended==0) found_next=0; b=[row-1,col-1]; %消除边框，与原图保持一致 A=[A;b]; B(row-1,col-1)=255; while(found_next==0) switch mod(direction,8) case 0 if (I0(row, col+1)~=R) row=row; %记录当前坐标 col=col+1; direction=5; %下一点由方向5开始查找 found_next=1; %找到边界点,准备跳出循环 end; case 1 if (I0(row-1, col+1)~=R) row=row-1; %记录当前坐标 col=col+1; direction=6; %下一点由方向6开始查找 found_next=1; %找到边界点 end; case 2 if (I0(row-1, col)~=R) row=row-1; %记录当前坐标 col=col; direction=7; %下一点由方向7开始查找 found_next=1; %找到边界点