How do you find the centroid of a particular Blob?

90 views (last 30 days)
I have the following code:
clc;
clear all;
% This part reads and shows the image selected
I = imread('115.2.jpg');
imshow(I)
% Ginput is used to select the Hinge position
[x,y] = ginput(1)
% drawpoint is used to mark the Hinge position
h = drawpoint('Position',[x y]);
% Convert the RGB image to a Grayscale image
G=im2gray(I);
level = graythresh(G)
% Binarize the grayscale image
BW = imbinarize(G,level);
% Display the Grayscale image
figure
title('Grayscale Image')
imshow(G)
% Display the Binarized image
figure
title('Binarized')
imshow(BW)
% Fill holes
BW2 = imfill(BW,'holes');
% Display the Image after filling holes
figure
imshow(BW2)
title('Filled Image')
% Mark the Hinge here as well
h = drawpoint('Position',[x y]);
imwrite(BW2, 'BW2.jpg')
% Find the centroid
measurements = regionprops(BW2, 'Centroid');
imshow(BW2);
centroids = cat(1,measurements.Centroid);
hold on;
plot(centroids(:,1), centroids(:,2), 'r+', 'MarkerSize', 30, 'LineWidth', 1);
hold off;
where the image 115.2 is
And on running that code, I get:
How do I find the coordinates of the centroid of the rectangular blob? How do i neglect all centroids except the centroid of the rectangular blob?

Accepted Answer

Gatech AE
Gatech AE on 26 Jan 2021
Since you're working with a binary image here, BW2, you can use the bwpropfilt function to clean up the small areas.
BW2 = bwpropfilt(BW2, 'Area', 1);
The last input argument determines the largest # of regions to keep, so it can be lowered to just the one. From there regionprops would only detect the single centroid.
You could also pass area as an input to regionprops as it is and use the centroid of the largest area.
measurements = regionprops(BW2, 'Centroid', 'Area');
[~,ind] = max(cat(1,measurements.Area));
centroid = measurements(ind).Centroid;
  1 Comment
KLETECH MOTORSPORTS
KLETECH MOTORSPORTS on 29 Jan 2021
Thanks, I made some edits and it worked.
I wanted to know how to save the centroid coordinates in the (x1,y1) format?
This is what i've done so far:
clc;
clear all;
% This part reads and shows the image selected
I = imread('115.2.png');
imshow(I)
% Ginput is used to select the Hinge position
[x,y] = ginput(1)
% drawpoint is used to mark the Hinge position
h = drawpoint('Position',[x y]);
% Convert the RGB image to a Grayscale image
G=im2gray(I);
level = graythresh(G)
% Binarize the grayscale image
BW = imbinarize(G,level);
% Display the Grayscale image
figure
title('Grayscale Image')
imshow(G)
% Display the Binarized image
figure
title('Binarized')
imshow(BW)
% Fill holes
BWfill = imfill(BW,'holes');
% Display the Image after filling holes
figure
imshow(BWfill)
title('Filled Image')
% Mark the Hinge here as well
h = drawpoint('Position',[x y]);
imwrite(BWfill, 'BWfill.png')
%Filter the rectangular bar from the image by Area
BWarea = imread('BWfill.png');
BW2 = bwareafilt(BWarea,[83800 83900]);
imwrite(BW2, 'BW2.png')
imshow("BW2.png")
% Find the centroid
measurements = regionprops(BW2, 'Centroid');
imshow("BW2.png");
centroids = cat(1,measurements.Centroid);
hold on;
plot(centroids(:,1), centroids(:,2), 'r+', 'MarkerSize', 30, 'LineWidth', 1);
h = drawpoint('Position',[x y]);
[A] = centroids %how do i save the centroid position as coordinates, and not a variable?
hold off;

Sign in to comment.

More Answers (1)

Image Analyst
Image Analyst on 26 Jan 2021
Edited: Image Analyst on 26 Jan 2021
I think I already did this swinging pendulum panel thing before when you had a video. Not sure why you're posting it again. But anyway to get the centroid of the largest blob
% Get largest blob
bw2 = bwareafilt(bw2, 1);
% Find the centroid
measurements = regionprops(BW2, 'Centroid');
x = measurements.Centroid(1);
y = measurements.Centroid(2);
hold on;
plot(x, y, 'r+', 'LineWidth', 2, 'MarkerSize', 50)

Community Treasure Hunt

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

Start Hunting!

Translated by