我在地面上标出了一些地块，然后在这些地块上收集了无人机镜头。 在QGIS中，我绘制了一个多边形以标记投影正射影像上每个图的边界。当然，这些多边形远非完美的南北方向。 现在，我想找到该多边形的角度或旋转角度，以便可以将网格对齐到样图边界内的子样本，类似于以下问题： Creating a regular polygon grid over a spatial extent, rotated by a given angle
但是，由于我有100多个这样的地块，所以我不想尝试反复尝试来找到角度。我可以使用“定向的最小边界框”在QGIS中实现此功能，该框可以自动计算角度并将其输出到属性表。
在R中找不到有关如何执行此操作的任何信息。 任何帮助将不胜感激。
这里有一些代码可以绘制示例多边形
poly = Polygon(cbind(c(2,4,3,1,2),c(2,3,5,4,2)))
polygons = Polygons(list(poly), "s3")
SpP = SpatialPolygons(list(polygons))
plot(SpP)
