本发明专利技术公开一种基于张量近似的海量数据多分辨率体绘制方法，首先对原始数据进行分块处理，得到若干数据块，然后对每一个数据块进行张量分解和多分辨率处理，最后对经张量分解和多分辨率处理后的每一个数据块进行重构处理，并创建二维纹理，完成地震数据的绘制。通过采用秩截断的方法有效的过滤原始数据中的噪声，并采用因子矩阵和核心张量秩截断的方式来替代对每一次秩试探所做的张量分解，准确的对每一个数据块确定秩的大小，节省数据块最佳秩选取的时间，并根据秩截断后的得到的数据，选择每个数据块的细节水平，快速有效地降低数据的整体分辨率大小，减少处理时间，从而实现多分辨率处理，并得到比传统基于信息熵的多分辨率更好的绘制效果。
【技术实现步骤摘要】
本专利技术属于图像处理领域，具体涉及一种体绘制技术。
技术介绍
对体数据的可视化技术是一项非常常用的技术，能够被广泛应用于许多领域中， 例如：医学领域、流体物理领域、气象领域、地质勘探领域等等。由于人类对于视觉信号的刺 激最为敏感，因此可视化技术可以将含有大量信息的文件、图片、表格等转换为三维图像， 便于研究人员直观地去观察和分析它们。 由于在现实生活中，常见的三维物体都是其表面，因此，在基于信息熵的模型中， 常采用表面表示的方式去绘制一个三维体。然而很多时候，人们所关心的恰恰是物体的内 部结构。例如在医学领域，人们更希望通过可视化技术观察器官内部的组织是否有病变；在 地质勘探领域，人们更希望通过可视化技术观察地层以下的各种地质结构。因此此时，表面 绘制就有了很大的局限性，它并不能帮助研究人员获得他们需要的信息。所以，人们需要一 个能够观测物体内部结构信息的可视化技术--体绘制技术。 体绘制技术的本质，就是通过技术手段，用一个三维的标量数据生成一张二维的 图片，并将其显示在电脑屏幕上。体绘制能够对物体的内部结构进行透视，使观察者能够一 览物体的整体而不仅仅是表面，因此具有广泛的应用价值。 然而，随着数据探测技术的发展，人们获得的数据量成几何倍数增长，加之现在全 球进入了互联网时代，互联网上的数据更是呈爆炸式增长。因此，用于体绘制的体数据量越 来越大。由于计算机寻址空间的限制，体数据的大小很容易就超过了计算机的显存甚至内 存大小。 于是，人们开始使用并行绘制来解决大数据量的问题。通过对体数据的不同部分 进行分布式计算的方法来完成对整个体数据的绘制。但是，由于在体数据中，我们所关心的 目标数据往往很小，采用并行体绘制显得得不偿失。 在这种情况下，产生了多分辨率体绘制技术。多分辨率体绘制技术是通过将数据 体划分为不同的分块，每一个分块赋予不同的分辨率来实现压缩数据量和减少绘制点数。 每个分块的分辨率又被称为分块的LOD(level of detail，细节水平）。因此，多分辨率技 术的关键即为如何确定每一个分块的细节水平。 通过计算分块的信息熵（方差）来确定一个分块的细节水平是一种十分常用的方 法。它通过计算分块的方差来得到分块的同质性。通常认为同质性高的分块，所包含的信 息量小，因此不需要较高的分辨率。因而，方差越大的分块，同质性越低，说明其蕴含的信息 量越大，因此其需要较高的分辨率。 对于某些数据，传统的依靠计算分块的数值特征并不能很好的实现多分辨率处 理。例如，在地质勘探中常用的地震数据。由于地震数据具有信噪比低，变化剧烈，同质区域 较少的特点，因此每个分块的方差都很高，多分辨率处理后的分块都普遍具有较高的分辨 率。同时由于地震数据中存在大量的噪声和其他无用信息，因而高的方差并不能代表其具 有较高的信息量。所以，单纯的多分辨率处理并不能有效的降低地震数据的数据量。另一 方面，理论上的地震数据具有十分显著的微结构特征，不同的结构值之间的轮廓十分清晰， 例如层位、断层。但是，由于地震数据中存在着大量的噪声和无用信息，导致实际应用中通 过传感器探测到的数据在局部是混乱的，结构之间界限是模糊的，从而导致研究人员不能 方便地从中分辨出他们所关心的结构。因此需要将地震数据中人们所关心的结构特征提取 出来，过滤掉噪声和无用信息，然后将其可视化，提供直观的清晰的显示。此时，基于信息熵 的多分辨率处理已经不能满足有效降低数据量和提取结构特征的要求了。
技术实现思路
本专利技术为解决的上述技术问题，提出一种基于张量近似的海量数据多分辨率体绘 制方法，通过采用秩截断的方法有效的过滤原始数据中的噪声，并通过自适应秩截断的方 法，准确的对每一个分块确定秩的大小，根据秩截断后的数据，选择分块的细节水平，从而 实现多分辨率处理，并得到比传统基于信息熵的多分辨率更好的绘制效果。 本专利技术采用的技术方案是：， 首先对原始数据进行分块处理，得到若干数据块，然后对每一个数据块进行张量分解和多 分辨率处理，最后对经张量分解和多分辨率处理后的每一个数据块进行重构处理，并创建 二维纹理，完成地震数据的绘制。 进一步地，具体包括以下步骤： Sl :对原始数据进行分块，得到若干数据块； S2 :对步骤Sl得到的每一个数据块进行张量分解； S3 :对步骤S2得到的每一个数据块进行多分辨率处理； S4 :对步骤S3得到的每一个数据块进行重构，并创建二维纹理，根据重构得到的 数据块进行地震数据绘制。 进一步地，所述步骤S2包括以下分步骤： S21 :对每个数据块根据各自的初始秩进行张量分解，得到因子矩阵和对应核心张 量，并设置初始秩截断参数秩为1 ; S22:对数据块的因子矩阵和对应的核心张量根据当前秩截断参数秩进行秩截断， 得到秩截断后的因子矩阵和对应的核心张量； S23 :根据步骤S22得到的因子矩阵和对应的核心张量，进行该数据块的重构，并 计算重构误差； S24 :判断步骤S23得到的重构误差是否满足收敛条件，若是，则进行步骤S25,否 贝1J，进行步骤S26 ; S25 :输出数据块在当前秩截断参数秩得到的因子矩阵和对应的核心张量； S26 :当前秩截断参数秩进行自加一操作，重复步骤S22至S25,得到每个数据块在 各自对应当前秩截断参数秩下的因子矩阵和对应的核心张量。 更进一步地，步骤S21所述的对每个数据块的秩进行初始化，该初始秩根据每个 数据块的分块尺寸确定。 更进一步地，所述步骤S22计算重构误差，具体为： 其中，e表示重构误差，A表示原始张量，2表示重构后的近似张量，I I I |F表示矩 阵的范书。 更进一步地，步骤S23所述的判断步骤S22得到的重构误差是满足收敛条件具体 为：判断当前秩截断参数秩是否小于或等于数据块的初始化秩R ;或者，判断当前秩截断参 数秩的重构误差e是否小于或等于当前秩截断参数秩的归一化重构误差?；;或者，判断当 前秩截断参数秩的重构误差e是否满足下式： 其中，e'表示上一秩截断参数秩的重构误差，Tp表示当前秩截断参数秩的重构误 差的提升值。 进一步地，所述步骤S3具体包括以下分步骤： S31 :根据数据块的分块个数为2的k次幂，得到分辨率为k个级别； S32 :根据初始的秩截断参数秩为1以及步骤S21中的每个数据块的秩的初始化值 R，得到秩的取值范围为，共有R个取值； S33 :根据（R_l+lV(k+l)，得到商为a，余数为b，则秩的分组为b组包括a+Ι个秩， 剩下的k+1-b组包括a个秩，随机组合得到分组排列； S34 :将秩的取值按照从大到小进行排序，并根据步骤S33得到的分组排列进行分 组。 进一步地，所述步骤S4根据下式对数据块进行重构： 其中，表示核心张量位于（Γι，r2, r3)位置的值，表示矩阵U(n)的第r n列的 列向量。 本专利技术的有益效果：本专利技术的一种基于张量近似的海量数据多分辨率体绘制方 法，通过采用秩截断的方法有效的过滤原始数据中的噪声并通过自适应秩截断的方法，采 用因子矩阵和核心张量秩截断的方式来替代对每一次秩试探所做的张量分解，准确的对每 一个数据块确定秩的大小，节省数据块最佳秩选取的时间，并根据秩截断后的得到的数据， 选择每个数据块的细节水平，快速有效地降低本文档来自技高网...
【技术保护点】
一种基于张量近似的海量数据多分辨率体绘制方法，其特征在于，首先对原始数据进行分块处理，得到若干数据块，然后对每一个数据块进行张量分解和多分辨率处理，最后对经张量分解和多分辨率处理后的每一个数据块进行重构处理，并创建二维纹理，完成地震数据的绘制。
【技术特征摘要】
【专利技术属性】
技术研发人员：鲁才，张力彬，曹琛，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。