我将一个shp文件作为SF表导入到R中,我想将多边形与数据连接起来。 仅在2种情况下,数据就涉及一个包含多个多边形的城市区域。例如,数据所指的城市是一个以6结尾的6位数字的代码。但是我有一张地图,仅在两个城镇的情况下有城镇地区,代码以2到6结尾。
我想加入这些多边形。 首先,我替换城市区域的ID代码以在我的数据中仅代表一个城镇的代码,以1结尾 然后我想我应该按代码分组并按组加入多边形。我该怎么办?
map_gminy <-
map_gminy %>%
mutate(
gmina =
if_else(str_detect(kod_teryt,"10610"),
"Łódz",
gmina),
kod_teryt =
if_else(str_detect(kod_teryt,"10610"),
"106101",
kod_teryt)
) %>%
mutate(
gmina =
if_else(str_detect(kod_teryt,"12610"),
"Kraków",
gmina),
kod_teryt =
if_else(str_detect(kod_teryt,"12610"),
"126101",
kod_teryt)
) %>%
group_by(kod_teryt) %>%
sf_???????(gemometry)