How to perform shape factor analysis (circularity) on list of x, y coordinates

Question

Suhaas Garre 2020-1-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/499135-how-to-perform-shape-factor-analysis-circularity-on-list-of-x-y-coordinates

评论： Adam Danz 2020-1-9

I have a list of x-y coordinates corresponding to the points along the perimeter of a shape. MATLAB's Image Processing Toolbox seems to allow for shape factor analysis of closed shapes within binary images, but I don't have enough points to fully close the shape.

So far, I've just plotted the x-y coordinates, but I'm having difficulty analyzing further. Are there built in mathematical methods for calculating max/min feret diameter and circularity using just a set of points?

Alternatively, are there other methods in image processing that are generally used to analyze shapes described by individual points?

7 个评论
显示 5更早的评论隐藏 5更早的评论

Suhaas Garre 2020-1-6

Yes, I certainly see what you mean. I’ll post some scatter plots tomorrow with some examples of what I’m referring to.

Thanks for your help so far!

Suhaas Garre 2020-1-7

编辑：Suhaas Garre 2020-1-7

I've updated the original question with an example of what I'm working with. I managed to close the shape yesterday, but as you can see, it is highly irregular. I've attached another example onto this comment as well.

请先登录，再进行评论。

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

Answer 1

Adam Danz 2020-1-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/499135-how-to-perform-shape-factor-analysis-circularity-on-list-of-x-y-coordinates#answer_409057

编辑：Adam Danz 2020-1-8

在 MATLAB Online 中打开

It looks like you're working with binary images in which case, you could use regionprops() along with the 'MajorAxisLength' and 'MinorAxisLength' properties to get the 'length' and 'width' of each blob. If those values are close to each other they can be categorized as square-ish. You'll have to choose a threshold such as

isSquareish = MinorAxisLength/MajorAxisLength > 0.95;

If you decide to go this route and get stuck, feel free to update us with your code so we can stay involved.

2 个评论
显示无隐藏无

Suhaas Garre 2020-1-8

I will give that a try, thank you for all your help!

I'll update this post as I make progress.

Adam Danz 2020-1-9

Sure, good luck!

请先登录，再进行评论。

Answer 2

Image Analyst 2020-1-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/499135-how-to-perform-shape-factor-analysis-circularity-on-list-of-x-y-coordinates#answer_409422

There is a new function bwferet() that you might want to look into. You might want to look into edge linking algorithms to close the shapes.

A circle could have the same major and minor axis lengths. I've found that looking at circularities (perimeter squared divided by 4*pi*area) or feret diameters is not that great a way but it could be one factor in the decision if a blob is square. One of the more reliable ways to find a square is to find the centroid of the points and then get the distances of all perimeter points to the centroid, then plot that and look for 4 peaks. Attached is my demo for shape recognition where I do that for a variety of polygons.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Adam Danz 2020-1-9

This is great.

BTW, I was assuming the OP wanted to classify elongated blobs from blobs without an obviously longer axis (hense the term square-ish). If the task is to identify shapes that are nearly square, then my simple minor/major axis idea isn't robust enough.

请先登录，再进行评论。

Answer 3

Meg Noah 2020-1-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/499135-how-to-perform-shape-factor-analysis-circularity-on-list-of-x-y-coordinates#answer_408917

在 MATLAB Online 中打开

For the final part of your question, 'Alternatively, are there other methods in image processing that are generally used to analyze shapes described by individual points?' one way is to raserize the data. This isn't a great way, but it gives ok results depending on the X, Y range and the accuracy you need:

% let inY and inX be your vectors <shapename>.Y, <shapename>.X
itheta = 0:360;
% fake data is a circle of radius R
R = 10.0;
inY = R.*sind(itheta); inX = R.*cosd(itheta);
h1 = figure(1);
gridDensity = 0.05; X1D = -2*R:gridDensity:2*R; Y1D = -2*R:gridDensity:2*R';
nX = numel(X1D); nY = numel(Y1D);
imagesc(X1D,Y1D,zeros(nY,nX),[0 255]);
set(h1,'position',[1 1 1000 800]);
set(gca,'units','pixels','position',[5 5 nX nY],'visible','off');
hold on;
colormap(gray);
axis equal
axis tight
plot(inY,inX,'color','w');
tim = getframe(gca); 
imorig = tim.cdata(:,:,1);
imorig(imorig>0) = 1;
% Note: If your data are in latitude/longitude values
%       then use vec2mtx to rasterize it
% rasterize the data values at grid size = 0.25
% gridDensity = 1/2; X1D = -2*R:0.25:2*R; Y1D = -2*R:0.25:2*R';
% [X2D,Y2D] = meshgrid(X1D,Y1D);
% [img, inRefVec] = vec2mtx(inY, inX, gridDensity);
% verify
h2 = figure(2);
imagesc(X1D,Y1D,imorig);
% apply standard blob detection and characterization
imblob = bwlabel(imorig,8);
stats = regionprops(imblob,'ALL');
diameterInPixels = norm(stats.BoundingBox([1 3]) - stats.Centroid);
diameterInDistance = diameterInPixels*gridDensity;
fprintf(1,'Diameter in Distance: %f\n',diameterInDistance);
 