粒度计算从给定的图像

37视图(30天)

显示旧的评论

Vikrant普拉塔普 2020年5月10日

1
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/524362-particle-size-calculation-from-given-image

回答: 图像分析 2020年5月10日

我要计算每个粒子的大小从输入图像使用图像处理技术。

答案(2)

蒂亚戈恩里克戈麦斯Lobato 2020年5月10日

0
链接

直接链接到这个答案

//www.tatmou.com/matlabcentral/answers/524362-particle-size-calculation-from-given-image answer_431546

的大小每一个粒子是一个任务自动化非常困难,你总会得到一些错误由于不同光线条件下,粒子分布,粒子重叠等。更容易处理和现实的一个问题是计算岩石的大小分布, 这将给你尺寸是如何统计的分布在图像。我拍了一些灵感从图像分析教程 ,使一个代码,你可以作为一个开始微调。这样你至少可以得到粒度分布的估计:

                             我= imread (“image.jpeg”);
                            
                             Ibinary = rgb2gray (I) > 150;你的问题尝试不同的百分比值
                            
                             Ibinary = imfill (Ibinary,“黑洞”);
                            
                             C = 0.008;%转换因子m /像素。这是一个猜测,您应该能够找到真正的一个
                            
                             [labeledImage, numberOfBlobs] = bwlabel (Ibinary);
                            
                             blobMeasurements = regionprops (labeledImage,“重心”,“EquivDiameter”);%我相信当量直径是一个很好的方法来测量尺寸,
                            
                             %你还可以检查不同的指标
                            
                             %得到当量直径
                            
                             EquivDiameter = [blobMeasurements.EquivDiameter];
                            
                             ValidDia =找到(EquivDiameter > 10);%的一些经验阈值
                            
                             blobMeasurements = blobMeasurements (ValidDia);
                            
                             EquivDiameter = EquivDiameter (ValidDia);
                            
                             %的概率估计的大小
                            
                             图,柱状图(EquivDiameter * C,“归一化”,“概率”),包含(的大小(m -未校准)),ylabel (“概率”)
                            
                             %显示一些粒子
                            
                             图中,imshow (Ibinary);持有在
                            
                             %显示区域形象
                            
                             为idx = 1:2:长度(ValidDia)%遍历所有气泡。
                            
                             质心= [blobMeasurements (idx) .Centroid (1) blobMeasurements (idx) .Centroid (2)];
                            
                             DiamSize = [“直径= '。,num2str ((EquivDiameter (idx) * C)));
                            
                             文本(质心(1),重心(2),DiamSize,“颜色”,“c”);
                            
                             viscircles(质心EquivDiameter (idx) / 2);
                            
                             结束

3评论
显示隐藏特殊照顾年长的评论

Vikrant普拉塔普 2020年5月10日

谢谢你的回答的蒂亚戈。

粒度分布也将我的目的服务。

我用下面的代码。产生的气泡是不正确的,不包括所有的粒子。

请建议更改的代码。

你也可以建议一种方法来衡量粒子而不是最长的一面直径。

                                  清晰的所有
                                 
                                  关闭所有
                                 
                                  rgb = imread (“image.jpg”);
                                 
                                  imshow (rgb)
                                 
                                  I = rgb2gray (rgb);
                                 
                                  imshow(我)
                                 
                                  I = adapthisteq(我);
                                 
                                  hy = fspecial (“索贝尔”);
                                 
                                  hx = hy”;
                                 
                                  Iy = imfilter(双(I),衔接,“复制”);
                                 
                                  第九= imfilter (hx双(I),“复制”);
                                 
                                  第九gradmag =√。^ 2 + Iy。^ 2);
                                 
                                  %,imshow (gradmag[]),标题(“梯度大小(gradmag)”)
                                 
                                  L =分水岭(gradmag);
                                 
                                  Lrgb = label2rgb (L);
                                 
                                  %,imshow (Lrgb)、标题(“梯度的分水岭变换级(Lrgb)”)
                                 
                                  se = strel (“磁盘”2);
                                 
                                  Io = imopen(我,se);
                                 
                                  %,imshow (Io)、标题(“开放(Io)”)
                                 
                                  Ie = imerode(我,se);
                                 
                                  Iobr = imreconstruct(即我);
                                 
                                  %,imshow (Iobr)、标题(“Opening-by-reconstruction (Iobr) ')
                                 
                                  国际奥委会= imclose (Io, se);
                                 
                                  %,imshow (Ioc)、标题(“开合(Ioc)”)
                                 
                                  Iobrd = imdilate (Iobr, se);
                                 
                                  Iobrcbr = imreconstruct (imcomplement (Iobrd) imcomplement (Iobr));
                                 
                                  Iobrcbr = imcomplement (Iobrcbr);
                                 
                                  %,imshow (Iobrcbr)、标题(“开合的重建(Iobrcbr)”)
                                 
                                  女性生殖器切割= imregionalmax (Iobrcbr);
                                 
                                  %,imshow(女性生殖器切割)、标题(“区域最大值开合的重建(女性生殖器切割)”)
                                 
                                  I2 =我;
                                 
                                  I2(女性生殖器切割)= 255;
                                 
                                  %,imshow (I2)、标题(“区域极大值叠加在原始图像(I2) ')
                                 
                                  se2 = strel(的(5,5));
                                 
                                  fgm2 = imclose(女性生殖器切割,se2);
                                 
                                  fgm3 = imerode (fgm2 se2);
                                 
                                  fgm4 = bwareaopen (fgm3 20);
                                 
                                  I3 =我;
                                 
                                  I3 (fgm4) = 255;
                                 
                                  %,imshow (I3)
                                 
                                  标题(“修改区域极大值叠加在原始图像(fgm4) ')
                                 
                                  bw = im2bw (Iobrcbr graythresh (Iobrcbr));
                                 
                                  图,imshow (bw)、标题(的阈值开合重建(bw) ')
                                 
                                  [labeledImage, numberOfBlobs] = bwlabel (bw);
                                 
                                  blobMeasurements = regionprops (labeledImage,“重心”,“EquivDiameter”);
                                 
                                  EquivDiameter = [blobMeasurements.EquivDiameter];
                                 
                                  ValidDia =找到(EquivDiameter > 10);%的一些经验阈值
                                 
                                  blobMeasurements = blobMeasurements (ValidDia);
                                 
                                  EquivDiameter = EquivDiameter (ValidDia);
                                 
                                  %的概率估计的大小
                                 
                                  图,柱状图(EquivDiameter“归一化”,“概率”),包含(的大小(m -未校准)),ylabel (“概率”)
                                 
                                  %显示一些粒子
                                 
                                  图中,imshow (bw);持有在
                                 
                                  %显示区域形象
                                 
                                  为idx = 1:2:长度(ValidDia)%遍历所有气泡。
                                 
                                  质心= [blobMeasurements (idx) .Centroid (1) blobMeasurements (idx) .Centroid (2)];
                                 
                                  DiamSize = [“直径= '。,num2str ((EquivDiameter (idx))));
                                 
                                  文本(质心(1),重心(2),DiamSize,“颜色”,“c”);
                                 
                                  viscircles(质心EquivDiameter (idx) / 2);
                                 
                                  结束

蒂亚戈恩里克戈麦斯Lobato 2020年5月10日

编辑:蒂亚戈恩里克戈麦斯Lobato 2020年5月10日

最长的粒子的大小可以实现MajorAxisLength属性

                                  regionprops (labeledImage“重心”,“EquivDiameter”,“MajorAxisLength”);
                                 

获得分配一个好处可以考虑的是,你不需要所有粒子对,而是只有一个子集。我说在最初的回答将很难找到所有和我不相信这不该是你的目标。

为了提高结果您可以考虑,岩石的主要问题是它们重叠太多。所以像大津阈值方法( graythresh (Iobrcbr)) 最有可能不会工作,试图关闭太多的漏洞。我认为最好的方法,在一个二元分割模式,是使用一个非常高的二进制阈值或高数量的侵蚀,以尽可能独立的十字路口,他们依赖的子集的假设,所以即使你可以删除一些岩石结果将仍然是合理的。也请注意 idx = 1: 2 :长度(ValidDia)。情节所有斑点去除这2,我最初让它所以图像不会有太多了。