话不多说,直接上代码
for num_graph=0:5:99 %每隔5张图片画一个边界图,这样看得更清楚
s1='C:\Users\86135\Desktop\competition\2001\2001aphoto\'
s2=int2str(num_graph)
s3='.bmp'
s=strcat(s1,s2,s3)
graphN=imread(s) %读入位图
edgeN=edge(graphN,'log') %求出边界点
for i =1:512
for j=1:512
if edgeN(i,j)==1
plot3(j,512-i,num_graph,'.')
hold on
end
end
end
end
axis equal
box on
tic %记录运行时间
mat_maxR=[]
mat_meanC0=[]
for l=0:99
graphN=[]
edgeN=[]
x=[]
y=[]
co=[]
s1='C:\Users\86135\Desktop\competition\2001\2001aphoto\'
s2=int2str(l)
s3='.bmp'
s=strcat(s1,s2,s3)
graphN=imread(s)
edgeN=edge(graphN,'log')
[x,y]=find(edgeN==1)
leftLine=min(y)+28
topLine=min(x)+