Карты с различными слоями в R

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

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

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

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

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

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

Я был бы очень признателен за вашу помощь!

Спасибо!

источник

2014-09-11 coolcats

Примечание: не используйте 'readShapeSpatial', использовать' rgdal 'package и' readOGR' - это будет считываться в любой проекционной информации в шейп-файле, который игнорируется 'readShapeSpatial'. – Spacedman

К оверлоту вы можете использовать аргумент sp.layout, см. this relevant question. Так что в вашем случае необходимо добавить следующие строки в ваш код:

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

который даст вам следующее choropleth:

enter image description here

источник

2014-09-17 10:43:46 BlankUsername

Отлично! Это работает!!! Спасибо! – coolcats

Карты с различными слоями в R

ответ

Смежные вопросы