在小册子R

Question

我需要将丹麦的地图划分为区域（有5个地区：Nordjylland，Midtjylland，Sydjylland，Sjælland和Hovedstaden）然后为这些地区着色，以便不同地区能够清晰地突出。我必须使用传单，因为它有其他功能，我将在以后使用。我在naturalearthdata.com上找到了一张地图，我认为我可以使用，但我无法弄清楚如何着色（甚至指示）这些区域。我试过的代码在

下面

library(rgdal)
library(leaflet)

download.file(file.path('http://www.naturalearthdata.com/http/',
                        'www.naturalearthdata.com/download/50m/cultural',
                        'ne_50m_admin_1_states_provinces_lakes.zip'), 
              f <- tempfile())
unzip(f, exdir=tempdir())

world <- readOGR(tempdir(), 'ne_50m_admin_1_states_provinces_lakes', encoding='UTF-8')
DK <- subset(world, name=="Denmark")
leaflet() %>% addTiles() %>% addTopoJSON(DK, weight = 1, color = "#444444", fill = TRUE)

如何使用naturalearthdata.com数据绘制不同国家的地区/州/省？我见过一个非常好的例子 http://www.56n.dk/kort/dk2050kort_age.html 但是没有可用的示例代码。

我在这里也找到了一个非常好的例子：https://rpubs.com/walkerke/leaflet_choropleth - 但我需要一张丹麦地图。

更新：我在http://www.kortforsyningen.dk找到了一个shapefile，可以解决问题。所以现在我的问题是如何将自己的数据与shapefile结合起来并将其绘制在传单中？如果我只是把

DK <- readOGR(".../shape", layer="REGION")
leaflet(data=DK)

我得到一个空白的屏幕......