本发明属于地图设计技术领域，涉及一种地图色彩的设计与配置方法，更具体涉及一种行政区划图自动着色的地图色彩效果优选方法。
背景技术：
色彩是重要的地图语言，在地图设计中作用不可忽视。行政区划图上行政区域需要用颜色进行区分，这些大面积的地图色彩往往决定了一幅地图的色彩风格。因此，行政区划图的色彩配置需要格外关注。
国内外学者目前主要研究行政区划图的“四色着色”问题，强调“行政区划图相邻色块颜色不相同”。常用的着色区域计算方法只是找到一个满足相邻政区色彩不相同条件下的解，没有较好的顾及着色区域间的面积均衡。另外，着色区域计算只顾及了地图要素间的空间关系，而色彩设计只顾及了色彩间的搭配原则。行政区划着色时，还需要考虑色彩在空间上的配置，不同的空间配置方案会引起人不同的地图色彩视觉感知。在实际地图制图过程中，制图员手工配置行政区划图色彩时往往耗时费力，且色彩配置有时难于符合高标准的制图要求。
经典的四色问题，传统的递归、回溯、贪心算法可以用来解决该问题；另外，随着计算机技术的发展，智能化算法也被应用于解决这个问题，如遗传算法等。地图色彩设计有相应的约束规则，为了使得地图色彩与所表达的信息一致，brewer等设计了适宜于地图着色的多种色彩模板集合，christophe等研发了相应的地图色彩设计工具。另外，色彩选择的方法和工具有很多，如hcsgenerator、colorschemeronline等。
技术实现要素：
本发明的目的在于，针对现有技术的空白，提出一种行政区划图自动着色的地图色彩效果优选方法，提高行政区划图着色时的效率与效果。
本发明所采用的技术方案是：一种行政区划图自动着色的地图色彩效果优选方法，其特征在于，包括以下步骤：
步骤1：利用回溯方法计算得到行政区划着色区域集合a＝{a1,a2,…,am}，其中，am表示有相同颜色的一组着色区域，并获取色彩模板集合c＝{c1,c2,…,cm}；
步骤2：对于着色区域集合a＝{a1,a2,…,am}，将色彩模板集合c＝{c1,c2,…,cm}分配到a中的指定区域共存在
具体实现包括以下子步骤：
步骤2.1：计算图面色彩协调程度sc；
步骤2.2：计算图面色彩对比程度，包括图面局部区域色彩对比程度sp与图面背景色彩对比程度sb；
步骤2.3：计算图面视觉平衡程度pb；
步骤2.4：综合考虑三个视觉因素对每一个色彩配置方案进行图面着色效果评价，选择其中最优的方案用于当前行政区划图的色彩配置。
本发明从地图色彩视知觉出发，详细研究了行政区划图着色过程中需要考虑的视知觉规律，并建立了相应的着色效果评价模型；在着色区域计算中，添加了面积均衡约束，在给定色彩模板集合前提下，设计并实现了顾及面积均衡的行政区划图自动着色算法。满足行政区划图的高标准的着色要求。
本发明量化研究了行政区划图色彩着色过程需要考虑的视觉因素，并建立相应的着色效果评价模型，能有效评价图面着色效果优劣，辅助提高地图图面着色效果。另外，实现了顾及面积均衡的地图着色区域计算，在给定色彩模板集合前提下，能快速实现行政区划图的色彩配置，提高工作效率，满足制图要求。
附图说明：
图1为本发明实施例的流程图；
图2为本发明实施例中顾及面积均衡的着色区域计算回溯方法数据处理流程图。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施示例对本发明作进一步的详细描述，应当理解，此处所描述的实施示例仅用于说明和解释本发明，并不用于限定本发明。
行政区划图的自动着色在给定色彩模板集合前提下，需要计算地图着色区域，并将色彩在着色区域中进行配置。地图着色区域计算需要顾及各色块的面积均衡，色彩的空间配置需要顾及人的视知觉规律，对相关的视知觉因素进行量化，并建立相应的评价模型，能实现色彩在着色区域中的配置。
请见图1，本发明提供的一种行政区划图自动着色的地图色彩效果优选方法，包括以下步骤：
步骤1：利用回溯方法计算得到地图着色区域集合a＝{a1,a2,…,am}，m为颜色数，am表示有相同颜色的一组着色区域。获取色彩模板集合c＝{c1,c2,…,cm}。其中，利用回溯方法计算地图着色区域时，添加了相应的面积均衡约束，顾及面积均衡求解m着色的着色区域计算回溯方法数据处理流程见图2，将行政区划图着色区域计算抽象为图着色问题，具体实现包括以下子步骤：
步骤1.1：将行政区划图转化为无环图g＝(v,e)，其中，将行政区划图中每个政区视为一个图节点，用vi表示；行政区划图中相邻区域对应两个图节点连接，连接边用en表示；v＝{v1,v2,…,vn}为图顶点集，e＝{e1,e2,…,en}为图边集；从色彩模板集合c中用尽可能少的颜色对顶点集v进行着色，且相邻节点颜色不相同；
步骤1.2：行政区划着色通常使用3-7个颜色为宜，设定初始颜色数m＝3。利用回溯法计算着色区域时，对每一个色彩赋权重，权重为该色彩所对应行政区划vi的图面面积avi。vi←cj表示在色彩模板集合中选择一个色彩cj对顶点vi着色，色彩选择时依次在剩余颜色中选择面积和最小的颜色，若顶点vi着色合法，且当前已经使用的色彩总数小于等于m，则对该颜色面积做累加，并移至下一个顶点继续着色。否则，m＝m+1，若m>7，则结束；否则重新选择有m个色彩的另一个色彩模板集合，回到步骤1。所述合法指对顶点vi着色时，相邻的顶点颜色不相同。
步骤2：对于着色区域集合a＝{a1,a2,…,am}，将色彩模板集合c＝{c1,c2,…,cm}分配到a中的指定区域共存在
具体实现包括以下子步骤：
步骤2.1：计算图面协调程度sc，其表达式为：
其中，m为配色区域个数，c(ai,aj)＝wij*η(ci,cj)，ci为区域ai颜色，cj为区域aj颜色，η(ci,cj)为利用色度学中的η模型计算出的两色彩间协调程度。wij为两区域相接部分长度，表示若两区域色彩越协调，且相接区域越长，则两区域色彩协调度大。
步骤2.2：计算图面色彩对比程度，包括图面局部区域色彩对比程度sp与图面背景色彩对比程度sb；图面局部区域色彩对比程度sp，其表达式为：
其中，m为配色区域个数，p(ai,aj)＝wij*d(ci,cj)，ci为区域ai颜色，cj为区域aj颜色，d(ci,cj)为两色彩在cielab色彩空间中的明暗度数值差异，wij为两区域相接部分长度，表示若两区域色彩间明暗度差异越大，且相接区域越长，则两邻近区域色彩对比度大。
计算图面背景色彩对比程度sb，其表达式为：
其中，n为配色区域个数，ci为区域si的颜色在cielab色彩空间中的明暗度，即若面积越大的区域明度越低，则全局色彩对比度越好。
步骤2.3：计算图面视觉平衡程度pb，其表达式为：
其中，其中，tx与ty为图面理论视觉中心的x，y坐标；px，py表示图面视觉中心，计算公式为：
其中，n为配色区域个数，san是区域an面积，xn和yn分别是对应区域an重心的x，y坐标，
步骤2.4：综合考虑三个视觉因素对色彩配置方案进行图面着色效果评价时采用的是多指标综合分析方法。计算
其中，fmin表示f的最小值，fmax表示f的最大值，f取值为sc、sp、sb和pb；
考虑到色彩对比有一定的上下阈值，并不是越大越好，或者越小越好，对于sp与sb，还需要进行阈值最优处理。计算公式为分别为：
其中，k1为图面局部区域色彩对比程度的最佳阈值，k2为图面背景色彩对比程度的最佳阈值，k1,k2∈(0,1)。
图面底色着色效果的优秀程度f(δ)计算公式为：
其中，w1，w2，w3分别为三个视觉因素对应的权重，且w1+w2+w3＝1。
本发明在充分研究地图区域面色色彩配置需要顾及的视觉规律基础上，对相关视知觉规律进行了量化，建立了相应的地图区域面色着色效果评价模型，能有效评价当前地图区域面色着色效果的优劣，自动提高图面着色效果，易于投入实用。
应当理解的是，本说明书未详细阐述的部分均属于现有技术。
应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。