clear,close all
I=imread('Image/mili.bmp');
imshow(I);
%%%估计背景图象
background=imopen(I,strel('disk',15));
figure,surf(double(background(1:8:end,1:8:end))),zlim([0 255]);
set(gca,'ydir','reverse');
%%从原始图片减去背景图象
I2=imsubtract(I,background);
figure,imshow(I2);
%%% 调节图象对比度
I3=imadjust(I2,stretchlim(I2),[0 1]);
figure,imshow(I3);
%%%使用阀值操作将图象转化为二进值图象
level=graythresh(I3);
bw=im2bw(I3,level);
figure,imshow(bw);
%%%对象个数
[labeled,numObjects]=bwlabel(bw,4);
%%%检查标记矩阵(imcrop)
grain=imcrop(labeled)
RGB_label=label2rgb(labeled