Trying to find the centroid of Polygon with Shapely but I do not understand how to tell the library that the vertices are (lon, lat). Perhaps I need to set a projection?
Here is my code:
from shapely.geometry import Polygon, Point # LON, LAT vertices = [ Point(-79.8726944444444, 8.68505555555556), Point(-79.8733888888889, 8.50419444444444), Point(-79.54552777777779, 8.68386111111111), Point(-79.54622222222221, 8.503), Point(-79.8726944444444, 8.68505555555556), ] p1 = Polygon(vertices) centroid = p1.centroid print(centroid) # POINT (-1805804163717.8823 6592764402.930745)
The result is clearly wrong.
It’s a degenerate shape. It’s bowtie that crosses itself. Shapely can’t handle that.