从R中的形状文件中找到海岸和国际寄宿生

时间:2017-04-04 03:40:38

标签: r gis r-raster sp maptools

我想从形状文件中自动创建两个变量:1。如果某个区域具有国际边界,则为虚拟指示符; 2.如果某个区域具有沿海边界,则为虚拟指示符。

例如,几内亚变量1将是下面带有红点的区域,变量2将带有蓝点(我用眼睛做了这些)。

library(raster)
sd0 <- getData(name = "GADM", country = "GIN", level = 2)
plot(sd0)

enter image description here

对于这些类型的特征,@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.指示一个区域是否没有外边界(即上面地图中没有点的所有那些)?

0 个答案:

没有答案