r - 在 R 中绘制 shapefile 中的元素

Question

stackoverflow 社区，

我有两个日本城市的 shapefile。我正在使用 R 为每个城市创建单独的绘图。我可以使用一个形状文件来完成这项工作，但相同的语法却无法使用另一个形状文件。以下是数据的代码和 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，因此我无法绘制各个元素。

非常感谢您的帮助:在这个问题上我已经用头撞墙有一段时间了!

Answer 1

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

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

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

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

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

#FALSE  TRUE 
# 1894     8 

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

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

#FALSE  TRUE  <NA> 
# 1893     1     8 

包装在中解决了这个问题:

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),])