如何在 Python 中计算多边形的质心答案

【问题标题】：How to compute the centroid of a polygon in Python如何在 Python 中计算多边形的质心
【发布时间】：2014-06-28 07:36:58
【问题描述】：

我想知道是否有任何算法可以在 OSM 中计算多边形的中心，因为我发现每个多边形都有不同的参数表达式：

"POLYGON((-171379.35 5388068.23,-171378.8 5388077.59,-171368.44 5388076.82,-171368.89 5388067.46,-171379.35 5388068.23))"

"POLYGON((-171379.3 5374226.42,-171375.96 5374227.32,-171378.95 5374239.82,-171365.69 5374243.03,-171364.16 5374237.14,-171365.92 5374236.76,-171364.26 5374230.26,-171362.67 5374230.63,-171360.11 5374220.19,-171376.6 5374216.13,-171379.3 5374226.42))"

在这个表达式中，每个表达式都有不同数量的参数，那么我如何计算它的质心呢？

我正在使用 postgis 从空间信息中获取多边形信息。

【问题讨论】：

我使用 numpy 找到了这段代码：oles-tutorials.googlecode.com/svn-history/r42/trunk/scipy2011/…
ST_Centroid，由 MaM 发布，会给你答案，但如果你想了解公式/理论，请查看底部的 en.wikipedia.org/wiki/Centroid。

标签： python openstreetmap postgis

【解决方案1】：

你可以使用PostGIS functionality或者例如python库shapely来实现它。

【讨论】：

【解决方案2】：

使用WKT 描述多边形。

使用 Python？使用Shapely。

from shapely import wkt
p1 = wkt.loads("POLYGON((-171379.35 5388068.23,-171378.8 5388077.59,
                -171368.44 5388076.82,-171368.89 5388067.46,
                -171379.35 5388068.23))")
p2 = wkt.loads("POLYGON((-171379.3 5374226.42,-171375.96 5374227.32,
                -171378.95 5374239.82,-171365.69 5374243.03,
                -171364.16 5374237.14,-171365.92 5374236.76,
                -171364.26 5374230.26,-171362.67 5374230.63,
                -171360.11 5374220.19,-171376.6 5374216.13,
                -171379.3 5374226.42))")
print(p1.centroid.wkt) # POINT (-171373.8710815240337979 5388072.5175872016698122)
print(p2.centroid.wkt) # POINT (-171370.4190352254081517 5374229.0108996396884322)

【讨论】：