r - 找到多边形的视觉中心

Question

跟进这个问题:gCentroid shifting centroid towards concentration of points ，但形状更复杂。

使用这个 SpatialPointsDataFrame我使用 gCentroid 计算了中心，但它位于形状之外。

# c. = [dput'd data from link]
poly <- SpatialPolygons(Srl = list(Polygons( srl = list(Polygon(coords = coordinates(c.))), ID = 1)))
plot(poly)
axis(1)
axis(2)
points(rgeos::gCentroid(poly)) 

如何返回落在指定形状内部的中心点的坐标？

正如这些形状的评论中所指出的，中心位于边界之外。我认为我正在寻找的东西最好描述为“质心”，即距离边界最远的多边形内点。

Answer 1

polylabelR::poi (多边形的不可访问性极点(视觉中心))似乎可以解决问题。相关方法:https://blog.mapbox.com/a-new-algorithm-for-finding-a-visual-center-of-a-polygon-7c77e6492fbc

poly <- SpatialPolygons(Srl = list(Polygons( srl = list(Polygon(coords = coordinates(c.))), ID = 1)))
plot(poly)
axis(1)
axis(2)
points(rgeos::gCentroid(c.)) 

p <- polylabelr::poi(xx$long, xx$lat, precision = 0.01)
points(p)