CN102364522A - 一种基于全球划分的大规模气象体数据绘制方法 - Google Patents
一种基于全球划分的大规模气象体数据绘制方法 Download PDFInfo
- Publication number
- CN102364522A CN102364522A CN201110325178XA CN201110325178A CN102364522A CN 102364522 A CN102364522 A CN 102364522A CN 201110325178X A CN201110325178X A CN 201110325178XA CN 201110325178 A CN201110325178 A CN 201110325178A CN 102364522 A CN102364522 A CN 102364522A
- Authority
- CN
- China
- Prior art keywords
- volume data
- meteorological
- piecemeal
- initial
- subdivision
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 abstract 1
- 238000007781 pre-processing Methods 0.000 abstract 1
- 238000012163 sequencing technique Methods 0.000 abstract 1
- 238000009877 rendering Methods 0.000 description 6
- 230000000007 visual effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
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)中,采用八叉树剖分法根据与每一可见初始分块对应的再次剖分次数n对该可见初始分块进行剖分。
与现有技术相比,本发明提供的基于全球划分的大规模气象体数据绘制方法具有以下优点:基于全球划分的策略对规则气象体数据进行绘制,剥离了与数据的相关性,可以将不同时区下、不同投影下的气象体数据集成绘制,能处理大规模的气象体数据,并可以在保证图像质量的同时大大地提高了绘制的速度;根据视点计算球面下视域内的可见分块,然后只对可见分块进行绘制,降低了对图形硬件资源的占用,同时提高了绘制效率;若视点移动,可动态地得到视点与初始分块的距离,并实时计算出分块剖分的级数,达到实时绘制的目的;利用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)中,采用八叉树剖分法根据与每一可见初始分块对应的再次剖分次数n对该可见初始分块进行剖分。
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 true CN102364522A (zh) | 2012-02-29 |
CN102364522B CN102364522B (zh) | 2013-02-13 |
Family
ID=45691088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110325178 Expired - Fee Related CN102364522B (zh) | 2011-10-24 | 2011-10-24 | 一种基于全球划分的大规模气象体数据绘制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102364522B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678657A (zh) * | 2013-12-24 | 2014-03-26 | 北京中科大洋科技发展股份有限公司 | 一种地形高程数据存储和读取的方法 |
CN103678658B (zh) * | 2013-12-24 | 2017-02-15 | 北京中科大洋科技发展股份有限公司 | 一种地形高程数据存储和读取的装置 |
CN111652961A (zh) * | 2020-07-10 | 2020-09-11 | 中国水利水电科学研究院 | 一种基于gpu的气象数据体绘制方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106845477B (zh) * | 2016-12-30 | 2020-07-28 | 武汉联影医疗科技有限公司 | 基于多个重建图像的感兴趣区域建立方法及其装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102074036A (zh) * | 2010-12-07 | 2011-05-25 | 中国地质大学(武汉) | 基于gpu加速的体数据动态剖切方法 |
-
2011
- 2011-10-24 CN CN 201110325178 patent/CN102364522B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102074036A (zh) * | 2010-12-07 | 2011-05-25 | 中国地质大学(武汉) | 基于gpu加速的体数据动态剖切方法 |
Non-Patent Citations (5)
Title |
---|
JOE KNISS ET AL: "Volume Rendering Multivariate Data to Visualize Meteorological Simulations: A Case Study", 《JOINT EUROGRAPHICS - IEEE TCVG SYMPOSIUM ON VISUALIZATION》 * |
YUYAN SONG ET AL: "An Atmospheric Visual Analysis and Exploration System", 《IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS》 * |
周琳娜 等: "大规模数据快速体绘制方法的设计与实现", 《计算机工程与设计》 * |
易会战: "体绘制算法研究与气象可视化系统实现", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
陈国良 等: "三维地质结构模型的切割分析技术及方法", 《计算机工程》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678657A (zh) * | 2013-12-24 | 2014-03-26 | 北京中科大洋科技发展股份有限公司 | 一种地形高程数据存储和读取的方法 |
CN103678657B (zh) * | 2013-12-24 | 2017-02-15 | 北京中科大洋科技发展股份有限公司 | 一种地形高程数据存储和读取的方法 |
CN103678658B (zh) * | 2013-12-24 | 2017-02-15 | 北京中科大洋科技发展股份有限公司 | 一种地形高程数据存储和读取的装置 |
CN111652961A (zh) * | 2020-07-10 | 2020-09-11 | 中国水利水电科学研究院 | 一种基于gpu的气象数据体绘制方法 |
CN111652961B (zh) * | 2020-07-10 | 2023-06-30 | 中国水利水电科学研究院 | 一种基于gpu的气象数据体绘制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102364522B (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6906707B2 (ja) | 地図における地理的領域のグリッドを識別するためのシステムおよび方法 | |
CN104090827B (zh) | 集成二、三维平台的管线数据同步及联动方法 | |
US20150123968A1 (en) | Occlusion render mechanism for point clouds | |
KR100956416B1 (ko) | 3차원 지리 정보 클라이언트 장치와 서버 장치 및 이를포함한 3차원 지리 정보 시스템 | |
CN106383828B (zh) | 一种绘图方法及电子设备 | |
CN103606192A (zh) | 一种基于三维虚拟地球的风场可视化展示方法 | |
CN104504742A (zh) | 一种基于电子沙盘系统的三维军标标绘方法 | |
CN102364522B (zh) | 一种基于全球划分的大规模气象体数据绘制方法 | |
CN110070616A (zh) | 存储器、基于gis平台的统计数据呈现方法、装置和设备 | |
CN103714213A (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 | |
Masood et al. | High‐performance virtual globe GPU terrain rendering using game engine | |
CN112597260A (zh) | 空气质量模式预报数据的可视化方法及装置 | |
CN115953541A (zh) | 一种四叉树lod地形生成方法、装置、设备及存储介质 | |
CN103675763B (zh) | 一种基于硬件加速的矢量化vts系统雷达视频绘制方法 | |
Deng et al. | Visualization of vector data on global scale terrain | |
Qiao et al. | A rapid visualization method of vector data over 3D terrain | |
CN112730743A (zh) | 空气质量模式预报数据的交互方法及装置 | |
CN106815359B (zh) | 基于移动gis平台的一种矢量地图立体边的生成方法 | |
KR20100060194A (ko) | 항법 장치에서 포인트 기반 렌더링 방법 및 장치 | |
CN115423917B (zh) | 一种全球三维风场实时绘制方法及系统 | |
CN103985150B (zh) | 一种图像处理方法及电子设备 | |
Sheng et al. | MCGIM-based model streaming for realtime progressive rendering |
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 ZHAOTU TECHNOLOGY CO.,LTD. Address before: 430074 No. 388 Lu Lu, Hongshan District, Hubei, Wuhan Patentee before: CHINA University OF GEOSCIENCES (WUHAN CITY) |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130213 |
|
CF01 | Termination of patent right due to non-payment of annual fee |