在MATLAB中删除背景和测量图像的特征

Question

我正在尝试测量此图像中显示的每个粒子的面积:

我设法使用这里显示的MSER获得每个粒子的一般形状:

但是我在删除背景时遇到了麻烦.我尝试使用MATLAB的imfill,但它没有填满所有的粒子,因为有些粒子在边缘被切断了.有关如何摆脱背景或以其他方式找到粒子区域的任何提示？干杯.

编辑:这是imfill看起来像:

编辑2:这是用于获取大纲的代码.我用这个用于MSER.

%Compute region seeds and elliptial frames.
%MinDiversity = how similar to its parent MSER the region is
%MaxVariation = stability of the region
%BrightOnDark is used as the  void is primarily dark. It also prevents dark
%patches in the void being detected.
[r,f] = vl_mser(I,'MinDiversity',0.7,...
                'MaxVariation',0.2,...
                'Delta',10,...
                'BrightOnDark',1,'DarkOnBright',0) ;

%Plot region frames, but not used right now
%f = vl_ertr(f) ;
%vl_plotframe(f) ;

%Plot MSERs
M = zeros(size(I)) ; %M = no of overlapping extremal regions

for x=r'
 s = vl_erfill(I,x) ;
 M(s) = M(s) + 1;
end

%Display region boundaries
figure(1) ;
clf ; imagesc(I) ; hold on ; axis equal off; colormap gray ;

%Create contour plot using the values
%0:max(M(:))+.5 is the no of contour levels. Only level 0 is needed so 
%[0 0] is used.

[c,h]=contour(M,[0 0]) ;;
set(h,'color','r','linewidth',1) ;

%Retrieve the image data from the contour image
f = getframe;
I2 = f.cdata;

%Convert the image into binary; the red outlines are while while the rest
%is black.
I2 = all(bsxfun(@eq,I2,reshape([255 0 0],[1 1 3])),3);
I2 = imcrop(I2,[20 1 395 343]);

imshow(~I2);

Answer 1

建议的解决方案/技巧和代码

看来你可以M在这里工作。您可以在此处使用的一种技巧是在图像的边界上填充零M，然后填充其孔。这将负责填充之前接触边界的斑点，因为现在由于零填充，不会有任何斑点接触边界。

因此，在有了之后M，您可以添加此代码 -

%// Get a binary version of M
M_bw = im2bw(M);

%// Pad zeros all across the grayscale image
padlen = 2; %// length of zeros padding
M_pad = padarray(M_bw,[padlen padlen],0);

%// Fill the holes
M_pad_filled = imfill(M_pad,'holes');

%// Get the background mask after the holes are gone
background_mask = ~M_pad_filled(padlen+1:end-padlen,padlen+1:end-padlen);

%// Overlay the background mask on the original image to show that you have
%// a working background mask for use
I(background_mask) = 0;
figure,imshow(I)

结果

输入图像 -

前景蒙版（这将是~background_mask）-

输出图像 -