CN111145302B - 顾及线宽一致性的矢量地图方形虚线符号绘制方法 - Google Patents

Abstract

本发明公开了一种顾及线宽一致性的矢量地图方形虚线符号绘制方法，包括以下步骤：对矢量线要素进行初始化，指定像素单位的线宽和虚实比例值，确定矢量线要素的世界坐标系和屏幕坐标系之间的转换规则。其次按照坐标转换规则将屏幕线宽转换为世界坐标系的线宽，将虚实比例值进行线宽映射，并在世界坐标系下对矢量线要素进行三角化操作。然后构建顾及线宽一致性的GPU状态协议，在CPU中获取初始最大线宽、地图缩放比例系数、虚实比例值等信息，利用状态协议传递到GPU中。最终，在GPU中获取状态协议后进行保持线宽一致的方形虚线符号绘制。本发明通过构建方形虚线符号绘制的状态协议，驱动CPU‑GPU协同工作，实现了高效地图矢量线要素绘制。

Description

顾及线宽一致性的矢量地图方形虚线符号绘制方法

技术领域

本发明涉及数字地图制图领域和地理信息***领域的矢量线要素绘制方法，尤其涉及一种顾及线宽一致性特征的矢量地图方形虚线符号绘制方法。

背景技术

地图是将真实地理世界中各种空间实体信息传递给读图者的有效媒介，地图符号化既是对空间实体的抽象，也是将空间实体进行视觉表达的绘制手段。数字地图制图是通过计算机来进行地图的可视化表达，点、线和面要素的符号化绘制是生成数字地图的基础。其中，矢量线要素的符号化表达既要考虑线宽、颜色、线型等绘图要素，还需要考虑随着地图缩放线宽是否保持一致的问题。方形虚线符号是表达矢量线要素的重要方法，通常用于表达铁路线、领海基线、未确定边界等。根据地图配图的需求，方形虚线地图符号在诸多情况下都需要随地图缩放保持线宽。

针对这种线宽一致性的绘图需求，现有的研究主要有两大类：

(1)基于纯软件绘制的方法。例如，论文《面向GIS的矢量线状符号填充算法研究》(张海泉,潘懋,吴焕萍,等.《地理与地理信息科学》,2004,20(04):11-14.)，给出了一种基于GDI+的纯软件绘制方法来表达矢量线要素。该种方法严格依赖于矢量线的屏幕坐标和屏幕宽度，当地图缩放时，需要不断的进行重新计算和绘制，严重降低了地图显示的效率。针对此效率低下的问题，后续的研究大多数集中于如何用硬件加速的方法，即GPU绘制的方法来实现矢量线符号的绘制。

(2)基于GPU绘制的方法。此类方法主要是研究如何通过GPU着色器(Shader)来绘制线符号。例如，论文《A function-based linear map symbol building and renderingmethod using shader language》(Yue SS,Yang JS,Chen M,等.《International Journalof Geographical Information Science》，2016,30(2):143-16)提出了一种利用OpenGLShader来构建线符号并用GPU绘制。该方法依赖于一个颜色二维数组来表达符号结构，虽然能够通过数组标记来表达和绘制虚线符号，但整个方法中只是针对固定地图比例尺的绘制，对于线宽的一致性问题完全没有涉及。而论文《一种基于屏幕的三维地图线状符号渲染方法》(刘君妍，陈雅茜，高亦远，等.《地球信息科学学报》，2018,20(8):1047-1054)提出了将矢量线要素映射成为纹理数据，通过纹理数据来绘制线；该方法主要是针对线要素与地形的融合，线要素的线宽处理明确是随着地图缩放而变化的。此外名称为“基于GPU的地图线形符号绘制方法及***”(申请号201310125110.6)的中国专利，名为“基于GPU的宽度渐变线型地图符号绘制方法”(申请号201610015703.0)的中国专利，都是在分析GPU绘制矢量线要素难点的基础上，提出了用Shader着色器来绘制线符号的方法。但是这些方法都是针对固定比例尺的线符号绘制，当地图比例尺变化时，为了保持线宽不变，前述方法都需要重新三角化和绘制。反复不断的三角化对绘制效率具有明显的影响，由此造成了地图的缩放交互仍然存在效率瓶颈。

因而，如何同时满足(1)随着地图比例尺缩放保持线宽一致，(2)矢量线要素不需要重新进行三角化，是数字地图制图领域一个亟待解决的问题。对于最经常使用的方形虚线符号而言，这种顾及线宽一致性的绘制需求在本质上依赖于CPU和GPU的协调工作。GPU的绘制是依赖于三角面片的，当一个矢量线被三角化为方形虚线符号后，GPU只能拿到三角形坐标和颜色，整个绘制是无状态的。如何将方形虚线符号的地图状态传递给GPU成为解决顾及线宽一致性绘制方形虚线符号的关键。

发明内容

本发明的目的是在GPU绘制环境下，针对方形虚线地图符号的线宽一致性绘制需求，提供一种能够在地图缩放过程中不改变线宽、不重新三角化的方形虚线地图符号绘制工艺。

技术方案：一种顾及线宽一致性的矢量地图方形虚线符号绘制方法，包括下列步骤：

步骤1、对矢量线要素的坐标序列进行初始化，设置以像素为单位的屏幕线宽，设置方形虚线的虚实比例值；

步骤2、构建矢量线要素的世界坐标系和屏幕坐标系之间的转换规则，通过此规则将屏幕线宽转换为世界坐标系的线宽，根据此世界坐标系线宽：顺序遍历矢量线要素的每个坐标，以半个世界坐标系的线宽为单位，垂直于线走向分别向两侧进行扩张，连接相邻的扩张点生成三角形集合，其中每个三角形的每个顶点都具有沿线走向的U值和垂直线走向的V值，并且还具有一个缩放方向向量D值；

步骤3、组装GPU着色器的状态协议：该状态协议包括两个层次，一个是顶点层，包括U值、V值、缩放方向D值，另一个是结构层，包括最大线宽M值、地图比例缩放S值、实部比例P1值、虚部比例P2值；

步骤4、在CPU中获取初始的世界坐标系线宽、地图比例尺缩放系数、虚实比例值，并分别按照状态协议中的M值、S值、P1值、P2值传递到GPU中，由此在GPU着色器中解析状态协议中的U、V、D、M、S、P1、P2值，动态计算当前线宽在屏幕上需要绘制的像素位置，并完成绘制。

进一步的，所述步骤2中沿线走向的U值是根据每个顶点沿着线要素的距离计算，其中起点扩张出来的两个顶点U值均为0，终点扩张出来的两个顶点U值均为整条线的长度除以线宽，其他点扩张出来的两个顶点U值均根据长度累积值除以线宽获得；垂直线走向的V值是根据每个顶点向两侧扩张的距离计算，线要素上的每个点，其扩张出来的两个顶点位于线走向左侧的V为1，位于右侧的0；顶点缩放方向的D值是一个二维向量，指明了扩张出来的每个顶点与原始点之间的连接方向，并且不是一个归一化的向量，D值二维向量的长度是根据扩张顶点与原始点之间的距离计算获得的。

进一步的，所述步骤3中最大线宽M值是根据步骤2中计算获得的世界坐标系的线宽；地图比例缩放S值是随着地图比例尺缩放动态更新的比例值；实部比例P1值可以是任意大于等于1.0的浮点数，代表方形虚线符号中需要填充颜色的长度，在传递给状态协议之前需要进行线宽映射，即将原始的P1值与线宽相除得到新的P1值；虚部比例P2值可以是任意大于0.0的浮点数，代表方形虚线符号中不需要填充颜色(即间隔)的长度，在传递给状态协议之前需要进行线宽映射，即将原始的P2值与线宽相除得到新的P2值；状态协议按照U值、V值、D值、M值、S值、P1值、P2值进行组织，其中U值和V值是一维浮点数、D值是二维浮点数、M值和S值是一维浮点数、P1值和P2值是一维浮点数，将U值、V值、D值按序组织成四维向量，将M值、S值、P1值、P2值按序组织成四维向量，最终形成的状态协议由两个向量构成：U-V-D四维向量和M-S-P1-P2四维向量，顾及线宽一致性的矢量地图方形虚线符号绘制必须严格遵从此协议，在CPU和GPU之间进行消息传递时也必须严格依照此协议进行消息组装。

和现有技术相比，本发明具有如下显著进步：本发明可以充分利用GPU硬件加速的绘制性能，突破GPU仅依靠坐标和颜色的无状态绘制局限，通过定制合理的状态协议，将矢量线的走向信息、线宽信息、虚实比例信息、地图缩放比例信息在CPU和GPU之间传递。通过两个四维向量的组合，能够有效的组装状态协议；按照U-V-D-M-S-P1-P2的结构能够方便的实现状态内容的获取，从而在GPU的着色器中进行保持线宽的绘制。相比于传统的方法，本发明能够在不重复三角化的基础上，定制清晰的状态协议，实现随着地图比例尺的缩放进行线宽一致的方形虚线绘制，从而提升地图的绘制效率和效果。

附图说明

图1是本发明方法的流程示意图；

图2是状态协议的基本结构；

图3是U值、V值和D值的说明，以及基于状态协议的矢量线要素三角化方法说明；(a)是对矢量线要素的U值、V值、D值的示例，(b)是矢量线要素三角化成三角形序列的示意；

图4是矢量线要素三角化后各顶点的UV值示例；

图5是M值、P1值、P2值的说明，以矢量线要素初始线宽的绘制为例；

图6是M值、S值、P1值、P2值的说明，以矢量线要素缩放后保持线宽的绘制为例，其中灰色填充部分为真实绘制的线宽，外层虚线部分为原始不处理前的绘制区域(即一开始三角化的结果)；

图7是利用本发明的方法绘制结果示意图，(a)所示为状态协议中虚实比例为[1.0:1.0]地图比例尺较小的线宽效果，(b)所示为状态协议中虚实比例为[1.0:1.0]地图比例尺变大后线宽的效果仍然一致，(c)所示为状态协议中虚实比例为[2.0:3.0]地图比例尺较小的线宽效果，(d)所示为状态协议中虚实比例为[2.0:3.0]地图比例尺变大后线宽的效果仍然一致。

具体实施方式

下面结合附图和具体实施例对本发明的技术方案做进一步说明。

图1为本发明方法的流程图。该方法包括以下步骤：首先对待绘制的矢量线要素进行初始化，指定需要在屏幕上呈现的线宽(像素单位)，设置方形虚线的虚实比例值，并确定矢量线要素的世界坐标系和屏幕坐标系之间的转换规则。接着，按照坐标转换规则将屏幕线宽转换为世界坐标系的线宽，将虚实比例值进行线宽映射，并在世界坐标系下对矢量线要素进行三角化操作。然后，按照所构建的顾及线宽一致性的GPU状态协议，在CPU中获取初始最大线宽、地图缩放比例系数、虚实比例值等信息并利用状态协议传递到GPU中。最终，在GPU中获取初始最大线宽和缩放系数等状态后进行保持线宽一致的方形虚线符号绘制。

在GPU绘制环境下，三角化是将矢量线要素进行方形虚线符号化绘制的关键，但是仅靠三角化并不能满足方形虚线符号的绘制需求。GPU着色器获取三角形序列，通过三角形的坐标和颜色来在屏幕上呈现待绘制对象。本质上，GPU的绘制流水线是一种无状态的过程，矢量线的线宽相关信息无法从CPU直接传递给GPU。而对于方形虚线符号而言，其需要沿着线的走向、按照确定的线宽进行绘制。基于此，本发明设计了一种面向方形虚线符号保持线宽一致性绘制的状态协议。CPU中需要按照此状态协议组织信息，GPU中需要按照此状态协议解析信息。

状态协议包含两个层面，一个是顶点层的U值，V值和D值，另一个是结构层的M值、S值、P1值和P2值，如图2所示；U值、V值、D值、M值、S值、P1值和P2值的内涵和计算方法在后续三个段落阐述。U值和V值分别是两个浮点型float数值，D值是一个二维向量，D值具有X和Y方向的分量D(x)和D(y)。M值和S值是两个浮点型的数值，P1值和P2值也是两个浮点型的数值。U-V-D三个值被组装成一个四维向量，M、S、P1、P2四个值被组装成一个四维向量。因此，面向方形虚线符号的线宽一致性绘制，CPU和GPU之间的信息传递就依赖于这两个四维向量的组合。之所以将U-V-D组装成一个四维向量、将M-S-P1-P2组装成一个四维向量，是为了在CPU-GPU信息传递过程中尽可能减少信息交互的次数，批量进行传递。

对于U-V-D值而言，在图3和图4中给出了相关的计算方法。U-V-D值是附加在每个扩展顶点上的。扩展顶点是指依据矢量线要素上的任一点，按照垂直与线走向的方向分别向两边进行扩张后得到的顶点；其中，扩张距离是线宽的一半(世界坐标系的线宽)。扩张顶点可以分为左侧点和右侧点(根据沿线走向来确定)。在图3(a)中给出线要素A-B-C的顶点扩张说明。以线段A-B为例，线段A-B扩张后形成左侧点A、右侧点A、左侧点B和右侧点B，这四个点虽然可以构造成两个三角形，但线段AB的结构特征被破坏了(A-B线段受到了拐角的影响)。因而，本发明在B点处分别对应于前后两个线段构造了垂向的左侧点B1和左侧点B2，用于区分线段区域和拐角区域。如图3(b)所示，拐角处分别构造左右两个三角形，线段A-B构造两个三角形，线段B-C构造两个三角形。如图4所示，所有右侧顶点的V值都是0，所有左侧顶点的V值都是1。任一个扩张顶点的U值都是通过计算其所关联的点与起始点之间的距离得到。起始点的两个扩张顶点的U值都是0，终止点的两个扩张顶点的U值都是整个矢量线的长度除以线宽。起始点与终止点中间的任一点都是按照累积距离除以线宽得到。D值是表明扩张顶点方向的二维向量，该向量是依据扩张顶点与关联点获得。需要注意的是，D向量不能进行归一化，其向量长度就是扩张顶点与关联点之间的长度。

对于M-S值而言，M值是在初始地图比例尺下显示的最大线宽的值，S值是相对于初始地图比例尺的缩放系数。如图5和图6所示，在GPU中解析状态协议里面的M-S值，利用以下公式重新计算任一扩张顶点的坐标：X＝D(x)*M*(1.0-S)/2.0，Y＝D(y)*M*(1.0-S)/2.0。D(x)和D(y)分别表示顶点方向向量的X分量和Y分量。在GPU中结合D值、M值和S值，动态计算每个点真正需要绘制的位置。在图6中外部虚线部分为初始比例尺下需要绘制的区域，中间灰色填充部分为缩放比例尺之后需要绘制的区域。可以看出初始区域大，大比例尺区域小；因为在大比例尺下，矢量线要素的扩张顶点直接转换为屏幕坐标后会比初始比例尺大，所以通过D向量和M-S值来控制扩张顶点向内收缩，从而达到一致线宽的绘制效果。

对于P1-P2值而言，在状态协议中是将虚实比例值进行线宽映射的，因而在GPU中利用公式U/(P1+P2)可以统一的获取当前像素是否在实部，在实部则绘制，在虚部则不绘制。如图6所示，随着地图比例尺的缩放，在GPU中进行U值、P1值和P2值与缩放系数S值的相乘，从而动态获取当前比例尺下U值、P1值和P2值，由此实现线宽一致的方形虚线绘制。

需要特别说明的是，基于本发明的绘制工艺，在线要素的拐角处能够识别是否处于实部区域；如果刚好处于实部与虚部的交界处，将会进行绘制，从而能够保留拐角特征，如图5中的拐角处所示。

在图7中给出了两个不同地图比例尺，方形虚线符号的绘制效果。左侧是小比例尺，(a)所示为状态协议中虚实比例为[1.0:1.0]地图比例尺较小的线宽效果，(c)所示为状态协议中虚实比例为[2.0:3.0]地图比例尺较小的线宽效果；右侧是大比例尺，(b)所示为状态协议中虚实比例为[1.0:1.0]地图比例尺变大后线宽的效果仍然一致，(d)所示为状态协议中虚实比例为[2.0:3.0]地图比例尺变大后线宽的效果仍然一致。可以看出小比例尺和大比例尺的绘制效果在线宽方面是一致的。

实施例：

相比于传统计算机可视化领域中的方形虚线绘制，地图绘制中的方形虚线符号需要考虑地图比例尺的变化。GPU绘制环境下，仅通过三角形集合的坐标和颜色来控制绘制流程，无法满足保持线宽一致性的需求。所以本发明设计了一种CPU-GPU状态协议，通过此协议将线宽相关的信息打包传递给GPU，在GPU中能够有序解析这些信息，从而达到顾及线宽一致性的矢量地图方形虚线符号绘制。

以下是利用OpenGL Shader给出了本发明的实施例：

(1)GPU状态协议结构

(2)CPU状态协议结构

Claims (2)

1.一种顾及线宽一致性的矢量地图方形虚线符号绘制方法，其特征在于，包括下列步骤：

步骤1、对矢量线要素的坐标序列进行初始化，设置以像素为单位的屏幕线宽，设置方形虚线的虚实比例值；

步骤2、构建矢量线要素的世界坐标系和屏幕坐标系之间的转换规则，通过此规则将屏幕线宽转换为世界坐标系的线宽，根据此世界坐标系线宽：顺序遍历矢量线要素的每个坐标，以半个世界坐标系的线宽为单位，垂直于线走向分别向两侧进行扩张，连接相邻的扩张点生成三角形集合，其中每个三角形的每个顶点都具有沿线走向的U值和垂直线走向的V值，并且还具有一个缩放方向D值；

步骤3、组装GPU着色器的状态协议：该状态协议包括两个层次，一个是顶点层，包括U值、V值、缩放方向D值，另一个是结构层，包括最大线宽M值、地图比例缩放S值、实部比例P1值、虚部比例P2值；

步骤4、在CPU中获取初始的世界坐标系线宽、地图比例尺缩放系数、虚实比例值，并分别按照状态协议中的M值、S值、P1值、P2值传递到GPU中，由此在GPU着色器中解析状态协议中的U值、V值、D值、M值、S值、P1值、P2值，动态计算当前线宽在屏幕上需要绘制的像素位置，并完成绘制；

所述步骤2中，沿线走向的U值根据每个顶点沿着线要素的距离计算获得，其中起点扩张出来的两个顶点U值均为0，终点扩张出来的两个顶点U值均为整条线的长度除以线宽，其它点扩张出来的两个顶点U值均为长度累积值除以线宽；垂直线走向的V值根据每个顶点向两侧扩张的距离计算获得，线要素上的每个点，其扩张出来的两个顶点位于线走向左侧的V值为1，位于右侧的V值为0；缩放方向D值是一个二维向量，二维向量的长度根据扩张顶点与原始点之间的距离计算获得；

所述步骤3中，最大线宽M值是世界坐标系的线宽；所述步骤3中，地图比例缩放S值是随着地图比例尺缩放动态更新的比例值；所述步骤3中，实部比例P1值是任意大于等于1.0的浮点数，代表方形虚线符号中需要填充颜色的长度，在传递给状态协议之前需要进行线宽映射，即将原始的P1值与线宽相除得到新的P1值；所述步骤3中，虚部比例P2值是任意大于0.0的浮点数，代表方形虚线符号中不需要填充颜色的长度，在传递给状态协议之前需要进行线宽映射，即将原始的P2值与线宽相除得到新的P2值。

2.根据权利要求1所述的顾及线宽一致性的矢量地图方形虚线符号绘制方法，其特征是，所述步骤3中，状态协议按照U值、V值、D值、M值、S值、P1值、P2值进行组织，其中U值和V值是一维浮点数、D值是二维浮点数、M值和S值是一维浮点数、P1值和P2值是一维浮点数，将U值、V值、D值按序组织成四维向量，将M值、S值、P1值、P2值按序组织成四维向量，最终形成的状态协议由两个向量构成：U-V-D四维向量和M-S-P1-P2四维向量。