Progressed with the neighbour region colouring problem, see the commented code in the end of visualization.R