1

我想从形状文件中自动创建两个变量：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

也许他们在其他地方（我对形状文件几乎没有经验）？是否有某个函数至少可以让我创建一个变量来指示一个区域是否没有外部边界（即上面地图中所有没有点的区域）？

4

0 回答 0