Freeman链码及其某些变化

%Freeman链码及其某些变化
f=imread('noisy_circular_stroke.tif');
subplot(3,2,1),imshow(f),title('原图')

h=fspecial('average',9);%掩模
g=imfilter(f,h,'replicate');%9*9平均掩模处理
subplot(3,2,2),imshow(g),title('用9*9平均掩模平滑后的图像')

g=im2bw(g,0.5)%阈值处理
subplot(3,2,3),imshow(g),title('经阈值处理后的图像')

B=boundaries(g)%边界坐标
d=cellfun('length',B);%B的长度
[max_d,k]=max(d);%最长边界
b=B{1};%最长边界所在区域
[M,N]=size(g);%g的大小
g=bound2im(b,M,N,min(b(:,1)),min(b(:,2)));%获取边界图像
subplot(3,2,4),imshow(g),title('二值图像边界')

[s,su]=bsubsamp(b,50);%二次取样
g2=bound2im(s,M,N,min(s(:,1)),min(s(:,2)));%二值图像
subplot(3,2,5),imshow(g2),title('二次取样后的边界')

cn=connectpoly(s(:,1),s(:,2));%重新连接
g2=bound2im(cn,M,N,min(cn(:,1)),min(cn(:,2)));%二值图像
subplot(3,2,6),imshow(g2),title('连接二次取样后中的点后的效果')

c=fchcode(su);
c.x0y0
c.fcc
c.mm
c.diff
c.diffmm