我想知道是否有任何算法我可以在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从空间信息中获取多边形信息。
答案 0 :(得分:4)
使用WKT描述多边形。
使用Python?使用匀称。
from shapely.wkt import loads as load_wkt
p1 = load_wkt("POLYGON((-171379.35 5388068.23,-171378.8 5388077.59,
-171368.44 5388076.82,-171368.89 5388067.46,
-171379.35 5388068.23))")
p2 = load_wkt("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)
答案 1 :(得分:1)
您可以使用PostGIS functionality或python lib shapely来实现它。