基于块内插值的多分辨率体绘制方法
【专利摘要】本发明公开了一种基于块内插值的多分辨率体绘制方法。其包括以下步骤：数据划分、分块细节水平选择、建立二级索引和体绘制。本发明的有益效果是：本发明的基于块内插值的多分辨率体绘制方法利用块内边界插值方法，避免了引入额外数据量和额外的数据预处理，消除了多分辨率体绘制的分块效应，并且利用GPU完成插值，使本发明的体绘制方法更加简单高效。
【专利说明】基于块内插值的多分辨率体绘制方法
【技术领域】
[0001]本发明属于体绘制方法【技术领域】，尤其涉及一种基于块内插值的多分辨率体绘制方法。
【背景技术】
[0002]体绘制是一种直接由三维数据场产生屏幕上二维图像的技术。数字图像对应的是描述数据元素的颜色和光强的二维阵列，这些元素成为像素，同理，一个三维数据场可以用一个具有相应值的三维阵列来描述，这些值称为体素。类似于数字图像的二维光栅，可以把体数据场看为一个三维光栅。一个典型的三维数据场是医学图像三维数据场，由CT(计算机断层成像)或MRI(核磁共振)扫描获得一系列的医学图像切片数据，把这些切片数据按照位置和角度信息进行规则化处理，然后就形成一个三维空间中由均匀网格组成的规则的数据场，网格上的每个节点为一个体素，描述了对象的密度等属性信息。体绘制以这种体素为基本操作单位，计算出每个体素对显示图像的影响。体绘制技术最大的优点是可以探索物体的内部结构，可以描述非常定形的物体，如肌肉，烟云等，而面绘制在这些方面比较弱。缺点是数据存储量大，计算时间较长。体绘制形成的图像一般是半透明的图像，颜色一般是人工指定的伪彩色。体绘制首先需要对数据进行分类处理，不同类别赋予不同的颜色和不透明度值，然后根据空间中视点和体数据的相对位置确定最终的成像效果。体绘制常用的算法有光线投射法，足迹表法，错切变形法，三维纹理贴图法等。体绘制技术是一种用于三维数据体显示的重要方法，通过对空间三维物体的直接展现，它能够准确表达出三维体数据中的数据分布、空间关系等特点。同时，通过旋转放大等一系列的交互手段，它能够形象直观的展示出数据内部的结构层次关系，因此，体绘制技术已经广泛应用于包括机械工程、地质勘探、医学以及流体物理在内的多个领域之中。随着数据探测和数据存储技术的发展，现在所用于体绘制的体数据的大小急速增加。由于体绘制技术需要预先将整个体数据读入显存，因此当数据量的大小超过显存容量的时候，体绘制就会受到限制。虽然可以采用将数据分块?弓Mt，然后分开读入内存分开绘制的方法解决内存不足的限制，但是该方法会导致频繁地硬盘到显存的读写操作，严重影响绘制的实时效率。在这种情况下，多分辨率体绘制技术应运而生。多分辨率体绘制技术采用了多分辨率表示的思想，按某种方式对数据体进行分块，为每个分块赋予不同的分辨率级别，并据此对其进行不同程度的抽样处理，以此减小绘制过程中载人显存的总体数据量。这种方式既能够在一定程度上降低用于绘制的整体数据量，又能最大限度的保证数据体内信息的完整，保证绘制质量。但是，由于在数据处理过程中采用了分块处理，同时，在生成数据纹理的时候，对不同大小的块进行了重新的整合1因此，生成的数据纹理和数据本身的空间位置关系不一样。在光线投影法中，沿着一条光线对数据进行采样时，就可能会导致采样点选取错误，从而在最终的绘制结果中出现边界效应。现有的解决边界效应的方法，大多是通过增加额外的数据量来保证采样过程的正确性。但是当数据量很大，内存限制极为严重的时候，通过增加额外数据量的方法就可能会导致最终的数据量超过内存容量限制，使得绘制变得困难。同时，该方法会增加数据预处理的复杂度。
【发明内容】
[0003]为了解决以上问题，本发明提出了一种基于块内插值的多分辨率体绘制方法。
[0004]本发明的技术方案是:一种基于块内插值的多分辨率体绘制方法，包括以下步骤:
[0005]S1.依据选取的块大小将数据体划分为多个大小相等的分块；
[0006]S2.对步骤SI中得到的每一个分块选择合适的分块细节水平，具体包括以下步骤:
[0007]S21.将步骤SI中得到的每一个分块进一步细分为多个子块；
[0008]S22.计算步骤S21中每个子块中所有数据点的均值，表示为:
【权利要求】
1.一种基于块内插值的多分辨率体绘制方法，其特征在于，包括以下步骤: S1.依据选取的块大小将数据体划分为多个大小相等的分块； S2.对步骤SI中得到的每一个分块选择合适的分块细节水平，具体包括以下步骤: S21.将步骤SI中得到的每一个分块进一步细分为多个子块； S22.计算步骤S21中每个子块中所有数据点的均值，表示为:
其中，BVgi表示第i个子块的均值，Sj为当前子块内第j个点的标量值，Ifflax为当前算法中所能使用的最高分块细节水平，I表示当前用于分析的细节水平； S23.根据步骤S22中得到的每个子块中所有数据点的均值，计算每个子块的均方差，表示为:
其中，Bi表示当前分块中第i个子块的均方差； S24.对当前分块的所有子块进行综合分析，计算在当前分块细节水平之下当前分块的平均均方差，表示为:
其中，V1表示在当前细节水平下当前分块的平均均方差，Smax表示当前分块数据的最大值，Smin表示当前分块数据的最小值； S25.定义满足Is厂avgi I >Bi的所有点为奇异点，计算在当前分块细节水平下每个子块内奇异点所占比例，并找出所有的奇异子块，计算奇异子块在当前分块内所占比例T1 ； S26.建立均方差和奇异子块比例列表，引入用户设定的门限值，计算当前分块的最终分辨率级别，表示为:
Ζ = {Ξ/|^<7:(}η{3/|7；<7；}
Im = min(L), 其中，Tu表示用户设定的均方差门限，Ts表示用户设定的奇异子块比例门限，L表示所有满足门限的分块细节水平，Im表示当前分块最终使用的分块细节水平； S3.创建二级索引，将分块数据整合为空间矩形，具体包括以下步骤: S31.对原始数据体组织进行处理，形成元祖体； S32.将所有低级分块组合为最高级分块，计算最高级分块的累计个数； S33.根据步骤S32中的最高级分块累计个数计算空间矩形大小； S34.根据步骤S33中的空间矩形大小，申请打包数据所需空间，并利用分级坐标方法，形成打包数据； S4.根据步骤S3中的空间矩形，在GPU片元上利用光线投射法进行体绘制，具体包括以下步骤:S41.按照采样间隔At对元祖体数据进行采样，并用该采样点处的光强度以及光衰减系数表征数据体中对应段的光强度和光衰减系数，则光学深度表示为:
其中，表示光线从O位置传递到t位置的累计衰减程度，k(i.At)表示光衰减程度，I表不米样点序号；S42.根据元祖体与打包数据的映射关系，计算采样点位置的颜色值Ci与不透明度Ai,表示为:
Ai = 1-θ^α,Δ?)Δ?
Ci = c (i.Δ t) Δ t ；S43.将所有点发出的光线进行融合，并在GPU片元着色器上利用块内边界插值方法得到绘制结果，表不为:
2.如权利要求1所述的基于块内插值的多分辨率体绘制方法，其特征在于:所述步骤S32中计算最高级分块的累计个数的公式具体为:
其中，N1表示将所有级别低于I的分块进行组合之后级别I的累积分块个数，L1则代表将I级及以下的所有分块全部进行组合之后I级别还剩余的分块个数。
3.如权利要求1所述的基于块内插值的多分辨率体绘制方法，其特征在于:所述步骤S33根据最高级分块累计个数计算空间矩形大小具体包括以下步骤: S331.设定空间矩形三个方向的边长分别为x、y、z，并将其初始值均设为I; S332.判断当前空间矩形是否已经满足空间要求； S333.若满足，则得到计算结果； S334.若不满足，则找到三条边长中的最小值，将其增加1，重复步骤S332。
4.如权利要求1所述的基于块内插值的多分辨率体绘制方法，其特征在于:所述步骤S34根据空间矩形大小，申请打包数据所需空间，并利用分级坐标方法，形成打包数据具体包括以下步骤: S341.初始化分级坐标，将其所有级别均设置为O; S5342.从当前的分辨率级别列表中的最高级别中取出一个分块，若最高级别列表为空，则从次闻级别列表中获取；S343.根据当前的分级坐标计算该分块存放的起始点位置； S344.将该分块写入步骤S343中得到的对应位置； S345.将分级坐标中对应级别的值增加一，并进行进位检查； S346.判断分辨率级别列表是否为空； S347.若为空，则操作结束； S348.若不为空，则继续步骤S342。
5.如权利要求4所述的基于块内插值的多分辨率体绘制方法，其特征在于:所述步骤S343根据当前的分级坐标计算该分块存放的起始点位置具体包括以下步骤: S3431.根据分级坐标中的最高级别坐标Flmax,计算当前分块应放置的大块位置ORG，计算公式表示为:
ORGz = Flfflax/(x.y).2lmax
ORGy = (Flmax% (x.y))/x.2lmax ；
ORGx = (Flmax% (x.y))% xi2lmax S3432.取出Γ级坐标F1,，其中Γ比步骤S3431中分块级别低一级； S3433.若步骤S3431中分块级别为O级，则操作结束； S3434.若步骤S3431中分块级别不为O级，则计算当前分块的起始点坐标，具体包括以下步骤: S34341.若F1, = O,则 ORG 不变； S34342.若匕,=1，则ORGx 增加 2〃 ； S34343.若F1, =2，则 ORGy 增加 21,； S34344.若F1, =3，则 0RGx，ORGy 增加 2〃 ；
S34345.若F1, = 4，则 ORGz 增加 21,；
S34346.若F1, =5，则 ORGx, ORGz 增加 21,；
S34347.若F1, =6，则 ORGy, ORGz 增加 21,；
S34348.若F1, =7，则 ORGx, ORGy, ORGz 增加 21,。
6.如权利要求1所述的基于块内插值的多分辨率体绘制方法，其特征在于:所述步骤S42中元祖体与打包数据的映射关系的确定方法具体包括以下步骤: S421.确定采样点所在分块位置，并计算该分块的起始点位置； S422.根据分块原始点位置与采样点位置，计算当前采样点的块内偏移； S423.根据当前分块在打包数据中的起始点位置及经过转换后的块内偏移，得到该点在打包数据中的实际位置。
7.如权利要求1所述的基于块内插值的多分辨率体绘制方法，其特征在于:所述步骤S43中块内边界插值方法具体包括以下步骤: S431.根据当前采样点的纹理坐标计算其块内相对坐标； S432.根据相对坐标判断其是否属于边界区域； S433.若采样点不属于边界区域，则直接取出其纹理值； S434.若采样点属于边界区域，则根据块内已知点的值估算该采样点的纹理值，具体包括以下步骤: S4341.设点A属于边界区域，记为V(A)； S4342.将点A沿Y方向分别平移相同距离，得到点B和点C，使点B处于边界上，点C处于边界内，分别记为V(B)和V (C)； S4343.根据点B和点C的值，计算得到点A的值，表示为:
V (A) = 2ν ⑶-V (c)。
【文档编号】G06T15/08GK104200511SQ201410427177
【公开日】2014年12月10日 申请日期:2014年8月27日 优先权日:2014年8月27日
【发明者】鲁才, 曹琛, 胡光岷 申请人:电子科技大学