2015-12-14 4 views
1

Я пытаюсь добавить слои из статистического источника, используя функцию addProviderTiles. Данные поступают из этого источника: http://www.cbs.nl/nl-NL/menu/themas/dossiers/nederland-regionaal/publicaties/geografische-data/archief/2015/wijk-en-buurtkaart-2014-art.htmR листовка Rstudio показывает пустую карту, используя addProviderTiles Листовка

Мой код выглядит следующим образом:

library(leaflet) 
library(rgdal) 

wijken2 <- readOGR("C://", layer="buurt_2014", verbose=FALSE) 
newijken <- subset(wijken2, wijken2$GM_NAAM %in% c("Amsterdam")) 
pal <- colorNumeric(palette = "Blues",domain = newijken$AUTO_TOT) 

leaflet(newijken) %>% addPolygons(stroke= FALSE, fillOpacity = 0.5, smoothFactor = 0.5, color= ~pal(AUTO_TOT)) %>% addProviderTiles("CartoDB.Positron") 

В результате основная пустая карта CartoDB.

Я пробовал аналогичный подход, используя данные от https://rstudio.github.io/leaflet/shapes.html. Это показывает добавленные слои.

Есть ли что-то, что мне нужно изменить в данных?

Спасибо!

+0

без данных это трудно воспроизвести вашу проблему. Можете ли вы включить свой набор данных? – MLavoie

+0

См. Ссылку: http://www.cbs.nl/nl-NL/menu/themas/dossiers/nederland-regionaal/publicaties/geografische-data/archief/2015/wijk-en-buurtkaart-2014-art.htm , Сайт находится на голландском языке, но данные можно найти по ссылке Wijk-en Buurtkaart 2014 (67Mb). – rickwol

ответ

0

Ваша проекция была неправильной. В longlat кажется, что это работает, и похоже, что это в правильной стране :) Этого хотят вы хотите?

так запустить

newijken_latlng <- spTransform(newijken, CRS("+proj=longlat +datum=WGS84")) 
leaflet() %>% 
    addPolygons(data = newijken_latlng, color= ~pal(AUTO_TOT)) %>% 
    addProviderTiles("CartoDB.Positron") 

enter image description here

+0

Спасибо, ты абсолютный герой! Это отлично работает! – rickwol

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