Постройте карты с разными слоями в R

Я рисую некоторые данные для немецких округов в R, используя spplot. -Это отлично работает.

Чтобы было легче понимать / читать карту или находить места, я хотел бы добавить границы федеральных земель. Карты взяты из немецких границ

Я использую readShapeSpatial для загрузки данных, а затем использую spplot со слоем для построения карты. Однако я не уверен

krs <- readShapeSpatial("VG250_Kreise.shp")
bdl <- readShapeSpatial("VG250_Bundeslaender.shp")

krs@data$n <- cut(rnorm(402, 0, 1), breaks=5)
seq.red <- rev(sequential_hcl(n=nlevels(krs@data$n), h=0))
spplot(krs, "n", col=grey(.9), col.regions=seq.red, pretty=TRUE)+
layer(sp.polygons(bdl))

Я не уверен, какую часть bdl мне нужно добавить, чтобы получить только границы.

Буду очень признателен за вашу помощь!

Спасибо!


r maps spatial
person coolcats    schedule 11.09.2014    source источник
comment
Примечание: не используйте readShapeSpatial, используйте пакет rgdal и функцию readOGR - это будет считываться в любой информации о проекции в шейп-файле, которая игнорируется readShapeSpatial.   -  person Spacedman    schedule 17.09.2014

Ответы (1)


arrow_upward
0
arrow_downward

Чтобы перекрыть график, вы можете использовать аргумент sp.layout, см. этот соответствующий вопрос. Итак, в вашем случае вы должны добавить в свой код следующее:

bunder<-list("sp.polygons", bdl)
spplot(krs, "n", col=grey(.9), col.regions=seq.red, pretty=TRUE,
       sp.layout=bunder)

Это даст вам следующий хороплет:

введите описание изображения здесь

person horseoftheyear    schedule 17.09.2014