Блог

Я строю шейп-файл полигона (SpatialPolygonDataFrame).Однако файл слишком большой, поэтому я попытался разделить его на основе названия места:

  library(maptools)
 mp3= readShapePoly('C:/Users/ghy_000c11a_e.shp')
 #subplot the shapefiles
 sub.shape = subset(mp1,mp1$NAME == "Fitz Hugh Sound")
 

Тем не менее, мне нужно добавить так много имен мест, поэтому мне было интересно, можно ли отобразить его на основе координат?

Я также пробовал 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")
 

Он отображает его в соответствии с указанным мной диапазоном координат, но я предполагаю, что он изменяет весь шейп-файл, поэтому, когда я пытаюсь его отобразить — все перепуталось.