Which one is best for calculating circularity ?

61 ビュー (過去 30 日間)
Muhammet
Muhammet 2013 年 8 月 22 日
コメント済み: Image Analyst 2021 年 6 月 15 日
Hi all, I need to calculate objects roundness.I am not sure about which of the below ones is the best choice for this purpose.
- R= 4*pi*A/P^2; where A is area in pixels, P is Perimeter.
-Eccentricity
Thank you for your opinions.particles I need to measure are like this :
  4 件のコメント
RAKESH KUCHANA
RAKESH KUCHANA 2021 年 6 月 15 日
Can you explain code for how to remove small white dot areas (in the red colour bounded areas) in the given image?
Image Analyst
Image Analyst 2021 年 6 月 15 日
@RAKESH KUCHANA, no, not here in Muhammet's 8 year old question (he's probably not interested in your problem), but we will in your new question if you post it.

サインインしてコメントする。

採用された回答

Image Analyst
Image Analyst 2013 年 8 月 22 日
編集済み: Image Analyst 2021 年 6 月 15 日
I use circularity quite a bit. The formula I use is the inverse of the one you use:
props = regionprops(mask, 'Area', 'Perimeter');
allAreas = [props.Area];
allPerimeters = [props.Perimeter];
circularity = (allPerimeters .^ 2) ./ (4 * pi * allAreas);
so for me the circularity of a circle is 1 and the higher it goes, the less circular it it. For you the lower it goes (closer to zero) the less circular the blob is. You can do it either way.
I don't often use eccentricity because you can get the same eccentricity for a perfectly smooth ellipse as for some tortuous-shaped blob because it fits the perimeter shape to an ellipse. For example an asterisk and a disc might both have the same eccentricity but there is a dramatic difference in shape that is not picked up by eccentricity. Don't use eccentricity unless you know that your blobs are fairly close to ellipses.
  5 件のコメント
Rai_313
Rai_313 2015 年 4 月 24 日
hi Image Analyst ,
Im would like to ask, why shouldnt we use the haralick circularity?
thank you.
Image Analyst
Image Analyst 2021 年 6 月 15 日
it looks like it's the inverse of the coefficient of variation or the mean radius divided by the standard deviation of the radii = [(mean R) / (stddev R)]. So for a circle, the value is infinity since the standard deviation would be zero. He claims it's more accurate for small, digitized polygons, and that may be true.
I prefer the other definition - what he calls the Rosenfield formula P^2/(4*pi*Area). I like that it's simpler to compute plus it's nicely intuitive in that the value for a circle is 1, not some big undefined number. How big does the number need to be to be circular for Haralick? 10? 100? 1000? I don't know. I guess you could create randomly shaped blobs and compute it both ways and see what values each way produces.
But it comes down to whatever one you want to use. I don't think it's like one is better than another, though Haralick claims his is better for small circles where it's so quantized that it's more octagonal than circular. I mean a circle with a radius of 1 pixel would essentially be a cross + shape, not a circle. So if you have very small blobs then you might use Haralick's formula, though I might counter that comparing blobs based on circularity when the blobs are so tiny they can't really be circular in the first place is not a proper comparison.

サインインしてコメントする。

その他の回答 (2 件)

Amith Kamath
Amith Kamath 2013 年 8 月 22 日
I would do one of the following:
1. Use the 'solidity' option in REGIONPROPS to estimate how 'curvy' the blob is, and solidity is defined as the area/convexArea.
2. Use the 'MajorAxisLength' and 'MinorAxisLength' properties in REGIONPROPS and compare the two. If they are close to equal, the object would be more circular (in most cases).
3. Use 'Area' and 'Perimeter' from REGIONPROPS to estimate circularity in my custom code. http://www.mathworks.com/matlabcentral/newsreader/view_thread/292521 may be useful to see.

Shashank Prasanna
Shashank Prasanna 2013 年 8 月 22 日
  1 件のコメント
Muhammet
Muhammet 2013 年 8 月 23 日
I am not looking for perfect circular objects almost with same radius.Objects are differing both in size,in circularity.Thank you for your concern.

サインインしてコメントする。

カテゴリ

Help Center および File ExchangeGeographic Plots についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by