clear?all;
for?kk=0:9
?????p1=ones（1616）;
????m=strcat（int2str（kk）‘.bmp‘）
????x=imread（m‘bmp‘）;
????bw=im2bw（x0.5）;
????[ij]=find（bw==0）;
????imin=min（i）;
????jmin=min（j）;
???imax=max（i）;
???jmax=max（j）;
???bw1=bw（imin:imaxjmin:jmax）;
???rate=16/max（size（bw1））;
???bw1=imresize（bw1rate）;
???[ij]=size（bw1）
???i1=round（（16-i）/2）;
???j1=round（（16-j）/2）;
???p1（i1+1:i1+ij1+1:j1+j）=bw1;
???p1=-1.*p1+ones（1616）;???
???for?m=0:15
???????p（m*16+1:（m+1）