CN102364522B - 一种基于全球划分的大规模气象体数据绘制方法 - Google Patents

一种基于全球划分的大规模气象体数据绘制方法 Download PDF

Info

Publication number
CN102364522B
CN102364522B CN 201110325178 CN201110325178A CN102364522B CN 102364522 B CN102364522 B CN 102364522B CN 201110325178 CN201110325178 CN 201110325178 CN 201110325178 A CN201110325178 A CN 201110325178A CN 102364522 B CN102364522 B CN 102364522B
Authority
CN
China
Prior art keywords
volume data
meteorological
subdivision
piecemeal
initial
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
CN 201110325178
Other languages
English (en)
Other versions
CN102364522A (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.)
Wuhan Zhao Tu Science and Technology Ltd.
Original Assignee
China University of Geosciences
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 China University of Geosciences filed Critical China University of Geosciences
Priority to CN 201110325178 priority Critical patent/CN102364522B/zh
Publication of CN102364522A publication Critical patent/CN102364522A/zh
Application granted granted Critical
Publication of CN102364522B publication Critical patent/CN102364522B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Generation (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明公开了一种基于全球划分的大规模气象体数据绘制方法，其包括如下步骤：对原始气象体数据进行预处理并对地球上空高度为h1至h2的空间进行划分；根据每一最终剖分子块的空间坐标范围获得体数据块，并对所有体数据块进行颜色映射，得到体数据块颜色索引块；依据各个体数据块与视点P1之间的距离对体数据块进行排序，得到体数据块的绘制列表；根据所得绘制列表，利用GPU调色板与所得体数据块颜色索引块对绘制列表中的每一体数据块进行体绘制；随机移动视点，重复上述步骤。本发明的主要优点在于基于全球划分的策略对规则气象体数据进行绘制，从而突破了GPU存储容量的瓶颈，达到在球面上高效、快速绘制大规模气象体数据的目的。

Description

一种基于全球划分的大规模气象体数据绘制方法
技术领域
本发明涉及一种基于全球划分的大规模气象体数据绘制方法。
背景技术
随着科技的发展，人们需要处理的数据也越来越复杂，这就需要计算机辅助可视化将复杂的数据信息转换为直观的图形图像并在屏幕上显示出来。海量数据只有通过可视化变成形象，人们才能从表面杂乱无章的海量数据中挖掘出其隐藏的规律，才能为科学发现、气象学、业务决策等领域提供依据。体绘制是一种有效的数据可视化手段，能够揭示体数据内部信息，而不仅仅是显示对象的表面信息，在气象领域取得了广泛的应用。
然而，目前体数据绘制技术远未成熟，在理论研究和实际应用中还有许多问题亟待解决。其存在的问题具体表现在以下几个方面：对于大规模的气象体数据而言，基于球面进行体绘制还没有成熟的方法；当需要对不同时区下、不同投影下的多个体数据绘制时，目前的体绘制方法没办法集成绘制，绘制效率不高；虽然随着计算机图形硬件的发展，基于GPU的体绘制算法与传统的基于CPU的体绘制算法相比，绘制速度有了极大的提高，但是当数据规模较大时，体绘制算法仍然难以达到实时绘制的速度；由于受GPU存储容量的限制，在绘制大规模的体数据时，例如大规模的气象体数据，绘制能力大大受限。
发明内容
为解决现有技术中的上述问题，本发明提供一种基于全球划分的大规模气象体数据绘制方法，本发明的主要优点在于基于全球划分的策略对规则气象体数据进行绘制，从而突破了GPU存储容量的瓶颈，达到在球面上高效、快速绘制大规模气象体数据的目的。
为达成上述目的，本发明采用如下技术方案：一种基于全球划分的大规模气象体数据绘制方法，其包括如下步骤：
（1）对原始气象体数据进行预处理并对地球上空高度为h1至h2的空间进行划分；
其中对原始气象体数据进行预处理是对原始气象体数据进行插值预处理，使原始非规则气象体数据变成规则气象体数据；
其中对地球上空高度为h1至h2的空间进行划分的步骤如下：
（I）对地球上空高度为h1至h2的空间进行全球初始划分，将此空间在经度上均匀划分成m个行、在纬度上均匀划分成n个列、在高度上均匀划分成k个层，得到m*n*k个初始分块，并对所得到的所有初始分块进行编号；
（II）确定视点在位置P1处时的可见视域范围，得到此视域范围内的所有q块可见初始分块的列表；
（III）在以位置P1为原点的视口坐标系下，计算步骤（II）中得到的所有可见初始分块中的每一可见初始分块与视点之间的距离di，1≤i≤q，并根据该距离di确定步骤（II）中得到的列表中每一可见初始分块所需的再次剖分次数ni，每一可见初始分块所需的再次剖分次数ni随着每一可见初始分块与视点之间的距离di的减小而增大，然后根据所得的每一可见初始分块的再次剖分次数ni对对应的可见初始分块进行剖分，得到最终剖分子块；
（2）根据每一最终剖分子块的空间坐标范围获取落在其坐标范围内的体数据信息，包含体数据信息的剖分子块形成体数据块，并对所有体数据块进行颜色映射，得到体数据块颜色索引块；
（3）依据各个体数据块与视点P1之间的距离、由近及远地对所有体数据块进行排序，得到体数据块的绘制列表；
（4）根据步骤（3）中得到的绘制列表，利用GPU调色板与步骤（2）中得到的体数据块颜色索引块的信息对步骤（3）中得到的绘制列表中的每一体数据块进行体绘制；
（5）将视点依次随机移动至位置P2、P3、P4、P5、P6……Pn，其中n为自然数，重复上述步骤（II）、（III）、（2）、（3）、（4）。
在上述步骤（1）中，用克里金插值方法对原始气象体数据进行插值预处理将原始非规则气象体数据变成规则气象体数据，其中原始气象体数据是非规则层状数据，其层与层之间是非等间距的，而规则气象体数据是层与层之间是等间距的层状数据。
在上述步骤（I）中，将位于第a行第b列第c层处的初始分块编号为（a，b，c），其中0≤a≤m-1、0≤b≤n-1、0≤c≤k-1。
在步骤（II）中，得到P1的视域范围内的所有q块可见初始分块（a1，b1，c1）、（a2，b2，c2）……的列表。
在步骤（III）中，采用八叉树剖分法根据与每一可见初始分块对应的再次剖分次数ni对该可见初始分块进行剖分。
与现有技术相比，本发明提供的基于全球划分的大规模气象体数据绘制方法具有以下优点：基于全球划分的策略对规则气象体数据进行绘制，剥离了与数据的相关性，可以将不同时区下、不同投影下的气象体数据集成绘制，能处理大规模的气象体数据，并可以在保证图像质量的同时大大地提高了绘制的速度；根据视点计算球面下视域内的可见分块，然后只对可见分块进行绘制，降低了对图形硬件资源的占用，同时提高了绘制效率；若视点移动，可动态地得到视点与初始分块的距离，并实时计算出分块剖分的级数，达到实时绘制的目的；利用GPU调色板记录颜色编号，体数据块的颜色值就根据颜色索引块找相应的颜色，有效地解决了GPU存储容量的瓶颈问题，充分利用了图形硬件的加速功能，提高了大规模数据的绘制速度；本发明实现了对大规模气象体数据在球面上的快速绘制，速度可以达到实时交互；本发明实现的基于球面的大规模气象体数据绘制效果形象逼真。
附图说明
图1是本发明的基于全球划分的大规模气象体数据绘制方法的流程图。
图2是本发明的初始划分的示意图。
具体实施方式
下面将结合具体实施例对本发明提供的基于全球划分的大规模气象体数据绘制方法做进一步详细说明。
步骤1：对原始气象体数据进行预处理：
用克里金插值方法对原始气象体数据进行插值预处理将原始非规则气象体数据变成规则气象体数据。其中原始气象体数据是非规则层状数据，其层与层之间是非等间距的，而规则气象体数据是层与层之间是等间距的层状数据。
步骤2：按照如下步骤对地球上空高度为h1至h2的空间进行划分：
（I）对地球上空高度为h1至h2的空间进行全球初始划分，将此空间在经度上均匀划分成m个行、在纬度上均匀划分成n个列、在高度上均匀划分成k个层，得到m*n*k个初始分块，并对所得到的所有初始分块进行编号，将位于第a行第b列第c层处的初始分块编号为（a，b，c），其中0≤a≤m-1、0≤b≤n-1、0≤c≤k-1。
（II）确定视点在位置P1处时的可见视域范围，得到此视域范围内的所有q块可见初始分块（a1，b1，c1）、（a2，b2，c2）……的列表。
（III）在以位置P1为原点的视口坐标系下，计算步骤（II）中得到的所有可见初始分块中的每一可见初始分块与视点之间的距离di，1≤i≤q，并根据该距离di确定步骤（II）中得到的列表中每一可见初始分块所需的再次剖分次数ni，每一可见初始分块所需的再次剖分次数ni随着每一可见初始分块与视点之间的距离di的减小而增大，初始分块距视点的距离越近所需的再次剖分次数越大，然后采用八叉树剖分法根据所得的每一可见初始分块的再次剖分次数ni对对应的可见初始分块进行剖分，得到最终剖分子块。
步骤3：根据每一最终剖分子块的空间坐标范围获取落在其坐标范围内的体数据信息，包含体数据信息的剖分子块形成体数据块，并对所有体数据块进行颜色映射，得到体数据块颜色索引块。
步骤4：依据各个体数据块与视点P1之间的距离、由近及远地对所有体数据块进行排序，得到体数据块的绘制列表。
步骤5：根据步骤4中得到的绘制列表，利用GPU调色板与步骤3中得到的体数据块颜色索引块的信息对步骤4中得到的绘制列表中的每一体数据块进行体绘制。
步骤6：将视点依次随机移动至位置P2、P3、P4、P5或P6……Pn，其中n为自然数，重复上述步骤（II）、（III）、3、4、5。视点随时间而移动，视点的移动会造成视域范围的变化，所得的可见分块列表也会产生变化，且可见分块与视点的距离也会发生变化，随机移动视点，实时计算这些参数，从而保证整个绘制过程是流畅、动态实时的。
本发明提出的全球划分的策略是基于球面的，全球划分虽然也是分块的思想，但划分的是空间范围，而不是对实际的体数据进行分块，可以将不同时区下、不同投影下的气象体数据集成绘制，能处理大规模的气象体数据，并可以在保证图像质量的同时大大地提高了绘制的速度。本发明解决的是球面上的大规模体数据的绘制方法，实际中，气象数据是在地球上空的，因此气象领域中气象体数据的绘制在实际展示中，更适用采用本发明方法。

Claims (5)

1.一种基于全球划分的大规模气象体数据绘制方法，其包括如下步骤：
（1）对原始气象体数据进行预处理并对地球上空高度为h1至h2的空间进行划分；
其中对原始气象体数据进行预处理是对原始气象体数据进行插值预处理，使原始非规则气象体数据变成规则气象体数据；
其中对地球上空高度为h1至h2的空间进行划分的步骤如下：
（I）对地球上空高度为h1至h2的空间进行全球初始划分，将此空间在经度上均匀划分成m个行、在纬度上均匀划分成n个列、在高度上均匀划分成k个层，得到m*n*k个初始分块，并对所得到的所有初始分块进行编号；
（II）确定视点在位置P1处时的可见视域范围，得到此视域范围内的所有q块可见初始分块的列表；
（III）在以位置P1为原点的视口坐标系下，计算步骤（II）中得到的所有可见初始分块中的每一可见初始分块与视点之间的距离di，1≤i≤q，并根据该距离di确定步骤（II）中得到的列表中每一可见初始分块所需的再次剖分次数ni，每一可见初始分块所需的再次剖分次数ni随着每一可见初始分块与视点之间的距离di的减小而增大，然后根据所得的每一可见初始分块的再次剖分次数ni对对应的可见初始分块进行剖分，得到最终剖分子块；
（2）根据每一最终剖分子块的空间坐标范围获取落在其坐标范围内的体数据信息，包含体数据信息的剖分子块形成体数据块，并对所有体数据块进行颜色映射，得到体数据块颜色索引块；
（3）依据各个体数据块与视点P1之间的距离、由近及远地对所有体数据块进行排序，得到体数据块的绘制列表；
（4）根据步骤（3）中得到的绘制列表，利用GPU调色板与步骤（2）中得到的体数据块颜色索引块的信息对步骤（3）中得到的绘制列表中的每一体数据块进行体绘制；
（5）将视点依次随机移动至位置P2、P3、P4、P5、P6……Pn，其中n为自然数，重复上述步骤（II）、（III）、（2）、（3）、（4）。
2.根据权利要求1所述的基于全球划分的大规模气象体数据绘制方法，其特征在于：在上述步骤（1）中，用克里金插值方法对原始气象体数据进行插值预处理将原始非规则气象体数据变成规则气象体数据，其中原始气象体数据是非规则层状数据，其层与层之间是非等间距的，而规则气象体数据是层与层之间是等间距的层状数据。
3.根据权利要求1所述的基于全球划分的大规模气象体数据绘制方法，其特征在于：在上述步骤（I）中，将位于第a行第b列第c层处的初始分块编号为（a，b，c），其中0≤a≤m-1、0≤b≤n-1、0≤c≤k-1。
4.根据权利要求1所述的基于全球划分的大规模气象体数据绘制方法，其特征在于：在步骤（II）中，得到P1的视域范围内的所有q块可见初始分块（a1，b1，c1）、（a2，b2，c2）……的列表。
5.根据权利要求1所述的基于全球划分的大规模气象体数据绘制方法，其特征在于：在步骤（III）中，采用八叉树剖分法根据与每一可见初始分块对应的再次剖分次数ni对该可见初始分块进行剖分。
CN 201110325178 2011-10-24 2011-10-24 一种基于全球划分的大规模气象体数据绘制方法 Active CN102364522B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110325178 CN102364522B (zh) 2011-10-24 2011-10-24 一种基于全球划分的大规模气象体数据绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110325178 CN102364522B (zh) 2011-10-24 2011-10-24 一种基于全球划分的大规模气象体数据绘制方法

Publications (2)

Publication Number Publication Date
CN102364522A CN102364522A (zh) 2012-02-29
CN102364522B true CN102364522B (zh) 2013-02-13

Family

ID=45691088

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110325178 Active CN102364522B (zh) 2011-10-24 2011-10-24 一种基于全球划分的大规模气象体数据绘制方法

Country Status (1)

Country Link
CN (1) CN102364522B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106845477A (zh) * 2016-12-30 2017-06-13 武汉联影医疗科技有限公司 基于多个重建图像的感兴趣区域建立方法及其装置

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103678658B (zh) * 2013-12-24 2017-02-15 北京中科大洋科技发展股份有限公司 一种地形高程数据存储和读取的装置
CN103678657B (zh) * 2013-12-24 2017-02-15 北京中科大洋科技发展股份有限公司 一种地形高程数据存储和读取的方法
CN111652961B (zh) * 2020-07-10 2023-06-30 中国水利水电科学研究院 一种基于gpu的气象数据体绘制方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102074036A (zh) * 2010-12-07 2011-05-25 中国地质大学（武汉） 基于gpu加速的体数据动态剖切方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102074036A (zh) * 2010-12-07 2011-05-25 中国地质大学（武汉） 基于gpu加速的体数据动态剖切方法

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
An Atmospheric Visual Analysis and Exploration System;Yuyan Song et al;《IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS》;20061031;第12卷(第5期);全文 *
Joe Kniss et al.Volume Rendering Multivariate Data to Visualize Meteorological Simulations: A Case Study.《Joint EUROGRAPHICS - IEEE TCVG Symposium on Visualization》.2002,全文.
Joe Kniss et al.Volume Rendering Multivariate Data to Visualize Meteorological Simulations: A Case Study.《Joint EUROGRAPHICS- IEEE TCVG Symposium on Visualization》.2002,全文. *
Yuyan Song et al.An Atmospheric Visual Analysis and Exploration System.《IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS》.2006,第12卷(第5期),全文.
三维地质结构模型的切割分析技术及方法;陈国良 等;《计算机工程》;20071031;第33卷(第20期);全文 *
周琳娜 等.大规模数据快速体绘制方法的设计与实现.《计算机工程与设计》.2009,第30卷(第23期),全文.
大规模数据快速体绘制方法的设计与实现;周琳娜 等;《计算机工程与设计》;20091231;第30卷(第23期);全文 *
易会战.体绘制算法研究与气象可视化系统实现.《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》.2003,(第1期),全文. *
陈国良 等.三维地质结构模型的切割分析技术及方法.《计算机工程》.2007,第33卷(第20期),全文.

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106845477A (zh) * 2016-12-30 2017-06-13 武汉联影医疗科技有限公司 基于多个重建图像的感兴趣区域建立方法及其装置
CN106845477B (zh) * 2016-12-30 2020-07-28 武汉联影医疗科技有限公司 基于多个重建图像的感兴趣区域建立方法及其装置

Also Published As

Publication number Publication date
CN102364522A (zh) 2012-02-29

Similar Documents

Publication Publication Date Title
CN104090827B (zh) 集成二、三维平台的管线数据同步及联动方法
CN110738721A (zh) 基于视频几何分析的三维场景渲染加速方法及系统
CN103606192A (zh) 一种基于三维虚拟地球的风场可视化展示方法
CN102364522B (zh) 一种基于全球划分的大规模气象体数据绘制方法
CN106383828A (zh) 一种绘图方法及电子设备
CN103714213B (zh) 一种复杂地理信息的并行绘制方法
Westerteiger et al. Spherical Terrain Rendering using the hierarchical HEALPix grid
CN111127603B (zh) 动画生成方法、装置、电子设备及计算机可读存储介质
Zhang et al. An efficient dynamic volume rendering for large-scale meteorological data in a virtual globe
GB2622666A (en) Three-dimensional virtual ecological environmental visualization integration and optimization system for large region
CN104657333A (zh) 基于gpu的动态二维矢量场流线可视化算法
CN115953541A (zh) 一种四叉树lod地形生成方法、装置、设备及存储介质
She et al. An efficient method for rendering linear symbols on 3D terrain using a shader language
CN103675763B (zh) 一种基于硬件加速的矢量化vts系统雷达视频绘制方法
CN116363324A (zh) 一种态势地图二三维一体化渲染方法
Li et al. Batched trajectory compression algorithm based on hierarchical grid coordinates
Deng et al. Visualization of vector data on global scale terrain
CN112730743A (zh) 空气质量模式预报数据的交互方法及装置
CN106815359B (zh) 基于移动gis平台的一种矢量地图立体边的生成方法
CN115423917B (zh) 一种全球三维风场实时绘制方法及系统
CN110738719A (zh) 一种基于视距分层优化的Web3D模型渲染方法
CN117237503B (zh) 一种地理要素数据加速渲染及装置
CN103985150B (zh) 一种图像处理方法及电子设备
Wang et al. 3D Scene Management Method Combined with Scene Graphs.
Yang et al. Large scale terrain real-time rendering on GPU using double layers tile quad tree and cuboids bounding error metric

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20190507

Address after: 430000 Wuhan Donghu Development Zone, Wuhan, Hubei Province, No. 11 Dongxin Road, Wuhan Overseas Students Pioneer Park, 2 floors A and F

Patentee after: Wuhan Zhao Tu Science and Technology Ltd.

Address before: 430074 No. 388 Lu Lu, Hongshan District, Hubei, Wuhan

Patentee before: China University of Geosciences (Wuhan)

TR01 Transfer of patent right