r - 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)

例如的上边缘。与北方有关。但是我如何通过脚本识别它？