Matlab中区域标签重映射：实现邻近区域差异化显示

问题背景 #

我有一个用于标识晶粒、虚构地图区域的标签矩阵，标签值范围可达数千。每个均质区域对应唯一标签值，地图不存在飞地，但可能存在仅有一个邻域的区域。示例矩阵如下：

[a a a c c d
 a a b c c d
 a b b d d d]

可视化痛点 #

尝试用现有配色方案（如lines、随机配色）可视化时，因区域尺寸差异较大，部分相邻区域颜色近乎一致，无法得到理想效果。

需求 #

是否存在一种算法，可对标签值进行重映射，使得相邻区域的标签值不同，同时保证重映射后的最大值尽可能小（类似政治地图的四色原理）？示例结果如下：

[1 1 1 3 3 1
 1 1 2 3 3 1
 1 2 2 1 1 1]

内容的提问来源于stack exchange，提问作者Crowley