How can I count the black particles without counting the small particles in the background?

Question

loris Ries 2023-4-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1945769-how-can-i-count-the-black-particles-without-counting-the-small-particles-in-the-background

编辑： Antoni Garcia-Herreros 2023-4-12

I want to automatically measure the particle size and number, but I can't because of the out-of-round shape and the contrast with the background.

rgb_im = imread('bild_1.jpg'); 
gray_im = adapthisteq(rgb2gray(rgb_im));
%  bw_im = imbinarize(gray_im,'adaptive','ForegroundPolarity','dark','Sensitivity',1); 
bw_im2 = ~imbinarize(gray_im,'adaptive','ForegroundPolarity','bright','Sensitivity',1); 
bw_im4 = imdilate(bw_im2,strel('disk',1)); 
BW1 = imfill(bw_im2,'holes'); 
BW2=255-BW1; 
figure, imshow(BW2)
imshowpair(BW1,BW2,'montage')

I have inverted the image and now I want to remove the small black dots from the background. However, I can't get that to work.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Antoni Garcia-Herreros 2023-4-12

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1945769-how-can-i-count-the-black-particles-without-counting-the-small-particles-in-the-background#answer_1214269

编辑：Antoni Garcia-Herreros 2023-4-12

在 MATLAB Online 中打开

Hello loris

Regionprops does exactly what you are looking for.

You can try something like this:

clear all
close all
rgb_im = imread('bild_1.jpg');  %Read RGB image
Error using imread>get_full_filename
File "bild_1.jpg" does not exist.

Error in imread (line 372)
        fullname = get_full_filename(filename);
gray_im=rgb2gray(rgb_im);   % RGB to grayscale
BW=~imbinarize(gray_im,'adaptive','ForegroundPolarity','bright','Sensitivity',1); %Binarize image
BWao = bwareaopen(BW,20); % Remove particles smaller than 20 pixels
BWaod = imdilate(BWao,strel('disk',2)); 
BWaodf=imfill(BWaod,'holes');
r=regionprops('table',BWaodf,'Centroid','Area','Circularity','EquivDiameter'); 
m=table2array(r); %Transform table to array
%Filter by cicularity
thr=0.6; % You may change this parameter
Points=m(m(:,4)>thr,:); % Array containing all the Particles 
% Points=[Area [pixels], Centroid X, Centroid Y, Circularity, Equiv Diameter ]
imshow(rgb_im)
hold on
viscircles([Points(:,2),Points(:,3)],Points(:,5)*0.5,'LineWidth',0.5);
Nparticles=size(Points,1); %Number of Particles
MeanD=mean(Points(:,5)); % Mean Diameter