R中多边形的罗盘方向（北，西......）

Question

我有一个GIS对象的多边形，例如。长方形。我想根据罗盘方向提取所有边缘。我怎么能这样做？

这是一个简单的例子

poly = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
polygons = Polygons(list(poly), "s3")
SpP = SpatialPolygons(list(polygons))
plot(SpP)

例如上面的edgs。与北方有关。但我如何通过脚本识别出来？