tm_shape(NLD_muni) + # Shape (group 1)
tm_polygons( # Layer: data-driven fill
fill = "income_high",
fill.scale = tm_scale_continuous(values = "brewer.purples"),
fill.legend = tm_legend(title = "Income")) +
tm_shape(NLD_prov) + # Shape (group 2)
tm_borders(col = "black", lwd = 2) +
tm_text(text = "name") +
tm_basemap("Esri.WorldTerrain") + # auxiliary layer
tm_compass() + # map component
tm_scalebar()
