data(Europe)
data(NLD_prov)
data(NLD_muni)
tm_shape(Europe) + tm_borders()
tm_shape(NLD_muni) +
tm_fill(col="population", convert2density=TRUE, style="kmeans") +
tm_borders("grey25", alpha=.5) +
tm_shape(NLD_prov) +
tm_borders("grey40", lwd=2) +
tm_layout_NLD(title="Population (per km2)", bg.color="white", draw.frame = FALSE)
Run the code above in your browser using DataLab