# How to find a figure's centroid inside a polyshape?

4 visualizaciones (últimos 30 días)
David Franco el 11 de Mzo. de 2020
Comentada: David Franco el 10 de Oct. de 2021
I have the following problem:
I am plotting the boundaries of several municipalities and I need to list each one.
As there are many municipalities, I created a script that does this numbering automatically.
I used the centroid function to find a better position automatically, but some municipalities, like the one shown below, have a concave shape and the centroi is outside the limits of the municipalities.
My question is: there is a automatic way to find a position as central as possible within the city limits?
This is the result I get when I use the centroid of each polyshape to put the number. But for some municipalities with irregular shapes the numbers may be outside the limits (not well centralized):
Ps.: The boundaries of the municipalities are polyshapes and I have all the information regarding the location of the borders (I attached the polyshape for this example).
##### 2 comentariosMostrar NingunoOcultar Ninguno
Matt J el 11 de Mzo. de 2020
Where would you consider the "center" of the above shape to be? Could you mark it on your posted figure?
David Franco el 11 de Mzo. de 2020
In fact Matt, anywhere within the limits is already viable for me.

Iniciar sesión para comentar.

Matt J el 11 de Mzo. de 2020
In fact Matt, anywhere within the limits is already viable for me.
If so, you could use nearestvertex command to project the centroid to the nearest point within city limits,
##### 5 comentariosMostrar 3 comentarios más antiguosOcultar 3 comentarios más antiguos
Matt J el 11 de Mzo. de 2020
Editada: Matt J el 11 de Mzo. de 2020
Perhaps as follows,
>> city=polyshape([0 0; 0 1; 1 0]);
>> A=area(city);
>> p = .8; %percentage
>> [d,res]=fzero(@(d)area(polybuffer(city,d))-p*A, 0)
d =
-0.0309
res =
8.0499e-09
>> innercity=polybuffer(city,d);
David Franco el 11 de Mzo. de 2020
Perfect! Thank you!

Iniciar sesión para comentar.

### Más respuestas (1)

Chad Greene el 5 de Oct. de 2021
I just ran into this problem when trying to place a text label in the middle of a crescent-shaped ice shelf. The centroid or the mean or median of the coordinates of the ice shelf polygon are all outside the bounds of the ice shelf. Here's the best solution I could come up with:
% Convert the outline to a polyshape:
P = polyshape(x,y);
% And get the delaunay triangulation of the polygon:
T = triangulation(P);
% Now find the center points of all the triangles:
[C,r] = circumcenter(T);
% Get the index of the centerpoint that has the largest radius from the boundary:
[~,ind] = max(r);
% These center coordinates are in the center of the fattest part of the polygon:
xc = C(ind,1);
yc = C(ind,2);
##### 5 comentariosMostrar 3 comentarios más antiguosOcultar 3 comentarios más antiguos
Chad Greene el 8 de Oct. de 2021
Alrighty, @David Franco, you inspired me to add a new polycenter function to the Climate Data Toolbox for Matlab.
[xc,yc] = polycenter(Snew);
mapshow(Snew,'FaceColor','c');
hold on
text(xc,yc,num2str((1:399)'),'fontsize',7,'horiz','center','vert','middle')
David Franco el 10 de Oct. de 2021

Iniciar sesión para comentar.

### Categorías

Más información sobre Elementary Polygons en Help Center y File Exchange.

R2019b

### Community Treasure Hunt

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

Start Hunting!

Translated by