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

Question

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

我设法使用此处显示的MSER获得每个粒子的一般形状：

但我在删除背景时遇到问题。我尝试使用MATLAB的imfill，但它没有填满所有的粒子，因为有些粒子在边缘被切断了。有关如何摆脱背景或以其他方式找到粒子区域的任何提示？ 欢呼声。

编辑：这是imfill的样子：

编辑2：这是用于获取大纲的代码。我使用this作为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） -

输出图像 -