CN114529633B - 一种支持gis线对象和面对象连续lod绘制的方法 - Google Patents

一种支持gis线对象和面对象连续lod绘制的方法 Download PDF

Info

Publication number
CN114529633B
CN114529633B CN202210426059.1A CN202210426059A CN114529633B CN 114529633 B CN114529633 B CN 114529633B CN 202210426059 A CN202210426059 A CN 202210426059A CN 114529633 B CN114529633 B CN 114529633B
Authority
CN
China
Prior art keywords
node
data
line
trapezoidal
nodes
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210426059.1A
Other languages
English (en)
Other versions
CN114529633A (zh
Inventor
吴明光
沈飞
闾国年
陈康正
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Normal University
Original Assignee
Nanjing Normal University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing Normal University filed Critical Nanjing Normal University
Priority to CN202210426059.1A priority Critical patent/CN114529633B/zh
Publication of CN114529633A publication Critical patent/CN114529633A/zh
Application granted granted Critical
Publication of CN114529633B publication Critical patent/CN114529633B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/20Drawing from basic elements, e.g. lines or circles
    • G06T11/206Drawing of charts or graphs
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/22Indexing; Data structures therefor; Storage structures
    • G06F16/2228Indexing structures
    • G06F16/2246Trees, e.g. B+trees
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical information databases
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/50Information retrieval; Database structures therefor; File system structures therefor of still image data
    • G06F16/51Indexing; Data structures therefor; Storage structures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/50Information retrieval; Database structures therefor; File system structures therefor of still image data
    • G06F16/56Information retrieval; Database structures therefor; File system structures therefor of still image data having vectorial format
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Remote Sensing (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了一种支持GIS线对象和面对象连续LOD绘制的方法,包括以点区域四叉树为原型,初始化空间索引,根据空间索引的层次性重载索引节点和数据的访问函数,通过范围相交查询构建层次编码;单独构建面要素的梯形格网,改进算法,迭代计算各节点的权重并作为增量数据录入数据库,构造包含权重值的要素数据结构;根据节点包围盒与可视化查询窗口的拓扑关系、像素比例尺等信息,从数据库中查询增量数据,对节点进行化简或累积简化;对于线,从线段主体和拐角样式两部分,剖分渲染线要素;对于面,以骨干梯形将梯形格网拆分为梯形链,简化并渲染梯形格网结构下的面要素。本发明可以实现矢量地图中线要素和面要素连续LOD的高性能绘制。

Description

一种支持GIS线对象和面对象连续LOD绘制的方法
技术领域
本发明涉及地理信息系统(GIS)、计算机辅助设计(CAD)、计算机地图制图(CAC),导航地图、电子地图、手机地图领域,特别涉及一种GIS线对象和面对象连续LOD绘制的方法。
背景技术
地图可视化一直是现代地图学的核心,地图可视化方法的性能和质量对高交互的制图和分析具有重要意义。大数据时代的到来,使得当前可视化方法难以应对矢量数据的海量性和复杂性带来的地图高交互性和数据查询与交互的挑战,矢量数据渲染所涉及的多尺度空间索引技术、LOD技术等方面的研究较为缺乏。
多尺度空间索引技术根据可视化尺度来获取用于快速显示的不同尺度的数据。常见的多尺度空间索引方法以Reactive-tree为代表,包括R树、R*树、SR树等。这些方法能够优先获取尺度较大、重要性较高的矢量数据,但破坏了空间索引的平衡性,改变了原有空间索引的存储方法,降低了数据查询效率。同时,当前多尺度索引多为面向数据快速获取,并未考虑数据调度成本与数据LOD成本。LOD技术分为静态LOD与连续LOD。静态LOD方法多围绕生成不同精度数据的质量和效率展开。常见的静态LOD方法包括顶点聚类算法、浮动栅格聚类方法、ICE等。静态LOD在可视化性能上有着优越表现,但需要较长的预处理时间,同时导致相邻比例尺下,不同层次细节之间难以平滑过渡。连续LOD的研究主要围绕矢量数据简化和累积简化。常见的简化算法多为从原始数据进行简化得到最终的简化数据,并未考虑简化后数据重复利用。目前还未形成一个通用的、高质量的GIS线对象和面对象连续LOD绘制方法。
发明内容
本发明的目的在于解决当前大规模矢量数据带来的地图可视化效率低、空间索引未考虑矢量数据LOD调度问题,提供一种基于GPU的要素连续LOD绘制方法,实现线要素和面要素连续LOD高性能渲染。
本发明的总体思路是:针对基于数据库的外部空间索引查询方法的数据查询与调度效率低问题,设计层次编码;针对矢量要素增量数据难以统一管理问题,设计矢量数据存储与组织结构;针对矢量数据LOD快速获取问题,设计基于层次编码索引结构的数据查询与调度方法;针对数据渲染效率低问题,设计基于GPU环境和要素样式的数据剖分与LOD算法。
本发明一种支持GIS线对象和面对象连续LOD绘制方法,包括以下几个步骤:
步骤1,以点区域四叉树初始化空间索引,重载索引节点和数据的访问函数,通过范围相交查询构建索引节点及其数据的层次编码;
步骤2,先单独构建面要素的梯形格网结构,再改进DP算法,迭代计算要素各节点的权重值,并将权重值作为增量信息录入数据库,再构建包含权重值的要素数据结构;
步骤3,根据节点包围盒与可视化查询窗口的相交关系,以及节点数据的累积简化状态,从数据库中查询增量数据,对节点进行化简或累积简化;
步骤4,在GPU中,去除相邻线段重复部分,绘制线段主体,再针对平角、斜角、圆角三种拐角样式,剖分绘制线段拐角;以骨干梯形将梯形格网拆分成梯形链,结合数据库中的增量信息,对细节梯形进行合并简化与累积简化,最终渲染面要素。
上述步骤1具体包括:(1)根据OGR图层的要素总量,计算索引层数,使得索引层数相对平衡,保证四叉树的高度适中;(2)以点区域四叉树为原型,生成基础的空间索引,并向树中插入要素的ID与范围;(3)对于点区域四叉树的节点访问和数据访问函数,进行自定义的重载,根据节点的层数生成相应的节点编码,再基于节点的编码,给节点下的各个数据赋予编码;(4)经过(1)、(2)和(3)可以得到反映空间索引层次性的层次编码索引结构。以层次编码为序的数据存储顺序与空间分布相关,即具有共同祖先节点的数据在线性空间中被存储到一起,能够大幅提高数据的查询和调度效率。
上述步骤2具体包括:(1)先单独对面要素构建梯形格网,逐一获取面要素外环和各内环上的点坐标与点个数,根据点的坐标数据与点的总数,利用水平扫描线,逐个扫描行将面要素切分成梯形面片,形成梯形连通图并基于梯形连通图,将梯形面片划分为骨干梯形与细节梯形两类;遍历梯形链,以梯形格网平行边为节点,梯形为节点连接线,构建平行边多段线图。(2)获取要素第一个和最后一个节点组成节点对,根据节点对的坐标数据,计算位于节点对之间各节点的权重值,找到权重值最大的节点,以该节点与节点对的两个节点分别形成新的两组节点对,再次计算节点对之间各顶点的权重值,以此不断迭代计算要素各节点的权重。(3)将节点权重值作为要素的增量信息,录入数据库。最终,构造包含权重值的要素数据结构。(4)经过(1)、(2)和(3)可以得到针对LOD、更完善的要素数据结构体系。
上述步骤3具体包括:(1)判断可视化窗口与空间索引范围是否相交,若相交,则取相交范围为实际查询范围,否则不进行渲染更新;再判断实际查询范围、比例尺、选择ID数,若任一项与前一次不同,则更新本次查询范围、比例尺、选择ID数,否则不进行渲染更新;(2)根据节点的包围盒与像素比例尺、可视化查询窗口与像素比例尺,结合节点数据的累积简化状态,对节点进行几何数据查询;(3)一方面,当内存中数据已满足可视化分辨率细节需求时,删除冗余的节点数据;另一方面,根据从数据库中查询得到的增量数据,对未满足可视化需求的节点进行展开;(4)经过(1)、(2)和(3)可以得到面向LOD可视化的渲染数据集。
上述步骤4具体包括:在GPU环境下,绘制线要素时,对于线的主体绘制,通过平行线推移去除相邻线段的重复部分,减少曲线的重叠绘制;对于线的拐角样式绘制,将线要素拐角的连接样式分为平角、斜角、圆角三种样式,连接推移顶点绘制平角样式,以交点将拐角缺口剖分成两三角形绘制斜角样式,采用顶点旋转方式绘制圆角样式。绘制面要素时,基于步骤2的梯形格网结构,根据梯形连通图,以骨干梯形将梯形格网拆分成多段由多个细节梯形组成的梯形链,同时从数据库中获取的增量信息。一方面,对相邻细节梯形基于权重进行合并简化;另一方面,以原始梯形和增量信息中的顶点,构成新的细节梯形进行累积简化。最终对简化结果进行渲染。
有益效果:本发明能够支持GIS线对象和面对象进行高性能的连续LOD绘制,满足地图可视化高交互的快速响应需求,同时直接保留原始几何数据的完整性,便于空间数据的查询与空间操作。与Reactive-tree相比,本方法的空间索引避免破环原有的空间索引结构,实现了对数据库中空间数据的范围查询,同时提高空间索引查询效率。与常见的静态、连续LOD方法相比,本方法顾及可视化中矢量数据量与数据LOD的关系,面向数据LOD渲染设计空间索引结构和数据调度方法,最终在GPU环境下完成线要素和面要素的剖分与渲染。
附图说明
图1本发明方法的流程图;
图2层次编码示意图;
图3数据组织结构与逻辑结果示意图,(a)要素组织结构,(b)数据与节点逻辑存储结构;
图4扫描剖分示意图;
图5数据查询流程图;
图6推移线重叠部分相对坐标示意图;
图7线要素剖分渲染示意图;
图8梯形合并示意图,(a)合并前的相邻梯形,(b)合并后梯形顶点的变化情况,(c)合并前的节点连接线,(d)合并后节点连接线的变化情况;
图9梯形增量构建示意图,(a)原始梯形,(b)根据增量信息简化得到的新梯形。
具体实施方式
下面结合附图和实施例对本发明做进一步详细说明。
本发明将以点区域四叉树为空间索引原型,构建考虑空间索引层次性的层次编码索引结构;单独构建面要素的梯形格网结构,并改进DP算法,计算要素节点权重并构建针对LOD的要素数据结构;根据可视化窗口和节点的包围盒与像素比例尺,结合节点的累积简化状态,对节点进行几何数据查询;最终在GPU环境下,基于线样式分线段主体与线段拐角两部分,剖分渲染线要素,根据骨干梯形将梯形格网拆分为梯形链,结合从数据库中获取的增量信息,对梯形格网结构下的面要素进行简化与渲染。
实施例1
如图1所示,我国某范围内路网地图数据连续LOD渲染包含以下步骤:
步骤1,以点区域四叉树初始化空间索引,重载索引节点和数据的访问函数,通过范围相交查询构建索引节点及其数据的层次编码;
步骤2,先单独构建面要素的梯形格网结构,再改进DP算法,迭代计算要素各节点的权重值,并将权重值作为增量信息录入数据库,再构建包含权重值的要素数据结构;
步骤3,根据节点包围盒与可视化查询窗口的相交关系,以及节点数据的累积简化状态,从数据库中查询增量数据,对节点进行化简或累积简化;
步骤4,在GPU中,去除相邻线段重复部分,绘制线段主体,再针对平角、斜角、圆角三种拐角样式,剖分绘制线段拐角;以骨干梯形将梯形格网拆分成梯形链,结合数据库中的增量信息,对细节梯形进行合并简化与累积简化,最终渲染面要素。
更具体的实施步骤如下:
四叉树初始化与层次编码索引构建
本发明改进点区域四叉树划分方法从上往下构建,但需要统计每个索引节点包含的点数据量(其他可视化参数)信息。本发明索引构建算法具体描述如下:
1)提取原始矢量数据集合D(D1、D2…)的地理范围,并以数据集包围盒中心点将地理范围拓展成正方形,基于数据的地理范围计算LevelMax等信息;
2)以当前数据集的统计信息(顶点数量)、节点编码,构建索引节点;
3)将地理范围均分为四子块空间,依据子块地理范围对数据集合进行分割;
4)遍历子块,判断子块节点深度是否等于LevelMax,若是,则将子块构建节点后,计算节点尺度S,并将子块中所有数据和节点信息存储至数据库中;
5)若否,则判断子块点数集合中是否存在点数量小于点数量阈值的Di,若是,则将子块构建节点后,计算节点尺度S,并将子块中Di和节点信息存储至数据库中;否则,以当前子块和子块中数据集合进行步骤2)。
与传统四叉树的构建算法相比,本发明提出的构建算法针对的是多要素集合,能够实现对多要素集进行索引构建。同时,本发明结合数据地理范围和可视化信息,能够自适应计算索引最大层级,避免树索引深度过高带来索引检索效率低的问题。本发明索引构建时,需要统计点的数量(可视化属性),同时需要计算节点的尺度S。节点的尺度S计算如下,假设Len为数据的拓展成方形的地理范围的宽或高,Z为节点在索引中的层级,节点尺度S的计算公式如下:
Figure 543219DEST_PATH_IMAGE001
本发明基于空间索引的层次性对索引节点和要素进行编码,使编码具有以下特征:一方面,父节点的编码小于孩子节点的编码;另一方面,同一层两节点A、B,若A的节点编码小于B的节点编码,则节点A的所有孩子节点编码同样小于B的节点编码。编码结果如图2所示。
图3为以编码为序的数据存储以及数据与节点逻辑存储示意图,图3的 (a)为要素 在数据库中的组织结构示意图,CODE字段为本发明的编码。从图3的 (b)可以看出以层次编 码为序的数据存储顺序与空间分布相关,即具有共同祖先节点的数据在线性空间中被存储 到一起,如要素3与6在存储上临近。同时,要素存储顺序符合基于索引深度查询,要素访问 顺序。基于编码结构能够实现对数据的范围查询,如当数据查询窗口包含索引节点2的地理 范围时,只需查询编码
Figure 802162DEST_PATH_IMAGE002
,即利用逻辑符号可快速查询、调度数据。
层次编码具有良好的拓展性,编码是基于树索引结构的层次性建立,并不会影响索引的构建方式和索引中数据存储位置,因此能应用于各类空间索引结构上。
梯形格网构建与节点权重计算
梯形格网是简化与剖分渲染一体化数据结构,能够对有效支撑面数据的连续LOD。本发明基于以梯形格网作为面要素的LOD结构,基于扫描线算法对面要素进行梯形剖分,对多边形原始顶点的纵坐标进行排序,依据纵坐标值生成x方向水平的扫描线,将多边形划分成按行分布的梯形面片,其空间位置由所在梯形四个顶点坐标值所决定。梯形剖分算法如下:
Step1:以排序后多边形原始节点y值构建扫描线数组。对多边形顶点进行y方向从大到小的唯一值(消除相同的y值)快速排序,依据排序结果建立扫描线数组S 1S 2、...、S n
Step2:以边节点的最小y值构建边链表。依据节点顺序构建边数组,以边数组的中每一条边的最小y值为Key值,按照从小到大的顺利快速排序,插入至边链表;
Step3:从扫描线S1开始,计算与当前扫描线Scur相交的所有边,对边与扫描线的交点依据x值进行快速排序,x值相同再依据边的斜率进行排序;
Step4:取下一条扫描线Snext,再次进行Step3,把两次计算的结果组成水平方向连续的梯形组;
Step5:接着更新ScurSnext进行下一个循环,直至Scur=Sn,扫描结束。
如图4,逐个提取梯形组中的梯形,即可实现多边形的梯形剖分。但剖分多边形的顶端和尾端可能会出现三角形,本发明通过存储两个相同的上顶点进行三角形的存储,同时也通过梯形具备两个相同的上顶点完成对梯形是否是三角形的判断。
几何数据节点权重算法决定了节点简化顺序和图形简化质量,本发明对于要素的LOD处理,基于DP算法设计本发明的数据节点权重算法。假设要素首节点Vs,要素尾部节点Ve,初始两端节点权重为无限大,则算法描述如下:
1)连接要素首尾节点形成直线VsVe
2)获取要素上距离直线VsVe距离最大的点Vi
3)判断Hi与节点两端顶点权重的大小,若小于,则以Hi作为点Vi的权重,若大于,则以两端节点中的权重最小值min(Hs,He)作为该点Vi的权重;
4)将要素以Vi节点为分割点形成新的两条要素,即Vs和Vi及其中间节点组成一段要素,Vi和Ve及其中间节点组成一段要素;
5)重复1)、2)、3)、4)步骤,直至要素中间不存在点,迭代结束。
基于本发明的算法,任意要素ViVj之间所有节点权重必须不高于Vi、Vj点权重,若高于,取两点权重最小值。如此通过权重调整,能够保证节点权重的单调性符合DP算法节点简化顺序,即优先简化的节点权重更小。
最终,以ID号、坐标数据和节点权重值作为成员变量,构建新的要素类,该类是针对LOD、更加完备的数据结构体系。
几何数据查询与LOD调度
本发明基于层次编码结构,设计针对数据库的空间数据范围查询方法。本发明的索引为可视化与数据调度一体化索引,结合索引查询包含求精和过滤两个过程,实现对数据的增量查询和化简查询,即基于查询需要,查询增量数据并对几何数据进行累积简化,也需要对其他数据进行简化。假设节点的像素比例尺为SN,当前地图可视化的像素比例尺为S,可视化查询窗口BV,节点包围盒为BN,本发明索引查询算法如图5所示。
1)判断子树根节点BN与BV是否相交,若不相交,则说明数据不在可视化范围内,对数据进行简化,并将节点所包含数据移除渲染数据集;
2)若相交,判断节点数据中数据是否已进行累积简化,若没有则以节点编码查询增量数据,进行数据累积简化;
3)若已进行累积简化,则判断当前节点SN是否大于S,若大于则说明内存中数据已满足可视化分辨率细节需求,将数据加载至渲染数据集;
4)若SN不大于S,则判断是否BV包含BN,若包含则说明当前节点所有数据都需进行渲染,基于节点编码进行增量信息范围查询,将累积简化后数据加入渲染数据集;若否,则继续遍历节点子树。
本发明方法查询的结果为有要素编码和索引节点编码构成的编码集。索引节点编码代表以该节点为祖先节点的所有要素都被查询窗口V包含。
在查询算法方面,本发明查询考虑节点与查询窗口的包含关系,当查询窗口包含节点时,则不再进行遍历子节点,有效避免了子节点的深度遍历。在外部索引查询结果集方面,本发明外部索引查询结构为编码集,而不是要素ID等信息。在基于结果集对数据库中数据进行查询调度方面,本发明以要素的编码进行离散查询,以索引节点编码进行范围查询。从而以索引节点编码进行范围查询,有效减少数据查询时的数据库扫描次数,降低数据库内部数据查询调度成本。同时本发明方法仅对空间索引进行编码和改变空间索引查询方法,并不会影响空间索引的构建方式和构建结果,因此本发明方法能拓展至其他空间索引结构上,具有适用性和拓展性。
基于要素样式与数据结构的剖分渲染
本发明针对线要素拐角的连接样式,采用平行线推移算法,在GPU中进行线要素的剖分绘制。以线段为单元,线要素的绘制分为线段主体(矩形区)绘制和拐角样式绘制。从而避免线的压盖绘制。
针对线段主体绘制需要去除与相邻线段重复部分。如图6所示。若以P2P3为X方向、 P2V3为Y方向,建立坐标系。P2P3长度为L,则矩形V3U3U4V4的坐标范围R2
Figure 87650DEST_PATH_IMAGE003
矩形V1U1U2V2的坐标范围R1
公共区域的范围
Figure 90241DEST_PATH_IMAGE004
。绘制矩形V1U1U2V2时,需要去除与矩形V3U3U4V4的 公共部分,即
Figure 879206DEST_PATH_IMAGE005
区域。在几何着色器阶段,计算顶点V1、U1、U2、V2 以 V3P2P3为坐标轴下的相对坐标;光栅化阶段会内插计算光栅化的像素点以V3P2P3为坐标轴下 的相对坐标;片段着色器阶段,判断是否属于R2,若属于则不处理当前像素点,从而避免当 前像素点的处理与绘制。
针对线的拐角样式剖分与绘制,如图7所示,从可视化的角度出发,线符号的连接 样式分为斜角、平角和圆角连接。①平角样式,连接推移顶点V2、P2、V3形成三角形V2 P2 V3即 完成平角样式的剖分,图7中(c)为基于剖分结果,线段V1 V2部分渲染结果示意图。②斜角样 式,计算线段V1 V2、V3 V4的交点,以交点将拐角缺口剖分成两三角形,图7中(g)为线段渲染 结果。③针对线的圆角样式,拟圆化处理完成拐角的平滑过渡,将圆弧部分均分成多个以线 宽
Figure 574629DEST_PATH_IMAGE006
为腰的等腰三角形。假设折线的宽度为w,为保证圆弧三角化变形较小,本发明定义 三角化变形最大变形长度为d,则三角形顶点角
Figure 347413DEST_PATH_IMAGE007
,采用顶点旋转 方式计算弧线V2 V3中剖分三角顶点坐标。顶点坐标计算公式如下:
Figure 888116DEST_PATH_IMAGE008
将剖分渲染方法与线数据LOD方法相结合,即可实现对矢量线LOD渲染。
本发明对面要素的简化分为合并简化与累计简化,以简化结果作为最终的渲染数据。
合并简化主要针对相邻细节梯形,对骨干梯形不进行简化。本发明以骨干节点对梯形连通图拆分,将图拆分成多段由多个细节梯形节点组成的梯形链。针对拆分后的细节梯形基于权重进行简化,最终简化为一个梯形。简化算法如下:
如图8中(a)、(b)所示,两相邻梯形T1、T2二者通过合并形成新的梯形,新梯形由原始梯形顶点V1V2(T1梯形的上边L1)和V5V6(T2梯形的下边组成)构成。梯形合并从形式上来看,移除V3V4线,产生变形。
若以构成梯形格网平行边为节点,梯形为节点连接线,如图8中(c)、(d)所示,则梯形T1为L1、L2,连接线,T2为L2为L3连接线,若简化L2则L1L3连接构成新的图形。
若从相邻梯形单侧线考虑,以梯形左侧顶点构成曲线V1V3V5,基于DP算法可知,当简化阈值大于顶点V3到V1V5距离时H3,移除顶点V3,同样当简化阈值大于顶点V4到V2V6距离H4时,则移除V4顶点。因此基于梯形格网的DP算法中,只有阈值大于L2节点的权重(max(H3、H4)),才能移除线段V3V4,进行梯形合并操作。
本发明基于以平行边为节点的梯形合并算法,设计了面数据基于细节梯形的DP算法。算法描述如下:
1)对面要素几何数据采用梯形剖分,构建梯形连通图;
2)基于梯形连通图划分梯形类型,基于骨干梯形,将连通图划分为多个由细节梯形构成的梯形链;
3)遍历梯形链,进行DP算法;
4)以梯形平行边为节点,以梯形为连接关系,构建平行边多段线图;
5)获取多段线图中最上最下节点Lt、Lb
6)获取节点上距离Lt、Lb距离最大的点Li和其权重Wi
7)判断Wi与多段线两端节点权重的大小,若小于,则以Wi作为节点Li的权重,若大于,则以两端节点中的权重最小值作为该节点的权重(min(Wt,Wb));
8)将多段线以Li点为分割点形成新的两条多段线,重复5)、6)、7)、8)步骤,直至多段线中间不存在节点,迭代结束,计算出节点(平行边)权重。
本发明的增量信息结构以及基于面数据累积简化算法如下:
面数据的增量信息以平行边为单位,由两个顶点、权重、更新梯形的ID和新生成梯形的ID组成。如图9中(a)所示,原始梯形为T1(V1、V2、V3、V4);如图9中(b)所示,增量信息为线段E(V5、V6;W;T1;T2)。获取增量信息E时,将T1分裂成T1和T2,新的梯形T2的坐标由原始梯形和增量信息中顶点构成,即T1变为(V1、V2、V5、V6)、T2变为(V5、V6、V3、V4),并将以堆栈存储当前增量信息。数据简化时,从堆栈中取出增量信息,将T2从渲染数据中移除,同时更新T1的坐标,移除增量信息。
本发明矢量数据符号化与剖分在GPU中完成,能发挥GPU的计算性能,减少CPU计算负载。

Claims (5)

1.一种支持GIS线对象和面对象连续LOD绘制的方法,其特征在于,包括如下步骤:
步骤1,以点区域四叉树初始化空间索引,重载索引节点和数据的访问函数,通过范围相交查询构建索引节点及其数据的层次编码;
步骤2,先单独构建面要素的梯形格网结构,再找到首尾节点之间权重值最大的节点,取该节点的权重值与首尾节点的权重值中的最小值,作该节点的权重,以此改进DP算法计算节点权重,再以该节点与首尾节点组成新的两组首尾节点,来迭代计算各节点的权重值,最终将权重值作为增量信息录入数据库,并且构建包含权重值的要素数据结构;
步骤3,根据节点包围盒与可视化查询窗口的相交关系,以及节点数据的累积简化状态,从数据库中查询增量数据,对节点进行化简或累积简化;
步骤4,在GPU中,去除相邻线段重复部分,绘制线段主体,再针对平角、斜角、圆角三种拐角样式,剖分绘制线段拐角;以骨干梯形将梯形格网拆分成梯形链,结合数据库中的增量信息,包括被更新的梯形的ID、新生成的梯形的ID,合并简化移除或累积简化生成的两个顶点的坐标与权重,对细节梯形进行合并简化与累积简化,最终渲染面要素。
2.根据权利要求1所述的一种支持GIS线对象和面对象连续LOD绘制的方法,其特征在于,所述步骤1具体包括:
步骤1.1,根据OGR图层的要素总量,计算索引层数,使得索引层数相对平衡,保证四叉树的高度适中;
步骤1.2,以点区域四叉树为原型,生成基础的空间索引,并向树中插入要素的ID与范围;
步骤1.3,对于点区域四叉树的节点访问和数据访问函数,进行自定义的重载,根据节点的层数生成相应的节点编码,再基于节点的编码,给节点下的各个数据赋予编码;
步骤1.4,经过上述步骤得到反映空间索引层次性的层次编码索引结构;以层次编码为序的数据存储顺序与空间分布相关,在线性空间中存储具有共同祖先节点的数据。
3.根据权利要求1所述的一种支持GIS线对象和面对象连续LOD绘制的方法,其特征在于,所述步骤2具体包括:
步骤2.1,单独对面要素构建梯形格网,逐一获取面要素外环和各内环上的点坐标与点个数,根据点的坐标数据与点的总数,利用水平扫描线,逐个扫描行将面要素切分成梯形面片,形成梯形连通图并基于梯形连通图,将梯形面片划分为骨干梯形与细节梯形两类;遍历梯形链,以梯形格网平行边为节点,梯形为节点连接线,构建平行边多段线图;
步骤2.2,获取要素第一个和最后一个节点组成节点对,根据节点对的坐标数据,计算位于节点对之间各节点的权重值,找到权重值最大的节点,并取该节点的权重值与节点对的权重值中的最小值,作为该节点的新权重值,再以该节点与节点对的两个节点分别形成新的两组节点对,再次计算节点对之间各顶点的权重值,以此不断迭代计算要素各节点的权重;
步骤2.3,将节点权重值作为要素的增量信息,录入数据库;构造包含权重值的要素数据结构;
步骤2.4,经过上述步骤得到针对LOD、完善的要素数据结构体系。
4.根据权利要求1所述的一种支持GIS线对象和面对象连续LOD绘制的方法,其特征在于,所述步骤3具体包括:
步骤3.1,判断可视化窗口与空间索引范围是否相交,若相交,则取相交范围为实际查询范围,否则不进行渲染更新;再判断实际查询范围、比例尺、选择的ID数,若任一项与前一次不同,则更新本次查询范围、比例尺、选择的ID数,否则不进行渲染更新;
步骤3.2,根据节点的包围盒与像素比例尺、可视化查询窗口与像素比例尺,结合节点数据的累积简化状态,对节点进行几何数据查询;
步骤3.3,当内存中数据已满足可视化分辨率细节需求时,删除冗余的节点数据;根据从数据库中查询得到的增量数据,对未满足可视化需求的节点进行展开;
步骤3.4,经过上述步骤得到面向LOD可视化的渲染数据集。
5.根据权利要求1所述的一种支持GIS线对象和面对象连续LOD绘制的方法,其特征在于,所述步骤4具体包括:在GPU环境下,绘制线要素时,对于线的主体绘制,通过平行线推移去除相邻线段的重复部分,减少曲线的重叠绘制;对于线的拐角样式绘制,将线要素拐角的连接样式分为平角、斜角、圆角三种样式,连接推移顶点绘制平角样式,以交点将拐角缺口剖分成两三角形绘制斜角样式,采用顶点旋转方式绘制圆角样式;绘制面要素时,基于步骤2的梯形格网结构,根据梯形连通图,以骨干梯形将梯形格网拆分成多段由多个细节梯形组成的梯形链,同时从数据库中获取的增量信息;一方面,对相邻细节梯形基于权重进行合并简化;另一方面,以原始梯形和增量信息中的顶点,构成新的细节梯形进行累积简化;最终对简化结果进行渲染。
CN202210426059.1A 2022-04-22 2022-04-22 一种支持gis线对象和面对象连续lod绘制的方法 Active CN114529633B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210426059.1A CN114529633B (zh) 2022-04-22 2022-04-22 一种支持gis线对象和面对象连续lod绘制的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210426059.1A CN114529633B (zh) 2022-04-22 2022-04-22 一种支持gis线对象和面对象连续lod绘制的方法

Publications (2)

Publication Number Publication Date
CN114529633A CN114529633A (zh) 2022-05-24
CN114529633B true CN114529633B (zh) 2022-07-19

Family

ID=81627783

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210426059.1A Active CN114529633B (zh) 2022-04-22 2022-04-22 一种支持gis线对象和面对象连续lod绘制的方法

Country Status (1)

Country Link
CN (1) CN114529633B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114708140B (zh) * 2022-06-07 2022-09-13 中国电建集团华东勘测设计研究院有限公司 一种基于单元的数值计算模型实时交互剖切可视化方法
CN115048008B (zh) * 2022-06-17 2023-08-15 浙江中控技术股份有限公司 一种hmi画面中对象的可视化方法和设备
CN115145930A (zh) * 2022-09-02 2022-10-04 南京师范大学 基于树状层次索引的gis矢量数据层次编码方法及装置
CN115168682B (zh) * 2022-09-05 2022-12-06 南京师范大学 一种大规模时空点数据lod绘制方法及装置
CN117237503B (zh) * 2023-09-18 2024-02-23 武汉含秋数据科技有限公司 一种地理要素数据加速渲染及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105336003A (zh) * 2015-09-28 2016-02-17 中国人民解放军空军航空大学 结合gpu技术实时流畅绘制出三维地形模型的方法
CN110910505A (zh) * 2019-11-29 2020-03-24 西安建筑科技大学 一种场景模型的加速渲染方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7940279B2 (en) * 2007-03-27 2011-05-10 Utah State University System and method for rendering of texel imagery
CN107564087B (zh) * 2017-09-11 2019-08-06 南京大学 一种基于屏幕的三维线状符号渲染方法
CN113268557B (zh) * 2021-06-10 2023-08-04 中国人民解放军国防科技大学 一种适应显示导向型可视化分析的快速的空间索引方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105336003A (zh) * 2015-09-28 2016-02-17 中国人民解放军空军航空大学 结合gpu技术实时流畅绘制出三维地形模型的方法
CN110910505A (zh) * 2019-11-29 2020-03-24 西安建筑科技大学 一种场景模型的加速渲染方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
An Efficient Visualization Method for Polygonal Data with Dynamic Simplification;Mingguang Wu等;《International Journal of Geo-Information》;20180402;全文 *
一种空间分布模式驱动的空间索引;吴明光;《测绘学报》;20150131;第44卷(第1期);全文 *
基于Hilbert曲线的STR索引改进算法;戴晶等;《武汉大学学报·信息科学版》;20140731;第39卷(第7期);全文 *

Also Published As

Publication number Publication date
CN114529633A (zh) 2022-05-24

Similar Documents

Publication Publication Date Title
CN114529633B (zh) 一种支持gis线对象和面对象连续lod绘制的方法
Garland et al. Fast polygonal approximation of terrains and height fields
Heckbert et al. Survey of polygonal surface simplification algorithms
De Floriani et al. Multiresolution models for topographic surface description
CN105405166B (zh) 一种基于线性四叉树的lod模型生成方法
CN109118588B (zh) 一种基于块分解的彩色lod模型自动生成方法
CN112181991B (zh) 基于快速构建kd树的地球模拟系统网格重映射方法
KR20100136604A (ko) 3차원 지형 영상 실시간 가시화 시스템 및 이의 방법
CN112115534A (zh) 一种将三维房屋模型转换带高度属性的二维矢量面的方法
Pascal et al. Fast adaptive quadtree mesh generation
CN115661374B (zh) 一种基于空间划分和模型体素化的快速检索方法
EP3563353A1 (en) Systems and methods for lightweight precise 3d visual format
Bartholdi III et al. Vertex‐labeling algorithms for the Hilbert spacefilling curve
CN106204719A (zh) 基于二维邻域检索的三维场景中海量模型实时调度方法
Azri et al. Review of spatial indexing techniques for large urban data management
De Floriani et al. Multiresolution modelling in geographical information systems
CN112102467A (zh) 一种基于gpu的并行八叉树生成、装置及电子设备
Zhang et al. A geometry and texture coupled flexible generalization of urban building models
CN113342999B (zh) 一种基于多层跳序树结构的变分辨率点云简化方法
KR20010046823A (ko) 정렬되지 않은 3차원 거리 데이터로부터 캐드모델 생성 방법
CN113849495A (zh) 一种点云动态哈希划分方法及设备
Prasher et al. Efficient Update and Retrieval of objects in a multiresolution geospatial database
Sivan Surface modeling using quadtrees
CN116310086A (zh) 一种基于四叉树的lod虚拟复杂场景简化方法
Damiand et al. Hierarchical representation for rasterized planar face complexes

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant