如何根据 R 中的属性删除 SpatialPolygonsDataFrame 的特定特征？答案

【问题标题】：How to remove specific features of SpatialPolygonsDataFrame based on attributes in R?如何根据 R 中的属性删除 SpatialPolygonsDataFrame 的特定特征？
【发布时间】：2021-10-03 15:42:42
【问题描述】：

由于不喜欢在已回答的已存在问题下提出问题，因此我的问题链接到此问题 URL：Simple way to subset SpatialPolygonsDataFrame (i.e. delete polygons) by attribute in R，您还可以在其中获取数据集。

有人知道如何根据值向量从 SpatialPolygonsDataFrame 中删除多个特征？

head(world.map@data)
#   FIPS ISO2 ISO3 UN                NAME   AREA  POP2005 REGION SUBREGION     LON     LAT
# 0   AC   AG  ATG 28 Antigua and Barbuda     44    83039     19        29 -61.783  17.078
# 1   AG   DZ  DZA 12             Algeria 238174 32854159      2        15   2.632  28.163
# 2   AJ   AZ  AZE 31          Azerbaijan   8260  8352021    142       145  47.395  40.430
# 3   AL   AL  ALB  8             Albania   2740  3153731    150        39  20.068  41.143
# 4   AM   AM  ARM 51             Armenia   2820  3017661    142       145  44.563  40.534
# 5   AO   AO  AGO 24              Angola 124670 16095214      2        17  17.544 -12.296

我想删除名称为“芬兰”、“挪威”、“瑞典”、“丹麦”的特征（无意冒犯，我只是没有斯堪的纳维亚半岛的数据；） )。

以下方法失败：

world.map <- world.map[world.map@data$NAME != c('Finland', 'Norway', 'Sweden', 'Denmark'), ]
# does not work at all

world.map <- dplyr::filter(world.map@data, NAMEN %in% c('Finland', 'Norway', 'Sweden', 'Denmark'))
# results in a dataframe, spatial information is lost

我猜想有一个使用base::which() 和%in% 的方法，但我不知道如何继续。 @Jeffrey Evans 在这篇文章中的答案参见 URL:https://gis.stackexchange.com/questions/57844/subset-a-spatialpolygonsdataframe-by-id-in-r 可以提供有用的提示...

【问题讨论】：

标签： r dataframe geospatial sf sp

【解决方案1】：

您可以按行 ID 子集 SpatialPolygonsDataFrame：

# row numbers of countries, exluding scandinavia
ids <- which(
!(world.map$NAME %in% c('Finland', 'Norway', 'Sweden', 'Denmark'))
)

# subset
not_scandinavia <- world.map[ids, ]

# plot
plot(not_scandinavia)

【讨论】：