stackoverflow社区,
我有两个日本自治市的形状文件.我正在使用R为每个城市创建单独的地块.我可以使用一个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, …