# 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
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?

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;
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
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
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;

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)
##### 2 CommentsShowHide 1 older comment
KLETECH MOTORSPORTS on 29 Jan 2021
Tried BW2, got the same error

R2020b

### Community Treasure Hunt

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

Start Hunting!