从R中的shapefile中提取50m方块的坐标

时间:2017-11-07 18:19:39

标签: r shapefile

我正在尝试从环境机构提供的数据(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?

0 个答案:

没有答案