som_cluster <- cutree(hclust(dist(som_model$codes)), 2)
my_palette <- c("#ffffff", '#dddddd', '#ccccff', '#aaaaaa')
plot(som_model, type = "mapping", labels = ud[,1], bgcol = my_palette[som_cluster])
add.cluster.boundaries(som_model, som_cluster)

