Как разместить данные PolySet в пакете PBSmapping?

Question

Это репост из gis.stackexchange .Вопрос довольно длинный и подробный , поэтому я не буду его повторять.Простая проблема в том, что я не могу найти способ поднабора данных PolySet, используемых в пакете PBSmapping.Одним из вариантов будет преобразование в SpatialPolygons с maptools , его поднабор и преобразование обратно в PolySet, но, очевидно, это всего лишь обходной путь.Любая помощь будет принята с благодарностью.

Answer 1

Вы ищете атрибут:

df[attr(df, "PolyData")$NAME_1 == 'Bumthang',]

EDIT:

Будет ли это работать?

Аргумент col обозначает цвет и связан со значением PID.

df$Z='white'
df[attr(df, "PolyData")$NAME_1 == 'Bumthang',]$Z ="red"

plotPolys(df,col=df$Z)

Вы также можете поиграть с аргументом border, который дает цвет границ. Если вы поставите «белый», то единственным регионом, который будет отображаться, будет Бумтханг.

Примерно так будет работать:

df$W = 'white'
#borders of the region
df[attr(df, "PolyData")$NAME_1 == 'Bumthang',]$W ="red"
plotPolys(df,col=df$Z,border=df$W)

Answer 2

Как насчет использования другой библиотеки для подмножества файла?Вы можете поместить данные с помощью rgdal, сохранить в файл, импортировать шейп-файл, содержащий подмножество, и построить его с помощью PBSMapping.

library(rgdal)
library(PBSmapping)

df <-  readOGR(".","df")
subset <- df[df$NAME_1=="Bumthang",]
writeOGR(subset, ".", "bumthang", driver="ESRI Shapefile")
bum = importShapefile("bumthang.shp")
plotPolys(bum, projection=TRUE)

Или вы можете вообще отказаться от PBSmapping.

 plot(subset, axes=TRUE)