我想从形状文件中自动创建两个变量:1。如果某个区域具有国际边界,则为虚拟指示符; 2.如果某个区域具有沿海边界,则为虚拟指示符。
例如,几内亚变量1将是下面带有红点的区域,变量2将带有蓝点(我用眼睛做了这些)。
library(raster)
sd0 <- getData(name = "GADM", country = "GIN", level = 2)
plot(sd0)
对于这些类型的特征,@data
插槽中似乎没有任何信息:
head(sd0@data)
# OBJECTID ID_0 ISO NAME_0 ID_1 NAME_1 ID_2 NAME_2 HASC_2 CCN_2 CCA_2 TYPE_2 ENGTYPE_2 NL_NAME_2
# 1 1 97 GIN Guinea 1 Boké 1 Boffa GN.BF NA Préfecture Prefecture
# 2 2 97 GIN Guinea 1 Boké 2 Boké GN.BK NA Préfecture Prefecture
# 3 3 97 GIN Guinea 1 Boké 3 Fria GN.FR NA Préfecture Prefecture
# 4 4 97 GIN Guinea 1 Boké 4 Gaoual GN.GA NA Préfecture Prefecture
# 5 5 97 GIN Guinea 1 Boké 5 Koundara GN.KD NA Préfecture Prefecture
# 6 6 97 GIN Guinea 2 Conakry 6 Conakry GN.CK NA Préfecture Prefecture
也许他们在其他地方(我对形状文件几乎没有经验)?是否有某个函数至少可以让我创建一个变量3.指示一个区域是否没有外边界(即上面地图中没有点的所有那些)?