我正在尝试从环境机构提供的数据(https://data.gov.uk/dataset/risk-of-flooding-from-rivers-and-sea1)中提取50平方米的坐标和与该方格相关的洪泛因子。当您从该页面点击到此页面时,我已经下载了shapefile格式(http://environment.data.gov.uk/ds/catalogue/#/8d57464f-d465-11e4-8790-f0def148f590)。
数据显示每50平方米的洪水因子。我对Shapefiles完全不熟悉,也是R的新手。从我读过的内容来看,我希望50m的正方形由'Polygons'表示并且使用
查看了多边形polys <- slot(data,"polygons")
然后使用
打印其中一些坐标for (i in 1:length(polys)) {
print(paste("Polygon #",i))
print(slot(slot(polys[[i]],"Polygons")[[1]],"coords" ))
}
我对输出感到困惑,因为我假设一个正方形将由四个点指定,但是,指定多边形的(对)坐标的数量变化很大。
这个假设是否正确?或者数据不是由它声称的50平方米组成?
如果它们确实是50平方米,是否有一种简单的方法来提取多边形中心的坐标及其ID?