How do I find the centre of gravity for an irregular shape?

30 次查看（过去 30 天）
Sudev Parthiban
Sudev Parthiban 2019-7-2
评论： Sudev Parthiban 2019-7-17
I have given the marked coordinates and region of interest in which I would like to find the centre of gravity for, please help me.
  2 个评论
Sudev Parthiban
Sudev Parthiban 2019-7-15
编辑：Jan 2019-7-15
It is an ultrasound image it is in .mat format... This is my code:
clc
clear
close all
cd('C:\Users\ADMIN\Desktop\Sudev\Objective 1')
load('ROI1.mat'); %Load data from file
figure(1); clf;
subplot(1,2,1);
imshow(paveikslas,[])
paveikslas(1:200,1:200,:) = 0;
figure(1); clf;
subplot(1,2,1);
imshow(paveikslas,[])
hold on
plot(x_poly, y_poly,'r','LineWidth',2)
title('Overlap of ROI on Bscan')
subplot(1,2,2);
plot(x_poly); hold on;
plot(y_poly,'-r');
title('Coordinates of ROI')
figure(2)
plot(x_poly, y_poly,'r','LineWidth',2)
axis equal
title('ROI to be calculated for Center and Axis')

请先登录，再进行评论。

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

采纳的回答

Bruno Luong
Bruno Luong 2019-7-15
You can get the coordinates of the centroid if you have a boundary coordinates of a polygonal region:
[xc,yc] = polygoncentroid(x_poly,y_poly)
using this function
function [xc,yc] = polygoncentroid(x,y)
xn = circshift(x,-1);
yn = circshift(y,-1);
A = x.*yn-xn.*y;
a = 3*sum(A);
xc = sum((x+xn).*A)/a;
yc = sum((y+yn).*A)/a;
end

更多回答（1 个）

Kaustav Bhattacharya
  4 个评论
显示 2更早的评论
Jan
Jan 2019-7-15
编辑：Jan 2019-7-15
The DICOM format is a file format. You have loaded the data from a MAT file. Therefore I think, that "DICOM format" is neither possible, nor does it explain, how the data are represented. All we know currently is that your data are not a binary image. Please provide any explicit details about the input. It would be inefficient if the readers guess, which data you have.
What do you exactly get for:
data = load('ROI1.mat')
Loading data directly into the workspace is prone to bugs. Prefer to catch the data in a variable like in my example.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 DICOM Format 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by