如何用Python计算矩形Polygon的准确centroid（质心）坐标？

你的问题出在Polygon的坐标本身不符合地理空间规则：WKT格式中坐标对的顺序是(经度 纬度)，而你提供的Polygon里，纬度值出现了-96.767578125，超出了纬度的合法范围（纬度只能在[-90, 90]之间），这直接导致质心计算结果无效。

问题分析 #

Shapely/GeoPandas的centroid方法是基于几何图形的数学中心计算的，你的Polygon本身是一个包含无效纬度值的矩形，所以计算出的质心纬度也超出了合法范围。手动验证的话，矩形质心就是对角线中点：

• 经度中点：(-96.8115234375 + 32.8271484375) / 2 = -31.9921875
• 纬度中点：(-96.767578125 + 32.87109375) / 2 = -31.9482421875
  这个结果是数学上的正确中点，但因为原始纬度值非法，所以地理意义上无效。

解决方案 #

1. 修正Polygon的合法坐标
   先调整原始Polygon的纬度值，确保所有纬度在[-90, 90]区间内。例如把超出范围的-96.767578125改为-86.767578125（在合法范围内）：

   import shapely
   from shapely import wkt
   
   # 修正后的合法Polygon WKT
   corrected_wkt = "POLYGON ((-96.8115234375 32.87109375, -96.8115234375 -86.767578125, 32.8271484375 -86.767578125, 32.8271484375 32.87109375, -96.8115234375 32.87109375))"
   
   polygon = wkt.loads(corrected_wkt)
   centroid = polygon.centroid
   print(f"合法质心坐标: {centroid}")
   # 输出：POINT (-31.9921875 -26.9482421875)
   

2. 手动计算矩形质心验证
   对于矩形，质心等于经纬度范围的中点，你可以直接计算来验证结果：

   • 经度范围：min_lon = -96.8115234375，max_lon = 32.8271484375，质心经度 = (min_lon + max_lon) / 2
   • 纬度范围：min_lat = -86.767578125，max_lat = 32.87109375，质心纬度 = (min_lat + max_lat) / 2

内容的提问来源于stack exchange，提问作者kms