5

stackoverflow 社区，

我有两个日本直辖市的 shapefile。我正在使用 R 为每个城市创建单独的地块。我可以使用一个 shapefile 来完成这项工作，但另一个 shapefile 的相同语法会失败。以下是数据的代码和 URL：

library(sp)
library(maptools)

# Map A - this one works
# Please download: http://www.filefactory.com/file/z26nirxoz53/n/JPN_adm_zip

# Enter your path for readShapePoly
japanMapA = readShapePoly("JPN_adm/JPN_adm2")
names(japanMapA) 

plot(japanMapA[japanMapA$ID_2 == 1199,])

# Map B - this one doesn't work
# Please download: http://geocommons.com/overlays/173340.zip

# Again, enter your path for readShapePoly    
japanMapB = readShapePoly("japan_ver71") 
names(japanMapB)

plot(japanMapB[japanMapB$JCODE == 45382,])

它抛出的错误是：

Error in plot(japanMapB[japanMapB$JCODE == 45382, ]) : 
error in evaluating the argument 'x' in selecting a method for function 'plot': Error in japanMapB[japanMapB$JCODE == 45382, ] : 
NAs not permitted in row index

在这种情况下，我不知道如何删除 NA，因此我无法绘制单个元素。

非常感谢您的帮助：在这个问题上我的头撞墙了一段时间！

4

2 回答 2

6

我认为这归结为适当的子设置问题。

您的japanMapB对象由一些元数据和一系列多边形组成，每个形状存储在japanMapB@polygons. 所以你有了：

> length(japanMapB$JCODE)
#[1] 1902
> length(japanMapB@polygons)
#[1] 1902

正如@PaulHiemstra 指出的那样，您的变量NA中有一些值JCODE

> table(is.na(japanMapB$JCODE))

#FALSE  TRUE 
# 1894     8

这意味着在尝试为要绘制的城市建立索引时，您会得到 NA 结果。

> table(japanMapB$JCODE==45382,useNA="always")

#FALSE  TRUE  <NA> 
# 1893     1     8

包装which解决了这个问题：

which(japanMapB$JCODE == 45382)
#[1] 1802
#You will now select to plot only the 1802th polygon stored in the data object
plot(japanMapB[which(japanMapB$JCODE == 45382),])
于 2013-08-19T05:28:00.163 回答
3

可能发生的情况是存在带有NAid 的城市/多边形。例如：

id = c(1,2,3,NA)
id == 2
# [1] FALSE  TRUE FALSE    NA

请注意，通过这种比较，NA中的值id不会导致FALSE，而是NA。使用id == 2for subsetting 会导致您看到的问题。@thelatemail 将调用包装在 中which，它将逻辑向量转换为一组索引，其中逻辑向量为TRUE：

which(id == 2)
# [1] 3

请注意，调用which消除了NA. 这个索引向量可以安全地用于子集化。通常，您可以使用na.omit删除NA's：

l[l == 2]
# [1]  2 NA
l[na.omit(l == 2)]
# [1] 2
于 2013-08-19T05:22:59.983 回答