CN105913475B - 一种时空变化过程动态可视化方法 - Google Patents
一种时空变化过程动态可视化方法 Download PDFInfo
- Publication number
- CN105913475B CN105913475B CN201610409366.3A CN201610409366A CN105913475B CN 105913475 B CN105913475 B CN 105913475B CN 201610409366 A CN201610409366 A CN 201610409366A CN 105913475 B CN105913475 B CN 105913475B
- Authority
- CN
- China
- Prior art keywords
- time
- space
- change
- intersected
- plane
- 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
Links
- 230000008859 change Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000004540 process dynamic Methods 0.000 title claims abstract description 12
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000009877 rendering Methods 0.000 claims abstract description 18
- 230000003247 decreasing effect Effects 0.000 claims abstract description 6
- 230000009466 transformation Effects 0.000 claims description 17
- 230000002123 temporal effect Effects 0.000 claims description 8
- 238000013519 translation Methods 0.000 claims description 8
- 230000008033 biological extinction Effects 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000005192 partition Methods 0.000 claims description 3
- 230000008520 organization Effects 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 230000000007 visual effect Effects 0.000 abstract description 13
- 238000012800 visualization Methods 0.000 abstract description 7
- 241000406668 Loxodonta cyclotis Species 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 150000001875 compounds Chemical class 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- PLXMOAALOJOTIY-FPTXNFDTSA-N Aesculin Natural products OC[C@@H]1[C@@H](O)[C@H](O)[C@@H](O)[C@H](O)[C@H]1Oc2cc3C=CC(=O)Oc3cc2O PLXMOAALOJOTIY-FPTXNFDTSA-N 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013523 data management Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 235000013350 formula milk Nutrition 0.000 description 1
- 238000004091 panning Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/10—Geometric effects
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N9/00—Details of colour television systems
- H04N9/12—Picture reproducers
- H04N9/31—Projection devices for colour picture display, e.g. using electronic spatial light modulators [ESLM]
- H04N9/3102—Projection devices for colour picture display, e.g. using electronic spatial light modulators [ESLM] using two-dimensional electronic spatial light modulators
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开了一种时空变化过程动态可视化方法，处于时空中的对象都有很强的时空依赖性和动态性，对时空变化过程进行可视化有助于人们更好地了解时空变化的机理。通过以x、y、t为坐标轴，构建时空点、线、面对象的时空立方体模型，在x、y、t构成的三维空间中，使用时间平面求取相交的时空立方体和时空立方体切面，通过颜色控制、显隐控制和时间平面的逐级递增递减，直观而有效地展现时空对象变化的整个过程，整个过程引入可视性计算和GPU顶点缓存，以提升时空对象的浏览体验和渲染效率。本方法直观而动态地描述时空变化的过程，帮助人们更好地理解时空变化的机理，并保证较好的可视化效率和体验，对于时空现象的研究和应用有着重要的意义。
Description
技术领域
本发明涉及GIS(Geographic Information System，地理信息系统)的时空数据可视化领域，具体涉及一种时空变化过程的动态可视化方法。
背景技术
地理时空对象总是处于连续变化之中，具有很强的时空依赖性和动态性特征。近年来，由于人们对时空资源的需求日益增长，时空变化愈加频繁和复杂。时空数据作为时空变化现象的数字化体现，其特征可以概括为两点：1)具有时空的动态性，2)数据日积月累。这些特征给时空变化的可视化带来一定的挑战，时空动态性表现力弱，绘制效率不高，用户体验差，都会给时空变化过程的可视化带来一定困扰。
时空立方体(space-time cube)模型，用几何立体图形描述二维空间沿着时间维的演变过程，在该模型中，空间实体演变历史对应其中一个时空体。可以根据需要在任何时间点上截取一个界面来确定空间实体在改时间点上的状态特征。时空立方体模型对地理变化描述更为形象简洁，适合作为时空变化动态可视化的数据模型。
四叉树索引结构简单清晰而且算法容易实现，是目前广泛应用的空间索引。通过四叉树空间索引，将数据所在的地理空间递归等分为四个子空间，存储于对应的的树结构中，直到树的深度到达某一值或满足某一要求后便停止分割。该四叉树中的根节点和中间节点不存储地理实体对象，叶子节点存储地理对象或对象的集合。在进行时空变化过程浏览时，尤其是空间数据规模较大时，可以根据坐标系姿态解算可见的空间范围，而根据需求调度相应的四叉树节点，能够减少不必要的对象的渲染，提高浏览效率。
OpenGL(Open Graphics Library)顶点缓存(VBO，Vertex Buffer object)直接在高性能显卡(GPU，Graphics Processing Unit)中开辟一个缓存区域来存储顶点数据，根据用户输入的“target”模式和“usage”模式，顶点缓冲区内存管理器将缓冲区对象放在储存器中最佳的位置。因为缓冲区只会在初始化和顶点数据发生变化时才会进行写操作，所以使用顶点缓存大大地减少了CPU-GPU之间的数据拷贝开销，能够显著地提升程序运行的效率。
发明内容
本发明的目的是克服现有技术存在的问题，提供一种时空变化动态可视化方法。以x、y、t为坐标轴，构建时空点、线、面对象的立方体模型，使用时间平面求取相交的立方体和立方体切面，通过颜色控制、显隐控制和时间平面的逐级递增递减，直观而有效地展现时空对象变化的整个过程，在时空对象的研究和应用中具有重要的意义。
本发明的目的是通过以下技术方案来实现的：一种时空变化过程动态可视化方法，包括如下步骤：
步骤1：记录时空对象的出生时间、消亡时间，建立“线性表+四叉树”的时空索引，以文件形式统一存储所要描述的时空对象信息；
步骤2：根据系统状态参数快速定位时空对象数据文件，将数据读入内存，进行平移和缩放操作以适应屏幕坐标系；
步骤3：构建时间平面，与内存中时空对象的时空立方体模型进行一一求交，记录当前相交对象和历史相交对象，将当前相交对象几何形状投影到时间平面；
步骤4：针对时间平面、历史相交对象、当前相交对象、当前相交对象投影结果构建顶点数组和顶点颜色数组，调用OpenGL顶点缓存进行场景渲染；
步骤5：上升或下降一级时间平面，将上一时间间隔的相交对象拷贝到历史相交对象，重复步骤3)和步骤4)；
步骤6：通过设置时间平面间隔参数，较为连续地对整个时空中的时空对象变化过程进行可视化展示。
进一步地，所述的步骤1中所述时空对象包括点、线、面三种形态的对象。
进一步地，所述步骤1具体如下：
(1.1)使用GIS开源库shapeLib从shapefile文件中读取时空对象的几何信息和属性信息；所述属性信息包括出生时间和消亡时间；
(1.2)根据出生时间以一定时间间隔对时空对象进行线性划分，对时空对象在出生时刻的空间几何图形创建四叉树索引。
(1.3)将时空对象以文件形式进行统一存储，存储结构如下：每一个时空对象的结构以一对能表明该时空对象类型的标记性语言作为开始和结束。记录时空对象的标识码、出生时间、消亡时间，以及时空对象在出生时刻的几何图形的边界约束，边界约束为边界点集合，边界点的标识码、坐标X和坐标Y同样记录在该时空对象的数据内容中。
进一步地，所述步骤2具体如下：
(2.1)根据系统状态参数计算整个时空中的可视空间范围；
所述系统状态参数包括当前用于求交的时间平面(iterTime)、世界坐标系(worldframe)的姿态、照相机坐标系(camera frame)的姿态和模型坐标系(model frame)的姿态；
所述可视空间范围的计算如下：相机坐标系与世界坐标系一致，将x,y,z方向的取值均为[-1,1]的世界坐标系根据当前的平移和缩放量进行平移和缩放的坐标逆变换，变换到模型坐标系；将模型坐标系根据初始设定的平移和缩放量进行平移和缩放的坐标逆变换，得到地理空间坐标系下的可视空间范围；
(2.2)每次进行模型平移或缩放变换操作时计算可视空间范围，变化前后的可视空间范围相差大于等于四叉树叶节点最小空间范围一半时重新读取数据；
(2.3)采用线性表的时间索引方式根据所分隔的年份命名根节点文件夹名称；
(2.4)采用四叉树的空间索引方式对空间进行划分，从划分空间的左上角象限开始，顺时针命名并创建文件夹节点，若该节点为叶子节点则创建data.txt文件，该文件统一存储该划分空间中的对象信息。
(2.5)在读取数据时，根据时间平面的值，定位根节点文件夹；再根据可视空间范围检索与该范围相交的四叉树中的父节点集合，对集合中的每个父节点进行递归得到父节点的四叉树检索路径列表。结合根节点文件夹和四叉树检索路径列表快速定位需要读入内存的时空对象数据文件。
进一步地，所述的步骤3包括：
(3.1)将时间平面与时空立方体投影到时间轴(t轴)上，分别得到点的投影和线段的投影，将点的投影和线段的投影进行相交计算，得到可视空间范围的当前相交对象；根据不断更新的当前相交对象得到历史相交对象，历史相交对象和当前相交对象不重复，且相对于当前时间平面处于消亡状态；
(3.2)计算当前相交对象在时间平面上的投影：首先取时空对象的时空立方体模型的任意一条时间线段与时间平面求交，计算得出相交点与时空对象边界约束点的时间距离td，然后将时空对象所有边界约束点在t轴上平移td，得到投影结果的边界约束点集合。
进一步地，所述的步骤4中采用glBufferData()为顶点缓存对象申请内存空间并进行初始化，glBufferData()的usage参数选取GL_DYNAMIC_DRAW。对时间平面、历史相交对象、当前相交对象、当前相交对象投影结果使用不同的颜色渲染方式，其中当前相交对象投影的渲染需要关闭深度缓存测试glDisable(GL_DEPTH_TEST)。
进一步地，所述的步骤6包括：
(6.1)根据时空数据在生存周期上的分布特征选取时间平面间隔参数；
(6.2)在x,y,t构成的三维空间中，通过逐级递增或递减时间平面，观察时空变化在不同阶段的特征、某一时空对象的变化历程，以及整个时空的演变过程；
(6.3)通过对点、线、面时空对象，以及当前时空对象、历史时空对象、当前对象投影结果进行显隐控制，以更好地观察整个时空变化的过程。
本发明与现有技术相比具有有益效果：
(1)本发明为GIS时空数据动态可视化提供了一种崭新的思考方式和实施方法。采用时间切面和时空立方体的表现形式，相比于现有的可视化方法，更加直观、简洁、形象。
(2)本发明采用线性表+四叉树的时空索引，结构清晰、算法容易实现，而且能够满足时空变化动态可视化的需求。所提出的结合四叉树的空间范围可视性计算，能够减少不必要对象的可视化并且较好地进行空间数据的调度，从而提升可视化的用户体验。
(3)本发明使用OpenGL顶点缓存技术进行场景渲染，大大提升了场景渲染的效率，增强了整个时空变化过程的可操作性。
附图说明
图1为本发明基于可视性计算的数据加载效果，(a)为系统初始状态数据加载情况，(b)为系统缩放操作后数据加载情况；
图2为本发明时空变化不同阶段的截图，(a)为时间平面处于1991年4月1日的地籍时空对象状态，(b)为时间平面处于1994年2月13日的地籍时空对象状态，(c)为时间平面处于1998年11月21日的地籍时空对象状态，(d)为时间平面处于2002年8月18日的地籍时空对象状态；
图3为本发明时空变化过程中不同时空对象的显隐控制，(a)为全部显示的时空对象情况，(b)为不显示当前相交对象投影结果的时空对象情况，(c)为不显示当前相交对象投影结果和历史对象的时空对象情况，(d)为只显示点对象的时空对象情况，(e)为只显示线对象的时空对象情况，(f)为只显示面对象时空对象情况，(g)为只显示当前相交对象投影结果的时空对象情况，(h)为当前相交对象投影结果的三维浏览。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明。
本发明提供的一种时空变化过程的动态可视化方法，包括如下步骤：
步骤1：记录时空对象的出生时间、消亡时间，建立“线性表+四叉树”的时空索引，以文件形式统一存储所要描述的时空对象信息；
步骤2：根据系统状态参数快速定位时空对象数据文件，将数据读入内存，进行平移和缩放操作以适应屏幕坐标系；
步骤3：构建时间平面，与内存中时空对象的时空立方体模型进行一一求交，记录当前相交对象m_interPrimitives和历史相交对象m_lsinterPrimitives，将当前相交对象m_interPrimitives几何形状投影到时间平面；
步骤4：针对时间平面、历史相交对象、当前相交对象、当前相交对象投影结果构建顶点数组和顶点颜色数组，调用OpenGL顶点缓存进行场景渲染；
步骤5：上升或下降一级时间平面，将上一时间间隔的相交对象m_interPrimitives拷贝到历史相交对象m_lsinterPrimitives，重复步骤3)和步骤4)；
步骤6：通过设置时间平面间隔参数，较为连续地对整个时空中的时空对象变化过程进行可视化展示。
所述的步骤1包括：
(1)所述时空对象包括点、线、面三种形态的对象。
(2)基于传统shapefile的空间数据管理模式的应用开发，受限于ArcGIS版权的约束，无法进行很好的扩展，所以采用文件形式进行时空对象的统一存储。
(3)使用GIS开源库shapeLib读取时空对象的几何信息和属性信息，分别是函数SHPOpen()和DBFOpen()；所述属性信息包括出生时间和消亡时间。其中padfX和padfY为指向对象几何图形端点的指针，通过padfX++和padfY++可以遍历所有节点获取相应坐标信息。而使用函数DBFGetFieldIndex()可以获取某一属性字段的索引号，进而使用函数DBFReadStringAttribute()读取属性字段具体值。
(4)近年来人们对于时空资源的需求日益增长，时空变化频繁，时空数据空间范围广且在时间维度上加速累积。为了提高读取和渲染的速度，根据出生时间以一定时间间隔对数据进行线性划分，在此基础上，对时空对象在出生时刻的空间几何图形创建四叉树索引。
(5)根据不同对象类型进行分层存储的传统形式，在时空变化动态可视化中没有优势，反而不利于统一的计算与渲染，和数据文件的快速检索。所述用于动态可视化的时空对象统一存储结构形式化表达如下：
＜ObjectType＞
[ObjectID,BeginTime,EndTime,Attr1,Attr2,…]＜ObjectConstraints＞
[Point,PointID,X,Y]
＜/ObjectType＞
所述的步骤2包括：
(1)整个过程中世界坐标系和相机坐标系重合。
(2)根据系统状态参数计算整个时空中的可视空间范围；
所述系统状态参数包括当前用于求交的时间平面(iterTime)，世界坐标系(worldframe)、照相机坐标系(camera frame)和模型坐标系(model frame)的姿态。
所述可视空间范围的计算如下：空间范围的计算只涉及x和y坐标信息，将三维坐标变换简化为二维坐标变换；相机坐标系与世界坐标系一致，将x,y,z方向的取值均为[-1,1]的世界坐标系根据当前的平移和缩放量进行平移和缩放的坐标逆变换，变换到模型坐标系；将模型坐标系根据初始设定的平移和缩放量进行平移和缩放的坐标逆变换，得到可视空间范围。变换公式如下：
二维坐标变换：
其中，S为缩放矩阵、T为平移矩阵、X为变换前的坐标矩阵、X′为变换后的坐标矩阵。
二维坐标逆变换：
其中，S-1为缩放逆矩阵、T-1为平移逆矩阵X为变换前的坐标矩阵、X′为变换后的坐标矩阵。
(3)如果每进行一次模型变换(平移、旋转、缩放)都重新进行数据调度，会对场景渲染、顶点缓存区更新造成一定压力，因此只考虑平移和缩放变换且当平移和缩放的累积量对应到空间范围的变化到达一定值时才重新进行数据读写，该值根据四叉树子结点划分的最小空间范围大小而定。若变化前后空间范围相差大于四叉树叶节点最小空间范围一半时，即可重新读取数据。变换累积量计算依赖于复合平移变换和复合缩放变换计算：
复合平移变换：S(s2x,s2y)·S(s1x,s1y)＝S(s2x·s1x,s2y·s1y)
复合缩放变换：T(t2x,t2y)·T(t1x,t1y)＝T(t2x+t1x,t2y+t1y)
(4)采用线性表的时间索引方式根据所分隔的年份命名根节点文件夹名称，例如1990，1995，2000。
(5)采用四叉树的空间索引方式对空间进行划分，从划分空间的左上角象限开始，顺时针命名并创建文件夹节点：0，1，2，3。若该节点为叶子节点则创建data.txt文件，该文件统一存储该划分空间中的对象信息。
(6)在读取数据时，根据时间平面的值，定位根节点文件夹；再根据可视空间范围检索与该范围相交的四叉树中的父节点集合，对集合中的每个父节点进行递归得到父节点的四叉树检索路径列表。对其中一项进行举例，例如“/1/2/1/3/”，使用strcat()函数在该字串前增加根节点年份字串，即可以快速定位需要读入内存的数据文件。
(7)父节点递归检索的函数如下：
所述的步骤3包括：
(1)将时间平面与时空立方体投影到时间轴(t轴)上，分别得到点的投影和线段的投影，将点的投影和线段的投影进行相交计算，即点的t坐标值与线段最小t坐标值和最大t坐标值的比较，得到可视空间范围的当前相交对象。
(2)根据不断更新的当前相交对象得到历史相交对象m_lsinterPrimitives，历史相交对象m_lsinterPrimitives和当前相交对象m_interPrimitives不重复，且相对于当前时间平面处于消亡状态；
(3)计算当前相交对象在时间平面上的投影：首先取时空对象的时空立方体模型的任意一条时间线段与时间平面求交，计算得出相交点与时空对象边界约束点的时间距离td，然后将时空对象所有边界约束点在t轴上平移td，得到投影结果的边界约束点集合。
所述的步骤4包括：
(1)传统glVertex()的方式是每次往GPU传顶点，CPU与GPU的数据传输消耗大，而顶点缓存是直接在GPU中开辟一个缓存区，一次性将一个场景所需的顶点数据传入GPU，减少了CPU与GPU的传输消耗，能够较大提升场景渲染的效率。
(2)具体实现主要包含以下顶点缓存函数：
glGenBuffers()创建顶点缓存对象；glBindBuffer()将顶点缓存对象设置为当前数组缓存对象(参数target取值GL_ARRAY_BUFFER)；glBufferData()为顶点缓存对象申请内存空间并进行初始化；glBufferSubData()初始化或更新顶点缓存对象；glDeleteBuffers()销毁缓存对象。
(3)随着时间平面的逐级递增或递减，缓存区会被周期性更改，因此glBufferData()的usage参数选取GL_DYNAMIC_DRAW。
(4)对时间平面、历史相交对象、当前相交对象、当前相交对象投影结果使用不同的颜色渲染方式：
时间平面顶点：透明灰色面，glColor4fm(0.7,0.7,0.7,0.5)；
历史相交对象顶点：灰色面，glColor3f(0.9,0.9,0.9)；
当前相交对象顶点：多色面渲染，针对类型分色，例如土地用途；
当前相交对象投影结果顶点：红色线框，glColorf3f(1.0,0.0,0.0)，需要关闭深度缓存测试glDisable(GL_DEPTH_TEST)；
(5)点、线、面时空立方体模型的绘制模式：
点：GL_POINTS；线：面(GL_POLYGON)+首尾相接线(GL_LINE_LOOP)双面绘制，通过控制法线方向实现；面：面(GL_POLYGON)+首尾相接线(GL_LINE_LOOP)，剔除背面。
所述的步骤6包括：
(1)根据时空数据在生存周期上的分布特征选取时间平面间隔参数。
(2)在x,y,t构成的三维空间中，通过逐级递增或递减时间平面，观察时空变化在不同阶段的特征、某一时空对象的变化历程，以及整个时空的演变过程。
(3)通过对点、线、面时空对象，以及现实时空对象、历史时空对象、现实对象投影结果进行显隐控制，以更好地观察整个时空变化的过程。
实施例：
以地籍时空数据的时空变化过程动态可视化为例。
第一步是使用shapeLib逐个读取现有地籍要素的出生时间、消亡时间等信息，读取过程中，根据要素的出生时间和空间包围盒坐标范围，建立“线性表+四叉树”的时空索引，并存储到相应路径下的data.txt中。
第二步是根据当前时间平面(默认19900000)在时间线性表中检索，定位父节点文件夹“1990”；根据模型坐标系姿态T(2,4)和S(1/10,1/10)，通过坐标逆变换计算模型可视空间范围(-30,-50)～(-10,-30)，根据初始设定的地理坐标往模型坐标转换时的固定初始缩放比例(1/100)和平移值(-430236.0,-3261083.0)，进行坐标先缩放和平移的逆变换，计算得到所要显示的地理空间范围(427236.0,3256083.0)～(429236.0,3258083.0)，在空间四叉树中进行递归检索，得到检索路径“/1/3/3”、“/1/3/4”、“/4/2/1”和“/4/2/2”，将四个路径下的data.txt数据读入内存，如图1所示；
第三步是构建时间平面，通过线段求交方式与内存中时空对象的时空立方体模型进行一一求交，记录当前相交对象m_interPrimitives和历史相交对象m_lsinterPrimitives，计算当前相交对象在时间平面上的几何形状投影；
第四步是针对时间平面、历史相交对象、当前相交对象、当前相交对象投影结果构建顶点数组和顶点颜色数组，调用OpenGL顶点缓存进行场景渲染，如图2所示；
第五步是上升或下降一级时间平面(按照间隔5000进行加或减)，将上一时间间隔的相交对象拷贝到历史相交对象m_lsinterPrimitives，重复步骤3)和步骤4)；
第六步是在连续调整时间平面的过程中(如图2所示)，在视窗中操作时空对象的显隐性，可选择显示点、线、面的时空对象中的一种或多种，可选择只显示或不显示当前相交对象在时间平面的投影结果，可选择只显示或不显示当前时间状态的历史时空对象。如图3所示。
Claims (7)
1.一种时空变化过程动态可视化方法，其特征在于，包括如下步骤：
步骤1：记录时空对象的出生时间、消亡时间，建立“线性表+四叉树”的时空索引，以文件形式统一存储所要描述的时空对象信息；
步骤2：根据系统状态参数快速定位时空对象数据文件，将数据读入内存，进行平移和缩放操作以适应屏幕坐标系；
步骤3：构建时间平面，与内存中时空对象的时空立方体模型进行一一求交，记录当前相交对象和历史相交对象，将当前相交对象几何形状投影到时间平面；
步骤4：针对时间平面、历史相交对象、当前相交对象、当前相交对象投影结果构建顶点数组和顶点颜色数组，调用OpenGL顶点缓存进行场景渲染；
步骤5：上升或下降一级时间平面，将上一时间间隔的相交对象拷贝到历史相交对象，重复步骤3和步骤4；
步骤6：通过设置时间平面间隔参数，较为连续地对整个时空中的时空对象变化过程进行可视化展示。
2.根据权利要求1所述的一种时空变化过程动态可视化方法，其特征在于，所述的步骤1中所述时空对象包括点、线、面三种形态的对象。
3.根据权利要求1所述的一种时空变化过程动态可视化方法，其特征在于，所述步骤1具体如下：
(1.1)使用GIS开源库shapeLib从shapefile文件中读取时空对象的几何信息和属性信息；所述属性信息包括出生时间和消亡时间；
(1.2)根据出生时间以一定时间间隔对时空对象进行线性划分，对时空对象在出生时刻的空间几何图形创建四叉树索引；
(1.3)将时空对象以文件形式进行统一存储，存储结构如下：每一个时空对象的结构以一对能表明该时空对象类型的标记性语言作为开始和结束；记录时空对象的标识码、出生时间、消亡时间，以及时空对象在出生时刻的几何图形的边界约束，边界约束为边界点集合，边界点的标识码、坐标X和坐标Y同样记录在该时空对象的数据内容中。
4.根据权利要求1所述的一种时空变化过程动态可视化方法，其特征在于，所述步骤2具体如下：
(2.1)根据系统状态参数计算整个时空中的可视空间范围；
所述系统状态参数包括当前用于求交的时间平面、世界坐标系的姿态、照相机坐标系的姿态和模型坐标系的姿态；
所述可视空间范围的计算如下：相机坐标系与世界坐标系一致，将x,y,z方向的取值均为[-1,1]的世界坐标系根据当前的平移和缩放量进行平移和缩放的坐标逆变换，变换到模型坐标系；将模型坐标系根据初始设定的平移和缩放量进行平移和缩放的坐标逆变换，得到地理空间坐标系下的可视空间范围；
(2.2)每次进行模型平移或缩放变换操作时计算可视空间范围，变化前后的可视空间范围相差大于等于四叉树叶节点最小空间范围一半时重新读取数据；
(2.3)采用线性表的时间索引方式根据所分隔的年份命名根节点文件夹名称；
(2.4)采用四叉树的空间索引方式对空间进行划分，从划分空间的左上角象限开始，顺时针命名并创建文件夹节点，若该节点为叶子节点则创建data.txt文件，该文件统一存储该划分空间中的对象信息；
(2.5)在读取数据时，根据时间平面的值，定位根节点文件夹；再根据可视空间范围检索与该范围相交的四叉树中的父节点集合，对集合中的每个父节点进行递归得到父节点的四叉树检索路径列表；结合根节点文件夹和四叉树检索路径列表快速定位需要读入内存的时空对象数据文件。
5.根据权利要求1所述的一种时空变化过程动态可视化方法，其特征在于，所述的步骤3包括：
(3.1)将时间平面与时空立方体投影到时间轴t轴上，分别得到点的投影和线段的投影，将点的投影和线段的投影进行相交计算，得到可视空间范围的当前相交对象；根据不断更新的当前相交对象得到历史相交对象，历史相交对象和当前相交对象不重复，且相对于当前时间平面处于消亡状态；
(3.2)计算当前相交对象在时间平面上的投影：首先取时空对象的时空立方体模型的任意一条时间线段与时间平面求交，计算得出相交点与时空对象边界约束点的时间距离td，然后将时空对象所有边界约束点在t轴上平移td，得到投影结果的边界约束点集合。
6.根据权利要求1所述的一种时空变化过程动态可视化方法，其特征在于，所述的步骤4中采用glBufferData()为顶点缓存对象申请内存空间并进行初始化，glBufferData()的usage参数选取GL_DYNAMIC_DRAW；对时间平面、历史相交对象、当前相交对象、当前相交对象投影结果使用不同的颜色渲染方式，其中当前相交对象投影的渲染需要关闭深度缓存测试glDisable(GL_DEPTH_TEST)。
7.根据权利要求1所述的一种时空变化过程动态可视化方法，其特征在于，所述的步骤6包括：
(6.1)根据时空数据在生存周期上的分布特征选取时间平面间隔参数；
(6.2)在x,y,t构成的三维空间中，通过逐级递增或递减时间平面，观察时空变化在不同阶段的特征、某一时空对象的变化历程，以及整个时空的演变过程；
(6.3)通过对点、线、面时空对象，以及当前时空对象、历史时空对象、当前对象投影结果进行显隐控制，以更好地观察整个时空变化的过程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610409366.3A CN105913475B (zh) | 2016-06-12 | 2016-06-12 | 一种时空变化过程动态可视化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610409366.3A CN105913475B (zh) | 2016-06-12 | 2016-06-12 | 一种时空变化过程动态可视化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105913475A CN105913475A (zh) | 2016-08-31 |
| CN105913475B true CN105913475B (zh) | 2019-01-25 |
Family
ID=56749931
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610409366.3A Active CN105913475B (zh) | 2016-06-12 | 2016-06-12 | 一种时空变化过程动态可视化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105913475B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109255004A (zh) * | 2018-09-18 | 2019-01-22 | 中国测绘科学研究院 | 面向地理空间数据的精度保证的在线交互式可视化方法 |
| CN109284291B (zh) * | 2018-09-28 | 2022-03-22 | 上海杰狮信息技术有限公司 | 一种基于自适应多基态修正模型的高效时空数据检索方法 |
| CN111783186B (zh) * | 2020-07-03 | 2021-03-16 | 北京安怀信科技股份有限公司 | 通用cad模型结果轻量可视化查看方法 |
| CN112634430B (zh) * | 2020-11-17 | 2022-06-21 | 青岛中海潮科技有限公司 | 体可视化方法、装置、设备和存储介质 |
| CN114742887B (zh) * | 2022-03-02 | 2023-04-18 | 广东工业大学 | 一种基于点线面特征融合的无人机位姿估计方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101231642A (zh) * | 2007-08-27 | 2008-07-30 | 中国测绘科学研究院 | 时空数据库管理方法及其系统 |
| CN102023983A (zh) * | 2009-09-11 | 2011-04-20 | 首都师范大学 | 统计时空数据库及其管理方法 |
| CN103714167A (zh) * | 2013-12-31 | 2014-04-09 | 上海师范大学 | 一种灾害信息时空数据在线三维显示方法 |
| CN103914558A (zh) * | 2014-04-16 | 2014-07-09 | 中南大学 | 一种基于时空统计的气象要素时空聚集模式挖掘方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090064051A1 (en) * | 2007-08-31 | 2009-03-05 | Florida Atlantic University | Interactive system for visualization and retrieval of video data |
-
2016
- 2016-06-12 CN CN201610409366.3A patent/CN105913475B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101231642A (zh) * | 2007-08-27 | 2008-07-30 | 中国测绘科学研究院 | 时空数据库管理方法及其系统 |
| CN102023983A (zh) * | 2009-09-11 | 2011-04-20 | 首都师范大学 | 统计时空数据库及其管理方法 |
| CN103714167A (zh) * | 2013-12-31 | 2014-04-09 | 上海师范大学 | 一种灾害信息时空数据在线三维显示方法 |
| CN103914558A (zh) * | 2014-04-16 | 2014-07-09 | 中南大学 | 一种基于时空统计的气象要素时空聚集模式挖掘方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105913475A (zh) | 2016-08-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105913475B (zh) | 一种时空变化过程动态可视化方法 | |
| Brath | 3D InfoVis is here to stay: Deal with it | |
| US8154544B1 (en) | User specified contact deformations for computer graphics | |
| KR100888528B1 (ko) | 3Ｄ 오브젝트의 셀프 셰도우(ｓｅｌｆ-ｓｈａｄｏｗ)를만들기 위해 데이터를 미리-저장할 수 있는 장치, 방법,응용프로그램 및 컴퓨터 판독가능한 매체 | |
| Li et al. | Visualizing dynamic geosciences phenomena using an octree-based view-dependent LOD strategy within virtual globes | |
| Hurter | Image-based visualization: Interactive multidimensional data exploration | |
| CN110503718B (zh) | 三维工程模型轻量化显示方法 | |
| CN108629834A (zh) | 一种基于单张图片的三维头发重建方法 | |
| CN109063279A (zh) | 基于粒子流轨迹追踪算法的三维空间风场动态模拟方法 | |
| Hastings et al. | Optimization of large-scale, real-time simulations by spatial hashing | |
| Fan et al. | An octree-based proxy for collision detection in large-scale particle systems | |
| CN110400366B (zh) | 一种基于OpenGL的实时洪水灾害可视化仿真方法 | |
| CN115859755A (zh) | 一种定常流场矢量数据的可视化方法、装置、设备及介质 | |
| Leist et al. | Interactive visualisation of spins and clusters in regular and small-world Ising models with CUDA on GPUs | |
| Ready et al. | Immersive visualisation of big data for river disaster management | |
| Boorboor et al. | Submerse: Visualizing storm surge flooding simulations in immersive display ecologies | |
| CN108875275A (zh) | 一种基于大规模流线的矢量场实时仿真方法及系统 | |
| She et al. | A building label placement method for 3D visualizations based on candidate label evaluation and selection | |
| US20030193501A1 (en) | Incremental resolution changes in multi-resolution meshes with update records | |
| Xuexian et al. | High-performance navigation and rendering of very-large scale landscape and seascape | |
| Xing et al. | MR environments constructed for a large indoor physical space | |
| Yan et al. | 3D visual design for mobile search result on 3G mobile phone | |
| Wang et al. | 3D Model Management and Dynamic Loading based on GIS Database | |
| JP2011248675A (ja) | 隠蔽処理プログラム、可視化処理方法及び装置 | |
| Wang et al. | A 3D Visualization Method of Global Ocean Surface Based on Discrete Global Grids |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |