Getting distance of centroids from each other.

Question

Sean Ivan Roxas 2022-4-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1692310-getting-distance-of-centroids-from-each-other

评论： Sean Ivan Roxas 2022-4-10

Is it possible to get the distances between each centroid? For example I have 5 centroids (centroid 1 - 5) in an image and I would like to get the distance between centroid 1 and centroid 2, then centroid 1 and centroid 3, then centroid 1 and centroid 4, then centroid 1 and centroid 5, then for centroid 2 and centroid 3, and so on. Then I will store these distances in an adjacency list or an adjacency matrix so that I could create a graph network.

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

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2022-4-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1692310-getting-distance-of-centroids-from-each-other#answer_938450

use pdist()

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

Sean Ivan Roxas 2022-4-9

在 MATLAB Online 中打开

I have this chunk of code already that prints a line between connected coins, but this only takes into account adjacent centroids not all. Can you please edit it such that it can also give me the distances between all the coins and store it in an adjacency matrix/list?

%% Blob Analysis
hBlobAnalysis = vision.BlobAnalysis('MinimumBlobArea', 10000,'MaximumBlobArea',150000);
[objArea,objCentroid,bboxOut] = step(hBlobAnalysis, IBWopen);
%% Print Centroid
CC = bwconncomp(IBWopen);
S = regionprops(CC, 'centroid');
centroids = cat(1, S.Centroid);
NodeCount = length(objCentroid)
figure
subplot(1,2,2)
imshow(pic)
hold on
plot(centroids(:, 1), centroids(:, 2), 'b*','LineWidth', 5)
hold off
%% Print Lines 
props = regionprops(IBWopen , 'Centroid', 'EquivDiameter');
meanDiam = mean([props.EquivDiameter]);
xy = vertcat(props.Centroid);
distances = pdist2(xy, xy);
hold on;
for k = 1 : size(distances, 1)
    thisRow = distances(k, :);
    closeIndexes = find(thisRow > 0 & thisRow < 1.12 * meanDiam);
   
    for k2 = 1 : length(closeIndexes)
        x1 = xy(k, 1);
        x2 = xy(closeIndexes(k2), 1);
        y1 = xy(k, 2);
        y2 = xy(closeIndexes(k2), 2);
        plot([x1, x2], [y1, y2], 'k-', 'LineWidth', 2)
    end
    
end
hold off;