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

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

Info

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
Application number
CN201110325178XA
Other languages
English (en)
Other versions
CN102364522B (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 Zhaotu Technology Co 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
Expired - Fee Related 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)中，采用八叉树剖分法根据与每一可见初始分块对应的再次剖分次数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对该可见初始分块进行剖分。
CN 201110325178 2011-10-24 2011-10-24 一种基于全球划分的大规模气象体数据绘制方法 Expired - Fee Related 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 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)

* Cited by examiner, † Cited by third party
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)

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

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 (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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