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