clc;
clear;
target='D:\Users\LIULIFEN310\Desktop\testresult\bw\'
image='D:\Users\LIULIFEN310\Desktop\testresult\test\'
%image=''
dst='D:\Users\LIULIFEN310\Desktop\testresult\analysis\';
list=dir(fullfile(image));
fileNum=size(list,1)-2;
for k=3:fileNum
filename=list(k).name;
inputImage=[image,filename];
pos=strfind(filename,'.');
relateFile=[filename(1:pos),'png']
relateTarget=[target,relateFile];
imgTarget=imread(relateTarget);
imgTarget=double(imgTarget)*255;
%find couter
% bw=im2bw(imgTarget);
% bw2=bwperim(bw,8);
Image=imread(inputImage);
B = bwboundaries(imgTarget,8);%get boundry x,y
x=B{1,1}(:,1);%cell數(shù)據(jù)采用{}獲取元素叨叙，array用[]獲取元素
y=B{1,1}(:,2);%x為縱坐標(biāo)匈挖，y為橫坐標(biāo)
l=length(x);
output=figure('Visible', 'off');%只寫出figure，不顯示
imshow(Image);
hold on;
for i=1:l
? ? plot(y(i),x(i),'b.');%
end
outputfile=[dst,filename(1:pos),'png'];
saveas(output,outputfile);
close(output);
end
