function?Touying（g）
imshow（g）
f=rgb2gray（g）;
[m?n]=size（f）;
[row?col]=size（f）;
T=graythresh（f）
T=T*255
for?i=1:row
????for?j=1:col
???????if（f（ij）>T）
???????????f（ij）=255;
???????else?f（ij）=0;
???????end
????end
end
f=im2bw（f）;
a=zeros（row1）;
b=zeros（col1）;
for?i=1:row
????for?j=1:col
????????a（i）=a（i）+f（ij）;
????end
end
for?j=1:col
????for?i=1:row
????????b（j）=b（j）+f（ij）;
????end
end
for?i=1:row
????if?（a（i）>50）
????????break;
????end
end
p=i;
for?i=row:-1:1
????if?（a（i）>50）
????????break;
????end
end
q=i;
f=f（p:q1:end）;
figureimshow（f）figureplot（1:rowa）figureplot（1:colb）
c=zeros（16）;
k=1;
if?（k<17）
????for?i=（c（k）+1）:col
????????if?（mod（k2））
????????????if?（b（i）>3）
????????????????c（k）=i;
????????????????k=k+1;
????????????end
????????else?if?（b（i）<3）
???????????????c（k）=i;
???????????????k=k+1;
????????????end
???