CN113487661B - 矢量场数据的可视化方法、装置、设备和可读存储介质 - Google Patents

矢量场数据的可视化方法、装置、设备和可读存储介质 Download PDF

Info

Publication number
CN113487661B
CN113487661B CN202110662577.9A CN202110662577A CN113487661B CN 113487661 B CN113487661 B CN 113487661B CN 202110662577 A CN202110662577 A CN 202110662577A CN 113487661 B CN113487661 B CN 113487661B
Authority
CN
China
Prior art keywords
vector field
field data
data
grid
streamline
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
CN202110662577.9A
Other languages
English (en)
Other versions
CN113487661A (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.)
Beijing Daoda Tianji Technology Co ltd
Original Assignee
Beijing Daoda Tianji Technology Co ltd
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 Beijing Daoda Tianji Technology Co ltd filed Critical Beijing Daoda Tianji Technology Co ltd
Priority to CN202110662577.9A priority Critical patent/CN113487661B/zh
Publication of CN113487661A publication Critical patent/CN113487661A/zh
Application granted granted Critical
Publication of CN113487661B publication Critical patent/CN113487661B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/40—Analysis of texture
    • G06T7/41—Analysis of texture based on statistical description of texture
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20081—Training; Learning
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20084—Artificial neural networks [ANN]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Probability & Statistics with Applications (AREA)
  • Multimedia (AREA)
  • Image Generation (AREA)

Abstract

本公开的实施例提供了一种矢量场数据的可视化方法、装置、设备和计算机可读存储介质。所述方法包括获取矢量场数据；将所述矢量场数据的拓扑结构转化为预设的拓扑结构；之前/之后/同时，对所述矢量场数据进行过滤，提取所述矢量场数据的特征信息；所述特征信息包括矢量的大小、方向和位置信息；根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据。以此方式，可以拓扑结构任意复杂的数据场，同时大幅提高了计算效率。

Description

矢量场数据的可视化方法、装置、设备和可读存储介质
技术领域
本公开的实施例一般涉及数据拓扑领域，并且更具体地，涉及一种矢量场数据的可视化方法、装置、设备和计算机可读存储介质。
背景技术
矢量是一种既有大小又有方向的量纲，矢量场的数据一般都具有复杂的拓扑关系和较高的维数，如温度、密度、压强等标量，速度、重力等矢量，还有涡流张力等张量。这些物理量往往定义在离散的不规则网格的结点、中心或表面上。根据空间数据处理的特点和可视化的基本技术，将矢量场数据的几何图形表示方法通常包括点场数据表示、线场数据表示和面场数据等。
但上述传统的矢量场数据的可视化只是针对较少数据量或较小范围数据内的可视化表现较好，针对大范围三维矢量场数据的情况，数据处理及可视化处理的时间会成几何倍的增加，无法满足实时可视化的需求。
发明内容
根据本公开的实施例，提供了一种矢量场数据的可视化方案。
在本公开的第一方面，提供了一种矢量场数据的可视化方法。该方法包括：
获取矢量场数据；
将所述矢量场数据的拓扑结构转化为预设的拓扑结构；之前/之后/同时，对所述矢量场数据进行过滤，提取所述矢量场数据的特征信息；所述特征信息包括矢量的大小、方向和位置信息；
根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据。
进一步地，所述将所述矢量场数据的拓扑结构转化为预设的拓扑结构包括：
通过数据组织方法将所述矢量场数据的拓扑结构转化为预设的拓扑结构；所述数据组织方法包括六面体体元或四面体体元的数据组织方法。
进一步地，所述根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据包括：
通过点图标法和/或矢量线法对所述已转化为预设拓扑结构的矢量场数据进行基于几何形状的映射，得到可视化的矢量和/或矢量场映射图；和/或
通过线积分卷积的方法对所述已转化为预设拓扑结构的矢量场数据进行基于纹理的映射，得到可视化的纹理映射图。
进一步地，所述通过线积分卷积的方法对所述已转化为预设拓扑结构的矢量场数据进行基于纹理的映射，得到可视化的纹理映射图包括：
根据所述矢量场数据的特征信息，生成纹理图像；
将所述预设形式的矢量场数据网格化，生成网格数据；所述网格数据包括n个网格点，其中，每个网格点分别对应所述纹理图像的一个像素点；所述n为大于等于1的正整数；
根据所述网格点和所述像素点的对应关系，将所述纹理图像映射到所述网格数据上；
对所述网格数据中的每个网格点分别进行流线跟踪，得到每个网格点的流线；
对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值；
生成可视化的纹理映射图。
进一步地，所述对所述网格数据中的每个网格点分别进行流线跟踪包括：
对所述网格数据中的每个网格点，分别通过变长的四阶Runge-Kutta方程进行流程跟踪。
进一步地，所述对所述网格数据中的每个网格点分别进行流线跟踪，得到每个网格点的流线；对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值包括：
通过并行的方式对所述网格数据中的每个网格点进行流线跟踪，得到每个网格点的流线；
通过并行的方式对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值，生成可视化的纹理映射图。
进一步地，所述对所述网格数据中的每个网格点分别进行流线跟踪，得到每个网格点的流线；对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值包括：
通过并行的方式分别对每个网格点进行流线跟踪和卷积，生成可视化的纹理映射图。
在本公开的第二方面，提供了一种矢量场数据的可视化装置。该装置包括：
获取模块，用于获取矢量场数据；
处理模块，用于将所述矢量场数据的拓扑结构转化为预设的拓扑结构；之前/之后/同时，对所述矢量场数据进行过滤，提取所述矢量场数据的特征信息；所述特征信息包括矢量的大小、方向和位置信息；
显示模块，用于根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据。
在本公开的第三方面，提供了一种电子设备。该电子设备包括：存储器和处理器，所述存储器上存储有计算机程序，所述处理器执行所述程序时实现如以上所述的方法。
在本公开的第四方面，提供了一种计算机可读存储介质，其上存储有计算机程序，所述程序被处理器执行时实现如根据本公开的第一方面的方法。
本申请实施例提供的矢量场数据的可视化方法，通过获取矢量场数据；将所述矢量场数据的拓扑结构转化为预设的拓扑结构；之前/之后/同时，对所述矢量场数据进行过滤，提取所述矢量场数据的特征信息；所述特征信息包括矢量的大小、方向和位置信息；根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据，实现了大范围三维矢量场的实时可视化，同时通过线积分卷积算法大幅提高了计算效率。
应当理解，发明内容部分中所描述的内容并非旨在限定本公开的实施例的关键或重要特征，亦非用于限制本公开的范围。本公开的其它特征将通过以下的描述变得容易理解。
附图说明
结合附图并参考以下详细说明，本公开各实施例的上述和其他特征、优点及方面将变得更加明显。在附图中，相同或相似的附图标记表示相同或相似的元素，其中：
图1示出了根据本公开的实施例的矢量场数据的可视化方法的流程图；
图2示出了根据本公开的实施例的基于箭头的矢量映射图；
图3示出了根据本公开的实施例的基于流线的矢量场映射图；
图4示出了根据本公开的实施例的矢量场数据的可视化装置的方框图；
图5示出了能够实施本公开的实施例的示例性电子设备的方框图。
具体实施方式
为使本公开实施例的目的、技术方案和优点更加清楚，下面将结合本公开实施例中的附图，对本公开实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本公开一部分实施例，而不是全部的实施例。基于本公开中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的全部其他实施例，都属于本公开保护的范围。
另外，本文中术语“和/或”，仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。
图1示出了根据本公开实施例的矢量场数据的可视化方法100的流程图。方法100包括：
S110，获取矢量场数据。
所述矢量场数据可以为本地预先存储的，也可以为通过与服务器（本公开的执行主体）连接的电子设备利用有线或无线方式获取的。
在一些实施例中，所述矢量场数据（原始数据）来自于数值计算或测量的结果。
进一步地，由于计算或测量结果的多样性，所以所述矢量场数据的类型通常比较复杂。即，所述矢量场数据具有的拓扑结构复杂和数据量庞大的特点。
S120，将所述矢量场数据的拓扑结构转化为预设的拓扑结构；之前/之后/同时，对所述矢量场数据进行过滤，提取所述矢量场数据的特征信息。
因为矢量场数据通常比较复杂，不利于后续步骤的进行，因此，需要先对获取的矢量场数据进行预处理。即，将所述矢量场数据的拓扑结构转化为预设的拓扑结构
在一些实施例中，采用多规则的有效过滤方法对所述矢量场数据进行预处理。
具体地，针对所述矢量场数据拓扑结构复杂的特点，采用六面体体元或四面体体元的数据组织方法，将复杂多样的拓扑结构转换为预设的拓扑结构。所述预设的拓扑结构为根据需求预先设定的便于后续处理的拓扑结构。
进一步地，所述六面体体元数据组织方法，用于将规则的和结构化的不规则矢量数据转换为预设的拓扑结构；
所述四面体体元的数据组织方法，用于将上述六面体体元数据组织方法不能转换的矢量数据转化为预设的拓扑结构。
进一步地，所述四面体体元的数据组织适用于拓扑结构任意复杂的数据场，但是对于规则的和结构化的不规则矢量数据的转换效果不如所述六面体体元数据组织方法。因此，对于规则的和结构化的不规则矢量数据的拓扑结构转换通常采用六面体体元数据组织方法。
针对数据量庞大的特点，采用特征过滤的方法对所述矢量场数据进行过滤，提取出所述矢量场数据中的特征信息，过滤掉无用信息（例如温度、密度、压强等标量和/或涡流张力等张量）。所述矢量场的特征特性信息包括后续映射（计算）步骤中需要用到的信息，包括矢量的大小、方向和/或位置等。
进一步地，所述特征过滤方法的规则，可根据应用场景进行预先设定，以便根据预先设定的规则对所述矢量场数据进行过滤，得到所需的矢量场特征信息。
在一些实施例中，所述特征过滤法包括：特征规则制定、特征学习、特征比对等，用于将对比结果中不符合特征的数据进行剔除。
需要说明的是，上述对所述矢量场数据进行过滤的步骤，可在步骤将所述矢量场数据的拓扑结构转化为预设的拓扑结构之前/之后/同时进行。
S130，根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据。
矢量场数据映射是矢量场可视化的核心，其目的是将已转化为预设拓扑结构的矢量数据转化为可通过图形或图像予以显示的几何数据。
在一些实施例中，通过点图标法和/或矢量线法对所述已转化为预设拓扑结构的矢量场数据进行基于几何形状（几何定位）的映射，得到可视化的矢量和/或矢量场映射图。
具体地，所述点图标法包括：
对于每一采样点，用具有大小和方向的图标映射采样点处矢量的大小和方向，常见的有箭头、锥体、有向线段等。所述采样点为用户根据需求预先设定的。
进一步地，运用箭头对矢量场进行映射时，如图2所示，用箭头的长短来表示矢量的大小，用箭头的方向来表示矢量的方向，并且可以用箭头的颜色来表示矢量场内的另一场值。
所述矢量线法包括：
在矢量场中，线上所有质点的瞬时速度都与之相切的线称为场线，在流体力学中，速度矢量场中的场线称为流线。流线与欧拉描述方法相联系。欧拉方法是设法在空间中的每一点上描述出流体随时间变化的状况，因为在不同时刻将有不同的流体质点经过空间某固定点，采用速度矢量来描写固定点上流体运动的变化规律，如图3所示。
流体质点的运动规律在数学上可表示为下列矢量形式：
Figure 847027DEST_PATH_IMAGE001
；
其中，所述V为速度矢量；
所述r为矢径；
所述t为时间，积分时作为常数。
在直角坐标系中为：
Figure 295326DEST_PATH_IMAGE002
其中，x、y、z、t为欧拉变数；
当x、y、z固定，t进行改变时，上式中的函数代表空间中某固定点速度随时间的变化规律。
流线是同一时刻不同质点所组成的曲线，用于描述该时刻不同流体质点的运动方向。通过数值积分的方法得到流线，流线的方程为：
Figure 606221DEST_PATH_IMAGE003
其中，d为速度矢量；
在一些实施例中，为了提高计算精度，采用二阶Runge-Kutta方法计算积分。
在一些实施例中，通过线积分卷积的方法对所述已转化为预设拓扑结构的矢量场数据进行基于纹理的映射，得到可视化的纹理映射图。
纹理空间在图像空间进行，具有图像空间的连续性，可生成具有图像空间分辨率的细致图形。基于纹理的可视化图像看表现出一定的几何形状，并且通过颜色的有序排列，表达出一定的方向信息，克服了传统的基于颜色的映射方法无法揭示方向的缺点。
在本公开中，采用线积分卷积的方法进行矢量映射。线积分卷积用滤波器沿流线卷积白噪声图像，能够可视化流线的方向性，反映整个矢量场的结构。同时克服了箭头、流线等可视化图像的混乱现象。即，能够更好的反应矢量场的特征。
具体地：
a，根据所述矢量场的特征信息绘制一幅纹理图像；
b，将所述矢量场网格化，生成网格数据；所述网格数据包括n个网格点，其中，每个网格点分别对应所述纹理图像的一个像素点；所述n为大于等于1的正整数；
c，根据所述网格点和所述像素点的对应关系，将所述纹理图像映射到所述网格数据上；
d，对所述网格数据中的每个网格点，分别通过变长的四阶Runge-Kutta方程进行流程跟踪，得到每个网格点的流线；
对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值；
具体地，
线积分卷积使用一维低通滤波器k（w）作为卷积核沿流线卷积自噪声图像合成纹理。输出纹理的每个像素值均通过积分卷积得到。
基于该像素沿矢量正、反方向对称积分得到流线，将流线上所有像素对应静输入噪声值按卷积核参与卷积，结果作为输出纹理的像素值（颜色值）。输入的白噪声图像各点不相关，而线积分卷积纹理在流线方向上相关，可以显示矢量场的方向信息。
每个像素点的纹理强度（颜色值），可通过如下公式计算得到：
Figure 798168DEST_PATH_IMAGE004
其中，所述
Figure 143699DEST_PATH_IMAGE005
表示流线；
所述s表示弧长；
所述
Figure 497320DEST_PATH_IMAGE006
表示弧长初始值；
所述
Figure 295512DEST_PATH_IMAGE007
表示积分增量；
线积分卷积纹理在点
Figure 25570DEST_PATH_IMAGE008
的强度为
Figure 225608DEST_PATH_IMAGE009
；
所述
Figure 484551DEST_PATH_IMAGE010
为归一化参数；
所述T为白噪声；
所述k为滤波器；
所述l为滤波器长度。
进一步地，采用盒式卷积核函数计算，则上式的离散形式为：
Figure 770038DEST_PATH_IMAGE011
其中，所述
Figure 38209DEST_PATH_IMAGE012
为流线上各离散点噪声纹理值；
所述n为流线上的离散点个数；
根据计算得到的每一个像素点的颜色值对图像进行显示。
在一些实施例中，由于每个网格点的流线跟踪运算是相互独立的，因此各个网格点的流线跟踪，可基于GPU的并行计算能力并行执行。即，步骤对所述网格数据中的每个网格点分别进行流线跟踪，得到每个网格点的流线，可并行执行。
在一些实施例中，由于网格点的流线卷积计算也是相互独立的，因此也可以基于GPU的并行计算能力并行处理。即，步骤对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值，可并行处理。
在一些实施例中，由于每个网格点的流线跟踪和卷积计算都是相互独立的，而且两个运算是连续的。即，对每个网格点都是先流线跟踪，然后根据流线跟踪得到的流线进行卷积计算。因此，可以将这两个运算看成一体，利用GPU的并行计算能力并行处理，加速大范围矢量场可视化的效果。
根据本公开的实施例，实现了以下技术效果：
本公开通过多规则的有效过滤方法对原始数据进行过滤，减少冗余或其它数据对可视化的干扰。同时，在本公开中每个网格点中的各种运算（流线卷积）相互独立，保证了流线跟踪和卷积计算之间互不干扰且连续进行，大幅提高了计算效率，可适用于拓扑结构任意复杂的数据场。
需要说明的是，对于前述的各方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本公开并不受所描述的动作顺序的限制，因为依据本公开，某些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于可选实施例，所涉及的动作和模块并不一定是本公开所必须的。
以上是关于方法实施例的介绍，以下通过装置实施例，对本公开所述方案进行进一步说明。
图4示出了根据本公开的实施例的一种矢量场数据的可视化装置400的方框图。如图4所示，装置400包括：
获取模块410，用于获取矢量场数据；
处理模块420，用于将所述矢量场数据的拓扑结构转化为预设的拓扑结构；之前/之后/同时，对所述矢量场数据进行过滤，提取所述矢量场数据的特征信息；所述特征信息包括矢量的大小、方向和位置信息；
显示模块430，用于根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，所述描述的模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
图5示出了可以用来实施本公开的实施例的电子设备500的示意性框图。如图所示，设备500包括中央处理单元（CPU）501，其可以根据存储在只读存储器（ROM）502中的计算机程序指令或者从存储单元508加载到随机访问存储器（RAM）503中的计算机程序指令，来执行各种适当的动作和处理。在RAM503中，还可以存储设备500操作所需的各种程序和数据。CPU 501、ROM 502以及RAM 503通过总线504彼此相连。输入/输出（I/O）接口505也连接至总线504。
设备500中的多个部件连接至I/O接口505，包括：输入单元506，例如键盘、鼠标等；输出单元507，例如各种类型的显示器、扬声器等；存储单元508，例如磁盘、光盘等；以及通信单元509，例如网卡、调制解调器、无线通信收发机等。通信单元509允许设备500通过诸如因特网的计算机网络和/或各种电信网络与其他设备交换信息/数据。
处理单元501执行上文所描述的各个方法和处理。例如，在一些实施例中，方法可被实现为计算机软件程序，其被有形地包含于机器可读介质，例如存储单元508。在一些实施例中，计算机程序的部分或者全部可以经由ROM 502和/或通信单元509而被载入和/或安装到设备500上。当计算机程序加载到RAM 503并由CPU 501执行时，可以执行上文描述的方法的一个或多个步骤。备选地，在其他实施例中，CPU 501可以通过其他任何适当的方式（例如，借助于固件）而被配置为执行方法。
本文中以上描述的功能可以至少部分地由一个或多个硬件逻辑部件来执行。例如，非限制性地，可以使用的示范类型的硬件逻辑部件包括：场可编程门阵列（FPGA）、专用集成电路（ASIC）、专用标准产品（ASSP）、芯片上系统的系统（SOC）、负载可编程逻辑设备（CPLD）等等。
用于实施本公开的方法的程序代码可以采用一个或多个编程语言的任何组合来编写。这些程序代码可以提供给通用计算机、专用计算机或其他可编程数据处理装置的处理器或控制器，使得程序代码当由处理器或控制器执行时使流程图和/或框图中所规定的功能/操作被实施。程序代码可以完全在机器上执行、部分地在机器上执行，作为独立软件包部分地在机器上执行且部分地在远程机器上执行或完全在远程机器或服务器上执行。
在本公开的上下文中，机器可读介质可以是有形的介质，其可以包含或存储以供指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合地使用的程序。机器可读介质可以是机器可读信号介质或机器可读储存介质。机器可读介质可以包括但不限于电子的、磁性的、光学的、电磁的、红外的、或半导体系统、装置或设备，或者上述内容的任何合适组合。机器可读存储介质的更具体示例会包括基于一个或多个线的电气连接、便携式计算机盘、硬盘、随机存取存储器（RAM）、只读存储器（ROM）、可擦除可编程只读存储器（EPROM或快闪存储器）、光纤、便捷式紧凑盘只读存储器（CD-ROM）、光学储存设备、磁储存设备、或上述内容的任何合适组合。
此外，虽然采用特定次序描绘了各操作，但是这应当理解为要求这样操作以所示出的特定次序或以顺序次序执行，或者要求所有图示的操作应被执行以取得期望的结果。在一定环境下，多任务和并行处理可能是有利的。同样地，虽然在上面论述中包含了若干具体实现细节，但是这些不应当被解释为对本公开的范围的限制。在单独的实施例的上下文中描述的某些特征还可以组合地实现在单个实现中。相反地，在单个实现的上下文中描述的各种特征也可以单独地或以任何合适的子组合的方式实现在多个实现中。
尽管已经采用特定于结构特征和/或方法逻辑动作的语言描述了本主题，但是应当理解所附权利要求书中所限定的主题未必局限于上面描述的特定特征或动作。相反，上面所描述的特定特征和动作仅仅是实现权利要求书的示例形式。
尽管已经采用特定于结构特征和/或方法逻辑动作的语言描述了本主题，但是应当理解所附权利要求书中所限定的主题未必局限于上面描述的特定特征或动作。相反，上面所描述的特定特征和动作仅仅是实现权利要求书的示例形式，本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。

Claims (7)

1.一种矢量场数据的可视化方法，其特征在于，包括：
获取矢量场数据；
将所述矢量场数据的拓扑结构转化为预设的拓扑结构；之前/之后/同时，对所述矢量场数据进行过滤，提取所述矢量场数据的特征信息；所述特征信息包括矢量的大小、方向和位置信息；
根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据；
所述将所述矢量场数据的拓扑结构转化为预设的拓扑结构包括：
通过数据组织方法将所述矢量场数据的拓扑结构转化为预设的拓扑结构；所述数据组织方法包括六面体体元或四面体体元的数据组织方法；
所述根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据包括：
通过点图标法和/或矢量线法对所述已转化为预设拓扑结构的矢量场数据进行基于几何形状的映射，得到可视化的矢量和/或矢量场映射图；和/或
通过线积分卷积的方法对所述已转化为预设拓扑结构的矢量场数据进行基于纹理的映射，得到可视化的纹理映射图；
所述通过线积分卷积的方法对所述已转化为预设拓扑结构的矢量场数据进行基于纹理的映射，得到可视化的纹理映射图包括：
根据所述矢量场数据的特征信息，生成纹理图像；
将预设形式的矢量场数据网格化，生成网格数据；所述网格数据包括n个网格点，其中，每个网格点分别对应所述纹理图像的一个像素点；所述n为大于等于1的正整数；
根据所述网格点和所述像素点的对应关系，将所述纹理图像映射到所述网格数据上；
对所述网格数据中的每个网格点分别进行流线跟踪，得到每个网格点的流线；
对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值；生成可视化的纹理映射图。
2.根据权利要求1所述的方法，其特征在于，所述对所述网格数据中的每个网格点分别进行流线跟踪包括：
对所述网格数据中的每个网格点，分别通过变长的四阶Runge-Kutta方程进行流线跟踪。
3.根据权利要求1所述的方法，其特征在于，所述对所述网格数据中的每个网格点分别进行流线跟踪，得到每个网格点的流线；对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值包括：
通过并行的方式对所述网格数据中的每个网格点进行流线跟踪，得到每个网格点的流线；
通过并行的方式对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值，生成可视化的纹理映射图。
4.根据权利要求1所述的方法，其特征在于，所述对所述网格数据中的每个网格点分别进行流线跟踪，得到每个网格点的流线；对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值包括：
通过并行的方式分别对每个网格点进行流线跟踪和卷积，生成可视化的纹理映射图。
5.一种矢量场数据的可视化装置，其特征在于，包括：
获取模块，用于获取矢量场数据；
处理模块，用于将所述矢量场数据的拓扑结构转化为预设的拓扑结构；之前/之后/同时，对所述矢量场数据进行过滤，提取所述矢量场数据的特征信息；所述特征信息包括矢量的大小、方向和位置信息；所述将所述矢量场数据的拓扑结构转化为预设的拓扑结构包括：通过数据组织方法将所述矢量场数据的拓扑结构转化为预设的拓扑结构；所述数据组织方法包括六面体体元或四面体体元的数据组织方法；
显示模块，用于根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据；所述根据所述矢量场数据的特征信息，将已转化为预设拓扑结构的矢量场数据转化为可视化的几何数据包括：通过点图标法和/或矢量线法对所述已转化为预设拓扑结构的矢量场数据进行基于几何形状的映射，得到可视化的矢量和/或矢量场映射图；和/或通过线积分卷积的方法对所述已转化为预设拓扑结构的矢量场数据进行基于纹理的映射，得到可视化的纹理映射图；所述通过线积分卷积的方法对所述已转化为预设拓扑结构的矢量场数据进行基于纹理的映射，得到可视化的纹理映射图包括：根据所述矢量场数据的特征信息，生成纹理图像；将预设形式的矢量场数据网格化，生成网格数据；所述网格数据包括n个网格点，其中，每个网格点分别对应所述纹理图像的一个像素点；所述n为大于等于1的正整数；根据所述网格点和所述像素点的对应关系，将所述纹理图像映射到所述网格数据上；对所述网格数据中的每个网格点分别进行流线跟踪，得到每个网格点的流线；对所述每个网格点的流线进行卷积，得到与所述流线对应的像素点的颜色值；生成可视化的纹理映射图。
6.一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1～4中任一项所述的方法。
7.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现如权利要求1～4中任一项所述的方法。
CN202110662577.9A 2021-06-15 2021-06-15 矢量场数据的可视化方法、装置、设备和可读存储介质 Active CN113487661B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110662577.9A CN113487661B (zh) 2021-06-15 2021-06-15 矢量场数据的可视化方法、装置、设备和可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110662577.9A CN113487661B (zh) 2021-06-15 2021-06-15 矢量场数据的可视化方法、装置、设备和可读存储介质

Publications (2)

Publication Number Publication Date
CN113487661A CN113487661A (zh) 2021-10-08
CN113487661B true CN113487661B (zh) 2022-12-23

Family

ID=77935381

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110662577.9A Active CN113487661B (zh) 2021-06-15 2021-06-15 矢量场数据的可视化方法、装置、设备和可读存储介质

Country Status (1)

Country Link
CN (1) CN113487661B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102831626A (zh) * 2012-06-18 2012-12-19 清华大学 极地投影模式下多变量时空数据的可视化方法
CN102999936A (zh) * 2012-11-19 2013-03-27 北京中海新图科技有限公司 一种海洋流场数据的三维流线体绘制算法
CN104658027A (zh) * 2015-02-11 2015-05-27 中国海洋大学 一种面向不规则海洋流场数据的三维流线动态可视化算法
CN107451958A (zh) * 2017-08-09 2017-12-08 中国科学院海洋研究所 一种面向ArcGIS极地矢量场可视化的角度转换算法
CN108171781A (zh) * 2018-01-05 2018-06-15 北京理工大学 一种基于图标的三维多变量矢量场数据可视化方法
CN109559369A (zh) * 2018-12-04 2019-04-02 中国科学院地理科学与资源研究所 一种面向流向特征的基于角度跟踪的海流矢量场可视化方法
CN111767679A (zh) * 2020-07-14 2020-10-13 中国科学院计算机网络信息中心 时变矢量场数据的处理方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9626565B2 (en) * 2014-07-01 2017-04-18 Digitalglobe, Inc. Automated conversion of two-dimensional hydrology vector models into valid three-dimensional hydrology vector models
CN108875275B (zh) * 2018-07-18 2023-02-17 成都信息工程大学 一种基于大规模流线的矢量场实时仿真方法及系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102831626A (zh) * 2012-06-18 2012-12-19 清华大学 极地投影模式下多变量时空数据的可视化方法
CN102999936A (zh) * 2012-11-19 2013-03-27 北京中海新图科技有限公司 一种海洋流场数据的三维流线体绘制算法
CN104658027A (zh) * 2015-02-11 2015-05-27 中国海洋大学 一种面向不规则海洋流场数据的三维流线动态可视化算法
CN107451958A (zh) * 2017-08-09 2017-12-08 中国科学院海洋研究所 一种面向ArcGIS极地矢量场可视化的角度转换算法
CN108171781A (zh) * 2018-01-05 2018-06-15 北京理工大学 一种基于图标的三维多变量矢量场数据可视化方法
CN109559369A (zh) * 2018-12-04 2019-04-02 中国科学院地理科学与资源研究所 一种面向流向特征的基于角度跟踪的海流矢量场可视化方法
CN111767679A (zh) * 2020-07-14 2020-10-13 中国科学院计算机网络信息中心 时变矢量场数据的处理方法及装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Streamline Variability Plots for Characterizing the Uncertainty in Vector Field Ensembles;Florian Ferstl 等;《IEEE TRANSACTIONS ON VISUALIZATION AND COMPUTER GRAPHICS》;20160131;第22卷(第1期);第767-776页 *
三维矢量场的流动拓扑结构抽取;王康健 等;《计算机辅助设计与图形学学报》;20081130;第20卷(第11期);第1514-1520页 *
矢量场流线自适应颜色映射增强算法;秦绪佳 等;《计算机科学》;20201130;第47卷(第11期);第236-240页 *

Also Published As

Publication number Publication date
CN113487661A (zh) 2021-10-08

Similar Documents

Publication Publication Date Title
CN112785674B (zh) 纹理贴图的生成方法、渲染方法、装置、设备及存储介质
US20210110599A1 (en) Depth camera-based three-dimensional reconstruction method and apparatus, device, and storage medium
CN108229654B (zh) 神经网络卷积运算装置及方法
Hlawatsch et al. Hierarchical line integration
Schoenenberger et al. Graph-based denoising for time-varying point clouds
US10540789B2 (en) Line stylization through graphics processor unit (GPU) textures
Ferstl et al. Interactive separating streak surfaces
CN107871337B (zh) 一种超音速二维流场数据的可视化方法
CN107203962B (zh) 一种利用2d图片制作伪3d图像的方法及电子设备
CN116704088A (zh) 三维模型渲染方法和装置、设备及存储介质
CN114139104A (zh) 基于偏微分方程处理流场数据的方法、装置及电子设备
Xi et al. Super resolution reconstruction algorithm of video image based on deep self encoding learning
CN114494564A (zh) 法线贴图生成方法、设备、电子设备及存储介质
CN118154752A (zh) 基于三维离散数据对物体着色的方法、介质及设备
CN114663597A (zh) 基于归一化扩展极线几何的实时结构光重建方法及装置
Bachthaler et al. Space-time Visualization of Dynamics in Lagrangian Coherent Structures of Time-dependent 2D Vector Fields.
CN113487661B (zh) 矢量场数据的可视化方法、装置、设备和可读存储介质
Vyatkin Method of binary search for image elements of functionally defined objects using graphics processing units
US8379023B2 (en) Calculating graphical vertices
CN115311403B (zh) 深度学习网络的训练方法、虚拟形象生成方法及装置
CN110738719A (zh) 一种基于视距分层优化的Web3D模型渲染方法
CN115409949A (zh) 模型训练方法、视角图像生成方法、装置、设备及介质
Mochurad et al. Comparison of the Efficiency of Parallel Algorithms KNN and NLM Based on CUDA for Large Image Processing.
US20180158194A1 (en) Determining Optical Flow
CN110246224B (zh) 网格模型的表面去噪方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 100085 room 703, 7 / F, block C, 8 malianwa North Road, Haidian District, Beijing

Applicant after: Beijing daoda Tianji Technology Co.,Ltd.

Address before: 100085 room 703, 7 / F, block C, 8 malianwa North Road, Haidian District, Beijing

Applicant before: Beijing daoda Tianji Technology Co.,Ltd.

GR01 Patent grant
GR01 Patent grant