我正在绘制多边形shapefile(SpatialPolygonDataFrame)。但是,文件太大,所以我尝试根据地点的名称对其进行子绘图:
library(maptools)
mp3= readShapePoly('C:/Users/ghy_000c11a_e.shp')
#subplot the shapefiles
sub.shape = subset(mp1,mp1$NAME == "Fitz Hugh Sound")
但是,我需要添加这么多的平面名称,所以我想知道它是否是基于坐标的可能的tosubplot?
我也尝试了PBSmapping:
#install.packages("PBSmapping")
df <- importShapefile('C:/Users/ghy_000h11a_e.shp')
sub <- clipPolys(df, xlim = c(-100 , 0) , ylim = c(51.3 , 70))
plot (sub, col="blue")
它按照我指定的坐标范围对其进行子绘图,但我猜它会改变所有的shapefile,所以当我尝试绘制它时 - 所有都搞砸了。