r - 使用 R 中的 shapefile 显示区域的子集

Question

我有一个英国的 shapefile:https://geoportal.statistics.gov.uk/Docs/Boundaries/Local_authority_district_(GB)_2014_Boundaries_(Generalised_Clipped).zip

我已将 shapefile 读入变量，UK

>UK <- readOGR(dsn = "....."
>England <- UK

我只想显示英语地方政府区域。它们在 LAD_DEC_2014_GB_BGC.dbf 中指定，其中 LAD14CD 以“E”开头

>UK@data

      LAD14CD                      LAD14NM                 LAD14NMW

0   E06000001                   Hartlepool                     <NA>
1   E06000002                Middlesbrough                     <NA>
2   E06000003         Redcar and Cleveland                     <NA>
371 W06000015                      Cardiff                 Caerdydd

>#filter UK@data and replace England@data with only English regions
>England@data <- UK@data$LAD14CD[c(grep("^E", UK$LAD14CD))]

>plot(England)

但是 grep 命令似乎将 shapefile 更改为一个因子，这意味着绘图如下所示:

Answer 1

使用此命令:

England <- UK@data$LAD14CD[c(grep("^E", UK$LAD14CD))]

...您仅对数据槽中的一列进行子集化，而不是对整个 shapefile 进行子集化，并将其分配给英格兰。

这应该可以完成工作:

England <- UK[grep("^E", UK@data$LAD14CD),]

注意，您需要在其中添加逗号!另外，您不需要将 grep 语句包装在 c() 中，但这并没有什么坏处，只是没有必要。